Low-degree melting of a metasomatized lithospheric mantle for the origin of Cenozoic Yulong monzogranite-porphyry,east Tibet: Geochemical and Sr–Nd–Pb–Hf isotopic constraints

SHRIMP zircon U–Pb dating, mineral chemical, element geochemical and Sr–Nd–Pb–Hf isotopic data have been determined for the Yulong monzogranite-porphyry in the eastern Tibet, China. The Yulong porphyry was emplaced into Triassic strata at about 39 Ma. The rocks are weakly peraluminous and show shoshonitic affinity, i.e., alkalis-rich, high K₂O contents with high K₂O / Na₂O ratios, enrichment in LREE and LILE. They also show some affinities with the adakite, e.g., high SiO₂ and Al₂O₃, and low MgO contents, depleted in Y and Yb, and enrichment in Sr with high Sr / Y and La / Yb ratios, and no Eu anomalies. The Yulong porphyry has radiogenic ⁸⁷Sr / ⁸⁶Sr (0.7063–0.7070) and unradiogenic ¹⁴³Nd / ¹⁴⁴Nd (ε_Nd = − 2.0 to − 3.0) ratios. The Pb isotopic compositions of feldspar phenocrysts separated from the Yulong porphyry show a narrow range of ²⁰⁶Pb / ²⁰⁴Pb ratios (18.71–18.82) and unusually radiogenic ²⁰⁷Pb / ²⁰⁴Pb (15.65–15.67) and ²⁰⁸Pb / ²⁰⁴Pb (38.87–39.00) ratios. In situ Hf isotopic composition of zircons that have been SHRIMP U–Pb dated is characterized by clearly positive initial ε_Hf values, ranging from + 3.1 to + 5.9, most between + 4 and + 5. Phenocryst clinopyroxene geothermometry of the Yulong porphyry indicates that the primary magmas had anomalously high temperature (> 1200 °C). The source depth for the Yulong porphyry is at least 100 km inferred by the metasomatic volatile phase (phlogopite–carbonate) relations. Detailed geochemical and Sr–Nd–Pb–Hf isotopic compositions not only rule out fractional crystallization or assimilation-fractional crystallization processes, but also deny the possibility of partial melting of subducted oceanic crust or basaltic lower crust. Instead, low degree (1–5%) partial melting of a metasomatized lithosphere (phlogopite–garnet clinopyroxenite) is compatible with the data. This example gives a case study that granite can be derived directly by partial melting of an enriched lithospheric mantle, which is important to understand the source and origin of diverse granites.     

Os–Hf–Sr–Nd isotope and PGE systematics of spinel peridotite xenoliths from Tok,SE Siberian craton: Effects of pervasive metasomatism in shallow refractory mantle

Os–Hf–Sr–Nd isotopes and PGE were determined in peridotite xenoliths carried to the surface by Quaternary alkali basaltic magmas in the Tokinsky Stanovik Range on the Aldan shield. These data constrain the timing and nature of partial melting and metasomatism in the lithospheric mantle beneath SE Siberian craton. The xenoliths range from the rare fertile spinel lherzolites to the more abundant, strongly metasomatised olivine-rich (70–84%) rocks. Hf–Sr–Nd isotope compositions of the xenoliths are mainly within the fields of oceanic basalts. Most metasomatised xenoliths have lower ¹⁴³Nd / ¹⁴⁴Nd and ¹⁷⁶Hf / ¹⁷⁷Hf and higher ⁸⁷Sr / ⁸⁶Sr than the host basalts indicating that the metasomatism is older and has distinct sources. A few xenoliths have elevated ¹⁷⁶Hf / ¹⁷⁷Hf (up to 0.2838) and plot above the Hf–Nd mantle array defined by oceanic basalts.¹⁸⁷Os / ¹⁸⁸Os in the poorly metasomatised, fertile to moderately refractory (Al₂O₃ ≥ 1.6%) Tok peridotites range from 0.1156 to 0.1282, with oldest rhenium depletion ages being about 2 Ga. The ¹⁸⁷Os / ¹⁸⁸Os in these rocks show good correlations with partial melting indices (e.g. Al₂O₃, modal cpx); the intercept of the Al–¹⁸⁷Os / ¹⁸⁸Os correlation with lowest Al₂O₃ estimates for melting residues (∼0.3–0.5%) has a ¹⁸⁷Os / ¹⁸⁸Os of ∼0.109 suggesting that these peridotites may have experienced melt extraction as early as 2.8 Gy ago. ¹⁸⁷Os / ¹⁸⁸Os in the strongly metasomatised, olivine-rich xenoliths (0.6–1.3% Al₂O₃) ranges from 0.1164 to 0.1275 and shows no apparent links to modal or chemical compositions. Convex-upward REE patterns and high abundances of heavy to middle REE in these refractory rocks indicate equilibration with evolved silicate melts at high melt / rock ratios, which may have also variably elevated their ¹⁸⁷Os / ¹⁸⁸Os. This inference is supported by enrichments in Pd and Pt on chondrite-normalised PGE abundance patterns in some of the rocks. The melt extraction ages for the Tok suite of 2.0 to 2.8 Ga are younger than oldest Os ages reported for central Siberian craton, but they must be considered minimum estimates because of the extensive metasomatism of the most refractory Tok peridotites. This metasomatism could have occurred in the late Mesozoic to early Cenozoic when the Tok region was close to the subduction-related Pacific margin of Siberia and experienced large-scale tectonic and magmatic activity. This study indicates that metasomatic effects on the Re–Os system in the shallow lithospheric mantle can be dramatic.     

Geochemical component relationships in MORB from the Mid-Atlantic Ridge, 22–35°N

