Photochemical transformation of residential wood combustion emissions: dependence of organic aerosol composition on OH exposure

Residential wood combustion (RWC) emits large amounts of gaseous and particulate organic aerosol (OA). In the atmosphere, the emission is transformed via oxidative reactions, which are under daylight conditions driven mainly by hydroxyl radicals (OH). This continuing oxidative ageing produces secondary OA and may change the healthand climate-related properties of the emission. However, it is not well known how the composition of RWCoriginated OA changes as the function of OH exposure. In this work, emissions from two modern residential logwood combustion appliances were photochemically aged in an oxidation flow reactor (OFR) with various OH exposure levels, reaching up to 6× 1011 s cm−3 (equivalent to 1 week in the atmosphere). Gaseous organic compounds were analysed by proton transfer reaction time-of-flight mass spectrometry (PTR-ToF-MS), while particulate OA was analysed online by a high-resolution soot particle aerosol mass spectrometer (SP-HR-ToF-AMS) and offline by in situ derivatization thermal desorption–gas chromatography–time-of-flight mass spectrometry (IDTD-GC-ToF-MS). Photochemical reactions increased the mass of particulate organic carbon by a factor of 1.3–3.9. The increase in mass took place during the first atmospheric equivalent day of ageing, after which the enhancement was independent of the extent of photochemical exposure. However, ageing increased the oxidation state of the particulate OA linearly throughout the assessed range, with 1H : C/1O : C slopes between −0.17 and −0.49 in van Krevelen space. Ageing led to an increase in acidic fragmentation products in both phases, as measured by the IDTD-GC-ToF-MS for the particulate and PTR-ToFMS for the gaseous phase. For the gaseous organic compounds, the formation of small carbonylic compounds combined with the rapid degradation of primary volatile organic compounds such as aromatic compounds led to a continuous increase in both the O : C and H : C ratios. Overall, the share of polycyclic aromatic compounds (PACs) in particles degraded rapidly during ageing, although some oxygensubstituted PACs, most notably naphthaldehydic acid, increased, in particular during relatively short exposures. Similarly, the concentrations of particulate nitrophenols rose extensively during the first atmospheric equivalent day. During continuous photochemical ageing, the dominant transformation mechanisms shifted from the initial gas-phase functionalization/condensation to the transformation of the particuPublished by Copernicus Publications on behalf of the European Geosciences Union. 6358 A. Hartikainen et al.: Photochemical transformation of residential wood combustion emissions late OA by further oxidation reactions and fragmentation. The observed continuous transformation of OA composition throughout a broad range of OH exposures indicates that the entire atmospheric lifetime of the emission needs to be explored to fully assess the potential climate and health effects of RWC emissions.

Abstract. Residential wood combustion (RWC) emits large amounts of gaseous and particulate organic aerosol (OA). In the atmosphere, the emission is transformed via oxidative reactions, which are under daylight conditions driven mainly by hydroxyl radicals (OH). This continuing oxidative ageing produces secondary OA and may change the health-and climate-related properties of the emission. However, it is not well known how the composition of RWCoriginated OA changes as the function of OH exposure. In this work, emissions from two modern residential logwood combustion appliances were photochemically aged in an oxidation flow reactor (OFR) with various OH exposure levels, reaching up to 6 × 10 11 s cm −3 (equivalent to 1 week in the atmosphere). Gaseous organic compounds were analysed by proton transfer reaction time-of-flight mass spectrometry (PTR-ToF-MS), while particulate OA was analysed online by a high-resolution soot particle aerosol mass spectrometer (SP-HR-ToF-AMS) and offline by in situ derivatization thermal desorption-gas chromatography-time-of-flight mass spectrometry (IDTD-GC-ToF-MS). Photochemical reactions increased the mass of particulate organic carbon by a factor of 1.3-3.9. The increase in mass took place during the first atmospheric equivalent day of ageing, after which the enhancement was independent of the extent of photochemical exposure. However, ageing increased the oxidation state of the particulate OA linearly throughout the assessed range, with H : C/ O : C slopes between −0.17 and −0.49 in van Krevelen space. Ageing led to an increase in acidic fragmentation products in both phases, as measured by the IDTD-GC-ToF-MS for the particulate and PTR-ToF-MS for the gaseous phase. For the gaseous organic compounds, the formation of small carbonylic compounds combined with the rapid degradation of primary volatile organic compounds such as aromatic compounds led to a continuous increase in both the O : C and H : C ratios. Overall, the share of polycyclic aromatic compounds (PACs) in particles degraded rapidly during ageing, although some oxygensubstituted PACs, most notably naphthaldehydic acid, increased, in particular during relatively short exposures. Similarly, the concentrations of particulate nitrophenols rose extensively during the first atmospheric equivalent day. During continuous photochemical ageing, the dominant transformation mechanisms shifted from the initial gas-phase functionalization/condensation to the transformation of the particu-

