Articles | Volume 15, issue 18
Research article
15 Sep 2022
Research article |  | 15 Sep 2022

Thermal–optical analysis of quartz fiber filters loaded with snow samples – determination of iron based on interferences caused by mineral dust

Daniela Kau, Marion Greilinger, Bernadette Kirchsteiger, Aron Göndör, Christopher Herzig, Andreas Limbeck, Elisabeth Eitenberger, and Anne Kasper-Giebl

The determination of mineral dust and elemental carbon in snow samples is of great interest, since both compounds are known to be light-absorbing snow impurities. Different analytical methods have to be used to quantify both compounds. The occurrence of mineral dust, which contains hematite, leads to a bias in the quantification of elemental carbon and organic carbon via thermal–optical analysis. Here we present an approach which utilizes this interference to determine the concentration of iron via thermal–optical analysis using a Lab OC / EC Aerosol Analyzer (Sunset Laboratory Inc.) and the EUSAAR2 protocol. For this, the temperature dependency of the transmittance signal determined during the calibration phase, i.e., when all carbonaceous compounds are already removed, is evaluated. Converting the transmittance signal into an attenuation, a linear relationship between this attenuation and the iron loading is obtained for loadings ranging from 10 to 100 µg Fe cm−2. Furthermore, evaluation of the transmittance signal during the calibration phase allows to identify samples which need to be re-evaluated, since the analysis of elemental carbon and organic carbon is biased by constituents of mineral dust. The method, which was initially designed for snow samples, can also be used to evaluate particulate matter samples collected within the same high alpine environment. When applying the method to a new set of samples it is crucial to check whether the composition of iron compounds and the sample matrix remain comparable. If other sources than mineral dust determine the iron concentration in particulate matter, these samples cannot be evaluated with thermal–optical analysis. This is shown exemplarily with data from particulate matter samples collected in a railway tunnel.

1 Introduction

Black carbon (BC) and mineral dust (MD) have been identified and evaluated as prominent light-absorbing snow impurities (LASI) (Doherty et al., 2018; Kaspari et al., 2014; Svensson et al., 2018; Tuzet et al., 2017; Warren and Wiscombe, 1980). Deposition of these compounds leads to a reduction in the surface albedo, influencing thawing and geochemical processes, accelerating snow and ice melt, and triggering albedo feedback (Ramanathan and Carmichael, 2008; Tuzet et al., 2020).

The presence of MD in snow can be determined via gravimetric measurements (e.g., Kuchiki et al., 2015) and particle counting methods (e.g., Di Mauro et al., 2015) or it is calculated based on the concentrations of Fe, Ca, Al, Si, and Ti (Li et al., 2017).

BC concentrations in snow and ice have been quantified through several methods that address either BC or elemental carbon (EC). Thermal–optical analysis (Birch and Cary, 1996; Chow et al., 1993), optical methods determining the spectral absorption of total particulate matter sampled on a filter (Grenfell et al., 2011), and laser-induced incandescence (Schwarz et al., 2012) are commonly applied. Regarding the differentiation between BC and EC, we refer the reader to the detailed discussion in Petzold et al. (2013). Reviews and comparisons of the advantages and disadvantages of the respective methods can be found in the literature (Kang et al., 2020; Lim et al., 2014; Qian et al., 2015). In this work we will focus on thermal–optical analysis (TOA), which makes it possible to differentiate between organic carbon (OC) and EC and represents the reference method for analysis of OC and EC in particulate matter samples (DIN e.V., 2017). The filter is heated in an inert (He) and subsequently in an oxidizing (He / O2) atmosphere, while the evolving carbon is converted to methane and determined with a flame ionization detector. To account for pyrolysis and to separate OC and EC precisely, the transmittance or reflectance of the filter is monitored. The time when the transmittance (or reflectance) reaches its initial value determines the split point for differentiating between OC and EC. Several temperature protocols were introduced to optimize the method (e.g., Cavalli et al., 2010; Chow et al., 2007). Liquid samples can be analyzed after filtration, when water-insoluble OC (WinsOC) and EC are retained on a quartz fiber filter. Previous works address uncertainties related to the filtration process such as an undercatch during the filtration, uneven filter loadings, or the loss of EC within vessels used for melting and filtration (Forsström et al., 2013; Kuchiki et al., 2015; Lim et al., 2014; Meinander et al., 2020; Torres et al., 2014; Wang et al., 2020).

Uncertainties related to TOA itself apply both to the analyses of filter residues of liquid samples and to the analyses of particulate matter samples, although there might be some basic differences between the two types of samples, such as a differing particle size distribution or a slightly modified chemical composition. Since we focus on the analytical procedure of TOA, we will avoid the term “WinsOC” throughout the paper and stick to “OC”. The concept of analyses and evaluation can be applied to the analysis of both sample types, i.e., particulate matter filters (OC is addressed) or filters loaded with liquid snow samples (WinsOC is addressed).

The OC / EC split and thus the quantification of both compounds, EC and OC, are influenced by the overall chemical composition. When elevated loadings of MD are present, the influence of carbonates (Karanasiou et al., 2011; Li et al., 2017) and metal oxides (Bladt et al., 2012, 2014) is important, as MD is rich in these compounds (e.g., Kandler et al., 2007). Here we focus on Fe oxides, which are the main light-absorbing minerals in MD and cause its reddish (hematite) or yellowish (goethite) color (Alfaro et al., 2004). Formenti et al. (2014) found Fe oxides (goethite and hematite) to account for 2 %–5 % of MD by mass, investigating the mineralogical composition of mineral dust in western Africa. At elevated temperatures, between 250 and 600 C, goethite changes to hematite (Liu et al., 2013). This temperature is already exceeded during the first part of TOA (He atmosphere), making hematite the main Fe oxide present on the filters during the time when the split point for OC / EC differentiation is set. Thus, Fe oxides do not only determine the properties of MD as a prominent component of LASI, but they also affect TOA. To exceed the method's limit of detection for OC and EC, filtration of a rather large amount of liquid snow or ice is desirable. By doing so, MD, containing Fe oxides, is enriched and interferences occur. Wang et al. (2012) used TOA for analyses and identified an extra decrease in optical reflectance during the 250 C heating stage in hematite and suggested shifting the reference value for the OC / EC split to this temperature step. Gul et al. (2018) suggested enhancing these adjustments even further, defining the reference value for the OC / EC split at 550 C for filters with higher MD loads.

