Quantifying uncertainty in the measurement of arsenic in suspended particulate matter by Atomic Absorption Spectrometry with hydride generator
Chemistry Central Journal volume 5, Article number: 17 (2011)
Arsenic is the toxic element, which creates several problems in human being specially when inhaled through air. So the accurate and precise measurement of arsenic in suspended particulate matter (SPM) is of prime importance as it gives information about the level of toxicity in the environment, and preventive measures could be taken in the effective areas. Quality assurance is equally important in the measurement of arsenic in SPM samples before making any decision. The quality and reliability of the data of such volatile elements depends upon the measurement of uncertainty of each step involved from sampling to analysis. The analytical results quantifying uncertainty gives a measure of the confidence level of the concerned laboratory. So the main objective of this study was to determine arsenic content in SPM samples with uncertainty budget and to find out various potential sources of uncertainty, which affects the results. Keeping these facts, we have selected seven diverse sites of Delhi (National Capital of India) for quantification of arsenic content in SPM samples with uncertainty budget following sampling by HVS to analysis by Atomic Absorption Spectrometer-Hydride Generator (AAS-HG). In the measurement of arsenic in SPM samples so many steps are involved from sampling to final result and we have considered various potential sources of uncertainties. The calculation of uncertainty is based on ISO/IEC17025: 2005 document and EURACHEM guideline. It has been found that the final results mostly depend on the uncertainty in measurement mainly due to repeatability, final volume prepared for analysis, weighing balance and sampling by HVS. After the analysis of data of seven diverse sites of Delhi, it has been concluded that during the period from 31st Jan. 2008 to 7th Feb. 2008 the arsenic concentration varies from 1.44 ± 0.25 to 5.58 ± 0.55 ng/m3 with 95% confidence level (k = 2).
In chemical metrology most of the important decisions are based on the quality and the reliability of results of quantitative analysis. In this context it is also important to have reliability of the results in the international trade and analytical results should be acceptable to all users within the country or outside the country. This can be achieved by the estimation of uncertainty in the measurements as per ISO/EURACHEM guidelines [1–5]. In this regard ISO has published "Guide to the expression of uncertainty in measurement (GUM)" in 1995 in collaboration with BIPM, IEC, IFCC, IUPAC, and OIML. Since last decades there is an increasing interest in arsenic study in the environment because of its toxicity to plants, animals and human beings. The major contamination sources of such toxic metals in the air are natural (volcanoes) as well as anthropogenic like industries [6, 7]. The industrial sectors like metal smelting, refining sector contributes major sources of contamination in the environment. The arsenic in atmosphere is generally found associated with fine particles (<2 μm), as arsenate (+5 oxidation state) and arsenite (+3 oxidation state). Recent study shows that lung cancer is significantly higher in the areas where arsenic concentration is 1.77 ng/m3 or more . Such studies become important when experimental results conclude that lung cancer is significantly higher in the areas where arsenic concentration is 1.77 ng/m3 or more. In the measurement of arsenic by AAS-HG there are several interference due to transition metal ions, which forms hydride during reaction with sodium borohydride, sodium hydroxide and hydrochloric acid [9–12]. So it is important to measure all the possible sources of uncertainty in measurement before make any conclusion. As the determination of arsenic with the uncertainty in an environmental sample is a tedious job as there are so many parameters from sampling to final value, which influences the result.
The evaluation of the uncertainty at every step of the experiment is one of the requirements of the standard ISO/IEC17025 for certain test methods to get accreditation. In the last several years uncertainty measurement have been carried in different materials by using different techniques [13–16], but there are hardly studies related to such complex matrix in environmental samples.
For the measurement of arsenic and uncertainty measurement in SPM samples by AAS-HG, we have collected seven samples of SPM on cellulose filter paper through HVS in diverse sites of Delhi. The locations of sampling site are marked in map of Delhi given in Figure 1. The Figure 2 represents the SPM concentration during the sampling period. There are several sources of uncertainty in chemical metrology like incomplete definition, sampling, environmental conditions, uncertainties of masses, volumes, equipment, reference values, measuring equipment approximation, assumptions incorporated in experimental methods, random variations, etc. The uncertainty is of two types 'Type A' (evaluated using statistical analysis of a series of observations) and 'Type B' (evaluated using other means than the statistical analysis of a series of observations). So in proposed study, an attempt has been made to choose such a complex matrix for the evaluation and expression of uncertainty in measurement of arsenic in SPM samples by AAS-HG, following documents published by ISO/IEC 17025:2005, EURACHEM guideline and GUM document.
2. Results and discussion
The arsenic concentration varies from country-to-country and depends mostly upon the industrial and transport emission. In Delhi the arsenic concentration were found to be from 1.44 ± 0.25 to 5.58 ± 0.55 ng/m3 at seven different sites in proposed period. The maximum concentration of arsenic was found at Ashram zone (S-1 site; 5.58 ± 0.55 ng/m3) and minimum concentration was found at Pitampura and NPL zone (S-4, S-7; 1.44 ± 0.25 ng/m3). The high concentrations of arsenic at Ashram is mainly due to transport, industry and thermal power plant near by ashram and in Loni Road area which is densely populated zone, the house hold activities, industrial and transport contribute towards more arsenic contamination in the air. In Spain a study has been carried out in the year 2001 and 2002 in PM2.5 and an average concentration of 6.4 ng/m3 has been reported which, is more than the concentration of arsenic observed from various sites of Delhi. In another study of Yoshikawa concluded that lung cancer is significantly higher in the areas where arsenic concentration is 1.77 ng/m3 or more. So it is evident from the proposed study that five site out of seven sites of Delhi have more arsenic then the reported value. So these are the highly risky zones due to high level of arsenic content in the air.
The precise and accurate determination of arsenic in SPM sample is a tedious job due to huge losses during sample digestion process and due to vaporizing nature. At the same time evaluation of combined uncertainty in such a complex matrix is very complicated as there are various parameters, which contribute uncertainty in the entire process. On determining the uncertainty for arsenic content measurement we were redundantly meticulous in identifying all uncertainty sources. In a routine determination of uncertainty we recommended that it is possible to be much less pedantic, and only reveal and realistically quantify the most significant components. The uncertainty components, which contribute a negligible value, could be ignored from the calculation of combined uncertainty value. The uncertainty components, which are found to contribute significant value during the subsequent computations they can be quantified more precisely and final combined uncertainty, can be recalculated. In the determination of arsenic concentration with uncertainty budget in environmental sample there are five major sources, which directly influence the final results those are [C As ], [V Test ], [W TSPM ], [V HVS ] and [W ANA ]. The repeatability under [C As ] is the major component, which contributes maximum uncertainty in arsenic determination. Repeatability depends upon various factors like volume of sample and acid taken for measurement in hydride generator. So we should take proper precautions during analysis. The reference standard stock solution is also other crucial parameters. So for the calibration of the instrument the standard stock solution must be off having low uncertainty. The Flow rate is also parameter which contributes uncertainty in HVS sampling, as flow rate get changed by several factors such as oily nature of the sample, presence of photochemical smog, wood smoke may block the filter paper and cause a rapid drop in air flow. Besides these high humidity and dense fog can cause the filter to become wet and thus airflow reduces through the filter paper. Sampling through HVS is another contributor for uncertainty in arsenic determination.
A generalized scheme for planning a measurement and a simple, practical approach to estimate and combining uncertainties has been demonstrated for the determination of arsenic in SPM samples. The results show that the arsenic concentration varies 1.44 ± 0.25 ng/m3 to 5.58 ± 0.55 ng/m3 with variation of sites with 95% confidence level. The experimental result shows that the heavy traffic site and industrial site are the major sources for the contamination of arsenic in environment.
4. Experimental Section
4.1 Instruments and Apparatus
Atomic Absorption Spectrometer (AAS) of Analytik Jena make Vario-6 with a hydride generator accessory was used for the analysis of arsenic. The arsenic measurements were carried out at 193.7 nm wavelength at optimum conditions of acids, temperature and reducing agent. The pipette of 5 mL and volumetric flask of 50 mL capacity used were of Borosil glass works India Limited. The pipettes and volumetric flasks were calibrated prior to analysis by National Physical Laboratory-India (National Metrology Institute of India) following international standard procedure and protocol. All the acid digestion and dilution work was carried out in cleaned laminar flow bench equipped with the proper exhaustive system. The sampling of SPM was carried out using High Volume Sampler (HVS-410, Make: Envirotech) and was collected on pre-desiccated cellulose filter papers. No chemical treatment of cellulose filter paper were carried out before sampling, but the filter paper were dried for 3 hrs in an oven at 105°C before sampling to remove moisture if any and to get constant weight. Flow rate of HVS was kept about 1.0 ± 0.2 m3/min. and the sampling time was eight hours.
Nitric acid (69%), Hydrochloric acid (35%) of GR grade (Guaranteed Reagent), which were further purified by sub boiling point distillation in a quartz glass device, Hydrogen Peroxide (50%), Sulphuric acid (98%) all E. Merck (India) make were used. De-ionized water (18 mega ohm resistivity) prepared from Millipore milli-Q element water purification system, USA was used throughout the process. The reference standards stock solution of 1.0 mg/litre (BND 301.03) of arsenic was used prepared by NPL-India. The stock solution was diluted upto working range of AAS-HG for arsenic and these solutions were used for the calibration of the AAS-HG instrument.
4.3 Sampling Site and Description
New Delhi is located 160 km south of the Himalayas at latitude 28°24' to 28°53'N and longitude 76°20'E to 77°20'E with an altitude of about 216 meter above mean sea level. Delhi's climate is mainly influenced by its inland position and the prevalence of continental air during major part of the year. Delhi has three distinct seasons: summer, monsoon and winter. For the proposed study seven different sites were selected namely Ashram (S-1), which is a heavy traffic zone and thousands of vehicles crossing this junction per hour, Azadpur (S-2) is a thickly populated residential area in North West Delhi district. It is India's largest wholesale market for vegetable and fruits called the Azad Market. Due to several industries in this area it is consider as industrial zone also. Loni Road (S-3) is densely populated area and also a traffic zone, Pitampura (S-4) is neither a heavy traffic zone nor a industrial zone, so we have considered this as an ambient zone, Highway no. 56 (S-5) is a traffic zone, Naraina (S-6) is an industrial zone, it also includes Mayapuri industrial area and NPL (S-7) is considered to be mixed zone. Two sides of NPL are surrounded by ridge area whereas other two sides are surrounded residential area, with traffic movement. The SPM samples were collected between 9 am to 6 pm because of the major operational activities carried out during daytime and in the night time the filter paper get moist due to dew and flow rate of HVS also decrease sharply after heavy loading of SPM. The sampling was carried out between 31st Jan. 2008 to 7th Feb. 2008.
4.4 Sample digestion procedure
For the determination of arsenic two replicates of 18 × 11.5 cm size (known weight at deposited area) of all the seven sites were taken into a separately cleaned air tight poly tetra fluoro ethylene (PTFE) vessel. To these vessels 15 mL of sub boiled nitric acid, 5 mL of hydrogen peroxide and 5 mL of de-ionized water were added and it was closed tightly with PTFE vessel lid. The PTFE vessels-containing samples were kept for overnight at room temperature to oxidized arsenic and released from SPM samples. These steps minimize the risk of losing of arsenic during digestion. After keeping the PTFE vessel for overnight, 4 to 5 drops of sulphuric acid was added to each vessel and PTFE vessel was kept on hot plate by covering with Teflon lid at 100°C for 15-20 minutes and then heated at 150°C -160°C. Further the lid of the PTFE vessel was opened and it was heated again after addition of 10 mL of hydrochloric acid. Further in syrupy condition the samples were heated again with 5 ml hydrochloric acid to ensure the complete removal of nitric acid. Finally the solutions were boiled once again with 15 ml 1:1 de-ionized water and hydrochloric mixture. After digestion the whole content in the vessel was transferred into 100 mL capacity quartz centrifuge tube and centrifuged at 5000 rotations per minute for 2-3 minutes to separate solid particles from the solutions. The supernatant liquid was transferred in to a glass beaker. The remaining residue in the vessel was washed several times with hot water to ensure complete transfer of digested sample from the vessel. Further the solutions were heated at 75°C for 45 minutes with 10 mL of 25% potassium iodide and 15 mL of 25% ascorbic acid solution to convert As (V) to As (III). The final solutions were made 50 mL by de-ionized water. In the final solution arsenic was determined by AAS-HG using respective standards and principal resonance line (193.7 nm) after reducing arsenic in the presence of sodium borohydride (3%), sodium hydroxide (1%) and hydrochloric acid (5%) using electro thermal heating at 900°C. Mean absorbance values of six readings of two replicates (three each) were taken into the consideration for calculation of the concentration of arsenic. Procedural blank for cellulose filter paper without sample was also done to check the blank levels for arsenic and correction applied wherever required. The details of the procedure are published elsewhere .
4.5 Determination of concentration of arsenic in SPM sample
In the determination of arsenic in SPM samples major sources of the uncertainty have been included in the calculation of combined uncertainty according to the EURACHEM/GUM guidelines. The concentration of arsenic C(M As ) in the SPM samples has been evaluated using AAS-HG technique by following equation.
Where; C(M As ) = Concentration of arsenic in SPM (ng/m3); C AS = Concentration of arsenic analyzed from AAS-HG in μg/litre; V TEST = Volume made of the test sample after wet digestion process in mL; W TSPM = Total weight of SPM deposit on 18 cm × 23 cm area; V HVS = volume of the air processed through HVS; W ANA = Weight of SPM deposit on 18 cm × 11.5 cm area taken for wet digestion process. The evaluated values for the above factors are given in table 1.
In this manuscript we have taken some important contributions that are stated above in equation-1. The parameters which contributed towards the uncertainty are also shown in 'Fish- bone or Cause and effect or Ishikawa diagram in Figure 3.
In accordance with GUM , the combined uncertainty for the mathematical model, which is a product or quotient form, is given by:
The sensitivity coefficient = piy/xi. Where pi is the power of the terms in the equation (1).
p1 for C AS = +1; p2 for V TEST = +1; p3 for W TSPM = +1; p4 for W ANA = -1 and p5 for V HVS = -1;
The combined uncertainty is given by:
The uncertainty evaluation in every step of the experiment has been discussed in the following sub-sections:
5.0 Uncertainty components in the measurement of concentration of arsenic obtained from AAS-HG [C AS ]
Referring to Figure 4, the major sources of the uncertainty for C AS are due to reference standard stock solution, repeatability, dilution in 50 mL volumetric flask from stock solution to working range and 5 mL capacity pipette used for dilution.
5.1 Reference standard stock solution
As per the certificate the strength of arsenic standard solution is 1.0 mg/L with 0.02 mL uncertainty. This stock solution was used for dilution up to working range. The details are given in table 2.
5.2 Uncertainty evaluation due to repeatability
The repeatability has also been considered for the determination of arsenic. The concentration of the arsenic was determined experimentally by analyzing six determinations of two replicates (three each). The standard deviation in values obtained by AAS-HG and repeatability has been calculated by using following EURACHEM guide. The details are given in table 3.
5.3 Uncertainty in 50 ml volumetric flask used for dilution from 1.0 mg/L (1000 μg/L) to working range
The final volume of standard arsenic solution made was 50 mL after several dilutions. The major uncertainties were due to temperature variation and calibration accuracy as per manufacturer's specification. The details are given in Figure 5 and table 4.
The uncertainty associated with the use of the volumetric flask of 50 mL capacity at a temperature other than the calibration temperature (temperature difference) is given by:
Where V is the volume of volumetric flask and ΔT is temperature difference
The same effect of temperature has been considered three times in the measurement of arsenic and the values have been calculated using equation (4): 50 mL volumetric flask used for dilution from stock solution to working range (refer to table 4), 5 mL pipette used for dilution (refer to table 5) and again 50 mL volumetric flask used for making volume after wet digestion of sample (refer to table 6).
5.4 Uncertainty due to 5 mL pipette used for dilution
5 mL aliquot was taken by 5 mL pipette for making dilution upto 50 mL. The main uncertainties were due to temperature variation and calibration accuracy as quoted by the manufacturer. The effect of temperature on 5 mL pipette has also been considered as per section 5.2. The details are given in Figure 6 and table 5.
Combined uncertainty due to section 5.0 [C AS ] has been calculated by the sum of subsection 5.1, 5.2, 5.3 and 5.4 as follows.
Similarly the combined uncertainty for sample S3-S7 is calculated.
6.0 Uncertainty evaluation in 50 mL volumetric flask used in making volume after wet chemical digestion [V TEST ]
The known weight of the sample was processed through wet chemical route and final volume was made 50 mL. Referring to Figure 7, the sources of uncertainty for V TEST are due to calibration uncertainty and effect of temperature variation on 50 mL volumetric flask. The uncertainty associated with temperature for 50 mL volumetric flask is carried out as per section 5.3. The details are given in Figure 7 and table 6.
Combined uncertainty of
7.0 Uncertainty evaluation due to volume of air processed through high Volume sampler [V HVS ]
Referring to Figure 8 the sources of the uncertainty for V HVS are due to the uncertainty due to flow rate of air through HVS, time recording by stop watch, ambient pressure measurement and measurement of temperature of the site.
7.1 Uncertainty due to flow rate of air through HVS
The volume of air sampled through HVS has been calculated on the basis of equation-5 and volume of air processed at standard temperature and pressure (STP) is calculated by using equation-6 given below.
Where, V1- Total volume of air processed through HVS, P1- Atmospheric pressure on the date of measurement (mm of mercury; table 7); T1- Standard temperature (273); T2- working temperature [temperature of the sampling site °C + standard temperature] and P2- Standard pressure (760 mm of mercury). Substituting value from table (7) in equation-6
Volume of air processed through HVS at
The detail of uncertainty due to flow rate of air is given in table 8.
7.2 Uncertainty due to ambient pressure measurement
The onsite pressure of the atmosphere was measured in millibar (mbar) by using ambient pressure measurement instrument and the value measure was converted in mm of mercury unit by multiplying mbar value by 0.75. The details are given in table 8.
7.3 Uncertainty in time recording due to Time Tracker Recorder
The total time was recording by using inbuilt time tracker recorder in the HVS used for processing the atmospheric air. The details are given in table 8.
7.4 Uncertainty in measuring temperature of the site
The onsite temperature of the open atmosphere was measured by calibrated thermometer three times in a day and average of three values were used for calculation of total volume of air processed at STP. The details are given in table 8.
Putting values in above from table 8 the combined uncertainty due to HVS;
8.0 Uncertainty evaluation due to balance used for weighing the blank and SPM deposited filter paper [W TSPM ] and [W ANA ]
The total weight of the SPM deposited on filter paper and weight of the filter paper used for wet digestion process was calculated as below. Taking example of sample-1.
Total area of the cellulose filter paper used for sampling = 20.3 cm × 25.4 cm = (515.6 cm2)
Total area of the filter paper used for deposition (18 cm × 23 cm) = 414 cm2
Total weight of the cellulose filter paper used for sampling (515.6 cm2) = 4.255 g
Empty weight of the filter paper of same deposit area (414 cm2) = 3.4164*g
Formula for calculation of weight on 414 cm2 area =
(It has been considered that the thickness of the cellulose filter used is uniform through out the area. On this basis we have calculated the weight of the deposit area (414 cm2).
Total suspended particulate matter deposit on the filter paper (414 cm2) = 3.6374
Weight of the SPM deposited (V-IV) = 0.2210 g
Weight of the SPM deposited area (18 cm × 11.5 cm) taken for analysis = 0.1105 g
The weight calculation for other samples (S-2 to S-7) were done in same way as given above for sample-1.
9. Calculation of concentration and associated uncertainty of arsenic in S-1 sample
The concentration of arsenic in the seven SPM has been determined by AAS-HG using equation-1. The uncertainty associated with the concentration of arsenic has been calculated using the uncertainty components. Substituting the values from table 1 into equation-1, the concentration of arsenic (ng/m3) in S-1 sample:
The combined uncertainties in the determination of arsenic in case of sample S-1 was carried out using equation-3. Putting values from summary table 10, combined uncertainty;
Therefore Concentration of arsenic for S-1 = 5. 58 ± 0.55 ng/m3
In the same way the concentration of arsenic in sample S-2 has been calculated. Putting values from table 1 concentration of arsenic is;
And the combined uncertainty in the determination of arsenic in case of sample S-2 is carried out using equation-3. Putting values from summary table 10, the combined uncertainty;
Therefore Concentration of arsenic for S-2 = 3.21 ± 0.39 ng/m3
Similarly combined uncertainties for concentration of arsenic in samples S-3 to S-7 have been calculated. Table 11 shows the analytical results of the arsenic with uncertainty at seven sites.
International Organization for Standardization: Guide to the expression of uncertainty in measurement. Geneva, Switzerland: ISO 98. 1995, 2
International Organization for Standardization: International vocabulary of basic and general terms in Metrology (VIM). Geneva, Switzerland: ISO 99. 2007, 3
International Organization for Standardization: General Requirements for the Competence of testing and calibration laboratories. Geneva, Switzerland: ISO 17025. 2005
International Organization for Standardization: Guide to the Use of Repeatability, Reproducibility and Trueness Estimate in Measurement Uncertainty Estimation. Geneva, Switzerland: ISO 21748: 2003
Ellison SLR, Williams A: EURACHEM/CITAC Guide: Use of uncertainty information in compliance assessment. 2007, Cooperation on the International Traceability in Analytical Chemistry: Switzerland, First
Baudo R, Giesy JP, Muntau H: Sediments: Chemistry and Toxicity of in-place pollutants. 1990, Lewis Publishers Inc.: Boca Raton, Florida
Moore W, Ramamoorthy S: Heavy metals in natural waters: Applied Monitoring and Impact Assessment. 1984, Springer, New York
Yoshikawa M, Aoki K, Ebine N, Kusunoki M, Okamoto A: Correlation between the arsenic concentrations in the air and the SMR of lung cancer. Environ Health Prev Med. 2008, 13: 207-218. 10.1007/s12199-008-0032-2.
Barth P, Krivan V, Hausbeck R: Cross-interferences of hydride-forming elements in hydride-generation atomic absorption spectrometry. Anal Chim Acta. 1992, 263: 111-118. 10.1016/0003-2670(92)85432-6.
Walcerz M, Bulska E, Hulanicki A: Study of some interfering process in the arsenic, antimony and selenium determination by hydride generation atomic absorption spectrometry. Fressenius J Anal Chem. 1993, 346: 622-626. 10.1007/BF00321258.
Erdem N, Henden E: Inter-element interferences in the determination of arsenic and antimony by hydride generation atomic absorption spectrometry with a quartz tube atomizer. Anal Chim Acta. 2004, 505: 59-65. 10.1016/S0003-2670(03)00199-5.
Singh Nahar, Ramchandran R, Sarkar AK: Determination of arsenic in Fly ash and sulphide ore by flame Atomic Absorption Spectrometry using Hydride generator. Asian J Chem. 2003, 15: 1593-1597.
Synek V, Subrt P, Marecek J: Uncertainties of mercury determinations in biological materials using an atomic absorption spectrometer: AMA 254. Accred Qual Assur. 2000, 5: 58-66. 10.1007/s007690050012.
Kayal N, Singh Nahar, Ojha VN, Gupta PK: Evaluation and expression of uncertainty in the determination of alumina in deodorants by complexometric method. Test and Eva (ASTM). 2009, 37: 321-328.
Kumar AR, Riyazuddin P: Underestimation of total arsenic concentration in groundwater samples determined by hydride generation quartz furnace atomic absorption spectrometry due to sample characteristics. Accred Qual Assur. 2007, 12: 455-458. 10.1007/s00769-007-0289-y.
Kayal N, Singh Nahar: The quantitative estimation of silica in rice husk ash by titrimetric method: A case study for uncertainty calculation. Mapan; A Journal of Metrology Society of India. 2010, 25 (2): 115-123.
Singh Nahar, Chauhan S, Singh K, Saud T, Saxena M, Soni D, Mandal TK, Bassin JK, Gupta PK: Determination of arsenic and mercury metals in suspended particulate matter by flame/flameless atomic absorption spectrometer. J Atmos Poll Res. 2010, 1: 112-117.
The authors are grateful to Prof. R. C. Budhani, Director, National Physical Laboratory for encouragement and permission to publish this work. The authors also duly thanks to Dr. J. K. Bassin, Director, NEERI for providing their sites for sampling.
The authors declare that they have no competing interests.
NS participated in the design of the study, carried out sampling, wet digestion of the sample and written work. VNO has contributed in calculation of various uncertainties, which have been identified in the experiments. NK has contributed in writing the manuscript. TA has contributed in analysis of the arsenic analysis by AASHG-PKG has given valuable suggestion during experiment as well written work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Singh, N., Ojha, V., Kayal, N. et al. Quantifying uncertainty in the measurement of arsenic in suspended particulate matter by Atomic Absorption Spectrometry with hydride generator. Chemistry Central Journal 5, 17 (2011). https://doi.org/10.1186/1752-153X-5-17