The main sources of uncertainty encountered during the analysis of the mass concentration of metals in ambient air as part of the operation of the UK Heavy Metals Monitoring Network are presented. It is observed that the uncertainty contribution from possible variation in the isotopic composition of the sample depends on the element in question, but can be significant (e.g., for Pb, Cd, and Hg). The working curve method for the ICP-MS analysis of metals in solution, with a low resolution, high throughput instrument measuring at one m/z ratio per element, relies on the relative abundance of the isotopes under consideration being the same in both the sample and the calibration solution. Calculation of the uncertainty in this analysis assumes that the isotopic composition variation within the sample and calibration solution is limited to a defined range. Therefore, in order to confirm the validity of this quantification methodology and its uncertainty budget, the isotopic composition of the calibration standards used for quantification has been determined. The results of this analysis are presented here.
1. Introduction
The general public
and the environment can be exposed to several classes of hazardous compounds
containing metallic elements which occur naturally or are released by domestic
or industrial processes [1].
The total concentration levels of Pb, Ni, As, and Cd allowable in the PM10 fraction of ambient air (particles with an aerodynamic diameter of 10 μm or
less) are now limited by European legislation [2–4].
In order to enforce this legislation, to measure human and environmental
exposure, and to show compliance with limit and target values, the total
concentration levels of ambient metals, at multiple sites on nationwide air
quality monitoring networks, need to be measured. To this end, nationwide
networks for the measurement of a wide range of particulate-borne and gaseous
pollutants are now well established in many developed countries around the
world. NPL currently manages and
operates the UK Heavy Metals Monitoring Network (the “network”) on behalf of
the UK Department for Environment, Food and Rural Affairs (Defra). The “network” consists of 24 monitoring sites
around the UK collecting PM10 particulate matter which is then sent
back to NPL for the analysis of the mass concentration of Ni, As, Cd, Pb, (as
required by European legislation) and also Hg, Cr, Cu, Fe, Mn, V, Zn, and Pt to
contribute to long-term UK data sets [5]. Whilst the data quality objectives laid down
by the air quality legislation are not especially exacting—the maximum
allowable expanded uncertainty for Pb determination is 25%, and for Ni, As, and
Cd is 40%—it is still necessary
to ensure that these objectives are routinely and consistently met.
Additionally NPL also sets a self-imposed maximum measurement uncertainty of
40% on the nonmandated metals.
To determine the mass concentration of particulate-phase metals in ambient
air, particles are collected onto air filters which are then digested in acid
before being analysed by inductively coupled plasma-mass spectrometry (ICP-MS).
The ICP-MS is calibrated by a working curve method using matrix matched
solutions prepared from commercially available elemental solutions certified
for metal mass fraction, and compared with NIST standard reference material
(SRMs) by NPL to ensure consistency and accuracy. Whilst NPL includes a
component of uncertainty to take account of the possibility of the isotopic
composition of metals in ambient air varying within natural limits, for a
working curve method this assumes that the isotopic composition of the
calibration standards being used also falls within this assumed range. If it
does not, then an additional uncertainty contribution, or a correction factor,
may need to be applied. In order to validate this hypothesis, the isotopic
composition of the standards used for these routine analyses requires measurement. This work
presents the results of this analysis and the effect of these findings on the
uncertainty budget of the measurement. We also present the contributions to
overall uncertainty budget from the possible variations in the isotopic
compositions of the ambient air samples.
2. Experimental
Particulate
samples were taken at all sites in the “network” using Partisol 2000
instruments (fitted with PM10 heads) operating at a calibrated flow
rate, nominally of 1m3⋅h−1, in accordance with European
standard method EN 12341 [6]. Samples were taken for a period of one week [7]
onto 47 mm diameter GN Metricel membrane filters. The analysis for
particulate-phase metals took place using a PerkinElmer Elan DRC II ICP-MS,
following NPL's UKAS accredited procedure, which is fully compliant with the
requirements of European standard method EN 14902 [8]
(the EU “reference method” for the analysis of metals in ambient air). Upon arrival at NPL, the filters sampled with
particulate matter were cut accurately in half, and each portion digested at
temperatures up to 220°C using a CEM Mars X microwave. The digestion mixtures
used were as follows.
Hg and Pt: 5 ml
of nitric acid and 5 ml hydrochloric acid.
All other metals:
8 ml of nitric acid and 2 ml hydrogen peroxide.
These digested solutions were then diluted with deionised water (Millipore, Milli Q, Mass, USA) prior to analysis. ICP-MS
analysis took place as previously described [9] using at least four-matrix-matched gravimetrically
prepared calibration solutions [10]
prepared from monoelemental standard solutions (VWR, checked for total
elemental composition against the NIST SRM 3100 series). A detector dead time
correction was applied [11]
and a full span dual detector linearity check was performed in order to
minimise any detector nonlinearity [9] since the concentration of different isotopes within
the samples may span several orders of magnitude. A quality control standard
was repeatedly analysed (after every two solutions), and the change in response
of the quality control standard was mathematically modelled to correct for the
long-term drift of the instrument. The short-term drift of the ICP-MS was corrected by the use of an internal
standards mixture (containing Y, In, Bi, Sc, Ga, and Rh) continuously added to
the all samples via a mixing block. Each sample was analysed in triplicate,
each analysis consisting of five replicates.
For each element, one isotope at one m/z value was chosen and
monitored. The mass of each metal in solution (and its uncertainty) was then
determined by a method of generalised least squares using XLGENLINE (an
NPL-developed programme [12])
to construct a calibration curve. The analysis of the isotopic ratios of the
calibration standards was performed by determining the blank-corrected
intensity at each appropriate m/z ratio.
The isotopic composition of the calibration standards was measured by
analysing the separate monoelemental standard solutions used to make up the
calibration solutions. All the usual corrections for isobaric and polyatomic
corrections were applied. The mass fraction of the sum of all isotopes of metal
analyte in each standard solution was approximately 1 μg/g. All isotopes of
the elements of interest were measured, not simply the ones used for
quantification. Whilst no additional effort was made to determine additional
corrections for isotopes not usually used for quantification by the NPL
procedure, additional type-B uncertainty components were included in the
uncertainty budget to account for unresolved inaccuracies owing to mass bias
and mass discrimination effects, based on conservative estimates from existing
literature data [13].
Given the low precision of the measurements, these factors are expected to have
a minimal contribution to the overall uncertainty [14]. Uncertainty contributions were
also added to account for residual dead time effects and detector nonlinearity.
Matrix effects were minimised by the matrix matching of all solutions prior to
analysis. Total expanded measurement uncertainties for each analysis are
calculated using a full GUM [15]
approach and are expressed with a coverage factor of k=2 representing
the 95% confidence interval.
3. Results and Discussion
The working curve
method for the ICP-MS analysis of metals in solution, with a low resolution,
high throughput instrument measuring at one m/z ratio per element relies
on the relative abundance of the isotopes under consideration being the same in
both the sample and the calibration solution. If this is not the case, a
multiplicative bias in the observed results for element X, δX,
will be observed, and is given by
δX=na,cal⋅∑ini,samna,sam⋅∑ini,cal,
where na,cal is the amount of isotope a of element X
in the calibration standard, ∑ini,cal is the amount of all isotopes of element X in
the calibration standard, na,sam is the amount of isotope a of element X
in the sample, and ∑ini,sam is the amount of all isotopes of element X in
the sample [16].
When the relative abundance of the isotope used for quantification in both the
sample and the calibration solution is the same, then δX=1 and no bias is observed. It is interesting to
note that the isotopic composition with respect to the other isotopes not used
for quantification has no effect on the measurement. (As and Mn are also
measured by the “network,” but are monoisotopic
and therefore not considered as part of this treatment).
Rather than assess the
isotopic composition of each individual sample, the uncertainty budget
developed for this measurement includes a component of uncertainty to recognise
that the isotopic composition of the sample may fall anywhere within the range
of natural variations, or the representative isotopic composition, whichever is
the larger range [17]. (The actual range of isotopic
composition in environmental samples may be considerably narrower [18].)
In practice, this assumption assigns δX=1 but imposes a relative uncertainty on this
value equal to the possible range of isotope abundances expected for the
isotopes used for quantification. Relatively little detail exists in the
literature on the isotopic composition of metals in ambient air particulates.
The vast majority of the work that has been published has been on Pb isotopic
composition, where the greatest variation is expected. Determination of isotope
ratios has been mostly used as a route to determining the origin of the Pb
sampled, particularly with regard to specific industrial processes or
long-range pollutant transport [14]. One study [19]
has examined the Pb isotopic composition in deposition in order to compare how
this changed before and after the closure of a local Pb mine. Others studies [20, 21]
have used the changing Pb isotope ratios in ambient particulate matter to
demonstrate the seasonal variation long-range transport of pollutants across
the Asian continent. Measured Pb isotope ratios have also been used as a route
to determine the changing origins of Pb emissions in an urban environment
during and after the phasing out of leaded petrol [22].
Cu and Zn isotopes ratios have been analysed near a large Zn refinery [23]
as a means of determining the origin of metallic ores, and Sr and Nd isotope
ratios have been used in the discrimination of emissions from various
industrial sources and traffic emissions, respectively [24].
In all cases, the observed ranges of the isotopic compositions fell well within
the natural ranges predicted [17] and these natural ranges have been used to construct
the uncertainty budget presented in this paper. Moreover, when the isotopic
composition of samples under consideration in this study has been measured
periodically, the abundance of the isotope used for quantitation has always
been well within these ranges as well [25].
The major contributions to the overall measurement uncertainty for the
determination of the mass concentration of metal in ambient air are shown in
Figure 1.
Relative contributions to the
standard uncertainty of the determination of metal mass concentration in
ambient air, as part of the UK Heavy Metals Monitoring Network. This example
shows a measurement with an overall expanded uncertainty at the 95% confidence
interval of approximately 20%. The changing uncertainty contribution from the
variation in the sample isotopic composition for the different metals measured
by the “network” is indicated by the additional lines and labelling on the
bottom bar.
As can be seen, the
contribution from the possible variation in the isotopic composition of the
sample is strongly dependent on the element being determined. This contribution
is very significant for Pb, significant for Hg, and Cd, less significant for
Cu, Zn and Pt, and negligible for Cr, Ni, Fe, and V. This is highlighted in Figure 2 which shows
the fraction of the overall measurement uncertainty contributed by uncertainty
in the isotopic composition of the sample being measured [25]. As expected from the data in Figure 1, this
contribution is very significant for Pb, notable for Hg, Cd and possibly Cu,
and insignificant for all other elements.
The fraction of the overall
measurement uncertainty for each metal contributed by the uncertainty in δX.
The summary uncertainty budget
presented in Figures 1 and 2 (and the measurement equation from which this has
been developed) is only valid if the abundance of the isotope used for
quantification in the calibration standards also falls within this range
allowed for the samples. If it does not,
then an additional uncertainty contribution, or a correction factor, may need
to be applied. This is not something that may be taken for granted since the
calibration standards may often have been prepared from isotopically enriched
pure materials, and must be measured in order to determine whether an increase
in the uncertainty estimate for the overall determination was required. The
results of the determination of the isotopic composition of the calibration
standards used are shown in Figures 3(a) and 3(b) which show that the majority
of isotopes demonstrated good agreement between the measured value and the
expected range of isotopic abundances. Table 1 highlights the level of this
agreement for the isotopes used for quantification.
Agreement between the measured abundances for the isotopes used for quantifications and the expected abundance ranges.
Isotope used for quantification
Agreement with predicted range
V50
No +0.01%
Cr52
No, −0.1%
Fe56
No, +0.4%
Ni60
No, +0.3%
Cu63
Yes
Zn66
Yes
Cd111
Yes
Pt194
Yes
Hg200
Yes
Pb208
Yes
(a) V, Cr, Fe, Ni, Cu, and Zn: comparison of
the measured relative isotopic abundance of the calibration standards (black
circles, with the grey bars representing the standard error of the mean)
against the expected range in natural, or representative, isotopic compositions
(whichever is the larger range) (black bars). The relative atomic mass number is
displayed for each isotope, with the boxed number being the isotope used for
the quantification of the samples. Values are normalised to the centre of the
natural (or representative) composition range for each isotope. The relative
abundance is displayed for each element in the separate plot beneath the main
chart. (b) Cd, Pt, Hg, and Pb: comparison
of the measured relative isotopic abundance of the calibration standards (black
circles, with the grey bars representing the standard error of the mean)
against the expected range in natural, or representative, isotopic compositions
(whichever is the larger range) (black bars). The relative atomic mass number
is displayed for each isotope, with the boxed number being the isotope used for
the quantification of the samples. Values are normalised to the centre of the
natural (or representative) composition range for each isotope. The relative
abundance is displayed for each element in the separate plot beneath the main
chart.
Where there is agreement with
the predicted range, we may assume that the uncertainty budget for the
measurement already covers the expected range of isotopic compositions for both
sample and calibration standard. Where Table 1 shows a lack of agreement, an
additional component of uncertainty equal to the discrepancy in agreement needs
to be added to justify the uncertainty statement. Given the large overall
measurement uncertainties reported for the measurement of metal mass
concentration in ambient air, these small additional uncertainty components are
unlikely to increase the overall uncertainty of the measurement significantly.
(An additional component of uncertainty has been added, rather than adding a
term into the uncertainty budget to correct for bias, since in all cases the
standard error of the mean of the isotopic composition measurements is greater
than the observed disagreement in all cases, and thus a bias correction term is
not justified.) Perhaps not surprisingly, Table 1 highlights that a disagreement
was found for the elements with the smallest predicted range of isotopic
compositions, where the accuracy of the isotopic measurement is more
critical.
A brief inspection of Figures 3(a) and 3(b) suggests that the bias between the measured isotopic abundance
and the centre of the expected range of isotopic abundances shows some tendency
to increase as the isotopic abundance decreases. This relationship is plotted
in Figure 4. Since the isotopes used for quantification are generally ones with
high abundances, little attention is usually paid to the less abundant
isotopes, and therefore, there may well be more bias in these measurements
owing to unresolved interferences, unresolved detector nonlinearity, and
instrument instabilities which have a larger proportional effect on the
measurement results than for high abundance species.
The relationship between the relative isotopic abundance of all the
isotopes considered and the relative difference between the measured isotopic
abundance and the centre of the expected range of isotopic abundance.
4. Conclusions
A summary of the main sources of
uncertainty encountered during the analysis of the mass concentration of metals
in ambient air as part of the operation of the UK Heavy Metals Monitoring
Network has been presented. It has been observed that the uncertainty
contribution from possible variations in the isotopic composition of the sample
depends on the element in question, but can be significant (as in the cases of
Pb, Cd, and Hg). The working curve method for the ICP-MS analysis of metals in
solution, with a low resolution, high throughput instrument measuring at one m/z ratio per element relies on the relative abundance of the isotopes under
consideration being the same in both the sample and the calibration solution.
Calculation of the uncertainty in this analysis assumes that the isotopic
composition variation within the sample and calibration solution is limited to
a defined range.
The results of the isotopic
analysis of these calibration standards have shown that the isotopic
composition of the calibration standards agrees with the expected range of
isotopic comparisons in the samples for all but four elements. In these cases,
additional uncertainty components were required to be added to the uncertainty
budget to account for this bias, although the increase in the overall uncertainty
of the measurement was not significant. It is interesting to note that the
isotopic composition with respect to the other isotopes not under consideration
not used for quantification has no effect, in theory, on the measurement.
The bias between the measured
isotopic abundance and the centre of the expected range of isotopic abundances
shows some tendency to increase as the isotopic abundance decreases. It has
been suggested that since the isotopes used for quantification are generally
ones with high abundance and little attention is paid to the lower abundance
isotopes, the bias in these measurements may well be due to unresolved
interferences, unresolved detector nonlinearity, and instrument instability
which have a larger proportional effect on the measurement results than for
high abundance species.
In future, it may be
expeditious to determine more rigorously the isotopic composition of samples
collected across the “network” so that the uncertainty contribution from δX for the elements where this is most
significant (in particular Pb) may be reduced, thereby reducing the overall
measurement uncertainty.
Acknowledgments
Useful discussions with Dr. Henrik Skov
from the National Environmental Research Institute (DMU, Denmark) are gratefully acknowledged. The UK Department for Innovation, University and Skills' funding of the National Measurement System Chemical and Biological Metrology Programme, and the UK Department for Environment, Food and Rural Affairs' funding of NPL's operation and management of the UK Heavy Metals Monitoring Network, are both gratefully acknowledged.
Ambient air pollution by As, Cd and Ni compounds, position paper. Brussels: working group on As, Cd and Ni compoundsEuropean Commission DG Environment, Brussels, Belgium, 2000Council Directive 1999/30/EC of 22 April 1999 relating to limit values for sulphur dioxide, nitrogen dioxide and oxides of nitrogen, particulate matter and lead in ambient air1999L1634160Council Directive 2004/107/EC of the European Parliament and of the Council of 15 December 2004 relating to arsenic, cadmium, mercury, nickel and polycyclic aromatic hydrocarbons in ambient air2005L023316BrownR. J. C.richard.brown@npl.co.ukThe use and abuse of limits of detection in environmental analytical chemistry20088796801BrownR. J. C.richard.brown@npl.co.ukYardleyR. E.MuhunthanD.Twenty-five years of nationwide ambient metals measurement in the United Kingdom: concentration levels and trends20081421–312714010.1007/s10661-007-9914-9EN 12341:1999 Air qualityDetermination of the PM10 fraction of suspended particulate matter. Reference method and field test procedure to demonstrate reference equivalence of measurement methodsCEN, Brussels, Belgium, 1999BrownR. J. C.richard.brown@npl.co.ukHoodD.daniel.hood@npl.co.ukBrownA. S.andrew.brown@npl.co.ukOn the optimum sampling time for the measurement of pollutants in ambient air20082008681471510.1155/2008/814715EN 14902:2005Ambient air quality-standard method for the measurement of Pb,Cd, As and Ni in the PM10 fraction of suspended particulate matterCEN, Brussels, Belgium, 2005BrownR. J. C.richard.brown@npl.co.ukYardleyR. E.BrownA. S.MiltonM. J. T.Sample matrix and critical interference effects on the recovery and accuracy of concentration measurements of arsenic in ambient particulate samples using ICP-MS200419570370510.1039/b316628gBrownA. S.andrew.brown@npl.co.ukBrownR. J. C.richard.brown@npl.co.ukA study of the critical uncertainty contributions in the analysis of PCBs in ambient air200820081417949810.1155/2008/179498NelmsS. M.QuételC. R.ProhaskaT.VoglJ.TaylorP. D. P.Evaluation of detector dead time calculation models for ICP-MS200116433333810.1039/b007913hSmithI. M.OnakunleF. O.XLGENLINE version 1.0, document CMSC/M/06/657NPL, Teddington, UK, 2007VoglJ.jochen.vogl@bam.deCharacterisation of reference materials by isotope dilution mass spectrometry200722547549210.1039/b614612kYipY.-C.ycyip@govtlab.gov.hkLamJ. C.-W.TongW.-F.Applications of lead isotope ratio measurements200827546048010.1016/j.trac.2008.02.011ISO Guide 98:1995Guide to the expression of uncertainty in measurement (GUM)ISO, Geneva, Switzerland, 1995MiltonM. J. T.martin.milton@npl.co.ukWielgoszR. I.Use of the international system of units (SI) in isotope ratio mass spectrometry200216232201220410.1002/rcm.836BöhlkeJ. K.de LaeterJ. R.de BièvreP.Isotopic compositions of the elements, 20012005341576710.1063/1.1836764CoplenT. B.tbcoplen@usgs.govBöhlkeJ. K.de BièvreP.Isotope-abundance variations of selected elements (IUPAC technical report)200274101987201710.1351/pac200274101987MunksgaardN. C.n_munksgaard@banks.ntu.edu.auParryD. L.Lead isotope ratios determined by ICP-MS: monitoring of mining-derived metal particulates in atmospheric fallout, Northern Territory, Australia19982171-211312510.1016/S0048-9697(98)00170-3HsuS.-C.schsu815@rcec.sinica.edu.twLiuS. C.JengW.-L.Lead isotope ratios in ambient aerosols from Taipei, Taiwan: identifying long-range transport of airborne Pb from the Yangtze Delta200640285393540410.1016/j.atmosenv.2006.05.003LeeC. S. L.LiX.-D.cexdli@polyu.edu.hkZhangG.LiJ.DingA.-J.WangT.Heavy metals and Pb isotopic composition of aerosols in urban and suburban areas of Hong Kong and Guangzhou, South China—evidence of the long-range transport of air contaminants200741243244710.1016/j.atmosenv.2006.07.035ZhengJ.jzheng@nirs.go.jpTanM.ShibataY.Characteristics of lead isotope ratios and elemental concentrations in PM10 fraction of airborne particulate matter in Shanghai after the phase-out of leaded gasoline20043881191120010.1016/j.atmosenv.2003.11.004MattielliN.RimetzJ.PetitJ.Zn-Cu isotopic study and speciation of airborne metal particles within a 5-km zone of a lead/zinc smelter19967018, supplement 1A40110.1016/j.gca.2006.06.808GeageaM. L.StilleP.pstille@illite.u-strasbg.frGauthier-LafayeF.MilletM.Tracing of industrial aerosol sources in an urban environment using Pb, Sr, and Nd isotopes200842369269810.1021/es071704cBrownR. J. C.ButterfieldD. M.GoddardS. L.MuhunthanD.WilliamsM.YardleyR. E.NPL Report AS20: Report to the Department of Environment, Food and Rural Affairs by the National Physical Laboratory: Annual Report for 2007 on the UK Heavy Metals Monitoring NetworkNPL, Teddington, UK, 2008