Articles | Volume 6, issue 3
https://doi.org/10.5194/gchron-6-465-2024
https://doi.org/10.5194/gchron-6-465-2024
Short communication/technical note
 | 
08 Aug 2024
Short communication/technical note |  | 08 Aug 2024

Technical note: RA138 calcite U–Pb LA-ICP-MS primary reference material

Marcel Guillong, Elias Samankassou, Inigo A. Müller, Dawid Szymanowski, Nathan Looser, Lorenzo Tavazzani, Óscar Merino-Tomé, Juan R. Bahamonde, Yannick Buret, and Maria Ovtcharova
Abstract

A promising primary reference material for U–Pb laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) carbonate dating is analysed and reported here. The new reference material (RM) is a botryoidal cement (C1) from sample RA138. The sample was collected in outcrop strata of mid-Carboniferous age in northern Spain near La Robla, and multiple aliquots have been meticulously prepared for distribution. RA138 is characterized by variable U/Pb ratios (from  1 to  19) and a relatively high and homogeneous U content ( 4 ppm). This material exhibits a low age uncertainty (0.2 %, 2 s; unanchored; ID-TIMS), allowing for the establishment of a well-defined isochron, particularly when anchored to the initial Pb ratio using LA-ICP-MS. Isotope dilution thermal ionization mass spectrometry (ID-TIMS) analyses of micro-drilled C1 cement (17 sub-samples) produce a lower-intercept age of 321.99 ± 0.65 Ma, an initial 207Pb/206Pb ratio of 0.8495 ± 0.0065, and a mean square weighted deviation (MSWD) of 5.1. The systematic uncertainty of 1.5 % observed in repeated LA-ICP-MS analyses challenges previous estimations of 2 %–2.5 % based on repeated analyses of ASH-15D and JT using WC-1 as primary reference material, underscoring the precision and reliability of RA138 for U–Pb dating applications.

1 Introduction

Recent years have seen a growing interest in the application of laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) to carbonate U–Pb geochronology. This method provided an accessible way to determine carbonate crystallization age, benefiting from the availability of the well-characterized, matrix-matched calcite primary reference material (RM) WC-1 (Roberts et al., 2017). The establishment of an accepted methodology for corrections and age calculations (Coogan et al., 2016; Li et al., 2014; Ring and Gerdes, 2016; Roberts et al., 2017) has resulted in an increasing number of publications (Hoareau et al., 2021; Kylander-Clark, 2020; Rasbury et al., 2023; Giorno et al., 2022) spanning a diverse range of applications including palaeoclimate (Chaldekas et al., 2022; Drake et al., 2017; Woodhead and Petrus, 2019; Gulbranson et al., 2022), hydrothermal (Brigaud et al., 2020; Burisch et al., 2017; MacDonald et al., 2019; Mottram et al., 2020; Piccione et al., 2019; Tavazzani et al., 2024), and tectonic (Hansman et al., 2018; Looser et al., 2021; Nuriel et al., 2017; Weinberger et al., 2020; Nuriel et al., 2019) activities and pedogenic and diagenesis processes (Elisha et al., 2021; Godeau et al., 2018; Mangenot et al., 2018; Methner et al., 2016; Scardia et al., 2019; Sindern et al., 2019).

Several additional RMs have been dated and distributed within the scientific community. However, their suitability as primary RMs is either limited due to their young age and low Pb count rates (Nuriel et al., 2021) or poses limitations due to their exhaustion (Guillong et al., 2020). Other published RMs, such as B6 (Pagel et al., 2018), lack independent age constraints from primary measurement techniques such as isotope dilution thermal ionization mass spectrometry (ID-TIMS). A comprehensive overview of available and proposed carbonate RMs can be found elsewhere (Wu et al., 2022).

While WC-1 is the most commonly used primary RM, it is not the most suitable due to its heterogeneity, and the scientific community is still in search of a better reference material (Roberts et al., 2017). The significant uncertainty in the ID isotope ratio mass spectrometry (ID-IRMS) age of WC-1 (254.4 ± 6.4 Ma) imposes a minimal uncertainty of > 2.5 % on any obtained carbonate LA-ICP-MS age. In addition, the very high mean square weighted deviation (MSWD) of 1069 (when the isochron is anchored to the LA-ICP-MS initial 207Pb/206Pb ratio of 0.85 ± 0.04), suggests that WC-1 most likely exhibits natural heterogeneities. Consequently, the use of WC-1 as a primary reference material imposes constraints on achievable precision and accuracy.

As part of a collaborative effort to use carbonate U–Pb LA-ICP-MS geochronology of syn-depositional botryoidal cements as a stratigraphic tool, we dated several samples from carbonate debris flows sourced from a series of steep and high-relief carbonate platforms that developed during the Carboniferous in the marine foreland basin of the Cantabrian Zone (CZ) in northern Spain. Some of these samples have yielded excellent isochrons, suggesting that they may be suitable as natural RMs for carbonate U–Pb dating by LA-ICP-MS. Here, we present a detailed characterization of one of these samples, RA138, which we identify as a promising RM for the growing LA-ICP-MS U–Pb geochronology community. Besides its suitability for producing well-defined isochrons, RA138 was chosen because of the available amount of material (> 6 kg), the ease of access to the location for any future sampling, and the abundance of botryoidal cement in the outcrop sample.

2 Sample

