Large Low Shear-wave Velocity Provinces Research Articles

Sluggish thermochemical basal mantle structures support their long-lived stability.

Large low shear-wave velocity provinces (LLSVPs) in the lowermost mantle are the largest geological structures on Earth, but their origin and age remain highly enigmatic. Geological constraints suggest the stability of the LLSVPs since at least 200 million years ago. Here, we conduct numerical modeling of mantle convection with plate-like behavior that yields a Pacific-like girdle of mantle downwelling which successfully forms two antipodal basal mantle structures similar to the LLSVPs. Our parameterized results optimized to reflect LLSVP features exhibit velocities for the basal mantle structures that are ~ 4 times slower than the ambient mantle if they are thermochemical, while the velocity is similar to the ambient mantle if purely thermal. The sluggish motion of the thermochemical basal mantle structures in our models permits the notion that geological data from hundreds of millions of years ago are related to modern LLSVPs as they are essentially stationary over such time scales.

Large-scale mantle heterogeneity as a legacy of plate tectonic supercycles

The Earth’s mantle is divided by the circum-Pacific subduction girdle into the African and Pacific domains, each featuring a large low-shear-wave-velocity province (LLSVP) in the lower mantle. However, how this hemispherical-scale mantle structure links to Earth’s plate tectonic evolution remains unclear. Previous geochemical work has suggested the presence of a north–south hemispheric subdivision, with large-scale mantle heterogeneities in the Southern Hemisphere, termed the DUPAL (Dupré and Allegre) anomaly. Here we compile elemental and isotopic data of both shallow-mantle-derived oceanic igneous rocks from mid-ocean ridges and deeper-mantle plume-related samples (ocean islands and oceanic plateaus) and analyse these using supervised machine learning classification methods. Data from both shallow- and deeper-mantle-sourced samples illustrate a consistent chemical dichotomy. Our results indicate that heterogeneities in the present-day shallow and deep mantle are not exclusively controlled by the north–south hemispheric DUPAL anomaly. Instead, they are consistent with a chemical dichotomy between the African and Pacific mantle domains and their associated LLSVPs. These observations can best be explained by tectonic supercycles over the past one billion years involving two supercontinents and two superoceans.

Periodicity in the Deccan Volcanism Modulated by Plume Perturbations at the Mid‐Mantle Transition Zone

AbstractIn peninsular India, the Deccan Traps record massive, continental‐scale volcanism in a sequence of magmatic events that corresponds with the timing of mass extinction at the Cretaceous‐Paleogene boundary. Although the Deccan volcanism is linked with the Réunion hotspot, the origin of its periodic magmatic pulses is still debated. We developed a numerical model replicating the geodynamic scenario of the African superplume underneath a moving Indian plate to explore the mechanism of magmatic pulse generation during the Deccan volcanism. Our model results revealed a connection between the Réunion hotspot and the African large low shear‐wave velocity province (LLSVP), suggesting that the pulses were produced from a thermochemical plume originated in the lower mantle. The ascending plume had stagnation at 660 km due to phase changes in the transition zone, and its head eventually underwent detachment from the tail under the influence of Indian plate movement to produce sequentially four major pulses (periodicity: 5–8 Ma), each giving rise to multiple secondary magmatic pulses at a time interval of ∼0.15–0.4 Ma. This study sheds a new light on the mechanism of periodic hotspot activities from a global perspective.

Heavy Mo isotope enrichment in the Pitcairn plume: Implications for the subduction cycle of anoxic sediments

