key: cord-0015403-p1ed577r
authors: Castro-Guijarro, Pablo Antonio; Álvarez-Vázquez, Eusebio Ramón; Fernández-Espinosa, Antonio José
title: A rapid Soxhlet and mini-SPE method for analysis of polycyclic aromatic hydrocarbons in atmospheric particles
date: 2021-02-16
journal: Anal Bioanal Chem
DOI: 10.1007/s00216-021-03188-9
sha: d76e399e019ff103bb9fe0e7fdd62516d6a219c5
doc_id: 15403
cord_uid: p1ed577r

An analytical method was validated with two reference materials of polycyclic aromatic hydrocarbons in atmospheric particles. Standard reference materials (SRMs) were incorporated into the matrix of unexposed cut quartz filters. The methodology was previously designed and extraction of polycyclic aromatic hydrocarbons (PAHs) from fortified filters was based on a rapid low-cost method, for a low consumption of volume and time. The optimisation combined a low-volume Soxhlet apparatus used in hot Soxhlet mode with a quick clean-up by solid-phase extraction with special cartridges. The quantification of target compounds was performed by gas chromatography/mass spectroscopy in SIM mode. Temperatures of injector and oven program of the GC-MS were also optimised. Experimental variables of both systems were successfully optimised and validated, achieving a robust analytical methodology. [Image: see text] SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00216-021-03188-9.

Polycyclic aromatic hydrocarbons (PAHs) are lipophilic compounds with important carcinogenic risk. These hazardous hydrophobic organic compounds (HOCs) led international organisations to establish policies for the environment [1] and indoor air [2] . The European Commission exposed a wide toxicological guidance since the 2001 position paper [3] to the Europe 2020 strategy [4] . Similar measures were taken by other European organizations [5] . In the atmospheric environment, PAH emissions are well known, as they originated from car exhausts [6] , domestic combustion [7] , industrial activities [8] , agriculture activities [9] and natural sources [10] , such as biomass burning [11] . Determination of low levels of PAHs in complex matrices such as atmospheric particles leads to a search for accurate analytical methodologies, which need to be optimised and validated with recognized material references.

Traditional extraction procedures for PAHs involve many techniques, such as direct extraction with organic solvents, sonication and Soxhlet extraction [12] . In order to speed up sample preparation, new methodologies such as accelerated solvent extraction [13] or microwave-assisted solvent extraction [14] or ultrasound-assisted solvent extraction [15] have been developed, saving processing time and solvent consumption. However, these new techniques are expensive. Soxhlet extraction represents an inexpensive method for solid samples, nevertheless is rather time-consuming and requires high solvent consumption. An improved Soxhlet technique, called hot Soxhlet, heats the extractor body at lower temperature than the boiling point of solvents to keep it in the liquid state [16] . Additionally, after extraction step, it needed purification steps by solid-phase extraction techniques like packed columns of silica, Florisil®, alumina or mixtures [17] .

Therefore, the aim of this paper is to optimise and validate a low-cost and rapid extraction methodology with low time and solvent consumption. The method was based on an improved and more rapid Soxhlet extraction plus an improved, short and more efficient SPE technique for the 16 PAHs included in the US EPA priority pollutant list [18] in airborne particles.

Regarding the validation process, suitable certified reference materials should be used, but it is much better with two references [19] to confirm the suitability of the method.

Additionally, the reference material was mixed with small pieces of unexposed filter, taking the mixture as a new reference for the same matrix of real samples of particles, which are collected on filters during the sampling. This spiking technique constituted an innovative idea in validation techniques. The relevance of the work is the fact that many modern techniques that reduce solvent volume and time consumption are costly when there is currently a global crisis due to SARS-CoV-2.

All solvents used in the present study were of analytical and chromatographic grade.

Optimisation of the new procedure: extraction, purification and quantification

