International Journal of Analytical Chemistry

Volume 2014 (2014), Article ID 389125, 7 pages

http://dx.doi.org/10.1155/2014/389125

Research Article

Estimation of Measurement Uncertainties for the DGT Passive Sampler Used for Determination of Copper in Water

Jesper Knutsson, Sebastien Rauch, and Gregory M. Morrison

Department of Water Environment Technology, Division of Civil and Environmental Engineering, Chalmers University of Technology, 412 96 Göteborg, Sweden

Received 20 May 2014; Accepted 24 July 2014; Published 1 September 2014

Academic Editor: Mohammad R. Pourjavid

Copyright © 2014 Jesper Knutsson et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Diffusion-based passive samplers are increasingly used for water quality monitoring. While the overall method robustness and reproducibility for passive samplers in water are widely reported, there has been a lack of a detailed description of uncertainty sources. In this paper an uncertainty budget for the determination of fully labile Cu in water using a DGT passive sampler is presented. Uncertainty from the estimation of effective cross-sectional diffusion area and the instrumental determination of accumulated mass of analyte are the most significant sources of uncertainty, while uncertainties from contamination and the estimation of diffusion coefficient are negligible. The results presented highlight issues with passive samplers which are important to address if overall method uncertainty is to be reduced and effective strategies to reduce overall method uncertainty are presented.

1. Introduction

The overall goal of environmental management programs is to provide a framework for assessing environmental status, identifying problem areas, and to continuously assess quality indicators to ensure that those are within established acceptable limits which ensure a “good and nondeteriorating status.” One of the indicators of environmental quality outlined by the Water Framework Directive of the European Union is heavy metal concentration in water bodies, including Cu, Pb, Cd, and Ni [1]. There is therefore a stated need to measure and assess the environmental concentration of these metals. This should be done using a method that is representative and that provides comparable results across EU member states, though the directive does not specify what level of uncertainty is considered sufficient.

A passive sampler is a device used to collect a target analyte in situ, both in gaseous and liquid media. Recently, passive samplers have found increasing use in the determination of metals and organic contaminants in water [2–4]. However, measurement uncertainty, relatively little investigated, is a perceived limitation of passive sampling in comparison to the more conventional grab and automated bottle sampling procedures. The work presented here aims at characterizing and assessing the uncertainty associated with the determination of time-weighted concentrations of labile metal ions in freshwater using passive sampling.

A passive sampler for metal sampling is typically composed of a membrane filter, a diffusion layer gel, and a receiving phase placed in a sampler housing, like the DGT (diffusive gradients in thin films) technique (Figure 1). The DGT passive sampler was first described by Allan et al. [4] and since then the technique has been used in a wide range of applications and is one the most widely used passive sampler techniques for quantification and speciation of metals in aquatic environments. The analyte accumulates on the receiving phase as a result of the chemical affinity of the analyte for the solid receiving phase. The amount of analyte accumulated is proportional to the average concentration of labile analyte in the water, the time the sampler is exposed, and other aquatic environmental factors such as temperature and turbulence. After sampler retrieval and determination of the collected amount of metal, the average bulk concentration of metal can be calculated (see [3–5])

See [6].

Figure 1: A schematic rendering of a DGT passive sampler showing its principal components.

In (1), denotes the bulk concentration of the analyte in the water body, is the mass of the analyte accumulated on the sampler, is the diffusion coefficient of the analyte in water at 20°C, is the thickness of material diffusion layer (MDL, consisting of membrane filter and diffusion layer gel), is the diffusion coefficient of the analyte in the MDL, is the effective thickness of the diffusion boundary layer that is formed at the water-sampler interface, is the time of exposure, and is the effective sectional area of diffusion.

Although there has been some consideration of overall uncertainty in passive sampler measurements [7], there is no published study evaluating the components of this uncertainty. The identification of key components contributing to overall uncertainty can support the improvement of procedures based on passive samplings, as well as reducing potential concerns about performance and reliability [7].

2. Materials and Methods

