A refined solution to Earth's hidden niobium: implications for evolution of continental crust and mode of core formation

New high-precision niobium (Nb) and tantalum (Ta) concentration data are presented for early Archaean metabasalts, metabasaltic komatiites and their erosion products (mafic metapelites) from SW Greenland and the Acasta gneiss complex, Canada. Individual datasets consistently show sub-chondritic Nb/Ta ratios averaging 15.1+/-11.6. This finding is discussed with regard to two competing models for the solution of the Nb-deficit that characterises the accessible Earth. Firstly, we test whether Nb could have sequestered into the core due to its slightly siderophile (or chalcophile) character under very reducing conditions, as recently proposed from experimental evidence. We demonstrate that troilite inclusions of the Canyon Diablo iron meteorite have Nb and V concentrations in excess of typical chondrites but that the metal phase of the Grant, Toluca and Canyon Diablo iron meteorites do not have significant concentrations of these lithophile elements. We find that if the entire accessible Earth Nb-deficit were explained by Nb in the core, only ca. 17% of the mantle could be depleted and that by 3.7 Ga, continental crust would have already achieved ca. 50% of its present mass. Nb/Ta systematics of late Archaean metabasalts compiled from the literature would further require that by 2.5 Ga, 90% of the present mass of continental crust was already in existence. As an alternative to this explanation, we propose that the average Nb/Ta ratio (15.1+/-11.6) of Earth's oldest mafic rocks is a valid approximation for bulk silicate Earth. This would require that ca. 13% of the terrestrial Nb resided in the Ta-free core. Since the partitioning of Nb between silicate and metal melts depends largely on oxygen fugacity and pressure, this finding could mean that metal/silicate segregation did not occur at the base of a deep magma ocean or that the early mantle was slightly less reducing than generally assumed. A bulk silicate Earth Nb/Ta ratio of 15.1 allows for depletion of up to 40% of the total mantle. This could indicate that in addition to the upper mantle, a portion of the lower mantle is depleted also, or if only the upper mantle were depleted, an additional hidden high Nb/Ta reservoir must exist. Comparison of Nb/Ta systematics between early and late Archaean metabasalts supports the latter idea and indicates deeply subducted high Nb/Ta eclogite slabs could reside in the mantle transition zone or the lower mantle. Accumulation of such slabs appears to have commenced between 2.5 and 2.0 Ga. Regardless of these complexities of terrestrial Nb/Ta systematics, it is shown that the depleted mantle Nb/Th ratio is a very robust proxy for the amount of extracted continental crust, because the temporal evolution of this ratio is dominated by Th-loss to the continents and not Nb-retention in the mantle. We present a new parameterisation of the continental crust volume versus age curve that specifically explores the possibility of lithophile element loss to the core and storage of eclogite slabs in the transition zone. (C) 2003 Elsevier Science B.V. All rights reserved.

Petrology and geochemistry of early cretaceous bimodal continental flood volcanism of the NW Etendeka, Namibia. Part 1: Introduction, mafic lavas and re-evaluation of mantle source components

The bimodal NW Etendeka province is located at the continental end of the Tristan plume trace in coastal Namibia. It comprises a high-Ti (Khumib type) and three low-Ti basalt (Tafelberg, Kuidas and Esmeralda types) suites, with, at stratigraphically higher level, interstratified high-Ti latites (three units) and quartz latites (five units), and one low-Ti quartz latite. Khumib basalts are enriched in high field strength elements and light rare earth elements relative to low-Ti types and exhibit trace element affinities with Tristan da Cunha lavas. The unradiogenic Pb-206/Pb-204 ratios of Khumib basalts are distinctive, most plotting to the left of the 132 Ma Geochron, together with elevated Pb-207/Pb-204 ratios, and Sr-Nd isotopic compositions plotting in the lower Nd-143/Nd-144 part of mantle array (EM1-like). The low-Ti basalts have less coherent trace element patterns and variable, radiogenic initial Sr (similar to0.707-0.717) and Pb isotope compositions, implying crustal contamination. Four samples, however, have less radiogenic Pb and Sr that we suggest approximate their uncontaminated source. All basalt types, but particularly the low-Ti types, contain samples with trace element characteristics (e.g. Nb/Nb-*) suggesting metasediment input, considered source-related. Radiogenic isotope compositions of these samples require long-term isolation of the source in the mantle and depletions (relative to unmodified sediment) in certain elements (e.g. Cs, Pb, U), which are possibly subduction-related. A geodynamic model is proposed in which the emerging Tristan plume entrained subducted material in the Transition Zone region, and further entrained asthenosphere during plume head expansion. Mixing calculations suggest that the main features of the Etendeka basalt types can be explained without sub-continental lithospheric mantle input. Crustal contamination is evident in most low-Ti basalts, but is distinct from the incorporation of a metasedimentary source component at mantle depths.

