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.     

Lithospheric contributions to high-MgO basanites from the Cumbre Vieja Volcano,La Palma,Canary Islands and evidence for temporal variation in plume influence

New geochemical and isotopic data are presented from the oldest part of the Cumbre Vieja volcano, La Palma (Canary Islands), located near the assumed emergence of the Canary mantle plume. The volcanics comprise a suite dominated by basanite flows with subordinate amounts of phono-tephrite, tephri-phonolite and phonolite flows and intrusives. Two compositionally different basanite groups have been identified, both with HIMU (high-μ)-type incompatible trace element characteristics: Primitive high-MgO basanites (10.7–12.1% MgO), found only at the base of a stratigraphic profile near Fuencaliente on the south coast, and intermediate-MgO basanites (6.0–7.3% MgO), exposed in the upper part of the profile and widespread on the east coast of La Palma. The high-MgO basanites are interpreted as near-primary mantle melts (primary composition 14–15% MgO) derived by progressive melting (2.9% to 4.5%) of a common lithospheric mantle source. Model calculations indicate that it is not possible to generate the intermediate-MgO basanites from the high-MgO group by crystal fractionation of observed phenocrysts. Relative to intermediate-MgO basanites, the high-MgO flows have lower concentrations of LIL and HFS elements, except for Ti, which is markedly enriched in the primitive rocks (3.7–4.7% TiO₂ vs 3.4–3.9% TiO₂). Fuencaliente volcanics display limited temporal isotopic variations suggested to be a result of mixing of melts originating from the rising plume and the metazomatized lithospheric mantle. ⁸⁷Sr / ⁸⁶Sr and ¹⁴³Nd / ¹⁴⁴Nd ratios range 0.70305–0.70311 and 0.51285–0.51291, respectively, while the corresponding ranges in Pb-isotope ratios are ²⁰⁶Pb / ²⁰⁴Pb = 19.46–19.64, ²⁰⁷Pb / ²⁰⁴Pb = 15.55–15.61, and ²⁰⁸Pb / ²⁰⁴Pb = 39.16–39.53. The overall variation of the Cumbre Vieja isotopic data can be accounted for by mixtures of three mantle components in the proportions 72–79% plume source (LVC = low velocity component), 9–16% depleted mantle (DM) and up to 12% enriched mantle (EMI). Negative Δ7 / 4 Pb (− 0.6 to − 5.4) in the Cumbre Vieja volcanics suggest derivation from a young HIMU mantle source. The relative abundance of plume source material increase in younger rocks in the Fuencaliente section, suggesting waning plume–lithosphere interaction during the emplacement of this part of the Cumbre Vieja volcano. The high-MgO volcanics define regular and systematic geochemical trends, interpreted as partial melting trends, when plotted against abundances of highly incompatible elements (P, Ce). Evaluation of minor and trace element variation in consecutive melts suggests control by residual amphibole, phlogopite, garnet and a Ti-bearing phase, possibly ilmenite. The melting mode changed gradually, allowing increasing input from residual phlogopite during partial melting. The residual mineralogy constrains the source region of the high-MgO basanites to the lowermost oceanic lithospheric mantle, presumably around 100 km depths.     

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.     

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.     

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.     

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.     

Relation between mantle shear zone deformation and metasomatism in spinel peridotite xenoliths of Jeju Island (South Korea): Evidence from olivine CPO and trace elements

This study presents geochemical and fabric analysis of coarse-grained, porphyroclastic and mylonitic spinel peridotite xenoliths derived from the shallow subcontinental lithospheric mantle of Jeju Island (South Korea). Fabric analysis of olivines in the studied peridotites shows activation of the high temperature (0 1 0)[1 0 0] slip system; however, in the porphyroclastic and mylonitic peridotites, activation of (0 k l)[1 0 0] slip system results in a weaker fabric. Formation of porphyroclastic and mylonitic fabrics are thought to form in a shear-dominated environment. The results of the trace element analysis reveal that the smaller the grain size and weaker the fabric, the more enriched in LREE and HFSE are the peridotites, which indicates a strong relationship between metasomatic agents and mantle shear zones.     

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.     

