# Determination of Selected Polycyclic Aromatic Compounds in Particulate Matter Samples with Low Mass Loading: An Approach to Test Method Accuracy

Authors: Susana García-Alonso; Rosa María Pérez-Pastor; Víctor Archilla-Prat; David Sanz-Rivera; Enrique Rojas-García; Gaizka Aragón-Urrengoetxea; Jesús Javier Rodríguez-Maroto
Journal: Journal of Chemistry
 (2017)

---

## Abstract

A miniaturized analytical procedure to determine selected polycyclic aromatic compounds (PACs) in low mass loadings (<10 mg) of particulate matter (PM) is evaluated. The proposed method is based on a simple sonication/agitation method using small amounts of solvent for extraction. The use of a reduced sample size of particulate matter is often limiting for allowing the quantification of analytes. This also leads to the need for changing analytical procedures and evaluating its performance. The trueness and precision of the proposed method were tested using ambient air samples. Analytical results from the proposed method were compared with those of pressurized liquid and microwave extractions. Selected PACs (polycyclic aromatic hydrocarbons (PAHs) and nitro polycyclic aromatic hydrocarbons (NPAHs)) were determined by liquid chromatography with fluorescence detection (HPLC/FD). Taking results from pressurized liquid extractions as reference values, recovery rates of sonication/agitation method were over 80% for the most abundant PAHs. Recovery rates of selected NPAHs were lower. Enhanced rates were obtained when methanol was used as a modifier. Intermediate precision was estimated by data comparison from two mathematical approaches: normalized difference data and pooled relative deviations. Intermediate precision was in the range of 10–20%. The effectiveness of the proposed method was evaluated in PM aerosol samples collected with very low mass loadings (<0.2 mg) during characterization studies from turbofan engine exhausts.


---

## Body

## 1. Introduction

There is an increasing concern about the chemical characterization of environmental samples of size-segregated samples [1–3]. This leads to the analysis of reduced sample sizes of particulate matter (PM) containing very low levels of pollutants. These samples require improving both analytical method extraction and quality. Among chemicals typically linked to PM, PACs, a well-known family of toxic contaminants, are of great concern. Conventional analytical procedures usually require samples with mass loadings in order of dozens of mg (i.e., ambient air samples collected by using high volume sampler during 24 hours). Among extraction techniques, soxhlet (SX), microwave assisted extraction (MAE), pressurized liquid extraction (PLE), and ultrasonic extractions (US) have been extensively used [4, 5]. Large amount of solvent (SX), high cost (PLE and MAE), and low extraction efficiency (US) are important limitations when large sets of samples with reduced mass loadings are considered for analysis.Extraction of PM samples with small mass loadings has been carried out using PLE [1, 2] and thermal desorption [6, 7]. These extractions are expensive in instrument and maintenance. Due to lower equipment cost and ease of operation, extraction by US stands as an attractive alternative [3, 8]. Hence, procedures based on reducing sample size and solvent requirements have been developed to achieve a miniaturization strategy for US extraction [9–11]. When the analyses of samples with low mass loadings are discussed, the validation of analytical method is an essential issue. Not enough attention is always paid to this matter.A validation study includes the estimation of the accuracy of analytical method, which reflects the closeness of a measured value to a true value. Therefore, two components of accuracy must be assessed: trueness and precision [12]. If intralaboratory tests have to be carried out, the performance of proposed method will be based on the use of recovery assessments, reference materials, or comparison with a reference method [13]. The comparison of results obtained by proposed method with those found with a comparison method of known accuracy is a useful tool for this kind of studies. Perhaps a conclusion of quantitative extraction after the analysis of certified reference materials could not be extended at real samples. Possible matrix effects coming from the sample material could be easily delimiting the quality of final results.As a part of aerosol characterization studies in gas turbine exhaust [12], PAC determination in PM samples with low mass loadings was required. In the present study, an affordable procedure based on an extraction by sonication and agitation with small volumes of solvent was evaluated. Quality of results was pursued by means of the evaluation of the extraction efficiency and intermediate precision from the analyses of PM samples. Keeping simplicity and low cost as advantages, our main goal was to evaluate analytical procedure for determination of selected PACs using real samples. This is rarely applied. Moreover, this work demonstrates the application of the proposed method.The applicability of the method was tested at the most disadvantageous case: PM samples with very low mass loadings (in order of micrograms) and PAC levels close to analytical detection limits (in order up to 1 ng per entire filter). These extremely low levels of PM and analyte require an assessment of the reliability of analytical measurements.


## 2. Materials and Methods

### 2.1. Reagents and Standards

Acetic acid (for analytical grade), dichloromethane (for organic trace analysis), and acetonitrile (HPLC grade) from SDS, Carlo Erba (Val de Reuil, France), were used for sample preparation and chromatographic analysis. Sodium borohydride (99%) and copper(II) chloride (97%) from Sigma Aldrich (Steinheim, Germany) were used for derivatization reaction to determine NPAHs.The PAH Mix 9 calibration check solutions (10µg mL−1 in acetonitrile) were supplied by Dr. Ehrenstorfer GmbH (Augsburg, Germany), containing naphthalene (Na), acenaphthene (Ace), fluorene (F), phenanthrene (Ph), anthracene (An), fluoranthene (Fl), pyrene (Pyr), benz[a]anthracene (BaA), chrysene (Cry), benz[b]fluoranthene (BbF), benz[k]fluoranthene (BkF), benz[a]pyrene (BaP), dibenz[a,h]anthracene (DBA), and benzo[g,h,i]perylene (BghiP). 1-Methyl-naphthalene (1-MeNa), 98.6%, and 2-methyl-naphthalene (2-MeNa), 97.2%, were purchased from Ultra Scientific (North Kingstown, USA).Individual nitro-PAH (NPAHs) calibration check solutions (methanol, 10µg/mL) were supplied by Dr. Ehrenstorfer Reference Materials (Augsburg, Germany). Selected NPAHs were 2-nitro-naphthalene (2-NNa), 3-nitro-phenanthrene (3-NPhe), 9-nitro-phenanthrene (9-NPhe), 1-nitro-pyrene (1-NPyr), and 3-nitro-fluoranthene (3-NFl). Stock solutions prepared in methanol were used to identify and quantify the aromatic compounds.The SRM 1649b, an urban dust, was obtained from the National Institute of Standards and Technology (NIST, USA).


### 2.2. Samples and Preparation

For the evaluation of the extraction efficiency and intermediate precision of results from the proposed method, twelve ambient air PM10 samples were taken and analyzed. Sampling site was the National Institute for Aerospace Technology (INTA), located at 20 km Northeast of Madrid city. A high volume MCV CAV-Amb sampler loaded with Whatman glass fiber filters (102 mm diameter) was used. Filter samples were wrapped separately in aluminum foil and stored in a freezer at −20°C until analysis, within 1 week from collection.The ambient air samples were taken with mass loadings of particulate matter between 3 and 40 mg as follows. Sampling was carried out during 24 hours (30 m3 h−1) for consecutive days between April and November of 2014. Only samples 10, 11, and 12 were collected during 30 min, 50 min, and 2 hours to achieve lower mass loadings of particulate matter (3 mg, 5 mg, and 11 mg, resp.). All filters were accurately weighed before and after sampling at constant humidity to a precision of 10 µg prior to sampling.Each sample was divided into four subsamples according to the flow diagram presented in Figures1s and 2s of the Electronic Supplementary Material (ESM), available online at https://doi.org/10.1155/2017/3517457. Hence, four subsamples of each PM10 sample were weighed, treated, and analyzed separately to evaluate four extraction procedures based on the following:(i)
MAE: PM subsamples were extracted using an Ethos Sel (Millestone) at 900 w (120°C), during 40 min. Dichloromethane was used as extraction solvent (15 mL).(ii)
PLE:PM subsamples were extracted using an ASE350 (Dionex Corporation) and dichloromethane (DCM) as extraction solvent (100°C, 10.34 psi and 1 static cycle). The total volume of extract which was obtained after extraction corresponded to 55 mL.(iii)
AG-D and AG-DM: PM subsamples were extracted in closed tubes containing 5 mL of extraction solvent (dichloromethane and dichloromethane/methanol 2 : 1, resp.). The extraction was performed in an ultrasonic bath for 15 min, using a Selecta Ultrasons-H (Barcelona, Spain). Subsequently the tubes were agitated for short time using an automatic shaker with a touch function (IKA vortex 1, Staufen, Germany). After the main extraction, additional extractions were carried out by adding 1 mL of extraction solvent (twice) and agitating for a short time. All extracts were then combined.The final extracts were filtered (0.22µm pore size PTFE) and slowly concentrated in a drying step by stream of nitrogen. Finally, dried samples were redissolved by adding 150–250 µL of methanol for PAH determination.Once selected PAHs were quantified; derivatization reaction of extracts was performed for NPAH determination [14]. The procedure was based on addition of 0.5 mL of 0.05% aqueous solution of CuCl2 and 0.05 g of NaBH4 to 0.5 mL of extract in methanol. This mixture was maintained at room temperature for 1 h. After this, 0.5 mL of water and 1 mL of dichloromethane were added and then agitated for extraction of amino-PAHs. Dichloromethane extract was collected with a syringe and two reextractions were performed to ensure recovery. Combined extracts were concentrated under nitrogen flow to change solvent (acetonitrile, 100 µL).The applicability of the method was then tested from the analyses of PM aerosols at very low mass loadings. Involved samples corresponded to the following:(i)
Mass loadings of 0.2 mg collected on quartz fiber filter (47 mm) by means of a Low Pressure Impactor(ii)
Mass loadings of 0.02 mg collected on aluminum (Al) foils (47 mm) by means of a holder filterThese samples were treated by AG-D method and analyzed as previously indicated.