In this work we investigate the influence of MD loads on TOA of snow samples collected in a high alpine environment. The main focus of our work is to investigate the temperature dependence of the light attenuation caused by Fe-containing compounds in MD. This led us to a new approach for approximating Fe concentrations via TOA. Therefore, we evaluate the transmittance signal during the calibration phase, i.e., when carbonaceous compounds are no longer present. The method is tested extensively for snow samples and evaluated briefly for particulate matter samples.

2 Collection of samples and sample preparation

2.1 Snow samples

Surface snow samples were collected on the upper platform of the global GAW (Global Atmosphere Watch) station Sonnblick Observatory, situated at 3106 m a.s.l., in the Austrian Alps. MD deposited on Mount Sonnblick often derives from the long-range transport of dust from desert regions (Greilinger et al., 2018). For sampling, the uppermost centimeters (approximately 5 cm) of snow were collected directly in Whirl-Pak™ bags (Nasco™). For this purpose, the bag was opened and the opening, which is reinforced by a wire handle, was used as a shovel to scoop in the top snow layer. Samples with and without a visible influence of MD were collected during multiple visits at the station between 2018 and 2020 and were kept frozen until further processing. In the laboratory, samples (approximately 600 g each) were melted gently in a glass beaker using a microwave (600 W). Aliquots of the liquid samples were filtrated onto precleaned (24 h at 600 C, cooling in a desiccator above distilled water) quartz fiber filters (Pallflex® Tissuquartz™; loaded area was circular with a diameter of 16 mm; half was used for analysis) applying a slight vacuum. After filtration, these snow sample filters were dried overnight in a desiccator above silica gel to facilitate further analysis.

Blank filters were prepared in a similar way, filtrating ultrapure water instead of the liquid snow sample.

2.2 PM10 samples

PM10, i.e., particles with an aerodynamic diameter equal to or less than 10 µm, was sampled onto quartz fiber filters (Pallflex® Tissuquartz™) in the same high alpine environment as the snow samples using a high-volume sampler (DIGITEL Elektronik AG, Switzerland; sampling duration: 1 week) and in tunnels with railway traffic using a low-volume sampler SEQ47/50 (Sven Leckel Ingenieurbüro GmbH, Germany; sampling duration: 4 h). No precleaning step of the filters was performed before sampling. Rectangular aliquots with an area of 1.5 cm2 were used for analysis.

2.3 Reference samples

Hematite (Fe2O3) was chosen as a reference substance for light-absorbing compounds in MD since goethite and hematite are the most abundant forms of Fe oxides in MD (Formenti et al., 2014) and goethite is converted to hematite at elevated temperatures starting at 250 C and being completed at 600 C (Liu et al., 2013), i.e., at temperatures which correspond to the conditions used in the inert phase of TOA already. Fe2O3 was suspended in ultrapure water (8.75 mg in 1 L flask) in an ultrasonic bath. Quartz fiber filters were loaded with different amounts of the suspension to cover the range of the Fe loadings found on filters loaded with snow samples containing MD (roughly 9–314 µg Fe cm−2). The procedure of suspending the particles and loading the filters was done on 2 different days to verify reproducibility.

A different set of samples was prepared with Standard Reference Material® 2709 (SRM 2709 San Joaquin Soil, National Institute of Standards and Technology, USA). SRM 2709 was suspended in ultrapure water (83.3 mg in 1 L flask) and quartz fiber filters were loaded covering the range of roughly 3–141 µg Fe cm−2.

3 Instrumentation

3.1 Thermal–optical analysis

Thermal-optical analysis (TOA) is the reference method for the determination of OC and EC in ambient aerosols (DIN e.V., 2017). Filter aliquots are heated in an inert and oxidizing atmosphere, the evolving carbonaceous compounds are converted to methane, which is quantified via a flame ionization detector. To correct for pyrolytic carbon, the transmittance of the filter aliquot is recorded. As explained in detail in the results section, the development of the transmittance signal is influenced by the Fe loading of the filter. Therefore, we solely investigate this transmittance signal in this paper. A Lab OC-EC Aerosol Analyzer (Sunset Laboratory Inc., USA) and the EUSAAR2 protocol (Cavalli et al., 2010; temperature steps are given in Table S1 in the Supplement) were used. The instrument logs the transmittance and reflectance of the sample at a wavelength of 660 nm during the measurement. We used the program OCEC834 for the measurement and the program Calc415 for the evaluation of the raw data (both Sunset Laboratory Inc., USA). Further processing of the transmittance data was conducted with an external program (Microsoft Excel, Microsoft Corporation, USA).

3.2 Inductively coupled plasma–optical emission spectroscopy and inductively coupled plasma–mass spectrometry

Inductively coupled plasma–optical emission spectroscopy (ICP-OES) and inductively coupled plasma–mass spectrometry (ICP-MS) were used to quantify Fe. For this, microwave-assisted digestion of the filter aliquots was carried out subsequent to TOA. Filters loaded with snow samples or reference material (Fe2O3, SRM 2709) samples and high alpine PM10 samples were digested using a microwave system (Multiwave 5000, Anton Paar, Austria) and METHOD 3052 (method for microwave-assisted acid digestion suitable for siliceous, organic, and other complex matrices; US EPA, 1996; HNO3 : HCl : HF 6 : 2 : 1, maximum temperature: 180 C) and were analyzed for Fe via ICP-OES (iCAP 6500 ICP-OES spectrometer, Thermo Scientific, USA). Filter aliquots of the particulate matter samples collected within the railway tunnel were digested using a microwave system (Multiwave 3000, Anton Paar, Austria/Start1500, MLS GmbH, Germany; aqua regia, maximum temperature: 220 C) and Fe was quantified via ICP-MS (iCap Q System instrument, Thermo Scientific, USA). The limit of detection was 0.4 µg cm−2 for ICP-OES and 0.1 µg cm−2 for ICP-MS.