Chemical characterization of earth's most ancient clastic metasediments from the Isua Greenstone Belt, southern West Greenland

New and published major and trace element abundances of elastic metasediments (mainly garnet-biotite-plagioclase schists) from the similar to 3.8 Ga Isua Greenstone Belt (IGB), southern West Greenland, are used in an attempt to identify the compositional characteristics of the protoliths of these sediments. Compositionally, the metasediments are heterogeneous with enrichment of LREE (La/Sm-chord = 1.1-3.9) and variable enrichment and depletion of HREE (Gd/Yb-chord = 0.8-4.3). Chondrite-normalized Eu is also variable, spanning a range from relative Eu depletion to enrichment (Eu/Eu* = 0.6-1.3). A series of geochemical and geological criteria provides conclusive evidence for a sedimentary origin, in disagreement with some previous studies that questioned the presence of genuine elastic metasediments. In particular, trace element systematics of IGB metasediments show strong resemblance to other well-documented Archaean clastic sediments, and are consistent with a provenance consisting of ultramafic, malic and felsic igneous rocks. Two schists, identified as metasomatized mafic igneous rocks from petrographic and field evidence, show distinct compositional differences to the metasediments. Major element systematics document incipient-to-moderate source weathering in the majority of metasediments, while signs of secondary K-addition are rare. Detailed inspection of Eu/Eu*, Fe2O3 and CIW (chemical index of weathering) relationships reveals that elevated iron contents (when compared to averages for continental crust) and strong relative enrichment in Eu may be due to precipitation of marine Fe-oxyhydroxides during deposition or diagenesis on the seafloor. Some of the IGB metasediments have yielded anomalous Nd-142 and W-182 isotopic compositions that were respectively interpreted in terms of early mantle differentiation processes and the presence of a meteorite component. Alternatively, W and possibly Nd isotopes could have been affected by thermal neutron capture on the Hadean surface. The latter process was tested in this study by analysis of Sm isotope compositions, which serve as an effective monitor for neutron capture effects. As no anomalous variation from terrestrial values was detected, we infer that isotope systematics (including W-182 and Nd-142) of IGB metasediments were not affected by neutron capture, but reflect decay of radioactive parent isotopes. Copyright (c) 2005 Elsevier Ltd.

An anisotropic viscous representation of Mohr-Coulomb failure for use is modelling coupled mantle-continent dynamics