Subduction redistributes elements between Earth's principal geochemical reservoirs, modifying the chemical composition of Earth's mantle, crust, atmosphere, and hydrosphere, and consequently having an impact on the evolution of life itself. Subduction of surface material that has been geochemically modified by low-temperature processes leads to mineralogical and chemical heterogeneities in mantle reservoirs over time and is recorded in modern ocean island basalts. One of the principal geochemical end members of the heterogeneous deep mantle, the enriched mantle 1 (EM-1) source of Pitcairn Island, has been attributed to the contribution of crustal material with vastly different chemical compositions and ages. The Mo isotope composition of lavas from Pitcairn Island constrains the nature of this recycled crustal component. Pitcairn lavas have elevated δ98/95Mo relative to the depleted mantle. The high δ98/95Mo is associated with high time-integrated 232Th/238U and 87Rb/86Sr, and low time-integrated 147Sm/144Nd and 238U/204Pb. These characteristics can be attributed to the recycling of nearly pristine pelagic sediments that were deposited in a Proterozoic anoxic deep-ocean into the sources of the Pitcairn Island lavas. The isotope composition of these lavas is similar to that of EM-1 hotspots from the South Atlantic, indicating the addition of reduced sediments in both of Earth's large low shear wave velocity provinces (LLSVPs). Consistent data from both locations imply that the subduction cycling of sedimentary redox-sensitive elements such as Mo, S, Se, and U into arc magmas was in these cases inefficient in the Precambrian and the chemical and isotopic signature of reduced sediments is preserved in the source of ocean island basalts bearing the EM-1 characteristics.

New evidence for the Ontong Java Nui hypothesis

The formation of the Ontong Java Nui super oceanic plateau (OJN), which is based on the model that the submarine Ontong Java Plateau (OJP), Manihiki Plateau (MP), and Hikurangi Plateau (HP) were once its contiguous fragments, could have been the largest globally consequential volcanic event in Earth’s history. This OJN hypothesis has been debated given the paucity of evidence, for example, the differences in crustal thickness, the compositional gap between MP and OJP basalts and the apparent older age of both plateaus relative to HP remain unresolved. Here we investigate the geochemical and 40Ar-39Ar ages of dredged rocks recovered from the OJP’s eastern margin. Volcanic rocks having compositions that match the low-Ti MP basalts are reported for the first time on the OJP and new ~ 96–116 Ma and 67–68 Ma 40Ar-39Ar age data bridge the temporal gap between OJP and HP. These results provide new evidence for the Ontong Java Nui hypothesis and a framework for an integrated tectonomagmatic evolution of the OJP, MP, and HP. The isotopic data imply four mantle components in the source of OJN that are also expressed in present-day Pacific hotspots sources, indicating origin from (and longevity of) the Pacific Large Low Shear-wave Velocity Province.

Formation of Thermochemical Heterogeneities by Core‐Mantle Interaction

AbstractEarth's core‐mantle boundary (CMB) shows a complex structure with various seismic anomalies such as the large low shear‐wave velocity provinces (LLSVPs) and ultra‐low velocity zones (ULVZs). As these structures are possibly induced by chemically distinct material forming a layer above the CMB, models of mantle convection made ad hoc assumptions to simulate the dynamics of this layer. In particular, density and mass were prescribed. Both conditions are critical for the dynamics but hardly constrained. Core‐mantle interaction is considered as one possible origin for this dense layer. For example, diffusion‐controlled enrichment of iron has been proposed. We here apply a chemical gradient between the mantle and the denser core to analyze the penetration of dense material into the mantle. As such, we employ 2D Cartesian models where a thermochemical layer at the base of the mantle develops self‐consistently by a diffusive chemical influx. Our simulations indicate that chemical diffusion is strongly affected by the convective mantle flow. This convection‐assisted diffusion yields a compositional influx mainly in the areas where slabs spread over the bottom boundary and sweep dense material aside to form accumulations with rising plumes atop. Like for a prescribed dense layer this process leads to chemically distinct piles, which are typically smaller (therefore more suited to explain ULVZs) but more persistent due to the constant chemical influx. Combining the influx scenario with the primordial layer can possibly explain the simultaneous existence of LLSVPs and ULVZs along with the observation of a core‐like isotopic composition in the mantle.

Contrasts in 2-D and 3-D system behaviour in the modelling of compositionally originating LLSVPs and a mantle featuring dynamically obtained plates

