Physical and chemical properties of black carbon and organic matter from different combustion and photochemical sources using aerodynamic aerosol classification

The physical and chemical properties of black carbon (BC) and organic aerosols are important for predicting their radiative forcing in the atmosphere. During the Soot Aerodynamic Size Selection for Optical properties (SASSO) project and a EUROCHAMP-2020 transnational access project, different types of light-absorbing carbon were studied, including BC from catalytically stripped diesel exhaust, an inverted flame burner, a colloidal graphite standard (Aquadag) and controlled flaming wood combustion. Brown carbon (BrC) was also investigated in the form of organic aerosol emissions from wood burning (pyrolysis and smouldering) and from the nitration of secondary organic aerosol (SOA) proxies produced in a photochemical reaction chamber. Here we present insights into the physical and chemical properties of the aerosols, with optical properties presented in subsequent publications. The dynamic shape factor (χ ) of BC particles and material density (ρm) of organic aerosols was investigated by coupling a charging-free Aerodynamic Aerosol Classifier (AAC) with a Centrifugal Particle Mass Analyzer (CPMA) and a Scanning Mobility Particle Sizer (SMPS). The morphology of BC particles was captured by transmission electron microscopy (TEM). For BC particles from the diesel engine and flame burner emissions, the primary spherule sizes were similar, around 20 nm. With increasing particle size, BC particles adopted more collapsed/compacted morphologies for the former source but tended to show more aggregated morphologies for the latter source. For particles emitted from the combustion of dry wood samples, the χ of BC particles and the ρm of organic aerosols were observed in the ranges 1.8–2.17 and 1.22– 1.32 g cm−3, respectively. Similarly, for wet wood samples, the χ and ρm ranges were 1.2–1.85 and 1.44–1.60 g cm−3, respectively. Aerosol mass spectrometry measurements show no clear difference in mass spectra of the organic aerosols in individual burn phases (pyrolysis or smouldering phase) with the moisture content of the wood samples. This suggests that the effect moisture has on the organic chemical profile of wood burning emissions is through changing the durations of the different phases of the burn cycle, not through the chemical modification of the individual phases. In this study, the incandescence signal of a Single Particle Soot Photometer (SP2) was calibrated with three differPublished by Copernicus Publications on behalf of the European Geosciences Union. 16162 D. Hu et al.: Physical and chemical properties of black carbon and organic matter from different sources ent types of BC particles and compared with that from an Aquadag standard that is commonly used to calibrate SP2 incandescence to a BC mass. A correction factor is defined as the ratio of the incandescence signal from an alternative BC source to that from the Aquadag standard and took values of 0.821± 0.002 (or 0.794± 0.005), 0.879± 0.003 and 0.843± 0.028 to 0.913± 0.009 for the BC particles emitted from the diesel engine running under hot (or cold idle) conditions, the flame burner and wood combustion, respectively. These correction factors account for differences in instrument response to BC from different sources compared to the standardised Aquadag calibration and are more appropriate than the common value of 0.75 recommended by Laborde et al. (2012b) when deriving the mass concentration of BC emitted from diesel engines. Quantifying the correction factor for many types of BC particles found commonly in the atmosphere may enable better constraints to be placed on this factor depending on the BC source being sampled and thus improve the accuracy of future SP2 measurements of BC mass concentrations.

