Tomographic imaging of subducted lithosphere below northwest Pacific island arcs

Extreme redox variations in a superdeep diamond from a subducted slab

The introduction of volatile-rich subducting slabs to the mantle may locally generate large redox gradients, affecting phase stability, element partitioning and volatile speciation¹. Here we investigate the redox conditions of the deep mantle recorded in inclusions in a diamond from Kankan, Guinea. Enstatite (former bridgmanite), ferropericlase and a uniquely Mg-rich olivine (Mg# 99.9) inclusion indicate formation in highly variable redox conditions near the 660 km seismic discontinuity. We propose a model involving dehydration, rehydration and dehydration in the underside of a warming slab at the transition zone-lower mantle boundary. Fluid liberated by dehydration in a crumpled slab, driven by heating from the lower mantle, ascends into the cooler interior of the slab, where the H₂O is sequestered in new hydrous minerals. Consequent fractionation of the remaining fluid produces extremely reducing conditions, forming Mg-end-member ringwoodite. This fractionating fluid also precipitates the host diamond. With continued heating, ringwoodite in the slab surrounding the diamond forms bridgmanite and ferropericlase, which is trapped as the diamond grows in hydrous fluids produced by dehydration of the warming slab.

Slab morphology and deformation beneath Izu-Bonin

Seismic tomography provides unique constraints on the morphology, the deformation, and (indirectly) the rheology of subducting slabs. We use teleseismic double-difference P-wave tomography to image with unprecedented clarity the structural complexity of the Izu-Bonin slab. We resolve a tear in the slab in the mantle transition zone (MTZ) between 26.5° N and 28° N. North of the tear, the slab is folded in the MTZ. Immediately above the fold hinge, a zone of reduced P-wavespeed may result from viscous dissipation within an incipient shear zone. To the south of the tear, the slab overturns and lies flat at the base of the MTZ. The ~680 km deep 2015 Bonin earthquake (Mw~7.9) is located at the northernmost edge of the overturning part of the slab. The localised tearing, shearing and buckling of the Izu-Bonin slab indicates that it remains highly viscous throughout the upper mantle and transition zone.

In the 1000 km long Izu-Bonin subduction zone to the south of Tokyo, the Pacific Plate descends beneath the Philippine Sea Plate. Here the authors use teleseismic double-difference tomography to image the complex morphology of the Izu-Bonin slab, especially in the mantle transition zone.

A Metastable Fo-III Wedge in Cold Slabs Subducted to the Lower Part of the Mantle Transition Zone: A Hypothesis Based on First-Principles Simulations

The metastable olivine (Ol) wedge hypothesis assumes that Ol may exist as a metastable phase at the P conditions of the mantle transition zone (MTZ) and even deeper regions due to inhibition of the phase transitions from Ol to wadsleyite and ringwoodite caused by low T in the cold subducting slabs. It is commonly invoked to account for the stagnation of the descending slabs, deep focus earthquakes and other geophysical observations. In the last few years, several new structures with the forsterite (Fo) composition, namely Fo-II, Fo-III and Fo-IV, were either experimentally observed or theoretically predicted at very low T conditions. They may have important impacts on the metastable Ol wedge hypothesis. By performing first-principles calculations, we have systematically examined their crystallographic characteristics, elastic properties and dynamic stabilities from 0 to 100 GPa, and identified the Fo-III phase as the most likely metastable phase to occur in the cold slabs subducted to the depths equivalent to the lower part of the MTZ (below the ~600 km depth) and even the lower mantle. As disclosed by our theoretical simulations, the Fo-III phase is a post-spinel phase (space group Cmc21), has all cations in sixfold coordination at P < ~60 GPa, and shows dynamic stability for the entire P range from 0 to 100 GPa. Further, our static enthalpy calculations have suggested that the Fo-III phase may directly form from the Fo material at ~22 GPa (0 K), and our high-T phase relation calculations have located the Fo/Fo-III phase boundary at ~23.75 GPa (room T) with an averaged Clapeyron slope of ~−1.1 MPa/K for the T interval from 300 to 1800 K. All these calculated phase transition pressures are likely overestimated by ~3 GPa because of the GGA method used in this study. The discrepancy between our predicted phase transition P and the experimental observation (~58 GPa at 300 K) can be explained by slow reaction rate and short experimental durations. Taking into account the P-T conditions in the cold downgoing slabs, we therefore propose that the Fo-III phase, rather than the Ol, highly possibly occurs as the metastable phase in the cold slabs subducted to the P conditions of the lower part of the MTZ (below the ~600 km depth) and even the lower mantle. In addition, our calculation has showed that the Fo-III phase has higher bulk seismic velocity, and thus may make important contributions to the high seismic speeds observed in the cold slabs stagnated near the upper mantle-lower mantle boundary. Future seismic studies may discriminate the effects of the Fo-III phase and the low T. Surprisingly, the Fo-III phase will speed up, rather than slow down, the subducting process of the cold slabs, if it metastably forms from the Ol. In general, the Fo-III phase has a higher density than the warm MTZ, but has a lower density than the lower mantle, as suggested by our calculations.

Tidal tomography constrains Earth's deep-mantle buoyancy

Earth's body tide-also known as the solid Earth tide, the displacement of the solid Earth's surface caused by gravitational forces from the Moon and the Sun-is sensitive to the density of the two Large Low Shear Velocity Provinces (LLSVPs) beneath Africa and the Pacific. These massive regions extend approximately 1,000 kilometres upward from the base of the mantle and their buoyancy remains actively debated within the geophysical community. Here we use tidal tomography to constrain Earth's deep-mantle buoyancy derived from Global Positioning System (GPS)-based measurements of semi-diurnal body tide deformation. Using a probabilistic approach, we show that across the bottom two-thirds of the two LLSVPs the mean density is about 0.5 per cent higher than the average mantle density across this depth range (that is, its mean buoyancy is minus 0.5 per cent), although this anomaly may be concentrated towards the very base of the mantle. We conclude that the buoyancy of these structures is dominated by the enrichment of high-density chemical components, probably related to subducted oceanic plates or primordial material associated with Earth's formation. Because the dynamics of the mantle is driven by density variations, our result has important dynamical implications for the stability of the LLSVPs and the long-term evolution of the Earth system.

A perovskitic lower mantle inferred from high-pressure, high-temperature sound velocity data

The determination of the chemical composition of Earth's lower mantle is a long-standing challenge in earth science. Accurate knowledge of sound velocities in the lower-mantle minerals under relevant high-pressure, high-temperature conditions is essential in constraining the mineralogy and chemical composition using seismological observations, but previous acoustic measurements were limited to a range of low pressures and temperatures. Here we determine the shear-wave velocities for silicate perovskite and ferropericlase under the pressure and temperature conditions of the deep lower mantle using Brillouin scattering spectroscopy. The mineralogical model that provides the best fit to a global seismic velocity profile indicates that perovskite constitutes more than 93 per cent by volume of the lower mantle, which is a much higher proportion than that predicted by the conventional peridotitic mantle model. It suggests that the lower mantle is enriched in silicon relative to the upper mantle, which is consistent with the chondritic Earth model. Such chemical stratification implies layered-mantle convection with limited mass transport between the upper and the lower mantle.

Mineral phases of the Earth´s mantle

Our knowledge of the structure of the Earth´s interior has been obtained by analysing seismic waves that travel in the Earth, and the reference Earth global models used by geophysicists are essentially seismological. Depth profiles of the seismic waves velocities reveal that the deep Earth is divided in several shells, separated by velocity and density discontinuities. The main discontinuity located at a depth of 2900 km corresponds to the transition between the mantle and the core. The Earth´s mantle can be further divided into the upper mantle and the lower mantle, with a transition zone characterised by two prominent increases in velocities observed at 410- and 660-km depths. This article will be focused on the mineral phases of the Earth´s mantle. The interpretation of seismological models in terms of chemical composition and temperature relies on the knowledge of the nature, structure and elastic properties of the candidate materials. We will describe to what extent recent advances in experimental mineral physics and X-ray diffraction have yielded essential knowledge on the structure and high-pressure high-temperature behaviour of pertinent materials, and major improvements in our understanding of the chemical and mineralogical composition of the Earth´s mantle.

Domain Graph of Arabidopsis Proteome by Comparative Analysis

The domain graph of domains and domain combinations of Arabidopsis thaliana is established based on pfam 14.0 database and analyzed via comparison with 10 eukaryotic, 30 bacterial, and 16 archaeal proteomes. The comparative analysis of the domain graphs provides a useful platform for revealing global insights on the evolution of plant kingdom. More importantly, it is a powerful tool for searching not only the possible new function of both plant-specific and nonspecific domains via specific domain combinations in Arabidopsis thaliana but also the functional role of unknown domains. As an example, we present the functional link between ubiquitin and Myb_DNA-binding domains via Bromodomain as the plant specific evidence for the association between transcription and ubiquitin. We further show that PentatricoPeptide Repeats (PPR) proteins have plant-specific links with a wide variety of domains responsible for RNA binding/metabolism, modulation of protein-protein interactions, ubiquitin-conjugation, cell growth/maintenance, catalysis, and others. This further supports the recently proposed association of PPR proteins with specific RNA transcripts and defined effector proteins. Moreover, the domain graph built from tissue-specific genes is frequently associated with DNA binding domains, suggesting that the differentiation of tissue cell types is contributed mostly by tissue-specific transcriptional process. DOGMA (DOmain Graph via coMparitive analysis for Arabidopsis thaliana) is available on-line with a variety of search tools at http://theory.med.buffalo.edu/DOGMA. The database, which allows user-specified search for plant specific domains and their combinations, will be useful as an additional tool for annotation of the proteins that play specific roles in plants and other organisms.

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.

The evolution of mantle mixing

We present a geochemical overview of the canonic model that suggests a two-layer mantle for most of the Earth's history. A change in the Rayleigh number may have modified the convection and now allows the subduction of oceanic plates into the lower mantle, which was not the case in the past. The measurement of stirring time in the source of mid-ocean-ridge basalt, together with Xe- and Pb-isotopic ratios in the mid-ocean-ridge-basalt source, suggests that the upper mantle is separated into two domains, one above the 400 km discontinuity (asthenosphere) with rapid mixing and short residence time, and another between the 400 and 670 km discontinuities with sluggish mixing and a residence time of ca. 1.5 x 10(9) yr.

Seismic evidence for catastrophic slab loss beneath Kamchatka

In the northwest Pacific Ocean, a sharp corner in the boundary between the Pacific plate and the North American plate joins a subduction zone running along the southern half of the Kamchatka peninsula with a region of transcurrent motion along the western Aleutian arc. Here we present images of the seismic structure beneath the Aleutian-Kamchatka junction and the surrounding region, indicating that: the subducting Pacific lithosphere terminates at the Aleutian-Kamchatka junction; no relict slab underlies the extinct northern Kamchatka volcanic arc; and the upper mantle beneath northern Kamchatka has unusually slow shear wavespeeds. From the tectonic and volcanic evolution of Kamchatka over the past 10 Myr (refs 3-5) we infer that at least two episodes of catastrophic slab loss have occurred. About 5 to 10 Myr ago, catastrophic slab loss shut down island-arc volcanic activity north of the Aleutian-Kamchatka junction. A later episode of slab loss, since about 2 Myr ago, seems to be related to the activity of the world's most productive island-arc volcano, Klyuchevskoy. Removal of lithospheric mantle is commonly discussed in the context of a continental collision, but our findings imply that episodes of slab detachment and loss are also important agents in the evolution of oceanic convergent margins.

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.

Subduction and slab detachment in the Mediterranean-Carpathian region

Seismic tomography models of the three-dimensional upper mantle velocity structure of the Mediterranean-Carpathian region provide a better understanding of the lithospheric processes governing its geodynamical evolution. Slab detachment, in particular lateral migration of this process along the plate boundary, is a key element in the lithospheric dynamics of the region during the last 20 to 30 million years. It strongly affects arc and trench migration, and causes along-strike variations in vertical motions, stress fields, and magmatism. In a terminal-stage subduction zone, involving collision and suturing, slab detachment is the natural last stage in the gravitational settling of subducted lithosphere.