Abstract

Thermal desorption (TD) coupled with gas chromatography/mass spectrometry (TD-GC/MS) is a simple alternative that overcomes the main drawbacks of the solvent extraction-based method: long extraction times, high sample manipulation, and large amounts of solvent waste. This work describes the optimization of TD-GC/MS for the measurement of airborne polycyclic aromatic hydrocarbons (PAHs) in particulate phase. The performance of the method was tested by Standard Reference Material (SRM) 1649b urban dust and compared with the conventional method (Soxhlet extraction-GC/MS), showing a better recovery (mean of 97%), precision (mean of 12%), and accuracy (±25%) for the determination of 14 EPA PAHs. Furthermore, other 15 nonpriority PAHs were identified and quantified using their relative response factors (RRFs). Finally, the proposed method was successfully applied for the quantification of PAHs in real 8 h-samples (PM10), demonstrating its capability for determination of these compounds in short-term monitoring.

1. Introduction

Polycyclic aromatic hydrocarbons (PAHs) are a class of persistent organic pollutants (POPs) comprising hundreds of individual substances. These compounds contain two or more fused aromatic rings (made up of carbon and hydrogen atoms) in linear, angular, or cluster arrangements [1]. They are semivolatile organic compounds (SVOC); thus, they are present in the atmosphere in both the gas and the particulate phases as well as dissolved or suspended in precipitation (fog or rain) [2].

PAHs are as by-products of incomplete combustion processes of organic matter [3], and primarily emitted from anthropogenic sources [4], being the mobile and major sources in urban areas [5, 6]. Their harmful health effects and persistence pose an environmental concern. Thus, these compounds were among the first atmospheric pollutants identified as suspected carcinogens [7]. Moreover, PAHs belong to the group of POPs included in the list of 16 POPs specified by the UNECE Convention on Long-range Transboundary Air Pollution Protocol on Persistent Organic Pollutants [8, 9]. Due to these features, the United States Environmental Protection Agency (US-EPA) has listed 16 of them as priority pollutants (16 EPA PAHs) [10]. The most toxic PAHs (5 and 6 rings) are linked to the particulate matter [11, 12]. Accordingly, many air pollution studies have been focused on PAHs bound to particulate matter, particularly PM10 and PM2.5 in order to assess their concentration, distribution, and sources.

Air monitoring for PAHs in urban areas is an important issue because the risk associated with human exposure is higher considering the population density [13, 14]. However, PAH data in urban air show large spatial and temporal uncertainties because of the complex sampling and analytical procedures required.

Sampling of particulate PAHs is mostly done by the collection of them on a filter (quartz or glass fiber), using high- or low-volume samplers [1517]. Once the PAHs have been collected, they have to be extracted for the final determination. The extraction of PAHs from multiple matrices is a difficult step. PAHs are found in the environment in very low concentrations; consequently, an effective extraction method, able to quantitatively separate the analytes from the matrix, is required. The widely used method is solvent desorption of sampling media (Soxhlet extraction, accelerated solvent extraction, microwave-assisted extraction, and ultrasonic-assisted extraction) followed by analysis of the compounds of interest by GC-MS (gas chromatography coupled to mass chromatography) or high-performance liquid chromatography coupled with florescence detection (HPLC-FLD) [18, 19], where the detection methods allow cutting most of analytical interferences.

The use of toxic organic solvents in the solvent-based extraction methods causes added difficulties with sample handling and generates large amounts of solvent waste, which is costly and can generate additional environmental problems. Furthermore, the sensitivity of the current analytical procedures limits time resolution of measurements; thus, most of the urban pollution studies rarely achieved temporal resolution measurements better than 24 h. Since the PAH composition of aerosols can vary according to the diurnal changes in the sources, meteorological conditions and atmospheric reactivity [20], the time resolution of 24 h seems not sufficient to comprehend their variability, fate, and behavior in the atmosphere [21]. For these reasons, the development of simpler and sensitive methods or the improvement of the existing ones is of great interest, for the detection, determination, and monitoring of PAHs.

In recent years, alternative analytical procedures for PAHs based on the use of solvent-free extraction methods have been studied [22, 23]. Thermal desorption (TD) involves heating sample materials or sorbents in a flow of inert carrier gas, so that retained organic volatiles and semivolatiles are released and transferred or injected into the analytical system (e.g., into the carrier gas stream of the GC column).