### 2.3. Chromatographic Analysis

Sample analyses were performed on an Agilent series 1200 liquid chromatograph coupled to an Agilent 1100 fluorescent detector (Waldbronn, Germany). Particular conditions were previously optimized and well established [14]. Briefly, a supelcosil C18 column (Supelco, 250∗4.6 mm, 5 µm) and an eclipse XDB C18 column (Agilent, 150∗4.6 mm, 5 µm) were used for PAH and NPAH determinations, respectively. The mobile phase rates were 1.5 mL and 1.0 mL for PAH and NPAH determinations, respectively. Operational conditions of HPLC/FD are compiled in Table 1s (ESM).Quantification was performed by external standards with five-point calibration. Linear response was established from 10 ng mL−1 to 1000 ng mL−1 for PAHs. Due to the lower detector sensitivity for NPAHs, its linear concentration range was from 100 ng mL−1 to 1000 ng mL−1. Regression coefficients were above 0.999. Analyte concentration should be in the middle of the calibration range of standards, but a response factor based on a lower/higher point of standard solutions was applied if this was not possible [15].


## 2.1. Reagents and Standards

Acetic acid (for analytical grade), dichloromethane (for organic trace analysis), and acetonitrile (HPLC grade) from SDS, Carlo Erba (Val de Reuil, France), were used for sample preparation and chromatographic analysis. Sodium borohydride (99%) and copper(II) chloride (97%) from Sigma Aldrich (Steinheim, Germany) were used for derivatization reaction to determine NPAHs.The PAH Mix 9 calibration check solutions (10µg mL−1 in acetonitrile) were supplied by Dr. Ehrenstorfer GmbH (Augsburg, Germany), containing naphthalene (Na), acenaphthene (Ace), fluorene (F), phenanthrene (Ph), anthracene (An), fluoranthene (Fl), pyrene (Pyr), benz[a]anthracene (BaA), chrysene (Cry), benz[b]fluoranthene (BbF), benz[k]fluoranthene (BkF), benz[a]pyrene (BaP), dibenz[a,h]anthracene (DBA), and benzo[g,h,i]perylene (BghiP). 1-Methyl-naphthalene (1-MeNa), 98.6%, and 2-methyl-naphthalene (2-MeNa), 97.2%, were purchased from Ultra Scientific (North Kingstown, USA).Individual nitro-PAH (NPAHs) calibration check solutions (methanol, 10µg/mL) were supplied by Dr. Ehrenstorfer Reference Materials (Augsburg, Germany). Selected NPAHs were 2-nitro-naphthalene (2-NNa), 3-nitro-phenanthrene (3-NPhe), 9-nitro-phenanthrene (9-NPhe), 1-nitro-pyrene (1-NPyr), and 3-nitro-fluoranthene (3-NFl). Stock solutions prepared in methanol were used to identify and quantify the aromatic compounds.The SRM 1649b, an urban dust, was obtained from the National Institute of Standards and Technology (NIST, USA).


## 2.2. Samples and Preparation

For the evaluation of the extraction efficiency and intermediate precision of results from the proposed method, twelve ambient air PM10 samples were taken and analyzed. Sampling site was the National Institute for Aerospace Technology (INTA), located at 20 km Northeast of Madrid city. A high volume MCV CAV-Amb sampler loaded with Whatman glass fiber filters (102 mm diameter) was used. Filter samples were wrapped separately in aluminum foil and stored in a freezer at −20°C until analysis, within 1 week from collection.The ambient air samples were taken with mass loadings of particulate matter between 3 and 40 mg as follows. Sampling was carried out during 24 hours (30 m3 h−1) for consecutive days between April and November of 2014. Only samples 10, 11, and 12 were collected during 30 min, 50 min, and 2 hours to achieve lower mass loadings of particulate matter (3 mg, 5 mg, and 11 mg, resp.). All filters were accurately weighed before and after sampling at constant humidity to a precision of 10 µg prior to sampling.Each sample was divided into four subsamples according to the flow diagram presented in Figures1s and 2s of the Electronic Supplementary Material (ESM), available online at https://doi.org/10.1155/2017/3517457. Hence, four subsamples of each PM10 sample were weighed, treated, and analyzed separately to evaluate four extraction procedures based on the following:(i)
MAE: PM subsamples were extracted using an Ethos Sel (Millestone) at 900 w (120°C), during 40 min. Dichloromethane was used as extraction solvent (15 mL).(ii)
PLE:PM subsamples were extracted using an ASE350 (Dionex Corporation) and dichloromethane (DCM) as extraction solvent (100°C, 10.34 psi and 1 static cycle). The total volume of extract which was obtained after extraction corresponded to 55 mL.(iii)
AG-D and AG-DM: PM subsamples were extracted in closed tubes containing 5 mL of extraction solvent (dichloromethane and dichloromethane/methanol 2 : 1, resp.). The extraction was performed in an ultrasonic bath for 15 min, using a Selecta Ultrasons-H (Barcelona, Spain). Subsequently the tubes were agitated for short time using an automatic shaker with a touch function (IKA vortex 1, Staufen, Germany). After the main extraction, additional extractions were carried out by adding 1 mL of extraction solvent (twice) and agitating for a short time. All extracts were then combined.The final extracts were filtered (0.22µm pore size PTFE) and slowly concentrated in a drying step by stream of nitrogen. Finally, dried samples were redissolved by adding 150–250 µL of methanol for PAH determination.Once selected PAHs were quantified; derivatization reaction of extracts was performed for NPAH determination [14]. The procedure was based on addition of 0.5 mL of 0.05% aqueous solution of CuCl2 and 0.05 g of NaBH4 to 0.5 mL of extract in methanol. This mixture was maintained at room temperature for 1 h. After this, 0.5 mL of water and 1 mL of dichloromethane were added and then agitated for extraction of amino-PAHs. Dichloromethane extract was collected with a syringe and two reextractions were performed to ensure recovery. Combined extracts were concentrated under nitrogen flow to change solvent (acetonitrile, 100 µL).The applicability of the method was then tested from the analyses of PM aerosols at very low mass loadings. Involved samples corresponded to the following:(i)
Mass loadings of 0.2 mg collected on quartz fiber filter (47 mm) by means of a Low Pressure Impactor(ii)
Mass loadings of 0.02 mg collected on aluminum (Al) foils (47 mm) by means of a holder filterThese samples were treated by AG-D method and analyzed as previously indicated.


## 2.3. Chromatographic Analysis

Sample analyses were performed on an Agilent series 1200 liquid chromatograph coupled to an Agilent 1100 fluorescent detector (Waldbronn, Germany). Particular conditions were previously optimized and well established [14]. Briefly, a supelcosil C18 column (Supelco, 250∗4.6 mm, 5 µm) and an eclipse XDB C18 column (Agilent, 150∗4.6 mm, 5 µm) were used for PAH and NPAH determinations, respectively. The mobile phase rates were 1.5 mL and 1.0 mL for PAH and NPAH determinations, respectively. Operational conditions of HPLC/FD are compiled in Table 1s (ESM).Quantification was performed by external standards with five-point calibration. Linear response was established from 10 ng mL−1 to 1000 ng mL−1 for PAHs. Due to the lower detector sensitivity for NPAHs, its linear concentration range was from 100 ng mL−1 to 1000 ng mL−1. Regression coefficients were above 0.999. Analyte concentration should be in the middle of the calibration range of standards, but a response factor based on a lower/higher point of standard solutions was applied if this was not possible [15].


## 3. Results and Discussion

### 3.1. Quality Control

Instrumental limit of detection of each PAC was calculated as twice the standard deviations deduced from analyses of three standard solutions containing 1 ng mL−1 and 10 ng mL−1 for selected PAHs and NPAH, respectively. Derived values were between 0.5 ng mL−1 and 5 ng mL−1 for selected PAHs and NPAHs, respectively.The analysis of blanks (i.e., chemical analysis excluding PM samples) indicated the potential cross contamination for the more volatile PAHs if extreme care is not taken. Due to ubiquitous character of PAH compounds in the air, laboratory inventory, and low volumes of extracts, there is a risk of contamination. Therefore, to determine blank contents for each set of sample analyses is recommended. The very low levels require a rigorous protocol to ensure minimal sample contamination from vials, syringes, solvents, and drying step.Blank filters of glass fiber, quartz, and aluminum were handled in the same way as the samples by agitation procedure and these trials concluded to detectable contents only for the more volatile PAHs. In detail, 1-methyl-naphthalene, 2-methyl-naphthalene, and phenanthrene were associated with increases more than 10 times the limit value considered instrumental. The detection limits were set as twice the standard deviation of the measured values.Following a more conventional protocol for testing accuracy of the proposed method, the reference material SRM 1649b from the National Institute of Standards and Technology of USA (NIST) was selected for initial comparisons. The reference material SRM 1649b provides typical urban dust matrix for evaluating the reliability of analytical methods. Four aliquots of 10 mg were subject to analytical protocol by agitation using dichloromethane as solvent. PAH mass fraction values were in the range of those certified (Table2s, ESM), except for benzo[a]pyrene and dibenzo[ah]anthracene, with lower recoveries (close to 70%). However, the recovery values were considered acceptable according to the published acceptable recovery percentages as a function of the analyte concentration [12].Regarding NPAH results, only 1-nitro pyrene could be compared and a recovery of 80% was reached.


