Subducted banded iron formations as a source of ultralow-velocity zones at the core-mantle boundary

Ultralow-velocity zones (ULVZs) are regions of the Earth's core-mantle boundary about 1-10 kilometres thick exhibiting seismic velocities that are lower than radial-Earth reference models by about 10-20 per cent for compressional waves and 10-30 per cent for shear waves. It is also thought that such regions have an increased density of about 0-20 per cent (ref. 1). A number of origins for ULVZs have been proposed, such as ponding of dense silicate melt, core-mantle reaction zones or underside sedimentation from the core. Here we suggest that ULVZs might instead be relics of banded iron formations subducted to the core-mantle boundary between 2.8 and 1.8 billion years ago. Consisting mainly of interbedded iron oxides and silica, such banded iron formations were deposited in the world's oceans during the late Archaean and early Proterozoic eras. We argue that these layers, as part of the ocean floor, would be recycled into the Earth's interior by subduction, sink to the bottom of the mantle and may explain all of the observed features of ULVZs.     

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.     

The influence of a chemical boundary layer on the fixity,spacing and lifetime of mantle plumes

Seismological observations provide evidence that the lowermost mantle contains superposed thermal and compositional boundary layers that are laterally heterogeneous. Whereas the thermal boundary layer forms as a consequence of the heat flux from the Earth's outer core, the origin of an (intrinsically dense) chemical boundary layer remains uncertain. Observed zones of 'ultra-low' seismic velocity suggest that this dense layer may contain metals or partial melt, and thus it is reasonable to expect the dense layer to have a relatively low viscosity. Also, it is thought that instabilities in the thermal boundary layer could lead to the intermittent formation and rise of mantle plumes. Flow into ascending plumes can deform the dense layer, leading, in turn, to its gradual entrainment. Here we use analogue experiments to show that the presence of a dense layer at the bottom of the mantle induces lateral variations in temperature and viscosity that, in turn, determine the location and dynamics of mantle plumes. A dense layer causes mantle plumes to become spatially fixed, and the entrainment of low-viscosity fluid enables plumes to persist within the Earth for hundreds of millions of years.     

Efficacy of the post-perovskite phase as an explanation for lowermost-mantle seismic properties