RA138 denotes that the sample was collected at metre 138 of section A, recorded at La Robla quarry near the village of La Robla in northern Spain. The precise geographical coordinates are 42.808355° N and 5.687950° W (EPSG:4258 (ETRS89)). A substantial sample exceeding 6 kg in weight was extracted from a microbial boundstone lithoclast characterized by the presence of exceptionally well-preserved botryoidal cements. This lithoclast was embedded within a matrix-supported carbonate debris flow interbedded within uppermost Serpukhovian limestone strata. Comprehensive discussions regarding the samples and the regional geology will be presented in a forthcoming publication. During the analytical process, several thin sections were prepared, and several blocks measuring approximately 30×20×5 mm were cut for subsequent polishing, cold-cathodoluminescence imaging, and the identification of botryoidal cement generations. Further details on selected individual blocks are available in Figs. S1–S5 in the Supplement.

3 Methods

LA-ICP-MS stands as the primary methodology for which RA138 serves as the designated RM. Furthermore, LA-ICP-MS was employed for comprehensive characterization, including analyses of trace element distribution and U–Pb age determination. Independently, two laboratories utilized ID-TIMS analyses as a supplementary method to establish the reference U–Pb isotopic ratios.

In tandem with these analytical techniques, cathodoluminescence (CL) imaging was acquired using an Olympus polarizing microscope equipped with a CITL 8200 Mk5-2 stage, operated at 15 kV and 130 µA, and an Olympus DP74 camera was used. Combined with microscope images of polished sections, this approach revealed growth zoning and the different phases present in the sample.

3.1 LA-ICP-MS

LA-ICP-MS analyses were conducted using a Resonetics (now Applied Spectra) RESOlution 155HR LA system coupled to a Thermo Element XR single-collector sector field ICP-MS. This setup included a high-capacity interface pump to improve sensitivity, along with the addition of small quantities of nitrogen (2 mL min−1) in the makeup argon gas flow (Wu et al., 2022). The same RA138 sample section was analysed during a series of individual sessions (n=23), aiming to demonstrate the repeatability and homogeneity of the material. The LA parameters employed for U–Pb dating were consistent with previous methodologies (Guillong et al., 2020) utilizing a 110 µm static spot ablation at 5 Hz and an energy density of 2 J cm−2. Key parameters are concisely outlined in Table S1 in the Supplement. Data reduction adhered to established procedures (Guillong et al., 2020; Roberts et al., 2017), with partial processing accomplished in Iolite 4 (Paton et al., 2011; Petrus and Kamber, 2012). This involved the selection of integration intervals, gas blank subtraction, and downhole fractionation correction (Paton et al., 2010), utilizing WC-1 and the UcomPbine (Chew et al., 2014) data reduction scheme, without implementing common Pb correction on the RM.

Further data treatment occurred in Excel, where the downhole-corrected 206Pb/238U ratio, 207Pb/206Pb ratio, raw counts per second, and error correlation were exported from Iolite. This treatment encompassed drift correction for the 206Pb/238U ratio, referencing NIST (National Institute of Standards and Technology) 614 and NIST 612; 207Pb/206Pb ratio calculation (ratio of mean); and normalization to the known NIST 614 ratio [0.8704 GeoReM preferred value (Jochum et al., 2005)]. Corrected 207Pb/206Pb ratios and intermediate 238U/206Pb ratios were subsequently plotted using IsoplotR (Vermeesch, 2018) in a Tera–Wasserburg (TW) diagram, enabling the calculation of a discordia model 1 anchored to a fixed initial Pb value (0.85 for WC-1 and 0.8495 ± 0.0065 for RA138) to determine a lower-intercept age.

To account for the known age of the primary reference material (RM), a correction factor was calculated and applied uniformly to all unknown 238U/206Pb ratios in each session, thus ensuring accuracy for the 238U/206Pb ratios of all spot analyses. The accuracy of this data reduction methodology was rigorously tested across all sessions using multiple secondary RMs, typically JT (Guillong et al., 2020), B6 (Pagel et al., 2018), and ASH-15D (Nuriel et al., 2021). Measurement uncertainties arising from the ratios, as well as uncertainties associated with drift correction and the 238U/206Pb correction factor, were quadratically propagated.

In addition to the dating analyses, we conducted an in-depth investigation of trace elements within the various phases using a laser spot diameter of 30 µm. NIST 612 served as the reference material, and stoichiometric calcium (40 wt % Ca) was employed as an internal standard for ablation yield correction. Trace element maps were acquired in scanning mode to examine the distribution of both trace elements (TEs) and U/Pb ratios.

For low-resolution mapping purposes, a square spot size of 91 µm, a repetition rate of 10 Hz, an energy density of 2 J cm−2, and a scanning rate of 50 µm s−1 were employed. A total of 32 isotopes were measured with a sweep time of 0.585 s, and data reduction was performed using Iolite 4.8.3. NIST 614 was utilized as the reference material, and a 40 % weight of calcium (Ca) served as the internal standard. All maps are shown in Figs. S6 and S7.

For high-resolution maps, a Teledyne Iridia LA system equipped with a Cobalt fast-washout ablation cell was coupled to an Agilent 8900 ICP-MS equipped with a quad-lock system (Norris et al., 2021) at the National History Museum in London. A square spot size of 20 µm, a repetition rate of 298.1 Hz, an energy density of 2 J cm−2, and a scanning rate of 424.3 µm s−1 were employed. Nine masses (24Mg, 43Ca, 47Ti, 55Mn, 57Fe, 86Sr, 89Y, 206Pb, and 238U) were analysed with a sweep time of 47 ms. Washout was tested to be < 5 ms to 10 % of the maximum signal from a single laser pulse. Data were processed and maps were constructed using Iolite 4 (Paton et al., 2011).

3.2 ID-TIMS