We validated in 2005 a five-stage analytical procedure for PAHs in airborne particles using the NIST 1649a (National Institute of Standard and Technology, Gaithersburg, MD, USA) [20] where extraction was performed by Soxhlet with 250 mL of solvent for 10 h. Also, we optimised in 2016 a fourstage procedure for fat in olive fruits [21] using 100 mL of solvent for 2.5 h. The proposed methodology reduced the method to only two stages before GC quantification, reducing time of extraction and volume consumption. Thus, optimisation of this new method was structured as follows:

Optimisation of PAH extraction was done by intercomparing three experiments using three Soxhlet apparatus, with solvent volumes of 250 mL [21] , 100 mL [18] and 25 mL. The lower the volume of the extractor body, the lower the time of extraction.

For these experiments, 100 mg of standard reference materials (SRM) 1649a was added to each extractor body over onequarter of unexposed QM/A quartz filters of 4 × 5 in. (Whatman International, Maidstone, England) cut in small pieces of 0.25 cm 2 (0.04 in. 2 ). Previously, quartz filters were heated in a muffle furnace at 500°C for 2 h to remove residual organic traces. The mixture employed was a 5:1 acetone/ methylene chloride (Merck, Darmstadt, Germany) mixture including pyrene-d10 (Dr. Ehrenstorfer GmbH, Augsburg, Germany) as deuterated surrogate standard. Volumes of organic extracts (around 250 mL and 100 mL) were reduced by rotary evaporation and then until 2 mL inside the test tubes by slow nitrogen concentration. Purification and quantification stages were carried out for the time being as in the 2005 validation [21] . Extracts of 25 mL extractor body were nitrogen-concentrated directly. After this volume optimisation, the time of extraction was then minimized on the best Soxhlet system obtained from 250, 100 and 25 mL.

The process of analyte isolation requires a sample clean-up stage to remove interferences. The extracts of filter from the extraction stage were purified comparing the OCC technique (open-column adsorption liquid chromatography) and the SPE technique in normal phase. In the OCC, alumina/silica (from Merck) was used as adsorbent [21] , and in the SPE, miniSpe-ed Plus silica gel cartridges (Applied Separations, Allentown, PA, USA) were proposed. SPE allows isolating PAHs from a sample through the same chemical principles of column chromatography but with low consumption of solvent and quickness.

The OCC technique was performed in glass columns (30 cm × 1 cm) filled with 1 g of each alumina/silica gel (top/bottom), conditioning it with 30 mL of n-hexane (from Merck). The extract (2 mL) was transferred to the top and the non-polar fraction was eluted with other 2 mL n-hexane, and 4 mL was discarded. The aromatic fraction was eluted with 7 mL of 20% methylene chloride in n-hexane and nitrogen concentration of < 1 mL-adding here the deuterated standard mixture-and making up to 2 mL into a chromatographic vial.

SPE technique with miniSpe-ed Plus silica gel cartridges was performed in a Varian vacuum manifold (Varian Inc., Scientific Instruments; Palo Alto, CA, USA). Each cartridge (450 mg/1 mL) was conditioned with 1.5 mL of n-hexane and 1.5 mL methylene chloride. Then, for the aliphatic fraction, 2 mL of extract plus a few microlitres of deuterated internal standards (ISs mix) were loaded and eluted with 3 mL of nhexane at a flow of 0.3-0.5 mL min −1 (higher elution speeds lead to low retention volume), and then it was discarded. Then, 3 mL of methylene chloride was used for eluting PAHs into a 2-mL vial. The chromatographic vial was located inside the vacuum manifold assisting it with a flow of nitrogen. So, a slow nitrogen concentration stage was not necessary.

The sixteen PAHs listed as priority air pollutants by the US EPA [18] were identified and quantified using gas chromatography (Agilent Series 6890A, Santa Clara, CA, USA) with mass spectrometry (Agilent Series 5973N). A capillary column of (5%-phenyl)-methylpolysiloxane (low polarity) was used (HP-5ms, 30 m × 0.25 mm i.d. × 0.25 μm film thickness) from Agilent Technologies. Helium (99.9995% purity) as carrier gas was operated at constant pressure.