The formation of volcanic centers at the Colorado Plateau as a result of the passage of aqueous fluid through the oceanic lithosphere and the subcontinental mantle: New implications for the planetary water cycle in the western United States

We provide new petrological evidence for the strong influence of water on the formation of the oceanic lithospheric mantle, the subcontinental mantle above, and the continental lithosphere. Our analysis throws new light on the hypothesis that new continental lithosphere was formed by the passage of silicate-rich aqueous fluid through the sub-continental mantle. In order to investigate this hypothesis, we analyzed a representative collection of lherzolite and harzburgite xenoliths from the sample volcano known as “The Thumb”, located in the center of the Colorado Plateau, western United States. The studied sample collection exhibits multi-stage water enrichment processes along point, line and planar defect structures in nominally anhydrous minerals and the subsequent formation of the serpentine polymorph antigorite along grain boundaries and in totally embedded annealed cracks. Planar defect structures act like monomineralic and interphase grain boundaries in the oceanic lithosphere and the subcontinental mantle beneath the North American plate, which was hydrated by the ancient oceanic Farallon plate during the Cenozoic and Mesozoic eras. We used microspectroscopical, petrological, and seismological techniques to confirm multi-stage hydration from a depth of ∼150 km to just below the Moho depth. High-resolution mapping of the water distribution over homogeneous areas and fully embedded point, line and planar defects in olivine crystals of lherzolitic and harzburgitic origin by synchrotron infrared microspectroscopy enabled us to resolve local wet spots and thus reconstruct the hydration process occurring at a depth of ∼150 km (T ≈ 1225 °C). These lherzolites originated from the middle part of the Farallon mantle slab; they were released during the break up of the Farallon mantle slab, caused by the instability of the dipping slab. The background hydration levels in homogeneous olivines reached ∼138 ppm wt H₂O, and the water concentration at the planar defects could reach up to ∼1000 ppm wt H₂O. However, the formation of antigorite in grain boundaries was found to be the primary hydration mechanism for harzburgitic samples originating from the subcontinental mantle (for hydration, T ≈ 600 °C). Additionally, the formation of antigorite in lherzolites could be found in annealed cracks. From these observations, we conclude that hydration induces multi-stage water enrichment of the mantle wedge by a process that is dominated by the growth and movement of ubiquitous cracks, which acts as planar defects. Cracks in the mantle seem to be the an important feature in both the water cycle of the subduction zone and the formation of the continental lithosphere.     

On the Vp/Vs–Mg# correlation in mantle peridotites: Implications for the identification of thermal and compositional anomalies in the upper mantle

We use thermodynamically self-consistent and hybrid methods to analyze the correlation of important physical parameters (e.g. bulk density, elastic moduli) with bulk Mg# and modal composition in mantle peridotites at upper mantle conditions. Temperature (anharmonic and anelastic), pressure and compositional derivatives for all these parameters are evaluated. The results show that the widely used correlations between Vp/Vs and Mg# in peridotites are strictly valid only for garnet-bearing assemblages at temperatures < 900 °C. The correlation breaks down when: i) spinel is the stable Al-rich phase in the assemblage and ii) when anelastic attenuation of seismic velocities becomes important (T ? 900 °C). This implies that the range of applicability of published Vp/Vs–Mg# correlations for the upper mantle is limited to a depth interval between the spinel–garnet phase transition and the 900 °C isotherm. We use numerical simulations to show that this depth interval is virtually nonexistent in lithospheres thinner than ～ 140 km and can comprise up to ～ 50% of the lithospheric mantle in thick (> 220 km) lithospheric domains. In addition, we show that for most of the upper mantle the expected Δ(Vp/Vs) values associated with compositional variations are smaller than the resolution limit of current seismological methods. All these considerations suggest that the Vp/Vs ratio is not a reliable measure of compositional variations and that for large parts of the upper mantle compositional anomalies cannot be separated from thermal anomalies on the basis of seismological studies only. We further confirm that the only reliable indicator of compositional anomalies in a peridotitic mantle is the ratio of density to shear wave velocities (ρ/Vs). Our results demonstrate that geophysical–petrological models (forward or inverse) that model these two fields (i.e. density and Vs) self-consistently within a robust thermodynamic framework are necessary for characterizing the small-scale thermal and compositional structure of the lithosphere and sublithospheric upper mantle.     

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.     