For the purpose of this study, a simple case was assumed; a DGT passive sampler with characteristics listed in Table 1 was used to determine dissolved Cu concentration in water ( =  gl⁻¹). We note that the estimation of uncertainty resulting from metal-ligand interactions is out of the scope of this paper and Cu is therefore considered fully labile and present as Cu²⁺. In the absence of metal complexes, the time weighted average concentration can be derived using (1).

Table 1: Predefined passive sampler characteristics and environmental conditions used as a basis in the uncertainty calculations.

Parameter
Property/Value

Passive sampler

    Diameter
2 cm

    Diffusion layer
Acrylamide gel with APA cross-linker (APA2) [23]

    Cellulose nitrate membrane
135 µm thickness and 0.45 µm pore size

    Receiving phase
Resin-gel containing Chelex resin

Environmental conditions

    pH
7.5

    Water temperature
25°C/298 K

    Turbulence
Estimated

The uncertainty budget presented here was estimated for a generic passive sampler under predefined environmental conditions (Table 1). The characteristics of the passive samplers were chosen based on the characteristics of existing commercially available samplers (DGT Research Ltd.) and the availability of data. Similarly, environmental conditions were selected based on the availability of data for specific samplers. Although a number of passive sampler technologies have been described in the literature [8–11], the general methodology presented in this work should be applicable to estimate measurement uncertainty for a broad range of passive samplers, even if the specific conclusions for the passive sampler system assessed here does not necessarily hold true for other types.

Uncertainty in passive sampling is expected from all steps in the analytical process, including preparation of the samplers, deployment, analyte extraction, analysis, and estimation of diffusion rates and pathways. Overall, the estimation of uncertainties and the propagation of uncertainties were based on standard methodology [12]. Input data for the calculation were obtained from the literature and our own results, depending on availability. A cause and effect diagram was created to visualize the sources of uncertainty in the analytical chain when using a passive sampler to determine time weighted average bulk concentration (Figure 2). A list of relevant parameters (see Table 2) was identified from the cause and effect diagram and the model equation as a basis for the construction of the uncertainty budget.

Table 2: Parameters for which uncertainty is determined and respective units.

Parameter
Unit
Definition

m²
Effective area of diffusional cross-section

m² s⁻¹
Diffusion coefficient of the Cu²⁺ ion in the MDL

m² s⁻¹
Diffusion coefficient of the Cu²⁺ ion in water

g
Accumulated amount of Cu²⁺ determined from sample

g
Contamination determined from field blank

Recovery during the extraction phase

K
Temperature in bulk water phase

hours
Exposure time

m
Diffusional boundary layer thickness

m
Diffusional pathway thickness of the MDL

Figure 2: Cause and effect diagram describing the uncertainties associated with the determination of bulk concentration , using a passive sampler. Dashed arrows indicate parameters whose uncertainty contribution was included in another parameter. The dashed box shows the uncertainty from instrument determination of analyte. Uncertainty analysis of the ICP-MS technique has been performed previously [19, 22] and was therefore not treated separately in this paper.

3. Results and Discussion

3.1. Uncertainties in Analyte Accumulation

3.1.1. Diffusional Pathway

When deploying a prepared passive sampler, the fully labile metal ion (Cu²⁺) accumulates on the receiving phase and the accumulation rate is governed by diffusion across a diffusion boundary layer (DBL, see Figure 3), a membrane filter of known thickness (0.135 mm), and a gel layer of a known thickness (0.80 mm). No assessment was found of the uncertainty of , but a low uncertainty level was assumed for the combined membrane filter and gel layer ( mm) based on the authors judgement [12].

Figure 3: Schematic representation of the concentration gradient that forms over the diffusional pathway.

The DBL is the water layer closest to the passive sampler-water interface that is not affected by the mixing conditions in the bulk water phase. This measure is a representation of the effective DBL as this is neither evenly distributed layer across the surface nor a true unmixed layer but rather a velocity gradient. The effective thickness of the DBL is subject to uncertainty. The uncertainty can be reduced by deploying several devices with varying , as described by Zhang et al. [13], but this procedure increases the scope and cost of measurement considerably. Therefore in the hypothetical scenario presented here, the DBL thickness was estimated to be  mm, covering a wide range of flow regimes, from fast flowing water to slow moving lake epilimnion [6]. The diffusion coefficient of the metal ion Me²⁺ depends in turn on the water temperature and on which media it is diffusing in. The total accumulated amount () depends on the accumulation rate and the length of the exposure in time ().