Chromatograms obtained at different temperatures of the injector and different optimised oven programs were compared. These temperatures were systematically optimised by using an orthogonal design approach. First, the inlet temperature was tested at 250, 260, 270 and 280°C with the rest of the instrumental conditions fixed. Second, the oven temperature program was optimised in order to obtain the best resolution and separation of chromatographic peaks. Three temperature programs were tested according the following conditions:

PrA: The initial temperature of 60°C was kept for 1 min, then it rose to 210°C at 15°C min −1 , was kept for 1 min, and finally to 280°C at 15°C min −1 , kept for 25 min. Total time: 41.7 min. PrB: The initial temperature of 60°C was kept for 2 min, then it rose to 200°C at 7°C min −1 , was kept for 2 min, and finally to 290°C at 10°C min −1 , kept for 35 min. Total time: 68.0 min. PrC: The initial temperature of 60°C was kept for 1 min, then it rose to 175°C at 20°C min −1 , was kept for 3 min, then it rose to 300°C at 5°C min −1 , kept for 20 min. Total time: 54.8 min.

Once temperatures were optimised, the instrument quantification method [22] required a calibration curve 4-1000 μg L −1 [ppb] for each. The chromatographic signal of each PAH was relative to a deuterated PAH of the acenaphtene-d 10 , phenanthrene-d 10 , chrysene-d 12 and perylene-d 10 internal standard mixture at 200 ppb (from Dr. Ehrenstorfer). Each peak was identified by the absolute and the relative retention times, and by comparison with the mass spectral library of the instrument [23] using a target ion, primary ion (T) and a qualifier molecular ion (Q) (see Supplementary Information (ESM), Table S1 ).

In order to study the accuracy (trueness and precision according ISO [24] ), two NIST-certified reference materials for PAHs were used, to cover two levels of concentrations: SRM 1649a-Urban Dust and SRM 1648a-Urban Particulate Matter (from NIST) [19, 25] . These studies were developed on SRMs 'with/without' filters. Studies to demonstrate the accuracy of the method included (see ESM) (a) recovery study, (b) t tests and F assays, (c) precision study by Horwitz ratios (HorRat), (d) linearity study, (e) sensitivity, (f) selectivity/ specificity determinations, (g) limits of detection and quantification, (h) ruggedness.

For (a), (b) and (c) studies, twelve replicates of 100 mg of the two SRMs were analysed. Other nine calibration series, which employed PAH standards and internal standards, constituted the base of the (d) to (g) validation studies. All experiments of ruggedness were performed with 100 mg of SRM. All statistics were done in agreement with guidelines for validation of the AOAC (Association of Official Analytical Chemists) [26] , the IUPAC (International Union of Pure and Applied Chemistry) [27] and the SANCO-DG (Directorate General for Health and Consumer Affairs at the European Commission) [28] .

(a) Recovery study: Recoveries were calculated for n = 12

replicates analysed in different days and weeks according Eq. (1; see ESM). Each recovery value was compared with the AOAC acceptable recoveries. (b) F assays and t tests: Hypothesis F tests for the precision and t tests for the trueness were done at p = 0.05 for the n = 12 replicates by comparing variances and mean values of the measurements with those certified by the NIST SRMs according Eqs. (2, 3 and 4; see ESM). (c) Precision study: The study of the intra-laboratory precision was done under reproducibility conditions (RSD R ), according to twelve replicates' experiment based on two daily sessions with duplicates (j = 1 and j = 2, morning/ afternoon sessions), for n = 1, 2 and 3 days of different weeks (between-days). The expected AOAC values of RSD for reproducibility (ERSD R ), the predicted RSD R values (PRSD R ) and the Horwitz ratio, HorRat, see Eq. (5; see ESM), were used as indicator of precision according the Horwitz or Thompson theories [29, 30] . (d) Linearity R 2 , r, L and CV: The linearity study was performed from data of different calibration curves during the validation processes. Coefficient of determination (R 2 ), the Pearson coefficient of correlation (r), the 'goodness' (t-significance), the percentage (L) of linearity [31] and the 'online linearity' (CV) were calculated according to Eqs. (6 to 10; see ESM). (e) Sensitivity γ and δ: The PAH sensitivity was determined as the 'sensitivity of calibration' (gamma, γ) and as the 'analytical sensitivity' (delta, δ). γ is defined as the slope of the regression curve (Eq. 11; see ESM). δ is the ratio of calibration sensitivity to standard deviation of the slope (Eq. 12; see ESM). (f) Selectivity: Interferences can be introduced through the sample matrix, the sampling system or the instrument system. To reduce interferences, we evaluated the use of pure solvents, the sample handling with laminar airflow cabinets and the use of blanks. Specifically, cleanup procedures were employed to remove most of these substances. In the case of co-eluting compounds, the mass spectrum can be easier interpreted working in SIM mode. (g) Limits of detection and quantification, criteria comparison: Limits of detection were essentials and they are evaluated by various criteria. First, both limits were determined based on 1978-IUPAC criteria (Eq. 13; see ESM) from the mean of procedural blanks (LOD Bl , ×3

