Geophysical and Geochemical Evidence for Deep Temperature Variations Beneath Mid-Ocean Ridges

Magmatism controls global oceanic transform fault topography

Oceanic transform faults play an essential role in plate tectonics. Yet to date, there is no unifying explanation for the global trend in broad-scale transform fault topography, ranging from deep valleys to shallow topographic highs. Using three-dimensional numerical models, we find that spreading-rate dependent magmatism within the transform domain exerts a first-order control on the observed spectrum of transform fault depths. Low-rate magmatism results in deep transform valleys caused by transform-parallel tectonic stretching; intermediate-rate magmatism fully accommodates far-field stretching, but strike-slip motion induces across-transform tension, producing transform strength dependent shallow valleys; high-rate magmatism produces elevated transform zones due to local compression. Our models also address the observation that fracture zones are consistently shallower than their adjacent transform fault zones. These results suggest that plate motion change is not a necessary condition for reproducing oceanic transform topography and that oceanic transform faults are not simple conservative strike-slip plate boundaries.

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.

Deep mantle origin of large igneous provinces and komatiites

Large igneous provinces (LIPs) resulted from intraplate magmatic events mobilizing volumes of magma up to several million cubic kilometers. LIPs and lavas with deep mantle sources have compositions ranging from komatiites found in Archean greenstone belts to basalts and picrites in Phanerozoic flood basalt and recent oceanic islands. In this study, we identify the mantle conditions appropriate to each type of lava based on an experimental study of the melting of pyrolite. The depth of the mantle source decreases from 600 to 700 km for the oldest komatiites to 100 to 300 km for picrites and basalts, and the extent of mantle melting ranges from 10 to 50%. We develop a geodynamical model that explains the origin of the hot mantle plumes capable of generating these melting P-T conditions. Within a superadiabatic temperature gradient persisting in the deep mantle, the ascent of hot mantle plumes creates excess temperatures up to 250 to 300 K by adiabatic decompression.

Discovery of distinct lithosphere-asthenosphere boundary and the Gutenberg discontinuity in the Atlantic Ocean

The plate tectonic theory requires a rigid lithosphere floating over a weak asthenosphere, separated by the lithosphere-asthenosphere boundary, which has been sometimes interpreted as the Gutenberg discontinuity. Using a deep seismic reflection technique, we report the presence of two continuous reflections covering 27 Ma to 58 Ma oceanic lithosphere in the Atlantic Ocean. We find that the upper reflection deepens with age and follows the ~1250°C isotherm, whereas the deeper reflection lies at a constant depth of ~75 km. We suggest that the upper reflection represents the thermally controlled lithosphere-asthenosphere boundary, whereas the lower reflection is the Gutenberg discontinuity, a frozen-in dehydration boundary separating the dry mantle melting region above from the hydrated mantle below formed at the ridge axis. We also find that thermal mantle anomalies rejuvenate the lithosphere, uplift the lithosphere-asthenosphere boundary, and destroy the Gutenberg discontinuity.

Oxidized primary arc magmas: Constraints from Cu/Zr systematics in global arc volcanics

Arc volcanics are more oxidized than mid-ocean ridge basalts (MORB), but it is debated whether this is a mantle feature or a result of magmatic evolution. Copper, a sulfur-loving element, has been used to trace the behavior of redox-sensitive sulfur during mantle melting and infer similar redox states of sub-arc and sub-ridge mantle. Previous studies, however, neglected elevated sulfur contents in the sub-arc mantle, leading to underestimation of oxygen fugacities, and did not recognize systematic Cu variations in arc volcanics. Here, we show that the Cu/Zr ratio is a sensitive indicator that responds to sulfur content, oxygen fugacity, and extent of melting of the mantle. Because of higher mantle S contents, Cu systematics of arc magmas require one log unit higher oxygen fugacities of sub-arc than sub-ridge mantle. Low Cu contents of thick-crusted arc volcanics result from low extents of melting of sulfur-rich mantle, obviating the need for deep crustal sulfide fractionation, with substantial implications for the origin of porphyry-Cu deposits.

On the relative temperatures of Earth's volcanic hotspots and mid-ocean ridges

