The Earth's mantle

Evidence of low velocity layers at the top and bottom of the Mantle Transition Zone (MTZ) below the Uttarakhand Himalaya, India

The Mantle Transition Zone (MTZ) beneath the Uttarakhand Himalaya has been modelled using Common Conversion Point (CCP) stacking and depth-migration of radial P-receiver functions. In the Uttarakhand Himalaya region, the depths of the 410-km discontinuity (d410) and the 660-km discontinuity (d660) are estimated to be approximately 406 ± 8 km and 659 ± 10 km, respectively. Additionally, the thickness of the mantle transition zone (MTZ) is modelled to be 255 ± 7 km. The average arrival times for d410 and d660 conversions are (44.47 ± 1.33) s and (71.08 ± 1.29) s, respectively, indicating an undisturbed slightly deeper d410 and a deformed noticeably deeper d660 in the area. The model identifies the characteristics of the d410 and d660 mantle discontinuities beneath the Lesser Himalayan region, revealing a thickening of the MTZ towards northeast, which could be due to gradual cooling or thickening of the Indian lithosphere towards its northern limit. We simulate a low-velocity layer (perhaps partially molten) above the d410 discontinuity at depths of 350 to 385 km, indicating the presence of a hydrated MTZ beneath the area. We also interpret a negative phase at d660 as a low-velocity layer between 590 and 640 km depths, which could be attributed to the accumulation of old subducted oceanic materials or increased water content at the bottom of the MTZ. Our results suggest the presence of residues from paleo-subducted lithospheric slabs in and below the mantle transition zone underlying the Uttarakhand Himalayas.

Vertical tearing of subducting plates controlled by geometry and rheology of oceanic plates

Lateral non-uniform subduction is impacted by continuous plate segmentation owing to vertical tearing of the subducting plate. However, the dynamics and physical controls of vertical tearing remain controversial. Here, we employed 3D numerical models to investigate the effects of trench geometry (offset by a transform boundary) and plate rheology (plate age and the magnitude of brittle/plastic strain weakening) on the evolution of shear stress-controlled vertical tearing within a homogenous subducting oceanic plate. Numerical results suggest that the trench offset geometry could result in self-sustained vertical tearing as a narrow shear zone within the intact subducting oceanic plate, and that this process of tearing could operate throughout the entire subduction process. Further, the critical trench offset length for the maturation of vertical tearing is impacted by plate rheology. Comparison between numerical modelling results and natural observations suggests that vertical tearing attributed to trench offset geometry is broadly developed in modern subduction and collision systems worldwide.

Seismic evidence for global basalt accumulation in the mantle transition zone

The mantle's compositional structure reflects the thermochemical evolution of Earth. Yet, even the radial average composition of the mantle remains debated. Here, we analyze a global dataset of shear and compressional waves reflecting off the 410- and 660-km discontinuities that is 10 times larger than any previous studies. Our array analysis retrieves globally averaged amplitude-distance trends in SS and PP precursor reflectivity from which we infer relative wavespeed and density contrasts and associated mantle composition. Our results are best matched by a basalt-enriched mantle transition zone, with higher basalt fractions near 660 (~40%) than 410 (~18-31%). These are consistent with mantle-convection/plate-recycling simulations, which predict that basaltic crust accumulates in the mantle transition zone, with basalt fractions peaking near the 660. Basalt segregation in the mantle transition zone also implies that the overall mantle is more silica enriched than the often-assumed pyrolitic mantle reference composition.

Seismic evidence for a 1000 km mantle discontinuity under the Pacific

Seismic discontinuities in the mantle are indicators of its thermo-chemical state and offer clues to its dynamics. Ray-based seismic methods, though limited by the approximations made, have mapped mantle transition zone discontinuities in detail, but have yet to offer definitive conclusions on the presence and nature of mid-mantle discontinuities. Here, we show how to use a wave-equation-based imaging method, reverse-time migration of precursors to surface-reflected seismic body waves, to uncover both mantle transition zone and mid-mantle discontinuities, and interpret their physical nature. We observe a thinned mantle transition zone southeast of Hawaii, and a reduction in impedance contrast around 410 km depth in the same area, suggesting a hotter-than-average mantle in the region. Here, we furthermore reveal a 4000-5000 km-wide reflector in new images of the mid mantle below the central Pacific, at 950-1050 km depth. This deep discontinuity exhibits strong topography and generates reflections with polarity opposite to those originating at the 660 km discontinuity, implying an impedance reversal near 1000 km. We link this mid-mantle discontinuity to the upper reaches of deflected mantle plumes upwelling in the region. Reverse-time migration full-waveform imaging is a powerful approach to imaging Earth's interior, capable of broadening our understanding of its structure and dynamics and shrinking modeling uncertainties.