Although BC and BrC are very important for climate, they are poorly represented in 80 atmospheric models (Zuidema et al., 2016). This is in part due to the complex microphysical properties of BC and the lack of accurate refractive index (RI) descriptions for both BC and BrC (Liu et al., 2020b). Fresh soot particles often exist in the form of aggregates composed of primary spherules with an irregular and highly fractal geometry (Xiong and Friedlander, 2001;Wentzel et al., 2003). The morphology of these aggregates change markedly during the 85 atmospheric aging process, influencing the corresponding particle size and optical properties (Zeng et al., 2019;Zhang et al., 2008). For example, after condensation of gaseous species such as sulfuric acid or water (under high relative humidity (RH) environments) on soot particles, or coagulation with the pre-existing particles, soot particles can experience restructuring and the shape of the soot particles becomes more similar to a spherical particle 90 (Zhang et al., 2008). The morphology of BC particles can be measured directly by using Scanning Electron Microscopy (SEM) or Transmission Electron Microscopy (TEM) (Fu et al., 2006;Chen et al., 2018;Ellis et al., 2016). However, the SEM/TEM approach only provides particle shape information in two dimensions and do not provide real time characterisation.
Alternatively, the particle morphology can be determined by measuring its size and mass with 95 different techniques (Chen et al., 2018;DeCarlo et al., 2004). A conventional approach is to classify particles (generally using a differential mobility analyser, DMA, to select monodisperse particles on their mobility size) and then measure particle mass using a particle mass analyser (Zhang et al., 2008;Park et al., 2003;Park et al., 2004a;Park et al., 2004b;Chen et al., 2018) (Wu et al., 2019). From the resulting information about particle mass for 100 different particle mobility sizes, the dynamic shape factor (χ, defined as the ratio of the drag force on the particle divided by the drag force on the particle's volume equivalent sphere) and fractal dimensions (D f ) can be retrieved (DeCarlo et al., 2004). Mobility-mass fractal dimension (D fm ) has been reported over a wide range of 2.2-2.8 for diesel exhaust particles (Park et al., 2004b). D fm has been reported as higher than the D f, -defined as the scaling 105 exponents between the radius of gyration of an aggregate and the radius of primary spherules composing the aggregate -but the two are not always directly equivalent, particularly in the transition regime (Sorensen, 2011).
In recent decades, the RI of BC or BrC was derived by measuring the optical properties for particles of controlled size, with studies commonly utilising the DMA to classify charged 110 aerosols on their electrical mobility diameter (Cotterell et al., 2020). However, the DMA approach suffers from transmitting a highly polydisperse and multi-modal distribution in terms of physical size. Since this technique relies on the selection of particles according to their electrical mobility from an aerosol ensemble with a distribution of charges, transmitting particles with a single electrical mobility can be achieved for multiple combinations of 115 particle charge and size. Thus, in addition to singly-charged particles of the desired size, larger particles with charge states greater than unity are selected which can impact optical measurements significantly. Careful consideration of the impacts of multiply charged particles on subsequent RI derivations can go some way to reducing uncertainty in the resultant RI, but this nevertheless remains a significant contributor to uncertainty (Cotterell et 120 al., 2020;Zarzana et al., 2014;Miles et al., 2011). Thus, the classification of particles without relying on electrical charge should reduce the uncertainty in refractive index retrievals from measured aerosol optical properties. Important additional considerations in the retrieval of refractive indices from optical spectroscopy data are the aerosol morphology (described 6 above) and mixing state. The mixing state can be probed using the Single Particle Soot 125 Photometer (SP2), which can measure the refractory BC (rBC) mass content and optical size of individual particles. However the SP2 needs an empirical calibration to retrieve the rBC mass from the incandescence signal (Laborde et al., 2012a). The conventional method to calibrate the incandescence channel of SP2 is using size selected Aquadag standards (Acheson Inc. USA) and then correcting to a calibration representative of ambient rBC by a 130 constant factor of 0.75 (Laborde et al., 2012b). However, few experiments since have independently verified this across various soot types.
To address the issues mentioned above, the Soot Aerodynamic Size Selection for Optical properties (SASSO) project utilised the Aerodynamic Aerosol Classifier (AAC) to classify particles according to aerodynamic diameter for size and mass distribution measurements and 135 optical evaluation (Tavakoli et al., 2014). Specifically, SASSO has used the AAC size selection of emissions from wood burning, diesel combustion and secondary organic aerosol (SOA) formation, prior to optical measurements using cavity ring-down and photoacoustic spectroscopy with the EXtinction, SCattering and Absorption of Light for AirBorne Aerosol Research (EXSCALABAR) instrumentation, custom-built by the Met Office (Cotterell et al., 140 2020;Cotterell et al., 2019a;Davies et al., 2018).
The purpose of this paper is to determine the physical and chemical properties of BC and organic aerosols from different combustion sources representing the combustion of gas (methane), liquid (diesel) and solid (wood), and to examine the variation in calibration SP2 constants for BC particles generated by these combustion sources to enable accurate 145 characterization of BC mass concentrations and mixing state in future studies. The key objectives of this work are:

7
(1) Develop reliable and repeatable methods of generating isolated BC and BrC using controlled combustion sources and a smog chamber (2) Derive the dynamic shape factor of BC particles and material density for organic 150 aerosols (3) Determine mass spectral profiles of organic aerosols produced from wood combustion from an aerosol mass spectrometer (AMS) (4) Explore the restructuring of BC particles in response to the controlled coating and humidification of aerosol samples 155 (5) Investigate the response of the SP2 incandescence signal to different types of the BC particles.
The optical data presented in this paper is limited to the qualitative characterisation of the materials under investigation for the sake of achieving objective (1). The quantitative parameterisation of the optical properties of the particles and associated development and 160 application of ambient optical property models will be the subject of future publications.

Experimental setup and methods
Figure 1 shows our instrument configuration during SASSO. We used the AAC to classify aerosols based on their aerodynamic size prior to characterizing the particle size distribution, chemical composition, aerosol mixing state and optical properties for the AAC-selected 165 aerosols. The Centrifugal Particle Mass Analyzer (CPMA) and Scanning Mobility Particle Sizer (SMPS) sampled aerosols downstream of the AAC to measure the mass and mobility size distributions of the AAC-selected particles. A Nafion humidifier and a custom-designed thermal denuder (TD) were used for the BC restructuring test.

8
The right panel in Fig. 1 shows the various instrument configurations used in this study to 170 target different measurements: Instrument set (1) was used for measurements of aerosol optical properties (and thereby enable the retrieval of refractive index for BC/BrC aerosols, subject of a future publication) and organic chemical measurement for wood combustion (setup (a)) and chamber (setup (c)) experiments. Instrument set (2) was used for dynamic shape factor measurements for BC and material density measurement for organic aerosols 175 (setup (a) and (b)). Instrument set (3) was used for SP2 incandescence signal calibration (setup (b)). Instruments set (4) was used for BC restructuring experiments (setup (c)).

Aerodynamic Aerosol Classifier (AAC)
The AAC (Cambustion Ltd, Cambridge, UK) is used to select aerosols within a narrow range 180 of aerodynamic diameters and does not suffer from the issue of multiple charges that affects selection using instruments such as the CPMA and DMA. The AAC uses a centrifugal force and sheath flow between two concentric rotating cylinders to produce an aerosol classified by aerodynamic diameter. The detailed information regarding the principle of AAC can be found in Tavakoli et al. (2014).

EXtinction, SCattering and Absorption of Light for AirBorne Aerosol Research (EXSCALABAR)
The EXSCALABAR instrument used in this work was developed by the Met Office (Exeter, UK), which can be operated in both the laboratory and from the UK atmospheric research aircraft (FAAM BAe-146). The operating principle of EXSCALABAR has been described in 190 detail in previous papers (Davies et al., 2018;Cotterell et al., 2020;Cotterell et al., 2019a;Cotterell et al., 2019b). In brief, the instrument uses cavity ring-down spectroscopy 9 (CRDS) to measure the dry aerosol extinction at 405 and 658 nm wavelengths and photoacoustic spectroscopy (PAS) to measure the dry aerosol absorption coefficient at the wavelengths of 405, 514 and 658 nm wavelengths. For deployment during SASSO, all cells 195 shared common sample conditioning; the aerosol sample relative humidity was reduced to <10% by passing through a Nafion drier (Perma Pure LLC) and Ozone and NOx were removed by an activated charcoal "honeycomb" scrubber (custom-built in-house). The sample passed through an impactor (Brechtel Manufacturing Inc., custom built for an 8 L min -1 volumetric flow rate) with aerodynamic cut-off (D50) of 1.3 µm before being drawn using Ozone before and after each set of experiments as well as at least every working week during the experimental periods, as described by Davies et al. (2018) and discussed further by Cotterell et al. (2019a).

225
The refractory black carbon (rBC) mass concentration was measured by a SP2 (Droplet Measurement Technologies, Colorado, USA). The SP2 uses the laser-induced incandescence to measure the rBC mass and optical size of individual BC particles with an intra cavity Nd:YAG laser operating at 1064 nm. The particle size can be determined by detecting the laser signal scattered by particles, with the scattering intensity maximum related to the optical 230 particle diameter through a calibration using polystyrene latex spheres. The optical size of BC-containing particles is determined by matching the measured scattering signal with calculations from light scattering calculations assuming a core-shell structure (core-shell Mie theory) (Moteki and Kondo, 2007). Because the scattering signal of the absorbing particle will be distorted during its transit through the laser beam due to the mass loss by laser heating, 235 the leading edge scattering signal before the onset of volatilisation is extrapolated to reconstruct the scattering signal of the absorbing particles (Gao et al., 2007;Liu et al., 2014).
This calibration was performed at the start of the measurement campaign. For those particles which contain absorbing materials such as the refractory BC, they will absorb 1064 nm light and then heat up and emit visible thermal radiation (incandesce). This incandescence signal is 240 directly proportional to the mass of rBC as determined by a calibration (Liu et al., 2010) with generated BC aerosols of known or independently-measured mass. In the atmosphere, besides BC, other materials (e.g., some metals and minerals) can incandesce at 1064 nm as well.
However, as boiling point temperatures of these materials are rather different to that of black carbon, it is easy to distinguish them in measurements made using the SP2 equipped with an 245 additional narrowband incandescence detector, such as the one used here .
Recently, Sedlacek et al. (2018) demonstrated that charring of light-absorbing organic particles at 1064 nm can produce the refractory black carbon and then overestimate the rBC mass concentration, however this does not affect pure BC particles and we saw no evidence for an incandescence signal associated with the pure organic particles measured in this study. 250 We investigated the effectiveness of this calibration from different types of BC particles in this study and report the outcomes of these investigations in Sect. 3.3.

High-Resolution Aerosol Mass Spectrometer (HR-AMS)
Non-refractory aerosol chemical compositions including sulfate, nitrate, ammonium, chloride, and organics were measured by a HR-AMS in real time. The HR-AMS was operated in fast-255 mode to capture the fast transition of the combustion phase during the wood combustion experiment (Kimmel et al., 2011). The instrument operation and data analysis of HR-AMS have been described in detail elsewhere (Alfarra et al., 2006;Allan et al., 2003Allan et al., , 2004DeCarlo et al., 2006). The ionization efficiency of the AMS was calibrated using monodisperse ammonium nitrate particles according to the method described by Jayne et al. 260 (2000) at various times during the experimental periods.

Centrifugal Particle Mass Analyzer (CPMA)
The CPMA (Cambustion Ltd., Cambridge, UK) uses opposing electrical and centrifugal fields to classify particles according to their mass-to-charge ratio. The ability to vary the 12 electrical field and rotation speed enables particle selection based on their mass. The 265 principles and operation of the CPMA have been described elsewhere (Olfert and Collings, 2005;Olfert et al., 2006). Combining the CPMA with a CPC (Model 3776, TSI, USA) and scanning across the mass range of interest provides the bulk aerosol mass size distribution.
As the CPMA uses an electrical classification method to select particles, it also suffers from an issue of multiple charging similar to a DMA, and it has problems associated with a 270 fraction of the uncharged particles that are also transmitted, particularly at the lower rotation speeds. In this study, an electrical ioniser (MSP Corp., USA) was used for wood combustion experiments and a 90 Sr radioactive ioniser was used for chamber experiments to neutralize particles before they were sampled by the CPMA. 275 Aerosol size distributions in the diameter range from 14.9 to 673.2 nm were measured by a commercial SMPS (TSI, USA, employing a model 3082 classifier unit, 3081 DMA and 3786 CPC) with sheath and sample flow rates of 3 L min -1 and 0.3 L min -1 , respectively. The DMA was operated in scanning mode with a scan time of 60 s and a retrace time of 4 s. The particle size distribution was corrected for multiple charge effects and diffusion loss using the 280 standard inversion algorithm in the SMPS software (AIM version 10). Before the experiment, the SMPS was calibrated using NIST certified polystyrene latex spheres (PSLs, Thermo Fisher Inc.).

Transmission electron microscopy (TEM) sampling and analysis
To investigate morphological properties, BC particles were collected on carbon coated copper 285 grids using an electrostatic sampler (ESPnano, DASH Inc., USA) (Miller et al., 2010) for transmission electron microscopy (TEM) analysis. The grids were analyzed using a FEI Titan 13 operated at 300 kV. Several areas of the grid were imaged to provide representative images for each sample since both the grid film and particles contained carbon. 290 We explored the impacts of the partitioning of coating materials on the structure of BC. For this BC restructuring experiment, a thermal denuder was utilised to remove the coating materials on coated BC particles. The home-built TD consisted of a stainless steel tube in a temperature controlled furnace . with additional temperature sensors on the outside of the tube. It is necessary to ensure flow through the TD whenever it is heated, even when bypassed to allow measurements of the unheated sample. A constant 2-2.5 L min -1 flow of the sample air through it. A vacuum line 300 maintained 2.0-2.5 L min -1 through whichever of the bypass or TD line was not in use. The residence time of the air sample in the heating section was 216-270 s. The TD was given 30 mins for its temperature to stabilize before sampling. The temperature of the TD was calibrated by measuring the temperature at the axial centre of the denuder's heating zone. In this study, the purpose of the TD is to remove as much of the organic coatings from the 305 combustion-generated particles as possible, rather than probing the volatility properties of the coating. Therefore, all heating zones of the TD were set to their maximum temperature of 180 °C. This upper temperature is lower than that achieved by other commercial TD units and minimises the risks of charring.

14
The multi-tube Nafion humidifier in the Manchester home-made Hygroscopicity Tandem Differential Mobility Analyser (HTDMA) system was used for the BC restructuring experiment. The principle and configuration of the HTDMA system including the Nafion humidifier flow system can be found in Good et al. (2010). The RH in the Nafion humidifier was controlled by adjusting the relative mixing ratio of dry air from a compressed air source 315 with humidified air generated by bubbling compressed air through a water-filled glass bulb.

The Manchester Aerosol Chamber
The Manchester Aerosol Chamber comprises an 18 m 3 collapsible Fluorinated ethylene propylene Teflon bag (3m (H) × 3m (L) × 2m (W)). It was used by Alfarra et al. (2012) to investigate the effect of photochemical aging and initial precursor concentration on the 320 composition and hygroscopic properties of secondary organic aerosol. The chamber was run as a batch reactor in which the composition of the gaseous precursors, oxidising environment, primary emissions or seed particles, relative humidity and temperature were controlled. Air was supplied to the chamber by a blower at a flow of 3 m 3 min -1 . The air was dried and filtered for gaseous impurities and particles using a series of Purafil (Purafil Inc., USA), 325 charcoal and HEPA filters (Donaldson Filtration, USA), prior to humidification with ultrapure deionised water. Halogen bulbs and two 6 kW Xenon arc lamps were mounted on the inside of the enclosure housing the bag, which was coated with a reflective space blanket (mylar) to maximise the irradiance in the bag and to ensure even illumination. The Xenon arc lamps were mounted on two opposite sides of the enclosure at different heights. The 330 combination of illumination was tuned and evaluated to mimic the atmospheric actinic spectrum over the wavelength range 290 -800 nm, and had a maximum total actinic flux of 1.4 × 10 14 (photon s -1 m -2 nm -1 ) over the region 460 -600 nm. The calculated j(O 1 D) value during the reported experiments was 1.23 × 10 -5 s -1 (290 -340 nm) and j(NO 2 ) was 1.5 × 10 -3 15 s -1 (290 -422 nm). The relatively large volume of the chamber allows the dilute sample to be 335 held for several hours without significant aerosol removal from wall losses, allowing the study of particles introduced directly or formed within the chamber over a period of several hours. The mean number and mass wall loss rates of particles inside the chamber were estimated as 9.17 ± 1.3 and 8.16 ± 1.5 × 10 -5 s -1 , respectively . Further details on the chamber are given by Alfarra et al. (2012).

Diesel Engine
The engine used in this study was a Volkswagen 1.9L SDI light duty diesel engine (EURO 4 car equivalent), mounted on a test rig (CM12; Armfield Ltd., Hampshire, UK) and coupled to an eddy current dynamometer. This setup was used for investigating the light absorption properties of black carbon aerosol by Liu et al. (2017) previously. The engine throttle and the 345 load on the dynamometer were controlled with dedicated software. The engine exhaust was passed through an oxidising catalytic converter (a retrofitted diesel oxidation catalyst consisting of a mixture of platinum and rhodium) followed by a standard Volkswagen silencer, then a computer operated pneumatic valve connected to the chamber, with 2 inch ID stainless steel tubing between each component. The whole length of the exhaust line was 4 m.

350
The fuel used was standard UK low sulphur diesel, obtained from a local fuel station.
For each run, the engine was started and warmed up at 2000 rpm with 30% load. Once the engine reached a steady temperature (after ~ 10 mins), a controlled amount of the exhaust was injected into the Manchester Aerosol Chamber by switching the valve for a set time (referred this condition as "hot engine"). For the cold start runs (in which no throttle or load 355 was applied, denoted as "cold idle"), the exhaust sample was injected within a minute of the engine starting.

Inverted flame burner
A miniature inverted soot generator (MISG, Argonaut Scientific) operated with propane was used to generate black carbon aerosol. The design, operation, and performance of the MISG 360 has been previously reported Moallemi et al., 2019). The MISG was operated under two conditions: with the flow rate of propane and air at 0.0625 L min -1 and 10 L min -1 respectively (Condition 1), and 0.0625 L min -1 and 7.5 L min -1 respectively (Condition 2). Under these two conditions, the MISG produces black carbon with average Elemental Carbon/Total Carbon of 95%, determined via thermal-optical analysis.

