Evidence of low velocity layers at the top and bottom of the Mantle Transition Zone (MTZ) below the Uttarakhand Himalaya, India

The Mantle Transition Zone (MTZ) beneath the Uttarakhand Himalaya has been modelled using Common Conversion Point (CCP) stacking and depth-migration of radial P-receiver functions. In the Uttarakhand Himalaya region, the depths of the 410-km discontinuity (d410) and the 660-km discontinuity (d660) are estimated to be approximately 406 ± 8 km and 659 ± 10 km, respectively. Additionally, the thickness of the mantle transition zone (MTZ) is modelled to be 255 ± 7 km. The average arrival times for d410 and d660 conversions are (44.47 ± 1.33) s and (71.08 ± 1.29) s, respectively, indicating an undisturbed slightly deeper d410 and a deformed noticeably deeper d660 in the area. The model identifies the characteristics of the d410 and d660 mantle discontinuities beneath the Lesser Himalayan region, revealing a thickening of the MTZ towards northeast, which could be due to gradual cooling or thickening of the Indian lithosphere towards its northern limit. We simulate a low-velocity layer (perhaps partially molten) above the d410 discontinuity at depths of 350 to 385 km, indicating the presence of a hydrated MTZ beneath the area. We also interpret a negative phase at d660 as a low-velocity layer between 590 and 640 km depths, which could be attributed to the accumulation of old subducted oceanic materials or increased water content at the bottom of the MTZ. Our results suggest the presence of residues from paleo-subducted lithospheric slabs in and below the mantle transition zone underlying the Uttarakhand Himalayas.

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.

Evidence of structural segmentation of the Uttarakhand Himalaya and its implications for earthquake hazard

The earthquake hazard associated with the Main Himalayan Thrust (MHT) is a critical issue for India and its neighbouring countries in the north. We used data from a dense seismic network in Uttarakhand, India, to model the lateral variations in the depths of MHT (2-6% drop in V_s at 12-21 km depths), Moho (a sharp increase in V_s (by ~ 0.5-0.7 km/s) at 39-50 km depths) and lithosphere (a marked decrease in V_s(~ 1-3%) at 136-178 km depths), across the Himalayan collisional front. Our joint inversion of radial PRFs and group velocity dispersion data of Rayleigh waves detects three NNE trending transverse lithospheric blocks segmenting the lithosphere in Uttarakhand Himalaya, which spatially correlate well with the northward extension of the Delhi -Haridwar Indian basement ridge, an inferred tectonic boundary and great boundary fault, respectively. Our radial receiver function imaging detects highly deformed and segmented crustal and lithospheric structures associated with three mapped transverse lithospheric blocks, suggesting a reduction in rupture lengths of future earthquakes, thereby, reducing earthquake hazards in Uttarakhand.

Evidence of secular variation in Archean crust formation in the Eastern Indian Shield

Understanding the dominant crustal accretion model in any Archean craton is the key to understanding the dominant geodynamic process responsible for early crust formation during the Hadean (> 4.0 Ga) and Archaean (4.0–2.5 Ga). The continental crust has been proposed to have formed through either horizontal/vertical accretion related to subduction or mantle plume tectonic processes. Here, the Moho depths and average crustal Vp/Vs ratios are modelled at 16 broadband stations in the Eastern Indian Shield (EIS) through HK stacking of radial P-receiver functions (PRFs). These modelled parameters are used to test both plume and subduction models, which might have played a key role in the crustal accretion of the EIS throughout the Archean. We observe a correlation between crustal age and composition within the ellipsoidal Paleoarchean cratonic domain in the Singhbhum-Odisha-Craton (SOC), which reveals an increase in age from the younger granitoid core of the SOC (with thinning of felsic crust) to the surrounding older greenstone belts (with thickening of felsic crust). A thinner mafic crust resulting from multiple magmatic events characterizes the neighbouring Meso-Proterozoic Chotanagpur Granitic Gneissic terrain (CGGT). The Common Conversion Point (CCP) image of radial PRFs reveals northward subduction of the Paleoarchean SOC below the Meso-Proterozoic CGGT.

Wedge tectonics in South China: constraints from new seismic data