SUMMARY More than two decades of systematic investigation has made steady progress towards generating plate-like surface behaviour in models of vigorous mantle convection. Accordingly, properties required to obtain dynamic plates from mantle convection have become widely recognized and used in both 2-D and 3-D geometries. Improving our understanding of the properties required to obtain durable (or replenishable) deep mantle features with LLSVP-like characteristics has received interest for a period with similar longevity. Investigation ultimately focuses on discovering the properties able to produce the presence of a detached pair of 3-D features, distinct from the ambient mantle. Here, we assume the large low shear-wave velocity provinces (LLSVPs) have a chemical origin by incorporating a compositionally anomalous and intrinsically dense (CAID) mantle component comprising 2–3.5 per cent of the total mantle volume. The feedback between plate formation and the presence of a CAID mantle component is investigated in both 2-D and 3-D spherical geometries. We explore the impact of both an intrinsic contrast in density and viscosity for the CAID component, with the objective of finding system parameter values that encourage the formation of a pair of LLSVP-like assemblages and a surface that exhibits the principle features of terrestrial plate tectonics; including recognizable and narrowly focused divergent, convergent and (in 3-D) transform plate boundaries that separate 8–16 distinct plate interiors. We present the results of nine 2-D and 11 3-D calculations and show that for some of the cases examined, a pair of CAID material provinces can be freely obtained in 2-D cases while maintaining a surface characterized by plate-like behaviour. However, specifying the same system parameters in the 3-D model does not readily yield a pair of enduring provinces for any values of the parameters investigated. Moreover, the inclusion of the CAID component in the mantle can affect the global geotherm so that in comparison to the surface behaviour obtained for the initial condition isochemical model, the surface behaviour of the cases incorporating the dense component are less exemplary of plate tectonics. In general, CAID material components that are 3.75–5 per cent denser than the surrounding mantle (at surface temperatures), and up to a factor of 100 times greater in intrinsic viscosity, form layers populated by voids, or nodes connected by tendril-like ridges that reach across the core–mantle boundary (CMB), rather than distinct piles resembling LLSVPs. Due to its inherently heavy and stiff character, in equilibrated systems, we find the CAID material becomes especially hot so that the temperature-dependence of its density and viscosity results in reduced distinction between the intrinsically dense assemblages and the ambient mantle. Accordingly, the CAID material forms masses on the CMB that are relatively less dense (0.625–1.5 per cent) and viscous than the adjacent mantle material, in comparison to the percentage differences obtained at common temperatures. We find that by adjusting our yield stress model to account for the influence of the CAID material on the geotherm, a highly satisfactory plate-like surface can be re-attained, however, the formation of a pair of LLSVP-shaped masses remains elusive.

D″ anisotropy inverted from shear wave splitting intensity

The D″ layer, located at the bottom of the mantle, is an active thermochemical boundary layer. The upwelling of mantle plumes, as well as possible plate subduction in the D″ layer, could lead to large-scale material transformation and mineral deformation, which could result in significant seismic anisotropy. However, owing to limited observations and immense computational cost, the anisotropic structures and geodynamic mechanisms in the D″ layer remain poorly understood. In this study, we proposed a new inversion method for the seismic anisotropy in the D″ layer quantitatively with shear wave splitting intensities. We first proved the linearity of the splitting intensities under the ray-theory assumption. The synthetic tests showed that, with horizontal axes of symmetry and ray incidences lower than 30º in the D″ layer (typical SKS phase), the anisotropy is well resolved. We applied the method to the measured dataset in Africa and Western Europe, and obtained strong D″ anisotropy in the margins of the large low shear-wave velocity provinces and subducting slabs. The new method makes it possible to obtain D″ anisotropy, which provides essential constraints on the geodynamical processes at the base of the mantle.

Spatial Characteristics of Recycled and Primordial Reservoirs in the Deep Mantle