The power and potential of TD allow configuring the technique in multiple adsorption-desorption stages, thus enhancing the concentration of the compounds of interest and detection limits. This higher sensitivity may provide shorter sampling times or lower sampling volumes. Another benefit of TD is that it is often possible to quantitatively retain target compounds during one or more of the trapping stages, while unwanted, for example, water and/or permanent gases, is selectively purged to vent. This avoids the entrance of unwanted compounds into the analytical system that could generate interferences during the analysis and/or damage to the equipment.

This work has tested and optimized different TD-GC/MS operation conditions in order to develop the best method that is able to sample and analyze airborne PAHs in particulate phase. The TD-GC/MS method was later validated by using a Standard Reference Material (SRM) 1649b urban dust and comparing with the conventional method based on solvent extraction (Soxhlet extraction-GC/MS). Moreover, the method was applied to measure PAH levels of 8 h PM10 samples in ambient air.

2. Materials and Methods

2.1. Reagents and Materials

A liquid certificated mixture of 16 EPA PAHs (2000 g·mL−1, SV Calibration Mix 5, Restek Corporation, USA) and a liquid deuterated mixture (200 g·mL−1, predeuterated internal standard PAH Mixture 6, Chiron AS, Norway) were used during the study. In Soxhlet extraction, decafluorobiphenyl, 4,4′-dibromooctafluorobiphenyl, 4,4′-dibromobiphenyl. (Restek, 2000 g·mL−1), and indeno [1,2,3-cd]pyrene-d12 (Chiron, 100 g·mL−1) were used as recovery standard for the assessment of extraction efficiency. Solutions were prepared by appropriate dilution in methanol, HPLC grade (99.9%, Lab-Scan Analytical Sciences, Poland).

The method was validated using the Standard Reference Material (SRM) 1649b urban dust, obtained from the National Institute of Standards and Technology (NIST, Gaithersburg, MD, USA).

2.2. TD-GC/MS Method

Sampling tubes (stainless steel tube of 5 mm outer diameter × 90 mm length) packed with filter were analyzed by using TD-GC/MS. The 16 EPA PAHs and deuterated PAHs were spiked on two one-eighth parts of a 47 mm quartz fiber filter (Whatman International Ltd., UK). These portions, suitably folded, were introduced into the sampling tubes.

Prior to use, the packed sampling tubes were conditioned by thermal cleaning under a helium flow rate of 100 mL·min−1 at 350°C for 30 min.

The NIST Standard Reference Material 1649b was handled in a similar way. Samples of the urban dust (10 mg) were weighed and placed on a one-eighth section of a 47 mm quartz fiber filter which was rolled and put into the sampling tube. Silanized glass wool (Supelco Inc., Bellefonte, USA) was introduced at the end and at the head of the desorption tubes in order to prevent system contamination.

Prior to use, filters and glass wool plugs were heated in a muffle furnace at 500°C for 24 h to remove trace organic compounds.

PAHs analysis was carried out using an automatic thermal desorber unit (Turbomatrix 150 ATD, Perkin Elmer S.L., USA) coupled by a fused silica capillary transfer line (5 m length × 0.32 mm I.D.) to a GC/MS detector (Clarus 500, Perkin Elmer S.L., USA). The chromatographic separation of PAHs was conducted on a Meta. X5 (silphenylene phase) capillary column: 30 m length × 0.25 mm I.D. × 0.25 mm film thickness (Teknokroma, Spain).

The helium gas carrier pressure employed in the GC/MS system was 145 kPa, and the column temperature was programmed as follows: initial temperature 100°C for 3 min, ramp of 10°C·min−1 until 250°C, ramp of 5°C·min−1 until 320°C, and finally temperature held at 320°C for 10 min. The total analysis time was 42 min per sample. The temperature of both the transfer lines (from TD to GC and from GC to MS) was held at 280°C, whereas the source temperature was 250°C. Simultaneous full scan (SCAN) and selective ion monitoring (SIM) modes were used for the identification and quantification of PAHs. Table 1 shows, according to their elution order, the PAHs determined in this study with their quantification ions. Supplementry Figures S1S3 show the representative SCAN chromatograms of the 16 EPA and deuterated PAHs.

2.3. Soxhlet Extraction-GC/MS Method

Between 300 and 500 mg of the NIST SRM 1649b urban dust was weighted and placed on one-eighth of a 150 mm prebaked (at 500°C for 24 h) quartz fiber filter (Whatman International Ltd.). Before folding the filter, it was spiked with the recovery standards.