Derivation of aerosol dynamic shape factor and material density
In this study, the SMPS and CPMA were arranged downstream of the AAC to measure the mobility diameter (D m ) and mass (M p ) of the particles within a narrow size range selected by the AAC based on their aerodynamic diameter (D a ). These parameters enabled the dynamic shape factor (χ) of the BC particles to be derived by using Eqn. (1) under the assumption that 370 the density ( ) of pure BC is 1.8 g cm -3 (DeCarlo et al. (2004)), in which C(D) is the Cunningham slip correction factor for a particle of diameter D using the parameters specified by Kim et al. (2005), and D ve is volume equivalent diameter which can be calculated from M p using Eqn. (2).
in which 0 is the reference density (1 g cm -3 ). This value D ve is then substituted into Eqn.
(1) 380 to calculate χ. We then iterate through a range of χ trial until χ and χ trial converge. The equations were solved iteratively using a standard Brent's method solver ('findroots' command, Igor Pro version 6.36, Wavemetrics).
In addition, the material density ( ) of organic aerosols can be calculated from Eqn. (4) and (5)  Exeter ( Fig. 1(a)). The main difference between the FPA and iCone approaches is that the there is a forced flow of air from beneath the sample in the FPA, whereas in the iCone the sample sits in ambient air conditions. Otherwise, both approaches rely on oxygen consumption calorimetry to measure the heat release rate from a burning object.