### 3.2. Measured PACs and Test for Assessing Trueness

The study for comparison among the investigated extraction methods was carried out to evaluate the performance of the developed procedure by agitation/sonication. In particular, ambient air samples were divided into four subsamples, weighed, and then considered as independent or parallel subsamples for analysis, as previously indicated.Measured PAC masses (in ng) per subsample obtained from the investigated extraction methods are compiled in Table1. The sum of masses for selected PAHs and NPAHs ranged from 24 ng to 800 ng and 5.5 ng to 240 ng, respectively. The individual quantities of most abundant PAHs varied between 1 and 150 ng and together they accounted for 75% of total mass. The most abundant PAHs were 2-MeNa, Ph, Fl, Pyr, Cry, BbF, and BghiP. Ace and F were not detected. Regarding selected NPAHs, 9-NPh and 1-NPyr were the most abundant, which ranged from 1.2 ng to 100 ng.Table 1
Measured mass of selected PACs (ng) obtained from analysis of subsamples using the investigated extraction procedures. PM (mg) andn correspond to PM mass loading expressed in mg and sample number, respectively. The abbreviations for PACs and extraction procedures are expanded in Reagents and Standards and Samples and Preparation, respectively.
n
PM (mg)
extr
1-MeNa
2-MeNa
Phe
An
Fl
Pyr
BaA
Cry
BbF
BkF
BaP
DBA
BghiP
3-NPh
9-NPh
1-NPyr
3-NFl
ΣPAH
ΣNPAH
1
3.9
MAE
7.1
75
46
0.18
34
47
23
77
113
58
40
2.9
95
2.3
7.8
13
9.1
618
32
3.7
PLE
8.5
89
34
3.2
41
39
19
51
85
30
37
1.5
70
8.1
35
21
13
508
77
4.0
AG-D
7.5
83
36
2.7
30
28
14
55
96
51
53
1.8
62
<1
19
11
8.0
520
38
3.8
AG-DM
20
154
61
3.5
51
46
5.0
30
51
9.7
9.6
0.72
22
9.7
34.2
8.4
11
464
63
2
10
MAE
3.7
40
37
3.4
40
63
17
53
71
19
7.2
2.5
55
4.8
13
14
<1
411
32
10
PLE
5.7
53
39
2.8
44
46
16
41
53
11
13
1.8
37
7.0
22
29
10
367
69
10
AG-D
6.5
67
35
2.0
33
32
9.0
34
51
12
16
1.8
26
1.4
6.2
14
6.5
326
28
10
AG-DM
13
99
37
2.8
39
36
8.8
43
83
28
30
1.8
48
7.1
45
16
9.9
470
78
3
8.1
MAE
3.6
46
23
2.8
32
44
8.5
24
29
5.0
5.3
1.0
33
2.3
7.8
17
2.6
256
30
8.1
PLE
8.9
53
24
1.8
36
33
11
26
37
7.6
8.8
1.2
28
6.1
14
25
7.8
277
52
8.1
AG-D
6.9
75
36
1.5
38
23
6.2
22
35
8.2
7.5
1.1
20
1.1
13
15
7.5
280
37
8.1
AG-DM
18
68
28
1.6
33
29
1.2
14
34
6.2
5.9
0.52
14
2.6
16
16
7.1
254
41
4
10
MAE
<1
23
72
0.79
79
80
19
74
83
27
<0.1
3.1
51
3.7
16
26
1.6
512
48
10
PLE
8.3
99
76
4.5
69
68
26
65
78
27
33
2.9
71
4.1
15
19
8.8
629
47
10
AG-D
7.9
71
64
3.4
53
47
9.7
34
56
16
29
3.7
31
22
25
14
7.4
424
69
10
AG-DM
13
101
78
6.4
71
65
12
48
60
15
21
1.2
26
19
60
25
13
518
117
5
8.5
MAE
<1
80
36
0.60
47
56
22
64
152
52
<0.1
4.2
120
10
48
23
3.5
634
84
8.5
PLE
4.2
89
92
4.9
62
63
26
67
148
55
60
3.7
121
13
44
53
33
795
142
8.5
AG-D
2.6
76
60
3.1
46
44
15
63
144
53
76
<0.1
91
5.6
59
46
30
673
141
8.5
AG-DM
2.0
60
51
3.0
43
40
8.1
44
115
37
52
2.4
64
27
106
65
41
520
238
6
7.9
MAE
2.7
11
9.3
1.8
24
27
16
30
26
12
13
1.2
28
<1
17
<1
<1
201
17
7.9
MAE
4.5
10
6.0
1.4
22
25
15
29
24
11
11
1.5
25
<1
12
<1
<1
188
12
7.9
AG-D
20
16
11
2.4
31
52
19
38
31
13
18
1.3
43
<1
24
<1
1.4
295
26
7.9
AG-D
16
21
13
2.6
32
37
22
44
34
16
19
1.8
50
<1
26
<1
1.4
307
27
7
5.6
MAE
2.5
12
11
2.3
36
42
12
21
19
10
12
0.58
15
<1
10
<1
<1
196
10
5.6
MAE
2.0
13
13
2.7
43
48
14
23
21
11
10
0.67
15
<1
13
1.2
<1
216
14
5.6
AG-D
19
24
23
4.4
51
64
17
28
24
12
14
0.63
24
<1
21
2.1
1.2
305
25
5.6
AG-D
16
19
17
4.0
49
61
16
27
24
12
15
0.70
24
<1
19
<1
1.2
285
20
8
6.0
MAE
11
5.9
4.6
1.3
25
25
11
20
19
10
6.4
1.0
14
<1
7.3
<1
<1
154
7.3
6.0
MAE
14
9.6
5.9
0.79
23
22
9.4
19
17
9.2
3.7
0.67
12
<1
6.1
<1
<1
145
6.1
6.0
AG-D
26
8.1
8.0
1.5
22
25
10
17
17
8.4
8.1
0.58
16
<1
11
<1
<1
167
11
6.0
AG-D
15
7.1
8.7
1.6
22
26
11
17
16
8.9
8.6
0.63
16
<1
14
<1
<1
160
14
9
6.4
MAE
4.5
7.3
11
2.3
34
37
11
35
18
8.2
7.2
0.40
9.9
<1
5.9
1.2
<1
186
7.1
6.4
MAE
2.2
12
9.8
2.1
33
38
10
29
17
7.6
6.2
0.69
8.7
<1
6.2
1.3
<1
176
7.5
6.4
AG-D
26
11
9.8
2.6
36
46
13
22
19
8.6
8.9
0.54
12
<1
11
1.3
<1
214
13
6.4
AG-D
16
19
13
3.1
42
53
14
24
21
9.3
10
0.53
14
<1
13
1.5
<1
239
15
10
0.73
PLE
<1
<1
15
0.29
3.6
4.8
1.6
1.4
2.7
1.0
1.5
<0.1
8.8
<1
<1
<1
<1
41
0.73
PLE
<1
<1
13
0.13
2.9
4.5
1.2
1.2
2.5
1.0
1.3
<0.1
9.6
<1
<1
<1
<1
37
0.73
AG-D
<1
<1
9.7
0.08
1.7
2.6
0.8
0.7
1.6
0.6
0.84
<0.1
5.6
<1
<1
<1
<1
24
0.73
AG-D
<1
<1
10
0.08
1.4
2.5
0.8
0.7
1.4
0.6
0.82
<0.1
5.1
<1
<1
<1
<1
24
11
1.2
PLE
<1
<1
2.8
<0.03
3.2
6.4
2.2
2.6
5.9
3.0
5.3
1.9
19
1.5
2.3
2.4
5.5
53
11.7
1.2
PLE
<1
<1
7.1
<0.03
2.7
5.4
1.9
1.9
4.6
2.3
4.7
0.65
13
2.8
3.6
2.0
8.1
45
16.5
1.2
AG-D
<1
<1
2.7
0.07
2.9
5.1
1.8
1.9
4.2
2.0
4.0
0.55
11
<1
1.4
1.4
2.6
37
5.5
1.2
AG-D
<1
<1
5.2
0.09
3.5
6.1
2.3
2.6
5.5
2.6
4.9
0.53
15
1.1
2.1
3.1
4.3
48
11
12
2.7
PLE
<1
<1
9.5
1.1
33
64
18
25
34
18
30
3.9
74
7.1
7.9
1.46
49
310
65
2.7
PLE
<1
<1
8.8
0.85
25
51
15
19
25
13
26
2.4
59
9.1
9.6
1.87
39
245
60
2.7
AG-D
<1
<1
7.5
0.52
27
50
15
23
28
14
28
2.7
61
11
14
3.8
39
257
67
2.7
AG-D
<1
<1
5.9
0.74
25
50
15
19
27
13
26
2.2
60
6.3
6.5
2.8
28
244
44
Max
26
101
92
6.4
79
80
26
74
152
55
76
4.2
121
27
106
65
49
795
238
Min
2.0
5.9
2.7
0.07
1.4
2.5
0.79
0.70
1.4
0.56
0.82
0.40
5.1
1.08
1.42
1.17
1.20
24
5.5
2
2
1
0.06
0.12
0.08
0.04
0.4
0.12
0.016
0.04
0.08
0.08For assessing trueness of AG-D method, ambient air samples 1–5 were divided and considered as independent subsamples for extraction by MAE, PLE, AG-D, and AG-DM methods. Mass aerosol loadings per subsample were about 40 mg (samples 2–5) and 16 mg (sample 1), so selected PACs were determined at around 10 mg per subsample and 4 mg per subsample after cutting it (Table1).The experiment was designed as follows: the first subsample was extracted by MAE, the second one by PLE, the third by AG-D, and the last one by AG-DM. A preliminary comparison was performed by grouping results by extraction technique and deducing mean of measured masses (Figure1). As can be seen, results were in good agreement. In contrast to results obtained from SRM 1649 analyses, no differences were observed among investigated methods for benzo[a]pyrene and dibenzo[ah]anthracene measurements. When applying the agitation method using a mix of dichloromethane/methanol as solvent, high levels of the more volatile PAHs were obtained, while heavier PAHs (BaA, Cry, BkF, BaP, and BghiP) provided lower response. This is coherent to the more polar character of solvent and explains the difficulty for extracting the heavier and less polar PAHs. Similarly, NPAHs (more polar) corresponded to higher levels using methanol as modifier for agitation extraction.Figure 1
Mean values (ng per subsample) of selected PACs measured studied extraction methods. The abbreviations for PACs and extraction procedures are expanded in Reagents and Standards and Samples and Preparation, respectively.Taking PLE results as reference values, we considered testing proposed method as the closeness of agreement between both AG-D and PLE results [16, 17]. In detail, the calculation of recoveries (Ri) was performed by computing the individual ratios for each PAC according to(1)Ri=xAGDxPLE,where xAGD and xPLE are the mass of quantified PAC (ng) after agitation/sonication and pressurized liquid extractions using dichloromethane as solvent, respectively. Results obtained from the analyses of samples 1–5 and 9–11 were used for calculation. Quantified mass of selected PAC were between 1 and 75 ng (Table 1).In order to detect outlier values, Grubb’s test was applied to each data set. To verify whether the minimum or the maximum value was acceptable, the ratioGi was calculated according to (2)Gi=Rm/M-R-isd,whereRm/M denotes the maximum or minimum ratio value, R-i is the mean ratio, and sd refers to the standard deviation of data set. The value from this equation was compared with two-tailed critical ttab, for (n-1) degrees of freedom at 95% confidence (where n is the number of values). WhenGi values were higher than ttab, they were rejected and a recalculation was carried out. Of the data set, 4 outliers were discarded. Final results are included in Table 2. Mean values were very acceptable for PAHs (>80%). But lower recoveries were obtained for the less abundant PAHs and NPAHs.Table 2
Max, min, and mean ratios estimated by taking PLE results as reference values. Relative standard deviation (sd) and statistict values (t) are also included for evaluating outlier data and whether the mean overall ratio was significantly different from 1.0. Tabulated two-side value (ttab) according to (n-1) freedom degrees is noted (ttab). Empty cells correspond to outlier data or measured values below detection limits. Mean recoveries established by taking PLE results as reference values are noted in bold.
Sample
Ratio AG-D/PLE
1-MeNa
2-MeNa
Phe
An
Fl
Pyr
BaA
Cry
BbF
BkF
BaP
DBA
BghiP
3-NPh
9-NPh
1-NPyr
3-NFl
1
R
1
0.89
0.93
1.07
0.84
0.73
0.71
0.72
1.07
1.13
0.88
0.53
0.54
0.60
2
R
2
1.14
1.25
0.89
0.72
0.75
0.70
0.55
0.83
0.95
1.07
1.21
1.00
0.70
0.21
0.28
0.48
0.63
3
R
3
0.77
0.70
0.55
0.83
0.95
1.07
0.86
0.97
0.70
0.18
0.95
0.61
0.96
4
R
4
0.96
0.72
0.84
0.75
0.77
0.68
0.37
0.52
0.72
0.57
0.87
1.25
0.44
0.73
0.85
9
R
5
0.62
0.85
0.65
0.63
0.73
0.71
0.60
0.95
0.97
0.96
1.25
0.75
0.43
1.35
0.87
0.92
10
R
6
0.71
0.40
0.48
0.55
0.57
0.52
0.57
0.56
0.60
0.58
11
R
7
0.80
1.06
0.93
1.01
0.99
0.92
0.87
0.90
0.42
0.80
0.51
0.59
1.02
0.51
12
R
8
0.73
0.64
0.90
0.88
0.91
0.95
0.92
0.88
0.95
0.78
0.91
1.16
0.76
Max
1.14
1.25
1.07
0.84
1.06
0.93
1.01
1.07
1.13
1.07
1.25
1.25
0.91
0.51
1.35
1.02
0.96
Min
0.62
0.72
0.65
0.40
0.48
0.55
0.37
0.52
0.57
0.56
0.60
0.42
0.44
0.18
0.28
0.48
0.51
Mean
0.87
0.94
0.83
0.67
0.75
0.71
0.62
0.82
0.89
0.85
0.95
0.91
0.69
0.33
0.74
0.71
0.75
sd
0.20
0.22
0.15
0.17
0.19
0.11
0.20
0.22
0.18
0.23
0.24
0.35
0.15
0.16
0.42
0.21
0.19
n
5
4
6
5
6
7
7
7
7
6
6
4
7
4
5
6
6
G
M
1.35
1.38
1.65
1.02
1.66
1.94
1.93
1.18
1.32
0.94
1.24
0.98
1.50
1.07
1.47
1.50
1.15
G
m
1.30
0.98
1.20
1.60
1.48
1.46
1.29
1.36
1.71
1.24
1.42
1.40
1.76
0.93
1.11
1.11
1.23
t
0.291
0.138
0.483
0.881
0.544
0.968
0.712
0.319
0.232
0.264
0.089
0.127
0.791
2.027
0.277
0.575
0.549
t
t
a
b
2.132
2.353
2.015
2.132
2.015
1.943
1.943
1.943
1.943
2.015
2.015
2.353
1.943
2.353
2.132
2.015
2.015A test of significance was then used for verifying if the ratio between the value of the mean of recoveriesR-i and its relative standard deviation (sd) was significantly different. The average recovery was tested for significance by using the Student’s t-test:(3)t=R-i-1srec/n.If the obtained t value is lower than the tabulated value for (n-1) degrees of freedom at 95% significance level, we consider the results from PLE and AG-D methods comparable for PAC measurements. Only Pyr and BghiP corresponded to slightly lower values from the proposed method. However, no significant differences were found (Table 2).Regarding NPAH results, lower ratios were found, especially for 3-NPhe. The low mass levels together with the low mass loading are critical factors to the relatively low sensitivity of the derivatization method proposed for HPLC/FD analysis. These limitations should be improved by extraction with a mix of dichloromethane/methanol as previously mentioned.Taking into account the fact that NPAHs are minor contributors and the higher interest in quantifying less volatile PAHs, dichloromethane was chosen.