New trace element and Hf, Nd, and Pb isotope data are reported for 22 basalts collected between 22°N and 35°N on the Mid-Atlantic Ridge. (La / Sm)_N ratios identify the presence of enriched (E)-MORB in the northernmost part of this area and normal (N)-MORB elsewhere. A negative correlation is observed when ¹⁴³Nd / ¹⁴⁴Nd is plotted against ²⁰⁶Pb / ²⁰⁴Pb, ²⁰⁷Pb / ²⁰⁴Pb, and ²⁰⁸Pb / ²⁰⁴Pb, whereas ¹⁷⁶Hf / ¹⁷⁷Hf appears not to correlate with any of the other isotopic ratios. The E-MORB samples are characterized by high ²⁰⁶Pb / ²⁰⁴Pb, ²⁰⁷Pb / ²⁰⁴Pb, ²⁰⁸Pb / ²⁰⁴Pb, and low ¹⁴³Nd / ¹⁴⁴Nd. Principal Component Analysis (PCA) of Pb isotopes alone identifies three, and only three, significant geochemical end-members (‘components’). Including Nd and Hf isotopic data in the PCA produces spurious components, partly because of curved mixing relationships, and partly because of fractionation during melting. Our preferred interpretation of why ¹⁷⁶Hf / ¹⁷⁷Hf is decoupled from the other isotopic ratios is, as inferred from recent experimental data, that the Hf isotopic compositions of the melt and the residue fail to equilibrate during melting. A strong correlation between (Sr / Nd)_N and (Eu / Eu*)_N indicates that plagioclase is a residual phase of N-MORB, but not of E-MORB melting. The three end-members identified in this study are the depleted mantle, a common-type component, and an enriched plume-type end-member. The common, or ‘C’-type, end-member is characteristic of E-MORB and may itself be a mixture containing recycled oceanic crust (the MORB suite, terrigenous sediments, and/or oceanic plateaus). The plume-type end-member is likely to represent the lower mantle and may involve some primordial material. It is shown that mantle isochrons in general and the Pb–Pb isochron in particular do not characterize a specific geodynamic process acting to create mantle heterogeneities.     

Enriched,HIMU-type peridotite and depleted recycled pyroxenite in the Canary plume: A mixed-up mantle

The Earth's mantle is chemically and isotopically heterogeneous, and a component of recycled oceanic crust is generally suspected in the convecting mantle [Hofmann and White, 1982. Mantle plumes from ancient oceanic crust. Earth Planet. Sci. Lett. 57, 421–436]. Indeed, the HIMU component (high µ = ²³⁸U/²⁰⁴Pb), one of four isotopically distinct end-members in the Earth's mantle, is generally attributed to relatively old (≥ 1–2 Ga) recycled oceanic crust in the form of eclogite/pyroxenite, e.g. [Zindler and Hart, 1986. Chemical geodynamics. Ann. Rev. Earth Planet. Sci. 14, 493–571]. Although the presence of the recycled component is generally supported by element and isotopic data, little is known about its physical state at mantle depths. Here we show that the concentrations of Ni, Mn and Ca in olivine from the Canarian shield stage lavas, which can be used to assess the physical nature of the source material (peridotite versus olivine-free pyroxenite) [Sobolev et al., 2007. The amount of recycled crust in sources of mantle-derived melts. Science 316, 412–417], correlate strongly with bulk rock Sr, Nd and Pb isotopic ratios. The most important result following from our data is that the enriched, HIMU-type (having higher ²⁰⁶Pb/²⁰⁴Pb than generally found in the other mantle end-members) signature of the Canarian hotspot magmas was not caused by a pyroxenite/eclogite constituent of the plume but appears to have been primarily hosted by peridotite. This implies that the old (older than ~ 1 Ga) ocean crust, which has more evolved radiogenic isotope compositions, was stirred into/reacted with the mantle so that there is not significant eclogite left, whereas younger recycled oceanic crust with depleted MORB isotopic signature (< 1 Ga) can be preserved as eclogite, which when melted can generate reaction pyroxenite.     

Hf and Nd isotopes in marine sediments: Constraints on global silicate weathering

The combined use of Lu–Hf and Sm–Nd isotope systems potentially offers a unique perspective for investigating continental erosion, but little is known about whether, and to what extent, the Hf–Nd isotope composition of sediments is related to silicate weathering intensity. In this study, Hf and Nd elemental and isotope data are reported for marine muds, leached Fe-oxide fractions and zircon-rich turbidite sands collected off the Congo River mouth, and from other parts of the SE Atlantic Ocean. All studied samples from the Congo fan (muds, Fe-hydroxides, sands) exhibit indistinguishable Nd isotopic composition (ε_Nd ~ ? 16), indicating that Fe-hydroxides leached from these sediments correspond to continental oxides precipitated within the Congo basin. In marked contrast, Hf isotope compositions for the same samples exhibit significant variations. Leached Fe-hydroxide fractions are characterized by ε_Hf values (from ? 1.1 to + 1.3) far more radiogenic than associated sediments (from ? 7.1 to ? 12.0) and turbidite sands (from ? 27.2 to ? 31.6). ε_Hf values for Congo fan sediments correlate very well with Al/K (i.e. a well-known index for the intensity of chemical weathering in Central Africa). Taken together, these results indicate that (1) silicate weathering on continents leads to erosion products having very distinctive Hf isotope signatures, and (2) a direct relationship exists between ε_Hf of secondary clay minerals and chemical weathering intensity.These results combined with data from the literature have global implications for understanding the Hf–Nd isotope variability in marine precipitates and sediments. Leached Fe-hydroxides from Congo fan sediments plot remarkably well on an extension of the ‘seawater array’ (i.e. the correlation defined by deep-sea Fe–Mn precipitates), providing additional support to the suggestion that the ocean Hf budget is dominated by continental inputs. Fine-grained sediments define a diffuse trend, between that for igneous rocks and the ‘seawater array’, which we refer to as the ‘zircon-free sediment array’ (ε_Hf = 0.91 ε_Nd + 3.10). Finally, we show that the Hf–Nd arrays for seawater, unweathered igneous rocks, zircon-free and zircon-bearing sediments (ε_Hf = 1.80 ε_Nd + 2.35) can all be reconciled, using Monte Carlo simulations, with a simple weathering model of the continental crust.     

Thermal history and origin of the Tanzanian Craton from Pb isotope thermochronology of feldspars from lower crustal xenoliths