3.3 X-ray powder diffraction

X-ray powder diffraction (PXRD) was used for specification of Fe present on the snow sample filters. Experiments were carried out in Bragg–Brentano geometry using an Empyrean diffractometer (Malvern PANalaytical B.V., Netherlands; scattering angle range of 5< 2θ< 135). A focusing mirror was used to provide Cu Kα1,2-radiation for the experiment. The beam divergence was defined by using a 1/4 fixed vertical entrance slit followed up by a 0.04 rad horizontal Soller slit and a 0.04 rad horizontal Soller slit on the secondary side in front of an open line detector (GaliPix detector). The detector to sample distance for this instrument was fixed to 240 mm.

The PXRD diagrams were evaluated using the Malvern PANalytical program suite HighScorePlus v4.6a (Degen et al., 2014). A background correction and a Kα2 strip were performed. Crystallographic phases were assigned based on the ICDD-PDF4+ database (Kabekkodu et al., 2002).

3.4 Scanning electron microscopy

Scanning electron microscopy (SEM) measurements allowed us to visualize particle sizes and their elemental composition. Thus, a comparison of filters loaded with snow samples and the reference substances Fe2O3 and SRM 2709 was possible. An FEI Quanta 200 (Thermo Fisher Scientific, USA) instrument equipped with an Octane Pro EDS System (EDAX, USA) was used for the analyses of loaded quartz fiber filters. Samples were coated with Au prior to analyses.

4 Results

4.1 Transmittance of filters loaded with snow samples during TOA

During TOA, the transmittance of the filters changes due to pyrolysis and removal processes of carbonaceous compounds, resulting in a clean filter after the analysis is completed. However, if MD containing colored Fe compounds, e.g., hematite, is present, the filters remain colored and the transmittance remains lower than for a clean filter, even after the complete removal of carbon. Thus, the transmittance, and the respective changes of this value observed throughout TOA, can either be due to the carbon load or due to constituents of MD. The influence of the various compounds becomes visible when the same filter is analyzed repeatedly, as illustrated in Fig. 1. The change in transmittance of a filter loaded with a snow sample without visible MD contamination during TOA is presented in Fig. 1a, while Fig. 1b shows the results obtained for a snow sample filter containing MD. For both samples, multiple measurements (three to five reruns) of the same filter were conducted to separate the change in transmittance deriving from carbonaceous compounds and from MD, which is left on the filter as a residue after the first analytical run. The course of the transmittance signal in the first run, when carbonaceous compounds are still present, is given as a black line. The average course of transmittance of the reruns is shown as a gray line, whereas the standard deviation of the reruns is plotted as a gray area. The reported temperature is the temperature inside the sample oven. The vertical lines denote the switch from an inert (He) to an oxidizing atmosphere (He / O2) and the start of the calibration phase (Cal). Before introducing O2 and starting the second temperature profile, the temperature is reduced temporarily. At the end, the thermogram includes a cooling phase. At that time, calibration is performed. The part of this cooling phase, which will be used for further evaluations, is highlighted in orange. The Fe loading on the filter loaded with MD was 30 µg Fe cm−2, while it was below 2.5 µg Fe cm−2 for the sample without visible MD contamination. This corresponds to Fe concentrations in the liquid snow sample of 1.1 mg Fe L−1 (sample containing MD, 54 mL of snow filtrated) and 51 µg Fe L−1 (sample without MD, 99 mL of snow filtrated). The total carbon loads on the filters were 66 µg C cm−2 (sample containing MD) and 28.5 µg C cm−2 (sample without MD), but for better clarity, the FID signal is not shown in any panel of Fig. 1.

Figure 1Transmittance during TOA of quartz fiber filters loaded with a snow sample without (a) and with (b) mineral dust contamination, a blank filter (c) and a filter loaded with Fe2O3 (d). Details of the filter loadings are given in the text.


Both samples, with and without MD, show pronounced changes in the transmittance signal during the first analysis represented by the black line. Changes are more pronounced in Fig. 1b, which can be expected due to the higher total carbon load. For the sample without MD (Fig. 1a), the split point could be set automatically and was reached during the He / O2 3 phase, i.e., the third temperature step during the oxidizing phase. For the sample containing MD (Fig. 1b), the initial transmittance was never reached, and thus no EC would be determined, while OC would be overestimated. Applying the approach described by Wang et al. (2012), i.e., using the transmittance at 250 C as a reference if MD is present, still no split point can be determined before the cooling phase starts. Using the approach of Gul et al. (2018) and taking the transmittance at 550 C as a reference, an identification of EC becomes possible. For completeness, we added a thermogram showing the FID signal of the sample containing MD in the Supplement (Fig. S1).

During the reruns (gray line) only minor changes in the transmittance signal are determined for the sample without MD (Fig. 1a), which is in good agreement with the recordings obtained for the blank filter (Fig. 1c). The sample containing MD (Fig. 1b) behaves differently. Pronounced changes of the transmittance are visible for the reruns as well. The transmittance decreases whenever the temperature increases, and vice versa. Obviously, these changes of transmittance are due to constituents of the remaining MD and not induced by pyrolysis or combustion of carbonaceous compounds. Since hematite has already been identified as the main light-absorbing compound in MD, several filters were loaded with hematite and analyzed as reference samples. To confirm the suitability of hematite as a reference compound, PXRD measurements were performed before and after TOA for a subset of the filters loaded with snow containing MD. Hematite was clearly identifiable in all samples, but due to the high background of the filter matrix and the rather low dust loadings, a more detailed specification of hematite could not be obtained. The changes in transmittance of Fe2O3 during TOA are shown exemplarily for a filter loaded with 19 µg Fe cm−2 in Fig. 1d. Again, replicate analyses (three reruns) were performed. The results of the first run are not shown, because marked differences in the transmittance signal were observed between the first measurement and the reruns, which could be attributed to the evolution of carbonaceous material attached to the hematite (total carbon approx. 5 µg C cm−2). The procedure to omit the first run is in accordance with Wang et al. (2012) and Yamanoi et al. (2009), who also preheated the hematite once before further analyses were conducted. The change in the transmittance signal of hematite compares well to the results of the snow sample containing MD (Fig. 1b). Again, an opposing trend of transmittance and temperature was observed. The variation spans an overall range of approximately 500 a.u. and the minimum values are found at the highest temperatures. The variation of transmittance which can be attributed exclusively to refractory and light-absorbing compounds in MD is lower than the variation observed during the first run of analyses, when carbonaceous compounds are present. Still, it is large enough to markedly affect the determination of the split point and explains why the initial value of transmittance cannot be reached at all, when MD loadings are high.