### 3.3. Estimation of the Intermediate Precision

From a more practical point of view, our goals were focused on reporting a range where results can be expected rather than detailed analyses of uncertainty sources. On this basis, intermediate precision was investigated by comparing results among the selected extraction methods and different PM mass loadings. In detail, seven ambient air samples were taken (samples 6–12) and each one was divided into four subsamples to analyze by pairs as follows (Figure2S, ESM).(i)
Comparison between MAE and AG-D procedures: Four samples with PM mass loadings between 5 and 8 mg were divided and analyzed by pairs in different days. In particular, results corresponded to the analyses of samples 5, 6, 7, and 8 (Table3).(ii)
Comparison between PLE and AG-D procedures: Analogously, three ambient air samples with lower PM mass loadings (0.7–2.8 mg) were divided and analyzed by pairs in different days. These low amounts of PM were chosen to investigate the scope of the developed method. In particular, results corresponded to the analyses of samples 9, 10, and 11 (Table3).Table 3
Values of intermediate precision deduced from results of the investigated extraction procedures using dichloromethane as solvent and different mass loadings of PM sample estimations were based on the pooled standard deviations (RSDpooled) and standard deviations of the normalized differences (RSDdiff/mean) calculations. Last column includes standard deviation obtained in previous work which was calculated as a standard deviation type A (sd/2) and obtained from the analysis of four subsamples split from a single sample.
PM mass loading (mg)
5–8 mg
0.7–2.8 mg
5 mg
RSD pooled
R
S
D
d
i
f
f
/
m
e
a
n
RSD pooled
R
S
D
d
i
f
f
/
m
e
a
n
sd/2
MAE
AG-D
MAE
AG-D
PLE
AG-D
PLE
AG-D
1-MeNa
0.32
0.27
0.52
0.17
2-MeNa
0.24
0.23
0.27
0.35
Phe
0.19
0.14
0.31
0.22
0.36
0.28
0.57
0.45
0.19
An
0.20
0.09
0.28
0.12
0.41
0.18
0.34
0.20
Fl
0.19
0.14
0.31
0.22
0.36
0.28
0.57
0.45
0.19
Pyr
0.20
0.09
0.28
0.12
0.18
0.34
0.20
0.18
BaA
0.07
0.06
0.11
0.08
0.15
0.10
0.05
0.17
0.15
Cry
0.07
0.14
0.11
0.21
0.12
0.08
0.08
0.12
0.16
BbF
0.08
0.06
0.12
0.07
0.15
0.10
0.08
0.16
0.15
BkF
0.09
0.07
0.12
0.09
0.18
0.15
0.08
0.26
0.15
BaP
0.06
0.06
0.09
0.08
0.16
0.12
0.11
0.19
0.33
DBA
0.06
0.06
0.09
0.06
0.16
0.11
0.11
0.17
BghiP
0.21
0.07
0.19
0.03
0.11
0.09
0.02
0.15
3-NPh
0.21
9-NPh
0.08
0.07
0.07
0.08
0.17
0.11
0.22
0.18
0.15
1-NPyr
0.28
3-NFl
0.20
0.57
0.28Once data compiled, intermediate precision was calculated based on two practical approaches for a better estimation:(1)
Percentage difference: by estimating a relative standard deviation of the normalized differences (Example A4, EURACHEM [18])(2)
RSDpooled: by calculating a pooled relative deviationNamely, the precision was estimated over a time period, using the standard deviation of pairs of results. The analyses corresponded to 7 samples analyzed 4 times over a period of time.


