Subcontinental lithospheric mantle discontinuities beneath the Eastern Himalayan Plate Boundary System, NE India

We use P-wave receiver function (P-RF) analysis and joint inversion with Rayleigh wave group velocity dispersion data to model the shear wave velocity (Vs) structure of subcontinental lithospheric mantle (SCLM) discontinuities beneath northeast (NE) India. The most prominent SCLM discontinuity is the Hales Discontinuity (H-D) observed beneath the Eastern Himalayan Foreland Basin (Brahmaputra Valley) and Shillong Plateau. The P-to-SV converted phase from the H-D (Phs) is a positive amplitude arrival at ∼10–12 s and has positive moveout with increasing ray-parameter. From joint inversion, the H-D is modelled at a depth range of 90–106 km, with 11–12 per cent Vs increase beneath the Brahmaputra Valley. Beneath the Shillong Plateau the H-D is at a depth range of 86–99 km, with 6–10 per cent Vs increase. An intralithospheric discontinuity (ILD) has been identified in the Shillong Plateau station P-RFs, as a positive amplitude PILDs phase, arriving at 8–8.5 s. This is modelled at a depth range of 66–75 km with Vs increase of 2–9 per cent. We construct 2-D profiles of depth-migrated common-conversion-point stack of P-RFs to distinguish the SCLM discontinuity arrivals from crustal phases. 3-D spline-interpolated surface of the H-D has been constructed to visualize its lateral variations. We use xenolith data from the Dharwar Craton, which has similar geological age, petrology and seismic structure as the Shillong Plateau, to petrologically model the SCLM H-D and ILD Vs structure in NE-India. From the calculated Vs structure we conjecture that the H-D is a petrological boundary between mantle peridotite and kyanite-eclogite, with its origin as metamorphosed paleosubducted oceanic slab, similar to other global observations. We further speculate that the shallower ILD could be formed as a contact between frozen asthenosphere-derived metasomatic melts within the SCLM.

Seismological evidence for the earliest global subduction network at 2 Ga ago

The earliest evidence for subduction, which could have been localized, does not signify when plate tectonics became a global phenomenon. To test the antiquity of global subduction, we investigated Paleoproterozoic time, for which seismic evidence is available from multiple continents. We used a new high-density seismic array in North China to image the crustal structure that exhibits a dipping Moho bearing close resemblance to that of the modern Himalaya. The relict collisional zone is Paleoproterozoic in age and implies subduction operating at least as early as ~2 billion years (Ga) ago. Seismic evidence of subduction from six continents at this age is interpreted as the oldest evidence of global plate tectonics. The sutures identified can be linked in a plate network that resulted in the assembly of Nuna, likely Earth's first supercontinent. Global subduction by ~2 Ga ago can explain why secular planetary cooling was not appreciable until Proterozoic time.

Seismic evidence for megathrust fault-valve behavior during episodic tremor and slip

Fault slip behavior during episodic tremor and slow slip (ETS) events, which occur at the deep extension of subduction zone megathrust faults, is believed to be related to cyclic fluid processes that necessitate fluctuations in pore-fluid pressures. In most subduction zones, a layer of anomalously low seismic wave velocities [low-velocity layer (LVL)] is observed in the vicinity of ETS and suggests high pore-fluid pressures that weaken the megathrust. Using repeated seismic scattering observations in the Cascadia subduction zone, we observe a change in the seismic velocity associated with the LVL after ETS events, which we interpret as a response to fluctuations in pore-fluid pressure. These results provide direct evidence of megathrust fault-valve processes during ETS.

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.

Constraints on the anisotropic contributions to velocity discontinuities at ∼60 km depth beneath the Pacific

Strong, sharp, negative seismic discontinuities, velocity decreases with depth, are observed beneath the Pacific seafloor at ∼60 km depth. It has been suggested that these are caused by an increase in radial anisotropy with depth, which occurs in global surface wave models. Here we test this hypothesis in two ways. We evaluate whether an increase in surface wave radial anisotropy with depth is robust with synthetic resolution tests. We do this by fitting an example surface wave data set near the East Pacific Rise. We also estimate the apparent isotropic seismic velocity discontinuities that could be caused by changes in radial anisotropy in S-to-P and P-to-S receiver functions and SS precursors using synthetic seismograms. We test one model where radial anisotropy is caused by olivine alignment and one model where it is caused by compositional layering. The result of our surface wave inversion suggests strong shallow azimuthal anisotropy beneath 0-10 Ma seafloor, which would also have a radial anisotropy signature. An increase in radial anisotropy with depth at 60 km depth is not well-resolved in surface wave models, and could be artificially observed. Shallow isotropy underlain by strong radial anisotropy could explain moderate apparent velocity drops (<6%) in SS precursor imaging, but not receiver functions. The effect is diminished if strong anisotropy also exists at 0-60 km depth as suggested by surface waves. Overall, an increase in radial anisotropy with depth may not exist at 60 km beneath the oceans and does not explain the scattered wave observations.

