skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: The Origin of Plagiogranites: Coupled SIMS O Isotope Ratios, U–Pb Dating and Trace Element Composition of Zircon from the Troodos Ophiolite, Cyprus

U–Pb ages, trace element content and oxygen isotope ratios of single zircons from five plagiogranite intrusions of the Troodos ophiolite were measured to determine their crystallization age and assess the importance of fractional crystallization versus crustal anatexis in their petrogenesis. The results indicate that oceanic magmatism in Troodos took place at 94·3 ± 0·5 Ma, about 3 Myr earlier than previously recognized. Later hydrothermal alteration has affected most of the Troodos plagiogranitic rocks, resulting in growth of new zircon and/or partial alteration of zircon domains, causing slightly younger apparent crystallization ages. The new age inferred for seafloor spreading and ocean crust accretion in Troodos nearly overlaps that of the Semail ophiolite in Oman (95–96 Ma), strengthening previous indications for simultaneous evolution of both ophiolites in similar tectonic settings. Average δ18O(Zrn) values in the Troodos plagiogranites range between 4·2 and 4·8 ‰. The lower values in this range are lower than those expected in equilibrium with mantle-derived melt (5·3 ± 0·6 ‰), indicating variable contribution from hydrothermally altered, deep-seated oceanic crust in most of the Troodos plagiogranite intrusions. The inferred substantial involvement of crustal component is consistent with the existence of a shallow axial magma chamber, typical of fast-spreading mid-ocean ridge settings, within the Troodos slow-spreading ridge environment. This apparent contradiction may be reconciled by episodically intense magmatism within an otherwise slow, magmatically deprived spreading axis.

