Metal saturation in the upper mantle

The oxygen fugacity f(O2)of the Earth's mantle is one of the fundamental variables in mantle petrology. Through ferric-ferrous iron and carbon-hydrogen-oxygen equilibria, f(O2) influences the pressure-temperature positions of mantle solidi and compositions of small-degree mantle melts. Among other parameters, f(O2) affects the water storage capacity and rheology of the mantle. The uppermost mantle, as represented by samples and partial melts, is sufficiently oxidized to sustain volatiles, such as H2O and CO2, as well as carbonatitic melts, but it is not known whether the shallow mantle is representative of the entire upper mantle. Using high-pressure experiments, we show here that large parts of the asthenosphere are likely to be metal-saturated. We found that pyroxene and garnet synthesized at >7 GPa in equilibrium with metallic Fe can incorporate sufficient ferric iron that the mantle at >250 km depth is so reduced that an (Fe,Ni)-metal phase may be stable. Our results indicate that the oxidized nature of the upper mantle can no longer be regarded as being representative for the Earth's upper mantle as a whole and instead that oxidation is a shallow phenomenon restricted to an upper veneer only about 250 km in thickness.     

Anisotropy of thermal diffusivity in the upper mantle.

Heat transfer in the mantle is a key process controlling the Earth's dynamics. Upper-mantle mineral phases, especially olivine, have been shown to display highly anisotropic thermal diffusivity at ambient conditions, and seismic anisotropy data show that preferred orientations of olivine induced by deformation are coherent at large scales (>50 km) in the upper mantle. Thus heat transport in the upper mantle should be anisotropic. But the thermal anisotropy of mantle minerals at high temperature and its relationship with deformation have not been well constrained. Here we present petrophysical modelling and laboratory measurements of thermal diffusivity in deformed mantle rocks between temperatures of 290 and 1,250 K that demonstrate that deformation may induce a significant anisotropy of thermal diffusivity in the uppermost mantle. We found that heat transport parallel to the flow direction is up to 30 per cent faster than that normal to the flow plane. Such a strain-induced thermal anisotropy implies that the upper-mantle temperature distribution, rheology and, consequently, its dynamics, will depend on deformation history. In oceans, resistive drag flow would result in lower vertical diffusivities in both the lithosphere and asthenosphere and hence in less effective heat transfer from the convective mantle. In continents, olivine orientations frozen in the lithosphere may induce anisotropic heating above mantle plumes, favouring the reactivation of pre-existing structures.     

Continental mantle signature of Bushveld magmas and coeval diamonds

The emplacement of the 2.05-billion-year-old Bushveld complex, the world's largest layered intrusion and platinum-group element (PGE) repository, is a singular event in the history of the Kaapvaal craton of southern Africa, one of Earth's earliest surviving continental nuclei. In the prevailing model for the complex's mineralization, the radiogenic strontium and osmium isotope signatures of Bushveld PGE ores are attributed to continental crustal contamination of the host magmas. The scale of the intrusion and lateral homogeneity of the PGE-enriched layers, however, have long been problematical for the crustal contamination model, given the typically heterogeneous nature of continental crust. Furthermore, the distribution of Bushveld magmatism matches that of seismically anomalous underlying mantle, implying significant interaction before emplacement in the crust. Mineral samples of the ancient 200-km-deep craton keel, encapsulated in macrodiamonds and entrained by proximal kimberlites, reveal the nature of continental mantle potentially incorporated by Bushveld magmas. Here we show that sulphide inclusions in approximately 2-billion-year-old diamonds from the 0.5-billion-year-old Venetia and 1.2-billion-year-old Premier kimberlites (on opposite sides of the complex) have initial osmium isotope ratios even more radiogenic than those of Bushveld sulphide ore minerals. Sulphide Re-Os and silicate Sm-Nd and Rb-Sr isotope compositions indicate that continental mantle harzburgite and eclogite components, in addition to the original convecting mantle magma, most probably contributed to the genesis of both the diamonds and the Bushveld complex. Coeval diamonds provide key evidence that the main source of Bushveld PGEs is the mantle rather than the crust.     

Resistance to mantle flow inferred from the electromagnetic strike of the Australian upper mantle