3.1.2. Diffusion Coefficients

The diffusion coefficients and are usually determined experimentally in a separate experiment. The determination itself is associated with uncertainty, and results are typically reported without associated uncertainty. A typical relative uncertainty of diffusion coefficients has been reported in the range 1.3–6.4% [14, 15] and for the purpose of this assessment we will use a value of 6.42 × 10−10 m² s⁻¹, which is the diffusion coefficient of Cu²⁺ in APA2 gel (a polyacrylamide hydrogel containing 15 % vol acrylamide and 0.3% agarose-derived cross-linker) and the upper value in the uncertainty interval, that is, 6.4% [15]. The diffusion coefficient of Cu²⁺ in water () is reported to be 1.14 times larger at 7.30 × 10−10 m² s⁻¹ [15]. For the purpose of this paper that same relative uncertainty was applied to both and . It should be noted that effective diffusion coefficients may also be significantly affected in low ionic strength solutions (<1 mM).

The diffusion coefficient depends on water temperature as described by the Stokes-Einstein equation:

where is the Boltzmann constant (m² kg s⁻² K⁻¹), is the temperature (K), is the viscosity of the medium (kg s−1 m⁻¹), and is the spherical diameter of the diffusing particle.

The uncertainty introduced from variability of was analysed (). The uncertainty in the experimental determination of was also estimated. The standard uncertainty in from uncertainty in water temperature was calculated to be 0.06 × 10⁻¹⁰ m² s⁻¹, and the combined standard uncertainty from the determination of and temperature was calculated through summation in quadrature to be × 10⁻¹⁰ m² s⁻¹. A similar treatment of resulted in  m² s⁻¹.

3.1.3. Effective Area

The effective area of the section through which diffusion occurs has been reported to be somehow larger than the nominal area due to lateral diffusion; that is, diffusion occurs in three dimensions [6, 7]. Warnken et al. report that the radius of the effective diffusion window is  cm and also note that the gel disc diameter had shrunk on average 0.12 cm () during drying prior to determination of the radius [16]. No estimate on uncertainty was given for this measure, so a 0.05 cm uncertainty was assumed based on the number of significant figures reported, and a rectangular distribution was selected due to the lack of information on the measurement.

Summation in quadrature was used to combine the uncertainties from the determination of effective radius and the estimation of the shrinkage in order to calculate the total uncertainty associated with the effective area [17]. The divisor was used to get the standard uncertainty of the shrinkage because of the assumed rectangular distribution, followed by summation in quadrature:

The combined uncertainty of the effective radius was calculated to be 0.0449 cm, making the effective radius of the sampler  cm. Using the derivative of the circle area function to calculate the uncertainty of the effective area, , gave the value  cm².

3.2. Uncertainties in Determination of Mass

3.2.1. Preparation and Handling

During preparation, transport, storage, and handling of the passive sampler devices there is a risk of contamination. The best assessment of the uncertainty from these sources comes from the use of field blanks [18]. The field blanks can be used to correct for contamination issues. We have estimated during field trials that the associated relative uncertainty resulting from contamination is typically in the order of 24% for passive sampler devices, with field blank values of  ng Cu²⁺ () (unpublished data).

3.2.2. Extraction

The analyte (Cu²⁺) is subsequently extracted from the receiving phase using a small volume of nitric acid. The recovery factor, , has been reported previously () [5]. The uncertainty was reported as an interval, and therefore a rectangular distribution was assumed.

3.2.3. Analysis/Determination

The resulting extract is diluted to a suitable volume concentration before analysis by a selected analytical technique. Inductively coupled plasma-mass spectrometry (ICP-MS) is widely used for the determination of trace metal concentrations in environmental samples and therefore, we estimate uncertainty for ICP-MS analysis in this paper. The ICP-MS instrument is calibrated using calibration standards prepared from certified standard solutions.

