Stable sulfur isotope measurements to trace the fate of SO 2 in the Athabasca oil sands region

. Concentrations and δ 34 S values for SO 2 and size-segregated sulfate aerosols were determined for air monitoring station 13 (AMS 13) at Fort MacKay in the Athabasca oil sands region, northeastern Alberta, Canada as part of the Joint Canada-Alberta Implementation Plan for Oil Sands Monitoring (JOSM) campaign from 13 August to 5 September 2013. Sulfate aerosols and SO 2 were collected on ﬁlters using a high-volume sampler, with 12 or 24 h time intervals. Sulfur dioxide (SO 2 ) enriched in 34 S was exhausted by a chemical ionization mass spectrometer (CIMS) operated at the measurement site and affected isotope samples for a portion of the sampling period. It was realized that this could be a useful tracer and samples collected were divided into two sets. The ﬁrst set includes periods when the CIMS was not running (CIMS-OFF) and no 34 SO 2 was emitted. The second set is for periods when the CIMS was running (CIMS-ON) and 34 SO 2 was expected

Abstract.Concentrations and δ 34 S values for SO 2 and sizesegregated sulfate aerosols were determined for air monitoring station 13 (AMS 13) at Fort MacKay in the Athabasca oil sands region, northeastern Alberta, Canada as part of the Joint Canada-Alberta Implementation Plan for Oil Sands Monitoring (JOSM) campaign from 13 August to 5 September 2013.Sulfate aerosols and SO 2 were collected on filters using a high-volume sampler, with 12 or 24 h time intervals.
Sulfur dioxide (SO 2 ) enriched in 34 S was exhausted by a chemical ionization mass spectrometer (CIMS) operated at the measurement site and affected isotope samples for a portion of the sampling period.It was realized that this could be a useful tracer and samples collected were divided into two sets.The first set includes periods when the CIMS was not running (CIMS-OFF) and no 34 SO 2 was emitted.The second set is for periods when the CIMS was running (CIMS-ON) and 34 SO 2 was expected to affect SO 2 and sulfate highvolume filter samples.
δ 34 S values for sulfate aerosols with diameter D > 0.49 µm during CIMS-OFF periods (no tracer 34 SO 2 present) indicate the sulfur isotope characteristics of secondary sulfate in the region.Such aerosols had δ 34 S values that were isotopically lighter (down to −5.3 ‰) than what was expected according to potential sulfur sources in the Athabasca oil sands region (+3.9 to +11.5 ‰).Lighter δ 34 S values for larger aerosol size fractions are contrary to expectations for primary unrefined sulfur from untreated oil sands (+6.4 ‰) mixed with secondary sulfate from SO 2 oxidation and ac-companied by isotope fractionation in gas phase reactions with OH or the aqueous phase by H 2 O 2 or O 3 .Furthermore, analysis of 34 S enhancements of sulfate and SO 2 during CIMS-ON periods indicated rapid oxidation of SO 2 from this local source at ground level on the surface of aerosols before reaching the high-volume sampler or on the collected aerosols on the filters in the high-volume sampler.Anticorrelations between δ 34 S values of dominantly secondary sulfate aerosols with D < 0.49 µm and the concentrations of Fe and Mn (r = −0.80 and r = −0.76,respectively) were observed, suggesting that SO 2 was oxidized by a transition metal ion (TMI) catalyzed pathway involving O 2 and Fe 3+ and/or Mn 2+ , an oxidation pathway known to favor lighter sulfur isotopes.
Correlations between SO 2 to sulfate conversion ratio (F (s)) and the concentrations of α-pinene (r = 0.85), βpinene (r = 0.87), and limonene (r = 0.82) during daytime suggests that SO 2 oxidation by Criegee biradicals may be a potential oxidation pathway in the study region.