Introduction
Biomass combustion is a major source of atmospheric particulate matter (PM) and is considered the main anthropogenic source of organic matter and the third largest contributor of black carbon (BC) emissions globally (Klimont et al., 2017). The use of wood fuels in small-scale residential settings is a main source for ambient organic aerosol (OA) in many parts of the world. For example, residential wood combustion (RWC) has been identified as a major source of ambient air fine particles in several European cities, where its relative contribution has been estimated to further increase in the future, while PM emissions from other sources, such as industry and traffic, are decreasing (Denier Van Der Gon et al., 2015;Klimont et al., 2017). The amount and contents of the RWC emissions depend greatly on combustion conditions, which are generally affected by the combustion procedure, fuel, and appliance technology (Bhattu et al., 2019;Nuutinen et al., 2014;Orasche et al., 2013;. In logwood-fired appliances, there is also a strong variation in the emissions during the different combustion phases of batches, with ignition producing the highest organic emissions (Bhattu et al., 2019;. However, the highest black carbon concentrations are emitted during the flaming phase, while the char burnout phase typically emits large amounts of carbon monoxide but low, mainly inorganic particulate emissions . Combustion conditions also affect the emissions of many toxic compounds, such as polycyclic aromatic compounds (PACs) (Kim et al., 2013;Orasche et al., 2013), and are consequently strongly linked with the adverse health effects of the emissions (Bølling et al., 2009;Kanashova et al., 2018;Kasurinen et al., 2018).
Many of the main organic species in fresh, not atmospherically aged wood smoke are connected to the composition of the wood fuel, such as levoglucosan, a common biomass burning marker, and lignin degradation products such as methoxyphenols and their derivatives (Elsasser et al., 2013;McDonald et al., 2000;Orasche et al., 2013). In addition, the gaseous organic emission from RWC contains hundreds of organic gaseous compounds (OGCs) (Bhattu et al., 2019;Bruns et al., 2017;Hartikainen et al., 2018;Hatch et al., 2017;McDonald et al., 2000), a majority of which can be classified as volatile organic compounds (VOCs) (Hatch et al., 2017), meaning that under ambient conditions they exist purely in gaseous phase. Notably, RWC is an impor-tant anthropogenic source for OGCs with high potential for secondary particulate organic aerosol (SOA) formation. The most potent SOA-precursor compounds include aromatic hydrocarbons and oxygenated species, such as phenolic and furanoic compounds (Bruns et al., 2016;Hartikainen et al., 2018), and removal of these compounds via either improved combustion conditions or for example catalytic cleaning has been shown to be efficient in lowering the SOA potential of RWC emissions (Czech et al., 2017;Pieber et al., 2018). Furthermore, the gaseous emissions contain high amounts of carbonyls, such as formaldehyde and acetaldehyde, with adverse health effects (Reda et al., 2015;USEPA, 2019).
OA has an atmospheric lifetime of approximately 1 week (Hodzic et al., 2016), during which its chemical composition and potential environmental and health effects are likely to transform extensively. In daylight conditions, hydroxyl radicals (OH) dominate this ageing process, where the oxidation of OGCs forms a variety of functionalized products. Some of these oxidized secondary organic species partition into the particulate phase, resulting in an enhancement of ambient air particulate organic matter concentrations (Robinson et al., 2007). While aromatic compounds are noted as the main SOA precursors from RWC, the complete pathway for SOA formation and the final SOA yields of complex OGC mixtures under different atmospheric conditions remain unclear (Bruns et al., 2016;Hartikainen et al., 2018;Hatch et al., 2017;McFiggans et al., 2019). Recent experiments on RWC exhaust estimate the mass of particulate OA to increase by a factor of 1.6-5.3 within approximately 1 d of photochemical ageing (Bertrand et al., 2017;Bruns et al., 2015b;Grieshop et al., 2009;Heringa et al., 2011;Tiitta et al., 2016). In addition, heterogeneous oxidation reactions of the particulate matter may be significant during the atmospheric transformation of RWC emissions: it has been reported that only a minority of the primary particulate OA (POA) from biomass burning remains unreacted after a few hours of atmospheric ageing (Hennigan et al., 2011;Tiitta et al., 2016). Photochemical ageing also decomposes hydrocarbonaceous PACs (i.e. polycyclic aromatic hydrocarbons, PAHs), which cause adverse health effects due to their carcinogenic properties. On the other hand, atmospheric ageing may lead to formation of even more toxic, oxygen-or nitrogen-substituted PACs (Finlayson-Pitts and Pitts, 2000;Keyte et al., 2013;Miersch et al., 2019). For instance, these substituted PACs have been reported to cause a substantial part of the particle-induced mutagenicity in Beijing, with a contribution of only 8 % relative to the concentration of hydrocarbon PACs (Wang et al., 2011). Similarly, oxidation in nitrogen oxide (NO x )-rich conditions can produce nitrophenols, which are known to be harmful to plant growth and human health (Harrison et al., 2005) and have been identified as an important constituent of light-absorbing organic matter ("brown carbon") (Moise et al., 2015;Zhang et al., 2011), thereby affecting atmospheric radiative forcing.
The photochemical ageing of RWC emissions has previously been studied mainly in smog chambers (Bertrand et al., 2017;Bruns et al., 2015a, b;Heringa et al., 2011;Tiitta et al., 2016), where the ageing was monitored as a batch process from fresh emission to up to 1 atmospheric equivalent day of exposure (eqv. d; assuming an ambient average OH concentration of 10 6 molec cm −3 ; Prinn et al., 2005). As an alternative, oxidation flow reactors (OFRs) with continuous sample flow have been increasingly utilized in combustion emission studies (Bruns et al., 2015a;Czech et al., 2017;Pieber et al., 2018). To achieve similar to or higher oxidant exposures than chambers in shorter residence times, OFRs have been used with high ozone concentrations together with high-intensity, low-wavelength UV lamps to generate OH-radical concentrations orders of magnitudes higher than those of smog chambers (Kang et al., 2007;Ihalainen et al., 2019;Simonen et al., 2017). Thus, OFRs enable measurements with better temporal resolution, which is a benefit when assessing the ageing of aerosols from dynamic sources, such as batchwise logwood combustion.
In this study, we aim to determine how atmospheric ageing changes the composition of organic aerosol, emitted by logwood-fired stoves, as the function of OH exposure. The atmospheric transformation of emissions from two RWC appliances fired with spruce and beech logwood was simulated using the photochemical emission ageing flow tube reactor (PEAR OFR; Ihalainen et al., 2019) at different exposure levels ranging up to a week of atmospheric age. The transformation of emissions and the formation of related secondary organic emissions were monitored with a comprehensive setup (Fig. 1), including analysis of gas-phase by proton transfer reactor time-of-flight mass spectrometry (PTR-ToF-MS) and investigation of the particulate phase online by aerosol mass spectrometry (SP-HR-ToF-AMS) and offline by targeted gas chromatography mass spectrometry and thermal-optical analyses. Together, these analyses enable the assessment of the OA chemical composition in both bulk and molecular level, from fresh emission to up to exposures representative of long-range transported smoke.

Experimental conditions
Experiments were conducted in the ILMARI laboratory of the University of Eastern Finland (http://www.uef.fi/ilmari, last access: 20 May 2020) with the experimental setup shown in Fig. 1. The two combustion appliances used were modern stoves with improved air intakes. First, a heat-storing masonry heater (Hiisi 4, Tulikivi Ltd., Finland) representing the typical modern logwood combustion technology utilized in northern Europe, was fired with spruce logwood. The combustion procedure in the masonry heater consisted of three 2.5 kg batches (35 min combustion time) of spruce logwood representing kiln-dried fuel (5 % moisture content, S-5 %), after which there was a 25 min char burning period prior to two 45 min batches with moist (22 % moisture content, S-22 %) spruce logwood. Second, a non-heat-retaining chimney stove (Aduro 9.3, Denmark) representing central European modern logwood stoves was fired with both beech and spruce logs. The chimney stove experiments consisted of five 2 kg batches (combustion time 40-55 min) of beech logwood (17 % H 2 O, B-17 %), followed by two 50 min batches of S-22 %. Beech was used in these experiments because it is the most common firewood used in central Europe, while spruce is used both in northern and central Europe. For ignition, 150 g of dry kindling was placed on the top of the first batch in the cold furnace. Each batch was divided into three parts: ignition, flaming, and burnout phase. The ignition phase was determined to last from the beginning of the batch to the moment of batch maximum flue gas CO 2 concentration, and the burnout phase began from the moment when the CO concentration started to elevate and remained at a high level until the end of the batch (Fig. 2). The modified combustion efficiency (MCE) as a function of time was calculated from primary flue-gas CO 2 and CO concentrations as CO 2 /( CO 2 + CO).
The exhaust was sampled from the stack with a 10 µm precut cyclone. The sample was diluted with a combination of a porous tube diluter and an ejector diluter (Dekati FPS ejector, Finland) and had a dilution ratio (DR) of 40-150 (Table 1) when fed to the PEAR OFR. In the chimney stove  experiments, an additional ejector diluter with a DR of 8 (Dekati DI-1000, Finland) was placed before the secondary online aerosol instruments. Compared to the no-ageing experiments, a higher DR was selected for the ageing experiments in response to the expected increase in particulate matter in the PEAR OFR during photochemical ageing. Measured concentrations were corrected for the dilution and normalized to stoichiometric dry flue-gas by multiplication with the stoichiometric correction factor (SCF) of Eq. (1) based on the background CO 2 concentration (CO 2,bg ), the CO 2 concentration in the diluted sample downstream of the PEAR OFR (CO 2,sec ), and the fact that the dry flue gas of wood combustion with no excess air contains 20.2 % CO 2 .
In addition, the diluted concentrations were normalized to 13 % flue-gas oxygen content.

Use of the PEAR OFR
The PEAR OFR  was used to continuously age the sample stream. In the setup, the extent of photon flux and consequential photochemical ageing were controlled by adjusting the voltage of the 254 nm UV lamps. The total flow through the PEAR OFR was 60 L min −1 , and the residence time within the PEAR was 139 s. In addition to the 55 L min −1 sample flow, humified and purified air was introduced to the PEAR OFR to obtain a relative humidity of 45 ± 5 %. In the PEAR OFR, the photolysis of the ozone forms hydroxyl radicals via Reactions (R1) and (R2). Thus, the extent of photochemical ageing in the PEAR OFR depends on the photon flux inside the reactor and the introduction of OH precursors, namely the externally fed O 3 and H 2 O vapours.
In the ageing experiments, ozone and butanol-d9 mixed with the main sample flow before the reactor inlet ( Fig. 1) resulted in concentrations of 2-11 ppm of O 3 (Table 1) and 80-200 ppb of butanol-d9 in the PEAR OFR before initiation of photochemical ageing. The OH exposure (OH exp ) was determined continuously based on the butanol-d9 concentrations downstream of the PEAR according to the method presented by Barmet et al. (2012). The OH exp for clean air in the PEAR prior to the exhaust input ranged from 8.2 × 10 10 to 1.6 × 10 12 molec cm −3 s, which is equivalent to 1-18 d in an ambient atmosphere (eqv. d) with an estimated 24 h average global OH concentration of 1 × 10 6 molec cm −3 (Prinn et al., 2005). The alternate reaction pathways in the PEAR OFR, namely exposure to photolysis (flux 254 nm,exp ), excited (O( 1 D) exp ), atomic oxygen (O( 3 P) exp ), and ozone (O 3,exp ), were assessed using the OFR exposure estimation equations of Peng et al. (2016) for OFRs with 254 nm lamps. Both the OH exp and the alternate reaction pathways were affected by the external OH reactivity (OHR ext ) (Li et al., 2015) of the gases in the sample fed to the PEAR OFR. The OHR ext was calculated with Eq. (2) from the concentrations (c i ) of the primary gaseous compounds (Table S1) and the butanol-d9 in the PEAR OFR, with their OH reaction rate constants k i .
Due to the differences in the emission concentrations during logwood combustion, the OHR ext and consequently OH exp vary also within a batch (Fig. S1). Average OHR ext was in the range 130-1300 s −1 with the highest contributions from CO, NO, and unsaturated hydrocarbons (Fig. S2). Particulate wall losses inside the PEAR OFR were minimized by conductive stainless-steel walls, laminar flow, and a relatively small surface-to-volume ratio (2.28 m 2 : 139 L) . The particulate losses were estimated to be approximately 6 % based on the loss of elemental carbon (EC) in the PEAR OFR, determined from the difference in the thermal-optical EC concentrations measured upstream and downstream of the reactor. The fate of the gaseous organic compounds capable of irreversible condensation under the present experimental conditions, that is, the low-volatility organic compounds (LVOCs), was estimated based on  by considering three possible depletion pathways: condensation onto particles, reactions with OH radicals, and condensation onto walls. To estimate the fraction of LVOCs condensing onto particle phase, condensation sinks ( Fig. S3) were calculated according to Lehtinen et al. (2003) using an average of the particle size distributions up-and downstream of the PEAR OFR. See Sect. S4.2 in the Supplement for further information on LVOC fate estimation.

Offline filter sampling and analysis
PM 1 filter samples were collected on Teflon (PTFE, Pall Corporation, P/N R2PJO47) and quartz fibre (QF; Pall Corporation, Tissuquartz) filters simultaneously from primary and secondary exhaust at a 10 L min −1 flow rate, following the methodology presented by Sippula et al. (2009). A pre-impactor (Dekati PM-10 impactor) was used to separate the particles with aerodynamic diameters of less than 1 µm (PM 1 ). For the masonry heater, samples were collected separately from the S-5 % combustion (100 min collection time) and S-22 % combustion (90 min collection time). From the chimney stove, two sample pairs were collected from the combustion of beech: full first two batches (85 min collection time) and fourth and fifth batches excluding the last 15 min (85 min collection time). The third chimney stove collection consisted of the combustion of two full batches of moist spruce (100 min collection time). The Teflon filters were weighted before and after sample collection to determine the total PM 1 mass of the emission. The amount of organic (OC) and elemental carbon (EC) in PM 1 was determined from the QF filters by using a thermal-optical carbon analyser (Sunset Laboratory Inc.) following the protocol NIOSH5040 (NIOSH, 1999). In situ derivatization thermal desorption-gas chromatography-time-of-flight mass spectrometry (IDTD-GC-ToF-MS) (Orasche et al., 2011), was applied for targeted analysis of organic compounds in the particulate phase from the QF filters. Nonpolar and polar compounds were identified and quantified using mixtures of isotopically labelled internal standards and calibration standards; see the Supplement, Sect. S5.1, for further information on the IDTD-GC-ToF-MS analysis procedure. The majority of the targeted compounds can be classified as either semior intermediately volatile, meaning that under the present experimental conditions they can be partitioned into both gaseous and particulate phases.

Online gas-phase measurements
A Fourier transform infrared spectrometer (FTIR DX4000, Gasmet Technologies Inc.) was implemented on the stack to measure the amount of NO x , CO 2 , CO, and 27 OGCs (Table S1) in the fresh exhaust. The measured primary nonmethane OGCs were grouped into four subgroups: alkanes, oxygenated compounds, and unsaturated and aromatic hy-drocarbons. The concentrations of CO 2 , O 3 , NO x , and SO 2 were monitored at the outlet of the PEAR OFR with a tracelevel single-gas analysers (ABB CO 2 analyser, Siemens), and OGCs were measured by PTR-ToF-MS (PTR-TOF 8000, Ionicon Analytik, Innsbruck, Austria), with H 3 O + as the reagent ion and an electric field to gas number density ratio (E / N) of 130. Mass calibration was done with H 3 O + (m/z 21.02) and 1,3-diiodobenzene (m/z 203.94), which was added as a calibrant for higher m/z. The processing of the PTR-ToF-MS data was done in a manner similar to that of earlier work . Reaction rates by Cappellin et al. (2012) were used when available; for other compounds, the reaction rate with H 3 O + was assumed to be 2 × 10 9 cm 3 s −1 (Table S2).

Online particulate-phase measurements
Particle concentrations and mobility size distributions were monitored with scanning mobility particle sizers before (SMPS 3082, TSI, size range 14.6-661.2 nm) and after the PEAR OFR (SMPS 3080, TSI, size range 15.1-661.2 nm). An electrical low-pressure impactor (ELPI; 10 L min −1 , Dekati) measured the particle aerodynamic size distribution and the concentration of primary particles in the size range 18.6-5950 nm. The total particulate mass after PEAR OFR was measured with a tapered element oscillating microbalance monitor (TEOM, Model 1405, Thermo Scientific).
The composition of submicron particulate matter after the PEAR OFR was measured by soot particle aerosol mass spectrometer (SP-HR-ToF-AMS, Aerodyne Research Inc). The SP-HR-ToF-AMS dual vaporizer mode was used, with the combination of the thermal vaporizer (600 • C) and the continuous wave laser vaporizer (1064 nm) enabling the study of both nonrefractory (NR-PM) and refractory lightabsorbing submicron aerosol particles (R-PM, e.g. refractory black carbon, rBC) (Onasch et al., 2012). Standard mass-based calibrations were performed for the ionization efficiency in the SP-HR-ToF-AMS using ammonium nitrate and Regal Black (Regal 400R Pigment Black, Cabot Corp.) particles (Jayne et al., 2000;Onasch et al., 2012). SP-HR-ToF-AMS was operated in V-mode from 12 to 555 m/z, and the two vaporizer configurations were alternated every 120 s, including the particle time-of-flight (PTOF) mode (duration 20 s). The SP-HR-ToF-AMS data were analysed using the standard analysis tools SQUIRREL v1.62A and PIKA v1.22D adapted in Igor Pro 8 (Wavemetrics). The data were corrected by time-dependent background gas-phase CO 2 subtraction using the online HEPA filter measurement technique. The interactions of inorganic salts with pre-deposited carbon on the tungsten vaporizer were corrected according to Pieber et al. (2016). The CO 2 -AMS interference calibration value was 0.3 % of the NO 3 concentration determined by NO 3 NH 4 calibration and corrected via the fragmentation table according to Pieber et al. (2016). Before correction, the relative effect of the interference on O : C ratio was typically under 1 % and 5 % at maximum. The applied relative ionization efficiency (RIE) was 1.4, which agrees with the OA mass concentrations of biomass-burning emissions within the stated ±38 % uncertainty in AMS (Bahreini et al., 2009;Xu et al., 2018). The elemental analysis of the OA was conducted using the improved-ambient method (Canagaratna et al., 2015), and the average carbon oxidation state (OS C ) of the OA was estimated as OS C ≈ 2× O : C − H : C (Kroll et al., 2011).
The SP-HR-ToF-AMS OA mass spectra was further examined by positive matrix factorization (PMF) using the method described in Lanz et al. (2007) and . Similar methods have been used previously for example for the assessment of RWC generated POA and SOA in a chamber (Bruns et al., 2015a;Tiitta et al., 2016) and in an OFR (Bruns et al., 2015a) or for time-resolved analysis of RWC OA emission constituents (Elsasser et al., 2013;. The PMF Evaluation Tool v.3.05 was applied, and the standard data pretreatment process was completed based on Ulbrich et al. (2009), including the application of minimum error criteria and down-weighting weak variables as well as m/z 44 (CO + 2 ) and water-related peaks. The final four-factor PMF solution covered 98 % of the OA spectra (2.2 % residual). Additional factors did not increase the realistic physical meaning of the solution, while fewer factors were insufficient for a meaningful presentation of the data. The factor identification was confirmed by comparing the time series and mass spectra of each factor with external tracers (nitrate, sulfate, ammonium, chloride, PAH, CO + 2 , C 2 H 3 O + , C 4 H + 9 , and C 2 H 4 O + 2 ). Furthermore, the factors were compared to logwood combustion mass spectra measured by Tiitta et al. (2016) from the ageing of spruce logwood exhaust in a smog chamber. The agreement of spectra was denoted with both a coefficient of determination (R 2 ) and the angle between two mass spectra vectors (Kostenidou et al., 2009), where an angle of less than 15 • indicates a good agreement between two mass spectra.

Results and discussion
This study comprehensively characterizes the chemical properties of RWC exhaust at different atmospheric ageing times by combining extensive information gathered from gas-phase and particulate-phase chemical analyses. In this section, we first discuss the dynamic combustion conditions and the characteristics of primary emissions from the logwood stoves utilized with different fuels, which define the starting point for the ageing experiments. Next, the ageing conditions in the PEAR OFR are evaluated in order to validate the atmospheric relevance of the results. Finally, we assess the changes in the gaseous and particulate OA during the ageing process under a variety of different oxidant concentrations. The observations of changes both in bulk-and molecular level chemical composition of the aerosol demonstrate that the major transformation pathway of OA changes from initial gas-phase functionalization followed by condensation to the transformation of the particulate OA by heterogeneous oxidation reactions and fragmentation. The study shows a linear dependency between OH exposure and organic aerosol oxidation state. Furthermore, OH-exposure-dependencies of specific OA constituents, such as nitrophenols, carboxylic acids, and PACs, are established.

Combustion conditions
The average modified combustion efficiency was greater than 0.97 for the masonry heater and greater than 0.95 for the chimney stove, with lower MCE occurring mainly during the char burnout periods (Fig. S4). These values are typical for modern batchwise-operated logwood appliances (Bhattu et al., 2019;Heringa et al., 2011). In addition to the variation within a combustion batch, there were also differences between the individual batches ( Fig. 2). Notably, the combustion conditions during the first batch were distinct from later batches, because of the ignition in a cold firebox. This caused a longer ignition period (determined by the rising CO 2 concentrations), which lasted for 24 ± 5 % (masonry heater) and 35 ± 4 % (chimney stove) of the total duration of the batch, whereas in the later batches the firebox temperatures are higher, and the ignition shortened to 9 ± 3 % (masonry heater) or 8 ± 4 % (chimney stove) of the total combustion time. Furthermore, the flaming phases of the first batches were shorter than those of the following batches of dry spruce or beech fuels. The emission profiles were affected by these batchwise differences, with ignition being the period for enhanced organic emissions, whereas the flaming phase was characterized by an increase in particulate emissions consisting mainly of black carbon, as expected based on previous work . The char burnout phases with these fuels were characterized by high CO concentrations, whereas in moist spruce combustion elevated CO concentrations were measured throughout the batch, thus making the burnout phases less distinguishable.

Gaseous organic emissions
The primary organic gaseous emissions in the undiluted flue gas were measured continuously by an FTIR multicomponent analyser. Additionally, the PTR-ToF-MS measured diluted, unaged emissions during the no-ageing experiments. These datasets complemented each other because, while FTIR was calibrated for 27 typical combustion-derived compounds (Table S1), a more detailed insight into the composition of the gaseous organic phase was acquired via the PTR-ToF-MS, with which 126 different molecular formulas for OGC in the primary emissions were identified in the m/z range of 40-180 (Table S2). However, PTR-ToF-MS is unable to detect compounds with proton affinities lower than that of water (691 kJ mol −1 ), such as alkanes. Furthermore, fragmentation in the PTR-ToF-MS limits the identification of compounds with similar mass-to-charge ratios as common fragment ions, including unsaturated aliphatic compounds, such as propene (C 3 H + 6 ) at m/z 41.04 or butene (C 4 H + 8 ) at m/z 57.08. However, these compounds were detected by the FTIR.
The major OGC groups measured by PTR-ToF-MS were carbonyls, aromatic hydrocarbons (ArHCs), furans, and phenols ( Fig. 3; Table S2). In addition, unsaturated aliphatic compounds constituted a substantial fraction of the total nonmethane OGC (NM-OGC) measured using the FTIR (Fig. 4). The NM-OGC emissions based on FTIR were 42.6 ± 9.2 and 101 ± 24.6 mg C m −3 for dry and moist spruce combustion in the masonry heater, respectively, and 89.4 ± 11.3 and 148.3 ± 23.0 mg C m −3 for beech and moist spruce in the chimney stove, respectively. Thus, the lowest NM-OGC concentrations were measured from the dry spruce combustion. Moist spruce combustion produced a factor of 1.7-3.6 higher emission than dry wood, which is well in line with studies by e.g. McDonald et al. (2000), where moist fuel produced 2-4 times more NM-OGC than dry fuel combustion. The difference between concentrations was highest for oxygenated compounds (factor of 4.2 ± 1.7) and for unsaturated compounds (2.7 ± 0.6), but statistically significant for all subgroups (paired t-test p values ≤ 0.02 for all groups for the consecutive dry/moist experiments; Table S4). Differences in NM-OGC emissions were also observed between the masonry heater and the chimney stove. Emissions from moist spruce combusted in the chimney stove were higher by a factor of 1.5 (p value of 0.05) compared to the masonry heater, with a statistically significant increase in unsaturated and aromatic hydrocarbons (factors of 1.6 and 1.7, respectively, p values < 0.01) based on FTIR measurements.
The NM-OGC emissions of the first batches exhibited a distinct time-dependent behaviour in comparison to the fol-  lowing batches: the emitted concentrations always increased both after ignition and at the end of flaming phase (Fig. 2), whereas in the following batches there was a sharp emission peak at ignition, after which the concentrations declined as soon as the flaming phase began. These findings agree with those by  and are influenced by the fact that cold ignition is performed from the top of a fuel batch, while the following batches are ignited from the bottom by the glowing embers. NM-OGC emissions were lowest during the burnout period when most of the fuel was already consumed (Fig. 2). Thus, NM-OGC emissions had a reverse time profile to the CO emissions, which peak during the char burnout phase. Furthermore, the contribution of the different organic compound groups to the total NM-OGC concentration differed in relation to time (Fig. 4). For example, the importance of aromatic species increased for dry fuels (S-5 % and B-17 %) for the flaming and burnout phases, while other species were pronouncedly emitted at ignition. In addition, the composition of aromatic species measured by PTR-ToF-MS in the unaged emission depended on the phase (Fig. S5). In the masonry heater, ArHCs had the highest con-tribution during the ignition phase, but their relative share decreased during flaming phase, while the share of furanoic and phenolic compounds increased. The relative importance of furanoic and phenolic compounds in the fresh exhaust of the flaming phase from a masonry heater has been also previously established . Overall, the share of ArHC in the fresh exhaust is higher and less phase dependent for beech combustion in a chimney stove than for other experiments. The share of N-containing aromatic compounds, namely nitrophenol and nitrocresol, also increased after ignition. These findings are important also when considering the potential of SOA formation, as aromatic compounds have been observed to be the major SOA precursors in RWC emissions, and removal of ArHCs from the flue gas either by improving combustion conditions or using e.g. catalytic cleaning has been noted to decrease also the resulting SOA formation (Bruns et al., 2016;Hartikainen et al., 2018;Pieber et al., 2018). Table 2. Primary PM 1 concentrations in dry, 13 % O 2 flue gas conditions. The number concentration and geometric mean mobility diameter (GMD) were derived from a scanning mobility particle sizer and PM 1 , organic carbon, and elemental carbon from filter samples.

Particulate emissions
The average primary PM 1 mass concentrations, 33-67 mg m −3 , and the number concentrations, 3.2×10 7 cm −3 -5.4 × 10 7 cm −3 (Table 2), were on a similar level with earlier studies reporting emissions from modern logwood stoves Nuutinen et al., 2014;. The combustion of dry spruce in the masonry heater emitted 1.5 and 2 times as much PM 1 mass compared to that of moist spruce and beech, respectively, mainly because of the higher elemental carbon emissions. Moist spruce generated similar PM 1 emissions with both combustion appliances. The ratio of organic carbon to elemental carbon (OC : EC) of dry spruce combustion was very low (0.07 ± 0.02, Table 2) compared to that of moist spruce combustion (0.31 ± 0.45 for the masonry heater and 0.25 ± 0.04 for the chimney stove) and beech combustion (0.15 ± 0.04). Such low OC : EC ratios have been previously reported for emissions from modern masonry heaters operating with dry logwood Miersch et al., 2019;Nuutinen et al., 2014). The particle size distribution (Fig. S3) from dry spruce combustion was clearly distinguishable from those of other wood fuels and had a geometric mean mobility diameter (GMD) of 95.5 nm, which was considerably larger than those of other fuels (GMD 52.8-68.4 nm). The soot-dominated composition of S-5 % exhaust likely increases the GMD, because soot particles are typically present as larger agglomerates than particles of inorganic origin (ash), which mainly form ultrafine particles (Tissari et al., 2008). The size distribution and number concentration of particles in an exhaust are important not only because of their link to potential health effects but also during ageing of the exhaust, because they affect the condensation sink (CS) (Lehtinen et al., 2003) of condensable vapours during the dilution and ageing process. Thus, the particle number concentrations and size distributions affect the fate of condensable vapours and the overall OA enhancement ratios.
The thermal-optically measured elemental carbon in both the primary and the secondary exhausts correlated well with the refractory black carbon (rBC) measured by SP-HR-ToF-AMS from the secondary exhaust (R 2 = 0.74 and 0.76, re-spectively; see Fig. S6). Analogous to the elemental carbon results, rBC emission was highest during the combustion of kiln-dried spruce (49.3 ± 13.7 mg m −3 ). Dry fuel has recently been found to increase the soot and PAC emissions also for combustion of birch logwood in sauna stoves . Considering the different combustion phases, rBC emissions were highest during high-temperature flaming combustion, as previously noted also by e.g. . For moist spruce combustion in the masonry heater, the rBC concentration dropped to 24.8 ± 12.3 mg m −3 , likely because of the lower temperature and consequentially slower burn rate. Similar rBC concentrations (24.0 ± 6.0 mg m −3 ) were also measured from spruce combustion in the chimney stove. Unlike the masonry heater where the rBC concentrations were similar throughout the three dry batches, the rBC concentration from the combustion of beech in a chimney stove decreased considerably after the first batch, from 57.7 ± 7.2 to 15.6 ± 5.2 mg m −3 . While ageing had no effect on the rBC mass, it plays an important role in the formation of SOA by acting as a seed for condensation during ageing. Furthermore, soot cores composed of elemental carbon are chemically active and may enhance the photooxidation of an OA condensed onto soot agglomerates through electron transfer (Li et al., 2018).

OFR conditions
During photochemical ageing, OHR ext is an important parameter affecting OH-radical consumption and reaction pathways of organic species. For the masonry heater experiments, the average OHR ext was 130-150 s −1 for dry spruce samples and 160-220 s −1 for moist spruce. In the chimney stove, the average OHR ext was 280-300 s −1 for beech and 320-370 s −1 for moist spruce samples, except for the low-DR experiment (Exp. 5, DR of 30) where the average OHR ext reached 820 and 1300 s −1 for the beech and moist spruce samples, respectively. The OFR conditions were divided into good, risky, and bad based on the ratio of the photon flux exposure to OH exp by the definitions of Peng and Jimenez (2017) (Fig. S7). Here, the ageing conditions were defined as mainly "risky" (4 × 10 5 cm s −1 < flux 254 nm,exp / OH exp < 10 7 cm s −1 ) during all the experiments. The OHR ext varied considerably during combustion cycles, but the limit for "bad" conditions was exceeded only during ignition phases of chimney stove experiments when the OHR ext peaks above 1000 s −1 for brief periods accounting for up to 4 % of total experiment time.
In terms of OA emission, the bad conditions accounted for 3 %-7 % of chimney stove OA emissions, excluding the S-22 % combustion for the low-DR experiment where 30 % of OA was emitted under bad conditions. Evaluation of the OHR ext was limited to the externally added butanol-d9 and the compounds measured from the primary exhaust with FTIR, with NO, CO, and unsaturated hydrocarbons as the main OHR ext producers (Fig. S2). In other words, the products of later-generation oxidation or from fragmentation from the particulate phase were not considered. The average NO x concentrations input to the PEAR OFR ranged from 150 to 420 ppb, except for the low-DR experiment with concentrations of 1140 and 751 ppb for beech and spruce combustion, respectively. NO is rapidly oxidized to NO 2 with the addition of O 3 and then partially converted to particulate nitrate. The subsequent low-NO conditions in the PEAR decreased the reactions of organic peroxy radicals (RO 2 ) with NO. Based on the RO 2 chemistry model by Peng et al. (2019) for 254 nm OFRs, the high photon fluxes and high concentrations of both OH and HO 2 in the PEAR OFR result in lower RO 2 lifetime and consequently lower RO 2 isomerization rate than in atmosphere, while the importance of RO 2 + OH reaction is enhanced due to lower-than-ambient HO 2 -to-OH-ratio (Fig. S8) (Peng et al., 2019).
The average batchwise OH exposure in the PEAR OFR during photochemical ageing experiments ranged from 0.5 to 7 eqv. d depending on the applied photon flux and the OHR ext of the sample. During the experiments, the photochemical exposure varied in line with the varying OHR ext during batchwise combustion, with the highest exposure occurring during the flaming phase and lower exposure during the ignition period (Fig. S1). The extent of alternative non-OH reaction pathways during photochemical ageing in the PEAR OFR were compared to the tropospheric conditions (Sect. S4.1), and the exposures to O( 1 D) and O( 3 P) in the PEAR OFR were estimated to be similar to those in ambient conditions (Fig. S9), excluding the ignition period during the low-DR experiment, where the importance of O( 3 P) exp briefly exceeded ambient conditions. The exposure to ozone in relation to OH radicals was estimated to be lower in the PEAR OFR than in the troposphere (O 3,exp /OH exp < 10 5 ). However, our estimations were based on the initial O 3 concentrations, whereas O 3 is expected to form in the PEAR OFR during photochemical ageing as a product of the OH + OGC reactions (Carter, 1994). This may have led to the O 3,exp /OH exp being slightly underestimated.
Several of the emitted OGCs were susceptible to photolysis at 254 nm, most importantly the ArHCs, which are among the main SOA precursors from RWC. Photolysis as a degradation pathway can exceed OH reactions also in the atmo-sphere for compounds with high photolysis rates (Hodzic et al., 2015). Here, the importance of photolysis was notable, e.g. for benzene, of which more than 40 % may have degraded via photolysis (Fig. S10). However, the importance of photolysis during photochemical ageing is inversely proportional to the ratio of the OH reaction rate to the photoabsorption cross section (σ abs ), and it can be considered a minor pathway for other main OGCs (Table S7), including other aromatic species such as toluene (< 10 % of total degradation during these experiments; Fig. S10).
The detailed results of the LVOC fate calculations are presented in Sect. S4.2. Briefly, the majority of LVOCs were estimated to condense into the particulate phase in all experimental conditions (Fig. S11). Similar to earlier OFR studies (Bruns et al., 2015a;Ihalainen et al., 2019;Simonen et al., 2017), the downstream particle number concentration and condensation sinks were influenced by new particle formation ( Fig. S12) with the extent of nucleation depending strongly on the ageing conditions and the concentrations of condensable vapours. However, the distinct formation of ultrafine particles does not necessarily represent the atmospheric fate of condensable vapours because of the fasterthan-ambient oxidation, which results in higher saturation ratios and faster gas-to-particle conversion in the PEAR OFR. Gas-to-particle conversion is also affected by the dilution ratio, which was in the PEAR OFR lower than what typically exists in atmosphere. Therefore, it has to be noted that the relative fraction of the measured semivolatile species in the condensed phase are likely lower at ambient dilution ratios than in the conditions of this study.

Transformation of gaseous phase
The composition of the OGCs changed throughout the studied exposure range, and the abundance of the secondary compounds increased with an increase in photochemical exposure. Namely, the OGCs measured with PTR-ToF-MS from the aged gas phase were governed by small carbonyls and fragmentation products such as acetic acid (C 2 H 4 O 2 , m/z 61.02). However, while the small oxidized organics became increasingly dominant in the gas phase, photochemical ageing particularly decreased the share of aromatic compounds (Fig. 3). The evolution of the total chemical composition is visualized in the van Krevelen diagram, which is often used for simplified characterization of particulate OA (Heald et al., 2010) and was here extended for the investigation of the gaseous organic phase measured by PTR-ToF-MS (Fig. 5). Photochemical ageing caused several simultaneous and subsequent functionalization reactions of the organic compounds, and the increase in the average H : C ratio together with an increasing O : C ratio led to linear slopes from +0.69 to +1.0 in the van Krevelen diagram for gaseous phase. Similar, positive slopes are derivable also from previous studies of ageing of RWC OGC in smog chambers Hartikainen et al., 2018). * Behaviour of RWC emission aged in a chamber (Hartikainen et al., 2018).
Most of the OGCs in primary RWC exhaust have relatively high OH reactivities, and thus their reactions are expected to take place within the first atmosphere equivalent hours of photochemical ageing Hartikainen et al., 2018). Therefore, the secondary OGC products, including the oxidized fragmentation products from the particulate phase, dominate the total OGC in the long-range transported smoke. In addition, the loss of compounds with low H : C and O : C ratios, such as aromatic species, is important for the change in the average composition. Furthermore, the oxidation of compounds undetected by the PTR-ToF-MS, such as alkanes, may introduce secondary products with higher proton affinities and their subsequent appearance in the spectra. Conversely, ageing may also lead to the growth of compounds outside the observed mass range (m/z 40-180). Compared to previously measured changes in RWC exhaust in a chamber , ageing in the PEAR OFR led to a slightly higher increase in the H : C ratio. This difference implies an increased fragmentation which is likely the result of faster OH exposure and differences in the RO 2 chemistry (Peng et al., 2019) in the PEAR OFR compared to ageing in a smog chamber. Aromatic compounds consisting of ArHCs, phenols, furans, N-containing aromatic compounds (N-aromatics), and other oxygenated aromatics were important constituents of the primary organic gas phase and formed 37 %-39 % of the fresh emission from the masonry heater and 33 %-34 % from the chimney stove, as measured by the PTR-ToF-MS (Fig. 3). Similar shares have been measured from fresh RWC exhaust by Bruns et al. (2017) (13 %-33 %) and Hartikainen et al. (2018) (33 %-36 %). However, after 2 eqv. d of exposure, their share decreased to less than 20 % of the identified OGCs, which agrees with previously reported conversions of aromatics during photochemical exposure Hartikainen et al., 2018). Overall, the photochemical reactions of aromatics are an important source of SOA because they form products that efficiently partition into the particulate phase. However, there are large differences between the conversion efficiencies of aromatic compound groups. While ArHCs comprise approximately half of the aromatic compounds in fresh gaseous exhaust, their share grows to over 70 % in 5-7 eqv. d, while the share of oxygenated aromatics decreases with ageing in line with their higher OH reactivity. Similar aromatic behaviour in RWC exhaust was observed earlier in a chamber with spruce exhaust , where the molar share of ArHCs in the total aromatic content increased from 45 % and 32 % to 63 % and 54 % during ageing of 0.6 and 0.8 eqv. d, respectively, while the share of furanoic and phenolic species decreased. Furthermore, N-aromatics were not detectable here with ages exceeding 1 eqv. d, although they have been observed to form with shorter exposures . The N-aromatics produced by the first stages of ageing may have partitioned to the particulate phase but are also degraded by subsequent reactions with OH. Simultaneously, the share of aliphatic nitrogen compounds (CHN and CHNO) to the total concentration was not substantially affected by photochemical exposure.
Carbonyls were divided into primary and secondary subgroups based on their behaviour during ageing. The primary carbonyl group, consisting mainly of acetaldehyde and to a smaller extent of compounds such as acrolein and butadiene, was prevalent in the fresh exhaust, but their share of identified compounds decreased from 13 %-27 % in unaged exhaust to 3 %-12 % in the highest exposures. This is the result of the degradation of these compounds and the introduction of high amounts of carbonyls in the secondary carbonyl group. The ratio between the two carbonyl groups increased linearly with age. The secondary carbonyl group was dominated by acetic acid, which was the most prevalent compound after extensive ageing in all experiments and covered over 30 % of the total measured OGC concentration from the highly aged S-5 % and B-17 % exhaust. The mainly small acidic compounds in the secondary carbonyl group were formed from the photochemical reactions of OGCs and from particulate OA, which is a consistent source of oxygenated OGCs such as acetic acid, formic acid, acetaldehyde, and acetone (Malecha and Nizkorodov, 2016). Of these, acetaldehyde is classified as a primary carbonyl because it reacts with OH 2 orders of magnitude faster than the other considered small carbonyls (Atkinson et al., 2006), and thus its concentration remained at a similar level throughout the ageing process. The photochemical ageing process increased the mass of particulate organic carbon measured with thermal-optical carbon analysis by factors of 1.3 and 3.9 for dry and moist spruce combustion in the masonry heater, respectively, and by factors of 2.0 and 3.0 for beech and moist spruce combustion in a chimney stove, respectively. This agrees well with the previously observed SP-HR-ToF-AMS-based OA enhancement factors (1.6-5.3) for RWC (Bertrand et al., 2017;Bruns et al., 2015b;Grieshop et al., 2009;Heringa et al., 2011;Tiitta et al., 2016) and with the thermaloptical analysis-based organic carbon enhancements previously measured for the same chimney stove (1.3-1.4 after 1.7-2.5 eqv d) (Miersch et al., 2019). The organic carbon concentrations after the PEAR OFR correlated well with the OA measured with SP-HR-ToF-AMS (R 2 = 0.85, Fig. S6). Ageing led to a linear increase in the SP-HR-ToF-AMSderived ratio of organic matter to organic carbon (OM : OC, Fig. S13), which rose from the initial average ratio of 1.8-2.2 to 2.7-3.0 during extended ageing. Similarly, the OS C of the OA increased as a function of the photochemical age throughout the tested exposure ranges (Fig. S13), indicating the existence of continuing reactions of the particulate phase after the rapid consumption of the majority of the primary gaseous SOA precursors. As expected, the organic carbon mass enhancement did not increase with continuous ageing because the major SOA precursors were already consumed by relatively short OH exposures. Continuing photochemical exposure may instead reduce the amount of particulate organic carbon , which acts as a source for volatile acidic compounds during photochemical ageing (Malecha and Nizkorodov, 2016;Paulot et al., 2011). See Table S8 for experiment-wise average concentrations and chemical compositions of the particulate matter measured by SP-HR-ToF-AMS.
The oxidation states measured after the PEAR OFR without oxidative ageing were highest for the chimney stove OA, with average unaged OS C 's of 0.22 and 0.41 for beech and spruce combustion, respectively, while dry and moist spruce in the masonry heater had average unaged OS C 's of −0.18 and 0.15, respectively. This indicates the existence of different combustion conditions in the studied appliances, with the masonry heater having lower emissions of highly oxygenated compounds and a higher share of unsaturated hydrocarbons compared to that of the chimney stove emissions. However, as a result of oxidative ageing, OS C surpassed 1.5 after 5 eqv. d regardless of the type of experiment, which exceeds the values typically observed in ambient aerosol (Kroll et al., 2011;Ng et al., 2011). This is likely because highly aged aerosol in ambient air is mixed continuously with fresh and less oxidized sources. However, the chemical evolution of OA in the PEAR OFR followed a similar trend to that ob- Figure 6. Van Krevelen diagram of the particulate organic aerosol measured by SP-HR-ToF-AMS. Numbers indicate batches combusted in the masonry heater (1-3: dry spruce; 4-5: moist spruce), whereas letters refer to chimney stove batches (A-E: beech; F-G: moist spruce) of this study. Ozonolysis experiments are marked with black and were not considered in slope calculations. Results from ageing of RWC exhaust with a PAM OFR by Bruns et al. (2015a) and Pieber et al. (2018), as well as behaviour of RWC OA aged in a chamber  are marked with scatter. served for typical ambient semivolatile oxygenated aerosol, with a van Krevelen slope of −0.5, which can be interpreted as a result of simultaneous fragmentation and acidgroup addition (Ng et al., 2011). Furthermore, the changes in OS C 's were similar to those previously observed for ambient aerosol aged extensively in a potential aerosol mass (PAM; Kang et al., 2007) OFR (OS C of up to 2; Ortega et al., 2016). The OS C 's of RWC exhaust aged in a PAM OFR have, however, been notably lower than observed here ( Fig. 6; Bruns et al., 2015a;Pieber et al., 2018).
A higher initial oxidation state in the chimney stove exhaust led to shallower slopes (−0.17 and −0.34) in the van Krevelen diagram than those in the masonry heater exhaust (−0.46 and −0.49, Fig. 6). The slope for the beech combustion aerosol was higher than for spruce, which may result from differences in fuel composition as has been observed also previously for biomass burning OA (BBOA) (Ortega et al., 2013). Previously, the ageing of spruce combustion exhaust from a masonry heater in a chamber has produced similar but slightly steeper van Krevelen slopes of −0.64 to −0.67 . The steepness may be affected by the shorter ageing (< 1 eqv. d) and consequent lower final OS C (maximum of +0.14) in the chamber, as the O : C ratio has been noted to level off at higher oxidation states (Ng et al., 2011). Slopes are also positively affected by fragmenta- tion, which may be enhanced in the PEAR OFR because of the more intensive UV radiation than in smog chambers or in ambient conditions. The van Krevelen slopes agreed well with those previously measured from single precursors aged in the PAM OFR, such as −0.48 for toluene or −0.46 for xylene (Lambe et al., 2011). Similar slopes have been observed also for ageing of OA from pellet combustion (−0.47 to −0.44;Czech et al. 2017) or open biomass combustion (−0.5;Ortega et al. 2013) with a PAM OFR.
The evolution of particulate organic aerosol was assessed also by the IDTD-GC-ToF-MS analysis of filter samples. When comparing the concentrations in the exhaust after the PEAR OFR, the concentrations of compounds with high oxidation states and low carbon number (n C ) increased during photochemical ageing (upper-right corner of Fig. 7). The locations of the measured organic compounds in the OS C : n C space are shown in Fig. S14, and their dilution-corrected concentrations downstream of the PEAR OFR in Table S11. In the OS C : n C space, the compounds which exhibited a major increase during photochemical ageing were located in or above the location of the low-volatility oxidized organic aerosol (LV-OOA) classified by Kroll et al. (2011). These compounds are products of the multistep oxidation process including both fragmentation and oxidative reactions.
OA formation is tied to the availability of organic precursors, and thus the formation of SOA was highest at the ignition phase of each batch (Fig. S15). In the masonry heater, the enhancement in particulate OA mass increased considerably with the introduction of moist logs, simultaneously with the increase in gaseous organics in the fresh emissions compared to the preceding dry spruce combustion. In contrast, the low organic emission by dry spruce combustion was reflected as a lower SOA formation. Another aspect related to the primary exhaust was the extent of OH exposure in the PEAR OFR, which was directly connected to the sample concentrations, as discussed in Sect. 3.3. Namely, the photochemical ageing was lower during periods of high emission, leading to lower oxidation states for the OA emitted during ignition.

