An empirical test of maximum depositional age through paired LA-ICP-MS and CA-ID-TIMS detrital zircon analysis with implications for the “Nevadan orogeny” in the Sierra Nevada foothills, California, U.S.A.
The interpretation of maximum depositional age (MDA) from U-Pb detrital zircon data acquired via laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) or secondary ion mass spectrometry (SIMS) is now routine; however, to date only a few studies have presented tests of such MDAs by subsequently analyzing a subset of the same detrital zircon grains with the more accurate and precise chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) method. We first generated LA-ICP-MS dates and MDAs from three turbidite sandstone samples of the Late Jurassic Mariposa Formation (Sierra Nevada foothills, California) containing large proportions of young detrital zircon grains (at the time of sediment deposition). We then removed 4–5 of the apparent youngest grains from the epoxy mount and analyzed them with CA-ID-TIMS. Despite a lack of low Th/U or high U (ppm) characteristics that might indicate U-Pb system disturbance, all LA-ICP-MS dates are younger than corresponding CA-ID-TIMS dates, on the same grain, by up to 11.8 m.y. when considered as point estimates. Only two out of 13 paired analyses overlap within error at 95% confidence. We interpret MDAs from our CA-ID-TIMS dates to be the youngest CA-ID-TIMS dates (YTDs). LA-ICP-MS MDAs based on fewer dates from the young tail of the youngest date distribution perform the worst (i.e., furthest from CA-ID-TIMS MDAs), whereas LA-ICP-MS MDA methods that incorporate more dates from the youngest date distribution perform better; the best performing method in all cases is the maximum likelihood algorithm–minimum. The performance of tested LA-ICP-MS MDA methods is improved by removing statistical outlier dates and by removing visually young “outlier” dates that drift away from the primary distribution of young dates. Our paired LA-ICP-MS/CA-ID-TIMS MDA workflow shows that the accuracy of MDA can be significantly improved by conducting CA-ID-TIMS on as few as four grains from a sample. The combination of our new CA-ID-TIMS MDAs with published CA-ID-TIMS analyses and petrochronology on the nearby Guadalupe igneous complex provides sub-million-year resolution of contemporaneous igneous and sedimentary systems during deformation and clarifies the timing of regional deformation that defines the local “Nevadan orogeny.” Rocks in the upper plate of the Bear Mountains fault zone have a CA-ID-TIMS YTD MDA of 151.71 ± 0.23 Ma, and rocks of the lower plate have a CA-ID-TIMS YTD MDA of 149.92 ± 0.11 Ma. These new MDAs, in conjunction with the observation of fabric-bearing Mariposa Formation xenoliths in the ca. 149.65 ± 0.10 Ma Guadalupe igneous complex, suggest that rocks in the upper plate of the Bear Mountains fault zone represent a slightly older (~2 m.y.) section of Mariposa Formation that was deformed and intruded prior to being juxtaposed against, and further deformed with, a slightly younger lower-plate section of Mariposa Formation in actively deforming, fault-bounded basins. Our observations are not consistent with traditional models that require that sedimentation of the Mariposa Formation ended by ca. 155 Ma. Instead, we interpret our data to be consistent with other evidence for a continuum of deformation in Late Jurassic to Early Cretaceous time and document that the regional “slatey cleavage” observed in the greater Mariposa Formation and used to define the “Nevadan orogeny” in our study area is largely younger than 149.92 ± 0.11 Ma.