Romi
Nambiar
*ab,
Wolfgang
Müller
ab and
David
Evans‡
ab
aInstitute of Geosciences, Goethe University Frankfurt, Altenhöferallee 1, 60438 Frankfurt am Main, Germany. E-mail: Nambiar@geo.uni-frankfurt.de
bFIERCE, Frankfurt Isotope & Element Research Center, Goethe University Frankfurt, 60438 Frankfurt am Main, Germany
First published on 16th April 2024
Non-matrix-matched calibration of laser ablation ICPMS (trace/major) element data is a common quantification strategy. However, LA sampling is associated with downhole elemental fractionation, potentially causing inaccuracies if the magnitude of fractionation between the sample and reference material (RM) differs. Here, we estimate fractionation factors (FFs) for different elements (El) in a range of RMs relative to NIST SRM610/612 (FFEl/Ca-NIST) and evaluate element-specific corrections for downhole fractionation using these measured FFEl/Ca-NIST. Significantly different mean El/Ca values were observed before and after correction, particularly for the alkali elements (all RMs), and B, Fe, and Zn (some RMs), notably improving accuracy, especially for the alkali elements. In cases where this methodology does not result in an accuracy improvement, this may help identify underlying issues in reported/reference values for RMs, given that this phenomenon should be accounted for. Overall, we recommend considering routine assessment of FFs and applying a FF correction to enhance data quality.
When utilising LA as a sampling methodology, spot analysis (depth-profiling) remains the approach with the highest (vertical) spatial resolution. This method is capable of revealing elemental heterogeneity at sub-μm resolution, providing insights into diverse mineral formation processes from elemental zoning in igneous minerals13 to biominerals formed by marine calcifying organisms.9,14–17 However, during LA spot analysis, ablation-induced fractionation as a function of time, relatable in broad terms to element volatility/condensation temperature,18 usually referred to as ‘downhole fractionation’, can be a significant source of inaccuracy, particularly when the magnitude of this effect differs between the sample and calibration reference material. These phenomena of downhole fractionations were already reported in the early days of the LA pioneer period in the 1990s.19 Downhole fractionation is usually influenced by an analyte's geochemical affinity with the chosen internal standard, with elements from the same grouping in Goldschmidt's20 geochemical classifications characterised by similar behaviour.18 However, there are some exceptions; for example, the alkali elements often do not fractionate in the same way as other lithophile elements.21,22
Efforts to minimise potential sources of inaccuracy due to downhole fractionation have focused to a large degree on hardware solutions and related analytical approaches, with the broad shift to 193 nm ArF lasers in the Earth and Environmental Sciences resulting in substantially lower downhole element fractionation compared to 213 nm lasers.23 In addition: (1) the ablation pit geometry (depth-diameter ratio) has been shown to exert an influence on downhole fractionation, with higher depth/width aspect ratios largely resulting in more pronounced fractionation,24,25 and (2) downhole fractionation is less pronounced in a helium atmosphere compared to argon at a similar depth/diameter ratio.24,26 To overcome issues related to downhole fractionation, the most ideal approach would be using matrix-matched reference materials. Where these are unavailable, calibration or characterisation of this issue using in-house standards has also been proposed, e.g., by co-precipitation of elements into a CaCO3 matrix27 or adding standard solutions to a powder base.28,29 However, these approaches, in turn, require the (time-consuming) comprehensive characterization of the composition and homogeneity of any standard materials before they are used for this purpose.
Correction for downhole fractionation is common practise when high precision (permil-level) accuracy and precision are required, for example, for U–Pb geochronology.30–32 However, despite this being a widely-known issue, to our knowledge, only one study has attempted to correct for downhole fractionation associated with trace element determinations following non-matrix matched standardisation,33 using a model derived from the compositionally-matching reference material applied to the sample (i.e., matrix matched reference material is still ideally required). Fractionation factors may be as large as ±30% (Jochum et al., 2014) and depend on the analyte and sample matrix, such that failing to account for this issue could lead to important inaccuracies given that low-% level accuracy and precision is desirable when making trace element measurements.
Here, we determine fractionation factors for various common analytes, including the alkali elements (Li, Na, K), alkaline earth metals (Mg, Sr, Ba), metalloid (B), and transition metals (Fe, Zn) in a range of reference materials including the silicate MPI-DING glasses, GOR128-G, GOR132-G, KL2-G, and ATHO-G as well as nano-pellets of the carbonate reference materials JCp-1NP, JCt-1NP and MACS-3NP.34 We compare this information to the respective fractionation factor values in NIST SRM610/612, the most commonly used primary standards, and evaluate the contribution of downhole fractionation to inaccuracy. Finally, we show how measured fractionation factors, inherently available in the data collected in any case, can be used to directly correct for element-specific fractionation.
RESOlution laser ablation (LA) system (now Applied Spectra) | |
Wavelength | 193 nm |
Sampling mode | Spot (depth profiling) analysis |
Fluence/repetition rate | ∼6.3 J cm−2/3 Hz |
Ablation spot size/ablation time | 50 μm/60 s |
He flow | 300–400 mL min−1 |
N2 flow | 2.5–4.0 mL min−1 |
![]() |
|
Element XR Inductively Coupled Plasma Mass Spectrometer (ICPMS) - (ThermoFisher Scientific) | |
Mass resolution | Medium mass resolution |
Torch RF power | 1300–1380 W |
Sample cone/skimmer cone | Ni Jet cone/Ni H cone |
Sample gas flow | 0.86–1.00 L min−1 |
Auxiliary gas flow | 0.65–0.90 L min−1 |
Monitored elements (m/z) | 7Li, 11B, 23Na, 25Mg, 39K, 43Ca, 56Fe, 66Zn, 88Sr, and 138Ba (43Ca used as internal standard) |
Sensitivity measured on NIST SRM612 (60 μm; 6 Hz) in low mass resolution (LR) mode | >4.5 million cps on 238U |
ThO+/Th+ (m/z 248/232) in LR mode | <1% |
Doubly charged production rate (m/z 22/44) in LR mode | <2% |
FF = (X/43Ca0.5t–1t)/(X/43Ca0–0.5t), |
![]() | ||
Fig. 1 (A) The relationship between the downhole fractionation factors of several elements in different reference materials with respect to (wrt) NIST SRM612 versus NIST SRM610 shows a tight correlation close to a 1![]() ![]() |
Averaging over all measurements, a tight correlation close to a 1:1 line was observed between the fractionation factor relative to NIST SRM610 versus NIST SRM612 for all elements in a range of reference materials, as expected, given their very similar matrix. Only the Ba FFs showed a poorer correlation between the seven RMs (R2 = 0.51, p = 0.07; Fig. S1†), although all values are within error of a 1:
1 line. Overall, this is in good agreement with the earlier finding that the fractionation of elements in NIST SRM610 relative to NIST SRM612 using a 193 nm ArF laser were non-resolvable with the exception of a few volatile elements such as V, Zn, and Pb.22 Given that overall FF relative to NIST SRM610 and NIST SRM612 are comparable, in the next section, we use the NIST SRM610 calibrated data to demonstrate the utility of downhole fractionation correction but note that our results imply that a correction based on either standard would function in the same way.
This exercise (further) highlights that the volatile elements (Li and Zn in particular) are substantially offset from values of 1, i.e., characterised by substantially differing downhole fractionation between RMs. Specifically, the percentage differences between the RMs and NIST SRM610 fractionation factors were −3% to 18% for Li, the element displaying the highest degree of downhole fractionation. The FFEl/Ca_percent of Na and K ranged between −3 to 4% and −3 to 5%, respectively. In contrast, the alkaline earth metals showed the lowest FFEl/Ca_percent with all values ranging from −1% to 2%, which is one reason that non-matrix matched standardisation produces accurate data in these cases.22 In the case of B, the carbonate reference materials were characterised by greater offsets compared to NIST SRM610 (JCp-1NP = −8%, JCt-1NP = −6%, MACS-3NP = −8%) than the MPI-DING reference materials. In addition, Fe and Zn fractionation factors were offset by −2% to 3% and −15% to 2%, respectively. In the case of Zn, relatively lower degrees of offset (<2%) were observed for one carbonate and three MPI-DING glasses, in contrast to ATHO-G (−7%), JCp-1NP (−15%) and JCt-1NP (−10%). However, the RSD (n = 148) of the measured Zn fractionation factor was ∼50% for JCp-1NP and JCt-1NP, possibly resulting from the low [Zn] (<0.3 μg g−1) of these materials or inhomogeneity.
(1) The molar El/Ca value was standardised using the in-house Matlab script as described in the methods section (here denoted as El/Cacal). The El/Cacal values represent the primary standardised data i.e., before the downhole fractionation correction is applied.
(2) The fractionation factor (FF) was determined for each individual spot analysis of the sample relative to the NIST primary calibration standard (averaged over each individual session/single instrument run given that the variability of FFNIST <0.5% for each session), denoted as FFEl/Ca-NIST:
![]() | (1) |
(3) A downhole fractionation correction was then applied to the (primary) standardised El/Cacal values using FFEl/Ca-NIST estimated for each individual spot analysis:
El/Cacorrected = El/Cacal/FFEl/Ca-NIST | (2) |
In practical terms, this is approximately equivalent to regressing the analyte/internal standard count ratios versus time back to the respective y-axis intercept (as proposed in the case of U–Pb analysis45), which we avoid here simply because this approach may be sensitive to outliers. The resulting change in the accuracy of the measurements i.e., percent offset from reported values for the seven geological RMs, before and after applying the fractionation factor correction using eqn (2), measured over an ∼18 month period, is shown in Fig. 3. In order to verify if the shifts in El/Ca values following the correction for downhole fractionation is statistically different from the primary standardized El/Ca data, we utilise Student's t-test at the 95% confidence interval.
As expected, we find that analytes characterised by a greater degree of downhole fractionation relative to the NIST glass (Fig. 2) are associated with statistically significant changes (El/Cacalversus El/Cacorrected). In particular, the alkali elements resulted in significantly different mean values after the correction described above (eqn (1) and (2)). Specifically, Li/Ca, Na/Ca, and K/Ca accuracy were broadly improved in almost all RMs utilised here (Fig. 3, Table S1†). In contrast, applying a correction to the alkaline earth metals, which are characterised by the lowest degree of downhole fractionation relative to the primary NIST standard (Fig. 2), resulted in statistically indistinguishable mean values except for Mg in KL2-G (however, the long-term accuracy of Mg/Ca before and after fractionation correction were within 4%). In the case of B/Ca, we find a significant shift in accuracy for the carbonate RMs characterised by higher FFEl/Ca-NIST (Fig. 2) but not the MPI-DING glasses (Fig. 3, Table S1†). This correction for B/Ca resulted in a substantial improvement in accuracy in the case of JCt-1NP but a worsening of accuracy in the case of JCp-1NP and MACS-3NP. In addition, significant shifts in Zn/Ca and Fe/Ca accuracy resulted from the procedure outlined here, resulting in a worse apparent accuracy in the case of Zn/Ca in ATHO-G but an improvement in Zn/Ca accuracy in JCp-1NP and Fe/Ca accuracy in KL2-G.
![]() | ||
Fig. 3 Percent offset in El/Ca from the reported value of different geological reference materials before (blue) and after (orange) a downhole fractionation correction (eqn (1) and (2)). Solid symbols represent analyte/standard combinations characterised by statistically different population means before and after downhole fractionation correction based on a Student's t-test (95% CI); asterisks depict combinations for which the correction made no significant difference. Where data points are missing, analytes were below the limit of detection, and/or there is no reported value. The grey shaded region represents the uncertainty in the reported values (MPI-DING reference materials: 95% CL; carbonate reference materials: 1SD). |
The above results demonstrate a substantial improvement in accuracy for several analytes in the reference materials studied here (while exerting no appreciable impact on precision). Given that the resultant El/Ca following the proposed correction are significantly different for elements characterised by a greater degree of downhole fractionation, our results highlight the potential of applying this correction in improving data quality while performing non-matrix matched calibration. At the very least, this exercise serves to highlight specific cases in which such a correction is warranted, namely, the non-matrix-matched calibration of the alkali elements, and even in cases in which accuracy is not improved we nonetheless suggest that the approach taken here should be routinely applied and built into data reduction software. The rationale for this is that the reference values for many analytes in many reference materials are not well-characterised, a well-known issue with LA-ICPMS trace element analysis,40 and these have often been determined with a degree of circularity (i.e., the reported/reference values are derived to a large extent from LA measurements22). Given that downhole fractionation is demonstrably a complication for some of the analyte/matrix combinations reported here (e.g., B in carbonates), this is an analytical phenomenon that should be accounted for. The worsening of accuracy that we observe in some cases (Fig. 3) may, therefore, point towards an issue with the reported values, previously masked to a degree by a coincident (analytical) offset in the same direction. In the case that the reported/reference values were derived (mainly) from LA measurements, an apparent worsening in accuracy following our procedure may also result from the previous lack of a correction along the lines of that suggested here (for example, the available values of [B] in JCp-1NP, JCt-1NP, and MACS-3NP are based on LA-ICPMS/LA-MC-ICPMS; GEOREM database, version 35 (ref. 40)).
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ja00018h |
‡ Present address: School of Ocean and Earth Science, University of Southampton, Southampton, SO14 3ZH, UK. |
This journal is © The Royal Society of Chemistry 2024 |