and LOQ Bl , ×10, n = 8), i.e. matrices containing only 1 μL of 50 pg μL −1 of the IS mix [26] .

Second, limits of detection and quantification were also estimated as instrument detection limit (LOD I , ×3 and LOQ I , ×10, n = 12) from the standard deviations (Eq. 14; see ESM) of calibration curves (ISO 11843, [32] ).

For chromatography, standards with concentration close to LOD I were required (n = 12). It is the noise detection limit, LOD N (×3) and LOQ N (×10) prepared from spiking blank samples: we spiked unexposed cut quartz filters with 1 μL of 100 pg μL −1 (100 ng of the 16 PAH mix standard and 50 ng of the IS mix). Calculation was from the signal-tonoise ratio (S/N) [33] according to Eq. 15 (see ESM).

Also, the limits were determined according to the method used by the ISO-17025-accredited laboratories. The limit of quantification was estimated from a standard at a near-zero concentration (NZC). Thus, three calibration curves were studied by quadruple at low concentrations < 25 μg L −1 [ppb]; preparing at high range (HR), 5-25 pg μL −1 ; medium range (MR), 0.5-2.5 pg μL −1 ; and low range (LR), 0.05-0.25 pg μL −1 . Then, the standard with the lowest concentration (determined with less than 3% error) is selected as the reference for the limit of quantification (LOQ ZC ). The limit of detection (LOD ZC ) was estimated as one-half of the quantification limit. Finally, we expressed the method detection limits (MDLs) and method quantification limits (MQLs) as those expressed in pg m −3 and ng g −1 .

(h) Ruggedness study: The method performance was evaluated using a 'ruggedness test' [34, 35] by the Youden method. It was based on the Plackett-Burman orthogonal design of eight experiments and seven factors (L8, 2 ). They included small changes in seven (A-G) chromatographic factors at two levels for each factor (overdefault or high-low) with regard to nominal conditions, affecting to the four following experimental variables: Ruggedness was determined by triplicate with 100 mg of both SRMs according the Hadamard matrix. We then evaluated the variables that best and worst adapted to the small 'accidental' changes.

The results for searching the best procedure are shown below. Statistical techniques were ANOVA for intercomparison of the three experiments of extraction, Snedecor-Fisher (F test) for precision and Student's (t test) for accuracy and intercomparison of the two purification techniques.

In the intercomparison of the three Soxhlet systems, the extraction time reached 10 h for the 250 mL extractor body, and 2.5 h for those of 100 mL. In the case of extraction at 25 mL, a special mini-Soxhlet system was required: the Quickfit® Soxhlet extractor body with only 20 mL of siphoning volume. The time of extraction was also optimised here. Procedure for 250 mL, 100 mL and 25 mL was satisfactorily compared (n = 5) of NIST SRM 1649a.