During TOA, transmittance is determined at a wavelength of 660 nm where Fe oxide minerals have a rather low absorbing potential compared to shorter wavelengths (Alfaro et al., 2004). Therefore, such a severe interference is not expected. However, these evaluations refer to ambient temperature, while high-temperature visible spectroscopy showed that the reflectance spectra and colors of hematite changed markedly with temperatures up to 800 C (Yamanoi and Nakashima, 2005; Yamanoi et al., 2009). The changes are especially pronounced at longer wavelengths (> 550 nm), meaning that the red color of hematite becomes black with temperature increase (Yamanoi et al., 2009). This phenomenon introduces a severe bias to the determination of the split point, as its position can be expected at a temperature above 500 C. At that point the darkening of the filter due to MD will be much more pronounced than at the initial temperature or even at the temperature ranges suggested earlier to correct for a potential bias of the split point (Gul et al., 2018; Wang et al., 2012) and a systematic error in the determination of EC and OC occurs. Using the transmittance at 250 C for the determination of the split point, the influence of MD still remains underestimated: at the time of the split point the temperature will definitely be higher and thus the reduction in the transmittance signal due to hematite more pronounced. It would be mere coincidence if this effect was offset by an early pyrolysis of OC. Using the transmittance at 550 C accounts better for this phenomenon; however, the difference of the black and the gray lines in Fig. 1b shows that marked amounts of pyrolytic carbon can already be formed at this point of the analysis. If this is the case it will lead to an overestimation of EC accompanied by an underestimation of OC.

The error in the determination of EC and OC cannot be quantified based on a single-run analysis and we suggest a second run to assess the relative influences of both the pyrolytic carbon and the MD. The need to perform a second run of analyses can be derived from the readings of the transmittance signal during the calibration phase. At this point of the analysis, all carbon is already removed, a wide temperature range is covered, and the gas flow through the main oven maintains oxidizing conditions, like in the second part of TOA when the split point is set. As expected, the transmittance signal is rather constant for the sample without MD contamination (orange part of the graph in Fig. 1a) and compares well to the transmittance signal of the blank filter (orange part of the graph in Fig. 1c). Samples containing MD show a marked increase in transmittance during the calibration phase, i.e., when the sample is cooling (Fig. 1b). The same effect is visible for the reference sample of hematite in Fig. 1d. The change in transmittance during the calibration phase is always visible, during the initial measurement of the sample and for the reruns. Thus, the decision for an additional run of TOA can be based on the calibration phase of the first run. If a marked increase in transmittance is visible, an interference has to be expected and we recommend reanalyzing the sample to set the split point more precisely. As explained in more detail later, a clear dependence of the transmittance signal on the temperature was observed for samples with Fe loadings above 10 µg Fe cm−2. A detailed description of a correction procedure is beyond the scope of this work, which focuses on the temperature dependence of the transmittance signal and its relationship to Fe concentrations.

For completeness, we would like to mention that high loads of MD might have other influences as well, since precombustion of EC might occur in the presence of inorganic compounds contained in MD, as was shown by Bladt et al. (2012) and Neri et al. (1997).

4.2 Using TOA for analysis of Fe in reference samples

The evaluation of the transmittance signal logged during the last 75 s of the calibration phase offers more possibilities than a qualitative judgment of whether MD might interfere with the OC / EC split point. The sample temperature in this interval changes between slightly above 700 C to slightly below 400 C. To consider the different loadings of the filters, the transmittance data were converted to an attenuation (ATN), defined in Eq. (1).

(1) ATN = 100 × ln I 0 I

This conversion follows the approach by Nicolosi et al. (2018), who defined I as the transmittance signal at a particular time step and I0 as the laser transmittance at the end of the TOA, which reflects an unloaded filter. In our case, I refers to a particular temperature rather than to a time step. Furthermore, I0 is more difficult to determine, as the initial transmittance of the unloaded filter is not reached at the end of the measurement, when the filter remains loaded with light-absorbing compounds like hematite. Using the transmittance signal at the end of the calibration phase as reference value I0, the calculated ATN becomes a relative value, expressing the changes relative to the transmittance at 400 C. This value will always be smaller than the actual ATN of the light-absorbing material. Calculating ATN700–400 based on the transmittance values at 700 C (I) and 400 C (I0), we approximate color changes of Fe-containing compounds, e.g., hematite, in this temperature interval. First, we evaluate the transmittance signal of filters loaded with two reference substances, Fe2O3 and SRM 2709. Plotting the Fe loading of the filter expressed as Fe (µg cm−2) versus ATN700–400 (Fig. 2) shows a similar trend for SRM 2709 and pure hematite. ATN700–400 values increase up to Fe loadings of approximately 150 µg Fe cm−2. At higher loadings, ATN700–400 starts to decrease again. This effect is due to the definition of ATN700–400 as a relative value, i.e., the difference between the actual ATN at 400 C (ATN400) and at 700 C (ATN700). ATN400 and ATN700 relate the transmittance at the respective temperatures to the transmittance value of the unloaded quartz fiber filter at room temperature (I0). These values could only be approximated, because we did not determine the transmittance values of the unloaded filters. The respective data are presented in the Supplement (Fig. S2). As may be expected for filter-based measurements, a linear relationship between Fe loadings and ATN could only be determined for low Fe values, but saturation occurs at higher loadings (Gundel et al., 1984). The onset of this effect is different for ATN400 and ATN700, i.e., it starts at lower Fe loadings for ATN700. Calculating ATN700–400 as the difference of ATN700 and ATN400 leads to the trend shown in Fig. 2, and to the unfavorable situation that a defined ATN700–400 can be attributed to two Fe loadings. However, this effect is of minor importance for actual analysis. Fe loadings above 150 µg Fe cm−2 are avoided, as these filters are prone to sample loss.