AbstractThe spatial distribution of the geochemical domains hosting recycled crust and primordial (high‐3He/4He) reservoirs, and how they are linked to mantle convection, are poorly understood. Two continent‐sized seismic anomalies located near the core‐mantle boundary—called the Large Low Shear Wave Velocity Provinces (LLSVPs)—are potential geochemical reservoir hosts. It has been suggested that high‐3He/4He hotspots are spatially confined to the LLSVPs, hotspots sampling recycled continental crust are associated with only one of the LLSVPs, and recycled continental crust shows no relationship with latitude. We reevaluate the links between LLSVPs and isotopic signatures of hotspot lavas using improved mantle flow models including plume conduit advection. While most hotspots with the highest‐3He/4He can indeed be traced to the LLSVP interiors, at least one high‐3He/4He hotspot, Yellowstone, is located outside of the LLSVPs. This suggests high‐3He/4He is not geographically confined to the LLSVPs. Instead, a positive correlation between hotspot buoyancy flux and maximum hotspot 3He/4He suggests that it is plume dynamics (i.e., buoyancy), not geography, which determines whether a dense, deep, and possibly widespread high‐3He/4He reservoir is entrained. We also show that plume‐fed EM hotspots (enriched mantle, with low‐143Nd/144Nd), signaling recycled continental crust, are spatially linked to both LLSVPs, and located primarily in the southern hemisphere. Lastly, we confirm that hotspots sampling HIMU (“high‐μ,” or high 238U/204Pb) domains are not spatially limited to the LLSVPs. These findings clarify and advance our understanding of deep mantle reservoir distributions, and we discuss how continental and oceanic crust subduction is consistent with the spatial decoupling of EM and HIMU.

The tungsten-182 record of kimberlites above the African superplume: Exploring links to the core-mantle boundary

Many volcanic hotspots are connected via ‘plume’ conduits to thermochemical structures with anomalously low seismic velocities at the core-mantle boundary. Basaltic lavas from some of these hotspots show anomalous daughter isotope abundances for the short-lived 129I-129Xe, 146Sm-142Nd, and 182Hf-182W radioactive decay systems, suggesting that their lower mantle sources contain material that dates back to Earth-forming events during the first 100 million years in solar system history. Survival of such ‘primordial’ remnants in Earth's mantle places important constraints on the evolution and inner workings of terrestrial planets.Here we report high-precision 182W/184W measurements for a large suite of kimberlite volcanic rocks from across the African tectonic plate, which for the past 250 million years has drifted over the most prominent thermochemical seismic anomaly at the core-mantle boundary. This so-called African LLSVP, or ‘large low shear-wave velocity province’, is widely suspected to store early Earth remnants and is implicated as the ultimate source of global Phanerozoic kimberlite magmatism. Our results show, however, that kimberlites from above the African LLSVP, including localities with lower mantle diamonds such as Letseng and Karowe Orapa A/K6, lack anomalous 182W signatures, with an average μ182W value of 0.0 ± 4.1 (2SD) for the 18 occurrences studied.If kimberlites are indeed sourced from the African LLSVP or superplume, then the extensive 182W evidence suggests that primordial or core-equilibrated mantle materials, which may contribute resolvable μ182W excesses or deficits, are only minor or locally concentrated components in the lowermost mantle, for example in the much smaller ‘ultra-low velocity zones’ or ULVZs. However, the lack of anomalous 182W may simply suggest that low-volume kimberlite magmas are not derived from hot lower mantle plumes. In this alternative scenario, kimberlite magmas originate from volatile-fluxed ambient convecting upper mantle domains beneath relatively thick and cold lithosphere from where previously ‘stranded’ lower mantle and transition zone diamonds can be plucked.

Distinct formation history for deep-mantle domains reflected in geochemical differences