Oxygen isotope constraints on the origin and differentiation of the Moon

We report new high-precision laser fluorination three-isotope oxygen data for lunar materials. Terrestrial silicates with a range of δ¹⁸O values (− 0.5 to 22.9‰) were analyzed to independently determine the slope of the terrestrial fractionation line (TFL; λ = 0.5259 ± 0.0008; 95% confidence level). This new TFL determination allows direct comparison of lunar oxygen isotope systematics with those of Earth. Values of Δ¹⁷O for Apollo 12, 15, and 17 basalts and Luna 24 soil samples average 0.01‰ and are indistinguishable from the TFL. The δ¹⁸O values of high- and low-Ti lunar basalts are distinct. Average whole-rock δ¹⁸O values for low-Ti lunar basalts from the Apollo 12 (5.72 ± 0.06‰) and Apollo 15 landing sites (5.65 ± 0.12‰) are identical within error and are markedly higher than Apollo 17 high-Ti basalts (5.46 ± 0.11‰). Evolved low-Ti LaPaz mare-basalt meteorite δ¹⁸O values (5.67 ± 0.05‰) are in close agreement with more primitive low-Ti Apollo 12 and 15 mare basalts. Modeling of lunar mare-basalt source composition indicates that the high- and low-Ti mare-basalt mantle reservoirs were in oxygen isotope equilibrium and that variations in δ¹⁸O do not result from fractional crystallization. Instead, these differences are consistent with mineralogically heterogeneous mantle sources for mare basalts, and with lunar magma ocean differentiation models that result in a thick feldspathic crust, an olivine–pyroxene-rich mantle, and late-stage ilmenite-rich zones that were convectively mixed into deeper portions of the lunar mantle. Higher average δ¹⁸O (WR) values of low-Ti basalts compared to terrestrial mid ocean ridge basalts (Δ=0.18‰) suggest a possible oxygen isotopic difference between the terrestrial and lunar mantles. However, calculations of the δ¹⁸O of lunar mantle olivine in this study are only 0.05‰ higher than terrestrial mantle olivine. These observations may have important implications for understanding the formation of the Earth–Moon system.     

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.     

The thermal state and strength of the lithosphere in the Spanish Central System and Tajo Basin from crustal heat production and thermal isostasy

In this work we have modeled the thermal structure of the lithosphere of the Spanish Central System and the Tajo Basin, and their implications for lithospheric strength. For this, we have used refined heat-producing elements (HPE) values to obtain new estimates of heat production rates in the Spanish Central System and Tajo Basin areas, which have been used joined to the relation between topography and thermal structure of the lithosphere to calculate the best-fit surface heat flows in the study area. Moreover, we have implemented a temperature-dependent thermal conductivity (appropriate for olivine) for the lithospheric mantle to improve the calculations of temperature profiles in the mantle. The geotherms so obtained, together with the implementation of a new rheological law for the upper lithospheric mantle, have been used to calculate refined estimations of the strength and effective elastic thickness of the lithosphere. We have obtained surface heat flow values of 84 mW m⁻² and ∼82 mW m⁻² for the Spanish Central System and the Tajo Basin, respectively. The thermal state of the lithosphere affects mantle temperatures, and hence may be playing an important role in the uplift and maintenance of the Spanish Central System.     

Ancient Os isotope signatures from the Ontong Java Plateau lithosphere: Tracing lithospheric accretion history