18
The aim of the wood combustion experiments was to produce different phases of combustion 400 that were as well separated as possible, rather than under natural burning conditions. This During the experiment, a wood sample (L×W×H (mm): 90×90×30) was first placed into a 420 custom-made stainless steel basket and then exposed to a radiant heat flux of 40 kW m -2 using the iCone Calorimeter. This radiant heating thermally decomposes the wood sample and generates volatile gases (pyrolysate) as the wood begins to pyrolyse. A continuously 19 operated spark igniter that was positioned in the released stream of pyrolysate acted as a source of ignition. Immediately following the initial placement of the wood sample under the 425 radiant heat source, pyrolysate was continuously released and the pyrolysate concentration immediately above the wood sample and in the vicinity of the spark igniter increased also.
The sample ignited once this pyrolysate reached sufficient concentrations and was well mixed with the surrounding air, at which point the igniter was switched off and removed from the air flow. Measurements were taken as soon as a sample was exposed to the radiant heat flux.

430
Each sample was allowed to flame for 5 minutes and the aerosols captured from this phase and then the flames were snuffed out manually leaving the fuel smouldering so that separate measurements of the smouldering phase could be taking alone for a further 5 minutes.
Throughout the experiment, the exhaust emissions from combustion were collected in a hood and entered into an exhaust duct. The concentration of oxygen, carbon dioxide and carbon 435 monoxide in the exhaust gas were measured by gas analysers integrated into the iCone system and used to calculate the heat release rate from the burning fuel. We installed an outshoot on this duct to allow part of the exhaust to be carried to a separate sampling system that allowed the measurement of the aerosol properties in the exhaust. Thus monitoring could occur in real time by the instruments (described above and shown in Fig. 1), after the exhaust 440 was diluted with the compressed air through a set of Dekati DI-1000 ejector diluters.