Volcanic hotspots are thought to be fed by hot, active upwellings from the deep mantle, with excess temperatures (T_ex) ~100° to 300°C higher than those of mid-ocean ridges. However, T_ex estimates are limited in geographical coverage and often inconsistent for individual hotspots. We infer the temperature of oceanic hotspots and ridges simultaneously by converting seismic velocity to temperature. We show that while ~45% of plume-fed hotspots are hot (T_ex ≥ 155°C), ~15% are cold (T_ex ≤ 36°C) and ~40% are not hot enough to actively upwell (50°C ≤ T_ex ≤ 136°C). Hot hotspots have an extremely high helium-3/helium-4 ratio and buoyancy flux, but cold hotspots do not. The latter may originate at upper mantle depths. Alternatively, the deep plumes that feed them may be entrained and cooled by small-scale convection.

Evolution of the Crustal and Upper Mantle Seismic Structure From 0-27 Ma in the Equatorial Atlantic Ocean at 2° 43'S

We present seismic tomographic results from a unique seismic refraction and wide-angle survey along a 600 km long flow-line corridor of oceanic lithosphere ranging in age from 0 to 27 Ma in the equatorial Atlantic Ocean at 2° 43'S. The velocities in the crust near the ridge axis rapidly increase in the first 6 Myr and then change gradually with age. The upper crust (Layer 2) thickness varies between 2 and 2.4 km with an average thickness of 2.2 km and the crustal thickness varies from 5.6 to 6 km along the profile with an average crustal thickness of 5.8 km. At some locations, we observe negative velocity anomalies (∼-0.3 km/s) in the lower crust which could be either due to chemical heterogeneity in gabbroic rocks and/or the effects of fault related deformation zones leading to an increase in porosities up to 1.6% depending on the pore/crack geometry. The existence of a low velocity anomaly beneath the ridge axis suggests the presence of partial melt (∼1.3%) in the lower crust. Upper mantle velocities also remain low (∼7.8 km/s) from ridge axis up to 5 Ma, indicating a high temperature regime associated with mantle melting zone underneath. These results suggest that the evolution of the crust and uppermost mantle at this location occur in the first 10 Ma of its formation and then remains unchanged. Most of the structures in the older crust and upper mantle are fossilized structures and could provide information about past processes at ocean spreading centers.

Global influence of mantle temperature and plate thickness on intraplate volcanism

The thermochemical structure of lithospheric and asthenospheric mantle exert primary controls on surface topography and volcanic activity. Volcanic rock compositions and mantle seismic velocities provide indirect observations of this structure. Here, we compile and analyze a global database of the distribution and composition of Neogene-Quaternary intraplate volcanic rocks. By integrating this database with seismic tomographic models, we show that intraplate volcanism is concentrated in regions characterized by slow upper mantle shear-wave velocities and by thin lithosphere (i.e. <100 km). We observe a negative correlation between shear-wave velocities at depths of 125–175 km and melt fractions inferred from volcanic rock compositions. Furthermore, mantle temperature and lithospheric thickness estimates obtained by geochemical modeling broadly agree with values determined from tomographic models that have been converted into temperature. Intraplate volcanism often occurs in regions where uplifted (but undeformed) marine sedimentary rocks are exposed. Regional elevation of these rocks can be generated by a combination of hotter asthenosphere and lithospheric thinning. Therefore, the distribution and composition of intraplate volcanic rocks through geologic time will help to probe past mantle conditions and surface processes.

Here, the authors compile a global geochemical database of Neogene-Quaternary intraplate volcanism. By comparing the distribution and composition of these rocks with tomographic models they show that intraplate volcanism can be used to constrain upper-mantle structure at the time of eruption.

A thin mantle transition zone beneath the equatorial Mid-Atlantic Ridge

The location and degree of material transfer between the upper and lower mantle are key to the Earth's thermal and chemical evolution. Sinking slabs and rising plumes are generally accepted as locations of transfer^1,2, whereas mid-ocean ridges are not typically assumed to have a role³. However, tight constraints from in situ measurements at ridges have proved to be challenging. Here we use receiver functions that reveal the conversion of primary to secondary seismic waves to image the discontinuities that bound the mantle transition zone, using ocean bottom seismic data from the equatorial Mid-Atlantic Ridge. Our images show that the seismic discontinuity at depths of about 660 kilometres is broadly uplifted by 10 ± 4 kilometres over a swath about 600 kilometres wide and that the 410-kilometre discontinuity is depressed by 5 ± 4 kilometres. This thinning of the mantle transition zone is coincident with slow shear-wave velocities in the mantle, from global seismic tomography^4-7. In addition, seismic velocities in the mantle transition zone beneath the Mid-Atlantic Ridge are on average slower than those beneath older Atlantic Ocean seafloor. The observations imply material transfer from the lower to the upper mantle-either continuous or punctuated-that is linked to the Mid-Atlantic Ridge. Given the length and longevity of the mid-ocean ridge system, this implies that whole-mantle convection may be more prevalent than previously thought, with ridge upwellings having a role in counterbalancing slab downwellings.