In order to better understand the nature and formation of oceanic lithosphere beneath the Early Cretaceous Ontong Java Plateau, Re–Os isotopes have been analysed in a suite of peridotite xenoliths from Malaita, Solomon Islands. Geological, thermobarometric and petrological evidence from previous studies reveal that the xenoliths represent virtually the entire thickness of the southern part of subplateau lithospheric mantle (< 120 km). This study demonstrates that vertical Os isotopic variations correlate with compositional variations in a stratified lithosphere. The shallowest plateau lithosphere (< 85 km) is dominated by fertile lherzolites showing a restricted range of ¹⁸⁷Os/¹⁸⁸Os (0.1222 to 0.1288), consistent with an origin from ~ 160 Ma Pacific lithosphere. In contrast, the basal section of subplateau lithospheric mantle (~ 95–120 km) is enriched in refractory harzburgites with highly unradiogenic ¹⁸⁷Os/¹⁸⁸Os ratios ranging from 0.1152 to 0.1196, which yield Proterozoic model ages of 0.9–1.7 Ga. Although the whole range of Os isotope compositions of Malaita peridotites is within the variations seen in modern abyssal peridotites, the contrasting isotopic compositions of shallow and deep plateau lithosphere suggest their derivation from different mantle reservoirs. We propose that the subplateau lithosphere forms a genetically unrelated two-layered structure, comprising shallower, typical oceanic lithosphere underpinned by deeper impinged material, which included a component of recycled Proterozoic lithosphere. The impingement of residual but chemically heterogeneous mantle, mechanically coupled to the recently formed, thin lithosphere, may have a bearing on the anomalous initial uplift and late subsidence history of the seismically anomalous plateau root.     

Tracking the lithium isotopic evolution of the mantle using carbonatites

Carbonatites are mantle-derived, intraplate magmas that provide a means of documenting isotopic variations of the Earth's mantle through time. To investigate the secular Li isotopic evolution of the mantle and to test whether Li isotopes document systematic recycling of material processed at or near the Earth's surface into the mantle, we analyzed the Li isotopic compositions of carbonatites and spatially associated mafic silicate rocks. The Li isotopic compositions of Archean (2.7 Ga) to Recent carbonatites (δ⁷Li = 4.1 ± 1.3 (n = 23, 1σ)) overlap the range typical for modern mantle-derived rocks, and do not change with time, despite ongoing crustal recycling. Thus, the average Li isotopic composition of recycled crustal components has not deviated greatly from the mantle value (~ + 4) and/or Li diffusion is sufficiently fast to attenuate significant heterogeneities over timescales of 10⁸ years. Modeling of Li diffusion at mantle temperatures suggests that limited δ⁷Li variation in the mantle through time reflects the more effective homogenization of Li in the mantle compared to radiogenic isotope systems. The real (but limited) variations in δ⁷Li that exist in modern mantle-derived magmas as well as carbonatites studied here may reflect isotopic fractionation associated with shallow-level processes, such as crustal assimilation and diffusive isotopic fractionation in magmatic systems, with some of the scatter possibly related to low-temperature alteration.     

Development of continental lithospheric mantle as reflected in the chemistry of the Mesozoic Appalachian Tholeiites, U.S.A.