The Earth’s mantle is currently divided into the African and Pacific domains, separated by the circum-Pacific subduction girdle, and each domain features a large low shear-wave velocity province (LLSVP) in the lower mantle. However, it remains controversial as to whether the LLSVPs have been stationary through time or dynamic, changing in response to changes in global subduction geometry. Here we compile radiogenic isotope data on plume-induced basalts from ocean islands and oceanic plateaus above the two LLSVPs that show distinct lead, neodymium and strontium isotopic compositions for the two mantle domains. The African domain shows enrichment by subducted continental material during the assembly and breakup of the supercontinent Pangaea, whereas no such feature is found in the Pacific domain. This deep-mantle geochemical dichotomy reflects the different evolutionary histories of the two domains during the Rodinia and Pangaea supercontinent cycles and thus supports a dynamic relationship between plate tectonics and deep-mantle structures. Earth’s deep-mantle domains are geochemically distinct. The African domain is enriched in subducted material, which suggests a different history from the Pacific domain and a dynamic relationship between plate tectonics and deep-mantle structures.

Plume interaction and mantle heterogeneity: A geochemical perspective

Mantle plumes originating from the Core-Mantle Boundary (CMB) or the Mantle Transition Zone (MTZ) play an important role in material transfer through Earth’s interior. The hotspot-related plumes originate through different mechanisms and have diverse processes of material transfer. Both the Morganian plumes and large low shear wave velocity provinces (LLSVPs) are derived from the D″ layer in the CMB, whereas the Andersonian plumes originate from the upper mantle. All plumes have a plume head at the Moho, although the LLSVPs have an additional plume head at the MTZ. We compare the geochemical characteristics of various plumes in an attempt to evaluate the material exchange between the plumes and mantle layers. The D″ layer, the LLSVPs and the Morganian plumes are consisted of subducted slab and post-perovskite from the lower mantle. Bridgmanite would crystallize during the upwelling process of the LLSVPs and the Morganian plumes in the lower mantle, and the residual is a basalt-trachyte suite. Unlike the Morganian plumes, the crystallization in the LLSVPs is insufficient that material accumulates beneath the MTZ to form a plume head. Typically, the secondary plumes above the plume head occur at the edge of the LLSVPs because it is easier for bridgmanite crystal separating from the plume head at the edge, and the residual material with low density upwells to form the secondary plumes. Meanwhile, Na and K are enriched during the long-term crystallization process, and then the basalt-phonolite suite appears in the LLSVPs. The geochemical characteristics of Andersonian plumes suggest that the basalt-rhyolite suite is the major component in the upper mantle. Meanwhile the basalt-rhyolite suite also appears in the LLSVPs and the Morganian plumes because of the assimilation and contamination in the plume head beneath the Moho.

The dynamics and impact of compositionally originating provinces in a mantle convection model featuring rheologically obtained plates

SUMMARY Previous geodynamic studies have indicated that the presence of a compositionally anomalous and intrinsically dense (CAID) mantle component can impact both core heat flux and surface features, such as plate velocity, number and size. Implementing spherical annulus geometry mantle convection models, we investigate the influence of intrinsically dense material in the lower mantle on core heat flux and the surface velocity field. The dense component is introduced into a system that features an established plate-like surface velocity field, and subsequently we analyse the evolution of the surface velocity as well as the interior thermal structure of the mantle. The distribution and mobility of the CAID material is investigated by varying its buoyancy ratio relative to the ambient mantle (ranging from 0.7 to 1.5), its total volume (3.5–10 per cent of the mantle volume) and its intrinsic viscosity (0.01–100 times the ambient mantle viscosity). We find at least three distinct distributions of the dense material can occur adjacent to the core–mantle boundary (CMB), including multiple piles of varying topography, a core enveloping layer and two diametrically opposed provinces (which can on occasion break into three distinct piles). The latter distribution mimics the morphology of the seismically observed large low shear wave velocity provinces (LLSVPs) and can occur over the entire range of CAID material viscosities. However, diametrically opposed provinces occur primarily in cases with CAID material buoyancy numbers of 0.7–0.85 (corresponding to contrasts in density between ambient and CAID material of 130 and 160 kg m−3, respectively) in our model (with an effective Rayleigh number of order 106). Steep and high topography piles are also obtained for cases featuring buoyancy ratios of 0.85 and viscosities 10–100 times that of the ambient mantle. An increase in relative density, as well as larger volumes of CAID material, lead to the development of a core enveloping layer. Our findings show that when two provinces are present core heat flux can be reduced by up to 50 per cent relative to cases in which CAID material is absent. Surface deformation quantified by Plateness is minimally influenced by variation of the properties of the dense material. Surface velocity is found to be reduced in general but mostly substantially in cases featuring high CAID material viscosities and large volumes (i.e. 10 per cent) or buoyancy ratios.