Figure 2Fe loading of quartz fiber filters loaded with Fe2O3 or SRM 2709 determined via ICP-OES plotted versus ATN700–400.


We attribute the differences between the filters loaded with SRM 2709 and Fe2O3 shown in Fig. 2 to differences in the overall loadings of the filters (i.e., the loading with other material than hematite) and the particle sizes. Such differences are confirmed via SEM analyses of selected filters loaded with both reference materials. While Fe2O3 filters only contained hematite, SRM filers showed different kinds of minerals, i.e., particles consist of Ca, K, Si, Al, Mg, Na, and O and just a minor contribution of Fe, which obviously leads to lower ATN values. Hematite particles on the Fe2O3 filters were generally larger in size than Fe-rich particulate matter in SRM 2709. Examples of the SEM analyses, elemental maps, and spectra are given in the Supplement (Fig. S3).

4.3 Applying TOA for the determination of Fe in snow and PM10 samples

Applying the approach described above to filters loaded with particulates in snow, a similar trend is obtained, although the snow sample filters generally yield smaller values of ATN700–400 (Fig. 3) compared to the measurements with pure hematite or SRM 2709 (Fig. 2). ATN700–400 values for snow sample filters are better comparable to the filters loaded with SRM 2709 than those loaded with Fe2O3. The chemical composition of particulates in snow is more similar to SRM 2709 than to pure hematite, as shown in Fig. S3 in the Supplement. Still, SEM images show differences in particle size of SRM and snow filters, which might induce the differences together with differences in the chemical composition, which cannot be identified via SEM.

Figure 3Fe loading of quartz fiber filters loaded with snow samples determined via ICP-OES plotted versus ATN700–400.


Restrictions have to be set to define a range which enables a definite allocation of ATN700–400 to an Fe loading. It is easy to distinguish between samples of the upper and lower part of the curve, as Fe loadings above 150 µg cm−2 always show transmittance values below 600 a.u. at a temperature of 400 C in the calibration phase. Regarding actual Fe loadings on snow sample filters, Fe loadings above 150 µg cm−2 are hardly found, as already mentioned. Additionally, an upper and lower limit for ATN700–400 values was set. The upper limit (i.e., ATN700–400< 29) excludes values with higher variation, while the lower limit (i.e., ATN700–400> 3.5) accounts for samples without MD, which do not show a distinct relationship between the transmittance and the sample temperature (see Fig. 1). Hence, these samples show very low or even negative values of ATN700–400. Based on the aforementioned restrictions (transmittance > 600 a.u. at 400 C in the calibration phase and 3.5 < ATN700–400< 29) Fe loadings from 10 to 100 µg cm−2 become accessible via TOA.

The selected data points are presented in Fig. 4. A fit was computed applying the least square approach to calculate Fe loadings based on the ATN700–400 values. The resulting function is given in Eq. (2) and features a coefficient of determination R2 of 0.9691. The residuals did not show distinctive features. The difference between the calculated and the measured Fe loadings divided by the measured Fe loading was below 10 % for 20 of 34 samples and above 25 % for seven samples, the maximum being 41 % and the median 8 %.

(2) Fe ( µ g cm - 2 ) = 2.9144 × ATN 700 400

Up to this point, we only addressed filters loaded with insoluble material filtrated from snow samples. Of course, MD also occurs in ambient aerosols and similar effects can be expected for these filters. To check whether the fit is also applicable for PM10 samples, ambient air filters obtained at the Sonnblick Observatory were analyzed and evaluated in the same way. Loadings above 10 µg Fe cm−2 were hardly achieved during the long-term data set, covering June 2017 to March 2022. Only six samples out of the data set showed Fe loadings in the appropriate range. The PM10 samples represent events with long-range transport of desert dust to the high alpine site and are in good agreement with the model obtained for the snow sample filters. These data points are also included in Fig. 4. Applying the fit calculated for the snow samples, the relative differences between the Fe loadings determined with TOA and with ICP-OES were between 3 % and 49 % with a median of 24 % for the PM10 samples.

Figure 4Linear fit for Fe loading and ATN700–400 of quartz fiber filters loaded with snow samples and comparison with high alpine PM10 samples showing sufficiently high Fe loadings.


For snow samples and PM10 samples from a high alpine environment, the evaluation of ATN700–400 enables the quantification of Fe between 10 and 100 µg cm−2 from TOA transmittance data using the same linear approach. Higher Fe loadings can be identified with TOA but should not be quantified. As the snow samples were collected over a period of 2 years, a variety of transport events is included in this evaluation. Still, we stress that different source regions might lead to differences in the composition and particle size distribution of MD, which might influence such a fit. Nevertheless, the applicability of the fit to at least a limited number of PM10 filters sampled in the same high alpine environment is a promising result.

Figure 5Fe loadings of all sample groups versus ATN700–400 (a) and Fe loading of PM10 samples from tunnel measurements versus ATN700–400 (b).