Firstly, F and t tests (p = 0.05) applied to both 250 mL and 100 mL series showed no significant difference between them for the 13 PAHs of the SRM (all F calc < 1.37 for F crit = 9.60 and all t calc < 1.02 for t crit = 2.31, p > 0.05). However, F and t tests applied to both 250 mL and 25 mL series and to both 100 mL and 25 mL series showed significant differences between 250 mL or 100 mL and 25 mL series (all F calc > 15.42 for F crit = 9.60 and t calc > 5.13 for t crit = 2.31, p < 0.05). Secondly, ANOVA showed significant differences along the three series of solvent volumes (F calc = 13.31 > F crit = 3.26, p = 0.008). In addition, recoveries of the 25 mL procedure were higher (+ 4.3-4.8%) than those of 250 mL and 100 mL (t calc < 0.96 for t crit = 2.07, p > 0.05) (Fig. 1) . As a result, the Quickfit® Soxhlet reduces the required volume of solvent to 25 mL and it allows a high number of cycles (67-83) and shorter cycles (65-80 s), reducing the total time of extraction to 60 min (see ESM, Fig. S1 ).

The low volume (24-25 mL) of the final extract implies a short duration of the NASE (nitrogen-assisted solvent evaporation), which, additionally, were done inside the vacuum manifold. On the other hand, to further increase the efficiency of the rapid extraction, this was performed under the 'hot Soxhlet' conditions [16] .

After extraction, 24-25 mL of the extract was concentrated to 2 mL inside the vacuum manifold before the purification stage. In the clean-up optimisation, both OCC and SPE series were compared (n = 5) using SRM 1649a (Fig. 2) , showing no significant difference between them (all F calc < 2.06 for F crit = 9.60 and t calc < 1.18 for t calc = 2.31, p > 0.05). However, recoveries obtained with SPE were higher (+ 2.0%) than for OCC (t calc < 0.81 for t crit = 2.07, p > 0.05).

In conclusion, according to the best recoveries obtained with Quickfit® Soxhlet and SPE, the methodology proposed combines the accelerated mini-Soxhlet extraction (AmSE) assisted by hot Soxhlet and the mini-solid-phase extraction with miniSpe-ed (mSPE) assisted with simultaneous nitrogen evaporation, as the rapid Soxhlet and solid-phase extraction (RSE-mSPE) method.

The final RSE-mSPE method was performed in two steps: First, accelerated mini-Soxhlet extraction on atmospheric filters for 60 min with 25 mL of solvent, using the Quickfit® Soxhlet, and second, a mini-solid-phase extraction with silica gel miniSpe-ed on filter extracts for the PAH isolation before GC injection.

Measurements in the chromatographic system were done using helium at a constant pressure of 20 psi, operating in pulsed splitless mode, and injecting 2 μL of all the samples, split opened after 30 s. Optimisation of the inlet temperature was tested at 250, 260, 270 and 280°C, fixing the other parameters. All PAHs of 1649a SRM were observed, studying especially the representative compound, benzo[a]pyrene [BaP] as the first marker of total PAHs in scientific reports [36] . When the injector temperature increases from 250 to 280°C (see ESM, Fig. S2 ), the maximum peak area of the BaP was given to 260°C; therefore, the temperature of 260°C was adopted as the optimal temperature of the injector. As for the PrA, PrB and PrC oven temperature programs, the three were obtained with good resolution between PAHs. The responses of chromatographic peaks were similar for the three programs tested. However, because the second is too long (68 min) and although the first is shorter (42 min), the third program has slightly higher peak resolution for some PAHs than the other two, so the oven temperature program selected was PrC (55 min): 60°C, 1 min, 20°C min −1 to 175°C, 3 min, 5°C min −1 to 300°C and hold for 20 min.

Finally, the mass selective detector was operated in electron ionization mode with electron energies of 70 eV, being the ion source temperature of 230°C. To improve sensitivity, quantitative analysis operated in selected ion monitoring (SIM) mode instead of full scan mode. The selected molecular ions of the different PAHs were shown in Table S1 (see ESM).