Constraining the chemical, rheological and electromagnetic properties of the lowermost mantle (D') is important to understand the formation and dynamics of the Earth's mantle and core. To explain the origin of the variety of characteristics of this layer observed with seismology, a number of theories have been proposed, including core-mantle interaction, the presence of remnants of subducted material and that D' is the site of a mineral phase transformation. This final possibility has been rejuvenated by recent evidence for a phase change in MgSiO3 perovskite (thought to be the most prevalent phase in the lower mantle) at near core-mantle boundary temperature and pressure conditions. Here we explore the efficacy of this 'post-perovskite' phase to explain the seismic properties of the lowermost mantle through coupled ab initio and seismic modelling of perovskite and post-perovskite polymorphs of MgSiO3, performed at lowermost-mantle temperatures and pressures. We show that a post-perovskite model can explain the topography and location of the D' discontinuity, apparent differences in compressional- and shear-wave models and the observation of a deeper, weaker discontinuity. Furthermore, our calculations show that the regional variations in lower-mantle shear-wave anisotropy are consistent with the proposed phase change in MgSiO3 perovskite.     

Iron-silica interaction at extreme conditions and the electrically conducting layer at the base of Earth's mantle

The boundary between the Earth's metallic core and its silicate mantle is characterized by strong lateral heterogeneity and sharp changes in density, seismic wave velocities, electrical conductivity and chemical composition. To investigate the composition and properties of the lowermost mantle, an understanding of the chemical reactions that take place between liquid iron and the complex Mg-Fe-Si-Al-oxides of the Earth's lower mantle is first required. Here we present a study of the interaction between iron and silica (SiO2) in electrically and laser-heated diamond anvil cells. In a multianvil apparatus at pressures up to 140 GPa and temperatures over 3,800 K we simulate conditions down to the core-mantle boundary. At high temperature and pressures below 40 GPa, iron and silica react to form iron oxide and an iron-silicon alloy, with up to 5 wt% silicon. At pressures of 85-140 GPa, however, iron and SiO2 do not react and iron-silicon alloys dissociate into almost pure iron and a CsCl-structured (B2) FeSi compound. Our experiments suggest that a metallic silicon-rich B2 phase, produced at the core-mantle boundary (owing to reactions between iron and silicate), could accumulate at the boundary between the mantle and core and explain the anomalously high electrical conductivity of this region.     

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'.     

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.     

Chemical interaction of Fe and Al(2)O3 as a source of heterogeneity at the Earth's core-mantle boundary

Seismological studies have revealed that a complex texture or heterogeneity exists in the Earth's inner core and at the boundary between core and mantle. These studies highlight the importance of understanding the properties of iron when modelling the composition and dynamics of the core and the interaction of the core with the lowermost mantle. One of the main problems in inferring the composition of the lowermost mantle is our lack of knowledge of the high-pressure and high-temperature chemical reactions that occur between iron and the complex Mg-Fe-Si-Al-oxides which are thought to form the bulk of the Earth's lower mantle. A number of studies have demonstrated that iron can react with MgSiO3-perovskite at high pressures and high temperatures, and it was proposed that the chemical nature of this process involves the reduction of silicon by the more electropositive iron. Here we present a study of the interaction between iron and corundum (Al(2)O3) in electrically- and laser-heated diamond anvil cells at 2,000-2,200 K and pressures up to 70 GPa, simulating conditions in the Earth's deep interior. We found that at pressures above 60 GPa and temperatures of 2,200 K, iron and corundum react to form iron oxide and an iron-aluminium alloy. Our results demonstrate that iron is able to reduce aluminium out of oxides at core-mantle boundary conditions, which could provide an additional source of light elements in the Earth's core and produce significant heterogeneity at the core-mantle boundary.     

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.     

Seismic detection of folded, subducted lithosphere at the core-mantle boundary

Seismic tomography has been used to infer that some descending slabs of oceanic lithosphere plunge deep into the Earth's lower mantle. The fate of these slabs has remained unresolved, but it has been postulated that their ultimate destination is the lowermost few hundred kilometres of the mantle, known as the D' region. Relatively cold slab material may account for high seismic velocities imaged in D' beneath areas of long-lived plate subduction, and for reflections from a seismic velocity discontinuity just above the anomalously high wave speed regions. The D' discontinuity itself is probably the result of a phase change in relatively low-temperature magnesium silicate perovskite. Here, we present images of the D' region beneath the Cocos plate using Kirchhoff migration of horizontally polarized shear waves, and find a 100-km vertical step occurring over less than 100 km laterally in an otherwise flat D' shear velocity discontinuity. Folding and piling of a cold slab that has reached the core-mantle boundary, as observed in numerical and experimental models, can account for the step by a 100-km elevation of the post-perovskite phase boundary due to a 700 degrees C lateral temperature reduction in the folded slab. We detect localized low velocities at the edge of the slab material, which may result from upwellings caused by the slab laterally displacing a thin hot thermal boundary layer.     

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.     

Slow differential rotation of the Earth's inner core indicated by temporal changes in scattering

The finding that the Earth's inner core might be rotating faster than the mantle has important implications for our understanding of core processes, including the generation of the Earth's magnetic field. But the reported signal is subtle--a change of about 0.01 s per year in the separation of two seismic waves with differing paths through the core. Subsequent studies of such data have generally supported the conclusion that differential rotation exists, but the difficulty of accurately locating historic earthquakes and possible biases induced by strong lateral variations in structure near the core-mantle boundary have raised doubt regarding the proposed inner-core motion. Also, a study of free oscillations constrained the motion to be relatively small compared to previous estimates and it has been proposed that the interaction of inner-core boundary topography and mantle heterogeneity might lock the inner core to the mantle. The recent detection of seismic waves scattered in the inner core suggests a simple test of inner-core motion. Here we compare scattered waves recorded in Montana, USA, from two closely located nuclear tests at Novaya Zemlya, USSR, in 1971 and 1974. The data show small but coherent changes in scattering which point toward an inner-core differential rotation rate of 0.15 degrees per year--consistent with constraints imposed by the free-oscillation data.     

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.     

Low-velocity zone atop the 410-km seismic discontinuity in the northwestern United States

The seismic discontinuity at 410 km depth in the Earth's mantle is generally attributed to the phase transition of (Mg,Fe)2SiO4 (refs 1, 2) from the olivine to wadsleyite structure. Variation in the depth of this discontinuity is often taken as a proxy for mantle temperature owing to its response to thermal perturbations. For example, a cold anomaly would elevate the 410-km discontinuity, because of its positive Clapeyron slope, whereas a warm anomaly would depress the discontinuity. But trade-offs between seismic wave-speed heterogeneity and discontinuity topography often inhibit detailed analysis of these discontinuities, and structure often appears very complicated. Here we simultaneously model seismic refracted waves and scattered waves from the 410-km discontinuity in the western United States to constrain structure in the region. We find a low-velocity zone, with a shear-wave velocity drop of 5%, on top of the 410-km discontinuity beneath the northwestern United States, extending from southwestern Oregon to the northern Basin and Range province. This low-velocity zone has a thickness that varies from 20 to 90 km with rapid lateral variations. Its spatial extent coincides with both an anomalous composition of overlying volcanism and seismic 'receiver-function' observations observed above the region. We interpret the low-velocity zone as a compositional anomaly, possibly due to a dense partial-melt layer, which may be linked to prior subduction of the Farallon plate and back-arc extension. The existence of such a layer could be indicative of high water content in the Earth's transition zone.     

Electromagnetic detection of a 410-km-deep melt layer in the southwestern United States

A deep-seated melt or fluid layer on top of the 410-km-deep seismic discontinuity in Earth's upper mantle, as proposed in the transition-zone 'water filter' hypothesis, may have significant bearing on mantle dynamics and chemical differentiation. The geophysical detection of such a layer has, however, proved difficult. Magnetotelluric and geomagnetic depth sounding are geophysical methods sensitive to mantle melt. Here we use these methods to search for a distinct structure near 410-km depth. We calculate one-dimensional forward models of the response of electrical conductivity depth profiles, based on mineral physics studies of the effect of incorporating hydrogen in upper-mantle and transition-zone minerals. These models indicate that a melt layer at 410-km depth is consistent with regional magnetotelluric and geomagnetic depth sounding data from the southwestern United States (Tucson). The 410-km-deep melt layer in this model has a conductance of 3.0 x 10(4) S and an estimated thickness of 5-30 km. This is the only regional data set that we have examined for which such a melt layer structure was found, consistent with regional seismic studies. We infer that the hypothesized transition-zone water filter occurs regionally, but that such a layer is unlikely to be a global feature.     

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.     

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.     

A sharp lithosphere-asthenosphere boundary imaged beneath eastern North America

Plate tectonic theory hinges on the concept of a relatively rigid lithosphere moving over a weaker asthenosphere, yet the nature of the lithosphere-asthenosphere boundary remains poorly understood. The gradient in seismic velocity that occurs at this boundary is central to constraining the physical and chemical properties that create differences in mechanical strength between the two layers. For example, if the lithosphere is simply a thermal boundary layer that is more rigid owing to colder temperatures, mantle flow models indicate that the velocity gradient at its base would occur over tens of kilometres. In contrast, if the asthenosphere is weak owing to volatile enrichment or the presence of partial melt, the lithosphere-asthenosphere boundary could occur over a much smaller depth range. Here we use converted seismic phases in eastern North America to image a very sharp seismic velocity gradient at the base of the lithosphere-a 3-11 per cent drop in shear-wave velocity over a depth range of 11 km or less at 90-110 km depth. Such a strong, sharp boundary cannot be reconciled with a purely thermal gradient, but could be explained by an asthenosphere that contains a few per cent partial melt or that is enriched in volatiles relative to the lithosphere.     

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.