Highly variable magmatic accretion at the ultraslow-spreading Gakkel Ridge

Crustal accretion at mid-ocean ridges governs the creation and evolution of the oceanic lithosphere. Generally accepted models1-4 of passive mantle upwelling and melting predict notably decreased crustal thickness at a spreading rate of less than 20 mm year-1. We conducted the first, to our knowledge, high-resolution ocean-bottom seismometer (OBS) experiment at the Gakkel Ridge in the Arctic Ocean and imaged the crustal structure of the slowest-spreading ridge on the Earth. Unexpectedly, we find that crustal thickness ranges between 3.3 km and 8.9 km along the ridge axis and it increased from about 4.5 km to about 7.5 km over the past 5 Myr in an across-axis profile. The highly variable crustal thickness and relatively large average value does not align with the prediction of passive mantle upwelling models. Instead, it can be explained by a model of buoyant active mantle flow driven by thermal and compositional density changes owing to melt extraction. The influence of active versus passive upwelling is predicted to increase with decreasing spreading rate. The process of active mantle upwelling is anticipated to be primarily influenced by mantle temperature and composition. This implies that the observed variability in crustal accretion, which includes notably varied crustal thickness, is probably an inherent characteristic of ultraslow-spreading ridges.

A role for subducting clays in the water transportation into the Earth's lower mantle

Subducting sedimentary layer typically contains water and hydrated clay minerals. The stability of clay minerals under such hydrous subduction environment would therefore constraint the lithology and physical properties of the subducting slab interface. Here we show that pyrophyllite (Al2Si4O10(OH)2), one of the representative clay minerals in the alumina-silica-water (Al2O3-SiO2-H2O, ASH) system, breakdowns to contain further hydrated minerals, gibbsite (Al(OH)3) and diaspore (AlO(OH)), when subducts along a water-saturated cold subduction geotherm. Such a hydration breakdown occurs at a depth of ~135 km to uptake water by ~1.8 wt%. Subsequently, dehydration breakdown occurs at ~185 km depth to release back the same amount of water, after which the net crystalline water content is preserved down to ~660 km depth, delivering a net amount of ~5.0 wt% H2O in a phase assemblage containing δ-AlOOH and phase Egg (AlSiO3(OH)). Our results thus demonstrate the importance of subducting clays to account the delivery of ~22% of water down to the lower mantle.

Tracking carbon from subduction to outgassing along the Aleutian-Alaska Volcanic Arc

Subduction transports volatiles between Earth's mantle, crust, and atmosphere, ultimately creating a habitable Earth. We use isotopes to track carbon from subduction to outgassing along the Aleutian-Alaska Arc. We find substantial along-strike variations in the isotopic composition of volcanic gases, explained by different recycling efficiencies of subducting carbon to the atmosphere via arc volcanism and modulated by subduction character. Fast and cool subduction facilitates recycling of ~43 to 61% sediment-derived organic carbon to the atmosphere through degassing of central Aleutian volcanoes, while slow and warm subduction favors forearc sediment removal, leading to recycling of ~6 to 9% altered oceanic crust carbon to the atmosphere through degassing of western Aleutian volcanoes. These results indicate that less carbon is returned to the deep mantle than previously thought and that subducting organic carbon is not a reliable atmospheric carbon sink over subduction time scales.

Ignition of the southern Atlantic seafloor spreading machine without hot-mantle booster

The source of massive magma production at volcanic rifted margins remains strongly disputed since the first observations of thick lava piles in the 1980s. However, volumes of extruded and intruded melt products within rifted continental crust are still not accurately resolved using geophysical methods. Here we investigate the magma budget alongside the South Atlantic margins, at the onset of seafloor spreading, using high-quality seismic reflection profiles to accurately estimate the oceanic crustal thickness. We show that, along ~ 75% of the length of the Early-Cretaceous initial spreading centre, the crustal thickness is similar to regular oceanic thickness with an age > 100 Ma away from hot spots. Thus, most of the southernmost Atlantic Ocean opened without anomalously hot mantle, high magma supply being restricted to the Walvis Ridge area. We suggest that alternative explanations other than a hotter mantle should be favoured to explain the thick magmatic layer of seaward dipping reflectors landward of the initial mid-oceanic ridge.