Recovery and RSD values obtained for both SRMs (n = 12) were shown in Table 1 . Average recovery for SRM1649a + filter was Rec = 97.9% (82.6-107.2%) and for SRM1648a + filter was Rec = 95.5% (75.7-104.8%), with RSD = 2.0% for SRM 1649a and 1.5% for SRM 1648a. All recoveries were 'acceptable' according to the AOAC ranges, showing excellent agreement between measured and certified values. Differences between experiences with/without pieces of filters were negligible, 97.8% without filters (SRM only) and 98.9% with filters (as real samples) for SRM 1649a and a similar difference for SRM 1648a.

The aromatic compound with the lowest recovery in SRM 1649a was fluorene (3 rings) with a value of 82.6% [37] . This was probably due to its volatility (MW = 166 uma and Pv = 0.09 Pa). Similar results in SRM 1648a were for acenaphthene and acenaphthylene (76%, 3 rings, MW = 152-154 uma and Pv = 0.3-0.9 Pa) and naphthalene (82%, 2 rings, MW = 128 uma and Pv = 8.6). On the other hand, the recovery of the extraction method from pyrene-d 10 was high with a value of 98.1 ± 2.8%.

Results of F tests ( Table 2) revealed that for all PAHs of both SRMs, differences between variances were not significant (F crit > F calc , p = 0.05), except for indeno[1,2,3-cd]pyrene and benzo[ghi]perylene in SRM 1649a; also, RSD R values for naphthalene and dibenzo[ah]anthracene were higher than those acceptable by the AOAC. Afterwards, t tests showed good results, but not for fluorene in SRM 1649a and acenaphthene, acenaphthylene, benzo[b]fluoranthene and dibenzo[a,h]anthracene in SRM 1648a, which showed significant differences against true values (t crit < t calc , p = 0.05). As a result, from the recoveries (a) and the t tests and F assays (b), we can affirm that the proposed new analytical methodology was traceable to both NIST SRMs without any important systematic error.

The intra-laboratory study ( Table 2) showed that values of RSD R were lower than AOAC-ERSD R values in all PAHs of SRM 1649a, so they were acceptable and similar than other works [38] . This fact occurred precisely for the two PAHs whose certified RSD values were high, indicating good precision for repeatability but not for reproducibility.

The linearity of the different calibration curves (4-1000 ppb) can be represented by the goodness or t significance (t calc ) and by other parameters (r, R 2 and L). In Table S2 (see ESM) , the parameters showed excellent results at p = 0.05. Coefficients of correlation r were over 0.997 (0.997-0.9997) and significantly different from zero (t calc > t crit ) and greater than the critical value 0.707 for a bad linearity. Coefficients of determination R 2 were over 0.993 (0.993-0.9994) [37] . Linearity L(%) were over 95% (96.1-98.9%) and coefficients of variation of the slope CV b (%) were less than 5% (1.12-3.91%).

The values obtained for gamma γ sensitivity (ESM Table S2 ) ranged from the less sensitive, such as 2. The selectivity study is necessary in complex matrices, such as filters of airborne particles; thus, it must reduce interferences from real samples and contamination from blanks. Current analysis used selective chemicals that removed interferences and contamination, such as specific cartridges (miniSpe-ed) and solvents of chromatographic grade, and selective instruments, such as gas chromatograph with capillary columns, and a mass spectrometry detector. In the case of coeluting compounds, the MS detector was set in the selected ion monitoring mode (SIM). The clean-up procedure used with miniSpe-ed removes aliphatic hydrocarbons and polar compounds. Besides, analysis of blanks also proved that analytical determinations were free from contaminants. In addition, possible contaminants were reduced by handling all formats of samples inside the laminar air-flow cabinet INDELAB®, Model IDL-48 V with a HEPA filter plus a charcoal layer, which ensure a clean air inside the cabinet. Consequently, the methodology resulted free from chemical interferences or at least with interferences controlled.