Geochemical analyses of dikes, sills, and volcanic rocks of the Mesozoic Appalachian Tholeiite (MAT) Province of the easternmost United States provide evidence that continental tholeiites are derived from continental lithospheric mantle sources that are genetically and geochronologically related to the overlying continental crust. Nineteen olivine tholeiites and sixteen quartz tholeiites from the length of this province, associated in space and time with the last opening of the Atlantic, display significant isotopic heterogeneity: initial ε_Nd = +3.8 to −5.7; initial ⁸⁷Sr/⁸⁶Sr= 0.7044−0.7072; ²⁰⁶Pb/²⁰⁴Pb= 17.49−19.14; ²⁰⁷Pb/²⁰⁴Pb= 15.55−15.65; ²⁰⁸Pb/²⁰⁴Pb= 37.24−39.11. In PbPb space, the MAT define a linear array displaced above the field for MORB and thus resemble oceanic basalts with DUPAL Pb isotopic traits. A regression of this array yields a secondary PbPb isochron age of ≈ 1000 Ma (μ₁ = 8.26), similar to Sm/Nd isochrons from the southern half of the province and to the radiometric age of the Grenville crust underlying easternmost North America. The MAT exhibit significant trace element ratio heterogeneity (e.g., Sm/Nd= 0.226−0.327) and have trace element traits similar to convergent margin magmas [e.g., depletions of Nb and Ti relative to the rare earth elements on normalized trace element incompatibility diagrams, Ba/Nb ratios (19–75) that are significantly greater than those of MORB, and low TiO₂ (0.39–0.69%)].Geochemical and geological considerations very strongly suggest that the MAT were not significantly contaminated during ascent through the continental crust. Further, isotope and trace element variations are not consistent with the involvement of contemporaneous MORB or OIB components. Rather, the materials that control the MAT incompatible element chemistry were derived from subcontinental lithospheric mantle. Thus: (1) the MAT/arc magma trace element similarities; (2) the PbPb and Sm/Nd isochron ages; and (3) the need for a method of introducing an ancient (> 2−3 Ga) Pb component into subcontinental mantle that cannot be much older than 1 Ga leads to a model whereby the MAT were generated by the melting of sediment-contaminated arc mantle that was incorporated into the continental lithosphere during arc activity preceding the Grenville Orogeny (≈ 1000 Ma).     

The Th/La and Sm/La conundrum of the Tethyan realm lamproites

The Oligocene–Pleistocene Tethyan Realm Lamproites, from the Mediterranean to Himalayas, have a clear subduction-related signature and represent one of the most exotic and rare ultrapotassic mantle-derived magmas. They share the major element and mineralogical characteristics of the lamproite clan rocks, but clearly define a distinct subgroup with respect to within-plate lamproites on the basis of a number of key trace element ratios and radiogenic isotopes. The most striking characteristic of the Tethyan Realm Lamproites is the positive correlation between Th/La (up to > 1.5) and Sm/La (up to > 0.3), which is opposite to what observed in subduction-related magmas worldwide and cannot be reconciled with typical slab recycling processes. The geochemical conundrum of the Tethyan Realm Lamproites requires a component with high Sm/La and Th/La (hence named SALATHO), in addition to a normal K₂O–rich sediment melt component (with low Th/La and Sm/La). The Pb isotope composition of the Tethyan Realm Lamproites also displays a mixing array from a normal sediment melt component to a high ²⁰⁸Pb/²⁰⁶Pb and low ²⁰⁶Pb/²⁰⁴Pb end-member that can be reconciled with the SALATHO component. This end-member requires a history of high time-integrated κ? (²³²Th/²³⁸U) and low time-integrated μ (²³⁸U/²⁰⁴Pb) and ought to be older than some hundreds of million years.We propose a multi-stage process for the formation of the mantle sources of the Tethyan Realm Lamproites related to the tectonic mélange domains (i.e. chaotic mixture of depleted peridotite, basalt, and sediment) accreted to the Eurasia plate during the collisional events of the northward drifting continental slivers from Gondwana. In a first stage, the mélange domains experienced high P and low T metamorphism with segregation and stabilisation of lawsonite and zoisite/epidote veins, which potentially match the geochemical characteristics required by the SALATHO component. Successively, the subduction of the Neotethys and Alpine Tethys oceanic plates produced normal K₂O-rich sediment melts that migrated through the mantle wedge and metasomatised the depleted lithospheric mantle blocks within the chaotic mélange, forming a clinopyroxene–phlogopite vein network. Eventually, the Tertiary orogenic belt collapses triggered the onset of low-degree melting of the low-solidus fractions within the mélange domains, producing lamproitic magmas.     