#### 3.3.1. Percentage Difference

The calculation was carried out by grouping results in pairs, according to each filter sample and extraction technique (Table3s, ESM). The normalized difference data was then calculated as the difference divided by mean value for each pair of data. These results really show a difference between two values of each pair as a percentage of both values. Therefore, they provide a measure of the overall variability between analyses. On this basis, calculating the standard deviation of the normalized differences gives an estimation of relative standard uncertainty due to run-to-run variation of the overall analytical process.


#### 3.3.2.RSDpooled

Results in pairs, according to the duplicate analyses of each sample and extraction technique, were grouped to deduce relative standard deviations. The mean relative standard deviation (RSDpooled) was then expressed according to (4)RSDpooled=∑rsdi2ni-1∑ni-1,where rsdi corresponds to relative standard deviation deduced from each ambient air sample, i, in duplicate.The results of intermediate precision obtained are compiled in Table3. Empty cells corresponded to values below detection limits. Also the relative deviations derived from previous work have been included [14]. It must be noted that these results were obtained from the analysis of a single sample after cutting four subsamples, which were treated and analyzed as independent samples. Deviation values deduced with this approach were slightly higher. Measures from a sole ambient air sample can lead to implying an overestimation of the analytical variability by mathematical calculation itself.In general, good concordance was found between results obtained from both mathematical approaches.As can be seen,RSDpooled corresponding to agitation results were generally below 10% for PAH and 9-NPhe with masses above 10 ng (subsamples with PM mass loadings between 5 and 8 mg, Table 3). Higher relative standard deviation was found (10–15%) when mass levels were between 1 and 5 ng (subsamples with PM mass loadings between 0.7 and 2.8 mg, Table 3). In the case of selected NPAHs, only 9-NPhe was quantified and no enough data were obtained for the rest of investigated NPAHs. For microwave method, lower values and increases of variability (up to 50%) were found for the isomers of methyl naphthalene, the most volatile. These observations likely correspond to losses due to the applied heat and the more volatile nature of these compounds.For NPAHs, intermediate deviation values were higher, around 20%, which is justified by the increase of experimental variability associated with additional derivatization step for chromatographic analysis. Calculations could be performed when higher NPAH mass levels were measured (subsamples with PM mass loadings between 5 and 8 mg, Table3). PLE shows a slightly lower variability of the results. This is in agreement with the higher efficiency of PLE extraction and difficulty of the extraction of several PACs from particulate matter [19].


### 3.4. Application to Real PM Samples with Very Low Mass Loadings

As a part of characterization studies from the turbofan engine exhaust [20], aerosols on quartz fiber filter and aluminum (Al) foils (47 mm) were collected by means of a Low Pressure Impactor and a filter holder, respectively. Measurements were focused on preliminary studies on particle size distributions of PM. Exhaust PM sampling was conducted to explore detection limits of analysis and optimize parameters such as time and flow at different engine conditions. This is often a limiting factor to obtain sufficient sample mass for analysis. In detail, the effectiveness of the proposed method was evaluated determining selected PACs in PM aerosol samples collected on quartz filters (at mass loadings of 0.2 mg) and Al foils (mass loadings of 0.02 mg).From an analytical point of view, these samples can be considered as a worst extreme case: analytical methods have to allow the quantification of analyte in PM mass loadings of a few micrograms. Figure2 corresponds to representative chromatograms obtained from the analyses of final extract volumes of 100 µL. It should be noted that the good resolution of chromatograms was obtained without the need of any additional step for clean-up. No selected NPAHs were detected in this kind of samples. The sum of PAH masses reached 75 ng in quartz filter samples, while those for Al foils were only up to 2 ng. Regarding PM collected on Al foils, very slight differences were obtained between samples, with most of the compounds being below 0.1 ng. All the measurements remained slightly above the quantified limits established, so we considered these levels as detectable, not for quantitative purposes. Table 4 includes selected PAH measurements (ng) from the analyses of quartz filters. As can be seen, measured levels reached orders of units of nanogram, expressed as mass per analyzed filter.Table 4
PAH measurements (ng) in PM mass loadings of about 0.2 mg collected in quartz filters.
Sample
1-MeNa
2-MeNa
Phe
An
Fl
Pyr
BaA
Cry
BbF
BkF
BaP
DBA
BghiP
ΣPAH
Blank 1
<1
<1
<1
<0.03
<0.1
<0.1
<0.3
<0.05
<0.1
<0.1
<0.1
<0.1
<0.1
Blank 2
<1
<1
<1
0.82
<0.1
<0.1
<0.3
<0.05
<0.1
<0.1
<0.1
<0.1
<0.1
F1
<1
3.5
5.0
0.98
11
18
1.9
4.5
2.0
0.63
0.35
<0.1
1.0
49
F2
<1
<1
<1
4.5
0.73
18
30
1.6
7.0
5.8
1.8
0.70
0.20
70
F3
40
<1
4.8
0.38
6.7
12
3.7
3.1
2.5
0.95
0.38
<0.1
1.5
76
F4
46
<1
3.8
0.16
4.6
8.1
2.4
6.1
1.6
0.51
0.27
<0.1
0.50
74
F5
5.8
<1
9.8
0.38
2.5
4.0
0.38
0.30
0.54
0.17
0.22
<0.1
0.85
25Figure 2
Chromatograms of samples collected on aluminum (a) and quartz (b) substrates at very low concentration levels. The PAH abbreviations are expanded in Reagents and Standards.
(a)
(b)Due to the lack of analytical assessments in this extremely low range of mass loadings, the real extent of this burden remains to be determined in future researches. Therefore, these preliminary results need additional studies on the characterization of size distribution from aircraft aerosol emissions. These goals are the main objectives to continue further studies from aircraft aerosol emissions.


## 3.1. Quality Control

Instrumental limit of detection of each PAC was calculated as twice the standard deviations deduced from analyses of three standard solutions containing 1 ng mL−1 and 10 ng mL−1 for selected PAHs and NPAH, respectively. Derived values were between 0.5 ng mL−1 and 5 ng mL−1 for selected PAHs and NPAHs, respectively.The analysis of blanks (i.e., chemical analysis excluding PM samples) indicated the potential cross contamination for the more volatile PAHs if extreme care is not taken. Due to ubiquitous character of PAH compounds in the air, laboratory inventory, and low volumes of extracts, there is a risk of contamination. Therefore, to determine blank contents for each set of sample analyses is recommended. The very low levels require a rigorous protocol to ensure minimal sample contamination from vials, syringes, solvents, and drying step.Blank filters of glass fiber, quartz, and aluminum were handled in the same way as the samples by agitation procedure and these trials concluded to detectable contents only for the more volatile PAHs. In detail, 1-methyl-naphthalene, 2-methyl-naphthalene, and phenanthrene were associated with increases more than 10 times the limit value considered instrumental. The detection limits were set as twice the standard deviation of the measured values.Following a more conventional protocol for testing accuracy of the proposed method, the reference material SRM 1649b from the National Institute of Standards and Technology of USA (NIST) was selected for initial comparisons. The reference material SRM 1649b provides typical urban dust matrix for evaluating the reliability of analytical methods. Four aliquots of 10 mg were subject to analytical protocol by agitation using dichloromethane as solvent. PAH mass fraction values were in the range of those certified (Table2s, ESM), except for benzo[a]pyrene and dibenzo[ah]anthracene, with lower recoveries (close to 70%). However, the recovery values were considered acceptable according to the published acceptable recovery percentages as a function of the analyte concentration [12].Regarding NPAH results, only 1-nitro pyrene could be compared and a recovery of 80% was reached.


## 3.2. Measured PACs and Test for Assessing Trueness