SP2 incandescence signal calibration and BC morphology investigation experiment
For the purposes of SP2 incandescence signal calibration and morphology investigation of BC particles, three types of the BC particles other than that emitted from wood combustion were investigated through chamber experiments ( Fig. 1(b)). For the BC particles generated by 445 the MISG: BC particles were introduced into the chamber through a Dekati DI-1000 ejector diluter. For the BC particles from the Aquadag standards: the Aquadag standards was first 20 generated by an aerosol atomizer (Topas ATM 226) and then passed through a diffusion drier (using silica desiccant) prior to entering the chamber. For the BC particles emitted from the engine: a hot engine running condition (2,000 r.p.m., 30% load and 10 min warm up) and a 450 cold idle condition were investigated. After the BC mass concentration had accumulated to adequate levels in the chamber, a catalytic stripper (operated at 350°C, Model CS10, Catalytic Instruments, Germany) was used to remove the coatings on the BC particles before they were sampled by the instruments.

455
For the purpose of testing models of BC and SOA mixing during the SASSO project, two procedures for generating SOA were used: one that yielded non-absorbing SOA and one that produced brown carbon. Experiments were conducted in the photochemical aerosol reaction chamber at the University of Manchester. During the experiment, NO 2 (10% v/v, with a balance gas of high purity N 2 (BOC, UK)) was injected directly into the bag from a custom-460 made gas cylinder via stainless steel tubing, and its concentration was measured using a chemi-luminescence gas analyzer (Model 42i, Thermo Scientific, MA, USA). After the NO x concentration in the bag reached the desired value, the precursor VOC was injected into a heated glass bulb (80°C) and flushed into the bag with high purity nitrogen. Subsequently, chamber lights were turned on to trigger SOA formation. The EXSCALABAR system was 465 used to determine the light absorption properties of the formed SOA particles in real time throughout the experiment.
In this study, 50 ppb NO x and 250 ppb α-pinene (Sigma-Aldrich) was used to form nonabsorbing SOA; this route to forming non-absorbing SOA is well characterized and has been reported previously (Nakayama et al., 2010). Here, we only focus on the methodology used 470 to produce the brown carbon SOA. Previous work has reported 'brown' SOA formation 21 through oxidation of aromatic precursors under high NO x conditions . In this study, we injected ~ 400 ppb cresol (Sigma-Aldrich) and ~ 15 ppb NO 2 (~ 60 ppb NO x in the chamber) into the chamber, and removed the UV filter from one of the Xenon arc lamps to increase photochemistry and accelerate and enhance the 'brown' SOA formation.

BC restructuring experiment
For the purpose of testing the influence of organic coatings and relative humidity on the structure of the BC particles, a series of experiments were conducted in this study. First, bare BC particles with minimal accompanying VOCs were injected into the chamber by using an oxidising catalytic converter, heated ejector dilutors, a catalytic stripper, Purafil, and 480 activated charcoal denuders between the exhaust line of the engine and the chamber inlet.
After the BC reached the desired concentration, BC injection was stopped. The particle size distribution of the dried bare BC particles was scanned by the AAC. In addition, the particle size distribution of the dried bare BC particles following 'humidity cycling', i.e. exposure to 90% RH for around 10 s and then dehydration to 10% RH, was measured by AAC for 485 comparison. Afterwards, we injected 50 ppb NO x and 250 ppb α-pinene into the chamber.
Photochemical reaction was initiated by turning the chamber lights on, and organic materials condense onto the BC particles. After the condensed organics equilibrate with the surrounding VOCs and the particles stabilised at a certain size (295 nm in aerodynamic size), the aerodynamic particle size distribution of the dried organic coated BC particles, and that of 490 the dried organic coated BC particles experienced the humidity cycling, was measured by AAC. Hereafter, a thermal denuder operated at 180 °C was added, and the organic coatings of the coated BC particles, that had either experienced the humidity cycling process or not, were removed by the TD, and then the size distribution of the BC core was measured by the AAC and SMPS. Immediately after ignition, the number (and mass) concentration of the OA decreased abruptly, while the rBC mass concentration increased sharply to 150 µg/m 3 (Fig. 2). We refer to this combustion phase following ignition as the flaming phase. The r OA at the onset of flaming combustion (the first 90s after ignition) is less than 0.05, indicating that most of the 515 emitted aerosols during this period are BC. As flaming combustion progressed, the extent and depth of charring increased and as the flammable gas flux slows the flame begins to diminish.