Generally, the analytical procedure using ICP-MS is subject to known and unknown interferences of which some can be compensated for, while others may persist, depending on specific instrument capabilities [18]. Furthermore, instrument drift, stability of stock solutions, and density of stock solutions will contribute to uncertainty [18] and the uncertainty budget of the instrumental analysis is a comprehensive topic in its own right. A simplified view is given in Figure 2 to highlight the importance of the analytical step. However, instrument performance and the typical uncertainty of the method have been addressed elsewhere [19] and are not repeated here. The reported standard relative uncertainty for Ni solutions containing 10 ng g⁻¹ or more was 7.5%, which was used for the calculations in this paper.

The estimated accumulated mass and mass on blank samples was determined using ICP-MS and then corrected for by the recovery factor according to

Using the rule for uncertainty propagation in quotients the estimate for becomes  g (see Table 3). A similar treatment of resulted in  g.

Table 3: Uncertainty budget for showing relative uncertainties for the variables and the combined standard uncertainty.

Symbol
Source of uncertainty
Type*
Standard uncertainty
Distribution
Divisor
Relative uncertainty

Estimated mass from ICP-MS analysis
A
 g
Normal
1
0.008

Recovery factor
B
0.0293
Rectangular

0.064

Uc (M)
Combined standard uncertainty
A
 g
Normal

0.038

Note: type of uncertainty refers to types A and B, using standard vocabulary for statistically evaluated uncertainty (A) and uncertainty evaluated by other methods (B).

3.3. Total Combined Uncertainty of the Passive Sampler Measurement

To estimate the combined standard uncertainty of the bulk concentration , the relation in the model equation (1) was used. Since it was a mixed expression, the rule of uncertainty propagation states that the combined uncertainty can be calculated using

See [20].

This means that the combined uncertainty is equal to the root square sum of the partial derivatives of the variables. However, it is also possible to derive a numerical solution as suggested by Kragten [21]. The approximation derived from this numerical method assumes linearity and small values of relative uncertainty, . While this is not always the case, the accuracy of the solution is still acceptable for most practical purposes [21].

A summary of the quantities and the associated standard uncertainties is presented in Table 4.

Table 4: Quantities, nominal values, and their associated uncertainty used in this work.

Quantity
Value
Standard uncertainty
Comment

3.66 cm²
0.30 cm²
See previous section and [16]

 m²/s
 m²/s
Empirical value [15]

 m²/s
 m²/s
Empirical value [15]

 g
 g
Observation

 g
 g
Observation

0.793
0.051
Observation [5]

168 h
0.3 h
Covers the time it takes to deploy and retrieves 5 passive samplers

25°C/298 K
4 K
Standard deviation of the measured temperature

 m
 m
Estimate [16]

 m
 m
Estimate

During calculations values were not rounded to avoid the introduction of additional uncertainty. The output of the numerical treatment of combined uncertainties can be seen in Table 5. The measurement output with associated uncertainty was =  μg l⁻¹. Using a coverage factor the result was instead  μg l⁻¹ (confidence interval ≈ 95%), or a relative uncertainty of 7.6% at .

Table 5: Uncertainty budget for determination of time weighted average concentration of Cu²⁺ in water using a DGT passive sampler.

Symbol
Source of uncertainty
Type
Standard uncertainty
Distribution
Divisor
(M) μg L⁻¹

Determination of accumulated mass
A
 g
Normal
1
0.49

Determination of contamination
A
 g
Normal
1
0.02

Diffusion coefficient in water
A
 m²/s
Normal
1
0.16

Thickness of material diffusion layer (MDL)
B
 m
Rectangular

0.33

Diffusion coefficient in MDL
A
 m²/s
Normal
1
0.16

Diffusion boundary layer
B
 m
Rectangular

0.29

Time
B
624 s
Rectangular

0.01

Effective area
A
 m²
normal
1
0.69

Uc ()
Combined standard uncertainty

Normal

0.98

Uc ()
Expanded standard uncertainty

Normal ()

1.95