Soxhlet extraction was performed by using Büchi extraction system B-811 (BÜCHI, Switzerland), an automated system that can be used to perform extraction according to the original Soxhlet principle. The samples were extracted with hexane using the Soxhlet Warm mode. This mode increases the solubility of the analytes, allowing an optimal extraction in 3 hours [24].

After the extraction process, the extracts of 5 mL were concentrated by a stream of dry nitrogen to a volume less than 0.5 mL. Finally, these extracts were diluted to 1.5 mL with methanol and spiked with deuterated PAHs.

Two-microliters  of aliquots from each extract was injected into the GC/MS with split mode. Table 2 collects the timed events and the oven program used in the GC/MS during the validation of the Soxhlet method.

Also, in this method, the GC/MS used simultaneously the SCAN and SIM mode for the identification and quantification of PAHs.

2.4. Sample Collection

Airborne particulate matter (PM10) samples were collected on preheated (at 500°C for 24 h) quartz fiber filters (150 mm diameter, Whatman International Ltd., United Kingdom) using a high-volume sampler (Digitel DHA-80, Digitel Elektronik AG, Switzerland) with a flow rate of 30 m3·h−1. DHA-80 stores 15 filters stretched in filter holders that are changed automatically at the preset time. DHA-80 has integrated temperature control in the filter storage section; in this way, the used filters can be stored at low temperatures (in this study at 4°C) after sampling.

Collected filters were put into individual Petri dishes, wrapped in aluminum foil, and kept in a 4°C freezer until analysis (<15 days) according to ISO 12884:2000 [25].

3. Results and Discussion

3.1. Optimization of Thermal Desorption Method

The thermal desorption process can be divided into two main stages: tube and trap desorption. In the first stage, target compounds are thermally desorbed from the sampling tube and transferred to the cold trap, where they are concentrated. After completing the primary desorption, the trapped compounds are released by quick heating of the trap and swept through the heated transfer line to the GC column.

To obtain the best analytical conditions in terms of sensitivity and reproducibility, different parameters in each desorption stage were tested.

For these tests, 1 L of the 16 EPA PAHs solutions of 20 ng·L−1 were spiked in sampling tubes packed with portions of quartz fiber.

3.1.1. Primary Desorption (Tube Desorption)

The conditions in the tube oven during this stage are key to guarantee an efficient desorption; thus, parameters such as the temperature, time, and flow in the tube oven were studied to optimize this process.

Different values of desorption temperatures, times, and flows were tested, considering factors such as the packing/sample matrix stability, the lability of the components of interest, and the temperature limitations of the system. Figure 1 shows the area of chromatographic peak for each of the 16 EPA PAHs (in %) obtained for each test.

The results demonstrated that an increase in the temperature of the oven tube enhances the desorption of particulate PAHs (Figure 1(a)). This improvement was remarkable for high molecular weight PAHs (IP, DBahA, and BghiP). A value of 320°C was selected as temperature in the first desorption stage. Regarding the time, the lowest value in the test (10 min) clearly showed significantly higher areas for most compounds (Figure 1(b)), indicating a more efficient desorption. This value was selected as desorption time in the tube. Finally, higher desorption flows enable better desorption of PAHs (Figure 1(c)). Flows higher than 150 mL·min−1 are not recommended, as they can generate problems in maintaining low temperatures in the trap zone during the first desorption stage [26]. Therefore, a flow of 150 mL·min−1 was selected as the optimal value.

3.1.2. Secondary Desorption (Trap Desorption)

To enhance PAH desorption from the trap (a quartz tube packed with glass wool), its high temperature has to be as high as possible. This temperature depends on the trap packing and equipment stability as well as on the target compounds. In this study, a value of 320°C (the value recommended by the manufacturer) was set, while its low temperature (values of −15, −10, and −5°C) and desorption time (values of 4, 6, and 12 min) were tested.

The area of chromatographic peak (in %) for low (2-3 rings: Naph, Ace, Acy, FL, Phe, and Ant), middle (4 rings: Ft, Pyr, BaA, and Chry), and high (5-6 rings: BbFt, BkFt, BaP, IP, DBahA, and BghiP) molecular weight PAHs obtained for each test are shown in Figure 2.

The temperature in the Peltier trap is a critical parameter in secondary desorption (Figure 2(a)), showing significant changes in the sensitivity for different values. The temperature of −10°C revealed the best results.