RA138 was analysed by isotope dilution thermal ionization mass spectrometry (ID-TIMS) in two sessions. All samples were purified at the University of Geneva and subsequently analysed by TIMS either in Geneva (session 1) or at both Geneva and ETH Zürich (session 2). The botryoidal cement C1 was drilled to produce 1–4 mg fine powder from polished rock surface using a Dremel drill with a drill bit diameter of 0.8 mm, allowing sample spots with approximately 1 mm diameter. The sample powders were collected in 12 mL PMP (polymethylpentene) beakers, dissolved in 100 µL 6 M HCl, and equilibrated with 5–9 mg of the EARTHTIME 205Pb–233U–235U(ET535) tracer solution (Condon et al., 2015; McLean et al., 2015) during 30 min on an 80 °C hotplate. The solutions were then dried down and re-dissolved in 1 M HBr. Pb was purified in an HBr ion exchange chemistry using 50 µL of AG1-X8 anion resin (Cl form; 200–400 µm mesh). The U cuts were converted to chloride form and then purified through a second pass in the AG1-X8 columns (rinsed in 3 and 6 M HCl; eluted in water). Due to the high cation concentration in carbonates, the U fractions were then further put through a separate set of 50 µL columns filled with RE resin (TrisKem B50-S; 50–100 µm). Pb and U fractions were dried down with trace H3PO4 (0.02 M) in separate 7 mL PFA (perfluoroalkoxy) beakers.

Uranium and Pb isotopic ratios were analysed on two thermal ionization mass spectrometers (TIMS), a Thermo Triton at the University of Geneva and a Thermo Triton Plus at ETH Zürich, both equipped with 1013Ω amplifiers. U and Pb aliquots were loaded at the respective labs on separate, zone-refined, outgassed Re filaments in 1 µL of silica gel emitter (modified after Gerstenberger and Haase, 1997). In both cases, Pb isotope measurements were done in static mode, either using Faraday cups for all Pb masses (ETH Zürich) or combining Faraday cups for 205–208Pb with an axial SEM for 204Pb (Geneva) (von Quadt et al., 2016). Instrumental mass fractionation of Pb isotopes was corrected with factors specific to each detector setup and derived from a compilation of mass fractionation factors measured in Pb isotopic standards and double-spiked unknowns. U was measured as UO2 in static mode, with oxide interferences corrected either using in-run 18O/16O determined with mass 272 (Wotzlaw et al., 2017) or assuming a value of 0.00205 ± 0.00004. Mass fractionation was corrected with the known spike 233U/235U value of 0.9950621 and assuming a sample 238U/235U of 137.818 ± 0.045 (2σ) (Hiess et al., 2012). The chosen value could introduce a bias if RA138 had a significantly different 238U/235U (Chen et al., 2021); however, the magnitude of the resulting age shift would be negligible.

Raw data were evaluated using the Tripoli software (Bowring et al., 2011) and reduced using the algorithms of Schmitz and Schoene (2007). Based on total procedural blank measurements at the University of Geneva, 3.5 pg of common Pb was assigned to laboratory blank and corrected with its long-term isotopic composition. Sample ages were evaluated using the 2D Tera–Wasserburg diagram in IsoplotR (Vermeesch, 2018), using the decay constants of Jaffey et al. (1971). The U/Pb ratios were not corrected for initial Th, Pa, or U disequilibrium.

https://gchron.copernicus.org/articles/6/465/2024/gchron-6-465-2024-f01

Figure 1Characterization of the different phases found in RA138. (a) Microscope image with four different phases, as well as spot locations for analyses shown in panel (c). (b) CL image of the square region with spot locations. (c) Tera–Wasserburg diagram of results obtained from the spots shown in panels (a) and (b) and matrix analyses. (d) Results of trace element concentration analysed in the four different phases.

Download

4 Results

4.1 Textures

Sample RA138 exhibits at least four distinct phases of precipitation, as illustrated in Fig. 1a (thin section image), Fig. 1b (CL-image), and Fig. 2 (LA-ICP-MS compositional maps):

  1. Matrix. Visually and under CL, the matrix has a patchy appearance and consists of clotted peloidal micrite with scarce bioclasts. Millimetre-sized internal porosities are lined by thin isopachous rims of fibrous cements (< 100 µm thick) and incipient botryoidal cement fans. Remaining porosity is occluded by blocky calcite cements.

  2. Cement C1 (botryoidal). This phase is characterized by a dark colour and a fibrous appearance observed in sections with concentric growth bands perpendicular to the direction of cement growth. Under CL, it appears dark (low or no luminescence) and relatively homogeneous.

  3. Replacement cement (C2). This phase replaces cement C1 and displays a fine-grained texture with a brighter appearance, both visually and under CL (high or medium luminescence).

  4. Radiaxial fibrous and sparry cement (RFSC3). This phase exhibits a white to transparent appearance and is very bright under CL (high luminescence). It is distinguished by the larger crystals exceeding 200 µm in size.

https://gchron.copernicus.org/articles/6/465/2024/gchron-6-465-2024-f02

Figure 2CL image and trace element maps of RA138. (a) CL image. (b) Mg concentration map showing the low Mg content of C1 compared to other regions. (c) Y concentration map. (d) U concentration map showing higher content in the C1 cement.

Download

4.2 LA-ICP-MS trace elements and distribution

Trace elements were quantified in the four distinct phases through segmentation (i.e. selection of region of interests) of a LA-ICP-MS map. Additionally, single-spot analyses were performed on the C1 botryoidal cement. A comprehensive presentation of all data is provided in Table S5. Selected elements are highlighted for comparative analysis among the phases in Fig. 1d.

Notably, the C1 botryoidal cement suggested as the potential new RM exhibits significantly lower concentrations of Mg, Mn, and Fe compared to all other phases and higher levels of U and Sr compared to C2 and RFSC3 phases. The C1 phase is identified best by this set of elements. In the matrix, U and Sr concentrations are comparable, but all elements show a more heterogeneous distribution. Uranium concentration in C1 is remarkably homogeneous, falling within the range of 4 µg g−1 with a precision as relative standard deviation (RSD) of less than 10 %.