Limits of detection from blanks (LOD Bl ) were lower (4-7 times) than the instrument detection limits (LOD I ) and the (Table 3 ) due to the purity of the chemicals used. Therefore, a more realistic alternative is a detection limit based on the signal-to-noise ratio (3 × S/N). As a result, values of the noise detection limits (LOD N ) were similar to the instrumental detection limits, that is LOD Bl << LOD I ≅ LOD N . Moreover, as the HR and MR ranges gave good linearity and all NZC standards were determined with accuracy (< 3% error), we focused on the low range (LR). Thus, at the lowest concentration obtained, we injected 12 replicates of the reference and quantification and detection limits (LOD ZC , onethird of LOQ ZC ) were estimated.

In brief, final comparison resulted in LOD Bl < LOD ZC < LOD N ≅ LOD I and the selected method detection limit (MDL) was the instrumental detection limit (LOD I ), expressed as pg μL −1 or as pg m −3 and ng g −1 . Values of Table 3 were lower than other studies [38] . However, Piñeiro-Iglesias et al. [19] reported limits similar for anthracene but not for fluoranthene and pyrene, which were higher than those reported here.

The ruggedness of a procedure must be established for 'inhouse' developed methods [33] . Seven factors (A-G) were tested along four chromatographic conditions: The factor variations are shown in Fig. 3 . The deviations of chromatographic parameters were below 12.5% for negative deviations and 5% for positive ones. Therefore, SRM recoveries were always within the range of AOAC percentages of -: Values of LOD Bl could not be estimated due to the absence of the compounds in the procedural blanks

Estimations for instrument detection limit (LOD I , n = 12), blank detection limit (LOD Bl , n = 8), noise detection limit (LOD N , n = 12) zero concentration detection limit (LOD ZC , n = 12) and method detection limit (MDL) expressed as concentration in the ambient air (MDL A , pg m −3 ) and in the solid particulate matter (MDL B , ng g −1 ) recovery. The only factor that most affects results is factor A (Soxhlet extraction time) when is less than 60 min. Therefore, there were no significant differences in performance of the method as a result of the deliberate changes, implying that our developed method is robust in terms of reliability. The complete optimised and validated method was (Fig. 4 ) 100 mg of one-half cut quartz filter of airborne particles was extracted in the AmSE system with 25 mL of 5:1 acetone/ methylene chloride for 60 min. After quick evaporation until 2 mL with NASE, the organic extract was purified by mSPE using a miniSpe-ed cartridge, eluting PAH simultaneously to nitrogen evaporation with 3 mL of methylene chloride.

The proposed RSE-mSPE method for extraction of PAHs in atmospheric particles helps in minimizing solvent volume and time extraction. The two main stages optimised AmSE+mSPE and improved the recovery and accuracy of PAH determination. Accelerated mini-Soxhlet extraction system (AmSE) reduced volume consumption in 90% (10 times), time consumption in 90% (10 times) and increased recoveries of SRM in 9%. Mini-solid-phase extraction (mSPE) also reduced time consumption in 90% (9 times) and increased recoveries of SRM in 3%. The other parameters of validation were successfully assessed with good agreement between certified and reference values, thereby no interferences being by the presence of quartz filters. Thus, it can be stated that the proposed RSE-mSPE methodology is a simple, fast and lowcost method that can be used, combining with GC/MS-SIM, for the determination of PAHs in airborne particulate matter.

The online version contains supplementary material available at https://doi.org/10.1007/s00216-021-03188-9.

Funding This work would not have been possible without the partial financial assistance of the Environmental Agency of the 'Junta de Andalucía', the Regional Government of Southern Spain, through the different Research Projects on Air Quality provided to our Research Group RNM-294 from 1995 to 2015. 

Regional Office for Europe. WHO regional publications, European series

Polynuclear aromatic hydrocarbons. In: WHO Guidelines for indoor air quality. Selected Pollutants

Ambient air pollution by polycyclic aromatic hydrocarbons (PAH). Position Paper

AUTH Chemical group leader. HBM4EU Project of the European Union's Horizon 2020 programme: Science and policy for a health future