In the study of the trap desorption time (Figure 2(b)), the results demonstrated that longer values do not implicate a higher efficiency and consequently a better detection, 6 min showed a better response than 12. This is especially significant with the lightest PAHs (LMW) which could be affected by the exposure to high temperatures, generating losses. By contrast, the heavier PAHs (MMW and HMW) showed higher concentrations after longer trap desorption times because they could need more time to be completely desorbed. Due to this, a trap duration of 6 min was selected as this value presented good desorption for 16 target PAHs.

3.1.3. Inlet and Outlet Split Flows

In order to enhance the process of two-stage thermal desorption, a double split mode was used. Therefore, the inlet (split flow as the tube is desorbed) and outlet (split flow as the trap is desorbed) split flows were also tested.

The inlet split flow plays an important role during primary desorption. This maintains a relatively high carrier gas flow through the sample tube, while at the same time establishes a reasonably low flow through the cold trap, aiding the complete removal from the sample tube and analyte retention. The deactivation of the inlet split (0 mL·min−1) generated a significant improvement in PAH desorption because the complete sample, without purge, arrived at the cold trap. With the increase of inlet split, the sensitivity decreased. Although in this study, an inlet split flow of 0 mL·min−1 showed the best results; it is recommended to activate this split in order to avoid the presence of unwanted compounds (permanent gases and water) in the trap. These could reduce the trap lifetime and interfere in the analysis. In order to find a compromise between these rules, an intermediate flow (23 mL·min−1) was considered as the optimal value.

The outlet split flow also plays an important function in the trap desorption: (1) adapting the effluent flow to a capillary column flow, it avoids the system saturation and (2) facilitating the release of the analytes, it guarantees a high enough flow through the trap during desorption. According to the manufacturer, at least 10 mL·min−1 of outlet split is necessary to minimize the air/water background on a mass spectrometer when atmospheric samples are analysed [26]. In this study, the results obtained for outlet split flows demonstrated that the increase of this parameter reduces the sensitivity of the technique, with losses becoming significant between 10 and 20 mL·min−1. Therefore, the manufacturer’s recommended flow (10 mL·min−1) was selected as the optimal value.

Finally, Table 3 summaries the optimized values for thermal desorption.

3.2. Desorption Efficiency

Once optimized, the efficiency of two-stage thermal desorption was studied. The efficiency was calculated by the following expression:where is the efficiency in , is the peak area of the analyte obtained from desorption of the sampling tube (previously loaded with PAHs), and is the peak area of the analyte obtained when the same sampling tube or the trap was desorbed the second time.

Figure 3 shows the tube and trap efficiencies obtained for each particulate PAH. The technique demonstrated a good efficiency with recoveries of the PAHs in the tube and trap higher than 94%.

3.3. TD-GC/MS Validation and Comparison with Soxhlet-GC/MS

In order to determine the performance of the method when applied to atmospheric PM samples, this was tested using the Standard Reference Material (SRM) 1649b urban dust.

The accuracy, repeatability, and recovery of the method were calculated by adding known amounts (approximately 10 mg) of the SRM 1649b to a one-eighth section of blank filters (). Before the analysis, filters were spiked with 1 L of the deuterated PAH internal standard solution (25 ng·L−1).

Table 4 shows the results obtained for each PAH by using the TD-GC/MS method, comparing the calculated concentration with the certified values.

Although the column used in this study demonstrated a good resolution for the 16 EPA PAHs, the presence of other PAHs in the urban dust can generate coelution problems with the target compounds [27]. Some PAH pairs such as BbFt and DBahA coeluted with the benzo[j]fluoranthene (BjFt) and dibenzo[a,c]anthracene (DBacA), respectively.

The TD-GC/MS method showed good precision with mean RSD values of 12.2. The accuracy of the TD-GC/MS method ranged from −22.8% to 25.1%, while the average recovery efficiency was 96.7. These performance parameters of the TD-GC/MS method accomplish the quality objectives for ambient air PAHs stated by ISO 12884:2000 [25], which establishes a precision of ±25%, an accuracy of ±20%, and a recovery efficiency between 75 and 125%. These requirements are accomplished for most PAH; however, there are some exceptions. The lowest molecular weight PAHs (Naph and Acy), with excessively high recoveries, confirmed the overestimation of these compounds when analyzed by using the TD-GC/MS method. These compounds could suffer losses during the sample preparation due to their high volatility. Besides, the presence of interfering compounds in the SRM and the low concentration of Acy could explain these overestimations. Therefore, this method was not applicable to the Naph and Acy analysis.