The influence of δ-(Al,Fe)OOH on seismic heterogeneities in Earth's lower mantle

The high-pressure phases of oxyhydroxides (δ-AlOOH, ε-FeOOH, and their solid solution), candidate components of subducted slabs, have wide stability fields, thus potentially influencing volatile circulation and dynamics in the Earth's lower mantle. Here, we report the elastic wave velocities of δ-(Al,Fe)OOH (Fe/(Al + Fe) = 0.13, δ-Fe13) to 79 GPa, determined by nuclear resonant inelastic X-ray scattering. At pressures below 20 GPa, a softening of the phonon spectra is observed. With increasing pressure up to the Fe3+ spin crossover (~ 45 GPa), the Debye sound velocity (vD) increases. At higher pressures, the low spin δ-Fe13 is characterized by a pressure-invariant vD. Using the equation of state for the same sample, the shear-, compressional-, and bulk-velocities (vS, vP, and vΦ) are calculated and extrapolated to deep mantle conditions. The obtained velocity data show that δ-(Al,Fe)OOH may cause low-vΦ and low-vP anomalies in the shallow lower mantle. At deeper depths, we find that this hydrous phase reproduces the anti-correlation between vS and vΦ reported for the large low seismic velocity provinces, thus serving as a potential seismic signature of hydrous circulation in the lower mantle.

New discovery of two seismite horizons challenges the Ries-Steinheim double-impact theory

The Nördlinger Ries and the Steinheim Basin are widely perceived as a Middle Miocene impact crater doublet. We discovered two independent earthquake-produced seismite horizons in North Alpine Foreland Basin deposits potentially related to both impacts. The older seismite horizon, demonstrated to be associated with the Ries impact, is overlain by distal impact ejecta in situ, forming a unique continental seismite-ejecta couplet within a distance of up to 180 km from the crater. The younger seismite unit, also produced by a major palaeo-earthquake, comprises clastic dikes that cut through the Ries seismite-ejecta couplet. The clastic dikes may have formed in response to the Steinheim impact, some kyr after the Ries impact, in line with paleontologic results that indicate a time gap of about 0.5 Myr between the Ries and Steinheim events. This interpretation suggests the Ries and Steinheim impacts represent two temporally separate events in Southern Germany that, thus, witnessed a double disaster in the Middle Miocene. The magnitude–distance relationship of seismite formation during large earthquakes suggests the seismic and destructive potential of impact-induced earthquakes may be underestimated.

Experimental evidence for silica-enriched Earth's lower mantle with ferrous iron dominant bridgmanite

Determination of the chemical composition of the Earth's mantle is of prime importance to understand the evolution, dynamics, and origin of the Earth. However, there is a lack of experimental data on sound velocity of iron-bearing Bridgmanite (Brd) under relevant high-pressure conditions of the whole mantle, which prevents constraints on the mineralogical model of the lower mantle. To uncover these issues, we have conducted sound-velocity measurement of iron-bearing Brd in a diamond-anvil cell (DAC) up to 124 GPa using Brillouin scattering spectroscopy. Here we show that the sound velocities of iron-bearing Brd throughout the whole pressure range of lower mantle exhibit an apparent linear reduction with the iron content. Our data fit remarkably with the seismic structure throughout the lower mantle with Fe²⁺-enriched Brd, indicating that the greater part of the lower mantle could be occupied by Fe²⁺-enriched Brd. Our lower-mantle model shows a distinctive Si-enriched composition with Mg/Si of 1.14 relative to the upper mantle (Mg/Si = 1.25), which implies that the mantle convection has been inefficient enough to chemically homogenize the Earth's whole mantle.

Combining Composition and Sonic Information Design in a New Electroacoustic Work

It might be suggested that composition and sonic information design are fundamentally different. However, some academic commentators and composers have explored the intersection between these disciplines. The authors presented one such work, Eonsounds: Fiamignano Gorge, at International Community of Auditory Display 2017. We argue that analysis of the aesthetic and informational choices in such hybrid works is essential to the development of sonic information design, with implications for the emergence of sonic information design as a subtype of human factors design. By acknowledging the relationship between aesthetics and information presentation, sound designers may develop designs that are safer and more user-friendly.