Conversely, the replacement cement C2 displays elevated concentrations of Mg, Mn, and Fe. Owing to its lower U content and higher initial Pb, spot analyses within this phase plot closer to the upper intercept in the Tera–Wasserburg (TW) diagram (Fig. 1c). This phase exhibits a broader scatter in U / Pb ratios than cement C1.

Additional high-resolution TE maps (Fig. 2) reveal the fibrous structure of C1 cement in the Y and U content and show some isolated, high-Mg spots in C1 cement. The higher spatial resolution imaging also captures details like the decreasing Y content towards the rim of the C1 cement, which are not detected on our low-resolution imaging (Fig. S6).

https://gchron.copernicus.org/articles/6/465/2024/gchron-6-465-2024-f03

Figure 3LA-ICP-MS U–Pb results for RA138 using WC-1 as RM. (a) Pooled isochron from 23 sessions. (b) Typical results from one analytical session. (c) Rank order plot of isochron ages and weighted mean calculation for 23 sessions showing excess scatter with a MSWD of 7.6. (d) Same data as panel (c) but including a systematic uncertainty of 1.5 %, resulting in a MSWD of 1.3.

Download

4.3 LA-ICP-MS U–Pb dating

Cement C1 was dated during 23 sessions (November 2022–July 2023), using WC-1 as the primary reference material for 238U/206Pb ratio correction. A detailed dataset is provided in Table S2a. Approximately 5 % of the RA138 analyses were defined as outliers based on their deviation from the isochron. This divergence is likely attributed to the ablation of different phases other than C1 or mixing of C1 with other phases (i.e. misplacing ablation spots or ablating other phases in lower parts of ablation pits), notably observed as bright spots within C1 (Fig. 1b). When analysing all different phases (C1, C2, RSFC3, and the matrix; Fig. 1c), the number of points deviating from the isochron in the TW space increases, some uncertainty ellipses are bigger due to low U concentration, the age becomes younger, and the initial Pb composition can be different.

The pooled TW isochron, comprising n=763 spot analyses (Fig. 3a and data in Table S2b), reveals a lower-intercept age of 319.25 ± 0.48 Ma (without propagation of excess uncertainty Ssys), with a mean square weighted deviation (MSWD) of 0.81, and an initial 207Pb/206Pb ratio of 0.8486 ± 0.0016 (unanchored). Figure 3b shows a representative session with n=30, producing a lower-intercept age of 322.69 ± 2.09 Ma (without propagation of excess uncertainty Ssys), a MSWD of 0.37, and an initial 207Pb/206Pb ratio of 0.8495 ± 0.0051 anchored to the ID-TIMS value of 0.8495 ± 0.0065.

The 23 lower-intercept ages span from 314.26 ± 2.25 to 323.83 ± 3.23 Ma, with MSWD values ranging from 0.12 to 2.0 (see Table S2a for details). The calculated weighted mean of the lower-intercept ages of all 23 sessions yields an age of 319.25 ± 1.33 Ma (0.4 %; without propagation of excess uncertainty Ssys) with a MSWD of 7.6 (Fig. 3c and data and calculations in Table S3). To address potential systematic uncertainties, a quadratic propagation of a 1.5 % uncertainty to each of the 23 individual intercept ages results in a comparable weighted mean age of 319.32 ± 1.1 (0.34 %; without propagation of excess uncertainty Ssys), with a more acceptable MSWD of 1.3 (Fig. 3d). All computations were conducted using IsoplotR and model 1 discordia with 95 % confidence interval uncertainties.

4.4 ID-TIMS U–Pb dating

We present a set of 21 ID-TIMS analyses conducted on botryoidal cement C1 during two distinct measurement campaigns. The initial campaign (n=9) took place at the University of Geneva, while the subsequent campaign involved a collaborative effort between the laboratories of ETH Zürich (n=8) and Geneva (n=4). Detailed ID-TIMS data are provided in Table S4, and the complete dataset of 21 analyses is shown in Fig. 4. The resulting isochron of all analyses exhibits a lower-intercept age of 320.76 ± 1.31 Ma, an initial 207Pb/206Pb ratio of 0.8378 ± 0.0102, and a significant overdispersion reflected in a MSWD of 76.

https://gchron.copernicus.org/articles/6/465/2024/gchron-6-465-2024-f04

Figure 4Tera–Wasserburg diagram of ID-TIMS results for RA138 cement C1. Reference isochron values are shown in blue (“preferred”), values for all data including four excluded points with highest 207Pb/206Pb ratios (“all”) are shown in black. The insert shows the CL image before ID-TIMS sampling and includes sampling locations.

Download

Notably, the four analyses with the highest 207Pb/206Pb ratio are distant from the isochron constructed with the remaining data, indicating potential mixing with non-C1-phase material which may be characterized by either a distinct common Pb composition or age (Fig. 1c). Examination of the sampling locations during the second campaign, as illustrated in the CL image (Fig. 4), suggests the possibility of mixtures of the micro-drilled sample due to the depth of pit required to obtain sufficient material for the ID-TIMS analyses. Excluding these four points from the age interpretation (Fig. 4) yields the preferred lower-intercept age of 321.99 ± 0.65, an initial 207Pb/206Pb ratio of 0.8495 ± 0.0065, and a MSWD of 5.1. The elevated MSWD value indicates some remaining sample heterogeneity; due to the large sampling volume, the mixed sampling of phases other than C1 seems possible. These refined results, including the overdispersion term (Vermeesch, 2018) for both intercepts, are considered the preferred values for utilization as a reference in LA-ICP-MS work.