Lithospheric layering in the North American craton

How cratons-extremely stable continental areas of the Earth's crust-formed and remained largely unchanged for more than 2,500 million years is much debated. Recent studies of seismic-wave receiver function data have detected a structural boundary under continental cratons at depths too shallow to be consistent with the lithosphere-asthenosphere boundary, as inferred from seismic tomography and other geophysical studies. Here we show that changes in the direction of azimuthal anisotropy with depth reveal the presence of two distinct lithospheric layers throughout the stable part of the North American continent. The top layer is thick ( approximately 150 km) under the Archaean core and tapers out on the surrounding Palaeozoic borders. Its thickness variations follow those of a highly depleted layer inferred from thermo-barometric analysis of xenoliths. The lithosphere-asthenosphere boundary is relatively flat (ranging from 180 to 240 km in depth), in agreement with the presence of a thermal conductive root that subsequently formed around the depleted chemical layer. Our findings tie together seismological, geochemical and geodynamical studies of the cratonic lithosphere in North America. They also suggest that the horizon detected in receiver function studies probably corresponds to the sharp mid-lithospheric boundary rather than to the more gradual lithosphere-asthenosphere boundary.

A sharp lithosphere–asthenosphere boundary imaged beneath eastern North America

Plate tectonic theory hinges on the concept of a relatively rigid lithosphere moving over a weaker asthenosphere, yet the nature of the lithosphere-asthenosphere boundary remains poorly understood. The gradient in seismic velocity that occurs at this boundary is central to constraining the physical and chemical properties that create differences in mechanical strength between the two layers. For example, if the lithosphere is simply a thermal boundary layer that is more rigid owing to colder temperatures, mantle flow models indicate that the velocity gradient at its base would occur over tens of kilometres. In contrast, if the asthenosphere is weak owing to volatile enrichment or the presence of partial melt, the lithosphere-asthenosphere boundary could occur over a much smaller depth range. Here we use converted seismic phases in eastern North America to image a very sharp seismic velocity gradient at the base of the lithosphere-a 3-11 per cent drop in shear-wave velocity over a depth range of 11 km or less at 90-110 km depth. Such a strong, sharp boundary cannot be reconciled with a purely thermal gradient, but could be explained by an asthenosphere that contains a few per cent partial melt or that is enriched in volatiles relative to the lithosphere.

Seismic anisotropy: tracing plate dynamics in the mantle

Elastic anisotropy is present where the speed of a seismic wave depends on its direction. In Earth's mantle, elastic anisotropy is induced by minerals that are preferentially oriented in a directional flow or deformation. Earthquakes generate two seismic wave types: compressional (P) and shear (S) waves, whose coupling in anisotropic rocks leads to scattering, birefringence, and waves with hybrid polarizations. This varied behavior is helping geophysicists explore rock textures within Earth's mantle and crust, map present-day upper-mantle convection, and study the formation of lithospheric plates and the accretion of continents in Earth history.

The electric Moho

Since Mohorovicić discovered a dramatic increase in compressional seismic velocity at a depth of 54 km beneath the Kulpa Valley in Croatia, the 'Moho' has become arguably the most important seismological horizon in Earth owing to its role in defining the crust-mantle boundary. It is now known to be a ubiquitous feature of the Earth, being found beneath both the continents and the oceans, and is commonly assumed to separate lower-crustal mafic rocks from upper-mantle ultramafic rocks. Electromagnetic experiments conducted to date, however, have failed to detect a corresponding change in electrical conductivity at the base of the crust, although one might be expected on the basis of laboratory measurements. Here we report electromagnetic data from the Slave craton, northern Canada, which show a step-change in conductivity at Moho depths. Such resolution is possible because the Slave craton is highly anomalous, exhibiting a total crustal conductance of less than 1 Siemens--more than an order of magnitude smaller than other Archaean cratons. We also found that the conductivity of the uppermost continental mantle directly beneath the Moho is two orders of magnitude more conducting than laboratory studies on olivine would suggest, inferring that there must be a connected conducting phase.