The study for comparison among the investigated extraction methods was carried out to evaluate the performance of the developed procedure by agitation/sonication. In particular, ambient air samples were divided into four subsamples, weighed, and then considered as independent or parallel subsamples for analysis, as previously indicated.Measured PAC masses (in ng) per subsample obtained from the investigated extraction methods are compiled in Table1. The sum of masses for selected PAHs and NPAHs ranged from 24 ng to 800 ng and 5.5 ng to 240 ng, respectively. The individual quantities of most abundant PAHs varied between 1 and 150 ng and together they accounted for 75% of total mass. The most abundant PAHs were 2-MeNa, Ph, Fl, Pyr, Cry, BbF, and BghiP. Ace and F were not detected. Regarding selected NPAHs, 9-NPh and 1-NPyr were the most abundant, which ranged from 1.2 ng to 100 ng.Table 1
Measured mass of selected PACs (ng) obtained from analysis of subsamples using the investigated extraction procedures. PM (mg) andn correspond to PM mass loading expressed in mg and sample number, respectively. The abbreviations for PACs and extraction procedures are expanded in Reagents and Standards and Samples and Preparation, respectively.
n
PM (mg)
extr
1-MeNa
2-MeNa
Phe
An
Fl
Pyr
BaA
Cry
BbF
BkF
BaP
DBA
BghiP
3-NPh
9-NPh
1-NPyr
3-NFl
ΣPAH
ΣNPAH
1
3.9
MAE
7.1
75
46
0.18
34
47
23
77
113
58
40
2.9
95
2.3
7.8
13
9.1
618
32
3.7
PLE
8.5
89
34
3.2
41
39
19
51
85
30
37
1.5
70
8.1
35
21
13
508
77
4.0
AG-D
7.5
83
36
2.7
30
28
14
55
96
51
53
1.8
62
<1
19
11
8.0
520
38
3.8
AG-DM
20
154
61
3.5
51
46
5.0
30
51
9.7
9.6
0.72
22
9.7
34.2
8.4
11
464
63
2
10
MAE
3.7
40
37
3.4
40
63
17
53
71
19
7.2
2.5
55
4.8
13
14
<1
411
32
10
PLE
5.7
53
39
2.8
44
46
16
41
53
11
13
1.8
37
7.0
22
29
10
367
69
10
AG-D
6.5
67
35
2.0
33
32
9.0
34
51
12
16
1.8
26
1.4
6.2
14
6.5
326
28
10
AG-DM
13
99
37
2.8
39
36
8.8
43
83
28
30
1.8
48
7.1
45
16
9.9
470
78
3
8.1
MAE
3.6
46
23
2.8
32
44
8.5
24
29
5.0
5.3
1.0
33
2.3
7.8
17
2.6
256
30
8.1
PLE
8.9
53
24
1.8
36
33
11
26
37
7.6
8.8
1.2
28
6.1
14
25
7.8
277
52
8.1
AG-D
6.9
75
36
1.5
38
23
6.2
22
35
8.2
7.5
1.1
20
1.1
13
15
7.5
280
37
8.1
AG-DM
18
68
28
1.6
33
29
1.2
14
34
6.2
5.9
0.52
14
2.6
16
16
7.1
254
41
4
10
MAE
<1
23
72
0.79
79
80
19
74
83
27
<0.1
3.1
51
3.7
16
26
1.6
512
48
10
PLE
8.3
99
76
4.5
69
68
26
65
78
27
33
2.9
71
4.1
15
19
8.8
629
47
10
AG-D
7.9
71
64
3.4
53
47
9.7
34
56
16
29
3.7
31
22
25
14
7.4
424
69
10
AG-DM
13
101
78
6.4
71
65
12
48
60
15
21
1.2
26
19
60
25
13
518
117
5
8.5
MAE
<1
80
36
0.60
47
56
22
64
152
52
<0.1
4.2
120
10
48
23
3.5
634
84
8.5
PLE
4.2
89
92
4.9
62
63
26
67
148
55
60
3.7
121
13
44
53
33
795
142
8.5
AG-D
2.6
76
60
3.1
46
44
15
63
144
53
76
<0.1
91
5.6
59
46
30
673
141
8.5
AG-DM
2.0
60
51
3.0
43
40
8.1
44
115
37
52
2.4
64
27
106
65
41
520
238
6
7.9
MAE
2.7
11
9.3
1.8
24
27
16
30
26
12
13
1.2
28
<1
17
<1
<1
201
17
7.9
MAE
4.5
10
6.0
1.4
22
25
15
29
24
11
11
1.5
25
<1
12
<1
<1
188
12
7.9
AG-D
20
16
11
2.4
31
52
19
38
31
13
18
1.3
43
<1
24
<1
1.4
295
26
7.9
AG-D
16
21
13
2.6
32
37
22
44
34
16
19
1.8
50
<1
26
<1
1.4
307
27
7
5.6
MAE
2.5
12
11
2.3
36
42
12
21
19
10
12
0.58
15
<1
10
<1
<1
196
10
5.6
MAE
2.0
13
13
2.7
43
48
14
23
21
11
10
0.67
15
<1
13
1.2
<1
216
14
5.6
AG-D
19
24
23
4.4
51
64
17
28
24
12
14
0.63
24
<1
21
2.1
1.2
305
25
5.6
AG-D
16
19
17
4.0
49
61
16
27
24
12
15
0.70
24
<1
19
<1
1.2
285
20
8
6.0
MAE
11
5.9
4.6
1.3
25
25
11
20
19
10
6.4
1.0
14
<1
7.3
<1
<1
154
7.3
6.0
MAE
14
9.6
5.9
0.79
23
22
9.4
19
17
9.2
3.7
0.67
12
<1
6.1
<1
<1
145
6.1
6.0
AG-D
26
8.1
8.0
1.5
22
25
10
17
17
8.4
8.1
0.58
16
<1
11
<1
<1
167
11
6.0
AG-D
15
7.1
8.7
1.6
22
26
11
17
16
8.9
8.6
0.63
16
<1
14
<1
<1
160
14
9
6.4
MAE
4.5
7.3
11
2.3
34
37
11
35
18
8.2
7.2
0.40
9.9
<1
5.9
1.2
<1
186
7.1
6.4
MAE
2.2
12
9.8
2.1
33
38
10
29
17
7.6
6.2
0.69
8.7
<1
6.2
1.3
<1
176
7.5
6.4
AG-D
26
11
9.8
2.6
36
46
13
22
19
8.6
8.9
0.54
12
<1
11
1.3
<1
214
13
6.4
AG-D
16
19
13
3.1
42
53
14
24
21
9.3
10
0.53
14
<1
13
1.5
<1
239
15
10
0.73
PLE
<1
<1
15
0.29
3.6
4.8
1.6
1.4
2.7
1.0
1.5
<0.1
8.8
<1
<1
<1
<1
41
0.73
PLE
<1
<1
13
0.13
2.9
4.5
1.2
1.2
2.5
1.0
1.3
<0.1
9.6
<1
<1
<1
<1
37
0.73
AG-D
<1
<1
9.7
0.08
1.7
2.6
0.8
0.7
1.6
0.6
0.84
<0.1
5.6
<1
<1
<1
<1
24
0.73
AG-D
<1
<1
10
0.08
1.4
2.5
0.8
0.7
1.4
0.6
0.82
<0.1
5.1
<1
<1
<1
<1
24
11
1.2
PLE
<1
<1
2.8
<0.03
3.2
6.4
2.2
2.6
5.9
3.0
5.3
1.9
19
1.5
2.3
2.4
5.5
53
11.7
1.2
PLE
<1
<1
7.1
<0.03
2.7
5.4
1.9
1.9
4.6
2.3
4.7
0.65
13
2.8
3.6
2.0
8.1
45
16.5
1.2
AG-D
<1
<1
2.7
0.07
2.9
5.1
1.8
1.9
4.2
2.0
4.0
0.55
11
<1
1.4
1.4
2.6
37
5.5
1.2
AG-D
<1
<1
5.2
0.09
3.5
6.1
2.3
2.6
5.5
2.6
4.9
0.53
15
1.1
2.1
3.1
4.3
48
11
12
2.7
PLE
<1
<1
9.5
1.1
33
64
18
25
34
18
30
3.9
74
7.1
7.9
1.46
49
310
65
2.7
PLE
<1
<1
8.8
0.85
25
51
15
19
25
13
26
2.4
59
9.1
9.6
1.87
39
245
60
2.7
AG-D
<1
<1
7.5
0.52
27
50
15
23
28
14
28
2.7
61
11
14
3.8
39
257
67
2.7
AG-D
<1
<1
5.9
0.74
25
50
15
19
27
13
26
2.2
60
6.3
6.5
2.8
28
244
44
Max
26
101
92
6.4
79
80
26
74
152
55
76
4.2
121
27
106
65
49
795
238
Min
2.0
5.9
2.7
0.07
1.4
2.5
0.79
0.70
1.4
0.56
0.82
0.40
5.1
1.08
1.42
1.17
1.20
24
5.5
2
2
1
0.06
0.12
0.08
0.04
0.4
0.12
0.016
0.04
0.08
0.08For assessing trueness of AG-D method, ambient air samples 1–5 were divided and considered as independent subsamples for extraction by MAE, PLE, AG-D, and AG-DM methods. Mass aerosol loadings per subsample were about 40 mg (samples 2–5) and 16 mg (sample 1), so selected PACs were determined at around 10 mg per subsample and 4 mg per subsample after cutting it (Table1).The experiment was designed as follows: the first subsample was extracted by MAE, the second one by PLE, the third by AG-D, and the last one by AG-DM. A preliminary comparison was performed by grouping results by extraction technique and deducing mean of measured masses (Figure1). As can be seen, results were in good agreement. In contrast to results obtained from SRM 1649 analyses, no differences were observed among investigated methods for benzo[a]pyrene and dibenzo[ah]anthracene measurements. When applying the agitation method using a mix of dichloromethane/methanol as solvent, high levels of the more volatile PAHs were obtained, while heavier PAHs (BaA, Cry, BkF, BaP, and BghiP) provided lower response. This is coherent to the more polar character of solvent and explains the difficulty for extracting the heavier and less polar PAHs. Similarly, NPAHs (more polar) corresponded to higher levels using methanol as modifier for agitation extraction.Figure 1
Mean values (ng per subsample) of selected PACs measured studied extraction methods. The abbreviations for PACs and extraction procedures are expanded in Reagents and Standards and Samples and Preparation, respectively.Taking PLE results as reference values, we considered testing proposed method as the closeness of agreement between both AG-D and PLE results [16, 17]. In detail, the calculation of recoveries (Ri) was performed by computing the individual ratios for each PAC according to(1)Ri=xAGDxPLE,where xAGD and xPLE are the mass of quantified PAC (ng) after agitation/sonication and pressurized liquid extractions using dichloromethane as solvent, respectively. Results obtained from the analyses of samples 1–5 and 9–11 were used for calculation. Quantified mass of selected PAC were between 1 and 75 ng (Table 1).In order to detect outlier values, Grubb’s test was applied to each data set. To verify whether the minimum or the maximum value was acceptable, the ratioGi was calculated according to (2)Gi=Rm/M-R-isd,whereRm/M denotes the maximum or minimum ratio value, R-i is the mean ratio, and sd refers to the standard deviation of data set. The value from this equation was compared with two-tailed critical ttab, for (n-1) degrees of freedom at 95% confidence (where n is the number of values). WhenGi values were higher than ttab, they were rejected and a recalculation was carried out. Of the data set, 4 outliers were discarded. Final results are included in Table 2. Mean values were very acceptable for PAHs (>80%). But lower recoveries were obtained for the less abundant PAHs and NPAHs.Table 2
Max, min, and mean ratios estimated by taking PLE results as reference values. Relative standard deviation (sd) and statistict values (t) are also included for evaluating outlier data and whether the mean overall ratio was significantly different from 1.0. Tabulated two-side value (ttab) according to (n-1) freedom degrees is noted (ttab). Empty cells correspond to outlier data or measured values below detection limits. Mean recoveries established by taking PLE results as reference values are noted in bold.
Sample
Ratio AG-D/PLE
1-MeNa
2-MeNa
Phe
An
Fl
Pyr
BaA
Cry
BbF
BkF
BaP
DBA
BghiP
3-NPh
9-NPh
1-NPyr
3-NFl
1
R
1
0.89
0.93
1.07
0.84
0.73
0.71
0.72
1.07
1.13
0.88
0.53
0.54
0.60
2
R
2
1.14
1.25
0.89
0.72
0.75
0.70
0.55
0.83
0.95
1.07
1.21
1.00
0.70
0.21
0.28
0.48
0.63
3
R
3
0.77
0.70
0.55
0.83
0.95
1.07
0.86
0.97
0.70
0.18
0.95
0.61
0.96
4
R
4
0.96
0.72
0.84
0.75
0.77
0.68
0.37
0.52
0.72
0.57
0.87
1.25
0.44
0.73
0.85
9
R
5
0.62
0.85
0.65
0.63
0.73
0.71
0.60
0.95
0.97
0.96
1.25
0.75
0.43
1.35
0.87
0.92
10
R
6
0.71
0.40
0.48
0.55
0.57
0.52
0.57
0.56
0.60
0.58
11
R
7
0.80
1.06
0.93
1.01
0.99
0.92
0.87
0.90
0.42
0.80
0.51
0.59
1.02
0.51
12
R
8
0.73
0.64
0.90
0.88
0.91
0.95
0.92
0.88
0.95
0.78
0.91
1.16
0.76
Max
1.14
1.25
1.07
0.84
1.06
0.93
1.01
1.07
1.13
1.07
1.25
1.25
0.91
0.51
1.35
1.02
0.96
Min
0.62
0.72
0.65
0.40
0.48
0.55
0.37
0.52
0.57
0.56
0.60
0.42
0.44
0.18
0.28
0.48
0.51
Mean
0.87
0.94
0.83
0.67
0.75
0.71
0.62
0.82
0.89
0.85
0.95
0.91
0.69
0.33
0.74
0.71
0.75
sd
0.20
0.22
0.15
0.17
0.19
0.11
0.20
0.22
0.18
0.23
0.24
0.35
0.15
0.16
0.42
0.21
0.19
n
5
4
6
5
6
7
7
7
7
6
6
4
7
4
5
6
6
G
M
1.35
1.38
1.65
1.02
1.66
1.94
1.93
1.18
1.32
0.94
1.24
0.98
1.50
1.07
1.47
1.50
1.15
G
m
1.30
0.98
1.20
1.60
1.48
1.46
1.29
1.36
1.71
1.24
1.42
1.40
1.76
0.93
1.11
1.11
1.23
t
0.291
0.138
0.483
0.881
0.544
0.968
0.712
0.319
0.232
0.264
0.089
0.127
0.791
2.027
0.277
0.575
0.549
t
t
a
b
2.132
2.353
2.015
2.132
2.015
1.943
1.943
1.943
1.943
2.015
2.015
2.353
1.943
2.353
2.132
2.015
2.015A test of significance was then used for verifying if the ratio between the value of the mean of recoveriesR-i and its relative standard deviation (sd) was significantly different. The average recovery was tested for significance by using the Student’s t-test:(3)t=R-i-1srec/n.If the obtained t value is lower than the tabulated value for (n-1) degrees of freedom at 95% significance level, we consider the results from PLE and AG-D methods comparable for PAC measurements. Only Pyr and BghiP corresponded to slightly lower values from the proposed method. However, no significant differences were found (Table 2).Regarding NPAH results, lower ratios were found, especially for 3-NPhe. The low mass levels together with the low mass loading are critical factors to the relatively low sensitivity of the derivatization method proposed for HPLC/FD analysis. These limitations should be improved by extraction with a mix of dichloromethane/methanol as previously mentioned.Taking into account the fact that NPAHs are minor contributors and the higher interest in quantifying less volatile PAHs, dichloromethane was chosen.