Collisional orogens form when tectonic forces amalgamte fragments of Earth's continental lithosphere. The sutures between individual fragments, or terranes, are potential sites of weakness that facilitate subsequent continental breakup. Therefore, the lithospheric architecture of collisional orogens provides key information for evaluating the long-term evolution of the continental interior: for example, the South China Block (SCB), where the tectonic history is severely obscured by extensive surface deformation, magmatism, and metamorphism. Using new passive-source seismic models, we show a contrasting seismic architecture across the SCB, with three prominent crustal dipping structures across the Jiangnan Orogen. Combined with constraints from multi-disciplinary regional geophysical datasets, these pronounced dipping patterns are interpreted as relict wedge-like lithospheric deformation zones initiated in the fossil collisions that assembled the Yangtze Block and the SCB. The overall trend of these tectonic wedges implies successive crustal growth along paleo-continental margins and is indicative of northward subduction and docking of accretional terranes. In contrast, no such dipping structures are preserved in the Cathaysia Block, indicating a weak and reorganized lithosphere. The variations in the deformation responses across the SCB reflect the long-term modifications of the lithosphere caused by prolonged collision and extension events throughout the tectonic history of the SCB. Our results demonstrate the critical roles that suture zones played in the successive growth and evolution of the continental lithosphere.

Improving Constraints on Planetary Interiors With PPs Receiver Functions

Seismological constraints obtained from receiver function (RF) analysis provide important information about the crust and mantle structure. Here, we explore the utility of the free-surface multiple of the P-wave (PP) and the corresponding conversions in RF analysis. Using earthquake records, we demonstrate the efficacy of PPs-RFs before illustrating how they become especially useful when limited data is available in typical planetary missions. Using a transdimensional hierarchical Bayesian deconvolution approach, we compute robust P-to-S (Ps)- and PPs-RFs with InSight recordings of five marsquakes. Our Ps-RF results verify the direct Ps converted phases reported by previous RF analyses with increased coherence and reveal other phases including the primary multiple reverberating within the uppermost layer of the Martian crust. Unlike the Ps-RFs, our PPs-RFs lack an arrival at 7.2 s lag time. Whereas Ps-RFs on Mars could be equally well fit by a two- or three-layer crust, synthetic modeling shows that the disappearance of the 7.2 s phase requires a three-layer crust, and is highly sensitive to velocity and thickness of intra-crustal layers. We show that a three-layer crust is also preferred by S-to-P (Sp)-RFs. While the deepest interface of the three-layer crust represents the crust-mantle interface beneath the InSight landing site, the other two interfaces at shallower depths could represent a sharp transition between either fractured and unfractured materials or thick basaltic flows and pre-existing crustal materials. PPs-RFs can provide complementary constraints and maximize the extraction of information about crustal structure in data-constrained circumstances such as planetary missions.

Iterative passive-source location estimation and velocity inversion using geometric-mean reverse-time migration and full-waveform inversion

Passive-seismic provides useful information for reservoir monitoring and structural imaging; for example, the locations of microseismic events are helpful to understand the extension of the hydraulic fracturing. However, passive-seismic imaging still faces some challenges. First, it is not easy to know where the passive-seismic events happened, which is known as passive-source locating. Additionally, the accuracy of the subsurface velocity model will influence the accuracy of the estimated passive-source locations and the quality of the structural imaging obtained from the passive-seismic data. Therefore the velocity inversion using the passive-seismic data is required to provide the velocity with higher accuracy. Focusing on these challenges, we develop an iterative passive-source location estimation and velocity inversion method using geometric-mean reverse-time migration (GmRTM) and full-waveform inversion (FWI). In each iteration, the source location is estimated using a high-resolution GmRTM method, which provides a better focusing of passive-source imaging compared to conventional wavefield scanning method. The passive-source FWI is then followed to optimize the velocity model using the estimated source location provided by GmRTM. The source location estimation and velocity inversion are implemented sequentially. We evaluate this iterative method using the Marmousi model data set. The experiment result and sensitivity analysis indicate that the proposed method is effective to locate the sources and optimize velocity model in the areas with complicated subsurface structures and noisy recordings.

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.

The non-cylindrical crustal architecture of the Pyrenees

We exploit the data from five seismic transects deployed across the Pyrenees to characterize the deep architecture of this collisional orogen. We map the main seismic interfaces beneath each transect by depth migration of P-to-S converted phases. The migrated sections, combined with the results of recent tomographic studies and with maps of Bouguer and isostatic anomalies, provide a coherent crustal-scale picture of the belt. In the Western Pyrenees, beneath the North Pyrenean Zone, a continuous band of high density/velocity material is found at a very shallow level (~10 km) beneath the Mauleon basin and near Saint-Gaudens. In the Western Pyrenees, we also find evidence for northward continental subduction of Iberian crust, down to 50-70 km depth. In the Eastern Pyrenees, these main structural features are not observed. The boundary between these two domains is near longitude 1.3 °E, where geological field studies document a major change in the structure of the Cretaceous rift system, and possibly a shift of its polarity, suggesting that the deep orogenic architecture of the Pyrenees is largely controlled by structural inheritance.