5 Discussion

Botryoidal cement C1 within sample RA138 displays a favourable level of age homogeneity and a large spread of U/Pb ratios, making it a suitable it a suitable primary RM for the 238U/206Pb correction in carbonate LA-ICP-MS U–Pb dating. The 0.2 % age uncertainty from the ID-TIMS analyses enhances the precision of the correction factor and consequently improves the age uncertainty in the unknown samples. This stands in contrast to WC-1 which has a higher uncertainty of 2.5 % (254.4 ± 6.4 Ma) (Roberts et al., 2017).

Furthermore, a re-evaluation of the excess uncertainty (Ssys) previously estimated for U–Pb analysis by LA-ICP-MS, particularly in zircons (Horstwood et al., 2016), indicates approximately 1.5 % excess uncertainty (Ssys) based on repeated analyses of RA138 (Fig. 3c and d), challenging the previously assumed range of 2 %–2.5 % (Guillong et al., 2020). RA138 thus provides an improved reference material with less scatter around the U/Pb isochron as a prerequisite for improving this method.

While RA138 contributes to enhanced precision, it is essential to acknowledge its lack of homogeneity. Apart from C1, other phases are present. Precise selection of analysis locations for calculating correction factors is crucial, avoiding less suitable locations such as the matrix or replacement cements C2 and RFSC3. C2, for instance, exhibits higher initial Pb than C1, with some analyses aligning with C1 but a substantial portion revealing open-system behaviour and younger crystallization ages (Fig. 1c). Radiaxial fibrous and sparry cement, RFSC3, containing less U, demonstrates an overall younger age with much higher variability (Fig. 1c).

The matrix, although visually heterogeneous and displaying a broad range of U–Pb ratios, is not recommended as RM due to potential open-system behaviour and some younger ages. Even when analysing only C1, up to 5 % of point analyses may be outliers, likely arising from the partial analysis of small parts of C2, identifiable as small bright spots in the CL image (Fig. 1b). Therefore, it is suggested to conduct approximately 30 analyses per session for a robust correction factor calculation. It is also recommended to image well-characterized RMs sections used in each laboratory (e.g. petrography and trace element screening) for better visualization and assessment of the different phases.

Finally, while cement C1 of RA138 exhibits a larger variation in U/Pb than other commonly used RMs, and unanchored data typically have an upper intercept within the uncertainty in the ID-TIMS initial 207Pb/206Pb ratio, we recommend anchoring the isochron to the ID-TIMS-determined initial 207Pb/206Pb ratio of 0.8495 ± 0.0065 for the purpose of accurate correction factor calculation.

6 Conclusions, availability, and outlook

Botryoidal cement C1 from RA138 provides a significant potential for reducing uncertainties related to primary reference material corrections during U–Pb LA-ICP-MS carbonate dating. Its variable U/Pb ratios, coupled with a relatively high and homogenous U content, contribute to a reasonably low uncertainty, resulting in a well-defined isochron, particularly when anchored to the initial Pb ratio. This material is sourced from an accessible outcrop in northern Spain, ensuring unlimited supply, and numerous aliquots have been prepared for distribution. Detailed information on aliquots is included in Figs. S1–S5, and interested parties can request additional details from the corresponding author. With the dissemination of this new RM, it is expected that achievable precision in LA-ICP-MS calcite dating will improve substantially. This is a crucial step forward, while the community is searching for RMs that would be homogeneous not only in age but also in U/Pb ratio – to be used directly for all corrections – making data reduction and downhole fractionation correction significantly easier.

Data availability

All data produced in this study are available in the Supplement.

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/gchron-6-465-2024-supplement.

Author contributions

MG, IM, ES, and MO planned the campaign. OMT and JRB collected the samples. ES, LT, NL, and MG prepared the samples. ES did the petrography. MG, LT, and NL did LA-ICP-MS U–Pb analyses. MG, LT, DS, and YB did LA-ICP-MS mapping analyses. IM and DS did ID-TIMS analyses. MG, LT, and NL did the CL imaging. MG wrote the draft. All authors reviewed and edited the paper.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

Acknowledgements

The authors acknowledge the critical reading and comments by Perach Nuriel and reviews by Sota Niki and Greg Ludvigson that helped improve the paper.

Financial support

This work has received funding from the Swiss National Science Foundation (grant nos. 200021_160019, 206021_133771, 200021_169849, and 200021_182556), the Spanish Government (grant nos. CGL2013-44458-P, PGC2018-099698-B25 I00, and PRX19/00423), the Principality of Asturias (grant no. AYUD/2021/51293), and European Regional Development Fund (FEDER). Dawid Szymanowski and Lorenzo Tavazzani have been partly supported by an ETH Zürich Career Seed Award.

Review statement

This paper was edited by Klaus Mezger and reviewed by Sota Niki and Greg Ludvigson.

References

Bowring, J. F., McLean, N. M., and Bowring, S. A.: Engineering cyber infrastructure for U-Pb geochronology: Tripoli and U-Pb_Redux, Geochem. Geophy. Geosy., 12, https://doi.org/10.1029/2010gc003479, 2011. 

Brigaud, B., Bonifacie, M., Pagel, M., Blaise, T., Calmels, D., Haurine, F., and Landrein, P.: Past hot fluid flows in limestones detected by Δ47-(U-Pb) and not recorded by other geothermometers, Geology, 48, 851–856, https://doi.org/10.1130/g47358.1, 2020. 

Burisch, M., Gerdes, A., Walter, B. F., Neumann, U., Fettel, M., and Markl, G.: Methane and the origin of five-element veins: Mineralogy, age, fluid inclusion chemistry and ore forming processes in the Odenwald, SW Germany, Ore Geol. Rev., 81, 42–61, https://doi.org/10.1016/j.oregeorev.2016.10.033, 2017. 

