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.     

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.     

The Earth's mantle

Seismological images of the Earth's mantle reveal three distinct changes in velocity structure, at depths of 410, 660 and 2,700 km. The first two are best explained by mineral phase transformations, whereas the third-the D" layer-probably reflects a change in chemical composition and thermal structure. Tomographic images of cold slabs in the lower mantle, the displacements of the 410-km and 660-km discontinuities around subduction zones, and the occurrence of small-scale heterogeneities in the lower mantle all indicate that subducted material penetrates the deep mantle, implying whole-mantle convection. In contrast, geochemical analyses of the basaltic products of mantle melting are frequently used to infer that mantle convection is layered, with the deeper mantle largely isolated from the upper mantle. We show that geochemical, seismological and heat-flow data are all consistent with whole-mantle convection provided that the observed heterogeneities are remnants of recycled oceanic and continental crust that make up about 16 and 0.3 per cent, respectively, of mantle volume.     

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.     

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.     

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.     

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.     

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.     

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.     

Evolution of helium isotopes in the Earth's mantle

Degassing of the Earth's mantle through magmatism results in the irreversible loss of helium to space, and high (3)He/(4)He ratios observed in oceanic basalts have been considered the main evidence for a 'primordial' undegassed deep mantle reservoir. Here we present a new global data compilation of ocean island basalts, representing upwelling 'plumes' from the deep mantle, and show that island groups with the highest primordial signal (high (3)He/(4)He ratios) have striking chemical and isotopic similarities to mid-ocean-ridge basalts. We interpret this as indicating a common history of mantle trace element depletion through magmatism. The high (3)He/(4)He in plumes may thus reflect incomplete degassing of the deep Earth during continent and ocean crust formation. We infer that differences between plumes and the upper-mantle source of ocean-ridge basalts reflect isolation of plume sources from the convecting mantle for approximately 1-2 Gyr. An undegassed, primordial reservoir in the mantle would therefore not be required, thus reconciling a long-standing contradiction in mantle dynamics.     

Solid-liquid iron partitioning in Earth's deep mantle

Melting processes in the deep mantle have important implications for the origin of the deep-derived plumes believed to feed hotspot volcanoes such as those in Hawaii. They also provide insight into how the mantle has evolved, geochemically and dynamically, since the formation of Earth. Melt production in the shallow mantle is quite well understood, but deeper melting near the core-mantle boundary remains controversial. Modelling the dynamic behaviour of deep, partially molten mantle requires knowledge of the density contrast between solid and melt fractions. Although both positive and negative melt buoyancies can produce major chemical segregation between different geochemical reservoirs, each type of buoyancy yields drastically different geodynamical models. Ascent or descent of liquids in a partially molten deep mantle should contribute to surface volcanism or production of a deep magma ocean, respectively. We investigated phase relations in a partially molten chondritic-type material under deep-mantle conditions. Here we show that the iron partition coefficient between aluminium-bearing (Mg,Fe)SiO(3) perovskite and liquid is between 0.45 and 0.6, so iron is not as incompatible with deep-mantle minerals as has been reported previously. Calculated solid and melt density contrasts suggest that melt generated at the core-mantle boundary should be buoyant, and hence should segregate upwards. In the framework of the magma oceans induced by large meteoritic impacts on early Earth, our results imply that the magma crystallization should push the liquids towards the surface and form a deep solid residue depleted in incompatible elements.     

Grain boundaries as reservoirs of incompatible elements in the Earth's mantle

The concentrations and locations of elements that strongly partition into the fluid phase in rocks provide essential constraints on geochemical and geodynamical processes in Earth's interior. A fundamental question remains, however, as to where these incompatible elements reside before formation of the fluid phase. Here we show that partitioning of calcium between the grain interiors and grain boundaries of olivine in natural and synthetic olivine-rich aggregates follows a thermodynamic model for equilibrium grain-boundary segregation. The model predicts that grain boundaries can be the primary storage sites for elements with large ionic radius--that is, incompatible elements in the Earth's mantle. This observation provides a mechanism for the selective extraction of these elements and gives a framework for interpreting geochemical signatures in mantle rocks.     

Effect of phase transitions on compressional-wave velocities in the Earth's mantle

The velocities of seismic waves in the Earth are governed by the response of the constituent mineral assemblage to perturbations in pressure and stress. The effective bulk modulus is significantly lowered if the pressure of the seismic wave drives a volume-reducing phase transformation. A comparison between the amount of time required by phase transitions to reach equilibrium and the sampling period thus becomes crucial in defining the softening and attenuation of compressional waves within such a two-phase zone. These phenomena are difficult to assess experimentally, however, because data at conditions appropriate to the Earth's deep interior are required. Here we present synchrotron-based experimental data that demonstrate softening of the bulk modulus within the two-phase loop of olivine-ringwoodite on a timescale of 100 s. If the amplitude of the pressure perturbation and the grain size are scaled to those expected in the Earth, the compressional-wave velocities within the discontinuities at 410, 520 and, possibly, 660 km are likely to be significantly lower than otherwise expected. The generalization of these observations to aluminium-controlled phase transitions raises the possibility of large velocity perturbations throughout the upper 1,000 km of the mantle.     

Melting of the Earth's lithospheric mantle inferred from protactinium-thorium-uranium isotopic data