Seismic evidence for uniform crustal accretion along slow-spreading ridges in the equatorial Atlantic Ocean

The crustal accretion along mid-ocean ridges is known to be spreading-rate dependent. Along fast-spreading ridges, two-dimensional sheet-like mantle upwelling creates relatively uniform crust. In contrast, the crust formed along slow-spreading ridges shows large along-axis thickness variations with thicker crust at segment centres, which is hypothesised to be due a three-dimensional plume-like mantle upwelling or due to focused melt migration to segment centres. Using wide-angle seismic data acquired from the equatorial Atlantic Ocean, here we show that the crustal thickness is nearly uniform (~5.5 km) across five crustal segments for crust formed at the slow-spreading Mid-Atlantic Ridge with age varying from 8 to 70 Ma. The crustal velocities indicate that this crust is predominantly of magmatic origin. We suggest that this uniform magmatic crustal accretion is due to a two-dimensional sheet-like mantle upwelling facilitated by the long-offset transform faults in the equatorial Atlantic region and the presence of a high concentration of volatiles in the primitive melt in the mantle.

Tectonic degassing drove global temperature trends since 20 Ma

The Miocene Climatic Optimum (MCO) from ~17 to 14 million years ago (Ma) represents an enigmatic reversal in Cenozoic cooling. A synthesis of marine paleotemperature records shows that the MCO was a local maximum in global sea surface temperature superimposed on a period from at least 19 Ma to 10 Ma, during which global temperatures were on the order of 10°C warmer than at present. Our high-resolution global reconstruction of ocean crustal production, a proxy for tectonic degassing of carbon, suggests that crustal production rates were ~35% higher than modern rates until ~14 Ma, when production began to decline steeply along with global temperatures. The magnitude and timing of the inferred changes in tectonic degassing can account for the majority of long-term ice sheet and global temperature evolution since 20 Ma.

3-D data of thermal regime, water content, and slab dehydration in Alaska

The data include the 3-D temperature field (degrees Celsius), water content (wt%), dehydration rate (wt%/km), and subduction velocity field (cm/yr) of the subducting plate, as well as the coastline and volcano distribution in Alaska. The data of the model region have dimensions of 800 × 1600 × 400 km (length × width × depth). The geometry of the subducted plate is well constrained by Slab2.0, and the plate ages are provided by EarthByte. The subduction velocities inside a prescribed 3-D constrained volume of the oceanic lithosphere are given based on the kinematic plate subduction modeling method and the MORVEL plate motion data. The observation of surface heat flow and Curie point depths are used to constrain the model thermal regime. The geophysical calculation is ensured after the subduction thermal regime reaches a steady state. Data are deposited in the TPDC repository, which has granted a persistent identifier https://data.tpdc.ac.cn/en/disallow/8b266d22-fea7-4259-9a5f-8ac0bd9e7869/. Data include (1) paraview_eq_USGS.vtk (earthquake catalog by IRIS, 2000-2010, Trabant et al., 2012), (2) paraview_slab.vtk (3-D thermal regime, slab water content and slab dehydration), (3) paraview_volcano.vtk (global volcanoes at NCEI, Siebert et al., 2010), and (4) paraview_map.vtk (coastline, GMT).

Separation of Gagua Rise from Great Benham Rise in the West Philippine Basin during the Middle Eocene

The West Philippine Basin (WPB) has started opening at ~ 58 Ma and ceased spreading at ~ 33 Ma, developing a fast spreading (~ 44 mm/yr half-spreading rate) magmatic episode between 58 and 41 Ma and the second amagmatic episode between 41 and 33 Ma. The occurrence of the first stage of spreading is closely related to the Oki-Daito mantle plume and related Benham Rise (BR) and Urdaneta Plateau (UP) activity. To the east of the Luzon–Okinawa Fracture Zone (LOFZ), BR was the most active volcanism from 48 to 41 Ma. The geomagnetic ages on both sides of the LOFZ have been determined; however, their causal relationship and evolution in the WPB remain unclear. In this study, we performed integrated analyses of multichannel seismic data and swath bathymetry data for the area to the west of the LOFZ. To the west of the LOFZ, the Gagua Rise (GR), is identified by a high residual free-air gravity anomaly, volcanic seamount chains and an overlapping spreading center. The GR is located at magnetic isochrons C20/C22 (50 to 44 Ma) and shows a thick oceanic crust of at least 12.7 km. We first propose an oceanic plateau named Great Benham Rise (GBR) which includes GR, UP and BR. We infer that the GR was a portion of the GBR since ~ 49 Ma and was separated from the GBR at ~ 41 Ma by the right-lateral LOFZ motion. Later, the relict GBR magmatism only continued in the area to the east of the LOFZ. Overall, the GBR dominates the spreading history of the WPB.

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.