Primitive Helium Is Sourced From Seismically Slow Regions in the Lowermost Mantle

AbstractA major goal in Earth Science has been to understand how geochemical characteristics of lavas at the Earth's surface relate to the location and formation history of specific regions in the Earth's interior. For example, some of the strongest evidence for the preservation of primitive material comes from low 4He/3He ratios in ocean island basalts, but the location of the primitive helium reservoir(s) remains unknown. Here we combine whole‐mantle seismic tomography, simulations of mantle flow, and a global compilation of new and existing measurements of the 4He/3He ratios in ocean island basalts to constrain the source location of primitive 4He/3He material. Our geodynamic simulations predict the present‐day surface expression of plumes to be laterally offset from their lower mantle source locations. When this lateral offset is accounted for, a strong relationship emerges between minimum 4He/3He ratios in oceanic basalts and seismically slow regions, which are generally located within the two large low shear‐wave velocity provinces (LLSVPs). Conversely, no significant relationship is observed between maximum 208Pb*/206Pb* ratios and seismically slow regions in the lowermost mantle. These results indicate that primitive materials are geographically restricted to LLSVPs, while recycled materials are more broadly distributed across the lower mantle. The primitive nature of the LLSVPs indicates these regions are not composed entirely of recycled slabs, while complementary xenon and tungsten isotopic anomalies require the primitive portion of the LLSVPs to have formed during Earth's accretion, survived the Moon‐forming giant impact, and remained relatively unmixed during the subsequent 4.5 billion years of mantle convection.

Radial thermo-chemical structure beneath Western and Northern Pacific from seismic waveform inversion

The Earth's deep mantle seismic structure is dominated by two large low shear-wave velocity provinces (LLSVPs) located beneath Africa and the Pacific. While the existence of these structures has been attested by many studies and data sets, their detailed nature, purely thermal or thermo-chemical, is still a matter of debate. Discriminating between these hypotheses requires constraints independent from seismic velocity structure. Seismic shear-wave attenuation, measured by the quality factor QS, strongly depends on temperature but not (or weakly) on composition. It may bring key information on temperature, resolving in turn the trade-off between temperature and composition. Here, we invert seismic waveform data for radial models of shear-wave velocity anomalies (dlnVS), and QS at two different locations in the Pacific, and from a depth of 2000 km down to the core-mantle boundary (CMB). We show that Western Pacific (WP) models, sampling the western tip of the Pacific LLSVP and the Caroline plume, cannot be explained by thermal anomalies alone, but require excess in iron of ∼4.0% from the CMB up to 2600 km. By contrast, Northern Pacific models (NP), if unaffected by seismic anisotropy, may have a purely thermal origin. Based on these observations, we build radial thermo-chemical models at WP and NP.

Geodynamic origin of carbonatites from the absolute paleotectonic reconstructions