Chaldekas, O., Vaks, A., Haviv, I., Gerdes, A., and Albert, R.: U-Pb speleothem geochronology reveals a major 6 Ma uplift phase along the western margin of Dead Sea Transform, Geol. Soc. Am. Bull., 134, 1571–1584, https://doi.org/10.1130/b36051.1, 2022. 

Chen, X., Tissot, F. L. H., Jansen, M. F., Bekker, A., Liu, C. X., Nie, N. X., Halverson, G. P., Veizer, J., and Dauphas, N.: The uranium isotopic record of shales and carbonates through geologic time, Geochim. Cosmochim. Ac., 300, 164–191, https://doi.org/10.1016/j.gca.2021.01.040, 2021. 

Chew, D. M., Petrus, J. A., and Kamber, B. S.: U-Pb LA-ICPMS dating using accessory mineral standards with variable common Pb, Chem. Geol., 363, 185–199, https://doi.org/10.1016/j.chemgeo.2013.11.006, 2014. 

Condon, D. J., Schoene, B., McLean, N. M., Bowring, S. A., and Parrish, R. R.: Metrology and traceability of U-Pb isotope dilution geochronology (EARTHTIME Tracer Calibration Part I), Geochim. Cosmochim. Ac., 164, 464–480, https://doi.org/10.1016/j.gca.2015.05.026, 2015. 

Coogan, L. A., Parrish, R. R., and Roberts, N. M. W.: Early hydrothermal carbon uptake by the upper oceanic crust: Insight from in situ U-Pb dating, Geology, 44, 147–150, https://doi.org/10.1130/g37212.1, 2016. 

Drake, H., Heim, C., Roberts, N. M. W., Zack, T., Tillberg, M., Broman, C., Ivarsson, M., Whitehouse, M. J., and Astrom, M. E.: Isotopic evidence for microbial production and consumption of methane in the upper continental crust throughout the Phanerozoic eon, Earth Planet. Sc. Lett., 470, 108–118, https://doi.org/10.1016/j.epsl.2017.04.034, 2017. 

Elisha, B., Nuriel, P., Kylander-Clark, A., and Weinberger, R.: Towards in situ U–Pb dating of dolomite, Geochronology, 3, 337–349, https://doi.org/10.5194/gchron-3-337-2021, 2021. 

Gerstenberger, H. and Haase, G.: A highly effective emitter substance for mass spectrometric Pb isotope ratio determinations, Chem. Geol., 136, 309–312, https://doi.org/10.1016/s0009-2541(96)00033-2, 1997. 

Giorno, M., Barale, L., Bertok, C., Frenzel, M., Looser, N., Guillong, M., Bernasconi, S. M., and Martire, L.: Sulfide-associated hydrothermal dolomite and calcite reveal a shallow burial depth for Alpine-type Zn-(Pb) deposits, Geology, 50, 853–858, https://doi.org/10.1130/g49812.1, 2022. 

Godeau, N., Deschamps, P., Guihou, A., Leonide, P., Tendil, A., Gerdes, A., Hamelin, B., and Girard, J. P.: U-Pb dating of calcite cement and diagenetic history in microporous carbonate reservoirs: Case of the Urgonian Limestone, France, Geology, 46, 247–250, https://doi.org/10.1130/g39905.1, 2018. 

Guillong, M., Wotzlaw, J.-F., Looser, N., and Laurent, O.: Evaluating the reliability of U–Pb laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) carbonate geochronology: matrix issues and a potential calcite validation reference material, Geochronology, 2, 155–167, https://doi.org/10.5194/gchron-2-155-2020, 2020. 

Gulbranson, E. L., Rasbury, E. T., Ludvigson, G. A., Möller, A., Henkes, G. A., Suarez, M. B., Northrup, P., Tappero, R. V., Maxson, J. A., Shapiro, R. S., and Wooton, K. M.: U-Pb Geochronology and Stable Isotope Geochemistry of Terrestrial Carbonates, Lower Cretaceous Cedar Mountain Formation, Utah: Implications for Synchronicity of Terrestrial and Marine Carbon Isotope Excursions, Geosciences, 12, 346, https://doi.org/10.3390/geosciences12090346, 2022. 

Hansman, R. J., Albert, R., Gerdes, A., and Ring, U.: Absolute ages of multiple generations of brittle structures by U-Pb dating of calcite, Geology, 46, 207–210, https://doi.org/10.1130/g39822.1, 2018. 

Hiess, J., Condon, D. J., McLean, N., and Noble, S. R.: 238U/235U Systematics in terrestrial uranium-bearing minerals, Science, 335, 1610–1614, https://doi.org/10.1126/science.1215507, 2012. 

Hoareau, G., Claverie, F., Pecheyran, C., Paroissin, C., Grignard, P.-A., Motte, G., Chailan, O., and Girard, J.-P.: Direct U–Pb dating of carbonates from micron-scale femtosecond laser ablation inductively coupled plasma mass spectrometry images using robust regression, Geochronology, 3, 67–87, https://doi.org/10.5194/gchron-3-67-2021, 2021. 

Horstwood, M. S., Košler, J., Gehrels, G., Jackson, S. E., McLean, N. M., Paton, C., Pearson, N. J., Sircombe, K., Sylvester, P., and Vermeesch, P.: Community-Derived Standards for LA-ICP-MS U-(Th-)Pb Geochronology–Uncertainty Propagation, Age Interpretation and Data Reporting, Geostand. Geoanal. Res., 40, 311–332, https://doi.org/10.1111/j.1751-908X.2016.00379.x, 2016. 