Seismic Crustal Structure and Morphotectonic Features Associated With the Chain Fracture Zone and Their Role in the Evolution of the Equatorial Atlantic Region

Oceanic transform faults and fracture zones (FZs) represent major bathymetric features that keep the records of past and present strike-slip motion along conservative plate boundaries. Although they play an important role in ridge segmentation and evolution of the lithosphere, their structural characteristics, and their variation in space and time, are poorly understood. To address some of the unknowns, we conducted interdisciplinary geophysical studies in the equatorial Atlantic Ocean, the region where some of the most prominent transform discontinuities have been developing. Here we present the results of the data analysis in the vicinity of the Chain FZ, on the South American Plate. The crustal structure across the Chain FZ, at the contact between ∼10 and 24 Ma oceanic lithosphere, is sampled along seismic reflection and refraction profiles. We observe that the crustal thickness within and across the Chain FZ ranges from ∼4.6-5.9 km, which compares with the observations reported for slow-slipping transform discontinuities globally. We attribute this presence of close to normal oceanic crustal thickness within FZs to the mechanism of lateral dike propagation, previously considered to be valid only in fast-slipping environments. Furthermore, the combination of our results with other data sets enabled us to extend the observations to morphotectonic characteristics on a regional scale. Our broader view suggests that the formation of the transverse ridge is closely associated with a global plate reorientation that was also responsible for the propagation and for shaping lower-order Mid-Atlantic Ridge segmentation around the equator.

Recycled arc mantle recovered from the Mid-Atlantic Ridge

Plate tectonics and mantle dynamics necessitate mantle recycling throughout Earth’s history, yet direct geochemical evidence for mantle reprocessing remains elusive. Here we present evidence of recycled supra-subduction zone mantle wedge peridotite dredged from the Mid-Atlantic Ridge near 16°30′N. Peridotite trace-element characteristics are inconsistent with fractional anhydrous melting typically associated with a mid-ocean ridge setting. Instead, the samples are best explained by hydrous flux melting which changed the melting reactions such that clinopyroxene was not exhausted at high degrees of melting and was retained in the residuum. Based on along-axis ridge depth variations, this buoyant refractory arc mantle is likely compensated at depth by denser, likely garnet-rich, lithologies within the mantle column. Our results suggest that highly refractory arc mantle relicts are entrained in the upper mantle and may constitute >60% of the upper mantle by volume. These highly refractory mantle domains, which contribute little to mantle melting, are under-represented in compilations of mantle composition that rely on inverted basalt compositions alone.

Plate tectonics necessitates mantle recycling throughout Earth’s history, yet direct geochemical evidence for mantle reprocessing remains elusive. Here, the authors present evidence of recycled supra-subduction zone mantle wedge peridotite dredged from the Mid-Atlantic Ridge near 16°30′N.

High H2O Content in Pyroxenes of Residual Mantle Peridotites at a Mid Atlantic Ridge Segment

Global correlations of mid-ocean-ridges basalt chemistry, axial depth and crustal thickness have been ascribed to mantle temperature variations affecting degree of melting. However, mantle H₂O content and elemental composition may also play a role. How H₂O is distributed in the oceanic upper mantle remains poorly constrained. We tackled this problem by determining the H₂O content of orthopyroxenes (opx) and clinopyroxenes (cpx) of peridotites from a continuous lithospheric section created during 26 Ma at a 11°N Mid-Atlantic Ridge segment, and exposed along the Vema Transform. The H₂O content of opx ranges from 119 ppm to 383 ppm; that of cpx from 407 ppm to 1072 ppm. We found anomalous H₂O-enriched peridotites with their H₂O content not correlating inversely with their degree of melting, although H₂O is assumed to be incompatible during melting. Inverse correlation of H₂O with Ce, another highly incompatible component, suggests post-melting H₂O enrichment. We attribute a major role to post-melting temperature-dependent diffusion of hydrogen occurring above the melting region, where water-rich melt flows faster than residual peridotites through dunitic conduits cross-cutting the uprising mantle. Accordingly, estimates of the H₂O content of the MORB mantle source based on H₂O in abyssal peridotites can be affected by strong uncertainties.