Evolution of BC and organics through wood combustion process
Char is a strong insulator and hence the char layer forms an obstacle to the conduction of heat into the lower uncharred layers of wood, reducing the production rate of the flammable gases 23 and decreasing the rate of heat release from combustion. As flaming subsides the number 520 (and mass) concentration of BC particles decreases (Fig. 2) and the r OA increased to around 0.2 after 300 s after the ignition.
After ~300 s (5 min) from the time of ignition the flames were snuffed manually to force the combustion from the flaming to the smouldering phase, thus transitioning the oxidation reaction from that on gas phase species to direct oxidation of the solid fuel. Fig. 2 shows that 525 the number (and mass) concentration of rBC in the smouldering phase decreased abruptly as flaming combustion ceased and the OA concentration increased sharply. The r OA was 1.0 during this period, indicating that almost all the emitted aerosols are composed of organic species only.

Brown carbon formation 530
As shown in Fig. 3, after the precursors (~ 400 ppb cresol and ~ 60 ppb NO x ) were introduced into the chamber and the lights were switched on, the SOA started to form and the particle size increased to ~ 200 nm after 40 minutes. Over this same time period, we measured increasing aerosol absorption and extinction coefficients, the single scattering albedo (SSA) was observed to increase sharply at the onset of the SOA formation and then 535 stabilized at around 0.81 (at 405 nm) eventually. Most of the evolution in SSA (i.e. increasing over time) is expected to be caused by the increase in particle size over a size range where SSA is very sensitive to particle size. We emphasise that we did not inject any primary ozone in this experiment, and the abrupt increase in apparent ozone concentration around the time of VOC injection is likely caused by the cresol (as an aromatic compound) 540 causing optical interference in the ozone instrument.

24
The most common way to calibrate the incandescence channel of SP2 is using monodisperse Aquadag standards (Acheson Inc. USA), which is then corrected for ambient rBC by applying a constant factor of 0.75 based on a previous laboratory comparison (Laborde et al.,545 2012b). As the incandescence signal is not independent for the different types of BC particles, using a constant factor of 0.75 may represent an uncertainty in the retrieved BC mass concentration. Thus, quantifying the correction factor for many types of BC particles found commonly in the atmosphere may enable better constraints to be placed on this factor depending on the BC source being sampled, thus improving the accuracy of future SP2 550 measurements of rBC mass concentrations. The broad range of BC aerosols generated in this study from different combustion sources served as an ideal platform to assess the variation in this correction factor to SP2-derived rBC mass.
In this study, the incandescence signal of the SP2 was measured for BC particles from catalytically stripped diesel engine exhaust emissions, an inverted flame burner, and 555 controlled flaming wood combustion, respectively, and compared with that measured from an Aquadag standard. The uncertainties here refer to precision of the fitted parameters reported by the Igor Pro fitting algorithm, based on analysis of residual data. As shown in Fig. 4, for the BC particles emitted from the diesel engine under hot engine and cold idle conditions ( Fig.   4(a)), the slopes of the incandescence signal with BC mass are 0.821 ± 0.002 and 0.794 ± 560 0.005 times of that measured from the Aquadag standard, respectively. Note that while some deviation from a perfect linear response is noted, this is small compared to the variation in slopes, so represents a minor source of uncertainty in comparison. These correction factors are 9.4% and 5.6% different with the common value of 0.75 (with the uncertainty less than 5%) recommended by Laborde et al. (2012b) when deriving the mass concentration of BC 565 emitted from diesel engines. For the BC particles generated from the flame burner ( Fig. 4(b)),  (Sedlacek et al., 2018), as shown in Fig. 2(c), the BC particles generated at the beginning of the flaming phase contained almost no organic species, with r OA values less than 0.05. Even if this OC were to be converted to EC with 100% 575 efficiency (which we consider to be highly unlikely), this would represent a very small error.
The differences in the correction factors derived in this study with the default value of 0.75 are 9.4% (5.6%), 17.2% and 12.4-21.7% for the BC particles emitted from engine with hot engine (or cold idle) condition, flame burner and wood combustion, respectively. We recommend that future studies utilizing the SP2 for rBC mass concentration measurements 580 use SP2 calibrations with the same type of BC as that to be studied.