Common and radiogenic Pb isotopic compositions of plagioclase and anti-perthitic feldspars from granulite-facies lower crustal xenoliths from the Labait Volcano on the eastern margin of the Tanzanian Craton have been measured via laser ablation MC-ICP-MS. Common Pb in plagioclase and a single stage Pb evolution model indicate that the lower crust of the Tanzanian Craton was extracted from mantle having a ²³⁸U/²⁰⁴Pb of 8.1 ± 0.3 and a ²³²Th/²³⁸U of 4.3 ± 0.1 at 2.71 ± 0.09 Ga (all uncertainties are 2σ). Since 2.4 Ga, some orthoclase domains within anti-perthites have evolved with a maximum ²³⁸U/²⁰⁴Pb of 6 and ²³²Th/²³⁸U of 4.3. The spread in Pb isotopic composition in the anti-perthitic feldspars yields single crystal Pb–Pb isochrons of ～ 2.4 Ga, within uncertainty of U–Pb zircon ages from the same sample suite. The Pb isotopic heterogeneities imply that these granulites resided at temperatures < 600 °C in the lower crust of the Tanzanian Craton from ca. 2.4 Ga to the present. In concert with the chemistry of surface samples, mantle xenoliths, and lower crustal xenoliths, our data imply that the cratonic lithosphere in Tanzania formed ca. ～ 2.7 Ga, in a convergent margin setting, and has remained undisturbed since 2.7 Ga.     

Source heterogeneity for the major components of ∼ 3.7 Ga Banded Iron Formations (Isua Greenstone Belt,Western Greenland): Tracing the nature of interacting water masses in BIF formation

We report trace element, samarium (Sm)–neodymium (Nd) and lead (Pb) isotopic data for individual micro-and mesobands of the Earth's oldest Banded Iron Formation (BIF) from the ∼ 3.7–3.8 Ga Isua Greenstone Belt (IGB, West Greenland) in an attempt to contribute to the characterization of the depositional environment and to the understanding of depositional mechanisms of these earliest chemical sediments. Rare earth element (REE)-yttrium (Y) patterns of the individual mesobands show features of modern seawater with diagnostic cerium (Ce/Ce^⁎), presodymium (Pr/Pr^⁎) and Y/holmium (Ho) anomalies. Very low high field strength elements (HFSE) concentrations indicate essentially detritus-free precipitation. Uranogenic Pb isotope data define a correlation line with a slope of 3691 ± 41 Ma, indicating that the uranium (U)–lead (Pb) system remained closed after the formation of this BIF. High ²⁰⁷Pb/²⁰⁴Pb relative to ²⁰⁶Pb/²⁰⁴Pb ratios compared to average mantle growth evolution models are a feature shared by BIF, penecontemporaneous basalts and clastic volcanogenic metasediments and are indicative of the ultimate high-μ (²³⁸U/²⁰⁴Pb) character of the source region, an essentially mafic Hadean protocrust. Sm–Nd isotopic relations on a layer-by-layer basis point to two REE sources controlling the back-arc basin depositional environment of the BIF, one being seafloor-vented hydrothermal fluids (εNd (3.7 Ga) ∼ + 3.1), the other being ambient surface seawater which reached its composition by erosion of parts of the protocrustal landmass (εNd(3.7 Ga) ∼ + 1.6). The validity of two different and periodically interacting water masses (an essentially two-component mixing system) in the deposition of alternating iron- and silica-rich layers is also reflected by systematic trends in germanium (Ge)/silicon (Si) ratios. These suggest that significant amounts of silica were derived from unexposed and/or destroyed mafic Hadean landmass, unlike iron which probably originated from oceanic crust following hydrothermal alteration by deep percolating seawater. Ge/Si distributional patterns in the early Archean Isua BIF are similar to those reported from the Paleoproterozoic Hamersley (Western Australia) BIF, but overall Ge concentrations are about one order of magnitude higher in the Archean BIF. This seems consistent with other lines of evidence that the ambient Archean seawater was enriched with iron relative to Proterozoic and recent seawater.     

New constraints on the HIMU mantle from neon and helium isotopic compositions of basalts from the Cook–Austral Islands

High ⁴He/³He ratios of 100 000 to 160 000 found at HIMU ocean islands (“high μ,” where μ is the U/Pb ratio) are usually attributed to the presence of recycled oceanic crust in the HIMU mantle source. However, significantly higher ⁴He/³He ratios are expected in recycled crust after residence in the mantle for periods greater than 1 Ga. In order to better understand the helium isotopic signatures in HIMU basalts, we have measured helium and neon isotopic compositions in a suite of geochemically well-characterized basalts from the Cook–Austral Islands. We observe ⁴He/³He ratios ranging from 56 000 to 141 000, suggesting the involvement of mantle reservoirs both more and less radiogenic than the mantle source for mid-ocean ridge basalts (MORBs). In addition, we find that the neon isotopic compositions of HIMU lavas extend from the MORB range to compositions less nucleogenic than MORBs. The Cook-Austral HIMU He–Ne isotopic compositions, along with Sr, Nd, Pb, and Os isotopic compositions, indicate that in addition to recycled crust, a relatively undegassed mantle end-member (e.g., FOZO) is involved in the genesis of these basalts. The association of relatively undegassed mantle material with recycled crust provides an explanation for the close geographical association between HIMU lavas and EM (enriched mantle)-type lavas from this island chain: EM-type signatures represent a higher mixing proportion of relatively undegassed mantle material. Mixing between recycled material and relatively undegassed mantle material may be a natural result of entrainment processes and convective stirring in deep mantle.     

Interaction of adakitic melt-peridotite: Implications for the high-Mg# signature of Mesozoic adakitic rocks in the eastern North China Craton