## 3.3. Estimation of the Intermediate Precision

From a more practical point of view, our goals were focused on reporting a range where results can be expected rather than detailed analyses of uncertainty sources. On this basis, intermediate precision was investigated by comparing results among the selected extraction methods and different PM mass loadings. In detail, seven ambient air samples were taken (samples 6–12) and each one was divided into four subsamples to analyze by pairs as follows (Figure2S, ESM).(i)
Comparison between MAE and AG-D procedures: Four samples with PM mass loadings between 5 and 8 mg were divided and analyzed by pairs in different days. In particular, results corresponded to the analyses of samples 5, 6, 7, and 8 (Table3).(ii)
Comparison between PLE and AG-D procedures: Analogously, three ambient air samples with lower PM mass loadings (0.7–2.8 mg) were divided and analyzed by pairs in different days. These low amounts of PM were chosen to investigate the scope of the developed method. In particular, results corresponded to the analyses of samples 9, 10, and 11 (Table3).Table 3
Values of intermediate precision deduced from results of the investigated extraction procedures using dichloromethane as solvent and different mass loadings of PM sample estimations were based on the pooled standard deviations (RSDpooled) and standard deviations of the normalized differences (RSDdiff/mean) calculations. Last column includes standard deviation obtained in previous work which was calculated as a standard deviation type A (sd/2) and obtained from the analysis of four subsamples split from a single sample.
PM mass loading (mg)
5–8 mg
0.7–2.8 mg
5 mg
RSD pooled
R
S
D
d
i
f
f
/
m
e
a
n
RSD pooled
R
S
D
d
i
f
f
/
m
e
a
n
sd/2
MAE
AG-D
MAE
AG-D
PLE
AG-D
PLE
AG-D
1-MeNa
0.32
0.27
0.52
0.17
2-MeNa
0.24
0.23
0.27
0.35
Phe
0.19
0.14
0.31
0.22
0.36
0.28
0.57
0.45
0.19
An
0.20
0.09
0.28
0.12
0.41
0.18
0.34
0.20
Fl
0.19
0.14
0.31
0.22
0.36
0.28
0.57
0.45
0.19
Pyr
0.20
0.09
0.28
0.12
0.18
0.34
0.20
0.18
BaA
0.07
0.06
0.11
0.08
0.15
0.10
0.05
0.17
0.15
Cry
0.07
0.14
0.11
0.21
0.12
0.08
0.08
0.12
0.16
BbF
0.08
0.06
0.12
0.07
0.15
0.10
0.08
0.16
0.15
BkF
0.09
0.07
0.12
0.09
0.18
0.15
0.08
0.26
0.15
BaP
0.06
0.06
0.09
0.08
0.16
0.12
0.11
0.19
0.33
DBA
0.06
0.06
0.09
0.06
0.16
0.11
0.11
0.17
BghiP
0.21
0.07
0.19
0.03
0.11
0.09
0.02
0.15
3-NPh
0.21
9-NPh
0.08
0.07
0.07
0.08
0.17
0.11
0.22
0.18
0.15
1-NPyr
0.28
3-NFl
0.20
0.57
0.28Once data compiled, intermediate precision was calculated based on two practical approaches for a better estimation:(1)
Percentage difference: by estimating a relative standard deviation of the normalized differences (Example A4, EURACHEM [18])(2)
RSDpooled: by calculating a pooled relative deviationNamely, the precision was estimated over a time period, using the standard deviation of pairs of results. The analyses corresponded to 7 samples analyzed 4 times over a period of time.