Seismic anisotropy is thought to result from the strain-induced lattice-preferred orientation of mantle minerals, especially olivine, owing to shear waves propagating faster along the a-axis of olivine crystals than along the other axes. This anisotropy results in birefringence, or 'shear-wave splitting', which has been investigated in numerous studies. Although olivine is also anisotropic with respect to electrical conductivity (with the a-axis being most conductive), few studies of the electrical anisotropy of the upper mantle have been undertaken, and these have been limited to relatively shallow depths in the lithospheric upper mantle. Theoretical models of mantle flow have been used to infer that, for progressive simple shear imparted by the motion of an overriding tectonic plate, the a-axes of olivine crystals should align themselves parallel to the direction of plate motion. Here, however, we show that a significant discrepancy exists between the electromagnetic strike of the mantle below Australia and the direction of present-day absolute plate motion. We infer from this discrepancy that the a-axes of olivine crystals are not aligned with the direction of the present-day plate motion of Australia, indicating resistance to deformation of the mantle by plate motion.     

Isotopic portrayal of the Earth's upper mantle flow field

It is now well established that oceanic plates sink into the lower mantle at subduction zones, but the reverse process of replacing lost upper-mantle material is not well constrained. Even whether the return flow is strongly localized as narrow upwellings or more broadly distributed remains uncertain. Here we show that the distribution of long-lived radiogenic isotopes along the world's mid-ocean ridges can be used to map geochemical domains, which reflect contrasting refilling modes of the upper mantle. New hafnium isotopic data along the Southwest Indian Ridge delineate a sharp transition between an Indian province with a strong lower-mantle isotopic flavour and a South Atlantic province contaminated by advection of upper-mantle material beneath the lithospheric roots of the Archaean African craton. The upper mantle of both domains appears to be refilled through the seismically defined anomaly underlying South Africa and the Afar plume. Because of the viscous drag exerted by the continental keels, refilling of the upper mantle in the Atlantic and Indian domains appears to be slow and confined to localized upwellings. By contrast, in the unencumbered Pacific domain, upwellings seem comparatively much wider and more rapid.     

The depth distribution of azimuthal anisotropy in the continental upper mantle

The most likely cause of seismic anisotropy in the Earth's upper mantle is the lattice preferred orientation of anisotropic minerals such as olivine. Its presence reflects dynamic processes related to formation of the lithosphere as well as to present-day tectonic motions. A powerful tool for detecting and characterizing upper-mantle anisotropy is the analysis of shear-wave splitting measurements. Because of the poor vertical resolution afforded by this type of data, however, it has remained controversial whether the splitting has a lithospheric origin that is 'frozen-in' at the time of formation of the craton, or whether the anisotropy originates primarily in the asthenosphere, and is induced by shear owing to present-day absolute plate motions. In addition, predictions from surface-wave-derived models are largely incompatible with shear-wave splitting observations. Here we show that this disagreement can be resolved by simultaneously inverting surface waveforms and shear-wave splitting data. We present evidence for the presence of two layers of anisotropy with different fast-axis orientations in the cratonic part of the North American upper mantle. At asthenospheric depths (200-400 km) the fast axis is sub-parallel to the absolute plate motion, confirming the presence of shear related to current tectonic processes, whereas in the lithosphere (80-200 km), the orientation is significantly more northerly. In the western, tectonically active, part of North America, the fast-axis direction is consistent with the absolute plate motion throughout the depth range considered, in agreement with a much thinner lithosphere.     

Melting in the Earth's deep upper mantle caused by carbon dioxide

The onset of partial melting beneath mid-ocean ridges governs the cycling of highly incompatible elements from the mantle to the crust, the flux of key volatiles (such as CO2, He and Ar) and the rheological properties of the upper mantle. Geophysical observations indicate that melting beneath ridges begins at depths approaching 300 km, but the cause of this melting has remained unclear. Here we determine the solidus of carbonated peridotite from 3 to 10 GPa and demonstrate that melting beneath ridges may occur at depths up to 330 km, producing 0.03-0.3% carbonatite liquid. We argue that these melts promote recrystallization and realignment of the mineral matrix, which may explain the geophysical observations. Extraction of incipient carbonatite melts from deep within the oceanic mantle produces an abundant source of metasomatic fluids and a vast mantle residue depleted in highly incompatible elements and fractionated in key parent-daughter elements. We infer that carbon, helium, argon and highly incompatible heat-producing elements (such as uranium, thorium and potassium) are efficiently scavenged from depths of approximately 200-330 km in the upper mantle.     

