A unified continental thickness from seismology and diamonds suggests a melt-defined plate

Discovery of flat seismic reflections in the mantle beneath the young Juan de Fuca Plate

Crustal properties of young oceanic lithosphere have been examined extensively, but the nature of the mantle lithosphere underneath remains elusive. Using a novel wide-angle seismic imaging technique, here we show the presence of two sub-horizontal reflections at ∼11 and ∼14.5 km below the seafloor over the 0.51-2.67 Ma old Juan de Fuca Plate. We find that the observed reflectors originate from 300-600-m-thick layers, with an ∼7-8% drop in P-wave velocity. They could be explained either by the presence of partially molten sills or frozen gabbroic sills. If partially molten, the shallower sill would define the base of a thin lithosphere with the constant thickness (11 km), requiring the presence of a mantle thermal anomaly extending up to 2.67 Ma. In contrast, if these reflections were frozen melt sills, they would imply the presence of thick young oceanic lithosphere (20-25 km), and extremely heterogeneous upper mantle.

Imaging Lithospheric Discontinuities Beneath the Northern East African Rift Using S-to-P Receiver Functions

Imaging the lithosphere is key to understand mechanisms of extension as rifting progresses. Continental rifting results in a combination of mechanical stretching and thinning of the lithosphere, decompression upwelling, heating, sometimes partial melting of the asthenosphere, and potentially partial melting of the mantle lithosphere. The northern East African Rift system is an ideal locale to study these processes as it exposes the transition from tectonically active continental rifting to incipient seafloor spreading. Here we use S-to-P receiver functions to image the lithospheric structure beneath the northernmost East African Rift system where it forms a triple junction between the Main Ethiopian rift, the Red Sea rift, and the Gulf of Aden rift. We image the Moho at 31 ± 6 km beneath the Ethiopian plateau. The crust is 28 ± 3 km thick beneath the Main Ethiopian rift and thins to 23 ± 2 km in northern Afar. We identify a negative phase, a velocity decrease with depth, at 67 ± 3 km depth beneath the Ethiopian plateau, likely associated with the lithosphere-asthenosphere boundary (LAB), and a lack of a LAB phase beneath the rift. Using observations and waveform modeling, we show that the LAB phase beneath the plateau is likely defined by a small amount of partial melt. The lack of a LAB phase beneath the rift suggests melt percolation through the base of the lithosphere beneath the northernmost East African Rift system.

Seismic Imaging of the North American Midcontinent Rift Using S-to-P Receiver Functions

North America's ~1.1-Ga failed Midcontinent Rift (MCR) is a striking feature of gravity and magnetic anomaly maps across the continent. However, how the rift affected the underlying lithosphere is not well understood. With data from the Superior Province Rifting Earthscope Experiment and the USArray Transportable Array, we constrain three-dimensional seismic velocity discontinuity structure in the lithosphere beneath the southwestward arm of the MCR using S-to-P receiver functions. We image a velocity increase with depth associated with the Moho at depths of 33-40 ± 4 km, generally deepening toward the east. The Moho amplitude decreases beneath the rift axis in Minnesota and Wisconsin, where the velocity gradient is more gradual, possibly due to crustal underplating. We see hints of a deeper velocity increase at 61 ± 4-km depth that may be the base of underplating. Beneath the rift axis further south in Iowa, we image two distinct positive phases at 34-39 ± 4 and 62-65 ± 4 km likely related to an altered Moho and an underplated layer. We image velocity decreases with depth at depths of 90-190 ± 7 km in some locations that do not geographically correlate with the rift. These include a discontinuity at depths of 90-120 ± 7 km with a northerly dip in the south that abruptly deepens to 150-190 ± 7 km across the Spirit Lake Tectonic Zone provincial suture. The negative phases may represent a patchy, frozen-in midlithosphere discontinuity feature that likely predates the MCR and/or be related to lithospheric thickness.

Predictions and Observations for the Oceanic Lithosphere From S-to-P Receiver Functions and SS Precursors

The ocean lithosphere is classically described by the thermal half-space cooling (HSC) or the plate models, both characterized by a gradual transition to the asthenosphere beneath. Scattered waves find sharp seismic discontinuities beneath the oceans, possibly from the base of the plate. Active source studies suggest sharp discontinuities from a melt channel. We calculate synthetic S-to-P receiver functions and SS precursors for the HSC and plate models and also for channels. We find that the HSC and plate model velocity gradients are too gradual to create interpretable scattered waves from the base of the plate. Subtle phases are predicted to follow a similar trend as observations, flattening at older ages. Therefore, the seismic discontinuities are probably caused by a thermally controlled process that can also explain their amplitude, such as melting. Melt may coalesce in channels, although channels >10 km thick should be resolvable by scattered wave imaging.

The growth of lithospheric diamonds

Natural diamonds contain mineral and fluid inclusions that record diamond growth conditions. Replicating the growth of inclusion-bearing diamonds in a laboratory is therefore a novel diagnostic tool to constrain the conditions of diamond formation in Earth's lithosphere. By determining the carbon isotopic fractionation during diamond growth in fluids or melts, our laboratory experiments revealed that lithospheric monocrystalline and fibrous and coated diamonds grow similarly from redox reactions at isotopic equilibrium in water and carbonate-rich fluids or melts, and not from native carbon. These new results explain why most of the lithospheric diamonds are characterized by a common carbon isotopic fingerprint, inherited from their common parent fluids and not from the mantle assemblage.

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.