In order to demonstrate the efficiency of the TD method as compared with other analytical methods, the conventional method (Soxhlet extraction-GC/MS) was also tested (Table 5). Between 300 and 500 mg of the NIST SRM 1649b urban dust was placed on one-eighth of a 150 mm prebaked quartz fiber filter, which was spiked with 1 L of a solution (0.5 ng·L−1) of the recovery standards. After the extraction process, the obtained extracts were spiked with 25 L of a deuterated PAH solution (20 ng·L−1).

The results of the Soxhlet extraction-GC/MS method showed a good recovery for 4-, 5-, and 6-ring PAHs with values between 72.8 and 131%; whereas the lightest PAHs (2- and 3-ring PAHs), except Phe, showed low recovery (<70%). The loss of these analytes during the extraction process in the Soxhlet Warm mode could be the main reason for these low recoveries. In the case of DBahA, although its coelution with DBacA was considered, its recovery continued to be high (>200%). This indicates an overestimation in the determination of this compound by the Soxhlet process. Regarding precision and accuracy, the Soxhlet extraction-GC/MS showed worse results, with an average precision of 34.9 and values of accuracy out of the limits ±20% for some PAHs.

Comparing both methods, the TD-CG/MS method demonstrated a better performance (good recovery, precision, and accuracy) for the determination of PAHs (except for Naph and Acy). By contrast, the manipulation of the samples in the Soxhlet process meant losses of the light PAHs (2- and 3-ring) and an overestimation of some PAHs, especially of the DBahA.

Regarding the detection limits, the average instrument detection limit (IDL) of the TD-GC/MS method was 0.04 ng and the average method detection limit (MDL), assuming a total sample volume 240 m3 (30 m3·h−1 for 8 h), was 2.89 × 10−3 ng·m−3 [28].

3.4. Extension of the Scope to Other PAHs

Although most environmental studies are focused on the analysis of 16 PAH listed by US-EPA, it could be interesting to determine other PAHs in order to have a better characterization of these compounds in terms of toxicity and sources. For this reason, besides the 16 EPA PAHs, other 15 PAHs were determinated by using TD-GC/MS. Table 1 shows, according to their elution order, the 16 EPA PAHs, the deuterated PAHs, and the 15 nonpriority PAHs.

Because SRM 1649b contains other compounds besides the 16 EPA PAHs, it was used to identify 15 nonpriority PAHs and to quantify them by relative response factors (RRFs). Supplementry Figures S4S7 show the SIM chromatograms of the target PAHs (16 EPA PAHs + 15 PAHs) obtained in the analysis of SRM 1649b. For quantification, the RRFs for each nonpriority PAH were calculated by the following equation:where is the peak area of nonpriority PAH, is the peak area of reference PAH compound, CPAH is the nonpriority PAH concentration in the NIST dust, and is the reference PAH concentration in the NIST dust.

Reference PAHs were selected according to the following criterion: the nearest of the 16 EPA PAHs to each new one, which provides the least variation in the RRF. Table 6 collects the reference PAHs, RRFs, and the relative standard deviations (RSDs) for each nonpriority PAH.

The nonpriority PAHs showed a range of RRFs between 0.31 and 4.74, with RSD of less than 15% for most compounds. In the case of Ret, the low chemical similarity between this compound and its reference PAH (BaA) could explain the poor precision in the determination of its RRF (>20%).

3.5. Application to Real Samples

After validation, the method described in this study was applied to extract and analyze samples collected in the city of Bilbao, Spain (longitude 2°56′56.24″W, latitude 3°15′44.86″N). Bilbao city is the most populated area in the Basque Country and the tenth largest in Spain (approximately 350,000 in the city and 1 million inhabitants in the metropolitan area). In this urban area, local traffic and stationary emissions from the surrounding industries are considered as the major sources of atmospheric pollutants [28].

During seven consecutive days per month, eight-hour PM10 samples were collected at a flow rate of 30 m3·h−1.

A total of 182 PM10 samples were collected over 9 months (between July 2013 and June 2014). Each sample was randomly cut into 8 portions of 1 cm2 and introduced into the sampling tube and analyzed using the optimized method. This was performed in the same way as other studies [29, 30], which demonstrated good homogeneity results when using small sections of the filters.

Table 7 shows the descriptive statistics (number of valid data, mean, standard deviation, minimum, maximum, 5th, and 95th percentiles) for individual PAH concentrations measured by using the TD-GC/MS method in the city of Bilbao (urban area). PAHs which showed overestimation in the SRM analysis (Naph and Acy) or poor precision in the RRF determination (Ret) were not measured in the real samples.