Stability of hydrous melt at the base of the Earth's upper mantle

Seismological observations have revealed the existence of low-velocity and high-attenuation zones above the discontinuity at 410 km depth, at the base of the Earth's upper mantle. It has been suggested that a small amount of melt could be responsible for such anomalies. The density of silicate melt under dry conditions has been measured at high pressure and found to be denser than the surrounding solid, thereby allowing the melt to remain at depth. But no experimental investigation of the density of hydrous melt has yet been carried out. Here we present data constraining the density of hydrous basaltic melt under pressure to examine the stability of melt above the 410-km discontinuity. We infer that hydrous magma formed by partial melting above the 410-km discontinuity may indeed be gravitationally stable, thereby supporting the idea that low-velocity or high-attentuation regions just above the mantle transition zone may result from the presence of melt.     

Vapour undersaturation in primitive mid-ocean-ridge basalt and the volatile content of Earth's upper mantle

The analysis of volatiles in magmatic systems can be used to constrain the volatile content of the Earth's mantle and the influence that magmatic degassing has on the chemistry of the oceans and the atmosphere. But most volatile elements have very low solubilities in magmas at atmospheric pressure, and therefore virtually all erupted lavas are degassed and do not retain their primary volatile signatures. Here we report the undersaturated pre-eruptive volatile content for a suite of mid-ocean-ridge basalts from the Siqueiros intra-transform spreading centre. The undersaturation leads to correlations between volatiles and refractory trace elements that provide new constraints on volatile abundances and their behaviour in the upper mantle. Our data generate improved limits on the abundances of carbon dioxide, water, fluorine, sulphur and chlorine in the source of normal mid-ocean-ridge basalt. The incompatible behaviour of carbon dioxide, together with the CO(2)/Nb and CO(2)/Cl ratios, permit estimates of primitive carbon dioxide and chlorine to be made for degassed and chlorine-contaminated mid-ocean-ridge basalt magmas, and hence constrain degassing and contamination histories of mid-ocean ridges.     

Density of hydrous silicate melt at the conditions of Earth's deep upper mantle

The chemical evolution of the Earth and the terrestrial planets is largely controlled by the density of silicate melts. If melt density is higher than that of the surrounding solid, incompatible elements dissolved in the melt will be sequestered in the deep mantle. Previous studies on dry (water-free) melts showed that the density of silicate melts can be higher than that of surrounding solids under deep mantle conditions. However, melts formed under deep mantle conditions are also likely to contain some water, which will reduce the melt density. Here we present data constraining the density of hydrous silicate melt at the conditions of approximately 410 km depth. We show that the water in the silicate melt is more compressible than the other components, and therefore the effect of water in reducing melt density is markedly diminished under high-pressure conditions. Our study indicates that there is a range of conditions under which a (hydrous) melt could be trapped at the 410-km boundary and hence incompatible elements could be sequestered in the deep mantle, although these conditions are sensitive to melt composition as well as the composition of the surrounding mantle.     

Re-Os isotopic evidence for long-lived heterogeneity and equilibration processes in the Earth's upper mantle

The geochemical composition of the Earth's upper mantle is thought to reflect 4.5 billion years of melt extraction, as well as the recycling of crustal materials. The fractionation of rhenium and osmium during partial melting in the upper mantle makes the Re-Os isotopic system well suited for tracing the extraction of melt and recycling of the resulting mid-ocean-ridge basalt. Here we report osmium isotope compositions of more than 700 osmium-rich platinum-group element alloys derived from the upper mantle. The osmium isotopic data form a wide, essentially gaussian distribution, demonstrating that, with respect to Re-Os isotope systematics, the upper mantle is extremely heterogeneous. As depleted and enriched domains can apparently remain unequilibrated on a timescale of billions of years, effective equilibration seems to require high degrees of partial melting, such as occur under mid-ocean ridges or in back-arc settings, where percolating melts enhance the mobility of both osmium and rhenium. We infer that the gaussian shape of the osmium isotope distribution is the signature of a random mixing process between depleted and enriched domains, resulting from a 'plum pudding' distribution in the upper mantle, rather than from individual melt depletion events.     

Cryptic striations in the upper mantle revealed by hafnium isotopes in southeast Indian ridge basalts

