A doubling of the post-perovskite phase boundary and structure of the Earth's lowermost mantle

The thermal structure of the Earth's lowermost mantle--the D' layer spanning depths of approximately 2,600-2,900 kilometres--is key to understanding the dynamical state and history of our planet. Earth's temperature profile (the geotherm) is mostly constrained by phase transitions, such as freezing at the inner-core boundary or changes in crystal structure within the solid mantle, that are detected as discontinuities in seismic wave speed and for which the pressure and temperature conditions can be constrained by experiment and theory. A recently discovered phase transition at pressures of the D' layer is ideally situated to reveal the thermal structure of the lowermost mantle, where no phase transitions were previously known to exist. Here we show that a pair of seismic discontinuities observed in some regions of D' can be explained by the same phase transition as the result of a double-crossing of the phase boundary by the geotherm at two different depths. This simple model can also explain why a seismic discontinuity is not observed in some other regions, and provides new constraints for the magnitude of temperature variations within D'.     

Anisotropy of Earth's D' layer and stacking faults in the MgSiO3 post-perovskite phase

The post-perovskite phase of (Mg,Fe)SiO3 is believed to be the main mineral phase of the Earth's lowermost mantle (the D' layer). Its properties explain numerous geophysical observations associated with this layer-for example, the D' discontinuity, its topography and seismic anisotropy within the layer. Here we use a novel simulation technique, first-principles metadynamics, to identify a family of low-energy polytypic stacking-fault structures intermediate between the perovskite and post-perovskite phases. Metadynamics trajectories identify plane sliding involving the formation of stacking faults as the most favourable pathway for the phase transition, and as a likely mechanism for plastic deformation of perovskite and post-perovskite. In particular, the predicted slip planes are {010} for perovskite (consistent with experiment) and {110} for post-perovskite (in contrast to the previously expected {010} slip planes). Dominant slip planes define the lattice preferred orientation and elastic anisotropy of the texture. The {110} slip planes in post-perovskite require a much smaller degree of lattice preferred orientation to explain geophysical observations of shear-wave anisotropy in the D' layer.     

Theoretical and experimental evidence for a post-perovskite phase of MgSiO3 in Earth's D" layer

The Earth's lower mantle is believed to be composed mainly of (Mg,Fe)SiO3 perovskite, with lesser amounts of (Mg,Fe)O and CaSiO3 (ref. 1). But it has not been possible to explain many unusual properties of the lowermost approximately 150 km of the mantle (the D" layer) with this mineralogy. Here, using ab initio simulations and high-pressure experiments, we show that at pressures and temperatures of the D" layer, MgSiO3 transforms from perovskite into a layered CaIrO3-type post-perovskite phase. The elastic properties of the post-perovskite phase and its stability field explain several observed puzzling properties of the D" layer: its seismic anisotropy, the strongly undulating shear-wave discontinuity at its top and possibly the anticorrelation between shear and bulk sound velocities.     

Dislocation creep in MgSiO3 perovskite at conditions of the Earth's uppermost lower mantle

Seismic anisotropy provides an important observational constraint on flow in the Earth's deep interior. The quantitative interpretation of anisotropy, however, requires knowledge of the slip geometry of the constitutive minerals that are responsible for producing rock fabrics. The Earth's lower mantle is mostly composed of (Mg, Fe)SiO3 perovskite, but as MgSiO3 perovskite is not stable at high temperature under ambient pressure, it has not been possible to investigate its mechanical behaviour with conventional laboratory deformation experiments. To overcome this limitation, several attempts were made to infer the mechanical properties of MgSiO3 perovskite on the basis of analogue materials. But perovskites do not constitute an analogue series for plastic deformation, and therefore the direct investigation of MgSiO3 perovskite is necessary. Here we have taken advantage of recent advances in experimental high-pressure rheology to perform deformation experiments on coarse-grained MgSiO3 polycrystals under pressure and temperature conditions of the uppermost lower mantle. We show that X-ray peak broadening measurements developed in metallurgy can be adapted to low-symmetry minerals to identify the elementary deformation mechanisms activated under these conditions. We conclude that, under uppermost lower-mantle conditions, MgSiO3 perovskite deforms by dislocation creep and may therefore contribute to producing seismic anisotropy in rocks at such depths.     

The elastic constants of MgSiO3 perovskite at pressures and temperatures of the Earth's mantle.

The temperature anomalies in the Earth's mantle associated with thermal convection can be inferred from seismic tomography, provided that the elastic properties of mantle minerals are known as a function of temperature at mantle pressures. At present, however, such information is difficult to obtain directly through laboratory experiments. We have therefore taken advantage of recent advances in computer technology, and have performed finite-temperature ab initio molecular dynamics simulations of the elastic properties of MgSiO3 perovskite, the major mineral of the lower mantle, at relevant thermodynamic conditions. When combined with the results from tomographic images of the mantle, our results indicate that the lower mantle is either significantly anelastic or compositionally heterogeneous on large scales. We found the temperature contrast between the coldest and hottest regions of the mantle, at a given depth, to be about 800 K at 1,000 km, 1,500 K at 2,000 km, and possibly over 2,000 K at the core-mantle boundary.     

MgSiO_3后钙钛矿热容和热膨胀的理论研究