Abundant dunite and harzbugite xenoliths are preserved in Early Cretaceous high-Mg# [63–67, where Mg# = molar 100 × Mg/(Mg + Fe_tot)] diorite intrusions from western Shandong in the North China Craton (NCC). Dunite and some harzburgite xenoliths typically preserve areas of orthopyroxenite (sometimes accompanied by phlogopite) either as veins or as zones surrounding chromite grains. Harzburgite is chiefly composed of olivine, orthopyroxene, minor clinopyroxene and chromian-spinel. High Mg#'s (averaging 91.4) and depletions in Al₂O₃ and CaO (averaging 0.52 wt.% and 0.29 wt.%, respectively) in harzburgite and dunite xenoliths suggest that they are residues formed by large degrees of polybaric melting. However, olivines and orthopyroxenes from dunite xenoliths spatially associated with orthopyroxenite display lower Mg#'s (i.e., 82–87 and 83–89, respectively), suggesting that an adakitic melt–peridotite reaction has taken place. This is consistent with the production of veined orthopyroxene or orthopyroxene + phlogopite in dunite and some harzburgite xenoliths in response to the introduction of adakitic melt into the previously depleted lithospheric mantle (i.e., harzburgite and dunite xenoliths). The presence of orthopyroxene in veins or as a zones surrounding chromite in peridotite xenoliths is thought to be representative of adakitic melt metasomatism. The dunite and harzbugite xenoliths are relatively rich in light rare earth elements (LREEs) and large ion lithophile elements (LILEs), poor in heavy rare earth elements (HREEs) and high field strength elements (HFSEs), and lack Eu anomalies on chondrite normalized trace element diagrams. The initial ⁸⁷Sr/⁸⁶Sr ratios and εNd(t) values for the xenoliths range from 0.7058 to 0.7212 and + 0.18 to ? 19.59, respectively. Taken together, these features, combined with the strong depletion in HFSE and the existence of Archean inherited zircons in the host rocks, suggest that the adakitic melt was derived from the partial melting of early Mesozoic delaminated lower continental crust. The interaction of the adakitic melt with peridotite is responsible for the high-Mg# character of the early Cretaceous diorites in western Shandong.     

Crustal structure of the Reykjanes Ridge near 62°N,on the basis of seismic refraction and gravity data

Explosion deep seismic sounding data sections of high quality had been obtained with RV Meteor in the Reykjanes Iceland Seismic Project (RRISP77 [Angenheister, G., Gebrande, H., Miller, H., Goldflam, P., Weigel, W., Jacoby, W.R., Pálmason, G., Björnsson, S., Einarsson, P., Pavlenkova, N.I., Zverev, S., Litvinenko, I.V., Loncarecic, B., Solomon, S., 1980. Reykjanes Ridge Iceland Seismic Experiment (RRISP 77). J. Geophys. 47, 228–238]) which close an information gap near 62°N. Preliminary results were presented by Weigel [Weigel, W., 1980. Aufbau des Reykjanes Rückens nach refraktionsseismischen Messungen. In: Weigel, W. (Ed.), Reykjanes Rücken, Island, Norwegischer Kontinentalrand. Abschlusskolloquium, Hamburg zur Meteor-Expedition, vol. 45. DFG, Bonn, pp. 53–61], and here we report on the data and results of interpretation. Clear refracted phases to 90 km distance permit crustal and uppermost mantle structure to be modelled by ray tracing. The apparent P-wave velocities are around 4.5, 6–6.5, 7–7.6 and 8.2–8.7 km/s, but no wide-angle reflections have been clearly seen. Accompanying sparker reflection data reveal thin sediment ponds in the axial zone and up to 400 m thick sediments at 10 Ma crustal age. Ray tracing reveals the following model below the sediments: (1) a distinct, 1–2 km thick upper crust (layer 2A) with Vp increasing with age (to 10 Ma) from <3.4 to 4.9 km/s and with a vertical gradient of 0.1–0.2 km/s/km, (2) a lower crust or layer 3 beginning at depths of 2 (axis) to 4 km (10 Ma age) below sea level with 6.1–6.8 km/s and similar vertical gradients as above, (3) the lower crust bottoms at 5.2–9.5 km depth below sea level (0–10 Ma) with a marked discontinuity, underneath which (4) Vp rises from about 7.5–7.8 km/s (0–10 Ma) with a positive vertical gradient of, again, 0.1–0.2 km/s/km such that 8 km/s would be reached at 12 km and deeper near the axis. Our preferred interpretation is that the mantle begins at the distinct discontinuity (“Moho”), but a deeper “Moho” of Vp ≈ 8 km/s cannot be excluded. From Iceland southward to 60°N several experiments show a decrease of crustal thickness from 14 to 8 km. Velocity trends with age across the ridge reflect cooling and filling of cracks, and thickness trends probably suggest volcanic productivity variations as previously suggested.Gravity inversion concentrates on a profile across the ridge with the above seismic a priori information; with 0.2–0.5 km depth uncertainty it leads to a good fit (±2.5 mGal where seismic data exist). Best fitting densities are (in kg/m³) for sediments, 2180; upper crust, 2450–2570; lower crust, 2850–2940; mantle lithosphere, 3215–3240 with a deficit for an asthenospheric wedge of no more than −100 kg/m³. The morphological ridges and troughs superimposed on the SE ridge flank are partly correlated, partly anti-correlated with the Bouguer anomaly and suggest that variable crustal density variations accompany the morphology variations.     

Balancing the global oceanic neodymium budget: Evaluating the role of groundwater

