Seismological constraints on a possible plume root at the core-mantle boundary

Recent seismological discoveries have indicated that the Earth's core-mantle boundary is far more complex than a simple boundary between the molten outer core and the silicate mantle. Instead, its structural complexities probably rival those of the Earth's crust. Some regions of the lowermost mantle have been observed to have seismic wave speed reductions of at least 10 per cent, which appear not to be global in extent. Here we present robust evidence for an 8.5-km-thick and approximately 50-km-wide pocket of dense, partially molten material at the core-mantle boundary east of Australia. Array analyses of an anomalous precursor to the reflected seismic wave ScP reveal compressional and shear-wave velocity reductions of 8 and 25 per cent, respectively, and a 10 per cent increase in density of the partially molten aggregate. Seismological data are incompatible with a basal layer composed of pure melt, and thus require a mechanism to prevent downward percolation of dense melt within the layer. This may be possible by trapping of melt by cumulus crystal growth following melt drainage from an anomalously hot overlying region of the lowermost mantle. This magmatic evolution and the resulting cumulate structure seem to be associated with overlying thermal instabilities, and thus may mark a root zone of an upwelling plume.     

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.     

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.     

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.     

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.     

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.     

Constraints on the composition of the Earth's core from ab initio calculations

Knowledge of the composition of the Earth's core is important for understanding its melting point and therefore the temperature at the inner-core boundary and the temperature profile of the core and mantle. In addition, the partitioning of light elements between solid and liquid, as the outer core freezes at the inner-core boundary, is believed to drive compositional convection, which in turn generates the Earth's magnetic field. It is generally accepted that the liquid outer core and the solid inner core consist mainly of iron. The outer core, however, is also thought to contain a significant fraction of light elements, because its density--as deduced from seismological data and other measurements--is 6-10 per cent less than that estimated for pure liquid iron. Similar evidence indicates a smaller but still appreciable fraction of light elements in the inner core. The leading candidates for the light elements present in the core are sulphur, oxygen and silicon. Here we obtain a constraint on core composition derived from ab initio calculation of the chemical potentials of light elements dissolved in solid and liquid iron. We present results for the case of sulphur, which provide strong evidence against the proposal that the outer core is close to being a binary iron-sulphur mixture.     

Ocean productivity before about 1.9 Gyr ago limited by phosphorus adsorption onto iron oxides

After the evolution of oxygen-producing cyanobacteria at some time before 2.7 billion years ago, oxygen production on Earth is thought to have depended on the availability of nutrients in the oceans, such as phosphorus (in the form of orthophosphate). In the modern oceans, a significant removal pathway for phosphorus occurs by way of its adsorption onto iron oxide deposits. Such deposits were thought to be more abundant in the past when, under low sulphate conditions, the formation of large amounts of iron oxides resulted in the deposition of banded iron formations. Under these circumstances, phosphorus removal by iron oxide adsorption could have been enhanced. Here we analyse the phosphorus and iron content of banded iron formations to show that ocean orthophosphate concentrations from 3.2 to 1.9 billion years ago (during the Archaean and early Proterozoic eras) were probably only approximately 10-25% of present-day concentrations. We suggest therefore that low phosphorus availability should have significantly reduced rates of photosynthesis and carbon burial, thereby reducing the long-term oxygen production on the early Earth--as previously speculated--and contributing to the low concentrations of atmospheric oxygen during the late Archaean and early Proterozoic.     

Mantle superplasticity and its self-made demise

The unusual capability of solid crystalline materials to deform plastically, known as superplasticity, has been found in metals and even in ceramics. Such superplastic behaviour has been speculated for decades to take place in geological materials, ranging from surface ice sheets to the Earth's lower mantle. In materials science, superplasticity is confirmed when the material deforms with large tensile strain without failure; however, no experimental studies have yet shown this characteristic in geomaterials. Here we show that polycrystalline forsterite + periclase (9:1) and forsterite + enstatite + diopside (7:2.5:0.5), which are good analogues for Earth's mantle, undergo homogeneous elongation of up to 500 per cent under subsolidus conditions. Such superplastic deformation is accompanied by strain hardening, which is well explained by the grain size sensitivity of superplasticity and grain growth under grain switching conditions (that is, grain boundary sliding); grain boundary sliding is the main deformation mechanism for superplasticity. We apply the observed strain-grain size-viscosity relationship to portions of the mantle where superplasticity has been presumed to take place, such as localized shear zones in the upper mantle and within subducting slabs penetrating into the transition zone and lower mantle after a phase transformation. Calculations show that superplastic flow in the mantle is inevitably accompanied by significant grain growth that can bring fine grained (≤1?μm) rocks to coarse-grained (1-10?mm) aggregates, resulting in increasing mantle viscosity and finally termination of superplastic flow.     

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.     

Inner-core shear-wave anisotropy and texture from an observation of PKJKP waves