Northward drift of the Azores plume in the Earth's mantle

Mantle plume fixity has long been a cornerstone assumption to reconstruct past tectonic plate motions. However, precise geochronological and paleomagnetic data along Pacific continuous hotspot tracks have revealed substantial drift of the Hawaiian plume. The question remains for evidence of drift for other mantle plumes. Here, we use plume-derived basalts from the Mid-Atlantic ridge to confirm that the upper-mantle thermal anomaly associated with the Azores plume is asymmetric, spreading over ~2,000 km southwards and ~600 km northwards. Using for the first time a 3D-spherical mantle convection where plumes, ridges and plates interact in a fully dynamic way, we suggest that the extent, shape and asymmetry of this anomaly is a consequence of the Azores plume moving northwards by 1-2 cm/yr during the past 85 Ma, independently from other Atlantic plumes. Our findings suggest redefining the Azores hotspot track and open the way for identifying how plumes drift within the mantle.

The causes of spatiotemporal variations in erupted fluxes and compositions along a volcanic arc

Decades of study on volcanic arcs have provided insight into the overarching processes that control magmatism, and how these processes manifest at individual volcanoes. However, the causes of ubiquitous and dramatic intra-arc variations in volcanic flux and composition remain largely unresolved. Investigating such arc-scale issues requires greater quantitative comparison of geophysical and geochemical data, linked through sets of common intensive variables. To work towards these goals, we use observed lava compositions to estimate the heat budget associated with Quaternary volcanism in the Cascades Arc and compare this to the heat required to produce the observed geophysical properties of the crust. Here we show that along-strike volcanic variability in the Quaternary Cascades Arc is primarily related to variations in the flux of basalt into the crust, rather than variations in their crustal storage history. This approach shows promise for studying other large-scale frontier geologic problems in volcanic arcs.

The primary causes of dramatic variations in volcanic flux and composition along strike in subduction zones remain largely unknown. Here we use a promising new approach to show that along-strike volcanic variability in the Quaternary Cascades Arc is primarily due to variations in the flux of basalt into the base of the crust, rather than crustal magma storage.

Petit-spot as definitive evidence for partial melting in the asthenosphere caused by CO2

The deep carbon cycle plays an important role on the chemical differentiation and physical properties of the Earth's mantle. Especially in the asthenosphere, seismic low-velocity and high electrical conductivity due to carbon dioxide (CO₂)-induced partial melting are expected but not directly observed. Here we discuss the experimental results relevant to the genesis of primitive CO₂-rich alkali magma forming petit-spot volcanoes at the deformation front of the outer rise of the northwestern Pacific plate. The results suggest that primitive melt last equilibrated with depleted peridotite at 1.8–2.1 GPa and 1,280–1,290 °C. Although the equilibration pressure corresponds to the pressure of the lower lithosphere, by considering an equilibration temperature higher than the solidus in the volatile–peridotite system along with the temperature of the lower lithosphere, we conclude that CO₂-rich silicate melt is always produced in the asthenosphere. The melt subsequently ascends into and equilibrates with the lower lithosphere before eruption.

The genesis of primitive CO₂-rich alkali magma at petit-spot volcanoes is still poorly understood. Here, the authors perform high pressure and temperature experiments to show that the CO₂-rich melt is derived from the asthenosphere.

Record of massive upwellings from the Pacific large low shear velocity province

Large igneous provinces, as the surface expression of deep mantle processes, play a key role in the evolution of the planet. Here we analyse the geochemical record and timing of the Pacific Ocean Large Igneous Provinces and preserved accreted terranes to reconstruct the history of pulses of mantle plume upwellings and their relation with a deep-rooted source like the Pacific large low-shear velocity Province during the Mid-Jurassic to Upper Cretaceous. Petrological modelling and geochemical data suggest the need of interaction between these deep-rooted upwellings and mid-ocean ridges in pulses separated by ∼10-20 Ma, to generate the massive volumes of melt preserved today as oceanic plateaus. These pulses impacted the marine biota resulting in episodes of anoxia and mass extinctions shortly after their eruption.