The distinctly different, ε_Nd(0) values of the Atlantic, Indian, and Pacific Oceans requires that the residence time of Nd in the ocean (i.e., τ_Nd) be on the order of, or less than, the ocean mixing time of ∼ 500–1500 yr. However, estimates of τ_Nd, based on river influxes, range from 4000 to 15,000 yr, thus exceeding the ocean mixing time. In order to reconcile the oceanic Nd budget and lower the residence time by roughly a factor of 10, an additional, as yet unidentified, and hence “missing Nd flux” to the ocean is necessary. Dissolution of materials deposited on continental margins has previously been proposed as a source of the missing flux. In this contribution, submarine groundwater discharge (SGD) is examined as a possible source of the missing Nd flux. Neodymium concentrations (n = 730) and ε_Nd(0) values (n = 58) for groundwaters were obtained from the literature in order to establish representative groundwater values. Mean groundwater Nd concentrations and ε_Nd(0) values were used along with recent estimates of the terrestrial (freshwater) component of SGD (6% of river discharge on a global basis) to test whether groundwater discharge to the coastal oceans could account for the missing flux. Employing mean Nd concentrations of the compiled data base (i.e., 31.8 nmol/kg for all 730 analyses and 11.3 nmol/kg for 141 groundwater samples from a coastal aquifer), the global, terrestrial-derived SGD flux of Nd is estimated to range between 2.9 × 10⁷ and 8.1 × 10⁷ mol/yr. These estimates are of the same order of magnitude, and within a factor of 2, of the missing Nd flux (i.e., 5.4 × 10⁷ mol/yr). Applying the SGD Nd flux estimates, the global average ε_Nd(0) of SGD is predicted to be − 9.1, which is similar to our estimate for the missing Nd flux (− 9.2), and in agreement with the mean (± S.D.) ε_Nd(0) measured in groundwaters (i.e., ε_Nd(0) = −8.9 ± 4.2). The similarities in the estimated SGD Nd flux and corresponding ε_Nd(0) values to the magnitude and isotope composition of the missing Nd flux are compelling, and suggest that discharge of groundwater to the oceans could account for the missing Nd flux. Future investigations should focus on quantifying the Nd concentrations and isotope compositions of groundwater from coastal aquifers from a variety of coastal settings, as well as the important geochemical reactions that effect Nd concentrations in subterranean estuaries in order to better constrain contributions of SGD to the oceanic Nd budget.     

Crustal seismic structure and depth distribution of earthquakes in the Archean Kuusamo region,Fennoscandian Shield

Two-dimensional crustal velocity models are derived from passive seismic observations for the Archean Karelian bedrock of north-eastern Finland. In addition, an updated Moho depth map is constructed by integrating the results of this study with previous data sets. The structural models image a typical three-layer Archean crust, with thickness varying between 40 and 52 km. P wave velocities within the 12–20 km thick upper crust range from 6.1 to 6.4 km/s. The relatively high velocities are related to layered mafic intrusive and volcanic rocks. The middle crust is a fairly homogeneous layer associated with velocities of 6.5–6.8 km/s. The boundary between middle and lower crust is located at depths between 28 and 38 km. The thickness of the lower crust increases from 5–15 km in the Archean part to 15–22 km in the Archean–Proterozoic transition zone. In the lower crust and uppermost mantle, P wave velocities vary between 6.9–7.3 km/s and 7.9–8.2 km/s. The average V_p/V_s ratio increases from 1.71 in the upper crust to 1.76 in the lower crust.The crust attains its maximum thickness in the south-east, where the Archean crust is both over- and underthrust by the Proterozoic crust. A crustal depression bulging out from that zone to the N–NE towards Kuusamo is linked to a collision between major Archean blocks. Further north, crustal thickening under the Salla and Kittilä greenstone belts is tentatively associated with a NW–SE-oriented collision zone or major shear zone. Elevated Moho beneath the Pudasjärvi block is primarily explained with rift-related extension and crustal thinning at ∼2.4–2.1 Ga.The new crustal velocity models and synthetic waveform modelling are used to outline the thickness of the seismogenic layer beneath the temporary Kuusamo seismic network. Lack of seismic activity within the mafic high-velocity body in the uppermost 8 km of crust and relative abundance of mid-crustal, i.e., 14–30 km deep earthquakes are characteristic features of the Kuusamo seismicity. The upper limit of seismicity is attributed to the excess of strong mafic material in the uppermost crust. Comparison with the rheological profiles of the lithosphere, calculated at nearby locations, indicates that the base of the seismogenic layer correlates best with the onset of brittle to ductile transition at about 30 km depth.We found no evidence on microearthquake activity in the lower crust beneath the Archean Karelian craton. However, a data set of relatively well-constrained events extracted from the regional earthquake catalogue implies a deeper cut-off depth for earthquakes in the Norrbotten tectonic province of northern Sweden.     

Age and nature of eclogites in the Huwan shear zone,and the multi-stage evolution of the Qinling-Dabie-Sulu orogen,central China

In situ LA-ICPMS U-Pb, trace element, and Hf isotope data in zircon demonstrate a Carboniferous age for eclogite-facies metamorphism in Siluro-Devonian protoliths in the Huwan shear zone, Dabie Mountains, Central China. This age contrasts with the more prevailing Triassic age for high- to ultrahigh pressure (HP to UHP) metamorphism in the Qinling-Dabie-Sulu orogen. Metamorphic zircon in two eclogite samples from Sujiahe is characterized by low Th/U ratios, small negative Eu anomalies, flat HREE patterns, and low ¹⁷⁶Lu/¹⁷⁷Hf ratios. These geochemical signatures suggest that the zircon crystallized in the presence of garnet and in the absence of plagioclase feldspar. Furthermore, temperatures of ~ 655 and ~ 638 °C, calculated using the Ti content of zircon, are consistent with their formation during eclogite-facies metamorphism. The weighted mean ²⁰⁶Pb/²³⁸U age of 309 ± 4 Ma (2δ) for this zircon improves previous age estimates for eclogite-facies metamorphism in the Huwan shear zone, ranging from 420 to 220 Ma. Metamorphic zircon from one eclogite sample from Hujiawan, most likely formed during prograde metamorphism, yields an equivalent age estimate of 312 ± 11 Ma. Magmatic zircon cores in the three samples yield ages for the magmatic protoliths of the eclogites ranging from 420 ± 7 to 406 ± 5 Ma, and post-dating the middle Paleozoic collision of the North China and the Qinling terrain. The zircon crystals in the three eclogite samples display a large variation of ε_Hf (t) values of ? 4.9 to 21.3. The metamorphic zircon overgrowths show the same range of ε_Hf (t) values as those of the inherited magmatic crystal interiors. This suggests that the metamorphic zircon overgrowths may have formed by dissolution-reprecipitation of pre-existing magmatic zircon thereby preserving their original Hf isotopic composition. The high ε_Hf (t) values suggest that the protoliths were derived from depleted mantle sources, most likely Paleotethyan oceanic crust; while the low ε_Hf (t) values are attributed to crustal contamination. Some eclogites in the Huwan shear zone have a distinctive signature of continental crust most probably derived from the Yangtze Craton. The coexistence of Paleozoic oceanic crust and Neoproterozoic continental crust with similar metamorphic ages in the Huwan shear zone implies that Paleozoic Paleotethyan oceanic crust was produced within a marginal basin of the northern Yangtze Craton. The opening of the Paleo-Tethyan ocean was slightly younger than the collision of the North China Craton and the Qinling terrain during the Late Paleozoic in the Qinling-Dabie-Sulu orogen. Subduction of the Paleo-Tethyan oceanic crust and associated continental basement resulted in the 309 ± 2 Ma (2σ) eclogite-facies metamorphism in the Huwan shear zone. The subsequent Triassic continent-continent collision led to the final coalescence of the Yangtze and Sino-Korean cratons. Amalgamation of the Yangtze and North China cratons was, therefore, a multistage process extending over at least 200 Ma.     