Magmatic evolution of a Cordilleran flare-up and its role in the creation of silicic crust

The role of magmatic processes as a significant mechanism for the generation of voluminous silicic crust and the development of Cordilleran plateaus remains a lingering question in part because of the inherent difficulty in quantifying plutonic volumes. Despite this difficulty, a growing body of independently measured plutonic-to-volcanic ratios suggests the volume of plutonic material in the crust related to Cordilleran magmatic systems is much larger than is previously expected. To better examine the role of crustal magmatic processes and its relationship to erupted material in Cordilleran systems, we present a continuous high-resolution crustal seismic velocity model for an ~800 km section of the active South American Cordillera (Puna Plateau). Although the plutonic-to-volcanic ratios we estimate vary along the length of the Puna Plateau, all ratios are larger than those previously reported (~30:1 compared to 5:1) implying that a significant volume of intermediate to silicic plutonic material is generated in the crust of the central South American Cordillera. Furthermore, as Cordilleran-type margins have been common since the onset of modern plate tectonics, our findings suggest that similar processes may have played a significant role in generating and/or modifying large volumes of continental crust, as observed in the continents today.

Seismological constraints on the crustal structures generated by continental rejuvenation in northeastern China

Crustal rejuvenation is a key process that has shaped the characteristics of current continental structures and components in tectonic active continental regions. Geological and geochemical observations have provided insights into crustal rejuvenation, although the crustal structural fabrics have not been well constrained. Here, we present a seismic image across the North China Craton (NCC) and Central Asian Orogenic Belt (CAOB) using a velocity structure imaging technique for receiver functions from a dense array. The crustal evolution of the eastern NCC was delineated during the Mesozoic by a dominant low seismic wave velocity with velocity inversion, a relatively shallow Moho discontinuity, and a Moho offset beneath the Tanlu Fault Zone. The imaged structures and geochemical evidence, including changes in the components and ages of continental crusts and significant continental crustal growth during the Mesozoic, provide insight into the rejuvenation processes of the evolving crust in the eastern NCC caused by structural, magmatic and metamorphic processes in an extensional setting. The fossil structural fabric of the convergent boundary in the eastern CAOB indicates that the back-arc action of the Paleo-Pacific Plate subduction did not reach the hinterland of Asia.

Subduction-driven recycling of continental margin lithosphere

Whereas subduction recycling of oceanic lithosphere is one of the central themes of plate tectonics, the recycling of continental lithosphere appears to be far more complicated and less well understood. Delamination and convective downwelling are two widely recognized processes invoked to explain the removal of lithospheric mantle under or adjacent to orogenic belts. Here we relate oceanic plate subduction to removal of adjacent continental lithosphere in certain plate tectonic settings. We have developed teleseismic body wave images from dense broadband seismic experiments that show higher than expected volumes of anomalously fast mantle associated with the subducted Atlantic slab under northeastern South America and the Alboran slab beneath the Gibraltar arc region; the anomalies are under, and are aligned with, the continental margins at depths greater than 200 kilometres. Rayleigh wave analysis finds that the lithospheric mantle under the continental margins is significantly thinner than expected, and that thin lithosphere extends from the orogens adjacent to the subduction zones inland to the edges of nearby cratonic cores. Taking these data together, here we describe a process that can lead to the loss of continental lithosphere adjacent to a subduction zone. Subducting oceanic plates can viscously entrain and remove the bottom of the continental thermal boundary layer lithosphere from adjacent continental margins. This drives surface tectonics and pre-conditions the margins for further deformation by creating topography along the lithosphere-asthenosphere boundary. This can lead to development of secondary downwellings under the continental interior, probably under both South America and the Gibraltar arc, and to delamination of the entire lithospheric mantle, as around the Gibraltar arc. This process reconciles numerous, sometimes mutually exclusive, geodynamic models proposed to explain the complex oceanic-continental tectonics of these subduction zones.

Active foundering of a continental arc root beneath the southern Sierra Nevada in California

Seismic data provide images of crust-mantle interactions during ongoing removal of the dense batholithic root beneath the southern Sierra Nevada mountains in California. The removal appears to have initiated between 10 and 3 Myr ago with a Rayleigh-Taylor-type instability, but with a pronounced asymmetric flow into a mantle downwelling (drip) beneath the adjacent Great Valley. A nearly horizontal shear zone accommodated the detachment of the ultramafic root from its granitoid batholith. With continuing flow into the mantle drip, viscous drag at the base of the remaining approximately 35-km-thick crust has thickened the crust by approximately 7 km in a narrow welt beneath the western flank of the range. Adjacent to the welt and at the top of the drip, a V-shaped cone of crust is being dragged down tens of kilometres into the core of the mantle drip, causing the disappearance of the Moho in the seismic images. Viscous coupling between the crust and mantle is therefore apparently driving present-day surface subsidence.