Physical properties of black carbon and organic aerosols
The morphology of BC particles change markedly during atmospheric aging, with associated impacts on the particle size and optical properties (Zeng et al., 2019;Zhang et al., 2008;Teoh et al., 2019). In this study, the dynamic shape factor of BC particles and the material density 585 of organic aerosols were derived by measuring the mass and/or mobility diameter of AACselected particles. For biomass burning aerosols generated through our controlled combustion of dry wood samples, the AAC was set to pass aerosols with an aerodynamic diameter D a =200 nm and we sampled BC or organic aerosols produced from the flaming or smouldering phases, respectively. Number concentration measurements from the CPC were 590 recorded as the CPMA was scanned to determine the aerosol mass distribution. Figure 5(a) shows that there are considerable differences in the mass distributions for the produced BC (during the flaming phase) and organic aerosols (during the smouldering phase). The peak mass of BC particles produced from the wood samples of the Scots pine, Poplar and Giant Redwood are 10.86, 12.47 and 7.00 fg, respectively. Meanwhile, the corresponding peak 595 mass for the organic aerosols are 1.72, 1.74 and 1.90 fg, respectively. This difference is likely due to the difference in density and morphology for BC and organic aerosols. BC normally has a large density and a more irregular morphology than organic aerosols. After the peak masses for the AAC-selected aerosols were ascertained, the CPMA was set to these values to further select the AAC-selected aerosols (D a =200 nm) by their peak mass for a subsequent 600 mobility size measurement using our SMPS. Figure 5(b) shows that the corresponding peak mobility size of the BC and organic particles after selection according to aerodynamic diameter are 394, 405, 292 nm and 139, 139, 140 nm for Scots pine, Poplar and Giant Redwood, respectively. From these aforementioned peak mobility diameters and aerosol masses, the χ of the BC particles and the of the organics can be calculated from the 605 methods described in Sect. 2.2. We highlight the discussion in Sect. 2.2 that stated that  could be calculated using either the aerodynamic or volume-equivalent diameter; as organic species may evaporate during the AAC selection (due to the high sheath flow rate) and thus bias the D a , only the mass and mobility size distributions were used in these calculations. of emitted particles from their combustion, in addition to any influence on the pre-ignition pyrolysis phase. The densities of organic aerosols reported/used in previous studies vary significantly, from 0.6 to 1.4 g cm -3 (Turpin and Lim, 2001;Nakao et al., 2013;Li et al., 2016).
For the biomass burning aerosols, Zhai et al. (2017) reported the effective density ranged from 1.35 to 1.51 g cm -3 for the carbonaceous aerosol produced from the agricultural residue 630 burning. Our results would expand the available data on the variation densities for organic aerosols from the wood combustion.
For the BC particles from the catalytically stripped Aquadag standard, diesel engine exhaust, and inverted flame burner, particles were first selected by the AAC prior to mobility size distribution measurements using the SMPS. The χ of the BC particles were derived from the 635 controlled aerodynamic size and the measured peak mobility diameter. Fig. 7(a) shows that the particle mobility diameter increases with the selected aerodynamic size for all three types of BC particles, while Fig. 7(b) shows the corresponding variations in . These figures demonstrate clearly that the dynamic shape factor varies with both the aerosol aerodynamic 28 size and with the BC source. For the BC particles generated from the Aquadag standard, the χ 640 increased sharply from 1.315 ± 0.004 to 1.460 ± 0.003 as D a increased from 75 nm to 125 nm, but then decreased with further increases in D a to a value of 1.255 ± 0.005 at D a =500 nm.
This result implies that the formed BC particles are more irregular as D a increases initially, but particles adopt more collapsed morphologies as D a increases further. These trends in  likely arise from the coagulation of the primary BC particles, which drives the fractal-like 645 morphology for small BC particles and relative compact shape for large particles. For BC particles generated from the diesel engine, the χ exhibits a monotonically decreasing trend with particle size over the reduced size range probed for this BC source, with  decreasing from 1.329 ± 0.003 at D a =150 nm to 1.243 ± 0.002 at D a =275 nm. The TEM images ( Fig.   7(c)) show that the primary spherules generated directly from the diesel engine are around 20 650 nm. As the BC particles selected in the experiment are much larger than the primary spherules, the coagulation process drives the large particles towards more compacted morphologies than the small particles and this leads to the decreasing trend of  with particle size. For the BC particles generated from the flame burner, the χ increases with particle size for both operating conditions and a maximum value was reached of 2.731 ± 0.021 for D a = 655 150 nm. The same increasing trend was reported by Slowik et al. (2004) previously, who reported an increase in χ from ~ 1.3 to ~ 3.0 when the mobility size increased from 100 nm to 300 nm for the BC particles generated with lower propane/O 2 ratio compared to the measurements we report here. The TEM images (Fig. 7(d)) shows that, similar to the BC particles emitted from diesel engine, the primary spherules generated from the flame burner 660 are also around 20 nm. As our experiments on BC particles from the flame burner selected particles at small sizes (D a in the range ~ 50 -150 nm), the increase of the dynamic shape factor with D a is probably due to the coagulation of primary generated particles. In addition, 29 by comparing Fig. 7(c) and Fig. 7(d) (more images are shown in Fig. S1), it is clear that the large BC particles generated from the flame burner tend towards aggregate morphology, in 665 contrast to the compact morphology of the large BC particles emitted from the diesel engine ( Fig. 9(c)). This leads to the flame burner generated BC particles exhibiting larger  than those from the diesel engine.

BC restructuring
Recent work has shown that soot particles maintain their structure after coating with SOA 670 and subsequent thermal denuding (Bhandari et al., 2017), but will restructure in response to humidification (Leung et al., 2017). We performed a series of experiments to better understand this restructuring. As shown in Fig. 8, there is no clear change in the measured particle aerodynamic diameter for the bare BC after they experience the humidity cycling process, implying that the bare BC particles retain their structure even after they experience 675 very large perturbations in RH, this result is similar to that reported by Bhandari et al. (2019).
In addition, no size change of the organic coated BC particles was observed after they experienced the humidity cycling process either. Based on these results, we cannot form conclusions on whether the BC core was restructured during the humidity cycling process.
There are two possibilities: (1) The BC cores retained their structure throughout humidity 680 cycling process; or (2) As shown in Figure 8, the coatings on BC particles are very thick and therefore dominate the particle, rendering our size measurement approach insensitive to any changes in BC core size from restructuring. The size of the bare BC particles is around 68 nm, but after coating by SOA, the size of the coated BC particles reaches up to around 300 nm.
Due to the very large coating thicknesses of the coated BC particles, even if the BC cores 685 were restructured during the humidity cycling process, any changes were not reflected in the overall particle size as this was dominated by the contribution from the coating organics. To 30 further examine these possibilities, the organic coatings of the coated BC particles, that had either experienced the humidity cycling process or not, were removed by the TD, and then the size distribution of the BC core was measured by the AAC and SMPS. As shown in Fig. 9, 690 a clear difference in BC core size was observed between the organic coated BC particles that had passed immediately through the TD ( Fig. 9(a)) and those that had first experienced the humidity cycling process ( Fig. 9(b)). The mobility diameter of the BC core size peaked at 79.3 nm for the former case but peaked at 23.6 nm for the latter one, implying that the BC core was restructured (becomes more compact) after coating with organics and then 695 experiencing the high humidity environment.