Primitive neon and helium isotopic compositions of high-MgO basalts from the Kerguelen Archipelago,Indian Ocean

The geochemical characteristics of mildly alkalic basalts (24–25 Ma) erupted in the southeastern Kerguelen Archipelago are considered to represent the best estimate for the composition of the enriched Kerguelen plume end-member. A recent study of picrites and high-MgO basalts from this part of the archipelago highlighted the Pb and Hf isotopic variations and suggested the presence of mantle heterogeneities within the Kerguelen plume itself. We present new helium and neon isotopic compositions for olivines from these picrites and high-MgO basalts (6–17 wt.% MgO) both to constrain the enriched composition of the Kerguelen plume and to determine the origin of isotopic heterogeneities involved in the genesis of Kerguelen plume-related basalts. The olivine phenocrysts have extremely variable ⁴He / ³He compositions between MORB and primitive values observed in OIB (∼90,000 to 40,000; i.e., R / R_a ∼8 to 18) and they show primitive neon isotopic ratios (average ²¹Ne / ²¹Ne_ext ∼0.044). The neon isotopic systematics and the ⁴He / ³He ratios that are lower than MORB values for the Kerguelen basalts clearly suggest that the Kerguelen hotspot belongs to the family of primitive hotspots, such as Iceland and Hawaii. The rare gas signature for the Kerguelen samples, intermediate between MORB and solar, is apparently inconsistent with mixing of a primitive component with a MORB-like source, but may result from sampling a heterogeneous part of the mantle with solar ³He / ²²Ne and with a higher (U, Th) / ³He ratio compared to typically high R / R_a hotspot basalts such as those from Iceland and Hawaii.     

Os isotope heterogeneity of the upper mantle: Evidence from the Mayarí–Baracoa ophiolite belt in eastern Cuba

In an attempt to quantify the extent of geochemical heterogeneity within a restricted and well dated portion of the upper mantle, 27 chromite separates from the 90 My old chromite deposits in the Mayarí–Baracoa ophiolite belt in eastern Cuba have been investigated for platinum group element (PGE) concentrations and Re–Os isotopic systematics. The samples are characterized by systematically subchondritic initial ¹⁸⁷Os/¹⁸⁸Os ratios and substantial heterogeneity. The initial ¹⁸⁷Os/¹⁸⁸Os ratios vary with chromite chemistry and with geographical distribution, reflecting differences in the Os isotopic evolution for the different upper mantle sections represented by the ophiolite. Accordingly, the Os isotope data might be divided into three groups. In the Moa–Baracoa district, where the chromite bodies are located in the mantle–crust transition zone, the calculated initial γOs values average − 0.97 ± 0.69 (n = 13). In the Sagua de Tanamo district, where chromite chemistry is highly variable and their location in relation the mantle sequence is less clear, the initial γOs values are intermediate, with an average of − 1.77 ± 0.80 (n = 7). In the Mayarí district, where the chromite bodies are located in the lower part of the mantle sequence, initial γOs values average − 2.66 ± 0.29 (n = 7). These subchondritic (i.e. negative) initial γOs values are most simply explained by Re depletion during ancient partial melting and/or melt percolation events.The Os isotope heterogeneity documented here indicates a high degree of geochemical complexity on small to intermediate length scales in the upper mantle. Our results, in combination with data on chromites from the literature, show that an “average present-day Os isotopic composition” for the hypothetical depleted MORB mantle (DMM) reservoir cannot be precisely established beyond the statement that it is “broadly chondritic”. Indeed, the upper mantle cannot be considered a sufficiently homogeneous geochemical “reservoir” to serve meaningfully as a baseline against which geochemical “anomalies” are evaluated. On the other hand, our findings are consistent with the “Statistical Upper Mantle Assemblage” or “SUMA”-concept, according to which a high level of geochemical heterogeneity is maintained in the upper mantle at all relevant length scales, as a result of the plate-tectonic cycle and intra-mantle processes such as melt-migration and metasomatism.