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.     

Experimental evidence for the existence of iron-rich metal in the Earth's lower mantle

The oxidation state recorded by rocks from the Earth's upper mantle can be calculated from measurements of the distribution of Fe3+ and Fe2+ between the constituent minerals. The capacity for minerals to incorporate Fe3+ may also be a significant factor controlling the oxidation state of the mantle, and high-pressure experimental measurements of this property might provide important insights into the redox state of the more inaccessible deeper mantle. Here we show experimentally that the Fe3+ content of aluminous silicate perovskite, the dominant lower-mantle mineral, is independent of oxygen fugacity. High levels of Fe3+ are present in perovskite even when it is in chemical equilibrium with metallic iron. Silicate perovskite in the lower mantle will, therefore, have an Fe3+/total Fe ratio of at least 0.6, resulting in a whole-rock ratio of over ten times that of the upper mantle. Consequently, the lower mantle must either be enriched in Fe3+ or Fe3+ must form by the disproportionation of Fe2+ to produce Fe3+ plus iron metal. We argue that the lower mantle contains approximately 1 wt% of a metallic iron-rich alloy. The mantle's oxidation state and siderophile element budget have probably been influenced by the presence of this alloy.     

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.     

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.     

Redox freezing and melting in the Earth's deep mantle resulting from carbon-iron redox coupling

Very low seismic velocity anomalies in the Earth's mantle may reflect small amounts of melt present in the peridotite matrix, and the onset of melting in the Earth's upper mantle is likely to be triggered by the presence of small amounts of carbonate. Such carbonates stem from subducted oceanic lithosphere in part buried to depths below the 660-kilometre discontinuity and remixed into the mantle. Here we demonstrate that carbonate-induced melting may occur in deeply subducted lithosphere at near-adiabatic temperatures in the Earth's transition zone and lower mantle. We show experimentally that these carbonatite melts are unstable when infiltrating ambient mantle and are reduced to immobile diamond when recycled at depths greater than ～250?kilometres, where mantle redox conditions are determined by the presence of an (Fe,Ni) metal phase. This 'redox freezing' process leads to diamond-enriched mantle domains in which the Fe(0), resulting from Fe(2+) disproportionation in perovskites and garnet, is consumed but the Fe(3+) preserved. When such carbon-enriched mantle heterogeneities become part of the upwelling mantle, diamond will inevitably react with the Fe(3+) leading to true carbonatite redox melting at ～660 and ～250 kilometres depth to form deep-seated melts in the Earth's mantle.     

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.     

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.     

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.     

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.     

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.     

Carbon solubility in olivine and the mode of carbon storage in the Earth's mantle

The total amount of carbon in the atmosphere, oceans and other near-surface reservoirs is thought to be negligible compared to that stored in the Earth's mantle. Although the mode of carbon storage in the mantle is largely unknown, observations of microbubbles on dislocations in minerals from mantle xenoliths has led to the suggestion that carbon may be soluble in silicates at high pressure. Here we report measurements of carbon solubility in olivine, the major constituent of the upper mantle, at pressures up to 3.5 GPa. We have found that, contrary to previous expectations, carbon solubility in olivine is exceedingly low--of the order of 0.1 to 1 parts per million by weight. Together with similar data for pyroxenes, garnet and spinel, we interpret this to imply that most carbon must be present as a separate phase in the deeper parts of the upper mantle, probably as a carbonate phase. Large-scale volcanic eruptions tapping such a carbonate-bearing mantle reservoir might therefore rapidly transfer large amounts of carbon dioxide into the atmosphere, consistent with models that link global mass extinctions to flood basalt eruptions via a sudden increase in atmospheric carbon dioxide levels.     

Pressure-induced changes in the compression mechanism of aluminous perovskite in the Earth's mantle

Although aluminium is the fifth most abundant element in the Earth's mantle, its effect on the physical properties of perovskite, the main mineral phase in the lower mantle, has largely been ignored. It is becoming clear, however, that many properties of MgSiO3 perovskites are remarkably sensitive to small amounts of aluminium. In particular, perovskite with only 5 wt% Al2O3 has a bulk modulus 10% lower than that of the pure magnesian end-member. The increased compressibility may be due to the high concentrations of oxygen vacancies required to balance the charge of the aluminium; if so, this would have important consequences for the mantle, as aluminous perovskites could be weaker, have lower seismic velocities and be hosts for water. To test whether oxygen vacancies exist in aluminous perovskites, I have calculated the compressibility of end-member defect-bearing perovskites using ab initio methods. The results show that perovskites with oxygen vacancies do have significantly greater compressibilities than those without such vacancies. But the results also suggest that oxygen vacancies become unfavourable at high pressures, in which case only the physical properties of the shallow lower mantle would be affected by aluminium-with the deeper mantle retaining properties similar to those of aluminium-free perovskite.     

Development of anisotropic structure in the Earth's lower mantle by solid-state convection

Seismological observations reveal highly anisotropic patches at the bottom of the Earth's lower mantle, whereas the bulk of the mantle has been observed to be largely isotropic. These patches have been interpreted to correspond to areas where subduction has taken place in the past or to areas where mantle plumes are upwelling, but the underlying cause for the anisotropy is unknown-both shape-preferred orientation of elastically heterogeneous materials and lattice-preferred orientation of a homogeneous material have been proposed. Both of these mechanisms imply that large-strain deformation occurs within the anisotropic regions, but the geodynamic implications of the mechanisms differ. Shape-preferred orientation would imply the presence of large elastic (and hence chemical) heterogeneity whereas lattice-preferred orientation requires deformation at high stresses. Here we show, on the basis of numerical modelling incorporating mineral physics of elasticity and development of lattice-preferred orientation, that slab deformation in the deep lower mantle can account for the presence of strong anisotropy in the circum-Pacific region. In this model-where development of the mineral fabric (the alignment of mineral grains) is caused solely by solid-state deformation of chemically homogeneous mantle material-anisotropy is caused by large-strain deformation at high stresses, due to the collision of subducted slabs with the core-mantle boundary.     

The elasticity of the MgSiO3 post-perovskite phase in the Earth's lowermost mantle

MgSiO3 perovskite has been assumed to be the dominant component of the Earth's lower mantle, although this phase alone cannot explain the discontinuity in seismic velocities observed 200-300 km above the core-mantle boundary (the D" discontinuity) or the polarization anisotropy observed in the lowermost mantle. Experimental and theoretical studies that have attempted to attribute these phenomena to a phase transition in the perovskite phase have tended to simply confirm the stability of the perovskite phase. However, recent in situ X-ray diffraction measurements have revealed a transition to a 'post-perovskite' phase above 125 GPa and 2,500 K--conditions close to those at the D" discontinuity. Here we show the results of first-principles calculations of the structure, stability and elasticity of both phases at zero temperature. We find that the post-perovskite phase becomes the stable phase above 98 GPa, and may be responsible for the observed seismic discontinuity and anisotropy in the lowermost mantle. Although our ground-state calculations of the unit cell do not include the effects of temperature and minor elements, they do provide a consistent explanation for a number of properties of the D" layer.     

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.     

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.     

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.