PMF analysis of particulate OA composition
PMF analysis applied to the exhaust produced a four-factor solution for the OA covering 98 % of the spectra. The spectra of the factors are shown in Fig. S16. Two of the factors were associated in particular with the primary OA from biomass combustion: pyrolysis-BBOA, formed especially during ignition, and flaming-BBOA, emitted pronouncedly during the flaming phase. The other two factors, semivolatile oxygenated OA (SV-OOA) and low-volatility oxygenated OA (LV-OOA), represent oxygenated organics with notably  higher OS C 's than those of the primary OA factors (Table 3).
Flaming-BBOA comprised 76 % and 55 % of the unaged OA from dry and moist spruce combustion in the masonry heater, respectively, but only 27 % and 23 % in the beech and moist spruce combustion OA from the chimney stove, respectively (Fig. 8), indicating a less-oxidizing higher-temperature flaming phase in the modern masonry heater. Flaming-BBOA is strongly related to the C 4 H + 9 ion (main peak of m/z 57; R 2 = 0.90, Fig. S17), which is used as a tracer for hydrocarbonlike compounds (Aiken et al., 2009). However, in contrast to the typical hydrocarbon-like OA factor characterized by a relatively low OS C (−1.7 to −1.6; Kroll et al., 2011), the flaming-BBOA contains more oxygen-containing functional groups and is similar to the primary biomass burning OA factor measured in the RWC exhaust in a chamber θ = 13.6 • ; R 2 = 0.95). Pyrolysis-BBOA, on the other hand, consisted of ions typical to the low-temperature pyrolysis products of wood combustion and its time series correlated well with the PACs (R 2 = 0.86, Fig. S17). Of the more-oxygenated factors, the LV-OOA was dominated by the CO + and CO + 2 ions and thus represented highly oxidized OA. The LV-OOA spectra corresponded well with the OH-induced SOA factor identified from the RWC exhaust aged in a smog chamber θ = 8.4 • ; R 2 = 0.98) and was comparable with that of ambient LV-OOA (Aiken et al., 2009) and the LV-OOA spectra of unaged wood combustion exhaust . LV-OOA was also present in the unaged exhaust of this study, excluding the dry spruce combustion in the masonry heater, which produced the least-oxidized primary exhaust. The SV-OOA, on the other hand, was related to the C 2 H 3 O + ion (R 2 = 0.78), which is indicative of carbonyl formation in the particulate OA (Ng et al., 2010). The SV-OOA was also comparable (θ = 11.5 • ; R 2 = 0.96) to a factor of SOA generated in a chamber in the previous work , where this factor was interpreted as ozonolysis-generated organic OA. Interestingly, SV-OOA was formed only during photochemical ageing and increased in line with higher exposures despite semivolatile compounds being products of the initial stages of OA oxidation. This further demonstrates the longcontinuing functionalization of OA during ageing alongside with the fragmentation processes.