Lamination and mixing in three fundamental flow sequences driven by electromagnetic body forces

This article pursues the idea that the degree of striations, called lamination, could be engineered to complement stretching and to design new sequential mixers. It explores lamination and mixing in three new mixing sequences experimentally driven by electromagnetic body forces. To generate these three mixing sequences, Lorentz body forces are dynamically controlled to vary the flow geometry produced by a pair of local jets. The first two sequences are inspired from the "tendril and whorl" and "blinking vortex" flows. The third novel sequence is called the "cat's eyes flip." These three mixing sequences exponentially stretch and laminate material lines representing the interface between two domains to be mixed. Moreover, the mixing coefficient (defined as 1-σ(2)/σ(0)(2) where σ(2)/σ(0)(2) is the rescaled variance) and its rate grow exponentially before saturation. This saturation of the mixing process is related to the departure of the mixing rate from an exponential growth when the striations' thicknesses reach the diffusive length scale of the measurements or species and dyes. Incidentally, in our experiments, for the same energy or forcing input, the cat's eyes flip sequence has higher lamination, stretching, and mixing rates than the tendril and whorl and the blinking vortex sequences. These features show that bakerlike in situ mixers can be conceived by dynamically controlling a pair of local jets and by integrating lamination during stirring stages with persistent geometries. Combined with novel insights provided by the quantification of the lamination, this paper should offer perspectives for the development of new sequential mixers, possibly on all scales.

Atomic-scale visualization of antisite defects in LiFePO4

We visualize the antisite exchange defects in LiFePO4 crystals with an ordered olivine structure by using annular dark-field scanning transmission electron microscopy (STEM). A recognizable bright contrast is observed in some of the Li columns of STEM images in a sample annealed at a lower temperature, which directly demonstrates the disordered occupations by Fe atoms. Furthermore, such exchange defects appear to be locally aggregated rather than homogeneously dispersed in the lattice, although their overall concentration is fairly low. The present study emphasizes the significance of atomic-level observations for the defect distribution that cannot be predicted by macroscopic analytical methods.

Splitting of the 520-kilometer seismic discontinuity and chemical heterogeneity in the mantle

Seismic studies indicate that beneath some regions the 520-kilometer seismic discontinuity in Earth's mantle splits into two separate discontinuities (at approximately 500 kilometers and approximately 560 kilometers). The discontinuity near 500 kilometers is most likely caused by the (Mg,Fe)2SiO4 beta-to-gamma phase transformation. We show that the formation of CaSiO3 perovskite from garnet can cause the deeper discontinuity, and by determining the temperature dependence for this reaction we demonstrate that regional variations in splitting of the discontinuity arise from variability in the calcium concentration of the mantle rather than from temperature changes. This discontinuity therefore is sensitive to large-scale chemical heterogeneity. Its occurrence and variability yield regional information on the fertility of the mantle or the proportion of recycled oceanic crust.

Helium isotopic evidence for episodic mantle melting and crustal growth

The timing of formation of the Earth's continental crust is the subject of a long-standing debate, with models ranging from early formation with little subsequent growth, to pulsed growth, to steadily increasing growth. But most models do agree that the continental crust was extracted from the mantle by partial melting. If so, such crustal extraction should have left a chemical fingerprint in the isotopic composition of the mantle. The subduction of oceanic crust and subsequent convective mixing, however, seems to have largely erased this record in most mantle isotopic systems (for example, strontium, neodymium and lead). In contrast, helium is not recycled into the mantle because it is volatile and degasses from erupted oceanic basalts. Therefore helium isotopes may potentially preserve a clearer record of mantle depletion than recycled isotopes. Here I show that the spectrum of 4He/3He ratios in ocean island basalts appears to preserve the mantle's depletion history, correlating closely with the ages of proposed continental growth pulses. The correlation independently predicts both the dominant 4He/3He peak found in modern mid-ocean-ridge basalts, as well as estimates of the initial 4He/3He ratio of the Earth. The correspondence between the ages of mantle depletion events and pulses of crustal production implies that the formation of the continental crust was indeed episodic and punctuated by large, potentially global, melting events. The proposed helium isotopic evolution model does not require a primitive, undegassed mantle reservoir, and therefore is consistent with whole mantle convection.