Since the discovery of the Earth's core a century ago, and the subsequent discovery of a solid inner core (postulated to have formed by the freezing of iron) seismologists have striven to understand this most remote part of the deep Earth. The most direct evidence for a solid inner core would be the observation of shear-mode body waves that traverse it, but these phases are extremely difficult to observe. Two reported observations in short-period data have proved controversial. Arguably more successful have been studies of longer-period data, but such averaging limits the usefulness of the observations to reported sightings. We present two observations of an inner-core shear-wave phase at higher frequencies in stacked data from the Japanese High-Sensitivity Array, Hi-Net. From an analysis of timing, amplitude and waveform of the 'PKJKP' phase we derive constraints on inner-core compressional-wave velocity and shear attenuation at about 0.3 Hz which differ from standard isotropic core models. We can explain waveform features and can partially reconcile the otherwise large differences between core wavespeed and attenuation models that our observations apparently suggest if we invoke shear-wave anisotropy in the inner core. A simple model of an inner core composed of hexagonal close-packed iron with its c axis aligned perpendicular to the rotation axis yields anisotropy that is compatible with both the shear-wave anisotropy that we observe and the well-established 3 per cent compressional-wave anisotropy.     

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.     

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.     

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.     

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.     

Experimental evidence that potassium is a substantial radioactive heat source in planetary cores

The hypothesis that (40)K may be a significant radioactive heat source in the Earth's core was proposed on theoretical grounds over three decades ago, but experiments have provided only ambiguous and contradictory evidence for the solubility of potassium in iron-rich alloys. The existence of such radioactive heat in the core would have important implications for our understanding of the thermal evolution of the Earth and global processes such as the generation of the geomagnetic field, the core-mantle boundary heat flux and the time of formation of the inner core. Here we provide experimental evidence to show that the ambiguous results obtained from earlier experiments are probably due to previously unrecognized experimental and analytical difficulties. The high-pressure, high-temperature data presented here show conclusively that potassium enters iron sulphide melts in a strongly temperature-dependent fashion and that (40)K can serve as a substantial heat source in the cores of the Earth and Mars.     

Hemispherical variations in seismic velocity at the top of the Earth's inner core

Knowledge of the seismic velocity structure at the top of the Earth's inner core is important for deciphering the physical processes responsible for inner-core growth. Previous global seismic studies have focused on structures found 100 km or deeper within the inner core, with results for the uppermost 100 km available for only isolated regions. Here we present constraints on seismic velocity variations just beneath the inner-core boundary, determined from the difference in travel time between waves reflected at the inner-core boundary and those transmitted through the inner core. We found that these travel-time residuals-observed on both global seismograph stations and several regional seismic networks-are systematically larger, by about 0.8 s, for waves that sample the 'eastern hemisphere' of the inner core (40 degrees E to 180 degrees E) compared to those that sample the 'western hemisphere' (180 degrees W to 40 degrees E). These residuals show no correlation with the angle at which the waves traverse the inner core; this indicates that seismic anisotropy is not strong in this region and that the isotropic seismic velocity of the eastern hemisphere is about 0.8% higher than that of the western hemisphere.     

Meteoritic dust from the atmospheric disintegration of a large meteoroid

Much of the mass of most meteoroids entering the Earth's atmosphere is consumed in the process of ablation. Larger meteoroids (> 10 cm), which in some cases reach the ground as meteorites, typically have survival fractions near 1-25 per cent of their initial mass. The fate of the remaining ablated material is unclear, but theory suggests that much of it should recondense through coagulation as nanometre-sized particles. No direct measurements of such meteoric 'smoke' have hitherto been made. Here we report the disintegration of one of the largest meteoroids to have entered the Earth's atmosphere during the past decade, and show that the dominant contribution to the mass of the residual atmospheric aerosol was in the form of micrometre-sized particles. This result is contrary to the usual view that most of the material in large meteoroids is efficiently converted to particles of much smaller size through ablation. Assuming that our observations are of a typical event, we suggest that large meteoroids provide the dominant source of micrometre-sized meteoritic dust at the Earth's surface over long timescales.     

Characterization of extrasolar terrestrial planets from diurnal photometric variability

The detection of massive planets orbiting nearby stars has become almost routine, but current techniques are as yet unable to detect terrestrial planets with masses comparable to the Earth's. Future space-based observatories to detect Earth-like planets are being planned. Terrestrial planets orbiting in the habitable zones of stars-where planetary surface conditions are compatible with the presence of liquid water-are of enormous interest because they might have global environments similar to Earth's and even harbour life. The light scattered by such a planet will vary in intensity and colour as the planet rotates; the resulting light curve will contain information about the planet's surface and atmospheric properties. Here we report a model that predicts features that should be discernible in the light curve obtained by low-precision photometry. For extrasolar planets similar to Earth, we expect daily flux variations of up to hundreds of per cent, depending sensitively on ice and cloud cover as well as seasonal variations. This suggests that the meteorological variability, composition of the surface (for example, ocean versus land fraction) and rotation period of an Earth-like planet could be derived from photometric observations. Even signatures of Earth-like plant life could be constrained or possibly, with further study, even uniquely determined.