Polycyclic aromatic compounds
PACs were overall more prominent in the primary exhaust of dry spruce combustion (3.2 % of total OA, Fig. 9) than from combustion of moist spruce (1.7 % or 2.0 %) or beech (1.7 %). As expected, the total PAC concentrations decreased because of ageing and contributed less than 0.5 % to OA for all cases after 3 eqv. d of ageing. Furthermore, ageing transformed the composition of the PACs assessed with the SP-HR-ToF-AMS HR-PAH analysis (P-MIP analysis, Sect. S5.2; Herring et al., 2015). While UnSubPAHs formed the most prominent PAC group of all the combustion experiments, ageing decreased their share from 60 % in unaged exhaust to 40 %-50 % after 3 eqv. d of ageing. Similarly to the PACs measured with SP-HR-ToF-AMS, ageing decreased the UnSubPAHs analysed with IDTD-GC-ToF-MS by 83 %-85 % in the dry spruce combustion in the masonry heater and by 90 %-98 % in the other situations. Of the most prominent UnSubPAHs, anthracene and fluoranthene degraded within the first eqv. d (Fig. 10). Also benzo[a]pyrene, which is used as the marker for total ambient PAHs (European Council, 2004), degraded by a factor of 5 as a result of photochemical ageing (Table S11).
The change in the PACs is likely to alter the potential health effects of the exhaust; although the total PAC concentration decayed, the simultaneous formation of oxy-and nitro-PAC derivatives known to be detrimental to health was observed. These substituted PACs have lower vapour pres- Figure 9. Relationship of photochemical ageing to the batchwise average ratios of (a) total HR-PAH concentration to OA concentration and of (b-f) HR-PAH subgroups to the total PAC concentration measured by the SP-HR-ToF-AMS. sures compared to those of parent PAHs and thus are more likely to condense on the particles (Shen et al., 2012). The share of oxygenated PACs to the total HR-PAH concentration measured by SP-HR-ToF-AMS increased from 15 %-19 % in unaged exhaust to 25 %-38 % in aged exhaust. The concentrations of both the SP-HR-ToF-AMS HR-OPAH and IDTD-GC-ToF-MS-derived OH-PAH also correlated with the SV-OOA PMF factor (Pearson r = 0.70 for OPAH and r = 0.88 for OH-PAH; see Table S12) pointing towards their continuous formation during ageing. Of the compounds measured by IDTD-GC-ToF-MS, the most notable increase was observed for naphthaldehydic acid, with high concentrations in aged aerosol (up to 100 µg m −3 in dilution-corrected flue gas). Interestingly, its concentration was highest at approximately 1 eqv. d, after which it decreased. Naphthaldehydic acid and other oxygenated PACs have also been previously found to form during photochemical ageing of RWC ex-haust (Bruns et al., 2015b;Miersch et al., 2019); however, we found that the photochemically enhanced naphthaldehydic acid concentration degraded after continuous ageing, although remaining considerably higher than in the unaged emissions.
In addition, the share of nitrogen-substituted PACs, including both NPAH and APAH, increased from a combined share of 5 % of HR-PAH tot measured by SP-HR-ToF-AMS in fresh exhaust to a maximum of 9 % in the aged exhaust. In general, particulate nitrogen-substituted PACs are formed in the atmosphere through the oxidation of gaseous PACs or via heterogenous reactions from UnSubPAHs and are also simultaneously degraded by photochemical reactions (Keyte et al., 2013). As UnSubPAHs in the present study were largely consumed after 3 eqv. d, the higher photochemical exposure times consequently led to situation where the precursors for 6372 A. Hartikainen et al.: Photochemical transformation of residential wood combustion emissions nitrogen-substituted PACs are not available, and therefore no further formation took place.