When plotting the relative standard uncertainties of the components graphically (Figure 4) it is obvious that the largest uncertainty was introduced from the effective cross-sectional area estimate (). The combined estimated uncertainties resulting from the determination of the lateral diffusion round the edges and the shrinkage of the gel resulted in an uncertainty that largely affects the end result, as it accounts for nearly 50% of the total uncertainty (see Figure 4). Uncertainties from the estimation of account for roughly 25% of the total uncertainty, with the most significant factor being the estimation of extraction recovery.

(a)

(b)

Figure 4: Relative standard uncertainty (a) and percentage of total uncertainty (b) for the variables in the model equation.

A sensitivity analysis shows that halving the uncertainty for the effective radius and shrinkage in the determination of would reduce the contribution of to the overall uncertainty to roughly 20% ( μg l⁻¹ or 6.0% relative overall uncertainty). Similarly, a reduction in the uncertainty in the recovery factor, , by 50%, would reduce the contribution from to overall uncertainty from 25% to approximately 9% ( μg l⁻¹ or 6.8% relative overall uncertainty). On the other hand, an increase in the uncertainty for diffusion layer thickness from 0.05 mm to 0.2 mm would result in  μg l⁻¹ or 12.2% relative overall uncertainty. This is a significant increase in overall uncertainty and illustrates the sensitivity of the method to inconsistencies in the gel-membrane layer interface. Furthermore, the effects of uncertainty changes in DBL and diffusive layer thickness are shown in Table 6.

Table 6: Results from sensitivity analysis, showing the effect on total uncertainty of the passive sampler measurement from reductions in uncertainty of selected parameters.

Parameter
Change in uncertainty
Result on total uncertainty

Effective area,
50% reduction
Reduction from 7.6% to 6.1% in overall relative uncertainty

Recovery factor,
50% reduction
Reduction from 7.6% to 6.9% in overall relative uncertainty

Diffusion boundary layer,
From 0.05 mm to 0.014 mm standard uncertainty
Reduction from 7.6% to 7.3% in overall relative uncertainty

Diffusion pathway thickness
50% reduction
Reduction from 7.6% to 7.3% in overall relative uncertainty

Diffusion pathway thickness
4 times increase
Increase from 7.6% to 12.2% in overall relative uncertainty

The sensitivity analysis shows that overall method uncertainty can be significantly reduced by addressing the proper sources of uncertainties and also that deterioration in diffusion layer consistency can have significant negative effects on overall method uncertainty.

4. Conclusion

An uncertainty analysis was performed for passive sampling of a metal ion in water to highlight critical steps in the method and to identify key factors for potential improvement. In the analysis performed here the uncertainty of the effective cross-sectional diffusion area was identified as the main contributor to overall uncertainty. Uncertainties in analyte recovery and material diffusion ranked second and third, respectively. An improvement in the estimation of was found to be an important step toward achieving a reduction in uncertainty in passive sampling. Optimization of the extraction procedure will provide a further reduction in overall uncertainty.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

European Commission, 2000/60/EC, 2000

H. Zhang and W. Davison, “Direct in situ measurements of labile inorganic and organically bound metal species in synthetic solutions and natural waters using diffusive gradients in thin films,” Analytical Chemistry, vol. 72, no. 18, pp. 4447–4457, 2000.

L. B. Blom, G. M. Morrison, J. Kingston et al., “Performance of an in situ passive sampling system for metals in stormwater,” Journal of Environmental Monitoring, vol. 4, no. 2, pp. 258–262, 2002.

I. J. Allan, J. Knutsson, N. Guigues, G. A. Mills, A. Fouillac, and R. Greenwood, “Evaluation of the chemcatcher and DGT passive samplers for monitoring metals with highly fluctuating water concentrations,” Journal of Environmental Monitoring, vol. 9, no. 7, pp. 672–681, 2007.

H. Zhang and W. Davison, “Performance characteristics of diffusion gradients in thin films for the in situ measurement of trace metals in aqueous solution,” Analytical Chemistry, vol. 67, no. 19, pp. 3391–3400, 1995.