To test whether the approach is generally applicable for aerosol samples with marked contributions of Fe compounds, we included samples collected in a railway tunnel. This represents a completely different environment since very high Fe contributions and different Fe compounds, such as magnetite, hematite, and Fe metal, can be expected (Eom et al., 2013), originating from abrasion of the tracks, wheels, and brakes of trains. For the tunnel samples included in Fig. 5a, Fe accounts for up to 73 % of particulate matter mass, but no specification of Fe compounds has been performed. This high share of Fe in particulate matter mass highlights the different origin of Fe compared to the high alpine samples, as Fe accounts for less than 5 % of mass in MD only. Overall, 79 filters collected in railway tunnels were evaluated. No correlation between ATN700–400 and Fe loadings was found for these samples, although loadings reach 325 µg Fe cm−2. The transmittance at 400 C during the calibration phase was never below 600 a.u., which was the case for the snow sample filters and the filters loaded with hematite at Fe loadings above 150 µg cm−2. Only 26 samples show an ATN700–400> 3.5 and thus meet the criteria for the evaluation based on ATN700–400. All other samples featured values for ATN700–400 below 3.5, partly even negative values. Even tunnel samples with Fe loadings up to 133 µg cm−2 showed negative values for ATN700–400. All samples are presented in Fig. 5b. “Tunnel PM10, Crit” is the subset of samples which meets the criteria defined for an unambiguous evaluation of Fe-containing compounds in snow. While snow samples that do not meet the criteria featured Fe loadings either below 10 or above 100 µg cm−2, tunnel samples that do not meet the criteria are found across the whole range of Fe loadings, i.e., 2–325 µg cm−2. Despite the high Fe loading on the tunnel filters, automatic split points could be set for most of the tunnel samples. The temperature-dependent influence of the Fe loads on the transmittance signal was smaller than the changes in transmittance deriving from carbonaceous compounds. This discrepancy between snow samples and tunnel samples is due to the different origin and form of Fe present and could also be noticed via visual inspection of the filters after TOA, as they were differently colored. Photos of two filters are presented in the Supplement (Fig. S4).

The comparison of all sample groups underlines that Fe loadings can only be deduced from the ATN700–400 signal when samples contain Fe in forms as present in MD, e.g., hematite and goethite. If other forms of Fe are expected, the evaluation cannot be applied. Additionally, it shows that samples from the same environment (snow and PM10 samples from high alpine environment) exhibit a similar behavior, while MD of a different origin (SRM 2709) might need an independent fit to estimate the Fe loading on the filters.

All data presented in this work were obtained using the EUSAAR2 protocol, but only the signals determined during the calibration phase are considered. Since this temperature range is also covered by the calibration phase of other commonly used thermal protocols, e.g., IMPROVE_A or NIOSH-like protocols, the evaluations remain valid for these protocols as well, as observed for a limited set of samples not discussed here. Note that the evaluation was conducted using the transmittance and not reflectance signal. Post-processing of data already available is possible, as no adjustment of the widely used thermal protocols is necessary.

5 Conclusion

Investigating the influence of MD, i.e., its light-absorbing constituent hematite, on TOA, we presented an easily applicable approach to determine Fe loadings on filters containing MD based on temperature changes of the transmittance signal recorded during the calibration phase of TOA. The method is based on the evaluation of quartz fiber filters loaded with snow samples containing different amounts of MD.

Considering the whole range of Fe loadings investigated, the relationship between the Fe loading and ATN700–400 was shown to be non-linear. By setting two straightforward criteria (ATN700–400 between 3.5 and 29 and transmittance signal at 400 C above 600 a.u.), it is possible to determine Fe loadings between 10 and 100 µg Fe cm−2 based on a linear approach. This approach allows to obtain another compound of LASI with the same method (TOA) as used for the determination of EC. This can be advantageous if the available amount of sample is limited and a proxy for MD is needed. However, a reliable conversion factor of Fe to MD has to be deduced by independent analysis or taken from the literature and is not presented here.

The method developed for the analyses of filters loaded with snow samples was successfully applied to the analyses of particulate matter filters originating from the same high alpine environment. We cannot give information about the general applicability to PM10 samples loaded with different types of MD, yet. We hypothesize that it should be possible to use the approach as long as MD and its constituent Fe2O3 remain the dominant source of Fe loadings. Still, an independent calibration might be needed to account for differences in MD composition. Fe loadings in particulate matter samples collected in a railway tunnel and thus having a strong influence of abrasion products from tracks and wheels and not from MD cannot be analyzed with the method presented here. This illustrates the limitation of the method.

Filters loaded heavily with MD (Fe loadings above 10 µg Fe cm−2) will experience a severe bias of the OC / EC split point. EC will be underestimated or cannot be determined at all, while OC will be overestimated.

We recommend rerunning the analyses, as soon as the evaluation of the transmittance signal in the calibration phase indicates an influence of MD. This second run of analyses will make it possible to set the split point more precisely.

Data availability

ATN700–400 values and Fe concentrations of all samples presented in this work are made available in Table S2 in the Supplement.


The supplement related to this article is available online at:

Author contributions

DK and AKG designed the experiments, performed the data evaluation and interpretation, and prepared the manuscript. Analyses were conducted by DK (preparation, TOA, and digestion of high alpine, Fe2O3 and SRM samples), BK (evaluation of digestion protocols, quality assurance), AG (TOA and digestion of tunnel samples), CH (ICP), and EE (SEM). MG performed formal analyses and was responsible for the collection of part of the snow samples. Resources were provided by AKG and AL. All authors contributed to reviewing and editing of the manuscript.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


Special thanks go to Klaudia Hradil for conducting PXRD analyses at the X-Ray Center at TU Wien, Austria, and for her valuable input.

Financial support

The publication of this paper has been supported by TU Wien Bibliothek through its Open Access Funding Program.

Review statement

This paper was edited by Mingjin Tang and reviewed by Jonas Svensson, Xin Wang, and one anonymous referee.


Alfaro, S. C., Lafon, S., Rajot, J. L., Formenti, P., Gaudichet, A., and Maille, M.: Iron oxides and light absorption by pure desert dust: An experimental study, J. Geophys. Res.-Atmos., 109, D08208,, 2004. 

Birch, M. E. and Cary, R. A.: Elemental carbon-based method for monitoring occupational exposures to particulate diesel exhaust, Aerosol Sci. Tech., 25, 221–241,, 1996. 

Bladt, H., Schmid, J., Kireeva, E. D., Popovicheva, O. B., Perseantseva, N. M., Timofeev, M. A., Heister, K., Uihlein, J., Ivleva, N. P., and Niessner, R.: Impact of Fe content in laboratory-produced soot aerosol on its composition, structure, and thermo-chemical properties, Aerosol Sci. Tech., 46, 1337–1348,, 2012. 

Bladt, H., Ivleva, N. P., and Niessner, R.: Internally mixed multicomponent soot: Impact of different salts on soot structure and thermo-chemical properties, J. Aerosol Sci., 70, 26–35,, 2014. 