Jaffey, A. H., Flynn, K. F., Glendenin, L. E., Bentley, W. C., and Essling, A. M.: Precision Measurement of Half-Lives and Specific Activities of U-235 and U-238, Phys. Rev. C, 4, 1889, https://doi.org/10.1103/PhysRevC.4.1889, 1971. 

Jochum, K. P., Nohl, L., Herwig, K., Lammel, E., Stoll, B., and Hofmann, A. W.: GeoReM: A new geochemical database for reference materials and isotopic standards, Geostand. Geoanalyt. Res., 29, 333–338, https://doi.org/10.1111/j.1751-908X.2005.tb00904.x, 2005. 

Kylander-Clark, A. R. C.: Expanding the limits of laser-ablation U–Pb calcite geochronology, Geochronology, 2, 343–354, https://doi.org/10.5194/gchron-2-343-2020, 2020. 

Li, Q., Parrish, R. R., Horstwood, M. S. A., and McArthur, J. M.: U-Pb dating of cements in Mesozoic ammonites, Chem. Geol., 376, 76–83, https://doi.org/10.1016/j.chemgeo.2014.03.020, 2014. 

Looser, N., Madritsch, H., Guillong, M., Laurent, O., Wohlwend, S., and Bernasconi, S. M.: Absolute Age and Temperature Constraints on Deformation Along the Basal Decollement of the Jura Fold-and-Thrust Belt From Carbonate U-Pb Dating and Clumped Isotopes, Tectonics, 40, e2020TC006439, https://doi.org/10.1029/2020TC006439, 2021. 

MacDonald, J. M., Faithfull, J. W., Roberts, N. M. W., Davies, A. J., Holdsworth, C. M., Newton, M., Williamson, S., Boyce, A., and John, C. M.: Clumped-isotope palaeothermometry and LA-ICP-MS U-Pb dating of lava-pile hydrothermal calcite veins, Contrib. Miner. Petrol., 174, 63, https://doi.org/10.1007/s00410-019-1599-x, 2019. 

Mangenot, X., Gasparrini, M., Gerdes, A., Bonifacie, M., and Rouchon, V.: An emerging thermochronometer for carbonate-bearing rocks: Δ47/(U-Pb), Geology, 46, 1067–1070, https://doi.org/10.1130/g45196.1, 2018. 

McLean, N. M., Condon, D. J., Schoene, B., and Bowring, S. A.: Evaluating uncertainties in the calibration of isotopic reference materials and multi-element isotopic tracers (EARTHTIME Tracer Calibration Part II), Geochim. Cosmochim. Ac., 164, 481–501, https://doi.org/10.1016/j.gca.2015.02.040, 2015. 

Methner, K., Mulch, A., Fiebig, J., Wacker, U., Gerdes, A., Graham, S. A., and Chamberlain, C. P.: Rapid Middle Eocene temperature change in western North America, Earth Planet. Sc. Lett., 450, 132–139, https://doi.org/10.1016/j.epsl.2016.05.053, 2016. 

Mottram, C. M., Kellett, D. A., Barresi, T., Zwingmann, H., Friend, M., Todd, A., and Percival, J. B.: Syncing fault rock clocks: Direct comparison of U-Pb carbonate and K-Ar illite fault dating methods, Geology, 48, 1179–1183, https://doi.org/10.1130/g47778.1, 2020. 

Norris, C. A., Danyushevsky, L., Olin, P., and West, N. R.: Elimination of aliasing in LA-ICP-MS by alignment of laser and mass spectrometer, J. Anal. Atom. Spectrom., 36, 733–739, https://doi.org/10.1039/d0ja00488j, 2021. 

Nuriel, P., Weinberger, R., Kylander-Clark, A. R. C., Hacker, B. R., and Craddock, J. P.: The onset of the Dead Sea transform based on calcite age-strain analyses, Geology, 45, 587–590, https://doi.org/10.1130/g38903.1, 2017. 

Nuriel, P., Craddock, J., Kylander-Clark, A. R. C., Uysal, I. T. S., Karabacak, V., Dirik, R. K., Hacker, B. R., and Weinberger, R.: Reactivation history of the North Anatolian fault zone based on calcite age-strain analyses, Geology, 47, 465–469, https://doi.org/10.1130/g45727.1, 2019. 

Nuriel, P., Wotzlaw, J.-F., Ovtcharova, M., Vaks, A., Stremtan, C., Šala, M., Roberts, N. M. W., and Kylander-Clark, A. R. C.: The use of ASH-15 flowstone as a matrix-matched reference material for laser-ablation U − Pb geochronology of calcite, Geochronology, 3, 35–47, https://doi.org/10.5194/gchron-3-35-2021, 2021. 

Pagel, M., Bonifacie, M., Schneider, D. A., Gautheron, C., Brigaud, B., Calmels, D., Cros, A., Saint-Bezar, B., Landrein, P., Sutcliffe, C., Davis, D., and Chaduteau, C.: Improving paleohydrological and diagenetic reconstructions in calcite veins and breccia of a sedimentary basin by combining Δ47 temperature, δ18Owater and U-Pb age, Chem. Geol., 481, 1–17, https://doi.org/10.1016/j.chemgeo.2017.12.026, 2018. 

Paton, C., Hellstrom, J., Paul, B., Woodhead, J., and Hergt, J.: Iolite: Freeware for the visualisation and processing of mass spectrometric data, J. Anal. Atom. Spectrom., 26, 2508–2518, https://doi.org/10.1039/c1ja10172b, 2011. 

Paton, C., Woodhead, J. D., Hellstrom, J. C., Hergt, J. M., Greig, A., and Maas, R.: Improved laser ablation U-Pb zircon geochronology through robust downhole fractionation correction, Geochem. Geophy. Geosy., 11, https://doi.org/10.1029/2009gc002618, 2010. 