W. Davison and H. Zhang, “Progress in understanding the use of diffusive gradients in thin films (DGT) back to basics,” Environmental Chemistry, vol. 9, no. 1, pp. 1–13, 2012.

C. Miège, S. Schiavone, A. Dabrin et al., “An in situ intercomparison exercise on passive samplers for monitoring metals, polycyclic aromatic hydrocarbons and pesticides in surface waters,” Trends in Analytical Chemistry, vol. 36, pp. 128–143, 2012.

W. Davison and H. Zhang, “In situ speciation measurements of trace components in natural waters using thin-film gels,” Nature, vol. 367, no. 6463, pp. 546–548, 1994.

W. G. Brumbaugh, J. D. Petty, T. W. May, and J. N. Huckins, “A passive integrative sampler for mercury vapor in air and neutral mercury species in water,” Chemosphere—Global Change Science, vol. 2, no. 1, pp. 1–9, 2000.

J. K. Kingston, R. Greenwood, G. A. Mills, G. M. Morrison, and L. B. Persson, “Development of a novel passive sampling system for the time-averaged measurement of a range of organic pollutants in aquatic environments,” Journal of Environmental Monitoring, vol. 2, no. 5, pp. 487–495, 2000.

J. Rozemeijer, Y. van der Velde, H. de Jonge, F. van Geer, H. P. Broers, and M. Bierkens, “Application and evaluation of a new passive sampler for measuring average solute concentrations in a catchment scale water quality monitoring study,” Environmental Science and Technology, vol. 44, no. 4, pp. 1353–1359, 2010.

S. L. R. Ellison and A. Williams, “Eurachem/CITAC guide: Quantifying Uncertainty in Analytical Measurement,” 2012

H. Zhang, W. Davison, R. Gadi, and T. Kobayashi, “In situ measurement of dissolved phosphorus in natural waters using DGT,” Analytica Chimica Acta, vol. 370, no. 1, pp. 29–38, 1998.

E. R. Unsworth, H. Zhang, and W. Davison, “Use of diffusive gradients in thin films to measure cadmium speciation in solutions with synthetic and natural ligands: comparison with model predictions,” Environmental Science and Technology, vol. 39, no. 2, pp. 624–630, 2005.

S. Scally, W. Davison, and H. Zhang, “Diffusion coefficients of metals and metal complexes in hydrogels used in diffusive gradients in thin films,” Analytica Chimica Acta, vol. 558, no. 1-2, pp. 222–229, 2006.

K. W. Warnken, H. Zhang, and W. Davison, “Accuracy of the diffusive gradients in thin-films technique: diffusive boundary layer and effective sampling area considerations,” Analytical Chemistry, vol. 78, no. 11, pp. 3780–3787, 2006.

I. J. Allan, J. Knutsson, N. Guigues, G. A. Mills, A. Fouillac, and R. Greenwood, “Chemcatcher and DGT passive sampling devices for regulatory monitoring of trace metals in surface water,” Journal of Environmental Monitoring, vol. 10, no. 7, pp. 821–829, 2008.

S. Nelms, Inductively Coupled Plasma Mass Spectrometry Handbook, Blackwell, 2005.

V. J. Barwick, S. L. R. Ellison, and B. Fairman, “Estimation of uncertainties in ICP-MS analysis: a practical methodology,” Analytica Chimica Acta, vol. 394, no. 2-3, pp. 281–291, 1999.

J. R. Taylor, An Introduction to Error Analysis, University Science books, Sausalito, Calif, USA, 1997.

J. Kragten, “Tutorial review. Calculating standard deviations and confidence intervals with a universally applicable spreadsheet technique,” The Analyst, vol. 119, no. 10, pp. 2161–2165, 1994.

P. Evans and B. Fairman, “High resolution ID-ICP-MS certification of an estuary water reference material (LGC 6016) and analysis of matrix induced polyatomic interferences,” Journal of Environmental Monitoring, vol. 3, no. 5, pp. 469–474, 2001.

H. Zhang and W. Davison, “Diffusional characteristics of hydrogels used in DGT and DET techniques,” Analytica Chimica Acta, vol. 398, no. 2-3, pp. 329–340, 1999.