The early drift of the Indian plate

Plate kinematic models propose that India and Sri Lanka (INDSRI) separated from Antarctica by extremely slow seafloor spreading that started in early Cretaceous times, and that a long-distance ridge jump left a continental fragment stranded off the Antarctic margin under the Southern Kerguelen Plateau ^1-3. Here, we present newly acquired magnetic and deep wide-angle seismic data that require a fundamental re-evaluation of these concepts. The new data clearly define the onset of oceanic crust in the Enderby Basin and off southern Sri Lanka, and date its formation with unprecedented confidence. The revised timing indicates that India and Sri Lanka detached from Antarctica earlier in the east than in the west. Furthermore, no compelling evidence for an extinct spreading axis is found in the Enderby Basin. A refined plate motion model indicates that India and Sri Lanka departed from Antarctica without major rift jumps, but by the action of three spreading ridges with different timings and velocities that must have been accommodated by significant intracontinental deformation.

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.

Caribbean plate tilted and actively dragged eastwards by low-viscosity asthenospheric flow

The importance of a low-viscosity asthenosphere underlying mobile plates has been highlighted since the earliest days of the plate tectonics revolution. However, absolute asthenospheric viscosities are still poorly constrained, with estimates spanning up to 3 orders of magnitude. Here we follow a new approach using analytic solutions for Poiseuille-Couette channel flow to compute asthenospheric viscosities under the Caribbean. We estimate Caribbean dynamic topography and the associated pressure gradient, which, combined with flow velocities estimated from geologic markers and tomographic structure, yield our best-estimate asthenospheric viscosity of (3.0 ± 1.5)*10¹⁸ Pa s. This value is consistent with independent estimates for non-cratonic and oceanic regions, and challenges the hypothesis that higher-viscosity asthenosphere inferred from postglacial rebound is globally-representative. The active flow driven by Galapagos plume overpressure shown here contradicts the traditional view that the asthenosphere is only a passive lubricating layer for Earth's tectonic plates.

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.

Oceanic mantle reflections in deep seismic profiles offshore Sumatra are faults or fakes

In the late 90’s, some faults identified within oceanic crust were demonstrated to be artifacts arising from out-of-plane scattering along linear sediment-buried fault scarps. Symmetrical mantle reflections observed southwest northern Sumatra on seismic reflection profiles have been identified as faults cutting through the upper mantle down to unprecedented depths reaching ~45 km. Seawater being conveyed along sub-vertical re-activated fracture zones (FZs) to the upper mantle, the mantle portions of FZs are serpentinized and act as mirrors for seismic rays. We suggest that the mantle features are not faults but artifacts resulting from out-of-plane reflections on these mirrors. Two perpendicular seismic profiles crossing the same FZ display two dipping features down to 30 km, which cannot be explained as faults from recent tectonic and structural constraints but merely as out-of-plane reflections on this FZ. This result confirms that most of mantle reflections observed southwest northern Sumatra are fakes rather than faults.

Melt movement through the Icelandic crust

We use both seismology and geobarometry to investigate the movement of melt through the volcanic crust of Iceland. We have captured melt in the act of moving within or through a series of sills ranging from the upper mantle to the shallow crust by the clusters of small earthquakes it produces as it forces its way upward. The melt is injected not just beneath the central volcanoes, but also at discrete locations along the rift zones and above the centre of the underlying mantle plume. We suggest that the high strain rates required to produce seismicity at depths of 10-25 km in a normally ductile part of the Icelandic crust are linked to the exsolution of carbon dioxide from the basaltic melts. The seismicity and geobarometry provide complementary information on the way that the melt moves through the crust, stalling and fractionating, and often freezing in one or more melt lenses on its way upwards: the seismicity shows what is happening instantaneously today, while the geobarometry gives constraints averaged over longer time scales on the depths of residence in the crust of melts prior to their eruption. This article is part of the Theo Murphy meeting issue 'Magma reservoir architecture and dynamics'.

