Widespread PREMA in the upper mantle indicated by low-degree basaltic melts.

Studies of ocean island basalts have identified a Prevalent Mantle (PREMA) component as a fundamental feature of mantle geochemical arrays; however, its origin and distribution are highly controversial, including its potential link to plumes sourced in low-shear-wave velocity provinces (LLSVPs) above the core-mantle boundary. In this study, we interrogate the compositional systematics of ~ 3500 Cenozoic oceanic and continental sodic basalts to provide insights into the origin and distribution of PREMA. We find that low-degree basaltic melts with high Nb concentrations located away from deep-mantle plumes have PREMA-like Sr-Nd-Hf isotopic signatures, implying that PREMA is highly fusible and not exclusively associated with LLSVPs. Geochemical modelling and mantle convection simulations indicate that PREMA could have been generated soon after Earth accretion, experiencing only minimal melting or enrichment, and then scattered throughout the upper mantle, rather than being the result of mixing between depleted and enriched mantle components.

Slab temperature controls on the Tonga double seismic zone and slab mantle dehydration.

Double seismic zones are two-layered distributions of intermediate-depth earthquakes that provide insight into the thermomechanical state of subducting slabs. We present new precise hypocenters of intermediate-depth earthquakes in the Tonga subduction zone obtained using data from local island-based, ocean-bottom, and global seismographs. The results show a downdip compressional upper plane and a downdip tensional lower plane with a separation of about 30 km. The double seismic zone in Tonga extends to a depth of about 300 km, deeper than in any other subduction system. This is due to the lower slab temperatures resulting from faster subduction, as indicated by a global trend toward deeper double seismic zones in colder slabs. In addition, a line of high seismicity in the upper plane is observed at a depth of 160 to 280 km, which shallows southward as the convergence rate decreases. Thermal modeling shows that the earthquakes in this "seismic belt" occur at various pressures but at a nearly constant temperature, highlighting the important role of temperature in triggering intermediate-depth earthquakes. This seismic belt may correspond to regions where the subducting mantle first reaches a temperature of ~500°C, implying that metamorphic dehydration of mantle minerals in the slab provides water to enhance faulting.

Trench-parallel flow and seismic anisotropy in the Mariana and Andean subduction systems.

Shear-wave splitting measurements above the mantle wedge of the Mariana and southern Andean subduction zones show trench-parallel seismically fast directions close to the trench and abrupt rotations to trench-perpendicular anisotropy in the back arc. These patterns of seismic anisotropy may be caused by three-dimensional flow associated with along-strike variations in slab geometry. The Mariana and Andean subduction systems are associated with the largest along-strike variations of slab geometry observed on Earth and are ideal for testing the link between slab geometry and solid-state creep processes in the mantle. Here we show, with fully three-dimensional non-newtonian subduction zone models, that the strong curvature of the Mariana slab and the transition to shallow slab dip in the Southern Andes give rise to strong trench-parallel stretching in the warm-arc and warm-back-arc mantle and to abrupt rotations in stretching directions that are accompanied by strong trench-parallel stretching. These models show that the patterns of shear-wave splitting observed in the Mariana and southern Andean systems may be caused by significant three-dimensional flow induced by along-strike variations in slab geometry.

Multiple volcanic episodes of flood basalts caused by thermochemical mantle plumes.

The hypothesis that a single mushroom-like mantle plume head can generate a large igneous province within a few million years has been widely accepted. The Siberian Traps at the Permian-Triassic boundary and the Deccan Traps at the Cretaceous-Tertiary boundary were probably erupted within one million years. These large eruptions have been linked to mass extinctions. But recent geochronological data reveal more than one pulse of major eruptions with diverse magma flux within several flood basalts extending over tens of million years. This observation indicates that the processes leading to large igneous provinces are more complicated than the purely thermal, single-stage plume model suggests. Here we present numerical experiments to demonstrate that the entrainment of a dense eclogite-derived material at the base of the mantle by thermal plumes can develop secondary instabilities due to the interaction between thermal and compositional buoyancy forces. The characteristic timescales of the development of the secondary instabilities and the variation of the plume strength are compatible with the observations. Such a process may contribute to multiple episodes of large igneous provinces.

Numerical models, geochemistry and the zero-paradox noble-gas mantle.

Numerical models of whole-mantle convection demonstrate that degassing of the mantle is an inefficient process, resulting in ca. 50% of the (40)Ar being degassed from the mantle system. In this sense the numerical simulations are consistent with the (40)Ar mass balance between the atmosphere and mantle reservoir. These models, however, are unable to preserve the large-scale heterogeneity predicted by models invoking geochemical layering of the mantle system. We show that the three most important noble-gas constraints on the geochemically layered mantle are entirely dependent on the (3)He concentration of the convecting mantle derived from the (3)He flux into the oceans and the average ocean-crust generation rate. A factor of 3.5 increase in the convecting-mantle noble-gas concentration removes all requirements for: a (3)He flux into the upper mantle from a deeper high (3)He source; a boundary in the mantle capable of separating heat from helium; and a substantial deep-mantle reservoir to contain a hidden (40)Ar rich reservoir. We call this model concentration for the convecting mantle the 'zero-paradox' concentration. The time-integrated flux of (3)He into the oceans is a robust observation, but only representative of the ocean-floor activity over the last 1000 years. In contrast, ocean-floor generation occurs over tens of millions of years. We argue that combining these two observations to obtain the (3)He concentration of the mantle beneath mid-ocean ridges is unsound. Other indicators of mantle (3)He concentration suggest that the real value may be at least a factor of two higher. As the zero-paradox concentration is approached, the noble-gas requirement for mantle layering is removed. We further consider the role that recycled material plays in ocean-island-basalt generation and show that a source with high (3)He and (3)He/(4)He must exist within the mantle. Nevertheless, only a small amount of this material is required to generate both the observed ocean-island (3)He/(4)He ratios and the concentrations inferred from basalt samples for this mantle source.

Development of anisotropic structure in the Earth's lower mantle by solid-state convection.

Seismological observations reveal highly anisotropic patches at the bottom of the Earth's lower mantle, whereas the bulk of the mantle has been observed to be largely isotropic. These patches have been interpreted to correspond to areas where subduction has taken place in the past or to areas where mantle plumes are upwelling, but the underlying cause for the anisotropy is unknown-both shape-preferred orientation of elastically heterogeneous materials and lattice-preferred orientation of a homogeneous material have been proposed. Both of these mechanisms imply that large-strain deformation occurs within the anisotropic regions, but the geodynamic implications of the mechanisms differ. Shape-preferred orientation would imply the presence of large elastic (and hence chemical) heterogeneity whereas lattice-preferred orientation requires deformation at high stresses. Here we show, on the basis of numerical modelling incorporating mineral physics of elasticity and development of lattice-preferred orientation, that slab deformation in the deep lower mantle can account for the presence of strong anisotropy in the circum-Pacific region. In this model-where development of the mineral fabric (the alignment of mineral grains) is caused solely by solid-state deformation of chemically homogeneous mantle material-anisotropy is caused by large-strain deformation at high stresses, due to the collision of subducted slabs with the core-mantle boundary.