skip to main content


This content will become publicly available on January 1, 2025

Title: Modeling apparent Pb loss in zircon U–Pb geochronology

Abstract. The loss of radiogenic Pb from zircon is known to be a major factor that can cause inaccuracy in the U–Pb geochronological system; hence, there is a need to better characterize the distribution of Pb loss in natural samples. Treatment of zircon by chemical abrasion (CA) has become standard practice in isotope dilution–thermal ionization mass spectrometry (ID-TIMS), but CA is much less commonly employed prior to in situ analysis via laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) or secondary ionization mass spectrometry (SIMS). Differentiating the effects of low levels of Pb loss in Phanerozoic zircon with relatively low-precision in situ U–Pb dates, where the degree of Pb loss is insufficient to cause discernible discordance, is challenging. We show that U–Pb isotopic ratios that have been perturbed by Pb loss may be modeled by convolving a Gaussian distribution that represents random variations from the true isotopic value stemming from analytical uncertainty with a distribution that characterizes Pb loss. We apply this mathematical framework to model the distribution of apparent Pb loss in 10 igneous samples that have both non-CA LA-ICP-MS or SIMS U–Pb dates and an estimate of the crystallization age, either through CA U–Pb or 40Ar/39Ar geochronology. All but one sample showed negative age offsets that were unlikely to have been drawn from an unperturbed U–Pb date distribution. Modeling apparent Pb loss using the logit–normal distribution produced good fits with all 10 samples and showed two contrasting patterns in apparent Pb loss; samples where most zircon U–Pb dates undergo a bulk shift and samples where most zircon U–Pb dates exhibited a low age offset but fewer dates had more significant offset. Our modeling framework allows comparison of relative degrees of apparent Pb loss between samples of different age, with the first and second Wasserstein distances providing useful estimates of the total magnitude of apparent Pb loss. Given that the large majority of in situ U–Pb dates are acquired without the CA treatment, this study highlights a pressing need for improved characterization of apparent Pb-loss distributions in natural samples to aid in interpreting non-CA in situ U–Pb data and to guide future data collection strategies.

 
more » « less
Award ID(s):
2243685
NSF-PAR ID:
10493102
Author(s) / Creator(s):
;
Publisher / Repository:
Geochronology
Date Published:
Journal Name:
Geochronology
Volume:
6
Issue:
1
ISSN:
2628-3719
Page Range / eLocation ID:
37 to 51
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Ascough, P. ; Dunai, T. ; King, G. ; Lang, A. ; Mezger, K. (Ed.)
    Detrital zircon geochronology by laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) is a widely used tool for determining maximum depositional ages and sediment provenance, as well as reconstructing sediment routing pathways. Although the accuracy and precision of U–Pb geochronology measurements have improved over the past 2 decades, Pb loss continues to impact the ability to resolve zircon age populations by biasing affected zircon toward younger apparent ages. Chemical abrasion (CA) has been shown to reduce or eliminate the effects of Pb loss in zircon U–Pb geochronology but has yet to be widely applied to large-n detrital zircon analyses. Here, we assess the efficacy of the chemical abrasion treatment on zircon prior to analysis by LA-ICP-MS and discuss the advantages and limitations of this technique in relation to detrital zircon geochronology. We show that (i) CA does not systematically bias LA-ICP-MS U–Pb dates for 13 reference materials that span a wide variety of crystallization dates and U concentrations, (ii) CA-LA-ICP-MS U–Pb zircon geochronology can reduce or eliminate Pb loss in samples that have experienced significant radiation damage, and (iii) bulk CA prior to detrital zircon U–Pb geochronology by LA-ICP-MS improves the resolution of age populations defined by 206Pb/238U dates (Neoproterozoic and younger) and increases the percentage of concordant analyses in age populations defined by 207Pb/206Pb dates (Mesoproterozoic and older). The selective dissolution of zircon that has experienced high degrees of radiation damage suggests that some detrital zircon age populations could be destroyed or have their abundance significantly modified during this process. However, we did not identify this effect in either of the detrital zircon samples that were analyzed as part of this study. We conclude that pre-treatment of detrital zircon by bulk CA may be useful for applications that require increased resolution of detrital zircon populations and increased confidence that 206Pb/238U dates are unaffected by Pb loss. 
    more » « less
  2. Abstract

    Detrital zircon (DZ) U–Pb geochronology has improved the way geologists approach questions of sediment provenance and stratigraphic age. However, there is debate about what constitutes an appropriate sample size (i.e., the number of dates in a DZ sample,n), which depends on project objectives, sample complexity, and, critically, analytical budget. Additionally, there is ongoing concern about bias introduced by zircon grain size. We tested a recently developed rapid (3 s/analysis) data acquisition method by multicollector laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) that incorporates an automated selection routine and calculates two‐dimensional grain geometry from polished sample surfaces. Eleven samples were analysed from below and above the Late Cretaceous (Campanian) basal Castlegate unconformity of the Book Cliffs, Utah, in a down‐depositional‐dip transect including Price, Horse, Tusher, and Thompson canyons. 12,448 new concordant dates were generated during two measurement sessions. Results are consistent with recent studies suggesting there is no major provenance change and little time (1–2 Myr) represented across the unconformity. Grain size and sample size both exert a strong control on sample dissimilarity. Age distributions constructed from subsamples of large grains are systematically less similar to whole samples; age distributions composed of small grains are overall more similar to whole samples. As such, North American sediment sources that produce large grains such as the Grenville and Yavapi‐Mazatzal belts can bias age distributions if only large grains are analysed. A sample size ofn = 100 is inadequate for characterizing age distributions as complex as those of the Book Cliffs, whereas a sample size ofn = 300 provides good characterization. Sample size ofn ≈ 1000 or more is unnecessary unless project objectives include scanning for subordinate age groups, such as when identifying the youngest grains for calculating a maximum depositional age (MDA). Dates used in MDA calculations acquired with rapid acquisition are best re‐analysed with longer LA‐ICP‐MS acquisition methods or isotope dilution thermal ionization mass spectrometry for increased accuracy and precision. We include new MATLAB code and open‐source software programs,DZpickandDZmda, for automated spot picking and calculating MDAs.

     
    more » « less
  3. Detrital zircon (DZ) U‐Pb laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) has revolutionised the way geologists approach many Earth science questions. Although recent research has focused on rapid sample throughput, acquisition rates are limited to 100–300 analyses h−1. We present a method to acquire zircon U‐Pb dates at rates of 120, 300, 600 and 1200 analyses h−1(30, 12, 6 and 3 s per analysis) by multi‐collector LA‐ICP‐MS. We demonstrate the efficacy of this method by analysing twelve zircon reference materials with dates from ~ 3465 to ~ 28 Ma. Mean offset from high‐precision dates increases with faster rates from 0.9% to 1.1%; mean random 1suncertainty increases from 0.6% to 1.3%. We tested this new method on a sandstone sample previously characterised by large‐nDZ geochronology. Quantitative comparison shows increased correspondence among age distributions comprising > 300 dates. This new method holds promise for DZ geochronology because (a) it requires no major changes to hardware, but rather modifications to software; (b) it yields robust age distributions well‐suited for quantitative analysis and maximum depositional age calculations; (c) there is only a minor sacrifice of accuracy and measurement uncertainty; and (d) there is less burden to researchers in terms of time investment and analytical cost.

     
    more » « less
  4. Abstract Despite being a prominent continental-scale feature, the late Mesoproterozoic North American Midcontinent Rift did not result in the break-up of Laurentia, and subsequently underwent structural inversion. The timing of inversion is critical for constraining far-field effects of orogenesis and processes associated with the rift's failure. The Keweenaw fault in northern Michigan (USA) is a major thrust structure associated with rift inversion; it places ca. 1093 Ma rift volcanic rocks atop the post-rift Jacobsville Formation, which is folded in its footwall. Previous detrital zircon (DZ) U-Pb geochronology conducted by laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) assigned a ca. 950 Ma maximum age to the Jacobsville Formation and led researchers to interpret its deposition and deformation as postdating the ca. 1090–980 Ma Grenvillian Orogeny. In this study, we reproduced similar DZ dates using LA-ICP-MS and then dated 19 of the youngest DZ grains using high-precision chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS). The youngest DZ dated by CA-ID-TIMS at 992.51 ± 0.64 Ma (2σ) redefines the maximum depositional age of the Jacobsville Formation and overlaps with a U-Pb LA-ICP-MS date of 985.5 ± 35.8 Ma (2σ) for late-kinematic calcite veins within the brecciated Keweenaw fault zone. Collectively, these data are interpreted to constrain deposition of the Jacobsville Formation and final rift inversion to have occurred during the 1010–980 Ma Rigolet Phase of the Grenvillian Orogeny, following an earlier phase of Ottawan inversion. Far-field deformation propagated >500 km into the continental interior during the Ottawan and Rigolet phases of the Grenvillian Orogeny. 
    more » « less
  5. We report U-Pb age determinations of carbonate nodules from an in situ paleosol horizon in the Upper Permian Balfour Formation and from several horizons of pedogenic nodule conglomerate (PNC) in the Triassic Katberg Formation, Karoo Basin, South Africa, using laser ablation−inductively coupled plasma−mass spectrometry (LA-ICP-MS). The paleosol sample yields an age of 252 ± 3 Ma, which overlaps with a previous high-precision U-Pb zircon date from a volcanic ash deposit 2 m above the paleosol. This relationship demonstrates the reliability of using LA-ICP-MS dating techniques on terrestrial pedogenic calcite. Two PNC samples collected at the base of the Katberg Formation within the same sandstone unit yield ages of 255 ± 3 Ma and 251 ± 3 Ma. The age of 251 ± 3 Ma overlaps with the high-precision U-Pb zircon date below the PNC and is a maximum age estimate of deposition for the base of the Katberg Formation. Our results show that reworked nodules in the same concentrated conglomerate lag can be of different ages, but that similarly aged nodules are spatially associated. In addition, two PNC samples collected higher in the section yield ages of 249 ± 3 Ma and 241 ± 3 Ma, providing maximum depositional ages for the lower to middle Katberg Formation for the first time. We demonstrate that pedogenic carbonate nodules can be dated with meaningful precision, providing another mechanism for constraining the age of sedimentary sequences and studying events associated with the Permian−Triassic transition in the central Karoo Basin, even though the extinction boundary may not be preserved in this area. 
    more » « less