The processes responsible for the generation of partial melt in the Earth's lithospheric mantle and the movement of this melt to the Earth's surface remain enigmatic, owing to the perceived difficulties in generating large-degree partial melts at depth and in transporting small-degree melts through a static lithosphere. Here we present a method of placing constraints on melting in the lithospheric mantle using 231Pa-235U data obtained from continental basalts in the southwestern United States and Mexico. Combined with 230Th-238U data, the 231Pa-235U data allow us to constrain the source mineralogy and thus the depth of melting of these basalts. Our analysis indicates that it is possible to transport small melt fractions--of the order of 0.1%--through the lithosphere, as might result from the coalescence of melt by compaction owing to melting-induced deformation. The large observed 231Pa excesses require that the timescale of melt generation and transport within the lithosphere is small compared to the half-life of 231Pa (approximately 32.7 kyr). The 231Pa-230Th data also constrain the thorium and uranium distribution coefficients for clinopyroxene in the source regions of these basalts to be within 2% of one another, indicating that in this setting 230Th excesses are not expected during melting at depths shallower than 85 km.     

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.     

The tungsten isotopic composition of the Earth's mantle before the terminal bombardment

Many precious, 'iron-loving' metals, such as gold, are surprisingly abundant in the accessible parts of the Earth, given the efficiency with which core formation should have removed them to the planet's deep interior. One explanation of their over-abundance is a 'late veneer'--a flux of meteorites added to the Earth after core formation as a 'terminal' bombardment that culminated in the cratering of the Moon. Some 3.8 billion-year-old rocks from Isua, Greenland, are derived from sources that retain an isotopic memory of events pre-dating this cataclysmic meteorite shower. These Isua samples thus provide a window on the composition of the Earth before such a late veneer and allow a direct test of its importance in modifying the composition of the planet. Using high-precision (less than 6 parts per million, 2 standard deviations) tungsten isotope analyses of these rocks, here we show that they have a isotopic tungsten ratio (182)W/(184)W that is significantly higher (about 13 parts per million) than modern terrestrial samples. This finding is in good agreement with the expected influence of a late veneer. We also show that alternative interpretations, such as partial remixing of a deep-mantle reservoir formed in the Hadean eon (more than four billion years ago) or core-mantle interaction, do not explain the W isotope data well. The decrease in mantle (182)W/(184)W occurs during the Archean eon (about four to three billion years ago), potentially on the same timescale as a notable decrease in (142)Nd/(144)Nd (refs 3 and 6). We speculate that both observations can be explained if late meteorite bombardment triggered the onset of the current style of mantle convection.     

A crystallizing dense magma ocean at the base of the Earth's mantle

The distribution of geochemical species in the Earth's interior is largely controlled by fractional melting and crystallization processes that are intimately linked to the thermal state and evolution of the mantle. The existence of patches of dense partial melt at the base of the Earth's mantle, together with estimates of melting temperatures for deep mantle phases and the amount of cooling of the underlying core required to maintain a geodynamo throughout much of the Earth's history, suggest that more extensive deep melting occurred in the past. Here we show that a stable layer of dense melt formed at the base of the mantle early in the Earth's history would have undergone slow fractional crystallization, and would be an ideal candidate for an unsampled geochemical reservoir hosting a variety of incompatible species (most notably the missing budget of heat-producing elements) for an initial basal magma ocean thickness of about 1,000 km. Differences in 142Nd/144Nd ratios between chondrites and terrestrial rocks can be explained by fractional crystallization with a decay timescale of the order of 1 Gyr. These combined constraints yield thermal evolution models in which radiogenic heat production and latent heat exchange prevent early cooling of the core and possibly delay the onset of the geodynamo to 3.4-4 Gyr ago.     

Thermochemical flows couple the Earth's inner core growth to mantle heterogeneity

Seismic waves sampling the top 100 km of the Earth's inner core reveal that the eastern hemisphere (40 degrees E-180 degrees E) is seismically faster, more isotropic and more attenuating than the western hemisphere. The origin of this hemispherical dichotomy is a challenging problem for our understanding of the Earth as a system of dynamically coupled layers. Previously, laboratory experiments have established that thermal control from the lower mantle can drastically affect fluid flow in the outer core, which in turn can induce textural heterogeneity on the inner core solidification front. The resulting texture should be consistent with other expected manifestations of thermal mantle control on the geodynamo, specifically magnetic flux concentrations in the time-average palaeomagnetic field over the past 5 Myr, and preferred eddy locations in flows imaged below the core-mantle boundary by the analysis of historical geomagnetic secular variation. Here we show that a single model of thermochemical convection and dynamo action can account for all these effects by producing a large-scale, long-term outer core flow that couples the heterogeneity of the inner core with that of the lower mantle. The main feature of this thermochemical 'wind' is a cyclonic circulation below Asia, which concentrates magnetic field on the core-mantle boundary at the observed location and locally agrees with core flow images. This wind also causes anomalously high rates of light element release in the eastern hemisphere of the inner core boundary, suggesting that lateral seismic anomalies at the top of the inner core result from mantle-induced variations in its freezing rate.     

Neon isotopes constrain convection and volatile origin in the Earth's mantle

Identifying the origin of primordial volatiles in the Earth's mantle provides a critical test between models that advocate magma-ocean equilibration with an early massive solar-nebula atmosphere and those that require subduction of volatiles implanted in late accreting material. Here we show that neon isotopes in the convecting mantle, resolved in magmatic CO2 well gases, are consistent with a volatile source related to solar corpuscular irradiation of accreting material. This contrasts with recent results that indicated a solar-nebula origin for neon in mantle plume material, which is thought to be sampling the deep mantle. Neon isotope heterogeneity in different mantle sources suggests that models in which the plume source supplies the convecting mantle with its volatile inventory require revision. Although higher than accepted noble gas concentrations in the convecting mantle may reduce the need for a deep mantle volatile flux, any such flux must be dominated by the neon (and helium) isotopic signature of late accreting material.