Cavalli, F., Viana, M., Yttri, K. E., Genberg, J., and Putaud, J.-P.: Toward a standardised thermal-optical protocol for measuring atmospheric organic and elemental carbon: the EUSAAR protocol, Atmos. Meas. Tech., 3, 79–89,, 2010. 

Chow, J. C., Watson, J. G., Pritchett, L. C., Pierson, W. R., Frazier, C. A., and Purcell, R. G.: The DRI thermal/optical reflectance carbon analysis system: description, evaluation and applications in US air quality studies, Atmos. Environ. A.-Gen., 27, 1185–1201,, 1993. 

Chow, J. C., Watson, J. G., Chen, L. W. A., Chang, M. O., Robinson, N. F., Trimble, D., and Kohl, S.: The IMPROVE_A temperature protocol for thermal/optical carbon analysis: maintaining consistency with a long-term database, JAPCA J. Air Waste Ma., 57, 1014–1023,, 2007. 

Degen, T., Sadki, M., Bron, E., König, U., and Nénert, G.: The highscore suite, Powder Diffr., 29, S13–S18,, 2014. 

Di Mauro, B., Fava, F., Ferrero, L., Garzonio, R., Baccolo, G., Delmonte, B., and Colombo, R.: Mineral dust impact on snow radiative properties in the European Alps combining ground, UAV, and satellite observations, J. Geophys. Res.-Atmos., 120, 6080–6097,, 2015. 

DIN e.V.: DIN EN 16909:2017, Außenluft – Messung von auf Filtern gesammeltem elementaren Kohlenstoff (EC) und organisch gebundenem Kohlenstoff (OC); Deutsche Fassung EN 16909:2017, Beuth Verlag GmbH, Berlin, 2017. 

Doherty, S. J., Hegg, D. A., Johnson, J. E., Quinn, P. K., Schwarz, J. P., Dang, C., and Warren, S. G.: Causes of variability in light absorption by particles in snow at sites in Idaho and Utah, J. Geophys. Res.-Atmos., 121, 4751–4768,, 2018. 

Eom, H. J., Jung, H. J., Sobanska, S., Chung, S. G., Son, Y. S., Kim, J. C., Sunwoo, Y., and Ro, C. U.: Iron speciation of airborne subway particles by the combined use of energy dispersive electron probe X-ray microanalysis and Raman microspectrometry, Anal. Chem., 85, 10424–10431,, 2013. 

Formenti, P., Caquineau, S., Desboeufs, K., Klaver, A., Chevaillier, S., Journet, E., and Rajot, J. L.: Mapping the physico-chemical properties of mineral dust in western Africa: mineralogical composition, Atmos. Chem. Phys., 14, 10663–10686,, 2014. 

Forsström, S., Isaksson, E., Skeie, R. B., Ström, J., Pedersen, C. A., Hudson, S. R., Berntsen, T. K., Lihavainen, H., Godtliebsen, F., and Gerland, S.: Elemental carbon measurements in European Arctic snow packs, J. Geophys. Res.-Atmos., 118, 13–614,, 2013. 

Greilinger, M., Schauer, G., Baumann-Stanzer, K., Skomorowski, P., Schöner, W., and Kasper-Giebl, A.: Contribution of saharan dust to ion deposition loads of high alpine snow packs in Austria (1987–2017), Front. Earth Sci., 6, 126,, 2018. 

Grenfell, T. C., Doherty, S. J., Clarke, A. D., and Warren, S. G.: Light absorption from particulate impurities in snow and ice determined by spectrophotometric analysis of filters, Appl. Optics, 50, 2037–2048,, 2011. 

Gul, C., Puppala, S. P., Kang, S., Adhikary, B., Zhang, Y., Ali, S., Li, Y., and Li, X.: Concentrations and source regions of light-absorbing particles in snow/ice in northern Pakistan and their impact on snow albedo, Atmos. Chem. Phys., 18, 4981–5000,, 2018. 

Gundel, L. A., Dod, R. L., Rosen, H., and Novakov, T.: The relationship between optical attenuation and black carbon concentration for ambient and source particles, Sci. Total Environ., 36, 197–202,, 1984. 

Kabekkodu, S. N., Faber, J., and Fawcett, T.: New Powder Diffraction File (PDF-4) in relational database format: advantages and data-mining capabilities, Acta Crystallogr. B, 58, 333–337,, 2002. 

Kandler, K., Benker, N., Bundke, U., Cuevas, E., Ebert, M., Knippertz, P., Rodríguez, S., Schütz, L., and Weinbruch, S.: Chemical composition and complex refractive index of Saharan Mineral Dust at Izaña, Tenerife (Spain) derived by electron microscopy, Atmos. Environ., 41, 8058–8074,, 2007. 

Kang, S., Zhang, Y., Qian, Y., and Wang, H.: A review of black carbon in snow and ice and its impact on the cryosphere, Earth-Sci. Rev., 210, 103346,, 2020. 

Karanasiou, A., Diapouli, E., Cavalli, F., Eleftheriadis, K., Viana, M., Alastuey, A., Querol, X., and Reche, C.: On the quantification of atmospheric carbonate carbon by thermal/optical analysis protocols, Atmos. Meas. Tech., 4, 2409–2419,, 2011. 

Kaspari, S., Painter, T. H., Gysel, M., Skiles, S. M., and Schwikowski, M.: Seasonal and elevational variations of black carbon and dust in snow and ice in the Solu-Khumbu, Nepal and estimated radiative forcings, Atmos. Chem. Phys., 14, 8089–8103,, 2014. 

Kuchiki, K., Aoki, T., Niwano, M., Matoba, S., Kodama, Y., and Adachi, K.: Elemental carbon, organic carbon, and dust concentrations in snow measured with thermal optical and gravimetric methods: Variations during the 2007–2013 winters at Sapporo, Japan, J. Geophys. Res.-Atmos., 120, 868–882,, 2015. 

Li, C., Yan, F., Kang, S., Chen, P., Han, X., Hu, Z., Zhang, G., Hong, Y., Gao, S., Qu, B., Zhu, Z., Li, J., Chen, B., and Sillanpää, M.: Re-evaluating black carbon in the Himalayas and the Tibetan Plateau: concentrations and deposition, Atmos. Chem. Phys., 17, 11899–11912,, 2017. 