Rhenium–osmium isotope and elemental behaviour during subduction of oceanic crust and the implications for mantle recycling

This study presents major-, trace-element, and rhenium–osmium (Re–Os) isotope and elemental data for basalts and gabbros from the Zermatt-Saas ophiolite, metamorphosed to eclogite-facies conditions during the Alpine orogeny. Igneous crystallisation of the gabbros occurred at 163.5 ± 1.8 Ma and both gabbro and basalt were subject to ‘peak’ pressure–temperature (P–T) conditions of > 2.0 GPa and ~ 600 °C at about 40.6 ± 2.6 Ma.Despite such extreme P–T conditions, Re–Os isotope and abundance data for gabbroic rocks suggest that there has been no significant loss of either of these elements during eclogite-facies metamorphism. Indeed, ¹⁸⁷Re–¹⁸⁷Os isotope data for both unaltered gabbros and gabbroic eclogites lie on the same best-fit line corresponding to an errorchron age of 160 ± 6 Ma, indistinguishable from the age of igneous crystallisation. In contrast, metamorphosed basalts do not yield age information; rather most possess ¹⁸⁷Re/¹⁸⁸Os ratios that cannot account for the measured ¹⁸⁷Os/¹⁸⁸Os ratios, given the time since igneous crystallisation. Taken with their low Re contents these data indicate that the basalts have experienced significant Re loss (∼ 50–60%), probably during high-pressure metamorphism. Barium, Rb and K are depleted in both gabbroic and basaltic eclogites. In contrast, there is no evident depletion of U in either lithology.Many ocean-island basalts (OIB) possess radiogenic Os and Pb isotope compositions that have been attributed to the presence of recycled oceanic crust in the mantle source. Published Re–Os data for high-P metabasaltic rocks alone (consistent with this study) have been taken to suggest that excessive amounts of oceanic crust are required to generate such signatures. However, this study shows that gabbro may exert a strong influence on the composition of recycled oceanic crust. Using both gabbro and basalt (i.e. a complete section of oceanic crust) calculations suggest that the presence of ≥ 40% of 2 Ga oceanic crust can generate the radiogenic Os compositions seen in some OIB. Furthermore, lower U/Pb ratios in gabbro (compared to basalt) serve to limit the ²⁰⁶Pb/²⁰⁴Pb ratios generated, while having a minimal effect on Os ratios. These results suggest that the incorporation of gabbro into recycling models provides a means of producing a range of OIB compositions having lower (and variable) ²⁰⁶Pb/²⁰⁴Pb ratios, but still preserving ¹⁸⁷Os/¹⁸⁸Os compositions comparable to HIMU-type OIB.     

Zircon U–Pb strain chronometry reveals deep impact-triggered flow

Large (> 100 km) meteorite impact cratering events play important roles in surface and biosphere evolution, however, their potential for widespread ductile modification of the lithosphere has been difficult to assess, due partly to our inability to isotopically age-correlate deep mineral fabrics with surface records. We have integrated benchmark U–Pb zircon dating methods (ID-TIMS, SHRIMP) with new microstructural techniques (EBSD, µXRD) to demonstrate that crystal–plastic deformation can cause rapid out-diffusion of radiogenic Pb and accompanying trace element alteration in crystalline zircon. We have used this phenomenon to directly date fabric in Archean zircons and xenoliths of the lower crust of South Africa at 2023 ± 15 million years, coeval with the 2020 ± 3 million year old Vredefort cratering event at surface, with extent ≥ 20,000 km². Our findings indicate that regional exogenic fabrics, similar to high-temperature tectonic fabrics, exist in ancient crust. Moreover, our results establish that crystal-plastic deformation in the lithosphere can now be directly dated and linked to planetary evolution by zircon U–Pb strain chronometry.     

Miocene to Late Quaternary Patagonian basalts (46–47°S): Geochronometric and geochemical evidence for slab tearing due to active spreading ridge subduction