Background document on polycyclic aromatic hydrocarbons (PAHs), update 2009. Hazardous Substances Series. The Executive Secretary

Environmental science for European refining. Division of the European Petroleum Refiners Association

Best environmental practice (BEP) for the reduction or prevention of emissions of polycyclic aromatic hydrocarbons (PAHs) from domestic combustion appliances. OSPAR Background Document. The Executive Secretary

Spatially high-resolved monitoring and risk assessment of polycyclic aromatic hydrocarbons in an industrial city

Sources, distribution and toxicity of polyaromatic hydrocarbons (PAHs) in particulate matter, air pollution, Vanda Villanyi

Sources, concentrations, fate and effects of polycyclic aromatic hydrocarbons (PAHs) in the environment, review. Part a: PAHs in air

Effect of northern boreal forest fires on PAH fluctuations across the artic

Trends in analytical techniques applied to particulate matter characterization: a critical review of fundaments and applications

Tenax-TA spiking approach of thermal desorption coupled to GC-MSMS for the quantification of PAHs in indoor air and dust

Rapid determination of 16 polycyclic aromatic hydrocarbons in PM2.5 by microwave assisted extraction-high performance liquid chromatography

Ultrasoundassisted solvent extraction (UASE) of total polycyclic aromatic hydrocarbons from mussels followed by spectrofluorimetric determination

Comparison of hot Soxhlet and accelerated solvent extractions with microwave and supercritical fluid extractions for the determination of polycyclic aromatic hydrocarbons and nitrated derivatives strongly adsorbed on soot collected inside a diesel particulate filter

H Oxygenated polycyclic aromatic hydrocarbons in atmospheric particulate matter: molecular characterization and occurrence

Priority pollutant list. From toxic and priority pollutants under the clean water act

Microwave assisted extraction of polycyclic aromatic hydrocarbons from atmospheric particulate samples

Particlesize distribution of polycyclic aromatic hydrocarbons in the urban air in southern Spain

Combining PLS regression with portable NIR spectroscopy to on-line monitor quality parameters in intact olives for determining optimal harvesting time

Method 8270E (SW-846): Semivolatile organic compounds by gas chromatography/ mass spectrometry (GC/MS), rev

Mass Spectrometry Data, Mass Spectrometry Data Center

Accuracy (trueness and precision) of measurement methods and results -part 1: General principles and definitions

Optimisation and validation of a low temperature microwave-assisted extraction method for analysis of polycyclic aromatic hydrocarbons in airborne particulate matter

Guidelines for standard method performance requirements. AOAC official methods of analysis

Harmonized guidelines for single-laboratory validation of methods of analysis

Technical active substance and plant protection products: guidance for generating and reporting methods of analysis in support of preand post-registration data requirements for Annex II (Section 4) of Regulation (EU) No 283/2013 and Annex (Section 5) of Regulation (EU) No 284/2013. Guidance document. Directorate General Health and Consumer Protection

The Horwitz function revisited

Precision parameters of methods of analysis required for nutrition labeling, part I. major nutrients

Statistical estimation of linear calibration range

ISO11843-2:2007 capability of detection-part 2: methodology in the linear calibration case

Title 40. Protection of environment. Chapter I. EPA. Subchapter D. Water programs. Part 136 -guidelines establishing test procedures for the analysis of pollutants. Appendix B to Part 136-Definition and procedure for the determination of the method detection limit

Eurachem guide: the fitness for purpose of analytical methods-a laboratory guide to method validation and related topics

Statistical manual of the AOAC-statistical techniques for collaborative tests (by)

Cancer risk assessment, indicators and guidelines for polycyclic aromatic hydrocarbons in the ambient air

Associations between inhaled doses of PM2.5-bound polycyclic aromatic hydrocarbons and fractional exhaled nitric oxide

Phase distribution of polycyclic aromatic hydrocarbons and their oxygenated and nitrated derivatives in the ambient air of a Brazilian urban area

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

Conflict of interest The authors declare no conflict of interest.