Organic acids
Photochemical ageing led to a considerable increase in particulate organic acids in exhaust aerosol analysed with IDTD-GC-ToF-MS. The amount of small multifunctional acids such as malic acid (C 4 H 4 O 4 ) and tartaric acid (C 4 H 6 O 6 ) increased by factors of up to over 200 (Fig. 10, Table S11). Increases in the amounts of also dicarboxylic acids such as succinic (C 4 H 6 O 4 ) and glutaric acid (C 5 H 8 O 4 ) were evident, although to a lesser factor. The comparison of IDTD-GC-ToF-MS measurements with the SP-HR-ToF-AMS measurements revealed that the concentrations of organic acids were connected to both LV-and SV-OOA factors (r > 0.7, Table S12), which also strongly increased during the ageing process. This was expected based on the association of SV-OOA with the C 2 H 3 O + ion indicating carbonyl formation.
In the n C -OS C space (Fig. 7), these compounds dominated the upper-right corner with the highest oxidation states and largest enhancement ratios and are close to the LV-OOA region specified for atmospheric OA by Kroll et al. (2011). The measured organic acids were mostly intermediately volatile products of the continuous fragmentation process. Thus, they may partition also to the gaseous phase where they are expected to participate in secondary gas-phase OH reactions. Overall, the concentrations of particulate organic acids were highest at approximately 3 eqv. d, after which they were degraded by further oxidation reactions.