Miocene to Quaternary large basaltic plateaus occur in the back-arc domain of the Andean chain in Patagonia. They are thought to result from the ascent of subslab asthenospheric magmas through slab windows generated from subducted segments of the South Chile Ridge (SCR). We have investigated three volcanic centres from the Lago General Carrera–Buenos Aires area (46–47°S) located above the inferred position of the slab window corresponding to a segment subducted 6 Ma ago. (1) The Quaternary Río Murta transitional basalts display major, trace elements, and Sr and Nd isotopic features similar to those of oceanic basalts from the SCR and from the Chile Triple Junction near Taitao Peninsula (e.g., (⁸⁷Sr/⁸⁶Sr)_o = 0.70396–0.70346 and εNd = + 5.5 − + 3.0). We consider them as derived from the melting of a Chile Ridge asthenospheric mantle source containing a weak subduction component. (2) The Plio-Quaternary (< 3.3 Ma) post-plateau basanites from Meseta del Lago Buenos Aires (MLBA), Argentina, likely derive from small degrees of melting of OIB-type mantle sources involving the subslab asthenosphere and the enriched subcontinental lithospheric mantle. (3) The main plateau basaltic volcanism in this region is represented by the 12.4–3.3-Ma-old MLBA basalts and the 8.2–4.4-Ma-old basalts from Meseta Chile Chico (MCC), Chile. Two groups can be distinguished among these main plateau basalts. The first group includes alkali basalts and trachybasalts displaying typical OIB signatures and thought to derive from predominantly asthenospheric mantle sources similar to those of the post-plateau MLBA basalts, but through slightly larger degrees of melting. The second one, although still dominantly alkalic, displays incompatible element signatures intermediate between those of OIB and arc magmas (e.g., La/Nb > 1 and TiO₂ < 2 wt.%). These intermediate basalts differ from their strictly alkalic equivalents by having lower High Field Strength Element (HFSE) and higher εNd (up to + 5.4). These features are consistent with their derivation from an enriched mantle source contaminated by ca. 10% rutile-bearing restite of altered oceanic crust. The petrogenesis of the studied Mio-Pliocene basalts from MLBA and MCC is consistent with contributions of the subslab asthenosphere, the South American subcontinental lithospheric mantle and the subducted Pacific oceanic crust to their sources. However, their chronology of emplacement is not consistent with an ascent through an asthenospheric window opened as a consequence of the subduction of segment SCR-1, which entered the trench at 6 Ma. Indeed, magmatic activity was already important between 12 and 8 Ma in MLBA and MCC as well as in southernmost plateaus, i.e., 6 Ma before the subduction of the SCR-1 segment. We propose a geodynamic model in which OIB and intermediate magmas derived from deep subslab asthenospheric mantle did uprise through a tear-in-the-slab, which formed when the southernmost segments of the SCR collided with the Chile Trench around 15 Ma. During their ascent, they interacted with the Patagonian supraslab mantle and, locally, with slivers of subducted Pacific oceanic crust that contributed to the geochemical signature of the intermediate basalts.     

Osmium isotope and highly siderophile element systematics of the lunar crust

Coupled ¹⁸⁷Os/¹⁸⁸Os and highly siderophile element (HSE: Os, Ir, Ru, Pt, Pd, and Re) abundance data are reported for pristine lunar crustal rocks 60025, 62255, 65315 (ferroan anorthosites, FAN) and 76535, 78235, 77215 and a norite clast in 15455 (magnesian-suite rocks, MGS). Osmium isotopes permit more refined discrimination than previously possible of samples that have been contaminated by meteoritic additions and the new results show that some rocks, previously identified as pristine, contain meteorite-derived HSE. Low HSE abundances in FAN and MGS rocks are consistent with derivation from a strongly HSE-depleted lunar mantle. At the time of formation, the lunar floatation crust, represented by FAN, had 1.4 ± 0.3 pg g^? 1 Os, 1.5 ± 0.6 pg g^? 1 Ir, 6.8 ± 2.7 pg g^? 1 Ru, 16 ± 15 pg g^? 1 Pt, 33 ± 30 pg g^? 1 Pd and 0.29 ± 0.10 pg g^? 1 Re (～ 0.00002 × CI) and Re/Os ratios that were modestly elevated (¹⁸⁷Re/¹⁸⁸Os = 0.6 to 1.7) relative to CI chondrites. MGS samples are, on average, characterised by more elevated HSE abundances (～ 0.00007 × CI) compared with FAN. This either reflects contrasting mantle-source HSE characteristics of FAN and MGS rocks, or different mantle–crust HSE fractionation behaviour during production of these lithologies. Previous studies of lunar impact-melt rocks have identified possible elevated Ru and Pd in lunar crustal target rocks. The new results provide no supporting evidence for such enrichments.If maximum estimates for HSE in the lunar mantle are compared with FAN and MGS averages, crust–mantle concentration ratios (D-values) must be ≤ 0.3. Such D-values are broadly similar to those estimated for partitioning between the terrestrial crust and upper mantle, with the notable exception of Re. Given the presumably completely different mode of origin for the primary lunar floatation crust and tertiary terrestrial continental crust, the potential similarities in crust–mantle HSE partitioning for the Earth and Moon are somewhat surprising. Low HSE abundances in the lunar crust, coupled with estimates of HSE concentrations in the lunar mantle implies there may be a ‘missing component’ of late-accreted materials (as much as 95%) to the Moon if the Earth/Moon mass-flux estimates are correct and terrestrial mantle HSE abundances were established by late accretion.     

The 3.26–3.24 Ga Barberton asteroid impact cluster: Tests of tectonic and magmatic consequences,Pilbara Craton,Western Australia