AMS Mass spectra of organic aerosols
Owing to the complex nature of biomass combustion in the natural environment, and the high sensitivity to the combustion environment, the organic mass spectral signatures from previous AMS experiments vary considerably. This variability makes it difficult to estimate the 700 contribution of biomass burning aerosol (BBA) to total PM through positive matrix factorisation (PMF) methods (Paglione et al., 2020). The mass spectra for the BBA generated from our repeatable and controlled laboratory experiments have the potential to tighten constraints on the AMS spectral signature(s) for BBA. Highly controlled and reproducible measurements of aerosol emissions from combustion of a common African biofuel source has 705 been achieved and reported by Haslett et al. (2018) previously. Here, we expand detail of the mass spectra observations reported by Haslett et al. (2018) by providing data sets on the combustion for different types of biomass to that studied in previous work.
The mass spectra of organic aerosols generated from the pyrolysis and smouldering phases (noise to signal ratio is less than 0.005), and the difference mass spectra between these two 710 phases, are shown in Fig. 10 and Fig. 11 for the combustion of dry and wet wood samples, 31 respectively. For all experiments, irrespective of the wood type and moisture content, the spectra for the pyrolysis and smouldering phases are dominated by hydrocarbon ion fragments, including C n H 2n-1 (m/z=41, 55), C n H 2n+1 (m/z=29, 43, 57), C n H 2n-3 (m/z=67) and C 5+n H 5+2n (m/z=77, 91). These peaks are associated with fragments of saturated alkenes, 715 alkanes, cycloalkanes and aromatic compounds, respectively. Our results are similar to the mass spectra reported by Haslett et al. (2018) for the controlled combustion of biomass in the laboratory. By comparing Fig. 10

725
The peaks at m/z 57, 60 and 73 are seen together in both combustion phases; these peaks are often used as markers for biomass burning organic aerosols (Alfarra et al., 2007), as very few other aerosol sources in the natural environment can contribute to these peaks. Interestingly, for the wood samples of Scots pine, Giant Redwood and Western red cedar, these peaks at m/z 57, 60 and 73 were more dominant in the smouldering phase. This is a similar 730 observation to that reported by Haslett et al. (2018). However, for the Poplar wood sample, the peaks are most dominant in the pyrolysis phase. This is particularly interesting because it suggests that the combustion of angiosperm wood (Poplar) and gymnosperm wood (Scots Pine, Giant Redwood and Western red cedar) will emit different chemical species depending on the combustion phase. This difference may be due to the resin content and composition (which is related to chemical volatility) within the wood. Future work should consider the impact of burning of different wood types and their respective relationship to combustion phase.
For all wood samples, the peak at m/z 44, which is mainly attributed to the CO 2 + fragment with a possible contribution from C 2 H 4 O + , is more prominent in pyrolysis phase than the 740 smouldering phase, which is in contrast with the results from Haslett et al. (2018). This could be a result of the different wood types. We note that the study of Haslett et al., (2018) used Rubber wood (Hevea brasiliensis), that produces latex sap. This hard tropical wood also has a bulk density much higher than the woods tested here (0.63 x 10 3 kg/m 3 ) and likely a higher lignin content where lignin is one the last compounds to be broken down by the combustion 745 process.

Conclusions
During the Soot Aerodynamic Size Selection for Optical properties (SASSO) project and a EUROCHAMP transnational access project, the physical and chemical properties of black carbon and organic matter from different combustion sources were investigated. For BC 750 particles from the diesel engine and flame burner emissions, the TEM images show that the primary spherules size are similar, around 20 nm. As particle size increases, BC particles adopted more collapsed/compacted morphologies from the diesel engine but tend to more aggregated morphologies from the flame burner source. For the particles emitted from the wood combustion, the χ of BC particles and the of organics were observed to range from 755  Processed data are available at https://doi.org/10.6084/m9.figshare.15172572.v1. Raw data is archived at the University of Manchester and is available on request.

Competing financial interests
The authors declare no competing financial interests.    (1), (2), (3) and (4) represent different measurement configurations used to enable characterisation of specific aerosol physiochemical parameters, as described in the main text. The instrument set (1) was used for measurements of aerosol optical properties (and thereby enable the retrieval of refractive index for BC/BrC aerosols, subject of a future publication) and organic chemical measurement for wood combustion (setup (a)) and chamber (setup (c)) experiments; instrument set (2) was used for dynamic shape factor measurements for BC and material density measurement for organic aerosols (setup (a) and (b)); the instrument set (3) was used for SP2 incandescence signal calibration (setup (b)); the instruments set (4) was used for BC restructuring experiments (setup (c)). organic particles (smouldering phase) emitted from wood combustion: (a) mass distribution measured by CPMA; (b) particle number distribution measured by the SMPS after the CPMA selected particles at the mass peak of the D a =200 nm particle mass distribution presented in (a). (green line) before and after they experience 'humidity cycling' (i.e. exposed to 90% RH around 10 s and then dehydrated to 10% RH) process. Da is the aerodynamic diameter.