Spin transition in magnesiowüstite in earth's lower mantle

Iron in the major lower mantle (LM) minerals undergoes a high spin (HS) to low spin (LS) transition at relevant pressures (23-135 GPa). Previous failures of standard first principles approaches to describe this phenomenon have hindered its investigation and the clarification of important consequences. Using a rotationally invariant formulation of LDA + U we report a successful study of this transition in low solute concentration magnesiowüstite, (Mg(1-x)Fe(x)(O), (x < 0.2) the second most abundant LM phase. We show that the HS-LS transition goes through an insulating (semiconducting) intermediate mixed spins (MS) state without discontinuous changes in properties, as seen experimentally. We show that the HS state crosses over smoothly to the LS state passing through an insulating MS state where properties change continuously, as seen experimentally.

The physical and chemical composition of the lower mantle

This article reviews some of the recent advances made within the field of mineral physics. In order to link the observed seismic and density structures of the lower mantle with a particular mineral composition, knowledge of the thermodynamic properties of the candidate materials is required. Determining which compositional model best matches the observed data is difficult because of the wide variety of possible mineral structures and compositions. State-of-the-art experimental and analytical techniques have pushed forward our knowledge of mineral physics, yet certain properties, such as the elastic properties of lower mantle minerals at high pressures and temperatures, are difficult to determine experimentally and remain elusive. Fortunately, computational techniques are now sufficiently advanced to enable the prediction of these properties in a self-consistent manner, but more results are required.A fundamental question is whether or not the upper and lower mantles are mixing. Traditional models that involve chemically separate upper and lower mantles cannot yet be ruled out despite recent conflicting seismological evidence showing that subducting slabs penetrate deep into the lower mantle and that chemically distinct layers are, therefore, unlikely.Recent seismic tomography studies giving three-dimensional models of the seismic wave velocities in the Earth also base their interpretations on the thermodynamic properties of minerals. These studies reveal heterogeneous velocity and density anomalies in the lower mantle, which are difficult to reconcile with mineral physics data.

Tungsten isotope evidence that mantle plumes contain no contribution from the Earth's core

Osmium isotope ratios provide important constraints on the sources of ocean-island basalts, but two very different models have been put forward to explain such data. One model interprets (187)Os-enrichments in terms of a component of recycled oceanic crust within the source material. The other model infers that interaction of the mantle with the Earth's outer core produces the isotope anomalies and, as a result of coupled (186)Os-(187)Os anomalies, put time constraints on inner-core formation. Like osmium, tungsten is a siderophile ('iron-loving') element that preferentially partitioned into the Earth's core during core formation but is also 'incompatible' during mantle melting (it preferentially enters the melt phase), which makes it further depleted in the mantle. Tungsten should therefore be a sensitive tracer of core contributions in the source of mantle melts. Here we present high-precision tungsten isotope data from the same set of Hawaiian rocks used to establish the previously interpreted (186)Os-(187)Os anomalies and on selected South African rocks, which have also been proposed to contain a core contribution. None of the samples that we have analysed have a negative tungsten isotope value, as predicted from the core-contribution model. This rules out a simple core-mantle mixing scenario and suggests that the radiogenic osmium in ocean-island basalts can better be explained by the source of such basalts containing a component of recycled crust.

Fractionation of elements by particle size of ashes ejected from Copahue Volcano, Argentina

The volcano Copahue, Neuquén province, Argentina has shown infrequent explosive eruptions since the 18th century. Recently, eruptive activity and seismicity were registered in the period July-October, 2000. As a consequence, ash clouds were dispersed by winds and affected Caviahue village located at about 9 km east of the volcano. Samples of deposited particles from this area were collected during this episode for their chemical analysis to determine elements of concern with respect to the health of the local population and its environment. Different techniques were used to evaluate the distribution of elements in four particle size ranges from 36 to 300 microm. X-ray powder diffraction (XRD) was selected to detect major components namely, minerals, silicate glass, fragments of rocks and sulfurs. Major and minor elements (Al, Ca, Cl, Fe, K, Mg, Mn, Na, S, Si and Ti), were detected by energy dispersive X ray analysis (EDAX). Trace element (As, Cd, Cr, Cu, Hg, Ni, Pb, Sb, U, V and Zn) content was quantified by inductively coupled plasma-mass spectrometry (ICP-MS). Nuclear activation analysis (NAA) was employed for the determination of Ce, Co, Cs, Eu, Hf, La, Lu, Rb, Sc, Sm, Ta and Yb. An enrichment was observed in the smallest size fraction of volcanic ashes for four elements (As, Cd, Cu and Sb) of particular interest from the environmental and human health point of view.

Chemical and seismological constraints on mantle heterogeneity

Recent seismological studies that use scattered waves to detect heterogeneities in the mantle reveal the presence of a small, distributed elastic heterogeneity in the lower mantle which does not appear to be thermal in nature. The characteristic size of these heterogeneities appears to be ca. 8 km, suggesting that they represent subducted recycled oceanic crust. With this stimulus, old ideas that the mantle is heterogeneous in structure, rather than stratified, are reinterpreted and a simple, end-member model for the heterogeneity structure is proposed. The volumetrically largest components in the model are recycled oceanic crust, which contains the heat-producing elements, and mantle depleted of these and other incompatible trace elements. About 10% of the mantle's mass is made up of recycled oceanic crust, which is associated with the observed small-scale seismic heterogeneity. The way this heterogeneity is distributed is in convectively stretched and thinned bodies ranging downwards in size from 8 km. With the present techniques to detect small bodies through scattering, only ca. 55% of the mantle's small-scale heterogeneities are detectable seismically.

Geodynamic and seismic constraints on the thermochemical structure and dynamics of convection in the deep mantle

We revisit a recent study by Forte & Mitrovica in which global geophysical observables associated with mantle convection were inverted and the existence of a strong increase in viscosity near a depth of 2000 km was inferred. Employing mineral-physics data and theory we also showed that, although there are chemical anomalies in the lowermost mantle, they are unable to inhibit the dominant thermal buoyancy of the deep-mantle mega-plumes below Africa and the Pacific Ocean. New Monte Carlo simulations are employed to explore the impact of uncertainties in current mineral-physics constraints on inferences of deep-mantle thermochemical structure. To explore the impact of the high-viscosity peak at a depth of 2000 km on the evolution of lower-mantle structure, we carried out time-dependent convection simulations. The latter show that the stability and longevity of the dominant long-wavelength heterogeneity in the lowermost mantle are controlled by this viscosity peak.

Zoned mantle convection

We review the present state of our understanding of mantle convection with respect to geochemical and geophysical evidence and we suggest a model for mantle convection and its evolution over the Earth's history that can reconcile this evidence. Whole-mantle convection, even with material segregated within the D" region just above the core-mantle boundary, is incompatible with the budget of argon and helium and with the inventory of heat sources required by the thermal evolution of the Earth. We show that the deep-mantle composition in lithophilic incompatible elements is inconsistent with the storage of old plates of ordinary oceanic lithosphere, i.e. with the concept of a plate graveyard. Isotopic inventories indicate that the deep-mantle composition is not correctly accounted for by continental debris, primitive material or subducted slabs containing normal oceanic crust. Seismological observations have begun to hint at compositional heterogeneity in the bottom 1000 km or so of the mantle, but there is no compelling evidence in support of an interface between deep and shallow mantle at mid-depth. We suggest that in a system of thermochemical convection, lithospheric plates subduct to a depth that depends - in a complicated fashion - on their composition and thermal structure. The thermal structure of the sinking plates is primarily determined by the direction and rate of convergence, the age of the lithosphere at the trench, the sinking rate and the variation of these parameters over time (i.e. plate-tectonic history) and is not the same for all subduction systems. The sinking rate in the mantle is determined by a combination of thermal (negative) and compositional buoyancy and as regards the latter we consider in particular the effect of the loading of plates with basaltic plateaux produced by plume heads. Barren oceanic plates are relatively buoyant and may be recycled preferentially in the shallow mantle. Oceanic plateau-laden plates have a more pronounced negative buoyancy and can more easily founder to the very base of the mantle. Plateau segregation remains statistical and no sharp compositional interface is expected from the multiple fate of the plates. We show that the variable depth subduction of heavily laden plates can prevent full vertical mixing and preserve a vertical concentration gradient in the mantle. In addition, it can account for the preservation of scattered remnants of primitive material in the deep mantle and therefore for the Ar and (3)He observations in ocean-island basalts.