Geodynamic origin of carbonatites is debated for several decades. One of hypotheses links their origin to large-volume mantle plumes rising from the core-mantle boundary (CMB). Some evidence exists for temporal and spatial relationships between the occurrences of carbonatites and large igneous provinces (LIPs), and both carbonatites and LIPs can be related to mantle plumes. A good example is the carbonatites of the Maymecha-Kotuy Province in the Polar Siberia, which were formed at the same time as the Siberian superplume event at ca. 250 Ma. In this study we use a recently published absolute plate kinematic modelling to reconstruct the position of 155 Phanerozoic carbonatites at the time of their emplacement. We demonstrate that 69% of carbonatites may be projected onto the central or peripheral parts of the large low shear-wave velocity provinces (LLSVPs) in the lowermost mantle. This correlation provides a strong evidence for the link between the carbonatite genesis and the locations of deep-mantle plumes. A large group of carbonatites (31%) has no obvious links to LLSVPs and, on the contrary, they plot above the “faster-than-average S-wave” zones in the deep mantle, currently located beneath North and Central America and China. We propose that their origin may be associated with remnants of subducted slabs in the mantle.

A Mantle Plume Origin for the Scandinavian Dyke Complex: A “Piercing Point” for 615 Ma Plate Reconstruction of Baltica?

AbstractThe origin of Large Igneous Provinces (LIPs) associated with continental breakup and the reconstruction of continents older than ca. 320 million years (pre‐Pangea) are contentious research problems. Here we study the petrology of a 615–590 Ma dolerite dyke complex that intruded rift basins of the magma‐rich margin of Baltica and now is exposed in the Scandinavian Caledonides. These dykes are part of the Central Iapetus Magmatic Province (CIMP), a LIP emplaced in Baltica and Laurentia during opening of the Caledonian Wilson Cycle. The >1,000‐km‐long dyke complex displays lateral geochemical zonation from enriched to depleted basaltic compositions from south to north. Geochemical modelling of major and trace elements shows these compositions are best explained by melting hot mantle 75–250 °C above ambient mantle. Although the trace element modelling solutions are nonunique, the best explanation involves melting a laterally zoned mantle plume with enriched and depleted peridotite lithologies, similar to present‐day Iceland and to the North Atlantic Igneous Province. The origin of CIMP appears to have involved several mantle plumes. This is best explained if rifting and breakup magmatism coincided with plume generation zones at the margins of a Large Low Shear‐wave Velocity Province (LLSVP) at the core mantle boundary. If the LLSVPs are quasi‐stationary back in time as suggested in recent geodynamic models, the CIMP provides a guide for reconstructing the paleogeography of Baltica and Laurentia 615 million years ago to the LLSVP now positioned under the Pacific Ocean. Our results provide a stimulus for using LIPs as piercing points for plate reconstructions.

Tungsten isotopes in mantle plumes: Heads it's positive, tails it's negative

The lowermost mantle is driven to Earth's surface by mantle plumes, providing a volcanic record of its structure and composition. Plumes comprise a head and tail, which melt to form large igneous provinces (LIPs) and ocean island basalts (OIBs), respectively. Recent analyses have shown that LIPs and OIBs exhibit tungsten (W) isotope heterogeneity that was created in the first ∼60 million years of our solar system's evolution. Moreover, the isotopic signature found in LIPs differs to that found in OIBs, revealing that the melt products of plume heads must be dominated by a different ancient mantle reservoir to that of plume tails. However, existing geodynamical studies suggest that plume heads and tails sample the same deep-mantle source region and, therefore, cannot account for any systematic differences in composition. Here, we present a suite of numerical simulations of thermo-chemical plumes and an isotopic model for W sources in the mantle. Our results demonstrate that the W isotope systematics of LIPs and OIBs can, under certain conditions, arise as a dynamical consequence of plumes forming in a heterogeneous, thermo-chemical boundary layer. We find that ultra low-velocity zones (ULVZs), which sit on the core–mantle boundary (CMB), likely contribute to the chemical diversity observed in OIBs but not LIPs, while any dense components residing inside large low shear-wave velocity provinces (LLSVPs) may contribute to both. This study places geochemical observations from Earth's surface in a geodynamically consistent framework and illuminates their relationship with seismically imaged features of the deep mantle.

Geochemistry and Distribution of Recycled Domains in the Mantle Inferred From Nd and Pb Isotopes in Oceanic Hot Spots: Implications for Storage in the Large Low Shear Wave Velocity Provinces