Nitrophenols
An increase in nitrophenols in the particulate phase was evident during photochemical ageing, which increased the concentrations of 4-nitrophenol (4-NP) and 4-nitrocatechol (4-NC) by respective factors of 2-30 and 30-3000 compared to non-aged exhaust (Fig. 10, Table S11) as measured by IDTD-GC-ToF-MS. The highest attained 4-NC concentration corresponded to 2 % of total SP-HR-ToF-AMS-based OA concentration, which is a notable fraction for a single compound in an ageing aerosol. These secondary nitrophenols were products of OH + phenolic-compound reactions and may have originated from both gas-phase and heterogenous reactions (Harrison et al., 2005). Also the ozonolysis experiments led to an extensive formation of 4-NP, whereas the amount of 4-NC did not increase. This discrepancy might be the result of absence of the photochemical production of catechol, which is the precursor for 4-NC formation (Finewax et al., 2018).
Nitrophenol concentrations were highest at relatively low OH exposures (1-2 eqv. d) and decreased with increased ageing. Similar trends with OH exposure were seen in the NPAH and APAH concentrations measured with SP-HR-ToF-AMS. Nitrophenols are reactive towards OH and photolysis in both gaseous and particulate phases, and their amount decreased from both phases after the first equivalent days of ageing. Thus, the measured nitrophenols can be considered to decay efficiently under high-OH exposures and, therefore, are not ideal biomass combustion markers for long-range transported smoke. On the hand, their concentrations still remained higher in the highly aged exhaust than in the fresh exhaust. Furthermore, in the particulate phase, they likely contributed to the formation of organic acids which are formed during the continuous photo-oxidation of the nitrophenols (Hems and Abbatt, 2018).