利用考虑了声子效应的准谐近似Debye模型结合从头算方法,研究了高温高压下MgSiO_3后钙钛矿的热力学性质.研究表明MgSiO_3后钙钛矿的状态方程与最新的实验结果吻合得很好.基于从头算和Debye模型,获得了压强上升至150 GPa,温度上升至2000 K时MgSiO_3后钙钛矿包括Debye温度、热容和热膨胀的热力学性质.     

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.     

MgSiO3后钙钛矿在三轴应力场中的第一原理模拟及地球物理学应用

MgSiO3后钙钛矿具有明显的弹性各向异性,可以解释地球核幔边界D″层中地震波的不连续性等现象,但目前的研究对晶格优选方位的机理一直没有好的解释.针对复杂的D″层应力环境中晶体特性与应力场和晶轴取向的密切关联,引入三轴应力场对矿物性质进行约束,模拟应力场状态下后钙钛矿相弹性性质,通过最小能量原理获得稳定的晶体空间取向类型.对具有最小能量类型的晶体弹性常数和地震波速研究发现,晶体a轴取向平行于最大压缩方向;以b轴为垂直对称轴的剪切波横向各向异性(横向各向同性平均)大于静水力场中的剪切波横向各向异性.结果支持以(010)为主滑移面的观点,为D″层各向异性和晶体取向研究提供新的数据和证据.     

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.     

Deformation of the lowermost mantle from seismic anisotropy

The lowermost part of the Earth's mantle-known as D″-shows significant seismic anisotropy, the variation of seismic wave speed with direction. This is probably due to deformation-induced alignment of MgSiO(3)-post-perovskite (ppv), which is believed to be the main mineral phase present in the region. If this is the case, then previous measurements of D″ anisotropy, which are generally made in one direction only, are insufficient to distinguish candidate mechanisms of slip in ppv because the mineral is orthorhombic. Here we measure anisotropy in D″ beneath North and Central America, where material from subducting oceanic slabs impinges on the core-mantle boundary, using shallow as well as deep earthquakes to increase the azimuthal coverage in D″. We make more than 700 individual measurements of shear wave splitting in D″ in three regions from two different azimuths in each case. We show that the previously assumed case of vertical transverse isotropy (where wave speed shows no azimuthal variation) is not possible, and that more complicated mechanisms must be involved. We test the fit of different MgSiO(3)-ppv deformation mechanisms to our results and find that shear on (001) is most consistent with observations and the expected shear above the core-mantle boundary beneath subduction zones. With new models of mantle flow, or improved experimental determination of the dominant ppv slip systems, this method will allow us to map deformation at the core-mantle boundary and link processes in D″, such as plume initiation, to the rest of the 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.     

Stability of magnesite and its high-pressure form in the lowermost mantle

Carbonates are important constituents of marine sediments and play a fundamental role in the recycling of carbon into the Earth's deep interior via subduction of oceanic crust and sediments. Study of the stability of carbonates under high pressure and temperature is thus important for modelling the carbon budget in the entire Earth system. Such studies, however, have rarely been performed under appropriate lower-mantle conditions and no experimental data exist at pressures greater than 80 GPa (refs 3-6). Here we report an in situ X-ray diffraction study of the stability of magnesite (MgCO(3)), which is the major component of subducted carbonates, at pressure and temperature conditions approaching those of the core-mantle boundary. We found that magnesite transforms to an unknown form at pressures above approximately 115 GPa and temperatures of 2,100-2,200 K (depths of approximately 2,600 km) without any dissociation, suggesting that magnesite and its high-pressure form may be the major hosts for carbon throughout most parts of the Earth's lower 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.     

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.     

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.     

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.     

Spin crossover and iron-rich silicate melt in the Earth's deep mantle

A melt has greater volume than a silicate solid of the same composition. But this difference diminishes at high pressure, and the possibility that a melt sufficiently enriched in the heavy element iron might then become more dense than solids at the pressures in the interior of the Earth (and other terrestrial bodies) has long been a source of considerable speculation. The occurrence of such dense silicate melts in the Earth's lowermost mantle would carry important consequences for its physical and chemical evolution and could provide a unifying model for explaining a variety of observed features in the core-mantle boundary region. Recent theoretical calculations combined with estimates of iron partitioning between (Mg,Fe)SiO(3) perovskite and melt at shallower mantle conditions suggest that melt is more dense than solids at pressures in the Earth's deepest mantle, consistent with analysis of shockwave experiments. Here we extend measurements of iron partitioning over the entire mantle pressure range, and find a precipitous change at pressures greater than ～76?GPa, resulting in strong iron enrichment in melts. Additional X-ray emission spectroscopy measurements on (Mg(0.95)Fe(0.05))SiO(3) glass indicate a spin collapse around 70?GPa, suggesting that the observed change in iron partitioning could be explained by a spin crossover of iron (from high-spin to low-spin) in silicate melt. These results imply that (Mg,Fe)SiO(3) liquid becomes more dense than coexisting solid at ～1,800?km depth in the lower mantle. Soon after the Earth's formation, the heat dissipated by accretion and internal differentiation could have produced a dense melt layer up to ～1,000?km in thickness underneath the solid mantle. We also infer that (Mg,Fe)SiO(3) perovskite is on the liquidus at deep mantle conditions, and predict that fractional crystallization of dense magma would have evolved towards an iron-rich and silicon-poor composition, consistent with seismic inferences of structures in the core-mantle boundary region.     

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.