Introduction
Sulfate aerosols are known to impact ecosystems and climate through their deposition and radiative effects.The deposition of sulfate aerosols can cause acidification of soils and lakes (Gerhardsson, 1994).Furthermore, their direct and indirect Published by Copernicus Publications on behalf of the European Geosciences Union.
i previously published under the name Roghayeh Ghahremaninezhad radiative effects can change the radiative budget at regional scales and alter climate (IPCC, 2001).
Sulfate aerosols can be primary or secondary.Primary particles are emitted directly from the surface to the atmosphere but secondary particles are formed in the atmosphere through gas to particle conversion.The majority of anthropogenic and natural sulfur is emitted as sulfur dioxide (SO 2 ) or oxidized to SO 2 in the atmosphere (Berresheim et al., 1995;Berresheim, 2002;Seinfeld and Pandis, 1998).Chin and Jacob (1996) and Chin et al. (2000) estimated that around 50 % of the globally emitted SO 2 is oxidized to form sulfate and the remainder is lost by dry and wet deposition.
Dry deposition is important and gives SO 2 a lifetime of about 3 days for a boundary layer with 1000 m depth (Hicks, 2006;Myles et al., 2007).Wet deposition is important intermittently for rainy days or days with fog.The lifetime of SO 2 in the atmosphere can vary greatly from hours to days depending on measurement location, season, time of day, etc.As an example, Hains (2007) measured SO 2 lifetime in the eastern US and found values of 19 ± 7 h.GEOS-Chem simulations suggest a value of 13 h during summer for the same location.
A detailed understanding of SO 2 oxidation pathways and their relative importance is critical for accurate representation of sulfate's spatial distribution as well as its impact on climate through aerosol radiative forcing.
The oil sands regions are of great interest because of the large quantities of SO 2 emissions (Fioletov et al., 2016;McLinden et al., 2012;Percy, 2013).Therefore, a comprehensive knowledge of SO 2 oxidation pathways important in this region is useful to identify where and how atmospheric sulfur species are transported and contribute to aerosol formation, growth, and acid deposition.
Oil sands extraction and upgrading processes can be a source of sulfate aerosols, SO 2 , and oxidants.The major sources of SO 2 emissions in the Athabasca oil sands region are upgrading and energy production operations (Kindzierski and Ranganathan, 2006).Simpson et al. (2010) observed SO 2 enhancements over the oil sands region with a maximum value of 39 parts per billion by volume (10 −9 ppb) relative to a background value of 102 parts per trillion by volume (10 −12 ppt).Howell et al. (2014) showed that both SO 2 and sulfate contributions from the Athabasca oil sands region are significant compared to estimates for potential background sources of sulfur such as annual forest fire emissions in Canada.Bardouki et al. (2003) by the use of positive matrix factorization (PMF) modeling suggested that secondary sulfate is the second most important contributor to PM 2.5 mass in Fort MacKay (31 %).
Sulfur dioxide is converted to sulfate in homogeneous and heterogeneous reactions.The oxidation pathway is a very important factor to determine the effects of the sulfate formed on the environment.Gas phase oxidation of SO 2 by hydroxyl radicals (OH) produces sulfuric acid (H 2 SO 4 ) gas, which can nucleate in the atmosphere to form new particles (Tanaka et al., 1994;Kulmala et al., 2004).These newly formed aerosol particles are buoyant and can be dispersed far from the emission source.Newly formed sulfate aerosols also impact direct radiative forcing by scattering sunlight back to space.These particles can grow by the addition of organics to create a large number of accumulation mode aerosols, which are more easily deposited on local surfaces, increasing the potential for acidification at regional to local scales.They also have the ability to form cloud condensation nuclei (CCN; Kulmala et al., 2004Kulmala et al., , 2007;;Benson et al., 2008).After forming CCN they can increase the albedo and lifetime of clouds (Twomey, 1991;Boucher and Lohmann, 1995).Homogeneous oxidation of SO 2 in the gas phase by OH is as follows (Burkholder et al., 2015): A range of 17 to 36 % of global sulfate production can be attributed to this pathway (Chin et al., 2000;Sofen et al., 2011;Berglen, 2004).Heterogeneous oxidation of SO 2 primarily occurs in cloud droplets, although oxidation on the surface of aerosols can be important regionally (Chin and Jacob, 1996).Heterogeneous oxidation prevents H 2 SO 4 gas production and new particle formation.Sulfate formed by this pathway can modify the aerosol size distribution, which affects both direct and indirect aerosol forcing.Scattering efficiency of the particle population can be increased, which is responsible for direct scattering (Hegg et al., 2004;Yuskiewicz et al., 1999).In addition, acidity of aerosols as well as their CCN activity of the particle population can be modified and affect the indirect radiative forcing (Mertes et al., 2005a, b).Eriksen et al. (1972) showed various steps in SO 2 dissolution before oxidation by major oxidants, these are H 2 O 2 , O 3 , and O 2 catalyzed by transition metal ions (TMIs) such as Fe 3+ or Mn 2+ in a radical chain reaction pathway (Herrmann et al., 2000).(Seinfeld and Pandis, 1998).
Field studies suggested that TMI-catalyzed oxidation is the dominant sulfate formation pathway in polluted environments in winter (Jacob et al., 1984(Jacob et al., , 1989;;Jacob and Hoffmann, 1983).Oxygen isotope measurements of sulfate aerosols collected at Alert, Canada (82.5 • N, 62.3 • W) showed that TMI-catalyzed SO 2 oxidation is significant during winter (McCabe et al., 2006).Recent studies have shown that the TMI-catalyzed oxidation pathway is underestimated (more than an order of magnitude) in all current atmospheric chemistry models (Harris et al., 2013a, b).For example, Harris et al. (2013a) measured the sulfur isotopic composition of SO 2 upwind and downwind of clouds and used the difference to calculate the fractionation that occurred for in-cloud SO 2 oxidation.They showed that SO 2 oxidation catalyzed by natural TMIs on mineral dust is the dominant in-cloud oxidation pathway and is underestimated by more than an order of magnitude in current atmospheric models.To the best of our knowledge there is no study to investigate the importance of the TMI-catalyzed pathway in SO 2 oxidation on the surface of aerosols in highly polluted areas such as the Alberta oil sands region during summer.
Until recently, OH-radical-initiated oxidation of SO 2 was considered the only gas phase oxidation pathway important in the atmosphere.However, recent measurements of the rate constants for oxidation of SO 2 by Criegee biradicals and model simulations of field observations have shown this pathway is more significant than previously thought (Berndt et al., 2012;Boy et al., 2013;Mauldin III et al., 2012;Sipilä et al., 2014).Criegee biradicals are formed through ozonolysis of unsaturated hydrocarbons such as biogenic terpenes (Boy et al., 2013;Welz et al., 2012).The rate constants of the reaction of Criegee biradicals and SO 2 are somewhat uncertain but researchers agree that the reaction is faster than what has been previously thought (e.g., 6 ×10 −13 and 8 ×10 −13 cm 3 molecule −1 s −1 for Criegee biradicals originating from the ozonolysis of α-pinene and limonene, respectively, Mauldin III et al., 2012).Several studies have linked biogenic volatile organic compound (BVOC) environments to an increase in SO 2 to sulfuric acid and/or sulfate conversion rates.For example, Mauldin III et al. (2012) reported the oxidation of SO 2 by Criegee biradicals faster than what has been thought before, during a field study in a boreal forest and confirmed the results by laboratory and theoretical studies.
In this study, we investigated the importance of the various SO 2 oxidation pathways, including Criegee biradicals in a polluted region with high volatile organic compound (VOC) emissions using measurements of sulfates, SO 2 concentrations, and isotopic composition.
Sulfur isotope analysis is a powerful tool to investigate SO 2 oxidation pathways in the atmosphere.As an example, Lin et al. (2017) used high-sensitivity measurements of cosmogenic 35 S in SO 2 and sulfate from the ambient boundary layer over coastal California and the Tibetan Plateau to identify oxidation of SO 2 to sulfate.The lifetime in summer ranged from 1 to 2 days suggesting that there might be oxidation pathways which are more important than previously thought.In this study, stable sulfur isotope values for SO 2 and sizesegregated sulfate aerosols were measured.δ 34 S values of potential sources in the region (Proemse et al., 2012a) and isotope fractionation data (Harris et al., 2012) were used to investigate atmospheric sulfur oxidation pathways in the Athabasca oil sands region.The sulfur dioxide to sulfate con- was also used as a tool to investigate the possible SO 2 oxidants in the region.Although the data represent a short period of time and do not reflect the variability on a seasonal timescale, Soares et al. (2018) showed that short-term measurements are more suitable for source identification.They mentioned that the source signals of NO 2 and SO 2 emissions are available in hourly to daily timescales and long-term observation may cause a loss in short term variation.

Study site
Sulfate aerosols and SO 2 measurements were made at a monitoring site next to the Wood Buffalo Environmental Association's (WBEA) air monitoring station 13 (AMS 13) site just south of Fort MacKay in the Athabasca oil sands region from 13 August to 5 September 2013 as part of the Joint Canada-Alberta Implementation Plan for Oil Sands Monitoring (JOSM) project (Liggio et al., 2016;Phillips-Smith et al., 2017).The location of AMS 13 is shown in Fig. 1.Stable sulfur isotopes can be used to investigate sulfur sources, transport, and chemistry such as the relative importance of oxidation pathways (Puig et al., 2008;Krouse and Grinenko, 1991).Sulfur has four stable isotopes: 32 S, 33 S, 34 S, and 36 S with relative abundances of ∼ 95, 0.75, 4.21 and 0.015 %, respectively.The isotopic composition is described using the delta notation: where n is the number of atoms, x S is the heavy isotope and V-CDT is the international sulfur isotope standard, Vienna Canyon Diablo Troilite, with the isotopic ratio of R 34 = 34 S 32 S = 0.044163, R 33 = 33 S 32 S = 0.007877 (Ding et al., 2001) and R 36 = 36 S 32 S = 1.05 × 10 −4 .For the purpose of this paper we only analyze δ 34 S values and use δ 33 S values to find enrichment of samples.
In addition to sulfur emissions from oil sands processing, aerosols can potentially be produced from vehicle exhaust.Combustion emissions from vehicles showed a δ 34 S of +5 ‰ for SO 2 from engine exhaust in Alberta and British Columbia (Norman et al., 2004;Norman, 2004).On average, diesel and gasoline contained very low amounts of sulfur (0.008 %, Norman, 2004) and combustion produces both primary sulfate as well as SO 2 .Other sulfur emissions in the region may result from anoxic conditions in the environment or the tailing ponds associated with sulfate-reducing bacteria.Biogenic emissions of hydrogen sulfide (H 2 S) have negative δ 34 S values which can be as negative as −30 ‰ (Wadleigh and Blake, 1999).H 2 S is oxidized to SO 2 with a lifetime of 1 day (Brimblecombe et al., 1989) and the sulfur isotopic composition is not expected to change during oxidation of H 2 S to SO 2 (Sanusi et al., 2006;Newman et al., 1991).
Differing isotopic contributions from sulfur sources can drive variations in aerosol sulfate δ 34 S values.Another reason for δ 34 S variation can be isotopic fractionation.The oxidation of SO 2 causes isotope fractionation between the products and reactants as long as the reaction is not complete.When the reactant is available as an infinite reservoir, the fractionation factor is calculated as where R = 34 S 32 S .Following the definition for α used by Harris et al. (2012) for both kinetic and equilibrium reactions, α < 1 means that the light isotopes react faster, so products are isotopically lighter than the reactant.During this study, minute quantities of 34 SO 2 were emitted from a chemical ionization mass spectrometer (CIMS) exhaust 50 m away from the high-volume sampler near the ground for special periods.Here we refer to these particular periods as CIMS-ON.The enrichment of 34 SO 2 was sufficiently large that isotopic fractionation can be neglected during CIMS-ON periods.However, sulfur sources and oxidation pathways can be examined using δ 34 S values for the periods when CIMS was not operational (CIMS-OFF).During SO 2 oxidation to sulfate, isotope fractionation occurs between reactants and products which is unique for each oxidation pathway.Note that sulfur isotope fractionation resulting from oxidation by Criegee biradicals is not currently known.Harris et al. (2012) reported temperature dependent fractionation factors for different SO 2 oxidation pathways as follows: SO 2 oxidation by OH radicals favors heavy isotopes and the fractionation decreases slightly with temperature (Eq.3).
Aqueous phase oxidation can occur by H 2 O 2 and O 3 , and fractionation during this pathway (Eq.4) also prefers heavy isotopes and decreases with temperature slightly.
The fractionation during the TMI-catalyzed oxidation pathway acts in the opposite direction to the other two pathways.TMI-catalysis is the only known oxidation pathway which favors lighter isotopes in the product sulfate and the fractionation strongly depends on temperature (Eq.5). (5) 4 Methods

Field measurements
Temperature, relative humidity, and wind speed and direction time series are shown in Fig. A1 in the Appendix.A diurnal cycle in relative humidity (RH) is evident for all days during the campaign except 25 August which was a rainy period.
A high-volume sampler placed at ground level with a flow rate of 0.99 ± 0.05 m 3 min −1 was used to collect aerosols and SO 2 .The high-volume sampler was fitted with a fivestage cascade impactor to collect size-segregated aerosols on glass fiber filters in five ranges of aerodynamic diameter as A (>7. 2 µm), B (3.0-7.2 µm), C (1.5-3.0 µm), D (0.95-1.5 µm), and E (0.49-0.95 µm).The final filter for fraction F < 0.49 µm was a 20.3 cm × 25.4 cm glass filter to collect aerosols with D < 0.49 µm.An SO 2 filter pretreated with potassium carbonate (K 2 CO 3 ) and glycerol solution was located beneath these six size-segregated aerosol filters (Norman, 2004).The sampling interval was 12 h (daytime 05:00 to 17:00 MDT (Mountain Daylight Time) and nighttime 17:00 to 05:00 the next day) for the first 12 days except 20 and 27 August after which samples were collected for 24 h (05:00 to 05:00).Field blanks were collected on three separate occasions at the start, in the middle, and at the end of the campaign.Filter blanks from the field were loaded and then unloaded, stored, and analyzed using the same protocols as samples.The high-volume sampler was turned off during field blank sampling.Filters were stored in ziplock bags and kept at temperatures less than 4 • C and transferred to the lab for analysis.
WBEA SO 2 data were used with a sampling interval of 5 min.Ozone and NO 2 mixing ratios were measured by UV absorption using a Thermo 49i O 3 monitor every 10 s and a blue diode laser cavity ring-down spectrometer every 1 s, respectively (data were averaged to 1 min; Odame-Ankrah, 2015; Paul and Osthoff, 2010).The slope uncertainties in these measurements were ±1 and ±10 %, respectively.Radiometer measurements using a pair of spectral radiometers (one facing the zenith, the other the nadir direction) were used to determine actinic flux and to calculate photolysis frequencies (j values;Osthoff et al., 2018).Iron (Fe) and Manganese (Mn) were measured by semi-continuous X-ray fluorescence measurements of metals taken every hour on a filter tape with a measurement uncertainty of ±10 % (Phillips-Smith et al., 2017).
Monoterpenes were measured hourly by gas chromatography ion-trap mass spectrometry (GC-IT-MS; Tokarek et al., 2017).VOCs and C 2 -C 12 were sampled in canisters over a period spanning 09:30 to 08:30 of the next day, and analyzed using gas chromatography mass spectrometery (GCMS).Detection limits for VOC measurements can be found in the online JOSM database (ftp://arqpftp:research@ftp.tor.ec.gc.ca/OS/AMS13, last access: 10 October 2017).
A chemical ionization mass spectrometer (CIMS) similar to the one described by Sjostedt et al. (2007) was used to measure OH reactivity at a distance of 10 m horizontally from the high-volume sampler.Enriched 34 SO 2 was emitted from an exhaust pipe at ground level less than 50 m to the east in an unused area containing shrubs.Enriched 34 SO 2 affected a portion of our samples during CIMS-ON periods; these periods were used to trace the fate of local 34 SO 2 emitted from the CIMS exhaust near the ground.

Analysis of high-volume filter samples
Filter papers were shredded and sonicated for 30 min in distilled deionized water in the laboratory (200 mL for SO 2 filters and filters to collect particles in the size range F < 0.49 µm and 75 mL for slotted filters to collect particles in sizes larger than 0.49 µm).For SO 2 filters, 1 mL of 30 % w/w hydrogen peroxide (from BDH) was added to oxidize the SO 2 to sulfate before sonication.Filter paper fibers were removed by 0.45 mm Millipore filtration, and 10 mL of the filtrate samples was analyzed using a Dionex ICS-1000 ion chromatography (IC) system with a Dionex IonPac AS14 column and electric conductivity detector to determine the concentration of sulfate with an uncertainty of 5 %.Prior to treatment, the pH of the remaining filtrate was measured and found to be ∼ 6.0.The remaining filtrate was treated with 0.5 mL of 10 % BaCl 2 (dihydrate 99 %, from EMD), and dilute (0.5 normal) OmniTrace HCl (34-37 %, from EMD) was added to samples until a pH of 3 was achieved.Approximately 100 µL of 0.5 normal HCl was used for aerosol filters.The samples were then heated to facilitate precipitation of BaSO 4 .Barium sulfate was isolated by Millipore filtration, and dried samples were packed into tin cups and analyzed with a PRISM II continuous-flow isotope ratio mass spectrometer (CF-IRMS) to obtain δ 34 S values (relative to V-CDT; Giesemann et al., 1994).The precision in measuring δ 34 S is ±0.3 ‰ which is determined as the standard deviation (1σ ) of δ 34 S for several standard runs.δ 34 S measurements were blank corrected using the sulfur concentration and δ 34 S values for field blanks.Insufficient sulfate was present for some samples after concentration blank correction.Although the concentration of sulfate was too small to perform blank correction for some samples, they displayed the same range for δ 34 S values as those which were blank corrected.This suggests little to no bias was introduced by blank correction.Therefore, δ 34 S values are reported from some samples which were not isotopically blank corrected.These samples are indicated with a * in Tables 1 and 2.
The PRISM II continuous flow isotope ratio mass spectrometer measures δ 34 S and δ 33 S simultaneously and the values for non-enriched samples were expected to be related according to the mass dependent fractionation (MDF) relation (δ 33 S ∼ 0.51δ 34 S).For this experiment, some of the samples were enriched in 34 S and they were identified by the use of the MDF relation between δ 34 S and δ 33 S of the standards for the same run.δ 33 S/δ 34 S was averaged for standards for each run and δ 33 S δ 34 S − 2σ was used as a cutoff criterion and data falling below this criterion were tagged as enriched.
Care was taken to analyze sufficient standards and blanks between enriched samples (CIMS-ON periods) to ensure carryover was minimal.Little to no deviation in standards and is also lost by wall loss in the instrument so the majority of the exhaust is in the form of 34 SO 2 .For the periods when the CIMS was operational (CIMS-ON), significant 34 S isotope enrichment was observed; therefore, samples were divided into two sets, CIMS-ON and CIMS-OFF.
The first set is for samples collected during the shutdown periods of the CIMS (CIMS-OFF).These CIMS-OFF periods were used to investigate the isotopic composition of sizesegregated sulfate aerosols and SO 2 in the region and the possible sources and formation pathways of sulfate aerosols.The second set (CIMS-ON) is for samples affected by enriched 34 S and is not used as an indicator of sulfur isotopic composition of sulfate aerosols in the region.Instead, the enriched 34 SO 2 is used as a natural tracer to follow the fate of SO 2 emitted from a local ground-based source and its oxidation.

Sulfur conversion ratio
In this paper we use the sulfur conversion ratio, which is defined as the portion of SO 2 which is converted to particulate sulfate: In this formula, [SO 4 ] is the concentration of sulfate aerosols with D < 0.49 µm.In this study the sulfate is dominantly secondary (Proemse et al., 2012a), corroborated here by the absence of soil indicators (Sect.5.1).F (s) can be affected by dry deposition.Since little is known about the appropriate dry deposition velocities in this region, potential variations between SO 2 and sulfate dry deposition rates are neglected in the analysis.Since F (s) is a measure of SO 2 to sulfate conversion, it is a measure of oxidant loading.Therefore, significant posi-tive correlation between F (s) and other compounds may be an indicator of the importance of that compound as a tracer for SO 2 oxidation.This formula can be used for both CIMS-ON and CIMS-OFF periods since the number of enriched molecules reaching the high-volume sampler is very small and cannot change F (s).The number of enriched molecules reaching the high volume sampler is calculated using equations described in Sect.4.3.3 and the fraction of enriched molecules in comparison to the total sulfur concentration is reported in Table A1.

Concentration of 34 S enriched molecules
The concentration of enriched molecules as 34 SO 2 and 34 SO 4 were calculated using the following equations during CIMS-ON periods.Isotope ratio (R) values show the ratio of sulfur isotopes to the most abundant isotope, which is 32 S for sulfur.

R 34
= n 34 S/n 32 S, ( 7) R 34 enriched = (n 34 S + n 34 S * )/n 32 S, ( 10) in which n 34 S * is the number of 34 S atoms reaching the filter from the CIMS exhaust and S total is the total number of sulfur atoms on the filter.The R 34 value is calculated as the average of R 34 values for samples without enrichment.There were R 33 data available from the IRMS but the uncertainty was high (±3 ‰) and we used the value for the international standard for sulfur V-CDT. 36S is included in calculations since the amount of 34 S * from the CIMS exhaust is on the same order of magnitude.R 34 enriched values were available for each sample.The concentration of sulfate for each sample was available from IC and the number of sulfur atoms as SO 2 or sulfate can be calculated.Then the number of 34 S from CIMS was calculated and divided by the volume of total sampled air and the number of 34 SO 2 * and 34 SO * 4 molecules cm −3 was calculated (Table A1).

Sulfur conversion ratio (F (s))
The sulfur conversion ratio (F (s), Eq. 6) was calculated for the smallest size fraction of measured sulfate (F < 0.49 µm ).Absence of Ca and Mg in this size fraction (all concentrations were below the IC detection limit of 0.1 mg L −1 ) indicates that primary soil particles were not present in this size fraction.Proemse et al. (2012a) suggested that less than 10 % of total sulfur emissions from two major stacks in the region in PM 2.5 were primary sulfate.Therefore, primary sulfate from stacks do not form a significant portion of sulfate aerosols in the D < 0.49 µm size range.Based on these two pieces of information, it is expected that sulfate particles on this size fraction are mostly (> 90 %) secondary.As a result, F (s) gives valuable information about which pathways dominate SO 2 oxidation and formation of sulfate aerosols.F (s) is not affected by enriched sulfate emissions during CIMS-ON periods (because the amount of 34 SO 2 emitted was relatively small, Table A1).Hence, F (s) reflects the conversion of SO 2 to sulfate for the entire measuring period.This implies negligible changes to F (s) values because of the CIMS emissions.
F (s) values were usually higher during the daytime in comparison to nighttime values (Tables 1 and 2), which was what we expected for OH-driven oxidation during daylight.In the troposphere, the OH radical is produced mainly from photolysis of O 3 to O( 1 D) and subsequent reaction with water vapor.If a steady state in O( 1 D) is assumed with respect to its production and loss, the (instantaneous) daytime OH production rate is proportional to A negative correlation was observed between F (s) and this (integrated) OH production rate during the daytime (r = −0.72,P value < 0.05; Fig. A2).However, two data points with the highest RH (25 and 26 August) drive this correlation, and no correlation was observed for the remainder of the samples.This suggests that there may be SO 2 oxidation pathways in addition to OH during the day in this region.
The time series for SO 2 during the campaign is shown in Fig. 3a.The time series was dominated by spikes in the SO 2 mixing ratio.Phillips-Smith et al. ( 2017) used PMF to determine concentration time series for five factors during the campaign.This analysis showed that on the 14, 23, and 24 August and 3 and 4 September were periods that the site was impacted by upgrader emissions.Concentrations of SO 2 , Fe, and Mn (measured in PM 2.5 ) were markedly higher during these periods (Fig. 3).
It is interesting to note that F (s) for daytime was higher than nighttime for all samples except periods when the site was impacted by plumes from major oil sands upgrading facilities (polluted periods; Phillips-Smith et al., 2017).A comparison between AM and PM values for F (s) for 23 and 24 August showed that nighttime values were almost double the daytime values.F (s) data were not available for the daytime of 14 August to compare with the nighttime value, but 14 August PM showed the highest value for F (s) (0.77) during the entire campaign (Tables 1 and 2).At night, aqueous phase oxidation is believed to be the dominant SO 2 transformation pathway as OH is absent (Chin and Jacob, 1996).No correlation was observed between the F (s) and O 3 mixing ratio for daytime, nighttime, or daily samples (Fig. A2).Therefore, it is expected that SO 2 oxidation occurs by the H 2 O 2 and/or the TMI-catalyzed pathways.
Fe and Mn concentrations in PM 2.5 aerosols, averaged over the nighttime high-volume sampling periods, are shown in Fig. A3 (Phillips-Smith et al., 2017).The averaged nighttime concentrations of Fe and Mn were higher during polluted periods (average values of 57 ± 20 ng m −3 and 1.5 ± 0.5 ng m −3 , respectively) in comparison to other periods (average values of 9 ± 3 ng m −3 and 0.13 ± 0.06 ng m −3 , respectively; Fig. A3).The data collected on 21 August PM were excluded from this analysis because the PMF analwww.atmos-chem-phys.net/18/1/2018/Atmos.Chem.Phys., 18, 1-24, 2018 Concentrations of Fe and Mn were associated with upgrader, soil, and haul road dust factors during polluted nighttime periods (14, 23, 24 August: Phillips-Smith et al., 2017).Although 21 August PM was not a polluted period, it showed high [Fe + Mn] × [H 2 O] values but F (s) was not high (Fig. 4b).For 21 August, the analysis by Phillips-Smith et al. (2017) showed that there was a peak for the soil factor but not upgrader and haul road dust.F (s) on this night was markedly lower than during periods when upgrader and haul road dust factors were high.F (s) for 25 August was also high since this was a rainy period (Sect.4.1).

δ 34 S values for size-segregated sulfate aerosols and SO 2 during CIMS-OFF periods
During CIMS-OFF periods 34 SO 2 emissions were absent, so δ 34 S values reflect the sulfur isotopic composition of the sulfur compounds in the region and/or fractionation as the SO 2 is oxidized and transported to the AMS 13 site.δ 34 S values during CIMS-OFF periods for SO 2 and size-segregated sulfate in size ranges F < 0.49 µm , E 0.49−0.95µm , D 0.95−1.5 µm , C 1.5−3.0µm, B 3.0−7.2µm , and A > 7.2 µm are shown in Table 1.
Possible oxidation pathways of SO 2 to sulfate were investigated using these δ 34 S values.Blank corrected δ 34 S values for SO 2 were +5.1 and +10.8 ‰.No negative δ 34 S values were observed for SO 2 .If it is assumed that no fractionation occurred during formation of primary sulfate in major stacks, then it is expected that δ 34 S values for SO 2 would be the same as primary sulfate (with an average of +7.3 ± 0.3 ‰ and +9.4 ± 2.0 ‰).The δ 34 S values of SO 2 ranged from +5.1 to +11.1 ‰ (Table 1) and are consistent with this assumption.The lowest value (+5.1 ‰) is consistent with a δ 34 S value for SO 2 from vehicle exhaust (Table 1).
δ 34 S values for size F < 0.49 µm particles ranged between +1.8 and +15.1 ‰ with an average of +7.4 ± 4.2 ‰.Although this average overlaps with values given by Proemse et al. (2012b) for primary sulfate from the stack emissions (+7.3 ± 0.3 and +9.4 ± 2.0 ‰), there were δ 34 S values lighter and heavier than what was expected from potential sulfur sources in the region in this size range.Therefore, δ 34 S of sulfate cannot be used as a quantitative indicator for industrial SO 2 emissions as isotope fractionation may have occurred as the stack emissions (SO 2 ) were transported to the AMS 13 site.As shown in Sect.5.1 sulfate particles in this size range are predominantly secondary; therefore, these data can be used to investigate the importance of different SO 2 oxidation pathways during transport.Particles in larger size ranges (E 0.49−0.95µm , D 0.95−1.5 µm , C 1.5−3.0µm , B 3.0−7.2µm , and A > 7.2 µm ) are expected to contain more primary sulfate and have lower δ 34 S values in comparison to the F < 0.49 µm size range.There were no nega-  tive values for sulfate particles in the size fraction F < 0.49 µm , but negative values were observed for the size fraction E 0.49−0.95µm .There was a tendency to lighter δ 34 S values for larger sulfate particles as shown in Fig. 5.

Correlation between Fe and Mn and sulfate concentration and δ 34 S values during CIMS-OFF periods
Sulfur dioxide can be oxidized in the aqueous phase by O 2 in the presence of TMIs predominantly by Fe 3+ and Mn 2+ (Herrmann et al., 2000).If this is an important oxidation pathway, more secondary sulfate is expected to be produced when the concentrations of catalysts are higher.Since concentrations of Fe and Mn were measured in PM 2.5 particles, the concentration of sulfate in impactor size fractions (F < 0.49 µm , E 0.49−0.95µm , D 0.95−1.5 µm , and C 1.5−3.0µm ) were added to find the concentration of sulfate in particles with D < 3 µm.The sulfate concentration in the impactor size range D < 3 µm is almost the same as the concentration in PM 2.5 since the concentration in size fraction C 1.5−3.0µm was very low (zero for all periods except polluted periods, which ranged between 0.58 and 1.76 µg m −3 ).The concentration for particles from the impactor with D < 3 µm is plotted versus Fe and Mn concentrations and the sum of Fe and Mn in Fig. 6.Positive correlations were observed for all three cases (r = 0.86, r = 0.89, r = 0.86, respectively; Fig. 6).Positive correlations were also observed when the concentrations of sulfate in the aerosol size fractions F < 0.49 µm and E 0.49−0.95µm were plotted against the concentrations of Fe and Mn and the sum of Fe and Mn (Fig. A4).There were not enough sulfate concentration data for size fractions C 1.5−3.0µm and D 0.95−1.5 µm to show the individual correlations with Fe, Mn, and the sum of Fe and Mn.When SO 2 is oxidized by the TMI-catalyzed pathway, the sulfur isotopic composition of the sulfate formed is lighter than the isotopic composition of the reactant SO 2 (Harris et al., 2012).Significant anti-correlations were apparent for sulfate δ 34 S values in the size fraction F < 0.49 µm when plotted against Fe and Mn concentrations (r = −0.80 and r = −0.76,respectively; Fig. 7).This suggests that lighter δ 34 S values occur in secondary sulfate in the presence of higher concentrations of Fe and Mn.Insufficient isotope data were available to create similar plots for other size fractions.
Positive correlations were also observed between concentrations of Fe and Mn and the concentration of SO 2 (r = 0.67 and r = 0.65, respectively; Fig. A5), which may indicate that they originate from the same source, or were transported together to the sampling site.

δ 34 S values of SO 2 and size-segregated sulfate aerosols during CIMS-ON periods
The release of 34 SO 2 from the CIMS allowed for an examination of SO 2 oxidation to sulfate under field conditions.An unexpected result was found: δ 34 S values for SO 2 and sulfate samples with D < 0.49 µm during the periods when the CIMS was operated (CIMS-ON) are shown in Table 2.The blank corrected data show that δ 34 S values for enriched SO 2 samples were only as high as +35.6 ‰, and there were values without enrichment ranging between +4.8 and +10.9 ‰ with an average value of +8.3 ± 1.8 ‰.All sulfate samples in the size range F < 0.49 µm representing SO 2 oxidation during the CIMS-ON periods were blank corrected, and all AM and PM samples were highly enriched in 34 S; the δ 34 S values were as high as +913 ‰ (Table 2).
A comparison between the isotopic composition of sulfate aerosols in the size range F < 0.49 µm and SO 2 samples (R SO 4 /R SO 2 ) showed that the sulfate particles with D < 0.49 µm were much more enriched in 34 S from the 34 SO 2 tracer released by the CIMS.The concentration of enriched sulfur as 34 SO 2 and 34 SO 4 molecules cm −3 is also calculated as described in Sect.4.3.3 and the data are reported in Table A1.
Since the CIMS exhaust was located to the southeast of the high-volume sampler, wind direction was considered as a potential factor in the analysis.No correlation (r = 0.16) was observed between the percent of time the high-volume sampler was downwind of the CIMS exhaust and the concentrations of 34 SO 2 or 34 SO 4 .

The role of Criegee biradicals in SO 2 oxidation
As mentioned in Sect.5.1, F (s) was higher during the daytime in comparison to nighttime except for polluted periods.No correlation (r = −0.36,excluding 25 and 26 August with the highest RH) was observed between F (s) and the inte- grated OH production rate, suggesting that another oxidation pathway for SO 2 was active during daytime.One likely pathway is oxidation of SO 2 by Criegee biradicals.Criegee biradicals are formed from ozonolysis of alkenes and may oxidize SO 2 to sulfate increasing F (s) (Mauldin III et al., 2012).Therefore, it is expected that correlations may exist between F (s) and precursors to Criegee biradicals.Positive correlations between F (s) and the concentration of α-pinene (r = 0.85), β-pinene (r = 0.87), and limonene (r = 0.82) were observed during daytime (Fig. 8).However, no correlations were observed between F (s) and monoterpenes during nighttime.
The concentration of monoterpenes showed a negative correlation with the mixing ratio of O 3 .There was a power law relationship between monoterpenes and O 3 mixing ratio during the daytime and a linear dependency at night (r = −0.60;Fig. 9).
Correlations with aromatic compounds generally fell into two categories.The first set includes compounds which were highly correlated with SO 2 and sulfate for impactor D < 0.49 µm (e.g., benzene).The second set contains the compounds which show no such correlations but were correlated with F (s) (Table A3).Styrene and p-cymene are two compounds with no correlation with SO 2 and sulfate but significant correlations with F (s) (r = 0.58 for both, and r = 0.66 and r = 0.71, respectively, when the rainy day data are omitted, P value < 0.05; Table A3, Fig. A6).They also show a positive correlation together (r = 0.71, P value < 0.05; Fig. A6).

Potential TMI-catalyzed SO 2 oxidation
Sulfur conversion ratios (F (s)) and sulfur isotope data for SO 2 and size-segregated sulfate aerosols were used to investigate the role of TMI-catalyzed SO 2 oxidation in the region.
Figure 2 exhibits the expected correlation between F (s) and OH in daytime but not at night when aqueous phase reactions are important.The similar slopes for nighttime and daytime F (s) versus RH plots suggests that SO 2 aqueous phase oxidation may be an important oxidation pathway for both day and night and the offset (intercept that is higher for daytime than nighttime) suggests there is additional gas phase SO 2 oxidation that takes place during the day.
Known aqueous phase oxidants for SO 2 are H 2 O 2 , O 3 , and O 2 in the presence of TMIs (Herrmann et al., 2000).No correlation was observed between F (s) and O 3 mixing ratios, which suggests that O 3 is of minor importance as an oxidant in the aqueous phase.The pH dependency of aqueous O 3 oxidation of SO 2 makes this reaction very slow at low pH (pH < 5.5).This reaction is also self-limiting and production of sulfate lowers the pH and slows down further reaction (Seinfeld and Pandis, 1998).Therefore, aqueous phase oxidation of SO 2 occurs mostly by H 2 O 2 and/or O 2 in the presence of TMIs.
The conversion ratio of SO 2 to sulfate (F (s)) was higher during the day than at night except during polluted periods (14,23,and 24 August).This is consistent with gas phase contributions to SO 2 oxidation in addition to aqueous phase oxidation that occurred both during the day and at night (Sect.5.1).On polluted nights, the SO 2 to sulfate conversion ratio was twice as high as during the day and on 14 August at night the highest (0.77) conversion ratio for the entire campaign was observed (Tables 1 and 2).Averaged Fe and Mn concentrations on these polluted nights coincided with the highest values for SO 2 to sulfate conversion (F (s); Fig. A3).
Whenever both RH and the sum of Fe and Mn concentrations were high at night, the proportion of SO 2 that was converted to sulfate (F (s)) was higher as well (Fig. 4).These conditions of coincident high RH and Fe + Mn concentrations was met on polluted nights during the campaign.On these nights the ratios of Fe / Mn were around 40 (38,40,42 for 14,23,24 August,respectively).This specific ratio may be a useful indicator for the source of Fe and Mn in aerosols.A particular night that was not classified as polluted (21 August) was identified as having Fe and Mn from soil (Phillips-Smith et al., 2017).The ratio of Fe / Mn on that night was 76, and the SO 2 to sulfate conversion ratio was indistinguishable from the remainder of the non-polluted nighttime samples (Fig. 4).Therefore, it is reasonable to suggest that a Fe / Mn value around 40 is associated with a non-soil source.The two remaining sources are upgrader emissions and haul road dust.This interpretation of Fe and Mn on polluted nights as originating from anthropogenic emissions (Fe / Mn ∼ 40) rather than soil is consistent with the higher solubility of anthropogenic TMIs relative to soil (Kumar et al., 2010).
Sulfur isotope measurements can provide the means to distinguish TMI from H 2 O 2 aqueous oxidation.Isotope fractionation will be evident in sulfate when a large reservoir of SO 2 (e.g., from stack emissions) mixes with oxidants during transport and produces accumulated sulfate product captured over 12 or 24 h.So long as the fraction of reaction is low (< 30 %) the difference in δ 34 S values for SO 2 and sulfate will reflect the magnitude and direction of the fractionation process.For the TMI-catalyzed pathway this direction is negative and produces lighter sulfate than SO 2 .This directly contrasts with fractionation for O 3 , H 2 O 2 , and OH oxidation pathways.Evidence that SO 2 released from tall stacks is transported high above the ground and mixes down toward the surface at AMS 13 has been demonstrated by Gordon et al. (2017) and should provide conditions meeting the requirement for the fraction of reaction less than 30 % described here.The observed δ 34 S values for size-segregated sulfate aerosols in this study were consistent with aqueous TMI rather than H 2 O 2 oxidation.Light δ 34 S values for sulfate aerosols were observed in the region in comparison to other potential atmospheric sulfur sources.An alternate explanation for isotopically light δ 34 S values in sulfate was proposed by Proemse et al. (2012b).Isotopically light δ 34 S values (−3.9 and +0.3 ‰) were reported by this group for sulfate from bulk and throughfall deposition (deposition of excess water onto the ground surface from wet leaves) in the Athabasca oil sands region, consistent with the observations in this study.Since these values were lighter than the potential sources in the region, this suggests a contribution of sulfate from a 34 S depleted source.Proemse et al. (2012b) suggested that the low δ 34 S values observed for atmospheric sulfate collected at two sites were due to H 2 S emitted from tailing ponds.Tailing ponds were in close proximity to the two sites where low δ 34 S values were found.Proemse et al. (2012b) suggested that H 2 S was oxidized to SO 2 and subsequently formed sulfate that then contributed to local sulfate deposition.The average value for δ 34 S of SO 2 during CIMS-ON and CIMS-OFF (non enriched values) periods was +7.9 ± 2.1 ‰.This value is in the range of δ 34 S of primary sulfate from two major stacks (Proemse et al., 2012a).No negative values were observed for δ 34 S of SO 2 .If H 2 S was the main source of atmospheric sulfur, the opposite pattern to that observed in Fig. 5, is expected.The reason is that isotopically light SO 2 from H 2 S oxidation is expected to produce secondary sulfate aerosols (from both homogeneous and heterogeneous reactions) in the smaller size fractions (F < 0.49 µm and E 0.49−0.95µm ) with isotopically light δ 34 S values.The larger A > 7.2 µm and B 3.0−7.2µm size aerosols contain primary sulfate from soil and would reflect δ 34 S values for untreated oil sand (+6.4 ‰; Proemse et al., 2012a) in addition to sulfate from H 2 S oxidation so they would have progressively more positive δ 34 S values.Therefore, a discernable contribution of H 2 S to isotopically light samples through an SO 2 oxidation pathway is ruled out.
Primary sulfate and SO 2 can originate from haul road dust or diesel exhaust.δ 34 S values for these two sources are +5 ‰ and higher (Norman, 2004;Norman et al., 2004).Therefore, if haul road dust and diesel primary sulfate were transported with Fe and Mn, then δ 34 S values should converge to +5 ‰ or higher.This should be particularly evident for the larger size aerosols (A > 7.2 µm and B 3.0−7.2µm ).In fact the opposite is observed in Fig. 5. Isotopically light δ 34 S values for sulfate aerosols in size ranges E 0.49−0.95µm , D 0.95−1.5 µm , C 1.5−3.0µm , B 3.0−7.2µm , and A > 7.2 µm were observed during CIMS-OFF periods.These values indicate that there was no, or only a very small, contribution of primary sulfate from major stacks.This leaves SO 2 from upgrader emissions as the most probable source of sulfate both for F < 0.49 µm size aerosols and for secondary sulfate formed on larger aerosol size fractions.δ 34 S values reflect isotope fractionation during oxidation of SO 2 rather than source signatures.This is supported by a positive correlation between the sum of sulfate in size fractions F < 0.49 µm , E 0.49−0.95µm , D 0.95−1.5 µm , and C 1.5−3.0µm and the concentrations of Fe and Mn and sum of Fe and Mn (r = 0.86, r = 0.89, and r = 0.86, respectively).The concentration of sulfate in size fractions F < 0.49 µm and E 0.49−0.95µm also showed positive correlations with the concentration of Fe and Mn.This indicates that when Fe and Mn were prevalent in aerosols, either more sulfate can be formed or Fe and Mn were transported to AMS 13 with SO 2 from a common emission source, likely upgrader emissions.There were also anti-correlations between δ 34 S values of sulfate in size fraction F < 0.49 µm and the concentrations of Fe and Mn.This shows that lighter δ 34 S values were associated with secondary sulfate formation and higher concentrations of Fe and Mn.One possible explanation for these observations may be the TMI-catalyzed SO 2 oxidation pathway during transport to the AMS 13 site.

CIMS-ON
Little 34 SO 2 reached the SO 2 filter in the high-volume sampler since high sulfur isotope enrichment was not observed for SO 2 samples (max δ 34 S = +35.6 ‰).Instead 34 SO 2 was oxidized to sulfate either as it moved in the atmosphere or on the filters in the high-volume sampler.This result was unexpected since previous studies of δ 34 S for sulfate and SO 2 showed no evidence of oxidation when SO 2 passed through the filters under marine or continental conditions (Ghahremaninezhad et al., 2016).The lack of 34 SO 2 and the predominance of 34 S molecules on sulfate aerosols demonstrates an oxidation pathway that is rapid and specific to the conditions at the ground level of the the AMS 13 site.

Potential oxidation of SO 2 by Criegee biradicals
The proportion of sulfate from SO 2 oxidation, F (s), during daytime is generally larger than F (s) at night (Tables 1  and 2).Greater vertical mixing is expected during the day than at night.Stack emissions high above ground (Gordon et al., 2017) undergo oxidation during transport to the AMS 13 site.Aloft, conventional oxidation pathways (i.e., OH-driven oxidation) are likely more important than near the surface.At the same time precursors to Criegee biradicals will be released and mixed upward.A larger F (s) during the day than at night suggests that during daytime gas phase SO 2 oxidation occurs in addition to aqueous phase oxidation.Typically, OH is expected to dominate gas phase SO 2 oxidation during the day.However, a correlation between F (s) and integrated OH production rate was not observed.Instead, positive correlations between F (s) and α-pinene, β-pinene and limonene were observed during the day but not at night (Fig. 8).This, combined with the loss of monoterpenes as daytime O 3 mixing ratio increased, suggests Criegee biradicals may be an important factor in SO 2 oxidation close to the surface during daytime.Monoterpenes are oxidized by O 3 to form Criegee biradicals which can be stabilized and oxidize SO 2 to form secondary sulfate.This pathway is potentially more important during the day but less so at night.At night, the emissions of monoterpenes continue into a shallow nocturnal boundary layer that is decoupled from the residual layer above it.The terpenes then titrate O 3 at the surface, leading to the observed anti-correlation and low surface O 3 mixing ratio which limits Criegee biradical production.
Reaction between O 3 and anthropogenic alkenes may also generate Criegee biradicals, potentially leading to higher SO 2 to sulfate conversion ratios (F (s)).Many anthropogenic alkenes and aromatics likely have sources in common with SO 2 since a correlation (P value < 0.05) was observed between them (Tables A2 and A3).Their emissions are likely injected into (and transported within) layers above the measurement site and only sporadically entrain to the surface during daytime.When this happens, relationships between F (s) and anthropogenic alkenes may be observed.As an examwww.atmos-chem-phys.net/18/1/2018/Atmos.Chem.Phys., 18, 1-24, 2018 ple, styrene and p-cymene did not correlate with SO 2 or secondary sulfate but they were correlated with F (s) (r = 0.66, r = 0.71, respectively).Styrene and p-cymene were also highly correlated with each other (r = 0.71) suggesting they originated from the same source or sources.It is likely that styrene and p-cymene are indicators of other anthropogenic alkenes that facilitate SO 2 oxidation (for instance, tetrachloroethene).

Conclusions
This is the first study to examine oxidation of SO 2 as it is transported above and within the boundary layer at AMS 13, a highly polluted environment, during summer.Sulfur dioxide (SO 2 ) and size-segregated sulfate aerosol concentrations and sulfur isotope compositions were measured during summer 2013 in the Athabasca oil sands region to investigate SO 2 oxidation pathways.δ 34 S values, F (s), and the relationship between secondary sulfate concentrations and Fe and Mn (in PM 2.5 ) show that there is the potential that a significant proportion of SO 2 is oxidized rapidly during both the day and at night.Aqueous phase oxidation by TMI catalysis is consistent with these results.The fraction of secondary sulfate was higher during the night than during the day for periods when the site was impacted by industrial plumes mixing downward from above.This, taken together with the high Fe and Mn concentrations in PM 2.5 at night, shows the importance of aqueous phase reactions, probably by the TMI pathway as SO 2 is transported from the stack to the site at night.In addition, a natural tracer experiment with enriched 34 S demonstrated that oxidation of SO 2 on the surface of aerosols is rapid.The results would be consistent with Criegee biradicals being an important daytime oxidation pathway for SO 2 at ground level, which was suggested in several recent high-profile papers (Mauldin III et al., 2012;Boy et al., 2013;Sipilä et al., 2014).

Figure 1 .
Figure 1.The Wood Buffalo Air Monitoring Station 13 (AMS 13) site, south of Fort MacKay (map data© 2018 Google)

4N.
Amiri et al.: Stable sulfur isotope measurements to trace the fate of SO 2 3 Sulfur isotopes

Figure 3 .
Figure 3. (a) SO 2 time series with a sampling interval of 5 min (there is a gap in 14 August data) and (b) hourly data for Mn (left axis) and Fe (right axis).Shaded areas indicate polluted periods.

FFigure 4
Figure 4. (a) ([Fe] + [Mn]) × [H 2 O] values for nighttime (PM) samples as an indicator of the TMI-catalyzed SO 2 oxidation pathway (Fe and Mn concentrations were averaged over the running periods of the high-volume sampler) and (b) F (s) values for the nighttime samples.Polluted periods and the soil episode are shown by gray and yellow shaded areas, respectively.

Figure 9 .
Figure 9. Sum of α-pinene, β-pinene, and Limonene versus ozone mixing ratio for daytime and nighttime and all data.

FigureFigure A2 .Figure A3 .Figure
Figure A1.(a) Temperature and relative humidity data with 1 min sampling interval.(b) Wind speed and wind direction with the sampling time interval of 5 min (data from WBEA meteorological station AMS 13).The gray shaded areas show the CIMS-ON periods.
Amiri et al.:Stable sulfur isotope measurements to trace the fate of SO 2 blanks was apparent after enriched δ 34 S values from CIMS-ON periods were analyzed.In this paper uncertainties are reported as 1σ standard deviation.=99and m / z = 49 in the negative ion spectrum of the mass spectrometer.An excess amount of 34 SO 2 compared to the required 34 SO 2 to complete titration of OH in the sample flow was used for ambient air OH reactivity measurements.Almost all of the flow entering was exhausted by the instrument which contained excess34SO 2 and formed H 34 2 SO 4 .In 1 min, n34 SO 2 = (7.4× 10 6 )n H 34 2 SO 4 .Some of the formed H 34 2 SO 4 www.atmos-chem-phys.net/18/1/2018/Atmos.Chem.Phys., 18, 1-24, 2018N.
* not blank corrected samples.

Table 2 .
δ 34 S values (‰) for SO 2 and sulfate with diameter D < 0.49 µm, ratio of sulfate to SO 2 isotope during CIMS-ON periods, F (s) values, and the error in F (s). Enriched samples were selected by comparing the mass-dependent fractionation relation between δ 34 S and δ 33 S for the sample and standards at the same run.Average uncertainty for δ 34 S values is ±0.5 ‰.

Table A3 .
Correlation coefficients (r) between selected aromatics (at least 15 out of 20 data points are above the detection limit), SO 2 , and sulfate.P values < 0.05 are indicated with a * .