more » « less
Award ID(s):
1658823 2004618
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Journal of Petrology
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present a new synthesis of oceanic crustal structure from two‐dimensional seismic profiles to explore differences related to spreading rate and age. Primary results are as follows: (1) Layer 2 has an average thickness of 1.84 km but is thicker for young slow‐spreading crust and thinner for young superfast‐spreading crust. At faster‐spreading rates the layer 2/3 boundary likely corresponds to the lithologic boundary between dikes and gabbros. At slow‐spreading centers, the layer 2/3 boundary is interpreted to mark a change in porosity with depth within the dikes. (2) Total crustal thickness averages 6.15 km and is similar across all spreading rates. (3) Velocities at the top of layer 2 increase rapidly from 3.0 km/s at 0 Ma to 4.6 km/s at 10.5 Ma, with a slower increase to 5.0 km/s at 170 Ma. The rapid increase in velocity at young ages is attributed to crack closure by precipitation of hydrothermal alteration products; the increase at older ages suggests that this process persists as the oceanic crust evolves. (4) There is a correlation between velocities at the top of layer 2 and sediment thickness, with velocities of 5.8–5.9 km/s associated with a sediment thickness of 4.0–4.3 km. The thick sediment may collapse large‐scale features such as lava tubes and fractures. (5) Average velocities at the top of layer 3 are lower for young slow‐spreading and intermediate‐spreading oceanic crust (6.1–6.2 km/s) than for older or faster‐spreading oceanic crust (6.5–6.7 km/s). These low velocities are likely associated with faults penetrating into the sheeted dikes.

    more » « less
  2. Garnet ages for eclogite and granulite from the Western Fiordland Orthogneiss (WFO) provide a precise age for high-grade metamorphism and partial melting of the lower crust in a Cretaceous magmatic arc currently exposed in Fiordland, New Zealand. U/Pb zircon ages and pluton areas indicate that a high magmatic flux event between 118 and 115 Ma added >3,000 km2 of mid- to lower-crustal plutons. The high flux event was followed by high temperature metamorphism and partial melting which resulted in pervasive leucosomes, and trondhjemite layers and veins. At least 1,800 km2 of the newly added crust was metamorphosed to garnet granulite facies orthogneiss. Thermobarometry and phase diagram models indicate that garnet grew at 850 to 1,000°C and 12 to 14 kbar in this monzodiorite and diorite gneiss of the Misty, Malaspina, and Breaksea plutons. Sm-Nd garnet-rock isochrons for these three plutons of the WFO (>700 km2of lower crust) indicate that peak temperatures were reached at 111.7±1.0 Ma (N=16). The isotopic and chemical composition of zircon indicate that the Cretaceous arc flare-up was most likely triggered by partial melting and hybridization of subducted oceanic crust and enriched subcontinental lithospheric mantle directly prior to cessation of arc magmatism. The driving mechanism for the terminal magmatic surge is inferred to be propagation of a discontinuous slab tear beneath the arc, or a ridge-trench collision event between 136 and 128 Ma. The lack of ca. 112 Ma plutons in the western part of Fiordland negates a magmatic heat source for garnet granulite metamorphism. Therefore, we infer that high heat flow associated with mantle advection at the base of the arc after the magmatic surge continued for several m.y., heating the lower crust to granulite facies temperatures 
    more » « less
  3. Abstract

    High‐precision dating of the metamorphic sole of ophiolites can provide insight into the tectonic evolution of ophiolites and subduction zone processes. To understand subduction initiation beneath a young, well‐preserved and well‐characterized ophiolite, we performed coupled zircon laser‐ablation inductively coupled mass spectrometry trace element analyses and high‐precision isotope dilution‐thermal ionization mass spectrometry U–Pb dating on 25 samples from the metamorphic sole of the Samail ophiolite (Oman‐United Arab Emirates). Zircon grains from amphibolite‐ to granulite‐facies (0.8–1.3 GPa, ~700–900°C), garnet‐ and clinopyroxene‐bearing amphibolite samples (n = 18) show systematic trends of decreasing heavy rare earth element slope (HREE; Yb/Dy) with decreasing Yb concentration, reflecting progressive depletion of the HREE during prograde garnet growth. For half of the garnet‐clinopyroxene amphibolite samples, Ti‐in‐zircon temperatures increase, and U–Pb dates young with decreasing HREE slope, consistent with coupled zircon and garnet growth during prograde metamorphism. In the remaining samples, there is no apparent variation in Ti‐in‐zircon temperature with decreasing HREE slope, and the combined U–Pb and geochemical data suggest zircon crystallization along either the prograde to peak or prograde to initial retrograde portions of the metamorphicP–T–tpath. The new data bracket the timing of prograde garnet and zircon growth in the highest grade rocks of the metamorphic sole between 96.698 ± 0.094 and 95.161 ± 0.064 Ma, in contrast with previously published geochronology suggesting prograde metamorphism at ~104 Ma. Garnet‐free amphibolites and leucocratic pods from lower grade (but still upper amphibolite facies) portions of the sole are uniformly HREE enriched (Yb/Dy > 5) and are ~0.5–1.3 Myr younger than the higher grade rocks from the same localities, constraining the temporal offset between the metamorphism and juxtaposition of the higher and lower grade units. Positive zircon εHf(+6.5 to +14.6) for all but one of the dated amphibolites are consistent with an oceanic basalt protolith for the sole. Our new data indicate that prograde sole metamorphism (96.7–95.2 Ma) immediately predated and overlapped growth of the overlying ophiolite crust (96.1–95.2 Ma). The ~600 ky offset between the onset of sole metamorphism in the northern portion of the ophiolite versus the start of ophiolite magmatism is an order of magnitude shorter than previously proposed (~8 Ma) and is consistent with either spontaneous subduction initiation or an abbreviated period of initial thrusting during induced subduction initiation. Taken together, the sole and ophiolite crust preserve a record of the first ~1.5 Myr of subduction. A gradient in the initiation of high‐grade metamorphism from the northwest (96.7 Ma) to southeast (96.0–95.7 Ma) may record propagation of the nascent subduction zone and/or variations in subduction rate along the length of the ophiolite.

    more » « less
  4. Abstract

    Lithospheric seismic anisotropy illuminates mid‐ocean ridge dynamics and the thermal evolution of oceanic plates. We utilize short‐period (5–7.5 s) ambient‐noise surface waves and 15‐ to 150‐s Rayleigh waves measured across the NoMelt ocean‐bottom array to invert for the complete radial and azimuthal anisotropy in the upper ∼35 km of ∼70‐Ma Pacific lithospheric mantle, and azimuthal anisotropy through the underlying asthenosphere. Strong azimuthal variations in Rayleigh‐ and Love‐wave velocity are observed, including the first clearly measured Love‐wave 2θand 4θvariations. Inversion of averaged dispersion requires radial anisotropy in the shallow mantle (2‐3%) and the lower crust (4‐5%), with horizontal velocities (VSH) faster than vertical velocities (VSV). Azimuthal anisotropy is strong in the mantle, with 4.5–6% 2θvariation inVSVwith fast propagation parallel to the fossil‐spreading direction (FSD), and 2–2.5% 4θvariation inVSHwith a fast direction 45° from FSD. The relative behavior of 2θ, 4θ, and radial anisotropy in the mantle are consistent with ophiolite petrofabrics, linking outcrop and surface‐wave length scales.VSVremains fast parallel to FSD to ∼80 km depth where the direction changes, suggesting spreading‐dominated deformation at the ridge. The transition at ∼80 km perhaps marks the dehydration boundary and base of the lithosphere. Azimuthal anisotropy strength increases from the Moho to ∼30 km depth, consistent with flow models of passive upwelling at the ridge. Strong azimuthal anisotropy suggests extremely coherent olivine fabric. Weaker radial anisotropy implies slightly nonhorizontal fabric or the presence of alternative (so‐called E‐type) peridotite fabric. Presence of radial anisotropy in the crust suggests subhorizontal layering and/or shearing during crustal accretion.

    more » « less
  5. Most known porphyry Cu±Au deposits are associated with moderately oxidized and sulfur-rich, calc-alkaline to mildly alkalic arc-related magmas in the Phanerozoic. In contrast, sodium-enriched tonalite–trondhjemite–granodiorite–diorite (TTG) magmas predominant in the Archean are hypothesized to be unoxidized and sulfur-poor, which together preclude porphyry Cu deposit formation. Here, we test this hypothesis by interrogating the causative magmas for the ∼2·7 Ga TTG-related Côté Gold, St-Jude, and Clifford porphyry-type Cu±Au deposit settings in the Neoarchean southern Abitibi subprovince. New and previously published geochronological results constrain the age of emplacement of the causative magmas at ∼2·74 Ga, ∼2·70 Ga, and∼2·69 Ga, respectively. The dioritic and trondhjemitic magmas associated with Côté Gold and St-Jude evolved along a plagioclase-dominated fractionation trend, in contrast to amphibole-dominated fractionation for tonalitic magma at Clifford. Analyses of zircon grains from the Côté Gold, St-Jude, and Clifford igneous rocks yielded εHf(t)±SD values of 4·5±0·3, 4·2±0·6, and 4·3±0·4, and δ18O±SD values of 5·40±0·11  , 3·91±0·13  , and 4·83±0·12  , respectively. These isotopic signatures indicate that, although these magmas are mantle-sourced with minimal crustal contamination, for the St- Jude and Clifford settings the magmas or their sources may have undergone variable alteration by heated seawater or meteoric fluids. Primary barometric minerals (i.e. zircon, amphibole, apatite, and magnetite–ilmenite) that survived variable alteration and metamorphism (up to greenschist facies) were used for estimating fO2 of the causative magmas. Estimation of magmatic fO2 values, reported relative to the fayalite–magnetite–quartz buffer as  FMQ, using zircon geochemistry indicates that the fO2 values of the St-Jude, Côté Gold, and Clifford magmas increase from  FMQ –0·3±0·6 to  FMQ +0·8±0·4 and to  FMQ +1·2±0·4, respectively. In contrast, amphibole chemistry yielded systematically higher fO2 values of  FMQ +1·6±0·3 and  FMQ +2·6±0·1 for Côté Gold and Clifford, respectively, which are consistent with previous studies that indicate that amphibole may overestimate the fO2 of intrusive rocks by up to 1 log unit. Micro X-ray absorption near edge structure (μ-XANES) spectrometric determination of sulfur (i.e. S6+/ S) in primary apatite yielded ≥ FMQ−0·3 and FMQ+1·4–1·8 for St-Jude and Clifford, respectively. The magnetite–ilmenite mineral pairs from the Clifford tonalite yielded  FMQ +3·3±1·3 at equilibrium temperatures of 634±21 ◦C, recording the redox state of the late stage of magma crystallization. Electron probe microanalyses revealed that apatite grains from Clifford are enriched in S (up to 0·1 wt%) relative to those of Côté Gold and St-Jude (below the detection limit), which is attributed to either relatively oxidized or sulfur-rich features of the Clifford tonalite. We interpret these results to indicate that the deposits at Côté Gold and Clifford formed from mildly (∼ FMQ +0·8±0·4) to moderately (∼ FMQ +1·5) oxidized magmas where voluminous early sulfide saturation was probably limited, whereas the St-Jude deposit represents a rare case whereby the ingress of externally derived hydrothermal fluids facilitated metal fertility in a relatively reduced magma chamber (∼ FMQ +0). Furthermore, we conclude that variable modes of formation for these deposits and, in addition, the apparent rarity of porphyry-type Cu–Au deposits in the Archean may be attributed to either local restriction of favorable metallogenic conditions, and/or preservation, or an exploration bias. 
    more » « less