The location in the Barberton Greenstone Belt (Kaapvaal Craton) of ∼3.26–3.24 Ga asteroid impact ejecta units at, and immediately above, a sharp break between a > 12 km-thick mafic–ultramafic volcanic crust (Onverwacht Group ∼3.55–3.26 Ga, including the ∼3.298 > 3.258 Ga Mendon Formation) and a turbidite–felsic volcanic rift-facies association (Fig Tree Group ∼3.258–3.225 Ga), potentially represents the first documented example of cause–effect relations between extraterrestrial bombardment and major tectonic and igneous events [D.R. Lowe, G.R. Byerly, F. Asaro, F.T. Kyte, Geological and geochemical record of 3400 Ma old terrestrial meteorite impacts, Science 245 (1989) 959–962; D.R. Lowe, G.R. Byerly, F.T. Kyte, A. Shukolyukov, F. Asaro, A. Krull, Spherule beds 3.47–3.34 Ga-old in the Barberton greenstone belt, South Africa: a record of large meteorite impacts and their influence on early crustal and biological evolution, Astrobiology 3 (2003) 7–48; A.Y. Glikson, The astronomical connection of terrestrial evolution: crustal effects of post-3.8 Ga mega-impact clusters and evidence for major 3.2 ± 0.1 Ga bombardment of the Earth–Moon system, J. Geodyn. 32 (2001) 205–229]. Here we correlate this boundary with a contemporaneous break and peak magmatic and faulting events in the Pilbara Craton, represented by the truncation of a 3.255–3.235 Ga-old volcanic sequence (Sulphur Springs Group—SSG) by a turbidite-banded iron formation–felsic volcanic association (Pincunah Hill Formation, basal Gorge Creek Group). These events are accompanied by ∼3.252–3.235 Ga granitoids (Cleland plutonic suite). The top of the komatiite–tholeiite–rhyolite sequence of the SSG is associated with a marker chert defined at 3.238 ± 3–3.235 ± 3 Ga, abruptly overlain by an olistostrome consisting of mega-clasts of felsic volcanics, chert and siltstone up to 250 × 150 m-large, intercalated with siliciclastic sedimentary rocks and felsic volcanics (Pincunah Hill Formation-basal Gorge Creek Group-GCG [R. M. Hill, Stratigraphy, structure and alteration of hanging wall sedimentary rocks at the Sulphur Springs volcanogenic massive sulphide (VMS) prospect, east Pilbara Craton, Western Australia. B.Sc Hon. Thesis, University of Western Australia (1997) 67 pp.; M.J. Van Kranendonk, A.H. Hickman, R.H. Smithies, D.R. Nelson, Geology and tectonic evolution of the Archaean north Pilbara terrain, Pilbara Craton, Western Australia, Econ. Geol. 97 (2002) 695–732; M.J. Van Kranendonk, Geology of the North Shaw 1 : 100 000 Sheet. Geological Survey Western Australia 1 : 100 000 Geological Series (2000) 86 pp., R. Buick, C.A.W. Brauhart, P. Morant, J.R. Thornett, J.G. Maniew, J.G. Archibald, M.G. Doepel, I.R. Fletcher, A.L. Pickard, J.B. Smith, M.B. Barley, N.J. McNaughton, D.I. Groves, Geochronology and stratigraphic relations of the Sulphur Springs Group and Strelley Granite: a temporally distinct igneous province in the Archaean Pilbara Craton, Australia, Precambrian Res. 114 (2002) 87–120]). The structure and scale of the olistostrome, not seen elsewhere in the Pilbara Craton, is interpreted in terms of intense faulting and rifting, supported by topographic relief represented by deep incision of overlying arenites (Corboy Formation) into underlying units [M.J. Van Kranendonk, Geology of the North Shaw 1 : 100 000 Sheet. Geological Survey Western Australia 1 : 100 000 Geological Series (2000) 86 pp.]. The age overlaps between (1) 3.255 ± 4–3.235 ± 3 Ga peak igneous activity represented by the SSG and the Cleland plutonic suite (Pilbara Craton) and the 3.258 ± 3 Ga S2 Barberton impact unit, and (2) 3.235 ± 3 Ga top SSG break and associated faulting and the 3.243 ± 4 S3–S4 Barberton impact units may not be accidental. Should correlations between the Barberton S2–S4 impact units and magmatic and tectonic events in the Pilbara Craton be confirmed, they would imply impact-triggered reactivation of mantle convection, crustal anatexis, faulting and strong vertical movements in Archaean granite–greenstone terrains associated with large asteroid impacts, culminating in transformation from sima-dominated crust to continental rift environments.     

The Lhasa Terrane: Record of a microcontinent and its histories of drift and growth

The Lhasa Terrane in southern Tibet has long been accepted as the last geological block accreted to Eurasia before its collision with the northward drifting Indian continent in the Cenozoic, but its lithospheric architecture, drift and growth histories and the nature of its northern suture with Eurasia via the Qiangtang Terrane remain enigmatic. Using zircon in situ U–Pb and Lu–Hf isotopic and bulk-rock geochemical data of Mesozoic–Early Tertiary magmatic rocks sampled along four north–south traverses across the Lhasa Terrane, we show that the Lhasa Terrane has ancient basement rocks of Proterozoic and Archean ages (up to 2870 Ma) in its centre with younger and juvenile crust (Phanerozoic) accreted towards its both northern and southern edges. This finding proves that the central Lhasa subterrane was once a microcontinent. This continent has survived from its long journey across the Paleo-Tethyan Ocean basins and has grown at the edges through magmatism resulting from oceanic lithosphere subduction towards beneath it during its journey and subsequent collisions with the Qiangtang Terrane to the north and with the Indian continent to the south. Zircon Hf isotope data indicate significant mantle source contributions to the generation of these granitoid rocks (e.g., ~ 50–90%, 0–70%, and 30–100% to the Mesozoic magmatism in the southern, central, and northern Lhasa subterranes, respectively). We suggest that much of the Mesozoic magmatism in the Lhasa Terrane may be associated with the southward Bangong–Nujiang Tethyan seafloor subduction beneath the Lhasa Terrane, which likely began in the Middle Permian (or earlier) and ceased in the late Early Cretaceous, and that the significant changes of zircon ε_Hf(t) at ~ 113 and ~ 52 Ma record tectonomagmatic activities as a result of slab break-off and related mantle melting events following the Qiangtang–Lhasa amalgamation and India–Lhasa amalgamation, respectively. These results manifest the efficacy of zircons as a chronometer (U–Pb dating) and a geochemical tracer (Hf isotopes) in understanding the origin and histories of lithospheric plates and in revealing the tectonic evolution of old orogenies in the context of plate tectonics.     

Evaluation of the 87Rb decay constant by age comparison against the U–Pb system

The Rb–Sr decay system is one of the most widely used geochronometers for obtaining ages and cooling rates of terrestrial magmatic, metamorphic, and hydrothermal events. It has also been extensively applied to date extraterrestrial, early solar system events. The accuracy of Rb–Sr ages, however, strongly depends on the accuracy of the ⁸⁷Rb decay constant (λ⁸⁷Rb). We determined λ⁸⁷Rb relative to the decay constants of ²³⁵U and ²³⁸U by comparing Rb–Sr ages of minerals with U–Pb ages obtained from the same intrusion. Comparison of U–Pb emplacement ages with high-precision Rb–Sr mineral ages from three rapidly cooled igneous rocks covering an age range of ca. 2.5 Ga yields an unweighted mean λ⁸⁷Rb of 1.393 ± 0.004 × 10^?11 yr^?1 (i.e., ± 0.3%), corresponding to a half-life of 49.76 × 10⁹ years. Because this decay constant is 2% lower than the presently recommended one, many previously published ages are 2% too young and the resulting geologic interpretations may need revision.