skip to main content


Title: A hidden Rodinian lithospheric keel beneath Zealandia, Earth's newly recognized continent
Abstract We present a data set of >1500 in situ O-Hf-U-Pb zircon isotope analyses that document the existence of a concealed Rodinian lithospheric keel beneath continental Zealandia. The new data reveal the presence of a distinct isotopic domain of Paleozoic–Mesozoic plutonic rocks that contain zircon characterized by anomalously low δ18O values (median = +4.1‰) and radiogenic εHf(t) (median = +6.1). The scale (>10,000 km2) and time span (>>250 m.y.) over which plutonic rocks with this anomalously low-δ18O signature were emplaced appear unique in a global context, especially for magmas generated and emplaced along a continental margin. Calculated crustal-residence ages (depleted mantle model, TDM) for this low-δ18O isotope domain range from 1300 to 500 Ma and are interpreted to represent melting of a Precambrian lithospheric keel that was formed and subsequently hydrothermally altered during Rodinian assembly and rifting. Recognition of a concealed Precambrian lithosphere beneath Zealandia and the uniqueness of the pervasive low-δ18O isotope domain link Zealandia to South China, providing a novel test of specific hypotheses of continental block arrangements within Rodinia.  more » « less
Award ID(s):
1650183 1655152
NSF-PAR ID:
10323892
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Geology
Volume:
49
Issue:
8
ISSN:
0091-7613
Page Range / eLocation ID:
1009 to 1014
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present >500 zircon δ18O and Lu-Hf isotope analyses on previously dated zircons to explore the interplay between spatial and temporal magmatic signals in Zealandia Cordillera. Our data cover ~8500 km2 of middle and lower crust in the Median Batholith (Fiordland segment of Zealandia Cordillera) where Mesozoic arc magmatism along the paleo-Pacific margin of Gondwana was focused along an ~100 km wide, arc-parallel zone. Our data reveal three spatially distinct isotope domains that we term the eastern, central, and western isotope domains. These domains parallel the Mesozoic arc-axis, and their boundaries are defined by major crustal-scale faults that were reactivated as ductile shear zones during the Early Cretaceous. The western isotope domain has homogenous, mantle-like δ 18O (Zrn) values of 5.8 ± 0.3‰ (2 St.dev.) and initial εHf (Zrn) values of +4.2 ± 1.0 (2 St.dev.). The eastern isotope domain is defined by isotopically low and homogenous δ18O (Zrn) values of 3.9 ± 0.2‰ and initial εHf values of +7.8 ± 0.6. The central isotope domain is characterized by transitional isotope values that display a strong E-W gradient with δ18O (Zrn) values rising from 4.6 to 5.9‰ and initial εHf values decreasing from +5.5 to +3.7. We find that the isotope architecture of the Median Batholith was in place before the initiation of Mesozoic arc magmatism and pre-dates Early Cretaceous contractional deformation and transpression. Our data show that Mesozoic pluton chemistry was controlled in part by long-lived, spatially distinct isotope domains that extend from the crust through to the upper mantle. Isotope differences between these domains are the result of the crustal architecture (an underthrusted low-δ18O source terrane) and a transient event beginning at ca. 129 Ma that primarily involved a depleted-mantle component contaminated by recycled trench sediments (10–20%). When data showing the temporal and spatial patterns of magmatism are integrated, we observe a pattern of decreasing crustal recycling of the low-δ18O source over time, which ultimately culminated in a mantle-controlled flare-up. Our data demonstrate that spatial and temporal signals are intimately linked, and when evaluated together they provide important insights into the crustal architecture and the role of both stable and transient arc magmatic trends in Cordilleran batholiths. 
    more » « less
  2. null (Ed.)
    Zircons widely occur in magmatic rocks and often display internal zonation finely recording the magmatic history. Here, we presented in situ high-precision (2SD <0.15‰ for δ 94 Zr) and high–spatial-resolution (20 µm) stable Zr isotope compositions of magmatic zircons in a suite of calc-alkaline plutonic rocks from the juvenile part of the Gangdese arc, southern Tibet. These zircon grains are internally zoned with Zr isotopically light cores and increasingly heavier rims. Our data suggest the preferential incorporation of lighter Zr isotopes in zircon from the melt, which would drive the residual melt to heavier values. The Rayleigh distillation model can well explain the observed internal zoning in single zircon grains, and the best-fit models gave average zircon–melt fractionation factors for each sample ranging from 0.99955 to 0.99988. The average fractionation factors are positively correlated with the median Ti-in-zircon temperatures, indicating a strong temperature dependence of Zr isotopic fractionation. The results demonstrate that in situ Zr isotope analyses would be another powerful contribution to the geochemical toolbox related to zircon. The findings of this study solve the fundamental issue on how zircon fractionates Zr isotopes in calc-alkaline magmas, the major type of magmas that led to forming continental crust over time. The results also show the great potential of stable Zr isotopes in tracing magmatic thermal and chemical evolution and thus possibly continental crustal differentiation. 
    more » « less
  3. Abstract The southern Coast Mountain batholith was episodically active from Jurassic to Eocene time and experienced four distinct high magmatic flux events during that period. Similar episodicity has been recognized in arcs worldwide, yet the mechanism(s) driving such punctuated magmatic behavior are debated. This study uses zircon Hf and O isotopes, with whole-rock and mineral geochemistry, to track spatiotemporal changes in southern Coast Mountains batholith melt sources and to evaluate models of flare-up behavior and crust formation in Cordilleran arc systems. Zircon Hf isotope analysis yielded consistently primitive values, with all zircon grains recording initial εHf between +6 and +16. The majority (97%) of zircons analyzed yielded δ18O values between 4.2‰ and 6.5‰, and only five grains recorded values of up to 8.3‰. These isotopic results are interpreted to reflect magmatism dominated by mantle melting during all time periods and across all areas of the southern batholith, which argues against the periodic input of more melt-fertile crustal materials as the driver of episodic arc magmatism. They also indicate that limited crustal recycling is needed to produce the large volumes of continental crust generated in the batholith. Although the isotopic character of intrusions is relatively invariant through time, magmas emplaced during flare-ups record higher Sr/Y and La/Yb(N) and lower zircon Ti and Yb concentrations, which is consistent with melting in thickened crust with garnet present as a fractionating phase. Flare-ups are also temporally associated with periods when the southern Coast Mountains batholith both widens and advances inboard. We suggest that the landward shift of the arc into more fertile lithospheric mantle domains triggers voluminous magmatism and is accompanied by magmatic and/or tectonic thickening. Overall, these results demonstrate that the magmatic growth of Cordilleran arcs can be spatially and temporally complex without requiring variability in the contributions of crust and/or mantle to the batholith. 
    more » « less
  4. Abstract This study addresses the question of how and where arc magmas obtain their chemical and isotopic characteristics. The Wooley Creek batholith and Slinkard pluton are a tilted, mid- to upper-crustal part of a vertically extensive, late-Jurassic, arc-related magmatic system in the Klamath Mountains, northern California. The main stage of the system is divided into an older lower zone (c. 159 Ma) emplaced as multiple sheet-like bodies, a younger upper zone (c. 158–156 Ma), which is gradationally zoned upward from mafic tonalite to granite, and a complex central zone, which represents the transition between the lower and upper zones. Xenoliths are common and locally abundant in the lower and central zones and preserve a ghost stratigraphy of the three host terranes. Bulk-rock Nd isotope data along with ages and Hf and oxygen isotope data on zircons were used to assess the location and timing of differentiation and assimilation. Xenoliths display a wide range of εNd (whole-rock) and εHf (zircon), ranges that correlate with rocks in the host terranes. Among individual pluton samples, zircon Hf and oxygen isotope data display ranges too large to represent uniform magma compositions, and very few data are consistent with uncontaminated mantle-derived magma. In addition, zoning of Zr and Hf in augite and hornblende indicates that zircon crystallized at temperatures near or below 800 °C; these temperatures are lower than emplacement temperatures. Therefore, the diversity of zircon isotope compositions reflects in situ crystallization from heterogeneous magmas. On the basis of these and published data, the system is interpreted to reflect initial MASH-zone differentiation, which resulted in elevated δ18O and lowered εHf in the magmas prior to zircon crystallization. Further differentiation, and particularly assimilation–fractional crystallization, occurred at the level of emplacement on a piecemeal (local) basis as individual magma batches interacted with partial melts from host-rock xenoliths. This piecemeal assimilation was accompanied by zircon crystallization, resulting in the heterogeneous isotopic signatures. Magmatism ended with late-stage emplacement of isotopically evolved granitic magmas (c. 156 Ma) whose compositions primarily reflect reworking of the deep-crustal MASH environment. 
    more » « less
  5. null (Ed.)
    Abstract The Duluth Complex (Minnesota, USA) is one of the largest mafic intrusive complexes on Earth. It was emplaced as the Midcontinent Rift developed in Laurentia’s interior during an interval of magmatism and extension from ca. 1109 to 1084 Ma. This duration of magmatic activity is more protracted than is typical for large igneous provinces interpreted to have formed from decompression melting of upwelling mantle plumes. While the overall duration was protracted, there were intervals of more voluminous magmatism. New 206Pb/238U zircon dates for the anorthositic and layered series of the Duluth Complex constrain these units to have been emplaced ca. 1096 Ma in <1 m.y. (duration of 500 ± 260 k.y.). Comparison of paleomagnetic data from these units with Laurentia’s apparent polar wander path supports this interpretation. This rapid emplacement bears similarities to the geologically short duration of well-dated large igneous provinces. These data support hypotheses that call upon the co-location of lithospheric extension and anomalously hot upwelling mantle. This rapid magmatic pulse occurred >10 m.y. after initial magmatism following >20° of latitudinal plate motion. A likely scenario is one in which upwelling mantle encountered the base of Laurentian lithosphere and flowed via “upside-down drainage” to locally thinned lithosphere of the Midcontinent Rift. 
    more » « less