### 3.3.1. Percentage Difference

The calculation was carried out by grouping results in pairs, according to each filter sample and extraction technique (Table3s, ESM). The normalized difference data was then calculated as the difference divided by mean value for each pair of data. These results really show a difference between two values of each pair as a percentage of both values. Therefore, they provide a measure of the overall variability between analyses. On this basis, calculating the standard deviation of the normalized differences gives an estimation of relative standard uncertainty due to run-to-run variation of the overall analytical process.


### 3.3.2.RSDpooled

Results in pairs, according to the duplicate analyses of each sample and extraction technique, were grouped to deduce relative standard deviations. The mean relative standard deviation (RSDpooled) was then expressed according to (4)RSDpooled=∑rsdi2ni-1∑ni-1,where rsdi corresponds to relative standard deviation deduced from each ambient air sample, i, in duplicate.The results of intermediate precision obtained are compiled in Table3. Empty cells corresponded to values below detection limits. Also the relative deviations derived from previous work have been included [14]. It must be noted that these results were obtained from the analysis of a single sample after cutting four subsamples, which were treated and analyzed as independent samples. Deviation values deduced with this approach were slightly higher. Measures from a sole ambient air sample can lead to implying an overestimation of the analytical variability by mathematical calculation itself.In general, good concordance was found between results obtained from both mathematical approaches.As can be seen,RSDpooled corresponding to agitation results were generally below 10% for PAH and 9-NPhe with masses above 10 ng (subsamples with PM mass loadings between 5 and 8 mg, Table 3). Higher relative standard deviation was found (10–15%) when mass levels were between 1 and 5 ng (subsamples with PM mass loadings between 0.7 and 2.8 mg, Table 3). In the case of selected NPAHs, only 9-NPhe was quantified and no enough data were obtained for the rest of investigated NPAHs. For microwave method, lower values and increases of variability (up to 50%) were found for the isomers of methyl naphthalene, the most volatile. These observations likely correspond to losses due to the applied heat and the more volatile nature of these compounds.For NPAHs, intermediate deviation values were higher, around 20%, which is justified by the increase of experimental variability associated with additional derivatization step for chromatographic analysis. Calculations could be performed when higher NPAH mass levels were measured (subsamples with PM mass loadings between 5 and 8 mg, Table3). PLE shows a slightly lower variability of the results. This is in agreement with the higher efficiency of PLE extraction and difficulty of the extraction of several PACs from particulate matter [19].


## 3.3.1. Percentage Difference

The calculation was carried out by grouping results in pairs, according to each filter sample and extraction technique (Table3s, ESM). The normalized difference data was then calculated as the difference divided by mean value for each pair of data. These results really show a difference between two values of each pair as a percentage of both values. Therefore, they provide a measure of the overall variability between analyses. On this basis, calculating the standard deviation of the normalized differences gives an estimation of relative standard uncertainty due to run-to-run variation of the overall analytical process.


## 3.3.2.RSDpooled

Results in pairs, according to the duplicate analyses of each sample and extraction technique, were grouped to deduce relative standard deviations. The mean relative standard deviation (RSDpooled) was then expressed according to (4)RSDpooled=∑rsdi2ni-1∑ni-1,where rsdi corresponds to relative standard deviation deduced from each ambient air sample, i, in duplicate.The results of intermediate precision obtained are compiled in Table3. Empty cells corresponded to values below detection limits. Also the relative deviations derived from previous work have been included [14]. It must be noted that these results were obtained from the analysis of a single sample after cutting four subsamples, which were treated and analyzed as independent samples. Deviation values deduced with this approach were slightly higher. Measures from a sole ambient air sample can lead to implying an overestimation of the analytical variability by mathematical calculation itself.In general, good concordance was found between results obtained from both mathematical approaches.As can be seen,RSDpooled corresponding to agitation results were generally below 10% for PAH and 9-NPhe with masses above 10 ng (subsamples with PM mass loadings between 5 and 8 mg, Table 3). Higher relative standard deviation was found (10–15%) when mass levels were between 1 and 5 ng (subsamples with PM mass loadings between 0.7 and 2.8 mg, Table 3). In the case of selected NPAHs, only 9-NPhe was quantified and no enough data were obtained for the rest of investigated NPAHs. For microwave method, lower values and increases of variability (up to 50%) were found for the isomers of methyl naphthalene, the most volatile. These observations likely correspond to losses due to the applied heat and the more volatile nature of these compounds.For NPAHs, intermediate deviation values were higher, around 20%, which is justified by the increase of experimental variability associated with additional derivatization step for chromatographic analysis. Calculations could be performed when higher NPAH mass levels were measured (subsamples with PM mass loadings between 5 and 8 mg, Table3). PLE shows a slightly lower variability of the results. This is in agreement with the higher efficiency of PLE extraction and difficulty of the extraction of several PACs from particulate matter [19].


## 3.4. Application to Real PM Samples with Very Low Mass Loadings

As a part of characterization studies from the turbofan engine exhaust [20], aerosols on quartz fiber filter and aluminum (Al) foils (47 mm) were collected by means of a Low Pressure Impactor and a filter holder, respectively. Measurements were focused on preliminary studies on particle size distributions of PM. Exhaust PM sampling was conducted to explore detection limits of analysis and optimize parameters such as time and flow at different engine conditions. This is often a limiting factor to obtain sufficient sample mass for analysis. In detail, the effectiveness of the proposed method was evaluated determining selected PACs in PM aerosol samples collected on quartz filters (at mass loadings of 0.2 mg) and Al foils (mass loadings of 0.02 mg).From an analytical point of view, these samples can be considered as a worst extreme case: analytical methods have to allow the quantification of analyte in PM mass loadings of a few micrograms. Figure2 corresponds to representative chromatograms obtained from the analyses of final extract volumes of 100 µL. It should be noted that the good resolution of chromatograms was obtained without the need of any additional step for clean-up. No selected NPAHs were detected in this kind of samples. The sum of PAH masses reached 75 ng in quartz filter samples, while those for Al foils were only up to 2 ng. Regarding PM collected on Al foils, very slight differences were obtained between samples, with most of the compounds being below 0.1 ng. All the measurements remained slightly above the quantified limits established, so we considered these levels as detectable, not for quantitative purposes. Table 4 includes selected PAH measurements (ng) from the analyses of quartz filters. As can be seen, measured levels reached orders of units of nanogram, expressed as mass per analyzed filter.Table 4
PAH measurements (ng) in PM mass loadings of about 0.2 mg collected in quartz filters.
Sample
1-MeNa
2-MeNa
Phe
An
Fl
Pyr
BaA
Cry
BbF
BkF
BaP
DBA
BghiP
ΣPAH
Blank 1
<1
<1
<1
<0.03
<0.1
<0.1
<0.3
<0.05
<0.1
<0.1
<0.1
<0.1
<0.1
Blank 2
<1
<1
<1
0.82
<0.1
<0.1
<0.3
<0.05
<0.1
<0.1
<0.1
<0.1
<0.1
F1
<1
3.5
5.0
0.98
11
18
1.9
4.5
2.0
0.63
0.35
<0.1
1.0
49
F2
<1
<1
<1
4.5
0.73
18
30
1.6
7.0
5.8
1.8
0.70
0.20
70
F3
40
<1
4.8
0.38
6.7
12
3.7
3.1
2.5
0.95
0.38
<0.1
1.5
76
F4
46
<1
3.8
0.16
4.6
8.1
2.4
6.1
1.6
0.51
0.27
<0.1
0.50
74
F5
5.8
<1
9.8
0.38
2.5
4.0
0.38
0.30
0.54
0.17
0.22
<0.1
0.85
25Figure 2
Chromatograms of samples collected on aluminum (a) and quartz (b) substrates at very low concentration levels. The PAH abbreviations are expanded in Reagents and Standards.
(a)
(b)Due to the lack of analytical assessments in this extremely low range of mass loadings, the real extent of this burden remains to be determined in future researches. Therefore, these preliminary results need additional studies on the characterization of size distribution from aircraft aerosol emissions. These goals are the main objectives to continue further studies from aircraft aerosol emissions.


## 4. Conclusions

We outlined the scope of a sonication/agitation extraction method with minimal handling and reagents for PAC determination in PM samples with low mass loadings (<10 mg). Quality assessment of results was investigated from analysis of real samples. Therefore, we can conclude the following:(i)
The proposed method has been tested by comparison with extractions performed by microwave and pressurized solvent from real samples with PM mass loadings of 10 mg. In general, ratios between results obtained from the proposed method and PLE were above 80% for the most abundant PAHs and dichloromethane as solvent. Determination of selected NPAHs requires the use of methanol as modifier in extraction solvent to improve extraction efficiency.(ii)
Evaluation of intermediate precision was carried out by using two practical approaches. Both were in agreement. From analyses of samples with mass loadings between 0.7 and 8 mg, precision of around 10% for the majority of investigated PAHs was found. Precision close to 20% was found in the case of the most volatile PAHs and NPAHs.(iii)
The effectiveness of the method was evaluated to analyze samples with both very low mass loadings (<0.2 mg) and concentration levels (a few ng in mass), that is, impactor substrates used to determine the nano- and submicron size distribution.



---

*Source: 3517457-2017-07-27.xml*