LIP formation and protracted lower mantle upwelling induced by rifting and delamination

Large Igneous Provinces (LIPs) are commonly attributed to mantle plumes, hot upwellings from the deep lower mantle, apparently unrelated to plate motions. However, LIPs often form in association with rifting and breakup. Using numerical modelling, we introduce a novel idea that explains plume-like mantle upwelling by plate tectonic processes. Our model indicates that rifting-induced delamination of orogenic lithosphere can perturb the thermochemical mantle stratification and induce lower mantle upwelling which causes syn-rift LIP formation followed by protracted and enhanced mid ocean ridge basalt (MORB) generation. Our model provides an explanation for the geographical correlation between the Caledonian suture, the North Atlantic Igneous Province (NAIP) and present-day Icelandic magmatism.

Scattered wave imaging of the oceanic plate in Cascadia

Fifty years after plate tectonic theory was developed, the defining mechanism of the plate is still widely debated. The relatively short, simple history of young ocean lithosphere makes it an ideal place to determine the property that defines a plate, yet the remoteness and harshness of the seafloor have made precise imaging challenging. We use S-to-P receiver functions to image discontinuities beneath newly formed lithosphere at the Juan de Fuca and Gorda Ridges. We image a strong negative discontinuity at the base of the plate increasing from 20 to 45 km depth beneath the 0- to 10-million-year-old seafloor and a positive discontinuity at the onset of melting at 90 to 130 km depth. Comparison with geodynamic models and experimental constraints indicates that the observed discontinuities cannot easily be reconciled with subsolidus mechanisms. Instead, partial melt may be required, which would decrease mantle viscosity and define the young oceanic plate.

Evidence for frozen melts in the mid-lithosphere detected from active-source seismic data

The interactions of the lithospheric plates that form the Earth's outer shell provide much of the evidentiary basis for modern plate tectonic theory. Seismic discontinuities in the lithosphere arising from mantle convection and plate motion provide constraints on the physical and chemical properties of the mantle that contribute to the processes of formation and evolution of tectonic plates. Seismological studies during the past two decades have detected seismic discontinuities within the oceanic lithosphere in addition to that at the lithosphere-asthenosphere boundary (LAB). However, the depth, distribution, and physical properties of these discontinuities are not well constrained, which makes it difficult to use seismological data to examine their origin. Here we present new active-source seismic data acquired along a 1,130 km profile across an old Pacific plate (148-128 Ma) that show oceanic mid-lithosphere discontinuities (oceanic MLDs) distributed 37-59 km below the seafloor. The presence of the oceanic MLDs suggests that frozen melts that accumulated at past LABs have been preserved as low-velocity layers within the current mature lithosphere. These observations show that long-offset, high-frequency, active-source seismic data can be used to image mid-lithospheric structure, which is fundamental to understanding the formation and evolution of tectonic plates.

Moho and basement depth in the NE Atlantic Ocean based on seismic refraction data and receiver functions

Seismic refraction data and results from receiver functions were used to compile the depth to the basement and Moho in the NE Atlantic Ocean. For interpolation between the unevenly spaced data points, the kriging technique was used. Free-air gravity data were used as constraints in the kriging process for the basement. That way, structures with little or no seismic coverage are still presented on the basement map, in particular the basins off East Greenland. The rift basins off NW Europe are mapped as a continuous zone with basement depths of between 5 and 15 km. Maximum basement depths off NE Greenland are 8 km, but these are probably underestimated. Plate reconstructions for Chron C24 (c. 54 Ma) suggest that the poorly known Ammassalik Basin off SE Greenland may correlate with the northern termination of the Hatton Basin at the conjugate margin. The most prominent feature on the Moho map is the Greenland–Iceland–Faroe Ridge, with Moho depths >28 km. Crustal thickness is compiled from the Moho and basement depths. The oceanic crust displays an increased thickness close to the volcanic margins affected by the Iceland plume.