AbstractSubduction of continental and oceanic crust is thought to give rise to geochemically distinct reservoirs in the mantle called EM (enriched mantle) and HIMU (high μ = 238U/204Pb), respectively. However, the locations of EM and HIMU domains in the Earth's interior are poorly constrained. We explore the geographic distribution of extreme EM (143Nd/144Nd ≤ 0.512630) and HIMU (206Pb/204Pb ≥ 20) geochemical signatures in ocean island basalts erupted at hot spots, highlighting three observations. First, hot spots geographically associated with the two large low shear wave velocity provinces (LLSVPs) have a similar range of EM compositions. If these hot spots are sourced by the LLSVPs via upwelling plumes, this observation is consistent with the hypothesis that the LLSVPs formed by similar processes and have similar geodynamic histories. Second, the EM and HIMU domains exhibit different latitudinal zonation: oceanic hot spots with the most extreme EM compositions (143Nd/144Nd < 0.5125) are concentrated in the Southern Hemisphere (14°S to 52°S latitude), while oceanic hot spots hosting HIMU compositions are found primarily near the tropical latitudes (38°N to 29°S). Third, all 13 oceanic hot spots with EM compositions (143Nd/144Nd ≤ 0.512630) are geographically associated with the LLSVPs; oceanic hot spots located far from the LLSVPs exhibit only non‐EM (143Nd/144Nd > 0.512630) compositions. In contrast, the HIMU domains do not show a clear geographic association with the LLSVPs. Therefore, EM and HIMU domains in the Earth's mantle exhibit different spatial distributions. This may reflect differences in subduction inputs of these two components, or differences in how they segregate or accumulate in the deep Earth.

Deformation, crystal preferred orientations, and seismic anisotropy in the Earth's D″ layer

We use a forward multiscale model that couples atomistic modeling of intracrystalline plasticity mechanisms (dislocation glide ± twinning) in MgSiO3 post-perovskite (PPv) and periclase (MgO) at lower mantle pressures and temperatures to polycrystal plasticity simulations to predict crystal preferred orientations (CPO) development and seismic anisotropy in D″. We model the CPO evolution in aggregates of 70% PPv and 30% MgO submitted to simple shear, axial shortening, and along corner-flow streamlines, which simulate changes in flow orientation similar to those expected at the transition between a downwelling and flow parallel to the core–mantle boundary (CMB) within D″ or between CMB-parallel flow and upwelling at the borders of the large low shear wave velocity provinces (LLSVP) in the lowermost mantle. Axial shortening results in alignment of PPv [010] axes with the shortening direction. Simple shear produces PPv CPO with a monoclinic symmetry that rapidly rotates towards parallelism between the dominant [100](010) slip system and the macroscopic shear. These predictions differ from MgSiO3 post-perovskite textures formed in diamond-anvil cell experiments, but agree with those obtained in simple shear and compression experiments using CaIrO3 post-perovskite. Development of CPO in PPv and MgO results in seismic anisotropy in D″. For shear parallel to the CMB, at low strain, the inclination of ScS, Sdiff, and SKKS fast polarizations and delay times vary depending on the propagation direction. At moderate and high shear strains, all S-waves are polarized nearly horizontally. Downwelling flow produces Sdiff, ScS, and SKKS fast polarization directions and birefringence that vary gradually as a function of the back-azimuth from nearly parallel to inclined by up to 70° to CMB and from null to ∼5%. Change in the flow to shear parallel to the CMB results in dispersion of the CPO, weakening of the anisotropy, and strong azimuthal variation of the S-wave splitting up to 250 km from the corner. Transition from horizontal shear to upwelling also produces weakening of the CPO and complex seismic anisotropy patterns, with dominantly inclined fast ScS and SKKS polarizations, over most of the upwelling path. Models that take into account twinning in PPv explain most observations of seismic anisotropy in D″, but heterogeneity of the flow at scales <1000 km is needed to comply with the seismological evidence for low apparent birefringence in D″.