Conclusions
The photochemical ageing of dynamically changing OA emitted from RWC was evaluated using the PEAR OFR to expose the exhaust to varying photochemical conditions for up to an equivalent of 1 week in the atmosphere. To evaluate typical northern and central European combustion emissions, two different appliances were used with regionally typical logwood fuels. While the primary concentrations of particulate OA were relatively similar for all of the assessed sources, the enhancement of the organic particulate carbon during ageing depended on the type of the fuel. In particular, the fuel moisture content affected the SOA production: dry fuel produced a lower organic mass and OA enhancement during ageing than moist wood, because of the significantly lower emission rates of the organic gases. However, a very low (5 %) logwood moisture content considerably increased the primary PM 1 emission because of the extensive soot formation. With current logwood combustion appliances, this presents a conflict in attempts to decrease emissions since the usually preferred dry logwood produces high BC emissions, while high moisture content of the fuel increases the amount of organic emissions, as shown here for spruce logwood. Similar phenomenon has, however, been recently noted also for hardwood fuels .
The particulate organic carbon mass in the RWC exhaust increased by a factor of 1.3-3.9 during photochemical ageing. Furthermore, photochemical ageing transformed the overall composition of the OAs. This was observed as a linear increase in the average carbon oxidation state of particulate OA throughout the investigated photochemical exposure range, while the ratio of organic carbon to total organic mass decreased. Photochemical ageing also caused multiform changes in the OA at the molecular level. Notably, small acidic reaction and fragmentation products became increasingly dominant in both particulate and gaseous phases with higher ageing. The concentrations of particulate nitrophenols were at their highest level after 1 eqv. d, after which they began to decay but remained higher than that in the primary exhaust. Similarly, nitrogen-containing aromatics were unobservable in the gaseous phase at longer exposures, although they have been observed to increase during exposures of less than 1 eqv. d . Ageing also enhanced the share of oxygen-and nitrogen-substituted polycyclic aromatic compounds in the PAC emissions. Of the oxygenated PACs, naphthaldehydic acid in particular increased considerably with concentrations peaking at approximately 1 eqv. d. However, PACs in total degraded almost completely after 3 eqv. d of ageing.
In general, several different oxidation mechanisms of organic aerosol are likely occurring simultaneously both in the atmosphere and in the OFR experiments of this study. However, based on the observed OA chemical compositions under a range of different OH exposures, different major transformation mechanisms for RWC exhaust under photochemical conditions can be roughly outlined. First, the results demonstrate that short OH exposures (∼ 1 eqv. d) are sufficient to functionalize the majority of gaseous SOA precursors of RWC exhaust and lead to their condensation into particulate phase. This mechanism dominates the overall OA transformation until the SOA precursors have been depleted. After this stage, the continuing OH exposure leads to further oxidation of particulate organic aerosol, which is likely explained by heterogeneous reactions between gasphase oxidants and particles. However, it is also possible that particulate-phase oxidation occurs via evaporation and homogeneous gas-phase oxidation followed by recondensation. It is important to note that partitioning of organic compounds is affected by dilution, which was lower in the PEAR OFR than in ambient air.
While several recent studies (Bertrand et al., 2017;Bruns et al., 2015b;Bruns et al., 2017;Grieshop et al., 2009;Hartikainen et al., 2018;Heringa et al., 2011;Pieber et al., 2018;Tiitta et al., 2016) have assessed the behaviour of RWC emissions in a relatively short timescale of less than 2 eqv. d, this study also highlights the importance of higher exposure levels towards chemical transformation of OA. Due to the potentially long atmospheric lifetimes of OA, long-term ageing is also important to consider in large-scale atmospheric models, which typically estimate SOA formation and characteristics based on short-term ageing experiments. The consideration of only the first stage of gas-phase functionalization and condensation may lead to underestimated oxygenation of the long-transported OA, while specific compound groups, such as nitrophenols or substituted-PACs, can be overestimated. In general, the potential health and climate effects of aerosols are to a large extent determined by their composition, which depends on their sources and the levels of atmospheric ageing. Thus, the characterization of aerosol emissions from different sources and their atmospheric transformation at different exposure levels would be crucial for assessment of the overall environmental effects of ambient air pollution.
Data availability. The data are available on request from the corresponding author.
Author contributions. AH, OS, PYP, MI, and HL designed the study. Measurements were performed by AH, PT, MI, PYP, MK, HL, HS, JT, and OS. AH, HC, MI, and OS made the assessments of the PEAR conditions. AH performed analyses of FTIR, SMPS, and PTR-ToF-MS data; PT and LH performed the AMS data analyses; JO performed IDTD-GC-ToF-MS analyses; and HK performed the thermal-optical analyses. OS, JT, RZ, and JJ supervised and acquired funding to the study. The paper was written by AH with contribution from all coauthors.