The Earth's mantle is isotopically heterogeneous on length scales ranging from centimetres to more than 10(4) kilometres. This heterogeneity originates from partial melt extraction and plate tectonic recycling, whereas stirring during mantle convection tends to reduce it. Here we show that mid-ocean ridge basalts from 2,000 km along the southeast Indian ridge (SEIR) display a bimodal hafnium isotopic distribution. This bimodality reveals the presence of ancient compositional striations (streaks) in the Indian Ocean upper mantle. The number density of the streaks is described by a Poisson distribution, with an average thickness of approximately 40 km. Such a distribution is anticipated for a well-stirred upper mantle, in which heterogeneity is continually introduced by plate tectonic recycling, and redistributed by viscous stretching and convective refolding.     

Pressure sensitivity of olivine slip systems and seismic anisotropy of Earth's upper mantle

The mineral olivine dominates the composition of the Earth's upper mantle and hence controls its mechanical behaviour and seismic anisotropy. Experiments at high temperature and moderate pressure, and extensive data on naturally deformed mantle rocks, have led to the conclusion that olivine at upper-mantle conditions deforms essentially by dislocation creep with dominant [100] slip. The resulting crystal preferred orientation has been used extensively to explain the strong seismic anisotropy observed down to 250 km depth. The rapid decrease of anisotropy below this depth has been interpreted as marking the transition from dislocation to diffusion creep in the upper mantle. But new high-pressure experiments suggest that dislocation creep also dominates in the lower part of the upper mantle, but with a different slip direction. Here we show that this high-pressure dislocation creep produces crystal preferred orientations resulting in extremely low seismic anisotropy, consistent with seismological observations below 250 km depth. These results raise new questions about the mechanical state of the lower part of the upper mantle and its coupling with layers both above and below.     

Contrasting origins of the upper mantle revealed by hafnium and lead isotopes from the Southeast Indian Ridge

The origin of the isotopic signature of Indian mid-ocean ridge basalts has remained enigmatic, because the geochemical composition of these basalts is consistent either with pollution from recycled, ancient altered oceanic crust and sediments, or with ancient continental crust or lithosphere. The radiogenic isotopic signature may therefore be the result of contamination of the upper mantle by plumes containing recycled altered ancient oceanic crust and sediments, detachment and dispersal of continental material into the shallow mantle during rifting and breakup of Gondwana, or contamination of the upper mantle by ancient subduction processes. The identification of a process operating on a scale large enough to affect major portions of the Indian mid-ocean ridge basalt source region has been a long-standing problem. Here we present hafnium and lead isotope data from across the Indian-Pacific mantle boundary at the Australian-Antarctic discordance region of the Southeast Indian Ridge, which demonstrate that the Pacific and Indian upper mantle basalt source domains were each affected by different mechanisms. We infer that the Indian upper-mantle isotope signature in this region is affected mainly by lower continental crust entrained during Gondwana rifting, whereas the isotope signature of the Pacific upper mantle is influenced predominantly by ocean floor subduction-related processes.     

Lateral heterogeneity in the Earth's mantle

The characteristics of seismic waves detected at the Large Aperture Seismic Array in Montana suggest that, beneath the Bonin Island Arc in particular, there are lateral inhomogeneities in the mantle.     

The essential signature of a massive starburst in a distant quasar

Observations of carbon monoxide emission in high-redshift (zeta > 2) galaxies indicate the presence of large amounts of molecular gas. Many of these galaxies contain an active galactic nucleus powered by accretion of gas onto a supermassive black hole, and a key question is whether their extremely high infrared luminosities result from the active galactic nucleus, from bursts of massive star formation (associated with the molecular gas), or both. In the Milky Way, high-mass stars form in the dense cores of interstellar molecular clouds, where gas densities are n(H2) > 10(5) cm(-3) (refs 1, 2). Recent surveys show that virtually all galactic sites of high-mass star formation have similarly high densities. The bulk of the cloud material traced by CO observations, however, is at a much lower density. For galaxies in the local Universe, the HCN molecule is an effective tracer of high-density molecular gas. Here we report observations of HCN emission from the infrared-luminous 'Cloverleaf' quasar (at a redshift zeta = 2.5579). The HCN line luminosity indicates the presence of 10 billion solar masses of very dense gas, an essential feature of an immense starburst, which contributes, together with the active galactic nucleus it harbours, to its high infrared luminosity.