Lim, S., Faïn, X., Zanatta, M., Cozic, J., Jaffrezo, J.-L., Ginot, P., and Laj, P.: Refractory black carbon mass concentrations in snow and ice: method evaluation and inter-comparison with elemental carbon measurement, Atmos. Meas. Tech., 7, 3307–3324,, 2014. 

Liu, H., Chen, T., Zou, X., Qing, C., and Frost, R. L.: Thermal treatment of natural goethite: Thermal transformation and physical properties, Thermochim. Acta, 568, 115–121,, 2013. 

Meinander, O., Heikkinen, E., Aurela, M., and Hyvärinen, A.: Sampling, Filtering, and Analysis Protocols to Detect Black Carbon, Organic Carbon, and Total Carbon in Seasonal Surface Snow in an Urban Background and Arctic Finland (> 60 N), Atmosphere-Basel, 11, 923,, 2020. 

Neri, G., Bonaccorsi, L., Donato, A., Milone, C., Musolino, M. G., and Visco, A. M.: Catalytic combustion of diesel soot over metal oxide catalysts, Appl. Catal. B-Environ., 11, 217–231,, 1997. 

Nicolosi, E. M. G., Quincey, P., Font, A., and Fuller, G. W.: Light attenuation versus evolved carbon (AVEC) – A new way to look at elemental and organic carbon analysis, Atmos. Environ., 175, 145–153,, 2018. 

Petzold, A., Ogren, J. A., Fiebig, M., Laj, P., Li, S.-M., Baltensperger, U., Holzer-Popp, T., Kinne, S., Pappalardo, G., Sugimoto, N., Wehrli, C., Wiedensohler, A., and Zhang, X.-Y.: Recommendations for reporting “black carbon” measurements, Atmos. Chem. Phys., 13, 8365–8379,, 2013. 

Qian, Y., Yasunari, T. J., Doherty, S. J., Flanner, M. G., Lau, W. K., Ming, J., Wang, H., Wang, M., Warren, S. G., and Zhang, R.: Light-absorbing particles in snow and ice: Measurement and modeling of climatic and hydrological impact, Adv. Atmos. Sci., 32, 64–91,, 2015. 

Ramanathan, V. and Carmichael, G.: Global and regional climate changes due to black carbon, Nat. Geosci., 1, 221–227,, 2008. 

Schwarz, J. P., Doherty, S. J., Li, F., Ruggiero, S. T., Tanner, C. E., Perring, A. E., Gao, R. S., and Fahey, D. W.: Assessing Single Particle Soot Photometer and Integrating Sphere/Integrating Sandwich Spectrophotometer measurement techniques for quantifying black carbon concentration in snow, Atmos. Meas. Tech., 5, 2581–2592,, 2012. 

Svensson, J., Ström, J., Kivekäs, N., Dkhar, N. B., Tayal, S., Sharma, V. P., Jutila, A., Backman, J., Virkkula, A., Ruppel, M., Hyvärinen, A., Kontu, A., Hannula, H.-R., Leppäranta, M., Hooda, R. K., Korhola, A., Asmi, E., and Lihavainen, H.: Light-absorption of dust and elemental carbon in snow in the Indian Himalayas and the Finnish Arctic, Atmos. Meas. Tech., 11, 1403–1416,, 2018.  

Torres, A., Bond, T. C., Lehmann, C. M., Subramanian, R., and Hadley, O. L.: Measuring organic carbon and black carbon in rainwater: evaluation of methods, Aerosol Sci. Tech., 48, 239–250,, 2014. 

Tuzet, F., Dumont, M., Lafaysse, M., Picard, G., Arnaud, L., Voisin, D., Lejeune, Y., Charrois, L., Nabat, P., and Morin, S.: A multilayer physically based snowpack model simulating direct and indirect radiative impacts of light-absorbing impurities in snow, The Cryosphere, 11, 2633–2653,, 2017. 

Tuzet, F., Dumont, M., Picard, G., Lamare, M., Voisin, D., Nabat, P., Lafaysse, M., Larue, F., Revuelto, J., and Arnaud, L.: Quantification of the radiative impact of light-absorbing particles during two contrasted snow seasons at Col du Lautaret (2058 m a.s.l., French Alps), The Cryosphere, 14, 4553–4579,, 2020. 

US EPA: Microwave Assisted Acid Digestion of Siliceous and Organically Based Matrices, OHW, Method 3052, U.S. Environmental Protection Agency, (last access: 13 September 2022), 1996. 

Wang, M., Xu, B., Zhao, H., Cao, J., Joswiak, D., Wu, G., and Lin, S.: The influence of dust on quantitative measurements of black carbon in ice and snow when using a thermal optical method. Aerosol Sci. Tech., 46, 60–69,, 2012. 

Wang, X., Zhang, X., and Di, W.: Development of an improved two-sphere integration technique for quantifying black carbon concentrations in the atmosphere and seasonal snow, Atmos. Meas. Tech., 13, 39–52,, 2020. 

Warren, S. G. and Wiscombe, W. J.: A model for the spectral albedo of snow. II: Snow containing atmospheric aerosols, J. Atmos. Sci., 37, 2734–2745,<2734:AMFTSA>2.0.CO;2, 1980. 

Yamanoi, Y. and Nakashima, S.: In situ high-temperature visible microspectroscopy for volcanic materials, Appl. Spectrosc., 59, 1415–1419,, 2005. 

Yamanoi, Y., Nakashima, S., and Katsura, M.: Temperature dependence of reflectance spectra and color values of hematite by in situ, high-temperature visible micro-spectroscopy, Am. Mineral., 94, 90–97,, 2009. 

Short summary
The use of thermal–optical analysis for the determination of elemental carbon (EC) and organic carbon (OC) can be biased by mineral dust (MD). We present a method that utilizes this interference to quantify iron contained in MD in snow samples. Possibilities and limitations of applying this method to particulate matter samples are presented. The influence of light-absorbing iron compounds in MD on the transmittance signal can be used to identify samples experiencing a bias of the OC / EC split.