A Rapid Preconcentration Method Using Modified GP-MSE for Sensitive Determination of Trace Semivolatile Organic Pollutants in the Gas Phase of Ambient Air
A Rapid Preconcentration Method Using Modified GP-MSE for Sensitive Determination of Trace Semivolatile Organic Pollutants in the Gas Phase of Ambient Air
Bulletin of the Korean Chemical Society. 2014. Oct, 35(10): 2995-3000
Copyright © 2014, Korea Chemical Society
  • Received : March 13, 2014
  • Accepted : June 16, 2014
  • Published : October 20, 2014
Export by style
Cited by
About the Authors
Miao, He
Qingjuan, Xu
Cui, Yang
Xiangfan, Piao
Narayanan, Kannan
Faculty of Environmental Studies, Universiti Putra Malaysia, 43400 UPM Serdang, Selangor Darul Ehsan, Malaysia
Donghao, Li

A sensitive concentration method utilising modified gas-purge microsyringe extraction (GP-MSE) was developed. Concentration (reduction in volume) to a microlitre volume was achieved. PAHs were utilised as semivolatile analytes to optimise the various parameters that affect the concentration efficiency. The injection rate and temperature were the key factors that affected the concentration efficiency. An efficient concentration (75.0-96.1%) of PAHs was obtained under the optimised conditions. The method exhibited good reproducibility (RSD values that ranged from 1.5 to 9.0%). The GP-MSE concentration method enhances the volume reduction (concentration factor), leading to a low method detection limit (0.5-15 ng L -1 ). Furthermore, this method offers the advantage of small-volume sampling, enabling even the detection of diurnal hourly changes in the concentration of PAHs in ambient air. Utilising this method in combination with GC-MS, the diurnal hourly flux of PAHs from the gas phase of ambient air was measured. Indeed, the proposed technique is a simple, fast, low-cost and environmentally friendly.
The development of sample pretreatment techniques for environmental analysis has been the subject of much researchers. 1-3 Although certain organic pollutants in the environment occur at trace levels, they can seriously threaten human health, either directly or indirectly, due to their toxicity and bioaccumulation. Therefore, the detection of trace environmental organic pollutants is necessary. The analysis of trace levels of organic pollutants in the environment by chromatographic techniques presents a great challenge due to their low detectability by instrumental analysis. Additionally, matrix interference can be a problem as well. 4-6 Conventional pretreatment techniques require large sample sizes/volumes, which in turn require large amounts of time and energy to process. The detection of PAHs in the gas phase in the atmosphere represents a classical problem in this context. 7-9 The analysis of trace semivolatile compounds usually requires the extraction of analytes from a sample matrix, followed by concentration (volume reduction) prior to analysis for better detectability in chromatographic techniques. 10,11 Many studies have been dedicated to the improvement of the extraction and cleanup of these compounds, but only a few have focused on the improvement of the concentration of analytes. 12,13 It is very important to enhance the concentration factor because doing so will improve the detectability of target compounds, allowing for small sample sizes and short analytical times.
Conventional concentration methods are based on volume reduction by the evaporation of the solvent that contains the analytes. Inert gas streams are routinely used in laboratories. 14,15 A common procedure is to blow an inert gas stream onto the surface of the solvent (extract). The more volatile solvents are removed, and the less volatile analytes are retained, which concentrates them. The rate of evaporation may be increased by heating, which accelerates the concentration process. Rotary evaporators are also used routinely in laboratories, utilising a water bath as the heat source. 16,17 These evaporation techniques can only achieve millilitre-level volume reductions; thus, further requires microconcentration is required to suit instrumental analysis.
In recent years, many reports on new concentration techniques such as the hollow-fibre membrane pervaporation technique have been published. 18,19 Pervaporation is an effective separation method utilising membrane permeation and evaporation. 20 The hollow-fibre membrane pervaporation technique shows a low efficiency for semivolatile compounds and could be “environmentally unfriendly” because the volatilised solvent escapes to the atmosphere. A purge-andtrap concentrator has also been commercialised, having been pioneered by Bellar and Lichtenberg in 1974. Many researchers have developed several methods based on the purge-and-trap technique. 21,22 Despite the method’s sensitivity and reliability, the operation of these methods is complex, and serious foaming and clogging of the gas easily occur. Moreover, the methods are more suitable for the detection of volatile compounds. 23,24
Gas-purge microsyringe extraction (GP-MSE) is a novel sample pretreatment technology developed by Yang et al ., 25 that has been patented. 26 The features of this technique are integration and miniaturisation. GP-MSE represents the most effective procedure for the extraction and preconcentration of analytes. This technology was developed from gasflow headspace liquid-phase microextraction (GF-HS-LPME), which is based on liquid-phase microextraction (LPME). 27 GP-MSE is an exhaustive extraction technique, and it is much faster (3 min) than existing methods, environmentally friendly, simple to operate, economical and demonstrated to have high extraction efficiency for volatile compounds and semivolatile chemicals. 27
An attempt was made in this study to develop a sensitive concentration method modelled after GP-MSE that can concentrate extracts to the microlitre level and greatly reduce the required sample volume. PAHs were selected as representative semivolatile analytes, and various parameters that affect the concentration efficiency, such as the temperature, the sample-injection rate and the nitrogen-flow rate, were optimised. Finally, an attempt was made to measure the diurnal hourly flux of PAHs in the gas phase of ambient air using this method.
Experimental Section
Material and Methods. The standard mixture used in this study contained sixteen PAHs: naphthalene (Nap), acenaphthylene (Acp), acenaphthene (Ace), fluorine (Flu), phenanthrene (Phe), anthracene (Anth), fluoranthene (Fluo), pyrene (Pyr), benzo[ a ]anthracene (B[ a ]A), chrysene (Chr), benzo[ b ]fluoranthene (B[ b ]F), benzo[ k ]fluoranthene (B[ k ]F), benzo[ a ]pyrene (B[ a ]P), indeno[1,2,3- cd ]pyrene (InP), dibenzo[ a,h ]anthracene (DBA), and benzo[ ghi ]perylene (BghiP). An internal standard ([ 2 H 10 ] phenanthrene) and a surrogate standard ([ 2 H 10 ] acenaphthene and [ 2 H 12 ] perylene) were purchased from Supelco (Bellefonte, PA, USA). The purity of the standards was generally higher than 99%. HPLC-grade dichloromethane was obtained from Caledon (Georgetown, Ont., Canada). Stock standard solutions of the PAHs (20 mg L –1 ) were prepared in methanol and were diluted to 20 μg L –1 with dichloromethane. Standard working solutions of various concentrations were prepared by diluting the stock solutions with dichloromethane. The internal standard was spiked into dichloromethane (400 μg L −1 ) when it was used as an extracting solvent. The surrogate standard was also spiked into dichloromethane (200 μg L –1 ). The standard solutions and the extracting solvent were stored in the dark at 0−4 °C until use. The purity of the nitrogen (N 2 ) used as a flow gas was 99.999%. A single-channel syringe pump (KDS100), made in the USA, was purchased from Shenzhen.
Sample Collection and Preparation. The level of PAHs in the gas phase in ambient air was determined by collecting air in the gas-flow liquid-phase enrichment apparatus (GFLPE) developed previously by our group (see Figure S1 in the supplementary information). The sampling site was Yanji City, where Yanbian University is located. The air sampling rate was 0.530 L min –1 , and the sampling time was 4 h. The gas was pumped into the gas absorber through a flow meter, and the organic compounds were absorbed in 100 mL dichloromethane. After sampling, the surrogate standard was spiked into the sample (20 ng L –1 ) and then concentrated to 1 mL by conventional techniques and further tested by GPMSE. These experiments have already been standardised in our laboratory.
PPT Slide
Lager Image
Schematic diagram of the modified GPMSE apparatus. A(1), A(2), A(3) represent the concentration process; B(1), B(2), B(3) represent the absorption process. ● denotes the target compound; the grey slash represents the organic solvent; the dotted line represents the gas-flow pathway; and the wavy curve represents glass wool. The numbers indicate the following: 1, syringe pump; 2, microsyringe (1 mL); 3, PTFE tube; 4, needle; 5, sample cell; 6, heat chamber; 7, T-value; 8, gas-flow controller; 9, digital monitor; 10, waste vial; 11, cooler; and 12, microsyringe barrel (100 μL).
GP-MSE Concentration Apparatus and Procedure. Please refer to Yang et al . 25,27 for a description of the GPMSE concentration apparatus. It should by noted, however, that the version used in this experiment is slightly different from the original, as explained in Figure 1 . The sampleinjection port and the solvent recovery port were enlarged. The sample-injection port included a syringe pump, a microsyringe (1 mL) and a PTFE tube (O.D.: 1.5 mm; I.D.: 0.75 mm; length: 10 cm). The solvent-recovering section featured a waste vial and a PTFE tube. The PTFE tubes in the two sections were connected with a single microsyringe needle and matched to each other. A microsyringe (1 mL) was used to inject the sample into the sample cell by the syringe pump, and a waste vial was used to collect the volatilised solvent.
The concentrating process is explained in Figure 1(a). The process consists of the following steps. (1) A sample cell was placed into the heat chamber, and then the heat chamber was sealed with a PTFE septum. The two PTFE tubes, one connected to the microsyringe (1 mL) and the other connected to the waste vial, were carefully inserted into the chamber via the septum. (2) The temperature of the chamber and the nitrogen flow were set to 80 °C and 2.0 mL min –1 , respectively. The 0.5 mL sample solution (solvent and analytes) loaded in the microsyringe (1 mL) was continuously injected into the sample cell at an injection rate of 167 μL min –1 , and then the concentration process was initiated. Upon injection, the solvent and the analytes in the sample were separated due to their dissimilar volatilities. It is noted that the sample cell should always be moist during this process. At a relatively low temperature, the solvent was volatilised easily due to its low boiling point and was carried away by the nitrogen flow into the waste vial. The analytes were retained in the sample cell during this process due to their relatively low volatility. (3) When the sample was completely injected, all of the parameters of the GP-MSE were shut down. In this concentration procedure, only the analytes and a few microlitres of solvent remained in the sample cell, and almost all of the evaporated solvent was recovered in the waste vial.
The extraction process is explained in Figure 1(b) . The process consists of the following steps. (1) The PTFE tubes were pulled out, followed by the insertion of a 100 μL microsyringe barrel (without the plunger) into the chamber through the septum. (2) The GP-MSE parameters were reset as follows. The temperature of the heat chamber was 280 °C, the extraction time was 3 min, the nitrogen-flow rate was 2.0 mL min –1 and the cooler temperature was –4 °C. A 5 μL volume of extracting solvent was added to the 100 μL microsyringe barrel, and then the extraction was initiated. A similar operating principle and practical result were observed by Yang et al . 25-27 : briefly, at a high extraction temperature, the analytes and a few microlitres of solvent in the sample cell were volatilised and carried away by the nitrogen flow, and they were absorbed by extraction of the solvent in the microsyringe barrel. (3) After extraction (3 min), the analytes were quantitatively absorbed by the extracting solvent. The microsyringe barrel was removed from the apparatus, and the extract was directly analysed by GC−MS.
Analytical Methods. Chromatographic analysis was performed by a Shimadzu GC 2010 with a DB5 fused-silica capillary column (30 m × 0.25 mm; thickness: 0.25 μm) and detected on a Shimadzu QPMS 2010 mass spectrometer with a quadrupole mass-spectrometer system. The injection temperature was 280 °C in the splitless mode, and the interface temperature was set to 280 °C. The column temperature was maintained at 80 °C for 2 min, and then it was programmed to increase by 20 °C min –1 to 100 °C and then by 10 °C min –1 to 280 °C; it was then held at this latter temperature for 21 min. The column flow was 1.20 mL min –1 . Helium (99.999% pure) was used as a carrier gas. The electron ionisation selected ion mode (EI-SIM) was used for quantification of the analytes. The solvent-cut time was 5.0 min. The ion-source temperature and the energy of the ionising electrons were set to 200 °C and 70 eV, respectively.
Calculation. To optimise the various parameters, the preconcentration efficiency and concentration factor were calculated.
The preconcentration efficiency (CE) is defined as follows:
CE = n i /n 0 = C i V i /C 0 V 0 ,
where n i and n 0 are the analyte masses in the final extract after GP-MSE concentration and in the sample, respectively; C i and C 0 are the analyte concentrations in the final extract after GP-MSE concentration and in the sample, respectively; and V i and V 0 are the volume of the extract after GP-MSE concentration and that of the original sample, respectively.
The concentration factor (CF) was defined as the ratio of the concentration of the analytes in the final extract after GPMSE concentration to that in the original sample and is expressed as follows: CF = C i /C 0 .
Results and Discussion
Optimization Parameters of GP-MSE Concentration. During the preconcentration process, the sample solution in the microsyringe was continuously injected into the sample cell by a syringe pump at a certain rate while the evaporated solvent was carried away by a gentle nitrogen flow when the temperature at which the concentration was performed (hereafter referred to as the “concentration temperature”) was relatively low. The analytes are retained in the sample cell during this process due to their relatively low volatility. Thus, the main parameters that affect the preconcentration efficiency are the sample-injection rate, the concentration temperature and the nitrogen-flow rate. To determine the effects of these parameters, 0.5 mL stock solutions of standard PAHs (4 μg L –1 ) were tested, and the results were calculated for six different PAHs.
The Evaporation Rate of the Solvent in the GP-MSE Concentration Process. It is well known that the target analytes will lose, if there is no solvent,. Maintaining moist conditions for the sample cell during the concentration process is very important for the preconcentration efficiency because the analytes must be absorbed in a solvent medium. First, we determined the evaporation rates of dichloromethane in the GP-MSE concentration process at various temperatures: 50 μL min –1 , 117 μL min –1 and 150 μL min –1 at 50 °C, 70 °C and 80 °C, respectively. Clearly, at high temperature, faster evaporation occurs (see Figure S2 in the supplementary information), but higher temperature not investigated in order to make the control is convenient in the experiment process, and ensure good repeatability.
PPT Slide
Lager Image
Concentration efficiency (CE) plotted as a function of concentration temperature. The experimental conditions were as follows. Initial volume: 0.5 mL; initial concentration: 4 μg L–1; final volume: 10 μL; and N2-flow rate: 2.0 mL min–1. At 50 °C, 70 °C and 80 °C, the injection rates were 67 μL min–1, 133 μL min–1 and 167 μL min–1, respectively.
Optimisation of the Concentration Temperature. We evaluated the effect of temperature on the preconcentration efficiency at 50 °C, 70 °C and 80 °C, respectively. The optimisation results are shown in Figure 2. The preconcentration efficiency (CE) is plotted as a function of the concentration temperature. All of the selected injection rates were slightly higher than the evaporation rate of the solvent during the GP-MSE concentration process. At 50 °C, 70 °C and 80 °C, we chose injection rates of 67 μL min –1 , 133 μL min –1 and 167 μL min –1 , respectively. The sample cell was kept moist during the entire concentration process. It was determined that the temperature variation had no significant effect on the concentration efficiency. The concentration time was reduced as the temperature was increased. At 80 °C, the concentration time was the shortest: only 3 min plus the extraction time was required. Thus, 80 °C was set as the operating temperature in the standardised procedure in order to achieve of rapid extraction.
Optimisation of the Sample-injection Rate. Next, we evaluated the preconcentration efficiency at various sampleinjection rates. The results are shown in Figure 3. We selected three injection rates, namely 133 μL min –1 , 167 μL min –1 and 200 μL min –1 , at 80 °C. At an excessively low injection rate the loss of the low-boiling-point analytes was high, but no effect was observed for the high-boiling-point compounds. The probable reason for this effect is that the evaporation of the solvent is rapid, making the cell dry and causing the lowboiling- point analytes to be lost at a slow injection rate. However, when the sample-injection rate was higher than the solvent-evaporation rate, then the sample cell was kept wet, and the low-boiling-point analytes were kept in the chamber due to the moist conditions. However, when the injection rate was too fast, the sample cell was overloaded, leading to a downward flow of sample into the gas pipe. As a result, all of the analytes showed significant loss. Therefore, the injection rate of the sample must be higher than the evaporation rate of the solvent; in other words, the sample cell must be kept moist throughout the entire process to prevent the loss of analytes.
PPT Slide
Lager Image
Concentration efficiency (CE) plotted as a function of the injection rate. The experimental conditions were as follows. Initial volume: 0.5 mL; initial concentration: 4 μg L–1; final volume: 10 μL; concentration temperature: 80 °C; and N2-flow rate: 2.0 mL min–1.
Optimisation of the N2-flow Rate. In the GP-MSE preconcentration technique, the solvent is evaporated from the extract (sample matrix), and then it is carried by the gentle inert-gas stream to the waste vial. To understand the effect of the gas-flow rate on the concentration efficiency, the gasflow rates were fixed at 0.5, 2.0 and 4.0 mL min –1 (0.5 and 4.0 mL min –1 were the minimum and the maximum rates of the gas mass-flow controller unit, respectively), and then the results were compared as shown in Figure 4. The concentration efficiency remained nearly unchanged at the varying gas-flow rates because the gentle gas flow did not change the moist state of the sample cell. A gas-flow rate of 2.0 mL min –1 was fixed as the standard flow rate in the remainder of the experiment.
PPT Slide
Lager Image
Concentration efficiency (CE) plotted as a function of the flow rate of N2. The experimental conditions were as follows. Initial volume: 0.5 mL; initial concentration: 4 μg L–1; final volume: 10 μL; concentration temperature: 80 °C; and injection rate: 167 μL min–1.
Evaluation of the Method Performance. A standard mixture of PAHs was used to evaluate the method detection limit (MDL), reproducibility and linearity of the proposed thechnique. The MDL ranged from 0.5 to 15 ng L –1 (three times the signal-to-noise ratio). The reproducibility was indicated by the relative standard deviation (RSD) for a standard mixture of PAHs at the 4 μg L –1 level. The proposed method exhibited good reproducibility, with RSD values ranging from 1.5 to 9.0% for three repetitions. To investigate the linearity of the method by experimentation, five different concentrations (0.2, 0.4, 1, 5 and 20 μg L –1 ) of the standard mixture of PAHs were used. Good linearity was obtained at all five concentrations, with the correlation coefficient (R 2 ) ranging from 0.9905 to 0.9995. The concentration efficiency was determined under the optimised conditions of 80 °C, 167 μL min –1 (injection rate) and 2.0 mL min –1 (gas-flow rate). The concentration efficiency for the 4 μg L –1 standard mixtures ranged from 75.0 to 96.1%. The results are summarised in Table 1 .
The analytical characteristics of GP-MSE concentration method
PPT Slide
Lager Image
Abbreviations: MDL = method detection limit (three times of the signalto- noise ratio); Nap = naphthalene; AcPy = acenaphthylene; AcP = acenaphthene; Flu = fluorene; Phe = phenanthrene; AnT = anthracene; FluA = fluoranthene; Pyr = pyrene; B[a]A= beazo[a]anthracee; Chr = chrysene; B[b]F = beazo[b]fluoranthene; B[k]F = beazo[k]fluoranthene; B[a]P = beazo[a]pyrene; IND = indeno[1,2,3-cd]pyrene; DBA= dibenz- [a,h]anthracene; B[ghi]P = benzo[ghi]perylene.
Application. The presence of PAHs in the gas phase of ambient air was determined by the proposed method. First, the results obtained from the GP-MSE preconcentration were compared with those obtained from the conventional concentration method (using a nitrogen blow-down device). Eight PAHs were detected in the gas phase of the atmospheric sample, and the concentrations were 445.1, 12.3, 63.5, 227.3, 383.8, 29.0, 38.7 and 25.1 ng m –3 for naphthalene, acenaphthylene, acenaphthene, fluorene, phenanthrene, anthracene, fluoranthene and pyrene, respectively, as determined by GP-MSE preconcentration (the recoveries of the surrogate standards ranged from 86.8% to 92.3%). The fact that the results of the GP-MSE preconcentration were similar to the results obtained using a conventional nitrogen blowdown device confirmed the reliability of the new method. The concentration factor for GP-MSE was much higher than that for the conventional techniques by up to six orders of magnitude ( Table 2 ).
Comparison of the detection of PAHs in the gas phase of ambient air by the proposed method compared with that by two other concentration methods
PPT Slide
Lager Image
aNitrogen blow-down device. bConcentration factor obtained from GP-MSE. cConcentration factor obtained from conventional concentration technique.
Additionally, we used the proposed method to analyse the diurnal hourly concentration flux of PAHs in the gas phase of ambient air. The PAH concentrations in the gas phase are listed in Table 3 . The eight PAHs naphthalene, acenaphthylene, acenaphthene, fluorene, phenanthrene, anthracene, fluoranthene and pyrene were detected, and their concentrations ranged from 4.3 to 177.7 ng m –3 (the recoveries of the surrogate standards ranged from 88.0% to 94.5%). The concentrations of naphthalene, fluorene and phenanthrene were higher than those of the other molecules. The PAH concentrations in the gas phase of the atmospheric sample were at the ultra-trace level; thus, traditional methods would require a large sample volume for instrumental detectability. The GP-MSE preconcentration method showed high sensitivity, by which it was possible to detect the diurnal hourly concentration of PAHs in the gaseous phase of the ambient air.
Diurnal hourly concentrations of PAHs in the gas phase of ambient air
PPT Slide
Lager Image
Diurnal hourly concentrations of PAHs in the gas phase of ambient air
GP-MSE preconcentration is a sensitive sample-pretreatment method that exhibits excellent analytical performance for the determination of trace semivolatile organic pollutants. A high concentration efficiency was achieved when the injection rate and the temperature were adjusted appropriately and the sample cell was kept moist during the concentration process. The GP-MSE preconcentration method achieved a microlitre-level reduction and greatly reduced the initial sampling volume. The diurnal hourly flux of PAHs in the gas phase of ambient air was determined by this method, which holds good prospects for the detection of trace- to ultra-trace-level semivolatile organic pollutants in the environment using relatively small sample volumes.
This study was supported by a grant from the National Natural Science Foundation of China (No. 21165020, 21265023).
Li Z. , Li D. , Oh J. R. , Je J. G. 2004 Chemosphere 56 611 -    DOI : 10.1016/j.chemosphere.2004.04.022
Wu N. , Huang H. , Zhang S. , Zhu Y. G. , Christie P. , Zhang Y. 2009 Environ. Pollut. 157 1613 -    DOI : 10.1016/j.envpol.2008.12.022
Kim S. K. 2012 Chemosphere 89 995 -    DOI : 10.1016/j.chemosphere.2012.07.016
Lee S. , Kim B. , Kim J. 2013 J. Chromatogr. A 1277 35 -    DOI : 10.1016/j.chroma.2012.12.047
Ferrer C. , Lozano A. , Agüera A. , Jiménez Girón A. , Fernández-Alba A. R. 2012 J. Chromatogr. A 1218 7634 -
Kannan N. , Petrick G. , Bruhn R. , Schulz-Bull D. E. 1998 Chemosphere 37 2385 -    DOI : 10.1016/S0045-6535(98)00296-3
Primbs T. , Piekarz A. , Wilson G. , Schmedding D. , Higginbotham C. , Field J. , Simonich S. M. 2008 Environ. Sci. Technol. 42 6385 -    DOI : 10.1021/es702160d
Hayward S. J. , Gouin T. , Wania F. 2010 Environ. Sci. Technol. 44 3410 -    DOI : 10.1021/es902512h
Arinaitwe K. , Kiremire B. T. , Muir D. C. G. , Fellin P. , Li H. , Teixeira C. , Mubiru D. N. 2012 Environ. Sci. Technol. 46 11524 -    DOI : 10.1021/es302238w
Zhang Z. L. , Zhou J. L. 2007 J. Chromatogr. A 1154 205 -    DOI : 10.1016/j.chroma.2007.03.105
Kannan N. , Petrick G. , Schultz-Bull D. E. , Duinker J. C. 1993 J. Chromatogr. A 642 425 -    DOI : 10.1016/0021-9673(93)80108-K
Herrero P. , Borrull F. , Marcé R. M. , Pocurull E. 2013 J. Chromatogr. A 1285 31 -    DOI : 10.1016/j.chroma.2013.02.034
Oellig C. , Schwack W. 2012 J. Chromatogr. A 1260 42 -    DOI : 10.1016/j.chroma.2012.08.088
Choi S. D. , Baek S. Y. , Chang Y. S. 2008 Atmos. Environ. 42 2479 -    DOI : 10.1016/j.atmosenv.2007.12.032
Pereira M. G. , Walker L. A. , Wright J. , Best J. , Shore R. F. 2009 Environ. Sci. Technol. 43 9010 -    DOI : 10.1021/es901805e
Liu S. , Tao S. , Liu W. , Liu Y. , Dou H. , Zhao J. , Wang L. , Wang J. , Tian Z. , Gao Y. 2007 Environ. Sci. Technol. 41 8256 -    DOI : 10.1021/es0716249
Lohmann R. , Dapsis M. , Morgan E. J. , Dekany V. , Luey P. J. 2011 Environ. Sci. Technol. 45 2655 -    DOI : 10.1021/es1025883
Sae-Khow O. , Mitra S. 2010 Anal. Chem. 82 5561 -    DOI : 10.1021/ac100426y
Saridara C. , Ragunath S. , Pu Y. , Mitra S. 2012 Anal. Chim. Acta 677 50 -
Sae-Khow O. , Mitra S. 2010 J. Chromatogr. A 1217 2736 -    DOI : 10.1016/j.chroma.2009.12.043
Deng X. , Xie P. , Qi M. , Liang G. , Chen J. , Ma Z. , Jiang Y. 2012 J. Chromatogr. A 1219 75 -    DOI : 10.1016/j.chroma.2011.11.031
Preter V. D. , Staeyen G. V. , Esser D. , Rutgeerts P. , Verbeke K. 2009 J. Chromatogr. A 1216 1476 -    DOI : 10.1016/j.chroma.2008.12.095
Soria A. C. , Martínez-Castro I. , Sanz J. 2009 J. Chromatogr. A 1216 3300 -    DOI : 10.1016/j.chroma.2009.01.065
Lara-Gonzalo A. , Sánchez-Uría J. E. , Segovia-García E. , Sanz- Medel A. 2008 Talanta 74 1455 -    DOI : 10.1016/j.talanta.2007.09.036
Yang C. , Piao X. , Qiu J. , Wang X. , Ren C. , Li D. 2011 J. Chromatogr. A 1218 1549 -    DOI : 10.1016/j.chroma.2011.01.018
Li D. , Yang C. , Ren C. 2008
Yang C. , Qiu J. , Ren C. , Piao X. , Li X. , Wu X. , Li D. 2009 J. Chromatogr. A 1216 7694 -    DOI : 10.1016/j.chroma.2009.09.039