In mantle convection models it has become common to make use of a modified (pressure sensitive, Boussinesq) von Mises yield criterion to limit the maximum stress the lithosphere can support. This approach allows the viscous, cool thermal boundary layer to deform in a relatively plate-like mode even in a fully Eulerian representation. In large-scale models with embedded continental crust where the mobile boundary layer represents the oceanic lithosphere, the von Mises yield criterion for the oceans ensures that the continents experience a realistic broad-scale stress regime. In detailed models of crustal deformation it is, however, more appropriate to choose a Mohr-Coulomb yield criterion based upon the idea that frictional slip occurs on whichever one of many randomly oriented planes happens to be favorably oriented with respect to the stress field. As coupled crust/mantle models become more sophisticated it is important to be able to use whichever failure model is appropriate to a given part of the system. We have therefore developed a way to represent Mohr-Coulomb failure within a code which is suited to mantle convection problems coupled to large-scale crustal deformation. Our approach uses an orthotropic viscous rheology (a different viscosity for pure shear to that for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation - neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane which is formally equivalent in any symmetric formulation. It is also important to allow strain-weakening of the material. The material should remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-Integration-point finite element code and will show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion, and comparisons between mantle convection models using the von Mises versus the Mohr-Coulomb yield criteria. The formulation itself is general and applies to 2D and 3D problems, although it is somewhat more complicated to identify the slip plane in 3D.

Anisotropy model for mantle convection

The paper presents a new theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including folding and kinking in multi-layered visco-elastic rock. The orientation of slip planes in the context of crystallographic slip is determined by the normal vector, the so-called director of these surfaces. The model is applied to simulate anisotropic natural mantle convection. We compare the evolution of the director and approximately steady states of isotropic and anisotropic convection. The isotropic case has a simple steady state solution, whereas the orthotropic convection model produces a continuously evolving patterning in tile core of the convection cell which makes only a near-steady condition possible, in which the thermal boundary layer appears to be well aligned with the flow and hence as observed in seismic tomomgraphy strong anistropic.

Investigation of Mantle Dynamics from Platinum Group Elements and Rhenium-Osmium Isotope Systematics of Mantle Xenoliths from Oahu, Hawaii

Intraplate volcanism that has created the Hawaiian-Emperor seamount chain is generally thought to be formed by a deep-seated mantle plume. While the idea of a Hawaiian plume has not met with substantial opposition, whether or not the Hawaiian plume shows any geochemical signal of receiving materials from the Earth’s Outer Core and how the plume may or may not be reacting with the overriding lithosphere remain debatable issues. In an effort to understand how the Hawaiian plume works I report on the first in-situ sulfides and bulk rock Platinum Group Element (PGE) concentrations, together with Os isotope ratios on well-characterized garnet pyroxenite xenoliths from the island of Oahu in Hawaii. The sulfides are Fe-Ni Monosulfide Solid Solution and show fractionated PGE patterns. Based on the major elements, Platinum Group Elements and experimental data I interpret the Hawaiian sulfides as an immiscible melt that separated from a melt similar to the Honolulu Volcanics (HV) alkali lavas at a pressure-temperature condition of 1530 ± 100OC and 3.1±0.6 GPa., i.e. near the base or slightly below the Pacific lithosphere. The 187Os/188Os ratios of the bulk rock vary from subchondritic to suprachondritic (0.123-0.164); and the 187Os/188Os ratio strongly correlates with major element, High Field Strength Element (HFSE), Rare Earth Element (REE) and PGE abundances. These correlations strongly suggest that PGE concentrations and Os isotope ratios reflect primary mantle processes. I interpret these correlations as the result of melt-mantle reaction at the base of the lithosphere: I suggest that the parental melt that crystallized the pyroxenites selectively picked up radiogenic Os from the grain boundary sulfides, while percolating through the Pacific lithosphere. Thus the sampled pyroxenites essentially represent crystallized melts from different stages of this melt-mantle reaction process at the base of the lithosphere. I further show that the relatively low Pt/Re ratios of the Hawaiian sulfides and the bulk rock pyroxenites suggest that, upon ageing, such pyroxenites plus their sulfides cannot generate the coupled 186Os- 187Os isotope enrichments observed in Hawaiian lavas. Therefore, recycling of mantle sulfides of pyroxenitic parentage is unlikely to explain the enriched Pt-Re-Os isotope systematics of plume-derived lavas.

Arc crust-magma interaction in the Andean Southern Volcanic Zone from thermobarometry, mineral composition, radiogenic isotope and rare earth element systematics of the Azufre-Planchon-Peteroa volcanic complex, Chile

The Andean Southern Volcanic Zone (SVZ) is a vast and complex continental arc that has been studied extensively to provide an understanding of arc-magma genesis, the origin and chemical evolution of the continental crust, and geochemical compositions of volcanic products. The present study focuses on distinguishing the magma/sub-arc crustal interaction of eruptive products from the Azufre-Planchon-Peteroa (APP 35°15'S) volcanic center and other major centers in the Central SVZ (CSVZ 37°S–42°S), Transitional SVZ (TSVZ 34.3–37.0°S), and Northern SVZ (NSVZ 33°S–34°30'S). New Hf and Nd isotopic and trace element data for SVZ centers are consistent with former studies that these magmas experienced variable depths of crystal fractionation, and that crustal assimilation is restricted to the lower crustal depths with an apparent role of garnet. Thermobarometric calculations applied to magma compositions constrain the depth of magma separation from mantle sources in all segments of the SVZ to(70-90 km). Magmatic separation at the APP complex occurs at an average depth of ~50 km which is confined to the mantle lithosphere and the base of the crust suggesting localized thermal abrasion both reservoirs. Thermobarometric calculations indicate that CSVZ primary magmas arise from a similar average depth of (~54 km) which confines magma separation to the asthenospheric mantle. The northwards along-arc Sr-Nd-Hf isotopic data and LREE enrichment accompanied with HREE depletion of SVZ mafic magmas correlates well with northward increasing crustal thickness and decreasing primary melt separation from mantle source regions indicating an increased involvement of lower crustal components in SVZ magma petrogenesis. ^ The study concludes that the development of mature subduction zones over millions of years of continuous magmatism requires that mafic arc derived melts stagnate at lower crustal levels due to density similarities and emplace at lower crustal depths. Basaltic underplating creates localized hot zone environments below major magmatic centers. These regions of high temperature/partial melting, and equilibration with underplated mafic rocks provides the mechanism that controls trace element and isotopic variability of primary magmas of the TSVZ and NSVZ from their baseline CSVZ-like precursors.^

Rare earth element concentrations and geochemistry of carbonate minerals from the Logatchev Hydrothermal Field and the Gakkel Ridge

The ultramafic-hosted Logatchev Hydrothermal Field (LHF) at 15°N on the Mid-Atlantic Ridge and the Arctic Gakkel Ridge (GR) feature carbonate precipitates (aragonite, calcite, and dolomite) in voids and fractures within different types of host rocks. We present chemical and Sr isotopic compositions of these different carbonates to examine the conditions that led to their formation. Our data reveal that different processes have led to the precipitation of carbonates in the various settings. Seawater-like 87Sr/86Sr ratios for aragonite in serpentinites (0.70909 to 0.70917) from the LHF are similar to those of aragonite from the GR (0.70912 to 0.70917) and indicate aragonite precipitation from seawater at ambient conditions at both sites. Aragonite veins in sulfide breccias from LHF also have seawater-like Sr isotope compositions (0.70909 to 0.70915), however, their rare earth element (REE) patterns show a clear positive europium (Eu) anomaly indicative of a small (< 1%) hydrothermal contribution. In contrast to aragonite, dolomite from the LHF has precipitated at much higher temperatures (~100 °C), and yet its 87Sr/86Sr ratios (0.70896 to 0.70907) are only slightly lower than those of aragonite. Even higher temperatures are calculated for the precipitation of deformed calcite veins in serpentine-talc fault schists form north of the LHF. These calcites show unradiogenic 87Sr/86Sr ratios (0.70460 to 0.70499) indicative of precipitation from evolved hydrothermal fluids. A simple mixing model based on Sr mass balance and enthalpy conservation indicates strongly variable conditions of fluid mixing and heat transfers involved in carbonate formation. Dolomite precipitated from a mixture of 97% seawater and 3% hydrothermal fluid that should have had a temperature of approximately 14 °C assuming that no heat was transferred. The much higher apparent precipitation temperatures based on oxygen isotopes (~ 100 °C) may be indicative of conductive heating, probably of seawater prior to mixing. The hydrothermal calcite in the fault schist has precipitated from a mixture of 67% hydrothermal fluid and 33% seawater, which should have had an isenthalpic mixing temperature of ~ 250 °C. The significantly lower temperatures calculated from oxygen isotopes are likely due to conductive cooling of hydrothermal fluid discharging along faults. Rare earth element patterns corroborate the results of the mixing model, since the hydrothermal calcite, which formed from waters with the greatest hydrothermal contribution, has REE patterns that closely resemble those of vent fluids from the LHF. Our results demonstrate, for the first time, that (1) precipitation from pure seawater, (2) conductive heating of seawater, and (3) conductive cooling of hydrothermal fluids in the sub-seafloor all can lead to carbonate precipitation within a single ultramafic-hosted hydrothermal system.

Seismic waveform inversion for core-mantle boundary topography

International audience

Savani: A variable resolution whole-mantle model of anisotropic shear velocity variations based on multiple data sets

International audience