Petrus, J. A. and Kamber, B. S.: VizualAge: A Novel Approach to Laser Ablation ICP-MS U-Pb Geochronology Data Reduction, Geostand. Geoanal. Res., 36, 247–270, https://doi.org/10.1111/j.1751-908X.2012.00158.x, 2012. 

Piccione, G., Rasbury, E. T., Elliott, B. A., Kyle, J. R., Jaret, S. J., Acerbo, A. S., Lanzirotti, A., Northrup, P., Wooton, K., and Parrish, R. R.: Vein fluorite U-Pb dating demonstrates post 6.2 Ma rare-earth element mobilization associated with Rio Grande rifting, Geosphere, 15, 1958–1972, https://doi.org/10.1130/ges02139.1, 2019. 

Rasbury, E. T., Piccione, G., Holt, W., and Ward, W. B.: Potential for constraining sequence stratigraphy and cycle stratigraphy with U-Pb dating of carbonates, Earth-Sci. Rev., 243, 104495, https://doi.org/10.1016/j.earscirev.2023.104495, 2023. 

Ring, U. and Gerdes, A.: Kinematics of the Alpenrhein-Bodensee graben system in the Central Alps: Oligocene/Miocene transtension due to formation of the Western Alps arc, Tectonics, 35, 1367–1391, https://doi.org/10.1002/2015tc004085, 2016. 

Roberts, N. M. W., Rasbury, E. T., Parrish, R. R., Smith, C. J., Horstwood, M. S. A., and Condon, D. J.: A calcite reference material for LA-ICP-MS U-Pb geochronology, Geochem. Geophy. Geosy., 18, 2807–2814, https://doi.org/10.1002/2016gc006784, 2017. 

Scardia, G., Parenti, F., Miggins, D. P., Gerdes, A., Araujo, A. G. M., and Neves, W. A.: Chronologic constraints on hominin dispersal outside Africa since 2.48 Ma from the Zarqa Valley, Jordan, Quaternary Sci. Rev., 219, 1–19, https://doi.org/10.1016/j.quascirev.2019.06.007, 2019. 

Schmitz, M. D. and Schoene, B.: Derivation of isotope ratios, errors, and error correlations for U-Pb geochronology using 205Pb-235U-(233U)-spiked isotope dilution thermal ionization mass spectrometric data, Geochem. Geophy. Geosy., 8, 2006GC001492, https://doi.org/10.1029/2006gc001492, 2007. 

Sindern, S., Havenith, V., Gerdes, A., Meyer, F. M., Adelmann, D., and Hellmann, A.: Dating of anatase-forming diagenetic reactions in Rotliegend sandstones of the North German Basin, Int. J. Earth Sci., 108, 1275–1292, https://doi.org/10.1007/s00531-019-01705-x, 2019. 

Tavazzani, L., Guillong, M., Giuliani, A., Fontboté, L., and Chelle-Michou, C.: Not so fast: Million-years of metal precipitation in Mississippi Valley type deposits inferred from in-situ petrochronology of hydrothermal carbonates, Earth Planet. Sc. Lett., 636, 118718, https://doi.org/10.1016/j.epsl.2024.118718, 2024.  

Vermeesch, P.: IsoplotR: A free and open toolbox for geochronology, Geosci. Front., 9, 1479–1493, https://doi.org/10.1016/j.gsf.2018.04.001, 2018. 

von Quadt, A., Wotzlaw, J. F., Buret, Y., Large, S. J. E., Peytcheva, I., and Trinquier, A.: High-precision zircon U/Pb geochronology by ID-TIMS using new 1013 ohm resistors, J. Anal. Atom. Spectrom., 31, 658–665, https://doi.org/10.1039/c5ja00457h, 2016. 

Weinberger, R., Nuriel, P., Kylander-Clark, A. R. C., and Craddock, J. P.: Temporal and spatial relations between large-scale fault systems: Evidence from the Sinai-Negev shear zone and the Dead Sea Fault, Earth-Sci. Rev., 211, 103377, https://doi.org/10.1016/j.earscirev.2020.103377, 2020. 

Woodhead, J. and Petrus, J.: Exploring the advantages and limitations of in situ U–Pb carbonate geochronology using speleothems, Geochronology, 1, 69–84, https://doi.org/10.5194/gchron-1-69-2019, 2019. 

Wotzlaw, J. F., Buret, Y., Large, S. J. E., Szymanowski, D., and von Quadt, A.: ID-TIMS U-Pb geochronology at the 0.1 ‰ level using 1013Ω resistors and simultaneous U and 18O/16O isotope ratio determination for accurate UO2 interference correction, J. Anal. At. Spectrom., 32, 579–586, https://doi.org/10.1039/c6ja00278a, 2017. 

Wu, S. T., Yang, Y. H., Roberts, N. M. W., Yang, M., Wang, H., Lan, Z. W., Xie, B. H., Li, T. Y., Xu, L., Huang, C., Xie, L. W., Yang, J. H., and Wu, F. Y.: In situ calcite U-Pb geochronology by high-sensitivity single-collector LA-SF-ICP-MS, Sci. China-Earth Sci., 65, 1146–1160, https://doi.org/10.1007/s11430-021-9907-1, 2022. 

Download
Short summary
RA138 is a new reference material for U–Pb dating of carbonate samples via laser ablation inductively coupled plasma mass spectrometry. RA138 exhibits variable U–Pb ratios and consistent U content, resulting in a precise isochron with low uncertainty. Isotope dilution thermal ionization mass spectrometry analyses fix a reference age of 321.99 ± 0.65 Ma. This research advances our ability to date carbonate samples accurately, providing insights into geological processes and historical timelines.