The average concentration of individual EPA PAHs in Bilbao ranged from 0.04 ± 0.05 to 0.50 ± 0.76 ng·m−3, whereas the nonpriority PAHs were between 0.01 ± 0.01 and 0.26 ± 0.32 ng·m−3. The EPA PAHs reported minimum values between and 0.03 ng·m−3 for most of the compounds, which are between 1.1 and 14.6 times the MDL, showing the suitability of the proposed method to determine particle-bound PAHs in real PM10 samples. Although the minimums of IP and DBahA were below their MDL, these values meant only the 5% or less of the measured samples.

Among compounds, BbFt was the major contributor to total PAHs (average concentration of 0.5 ± 0.76 ng·m−3), followed by Pyr (0.27 ± 0.23 ng·m−3), Ft (0.26 ± 0.22 ng·m−3), BeP (0.26 ± 0.32 ng·m−3), and Chry (0.22 ± 0.24 ng·m−3). The high presence of these compounds in PM10 fraction has been reported by previous studies [31, 32] in urban areas with traffic loads.

4. Conclusions

The method developed in this study, based on thermal desorption, showed a good efficiency for the determination of particle-bound PAHs. The use of a solvent-free extraction technique has showed numerous advantages (less sample manipulation and analysis time, reduced exposure risk, and higher sensitivity and reliability) that enable a better performance (good recovery, precision, and accuracy) for the determination of particle-bound PAHs; however, the lowest molecular weight PAHs (Naph and Acy) could be overestimated by this methodology.

Parameters such as tube and trap temperature, time, desorption, and split flows (inlet and outlet) were critical in the thermal desorption of PAHs. The optimized TD-GC/MS method showed an efficient desorption of PAHs with recoveries higher than 94%.

The results obtained in the validation of TD-GC/MS by standard reference material (urban dust) demonstrated that this is a reliable method to determine particulate PAHs in aerosol samples (good precision and accuracy), with average recovery efficiency of 96.67 and a mean RSD value of 12.18. Comparing with the conventional method Soxhlet-GC/MS, the TD-CG/MS method demonstrated a better performance for the determination of PAHs. Besides 16 EPA PAHs, the TD-GC/MS method demonstrated its ability to quantify other PAHs in aerosol samples.

Finally, the method was successfully applied for the quantification of PAHs in real PM10 samples collected with a time resolution of 8 h.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors gratefully thank the University of the Basque Country UPV/EHU (Ref.: GIU 13/25, GIU 16/03, and UFI 11/47) and the Spanish Ministry of Science and Innovation (MICINN) for financing the project PROMESHAP (Ref.: CTM 2010-20607). Iñaki Elorduy wants to thank the MICINN for his doctoral grant Ministerio de Ciencia e Innovació.

Supplementary Materials

The representative SCAN chromatograms of the 16 EPAs and deuterated PAHs and also the SIM chromatograms of the target PAHs (16 EPA PAH + 15 PAH) obtained in the SRM 1649b analysis are shown in the supplementary material. Figure S1: chromatogram of 16 EPA and deuterated PAHs in SCAN mode, from 0 to 15 min: (1) Naph-d8, (2) Naph, (3) Bph-d10, (4) Acy, (5) Ace, and (6) FL. Figure S2: chromatogram of 16 EPA and deuterated PAHs in SCAN mode, from 15 to 20.5 min: (7) Phe-d10, (8) Phe, (9) Ant, (10) Ft, (11) Pyr-d10, and (12) Pyr. Figure S3: chromatogram of 16 EPA and deuterated PAHs in SCAN mode, from 22.5 to 35.5 min: (13) BaA-d12, (14) BaA, (15) Chry, (16) BbFt, (17) BkFt, (18) BaP-d10, (19) BaP, (20) IP, (21) DBahA, (22) BghiP-d12, and (23) BghiP. Figure S4: PAHs and deuterated PAHs in SIM windows (m/z 226, 240, 228, and 234) in the analysis of NIST SRM 1649b dust. Figure S5: PAHs and deuterated PAHs in SIM windows (m/z 252 and 264) in the analysis of NIST SRM 1649b dust. Figure S6: PAHs and deuterated PAHs in SIM windows (m/z 276, 288, and 278) in the analysis of NIST SRM 1649b dust. Figure S7: PAHs and deuterated PAHs in m/z 300 SIM window in the analysis of NIST SRM 1649b dust. (Supplementary Materials)