Low-velocity zone atop the 410-km seismic discontinuity in the northwestern United States

The seismic discontinuity at 410 km depth in the Earth's mantle is generally attributed to the phase transition of (Mg,Fe)2SiO4 (refs 1, 2) from the olivine to wadsleyite structure. Variation in the depth of this discontinuity is often taken as a proxy for mantle temperature owing to its response to thermal perturbations. For example, a cold anomaly would elevate the 410-km discontinuity, because of its positive Clapeyron slope, whereas a warm anomaly would depress the discontinuity. But trade-offs between seismic wave-speed heterogeneity and discontinuity topography often inhibit detailed analysis of these discontinuities, and structure often appears very complicated. Here we simultaneously model seismic refracted waves and scattered waves from the 410-km discontinuity in the western United States to constrain structure in the region. We find a low-velocity zone, with a shear-wave velocity drop of 5%, on top of the 410-km discontinuity beneath the northwestern United States, extending from southwestern Oregon to the northern Basin and Range province. This low-velocity zone has a thickness that varies from 20 to 90 km with rapid lateral variations. Its spatial extent coincides with both an anomalous composition of overlying volcanism and seismic 'receiver-function' observations observed above the region. We interpret the low-velocity zone as a compositional anomaly, possibly due to a dense partial-melt layer, which may be linked to prior subduction of the Farallon plate and back-arc extension. The existence of such a layer could be indicative of high water content in the Earth's transition zone.

Formation and evolution of Archaean cratons: insights from southern Africa

Archaean cratons are the stable remnants of Earth’s early continental lithosphere, and their structure, composition and survival over geological time make them unique features of the Earth’s surface. The Kaapvaal Project of southern Africa was organized around a broadly diverse scientific collaboration to investigate fundamental questions of craton formation and mantle differentiation in the early Earth. The principal aim of the project was to characterize the physical and chemical nature of the crust and mantle of the cratons of southern Africa in geological detail, and to use the 3D seismic and geochemical images of crustal and mantle heterogeneity to reconstruct the assembly history of the cratons. Seismic results confirm that the structure of crust and tectospheric mantle of the cratons differs significantly from that of post-Archaean terranes. Three-dimensional body-wave tomographic images reveal that high-velocity mantle roots extend to depths of at least 200 km, and locally to depths of 250–300 km beneath cratonic terranes. No low-velocity channel has been identified beneath the cratonic root. The Kaapvaal Craton was modified approximately 2.05 Ga by the Bushveld magmatic event, and the mantle beneath the Bushveld Province is characterized by relatively low seismic velocities. The crust beneath undisturbed Archaean craton is relatively thin (c. 35–40 km), unlayered and characterized by a strong velocity contrast across a sharp Moho, whereas post-Archaean terranes and Archaean regions disrupted by large-scale Proterozoic magmatic or tectonic events are characterized by thicker crust, complex Moho structure and higher seismic velocities in the lower crust. A review of Re-Os depletion model age determinations confirms that the mantle root beneath the cratons is Archaean in age. The data show also that there is no apparent age progression with depth in the mantle keel, indicating that its thickness has not increased over geological time. Both laboratory experiments and geochemical results from eclogite xenoliths suggest that subduction processes played a central role in the formation of Archaean crust, the melt depletion of Archaean mantle and the assembly of early continental lithosphere. Co-ordinated geochronological studies of crustal and mantle xenoliths have revealed that both crust and mantle have experienced a multi-stage history. The lower crust in particular retains a comprehensive record of the tectonothermal evolution of the lithosphere. Analysis of lower-crustal xenoliths has shown that much of the deep craton experienced a dynamic and proteracted history of tectonothermal activity that is temporally associated with events seen in the surface record. Cratonization thus occurred not as a discrete event, but in stages, with final stabilization postdating crustal formation.

Mantle discontinuity structure beneath the southern east pacific rise from P-to-S converted phases

Receiver functions derived from teleseismic body waves recorded by ocean-bottom seismometers on the southern East Pacific Rise reveal shear waves converted from compressional waves at the mantle discontinuities near 410- and 660-kilometer depth. The thickness of the mantle transition zone between the two discontinuities is normal relative to the global average and indicates that upwelling beneath the southern East Pacific Rise is not associated with an excess temperature in the mantle transition zone.