Search

Earthquakes are frequently accompanied by electromagnetic (EM) anomalies. These anomalies are thought to be caused by earthquakes but the generation mechanism is still unclear. The piezoelectric effect has been proposed as a possible mechanism but the EM responses to earthquakes due to such an effect has not been well understood. In this article, we study the EM signals generated by an earthquake source due to the piezoelectric effect. We develop a semi-analytical method to solve the seismic and EM fields in a 3D layered model and conduct numerical simulations to investigate the characteristics of the EM fields. The results show that the earthquake can generate two kinds of EM signals. One is the early-EM signal which arrives earlier than the seismic wave. The other is the co-seismic EM signal accompanying the seismic wave. For an earthquake the co-seismic electric field can reach ~10 μV/m and the magnetic field can reach ~10-4 nT. We also study the sensitivity of the co-seismic EM fields to the rock conductivity. The results show that the co-seismic EM fields are mainly affected by the conductivity of the shallow layer, and they are also affected by the conductivity of the deep layer when the top layer is thin.

We study the Moho, the mid-lithospheric discontinuity (MLD), and the lithosphere-asthenosphere boundary (LAB) from southern Africa to northern Arabia, from Archean cratons to active rifts, at 1° resolution using our comprehensive new database of shear-wave receiver functions (SRFs). The good agreement between the Moho depth obtained from our SRFs and published P-wave receiver function (PRF) results provides confidence that our images of deeper lithospheric discontinuities are robust, including boundaries not normally visible on PRFs. We map the Moho and a deeper negative velocity gradient (NVG) almost everywhere we have data coverage. Our synthetic tests and comparisons of SRFs processed with and without deconvolution, and with varying filter parameters, indicate the observed NVG represents earth structure, not a processing artifact. Depth comparisons with seismic tomography and tectonothermal age studies suggest the NVG represents the MLD beneath Archean cratons but represents the LAB beneath non-cratonic regions. Both preserved crustal thickness and lithospheric thickness in the Nubia-Somalia-Arabia plates are statistically thinner for Phanerozoic and late Proterozoic terranes and older regions reactivated during these eras, than for cratons not reworked since the early Proterozoic or Archean. In contrast, NVG depth is uniform for all tectonothermal ages, though with a possible increase in amplitude with age. The equivalence of NVG depth and LAB depth in Phanerozoic lithosphere suggests that low-wavespeed compositions are frozen into the lithosphere as it thickens by cooling, forming our observed MLD at the present day.

SUMMARY The rate of occurrence of anomalous ultra-low frequency electromagnetic (ULFEM) pulses has been claimed to have increased days to weeks prior to the M5.4 2007 and M4.0 2010 Alum Rock earthquakes. We re-examine the previously reported ultra-low frequency (ULF: 0.01–10 Hz) magnetic data recorded at a QuakeFinder site located 9 km from the earthquake hypocentre, and compare to data from a nearby Stanford-USGS site located 42 km from the hypocentre, to analyse the characteristics of the pulses and assess their origin. Using pulse definitions and pulse-counting algorithms analogous to those previously reported, we corroborate the increase in pulse counts before the 2007 Alum Rock earthquake at the QuakeFinder station, but we note that the number of pulses depends on chosen temporal and amplitude detection thresholds. These thresholds are arbitrary because we lack a clear physical model or basis for their selection. We do not see the same increase in pulse counts before the 2010 Alum Rock earthquake at the QuakeFinder or Stanford-USGS stations. In addition, the majority of pulses in the QuakeFinder data and Stanford-USGS data do not match temporally, indicating the pulses lack a common origin and are not from lightning or solar-driven ionospheric/magnetospheric disturbances. Our assessment of the temporal distribution of pulse counts shows pulse counts increase during peak human activity hours, suggesting these pulses result from local cultural noise and are not tectonic in origin. The many unknowns about the character and even existence of precursory earthquake pulses means that standard numerical and statistical tests cannot easily be applied. Yet here we show that exhaustive investigation of many different aspects of ULFEM signals can be used to properly characterize their origin.

SUMMARY Very clear coseismic electric and magnetic signals accompanying seismic waves were observed during the 2017 Mw 6.5 Jiuzhaigou earthquake, which took place in western China. In order to understand the generation mechanism of these observed signals, we simulate electric and magnetic responses to this specific earthquake based on three mechanisms, namely, the electrokinetic effect, the motional induction effect and the rotation effect of the coil-type magnetometer. We conduct the simulations using a point source model and a realistic layered earth model and compare to the observed data in the frequency band 0.05–0.3 Hz. Our results show that the electrokinetic effect can explain the observed electric fields in both waveform and amplitude, but it cannot explain the magnetic signals accompanying the Rayleigh wave. The motional induction effect cannot explain either the coseismic electric or magnetic data because it predicts much weaker coseismic electric and magnetic fields than the observed data. The magnetic fields resulting from the rotation of the magnetometer agree with the observed data in the waveforms though their amplitudes are two to four times smaller than the observed data. Our simulations suggest that the electrokinetic effect is responsible for the generation of coseismic electric fields and that rotation of the coil magnetometer is likely the main cause of coseismic magnetic fields. The results improve our interpretation of the coseismic electromagnetic (EM) phenomenon and are useful for understanding other kinds of earthquake-associated EM phenomena.

SUMMARY Virtual deep seismic sounding (VDSS) uses the arrival time of post-critical SsPmp relative to the direct S wave to infer Moho depth at the Pmp reflection point. Due to the large offset between the virtual source and the receiver, SsPmp is more sensitive to lateral variations of structures than near-vertical phases such as Ps, which is used to construct conventional P receiver functions. However, the way post-critical SsPmp is affected by lateral variations in lithospheric structure is not well understood, and previous studies largely assumed a 1-D structure when analysing SsPmp waveforms. Here we present synthetic tests with various 2-D models to show that lateral variations in lithospheric structures, from the lithosphere–asthenosphere boundary (LAB) to sedimentary basins, profoundly affect traveltime, phase and amplitude of post-critical SsPmp, and that a 1-D approximation is usually inappropriate when analysing 2-D data. Despite these strong effects we show, with synthetic examples and the ChinArray data from the Ordos Block in northern China, that a simple ray-theory-based back-projection method can retrieve the geometry of the crust–mantle boundary (CMB) given array observations in cases with moderate lateral variations in the CMB and/or the LAB. The success of our back-projection method indicates that ray-theory approximations are sufficient in modelling SsPmp traveltimes in the presence of moderate lateral heterogeneity. In contrast, we show that the ray theory is generally insufficient in modelling SsPmp phase shifts in a strongly heterogeneous lithosphere due to non-planar downgoing P waves incident at the CMB. Nonetheless, our results demonstrate the feasibility of direct imaging of the CMB with post-critical SsPmp even in the presence of 2-D variations of lithospheric structure.

The Himalaya has experienced large damaging earthquakes over the past few centuries, most recently the damaging 25 April 2015 M 7.8 Gorkha earthquake in Nepal. Because of the continued earthquake risk presented by the continental collisional plate boundary at the Main Himalayan thrust and the high population densities in the region, collecting and processing data related to recent large earthquakes in this region is critically important for improving our understanding of the regional tectonics and earthquake hazard. Following the 2015 Gorkha earthquake, we deployed a National Science Foundation-funded rapid-response aftershock network known as the Nepal Array Measuring Aftershock Seismicity Trailing Earthquake network across the rupture area for 11 months beginning 7 weeks after the mainshock. The network consisted of 41 broadband and short-period seismometers, and 14 strong-motion sensors at 46 sites across eastern and central Nepal. The network spanned a region approximately 210 km along strike by 110 km across strike with a station spacing of 20–25 km. In this article, we report lessons learned from this deployment as well as details of the publicly accessible dataset including data recovery, data quality, and potential for future research.

During continent–continent collision, does the downgoing continental plate underplate far inboard of the collisional boundary or does it subduct steeply into the mantle, and how is this geometry manifested in the mantle flow field? We test conflicting models for these questions for Earth’s archetypal continental collision forming the Himalaya and Tibetan Plateau. Air-corrected helium isotope data (3He/4He) from 225 geothermal springs (196 from our group, 29 from the literature) delineate a boundary separating a Himalayan domain of only crustal helium from a Tibetan domain with significant mantle helium. This 1,000-km-long boundary is located close to the Yarlung-Zangbo Suture (YZS) in southern Tibet from 80 to 92°E and is interpreted to overlie the “mantle suture” where cold underplated Indian lithosphere is juxtaposed at80 km depth against a sub-Tibetan incipiently molten asthenospheric mantle wedge. In southeastern Tibet, the mantle suture lies 100 km south of the YZS, implying delamination of the mantle lithosphere from the Indian crust. This helium-isotopic boundary helps resolve multiple, mutually conflicting seismological interpretations. Our synthesis of the combined data locates the northern limit of Indian underplating beneath Tibet, where the Indian plate bends to steeper dips or breaks off beneath a (likely thin) asthenospheric wedge below Tibetan crust, thereby defining limited underthrusting for the Tibetan continental collision.

In April 2015, the lower locked portion of the Main Himalayan Thrust ruptured beneath Nepal, causing the disastrous Mw 7.8 Gorkha earthquake. Elucidating the enigmatic geometry of this plate boundary fault is important for understanding the nucleation and arrest of large earthquake ruptures as well as the seismic hazard, topography and tectonics of the Himalaya. Here we interpret the geometry of the Main Himalayan Thrust from the spatial distribution and rupture patterns of a dynamic sequence of aftershocks following the Gorkha earthquake, which were recorded by a rapidly deployed dense seismic network. We find that the thrust comprises two north-dipping subhorizontal planes that are connected by a system of bounded imbricate thrust faults; this structure is known as a duplex. We propose that this duplex acts as an impediment to plate convergence and accommodates tectonic stress along its complex system of faults. Such a prominent structure to the Main Himalayan Thrust is consistent with surface geological studies but challenges geophysically derived conventional models with simpler geometries. The Main Himalayan Thrust comprises two fault planes connected by imbricated faults, a structure that impedes convergence, according to an analysis of the distribution and orientation of aftershocks of the 2015 Gorkha earthquake in Nepal.

The occurrence of earthquakes in the continental upper mantle is highly debated, and bears directly on lithospheric rheology (e.g. “jelly-sandwich” vs. “creme-brulee” models). Because...

The deep structure of the continental collision between India and Asia and whether India's lower crust is underplated beneath Tibet or subducted into the mantle remain controversial. It is also unknown whether the active normal faults that facilitate orogen-parallel extension of Tibetan upper crust continue into the lower crust and upper mantle. Our receiver-function images collected parallel to the India-Tibet collision zone show the 20-km-thick Indian lower crust that underplates Tibet at 88.5-92°E beneath the Yarlung-Zangbo suture is essentially absent in the vicinity of the Cona-Sangri and Pumqu-Xainza grabens, demonstrating a clear link between upper-crustal and lower-crustal thinning. Satellite gravity data that covary with the thickness of Indian lower crust are consistent with the lower crust being only ∼30% eclogitized so gravitationally stable. Deep earthquakes coincide with Moho offsets and with lateral thinning of the Indian lower crust near the bottom of the partially eclogitized Indian lower crust, suggesting the Indian lower crust is locally foundering or stoping into the mantle. Loss of Indian lower crust by these means implies gravitational instability that can result from localized rapid eclogitization enabled by dehydration reactions in weakly hydrous mafic granulites or by volatile-rich asthenospheric upwelling directly beneath the two grabens. We propose that two competing processes, plateau formation by underplating and continental loss by foundering or stoping, are simultaneously operating beneath the collision zone.

The Tibetan Plateau's enigmatic collapse, which was indicated based on extensional faulting in the interior of the plateau and large eastward-directed strike-slip faults, has inspired intensive geologic inquiry and debate. Interaction of the subducting Indian plate with the overlying Asian lithosphere is one factor that may be affecting the plateau. A more detailed image of the subducting Indian plate can help constrain its role, but regional high-resolution seismic data have not been widely available at a relevant scale. Here we present an integrated interpretation based on two types of seismic datasets: (1) two N-S oriented deep seismic-reflection profiles that cross the Yarlung-Zangbo suture, and (2) two E-W receiver function profiles that cross the Xainza-Dingjye and Nyima-Tingri rifts, respectively. These images reveal different crustal-scale structures in the dominant collision zone between the western and central regions, as well as distinctly offset Moho discontinuities across the N-S trending rift grabens in southern Tibet. These crustal variations, combined with previous tomographic studies in Tibet, outline an easterly tilt of the subducting Indian slab, along which crust-mantle decoupling occurred. Together with the spatio-temporal distribution of synchronous potassium-rich volcanics exposed within the grabens in southern Tibet, this offset Moho structure beneath the grabens lends support to the hypothesis that eastward steepening of the subducting Indian slab was accompanied by slab tearing beneath each north-south striking rift. This overall crustal geometric variation suggests a stepwise flattening of the torn Indian slab from west to the east, a process that brought easterly migration of gravitational instability to the overriding Tibetan crust that drove collapse of the Tibetan Plateau once gravitational instability reached a critical level to the east.

We use earthquake geothermometry, measured heat flow, and structural constraints from P-wave receiver functions to model the thermal evolution of the lithosphere beneath Harrat Lunayyir. We suggest that the lithosphere thinned to its present 60-km thickness in a second stage of lithospheric thinning at 15–12 Ma following initial Red Sea extension at ∼27 Ma. Harrat Lunayyir is an active volcanic field located in the Arabian Shield >150 km east of the Red Sea rift axis. In the lithospheric mantle beneath Harrat Lunayyir we locate 64 high-frequency earthquakes at depths of 42–48 km, all with m L ≤ 2.5 . These brittle-failure earthquakes must have nucleated at relatively low temperatures, based upon global maximum nucleation depths and temperature-dependent-deformation experimental results. Therefore, the mantle earthquakes show that the upper mantle lithosphere is not in thermal equilibrium with the shallow (60 km) underlying asthenosphere. Our thermal modeling indicates that the lithosphere beneath Harrat Lunayyir thinned to its current 60-km thickness at 12 ± 2 Ma, as constrained by thermal modeling of: (1) surface heat-flow; (2) the depth to the mid-crustal brittle–ductile transition, and (3) the depth to the upper-mantle brittle–ductile transition.

At the boundary between the Western Transverse Ranges province and Inner Continental Borderland of Southern California, strain is partitioned across the sinistral-oblique Anacapa-Dume Fault system. Sycamore Knoll, a pop-up structure 20 km west of Point Dume in the hanging wall of the north-dipping Anacapa-Dume Fault, stands out as an anomalous bathymetric high along the Southern California coast. We synthesize new and existing geophysical datasets, including multibeam bathymetry/backscatter and single/multichannel seismic, to understand the morphology and tectonic implications of Sycamore Knoll. We identify a complete eustatic low-stand submerged marine-terrace package encircling Sycamore Knoll platform that corresponds to the Last Glacial Maximum (LGM). The marine-terrace deposits and wave-planed surface of Sycamore Knoll require a late Quaternary uplift rate of 0.55 (+0.5/−0.3) mm/yr relative to sea-level. By combining our uplift rate with previous estimates of basin subsidence and a revised structural model of the Anacapa-Dume Fault, we determine a maximum slip rate in latest Quaternary time of 1.50 (+0.5/−0.3) mm/yr and hence a post-LGM rate of convergence across the Anacapa-Dume Fault of 0.60 (+0.2/−0.1) mm/yr.

We discriminate upper-mantle earthquakes from crustal earthquakes based on the amplitude ratio of seismic waves Sn and Lg (‘ S n / L g ’), a prominent feature of regional seismograms that is visible to the naked eye. Crucially, our new method uses only the waveforms of the candidate earthquake, unlike previous methods to identify upper-mantle earthquakes that introduce potentially large errors by comparing hypocentral depths with independent measures of crustal thickness. Our synthetics show that the Love-wave higher modes that form individual Sn and Lg Airy phases on the transverse component are preferentially excited when the source is respectively below or above the Moho. We use three previously recognized mantle events from southern Tibet to validate our new approach and show that focal mechanism, intrinsic attenuation, geometrical spreading etc., can be ignored to first order. We then identify two new upper-mantle earthquakes in NW Tibet where, previously, only lower-crustal events had been reliably demonstrated, thereby showing this NW-Tibet lithosphere is seismogenic at all depths i.e. that upper-mantle and lower-crustal earthquakes co-exist. Our method has potential for expanding the global catalog of continental earthquakes reliably determined to be close above or close below the Moho.

Virtual Deep Seismic Sounding (VDSS) has recently emerged as a novel method to image the crust-mantle-boundary (CMB) and potentially other lithospheric boundaries. In Liu et al., 2018 (“Part 1”), we showed that the arrival time and waveform of post-critical SsPmp, the post-critical reflection phase at the CMB used in VDSS, is sensitive to multiple attributes of the crust and upper mantle. Here, we present a synthesis on the methodology of deriving Moho depth, crustal average Vp and uppermost mantle Vp from single-station observations of postcritical SsPmp under an 1D assumption. We first verify our methods on synthetic examples and then substantiate it with a case study using data collected by the Yellowknife and POLARIS array in SW Slave Province, which shows good agreement between crustal and upper-mantle properties derived with VDSS and the ones given by previous studies using active source experiments and conventional P receiver functions (PRF). Moreover, we show that PRF-VDSS joint analysis is capable of placing tighter constraint on average crustal composition, which cannot be achieved with either method alone. We demonstrate with our PRF-VDSS joint analysis that the SW Slave Province, Canada has an intermediate crust composition, which likely indicates a Neoarchean age

Gross morphological structures of the Main Himalayan Thrust control the rupture length of the 2015 Mw 7.8 Nepal earthquake., The Himalaya orogenic belt produces frequent large earthquakes that affect population centers along a length of over 2500 km. The 2015 Gorkha, Nepal earthquake (Mw 7.8) ruptured the Main Himalayan Thrust (MHT) and allows direct measurements of the behavior of the continental collision zone. We study the MHT using seismic waveforms recorded by local stations that completely cover the aftershock zone. The MHT exhibits clear lateral variation along geologic strike, with the Lesser Himalayan ramp having moderate dip on the MHT beneath the mainshock area and a flatter and deeper MHT beneath the eastern end of the aftershock zone. East of the aftershock zone, seismic wave speed increases at MHT depths, perhaps due to subduction of an Indian basement ridge. A similar magnitude wave speed change occurs at the western end of the aftershock zone. These gross morphological structures of the MHT controlled the rupture length of the Gorkha earthquake.

We present new shear-wave splitting data showing systematic lateral variations in upper-mantle anisotropy across the plate boundary in southernmost California (USA). Beneath the Peninsular Ranges batholith, fast polarization directions parallel the direction of former Farallon subduction, suggestive of a slab remnant. Near the eastern edge of the batholith, across the Elsinore fault, fast polarization directions change rapidly to align with the direction of San Andreas fault shear. We infer that the Elsinore fault penetrates the entire lithosphere and may represent a future localization of the plate boundary that is migrating west from the San Andreas fault. Beneath the Salton Trough and the Chocolate Mountains region, large splitting times, despite a very thin lithosphere, imply vertical melt pockets in the uppermost mantle aligned in the shear direction. Largest splitting times, ∼1.2 s, are seen closest to the Sand Hills fault that projects southeast from the San Andreas fault. Further east, in the southern Basin and Range province, fast directions align with North America absolute plate motion.

The fate of the Indian plate during continental collision with Asian terranes, and the proportion of the Indian crust that is underthrust or subducted beneath Tibet as opposed to transferred to the upper (Himalayan) plate, are much debated. The active geometry of low-angle underthrusting or subduction of the Indian plate beneath the Lesser and Greater Himalayan thrust sheets is well known from seismic imaging. Previously, only lower-resolution images of the Main Himalayan Thrust have been obtained beneath the Yarlung Zangbo suture that separates Indian and Asian rocks at the surface. It remains controversial whether the orogenic wedge between the Main Himalayan Thrust and the Yarlung Zangbo suture, formed of Indian crust transferred to the upper plate, is evolving by thrust-faulting in a critical-taper wedge or by southward extrusion of a ductile channel flow. Here we present a seismic reflection profile across the western Himalaya at 81.5° E, and show that the Main Himalayan Thrust dips ∼20° to ∼60 km depth beneath the Yarlung Zangbo suture, approaching a continuous Moho reflection at ∼70–75 km depth. The Indian crust being transported northwards beyond the Yarlung Zangbo suture is no more than ∼15 km thick, reduced from its original ∼40 km thickness by transfer of material from the lower plate to the upper plate through crustal-scale duplexing. The fate of the Indian plate during collision with Asia is debated. Seismic images of the western Himalaya reveal large-scale thrust faults that transfer Indian crust upwards, into the overriding Asian plate.

We perform a finite-frequency tomographic inversion to image 3D velocity structures beneath southern and central Tibet using teleseismic body-wave data recorded by the TIBET-31N passive seismic array as well as waveforms from previous temporary seismic arrays. High-velocity bodies dip ∼40° northward beneath the Himalaya and the Lhasa Terrane. We interpret these high-velocity anomalies as subducting Indian Continental Lithosphere (ICL). The ICL appears to extend further north in central Tibet than in eastern Tibet, reaching 350 km depth at ∼31°N along 85°E but at ∼30°N along 91°E. Low P- and S-wave velocity anomalies extend from the lower crust to ≥180 km depth beneath the Tangra Yum Co Rift, Yadong-Gulu Rift, and the Cona Rift, suggesting that rifting in southern Tibet may involve the entire lithosphere. The anomaly beneath Tangra Yum Co Rift extends down to about 180 km, whereas the anomalies west of the Yadong-Gulu Rift and east of the Cona Rift extend to more than 300 km depth. The low-velocity upper mantle west of the Yadong-Gulu Rift extends furthest north and appears to connect with the extensive upper-mantle low-velocity region beneath central Tibet. Thus the northward-subducting Indian Plate is fragmented along north–south breaks that permit or induce asthenospheric upwellings indistinguishable from the upper mantle of northern Tibet.

The Tibetan plateau is bounded in part by cratonic blocks, and in part by tectonothermally younger and weaker crust. How the Tibetan plateau interacts with bordering non-cratonic crust is a continuing debate, whether the plateau is currently enlarging its area, and if so, whether by thrusting or by lower-crustal flow, or both. The gently sloping northeastern margin of Tibet has been inferred to incorporate a lower-crustal channel flow outward towards the Ordos block that is a part of the North China Craton. We calculated receiver functions using teleseismic waveforms recorded by a dense nodal seismic array crossing from the Longxi basin of northeastern Tibet across a transition zone marked by the Liupan Shan (mountains) into the Ordos block. Our image shows the Longxi upper crust overthrusts onto the Ordos block building the Liupan Shan, and the Longxi lower crust underthrusts beneath the Ordos forming a 10-km thick crustal root. The Ordos crust forms a crustal-scale wedge inserted into the Longxi crust, and deformation is limited to a narrow zone of

We studied the crustal structure and tectonics in the north Tibetan Plateau from the Songpan-Ganzi terrane to the Qaidam Basin using teleseismic receiver-function imaging, across a major lithospheric boundary, the Kunlun-Qaidam boundary, where previous studies suggest a ~15–20-km change in crustal thickness from thicker crust in the Kunlun Mountains to thinner crust in the Qaidam Basin. We report P receiver functions for 70 stations, largely the International Deep Profiling of Tibet and the Himalaya (INDEPTH), phase IV, experiment. Our most dense station coverage is located along the roughly north-south INDEPTH-IV active-source seismic profile at approximately 95° E longitude. Azimuthal and geographical changes in the receiver functions reveal significant changes in crustal structure and Vp/Vs from across the study area. Receiver functions show strong converters that we interpret as the Moho at ~70 km depth beneath the Qiangtang, Songpan-Ganzi terranes and Kunlun Mountains and at ~50 km depth beneath the central Qaidam Basin. This large change in crustal thickness occurs >50 km north of the North Kunlun strike-slip fault, on which the 2001 M8.1 Kunlun earthquake occurred. Receiver functions for some of the stations north of the thickness change at the Kunlun-Qaidam boundary also show a deeper ~70-km bright converter in addition to the 50-km converter. The two converters appear to overlap by up to ~30 km in some locations along the south Qaidam Basin. We combine previous results with these new results to discuss implications for mechanisms for crustal thickening in the north Tibetan Plateau including crustal flow and crustal injection. At depths imaged here, shallower than ~100 km, we see no evidence of southward subduction of Eurasian lithosphere.

The relative contributions of tectonic and magmatic processes to continental rifting are highly variable. Magnetic, gravity, and seismic reflection data from Surprise Valley, California, in the northwest Basin and Range, reveal an intrabasin, fault-controlled, ∼10-m-thick dike at a depth of ∼150 m, providing an excellent example of the interplay between faulting and dike intrusion. The dike, likely a composite structure representing multiple successive intrusions, is inferred from modeling a positive magnetic anomaly that extends ∼35 km and parallels the basin-bounding Surprise Valley normal fault on the west side of the valley. A two-dimensional high-resolution seismic reflection profile acquired across the magnetic high images a normal fault dipping 56°E with ∼275 m of throw buried ∼60 m below the surface. Densely spaced gravity measurements reveal a

We use ambient-noise tomography to improve CVM-H11.9, a community velocity model of southern California. Our new 3D shear-velocity model with 0.05°x0.05° lateral and 1-km vertical blocks reveals new structure beneath the San Andreas Fault (SAF), Peninsular Ranges batholith (PRB), southern Sierra Nevada batholith (SNB) and the Salton Trough (ST). We use four years of data recorded on 849 broadband stations, vastly more than previous studies and including our own broadband Salton Seismic Imaging Project, a 40-station transect across the ST, as well as other campaign stations in both Mexico and the U.S.A. Mean lower-crust and upper-mantle wavespeeds (3.6 km/s @ 20km, 4.2 km/s @ 40km) are low by global standards. Across the SAF, southeast of San Gorgonio Pass, we observe vertical to steeply-dipping lateral velocity contrasts that extend beneath the Moho. Beneath the western PRB and westernmost southern SNB, we observe relatively high shear-velocities (≥ 3.8 km/s) in the lower crust that we interpret as the mafic roots of the overlying arc. Relatively high-velocity upper-mantle (up to ∼4.5 km/s) may be part of the intact arc, or possibly a remnant of the Farallon plate. Beneath the ST, we observe zones of low shear-velocity in the lower crust and upper mantle which permit up to ∼4.5% melt in the lower crust and up to ∼6% melt in the upper mantle, depending on the assumed composition and pore geometry. Our results preclude the existence of older continental crust beneath the ST and support the creation of new crust beneath the ST. This article is protected by copyright. All rights reserved.

a b s t r a c t West-east extension is a prominent tectonic feature of southern and central Tibet despite ongoing north- south (N-S) convergence between India and Eurasia. Knowledge of deep structure beneath the N-S trending rifts is key to evaluating models proposed for their origin, including gravitational collapse, oblique convergence along the arcuate plate boundary, and mantle upwelling. We model direct S and Moho-reflected SsPmp phases at teleseismic distances to constrain variations in crustal thickness across the major rifts crossed by a ∼900-km long, W-E broadband array in the Lhasa Terrane. Crustal thicknesses are ∼70-80 km. However, Moho depth decreases by ∼10 km within a horizontal distance of 100 km west of the Yadong-Gulu rift (YGR) and Nyainquentanghla mountains (NQTL). This Moho uplift, taken with deep, extensional focal mechanisms and reduced seismic velocity in the upper mantle, suggests that asthenospheric upwelling has significantly contributed to the pattern of extension across the YGR and NQTL. The ∼100-km separation between surface rift and Moho uplift is likely enabled by partial decoupling across a ductile middle crust. © 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

To understand the along-strike variation of crustal deformation and tectonic processes in the India–Tibet continental collision zone, we deployed a linear array of broadband seismic stations along 92° E to image lithospheric structure. Our receiver-function cross-section reveals a prominent negative converter dipping ∼20° north from ∼10–55 km depth below sea-level, almost through the whole crust, beneath the southern Lhasa terrane. We interpret it to be a manifestation of the Yarlung–Zangbo suture zone (YZS) separating the continental crust of the Indian and Eurasian plates. This implies the hypothesized channel-flows of Indian middle crust extruding southwards from Tibet are limited at this longitude to the southernmost portion of the Lhasa terrane. A positive converter, consistent with previous suggestions of eclogite formation, is seen about 10–15 km above the Moho and continuing 50 km north of the 20°-dipping YZS converter. We image this positive converter continuously from ∼60 km south of the surface trace of the YZS to the vicinity of the Jiali fault, supporting the interpretation of sub-horizontal underplating of Tibetan crust by Indian crust to ∼31° N at 85° E on the Hi-CLIMB transect. However, we also show a negative mantle converter sub-parallel to the crustal YZS converter, from the northern limit of the underplating Indian lower crust to at least 140 km depth, that we interpret as the base of Tibetan lithosphere overlying an asthenospheric mantle wedge. Based on the lithospheric structure observed in this and other studies, we infer that Indian mantle lithosphere currently detaches from Indian lower crust at the “mantle suture” that is nearly 50 km south of the surface trace of the YZS at 92° E, south of the mantle suture suggested by INDEPTH transect beneath the surface trace of the YZS at ∼90° E, and far south of the mantle suture suggested to be at the 31° N northern limit of underthrusting Indian lower crust suggested by Hi-CLIMB transect at ∼85° E. A change from underplating in the west to steep subduction in the east can reconcile all these observations.

Himalayan tectonic activity is triggered by downward penetration of the Indian plate beneath the Asian plate. The subsurface geometry of this interaction has not been fully investigated. This study presents novel constraints on this geometry provided by two newly obtained, deep seismic reflection profiles. The profiles cover 100- and 60-km transects across the Yarlung-Zangbo suture of the Himalaya-Tibet orogen at c. 88°E. Both profiles show a crustal-scale outline of the subducting Indian crust. This outline clearly shows Indian understhrusting southern Tibet, but only to a limited degree. When combined with a third seismic reflection profile of the western Himalayas, the new profiles reveal progressive, eastward steepening and shortening in the horizontal advance of the subducting Indian crust.

The Tibetan Plateau's enigmatic collapse, which was indicated based on extensional faulting in the interior of the plateau and large eastward-directed strike-slip faults, has inspired intensive geologic inquiry and debate. Interaction of the subducting Indian plate with the overlying Asian lithosphere is one factor that may be affecting the plateau. A more detailed image of the subducting Indian plate can help constrain its role, but regional high-resolution seismic data have not been widely available at a relevant scale. Here we present an integrated interpretation based on two types of seismic datasets: (1) two N-S oriented deep seismic-reflection profiles that cross the Yarlung-Zangbo suture, and (2) two E-W receiver function profiles that cross the Xainza-Dingjye and Nyima-Tingri rifts, respectively. These images reveal different crustal-scale structures in the dominant collision zone between the western and central regions, as well as distinctly offset Moho discontinuities across the N-S trending rift grabens in southern Tibet. These crustal variations, combined with previous tomographic studies in Tibet, outline an easterly tilt of the subducting Indian slab, along which crust-mantle decoupling occurred. Together with the spatio-temporal distribution of synchronous potassium-rich volcanics exposed within the grabens in southern Tibet, this offset Moho structure beneath the grabens lends support to the hypothesis that eastward steepening of the subducting Indian slab was accompanied by slab tearing beneath each north-south striking rift. This overall crustal geometric variation suggests a stepwise flattening of the torn Indian slab from west to the east, a process that brought easterly migration of gravitational instability to the overriding Tibetan crust that drove collapse of the Tibetan Plateau once gravitational instability reached a critical level to the east.

The northeastern Tibetan plateau results from the superposition of the Caledonian Qilian orogenic belt and Cenozoic reactivation following continental collision between the Indian and Asian plates. In order to identify the constraints on the Qilian orogenic mechanism and the expansion of the plateau, we carried out a 430-km-long wide-angle seismic experiment between Jingtai and Hezuo. We herein present the crustal P-wave velocity structure model resulting from the interpretation of this wide-angle seismic dataset. The principal characteristics of the crustal velocity model include: (1) thinning of the crust from south to north within a range of about 48–54 km with some undulation of the Moho beneath the Northern Qilian block; (2) the thickness of the sedimentary layer and its average P-wave velocity exhibit the obvious tectonic division of the Qaidam–Kunlun–West Qinling belt, the central and Northern Qilian, and the Alax blocks; (3) the lower crustal layer has a P-wave velocity of 6.6–6.8 km/s beneath the Qaidam–Kunlun–West Qinling belt, Northern Qilian orogenic belt and Alax block, but has a very low value of 6.4–6.5 km/s beneath the Central Qilian block; and (4) the Central Qilian block shares a similar crustal P-wave velocity–depth relationship with the Sierra Nevada arc or the accretionary crust. The Qaidam–Kunlun–West Qinling and the Northern Qilian orogenic belts have a crustal P-wave velocity–depth relationship characteristic of the crust of the global average continent. The particular Vp–depth relationship beneath the Central Qilian orogenic belt may result from a Caledonian accretionary orogeny or from the delamination of the ecologized lower crust during the arc–continental collisional orogeny. The lower velocity layer (LVL) exists at the bottom of the upper crust, to the north of the middle of the central Qilian block. We infer that this LVL plays an important role as an intracrustal decollement to expand the plateau.

The Tibetan Plateau is undergoing north–south shortening accompanied by west–east extension, as evidenced by the widespread development of north–south trending normal faults, grabens and rifts. While the mode of the north–south shortening has been the main focus of most international studies, knowledge of the deep structure beneath South Tibet is required for understanding the mechanism of the west–east extension. The onset of the north–south trending normal faulting is commonly taken as an indicator that the Tibetan Plateau was uplifted to a near-maximum elevation before entering a collapsing stage. Here we report on the receiver functions of a seismological experiment across the northern segment of the Yadong-Gulu Rift (YGR), one of the youngest rifts in South Tibet. The migrated receiver function images reveal that the YGR is a high-angle normal fault characterized by a 5-km Moho rise from its western to eastern flank, together with distinct differences in the crustal structure and intracrustal seismic conversion patterns between the two flanks. This highly asymmetric lithospheric structure suggests whole-crustal extension controlled by a simple/general shear rifting mechanism. This simple/general shear rifting in the YGR is attributed to an eastward (horizontal) shear at the base of the upper crust, as evidenced by the observed Tibetan GPS velocity field and our observation of shear wave splitting discrepancy among the upper crust, lower crust and lithospheric mantle. We propose that in the YGR, simple shear rifting accommodates the northward injection of the Indian lithosphere, which may suggest that the onset of the north–south normal faulting does not indicate gravitational collapse of the Tibetan lithosphere.

We use common conversion point (CCP) stacking of Ps receiver functions to image the crustal structure and Moho of the Garhwal Himalaya of India. Our seismic array of 21 broadband seismometers spanned the Himalayan thrust wedge at 79–80°E, between the Main Frontal Thrust and the South Tibet Detachment, in 2005–2006. Our CCP image shows the Main Himalayan Thrust (MHT), the detachment at the base of the Himalayan thrust wedge, with a flat-ramp-flat geometry. Seismic impedance contrasts inferred from geologic cross-sections in Garhwal imply a negative impedance contrast (velocity decreasing downward) for the upper flat, located beneath the Lower Himalaya, and a positive impedance contrast (velocity increasing downward) for the ramp, located beneath the surface trace of the Munsiari Thrust (or MCT-I). At the lower flat, located beneath the Higher Himalaya, spatially coincident measurements of very high electrical conductivities require the presence of free fluids, and we infer a negative impedance contrast on the MHT caused by ponding of these fluids beneath the detachment. Our seismic image indicates that the upper flat of the MHT is ∼ 10 km below sea level and dips north at ∼ 2 ° , connecting to a mid-crustal ramp which is ∼ 10 km high and dips at ∼ 16 ° . The lower flat is 20–25 km below sea level and dips at ∼ 4 ° . The Main Central Thrust (MCT) appears as a negative impedance contrast, dipping at ∼ 16 ° . The Moho is nearly horizontal at 35–45 km depth beneath the Sub-Himalaya and Lower Himalaya, deepening to 50 km or more beneath the Higher Himalaya. This depth is 10–25 km shallower than in the NW Indian Himalaya and 5–10 km shallower than in central Nepal, requiring significant along-strike variations in crustal thickness. The observed thickness of subducted Indian crust in Garwhal is 20–28 km.

[1] We use ambient-noise tomography to map regional differences in crustal Rayleigh-wave group velocities with periods of 8–40 s across north Tibet using the International Deep Profiling of Tibet and the Himalaya phase IV arrays (132 stations, deployed for 10–24 months). For periods of 8–24 s (sensitive to midcrustal depths of ~5–30 km), we observe striking velocity changes across the Bangong-Nujiang and Jinsha suture zones as well as the Kunlun-Qaidam boundary. From south to north, we see higher velocities beneath the Lhasa terrane, lower velocities beneath the Qiangtang, higher velocities in the Songpan-Ganzi and Kunlun terranes, and the lowest velocities beneath the Qaidam Basin. Maps at periods of 34 and 40 s (sensitive to the middle and lower crust at depths of ~30–60 km) do not show evidence of changes across those boundaries. Any differences between the Tibetan terrane lower crusts that were present at accretion have been erased or displaced by Cenozoic processes and replaced almost ubiquitously by uniformly low velocities.

The Karakoram fault (KKF) is the 1000 km-long strike-slip fault separating the western Himalaya from the Tibetan Plateau. From geologic and geodetic data, the KKF is argued either to be a lithospheric-scale fault with hundreds of km of offset at several cm/a, or to be almost inactive with cumulative offset of only a few tens of kilometers and to be just the upper-crustal localization of distributed deformation at depth. Here we show 3He/4He ratios in geothermal springs along a 500-km segment of the KKF are 3–100 times the normal ratio in continental crust, providing unequivocal evidence that a component of these hydrologic systems is derived from tectonically active mantle. Mantle enrichment is absent along the Indus–Yarlung suture zone (ISZ) just 35 km southwest of the KKF, suggesting that the mantle fluids flow only within the KKF. Within the last few Ma, the KKF must have accessed tectonically active Tibetan mantle northeast of the “mantle suture” which we therefore locate vertically beneath the KKF, very close to the surface trace of the ISZ. Hence, in southwestern Tibet, Indian crust may not now be underthrusting substantially north of the ISZ, even though Miocene underthrusting may have placed Indian crust north of the ISZ in the lower half of the Tibetan Plateau crust. This is in significant contrast to central and eastern Tibet where underthrust Indian material not only forms the lower half of the Tibetan crust but is also currently underthrusting for ∼200 km north of the ISZ. Our new constraint on KKF penetration to the mantle allows an improved description of the continuously evolving Karakoram fault, as a tectonically significant yet perhaps geologically ephemeral lithospheric structure.

SUMMARY From the S-wave data collected along a 270-km-long profile spanning the Kunlun mountains in NE Tibet, 14595 Sg phase arrivals and 21 SmS phase arrivals were utilized to derive a whole-crustal S velocity model and, together with a previously derived P velocity model, a Poisson's ratio (σ) model beneath the profile. The final tomogram for the upper 10–15 km of the crust reveals the lower velocities associated with the predominantly Neogene-Quaternary sediments of the Qaidam basin to the north and the higher velocities associated with the predominantly Palaeozoic and Mesozoic upper crustal sequences of the Songpan-Ganzi terrane and Kunlun mountains to the south. This study finds no evidence that the Kunlun mountains are involved in large-scale northward overriding of the Qaidam basin along a shallow south-dipping thrust. The σ in the upper 10–15 km of the crust are often lower than 0.25, indicating a preponderance of quartz-rich rocks in the upper crust beneath the profile. Below 10–15 km depth, the remainder of the crust down to the Moho has an average σ of 0.24 beneath the Songpan-Ganzi terrane and Kunlun mountains and 0.25 below the Qaidam basin. These low σ are similar to other low σ found along other profiles in the northeastern part of the plateau. Assuming an isotropic situation and no significant variation in σ between 10–15 km depth and the Moho, then the lower crust between 25–30 km depth below sea level and the Moho with P velocities varying from 6.6 km s−1 at the top to around 6.9 km s−1 at the base and σ of 0.24–0.25 should comprise intermediate granulites in the upper part transitioning to granulite facies metapelites in the lower part. As the pre-Cenozoic Qaidam basin crust has probably not lost any of its lower crust during the present Himalayan orogenic cycle in the Cenozoic and only has a σ of 0.245–0.25, then it appears that the pre-Cenozoic Qaidam basin crust involved in the collision is more felsic and thus weaker and more easily deformable than normal continental crust with a global average σ of 0.265–0.27 and the Tarim and Sichuan basin crusts. This situation then probably facilitates the collision and promotes the formation of new high plateau crust at the NE margin of Tibet. South of the Qaidam basin, the crust of the Songpan-Ganzi terrane and Kunlun mountains has an even lower average crustal σ of 0.23–0.24 and is thus presumably even weaker and more easily deformable than the crust beneath the Qaidam basin. This then supports the hypothesis of Karplus et al. that ‘the high Tibetan Plateau may be thickening northward into south Qaidam as its weak, thickened lower crust is injected beneath stronger Qaidam crust'.

Most compressional orogens include salients and recesses along their strike, the age and origin of which can be hard to ascertain. In the Kangra recess in the trace of the Main Boundary Thrust (MBT), the largest such recess within the active Himalayan orogen, the Sub-Himalayan sedimentary fold-thrust belt increases in width to as much as 90 km (the Kangra Basin), but narrows to as little as 10 km in the adjoining Nahan salient of the MBT (the Subathu Basin) to the southeast. New seismic reflection profiling places the Himalayan decollement at 6–8 km depth above thin but reflective Meso- to Neoproterozoic Vindhyan (Lesser Himalayan Series-equivalent) strata. These data show that the Vindhyan sedimentary rocks are thinner in the Kangra recess than further southeast, allowing the hypothesis that the width of the Lesser Himalayan thrust belt, and the existence of the Kangra recess, could be related to the pre-deformation basin thickness. This hypothesis obviates the need for control of the Kangra recess by a lateral ramp in the Main Himalayan Thrust, so making it more likely that the Kangra segment could rupture as part of an earthquake far larger than the devastating 1905 M = 7.8 Kangra earthquake. Below the Proterozoic sedimentary rocks, our reflection data show a west–southwest-dipping reflective fabric spanning a 30 km-crustal thickness, which we infer corresponds to a widespread “Ulleri-Wangtu” orogenic event at c. 1850 Ma affecting a pre-Tethyan Indian continental margin, thickening the basement by c. 20%. The deepest 10 km of this ~ 50 km-thick crust shows a more horizontal, arguably younger, reflectivity, though the Moho is not clearly marked by strong reflectors.

We herein present a new seismic refraction/wide-angle reflection profile that crosses the Songpan–Ganzi terrane, the Animaqing suture zone and the eastern Kunlun mountains (comprised of the South Kunlun and Middle Kunlun blocks separated by the Middle Kunlun fault). The profile is 380 km long and extends from Moba to Guide in eastern Tibet. The crustal thickness is about 62 km under the Songpan–Ganzi terrane, 62–64 km under the South Kunlun, and 60 km under the Middle Kunlun block. The Songpan–Ganzi flysch seems to be present up to a depth of 15 km south of the Animaqing suture zone, and up to a depth of 10 km in the Middle Kunlun block, with thicknesses elsewhere that depend on assumptions about the likely lithologies. The profile exhibits clear lateral variations both in the upper and lower crust, which are indicative of different crustal blocks juxtaposed by the Kunlun fault system. Whether or not the Songpan–Ganzi flysch was originally deposited on oceanic crust, at the longitude of our profile (100°E) it is now underlain by continental crust, and the presence of continental crust beneath the Songpan–Ganzi terrane and of a continental arc under the South Kunlun block suggest Paleozoic continent–continent arc collision in the eastern Kunlun Mountains. Comparison of crustal velocity columns from all wide-angle seismic profiles across the eastern Kunlun mountains indicates a remarkable west-to-east change in the Moho topography across the Kunlun fault system (15–20 km Moho step at 95°E, but only 2–5 km along our profile at 100°E). Lower-crustal thickness of the Kunlun terranes is rather uniform, about 35 km, from 80°–95°E, which suggests that similar thrust-thickening processes have played a role where the Qaidam Basin abuts the Kunlun fault, but thins to 20–25 km at 100°E, east of the Qaidam Basin. The increased crustal thickness from 93° to 98°E compared to that at 100°E may be due to the differences in the thickness of the crust of the two plates before their collision, and/or largely achieved by thickening of the lower crust, perhaps indicating a crustal flow mechanism operating more strongly in the western region.

Ground magnetic data contain information, not pre-sent in aeromagnetic data, which may be useful for precisely mapping near-surface faults and contacts, as well as constraining or aiding interpretation of other geophysical methods. However, collecting ground magnetic data on foot is labor-intensive and is therefore limited to small surveys. In this article, we present two newly developed all-terrain vehicle (ATV) magnetometer systems that significantly expand the survey area that is possible in a ground magnetic survey without greatly reducing the quality of data.

Seismic reflection profiling demonstrates that the active, significant-offset Surprise Valley Fault that marks the western boundary of the Basin-and-Range Province in northernmost California dips at a moderate angle, only ~ 30° to ~ 2 km depth. A nearby seismic refraction/wide-angle reflection profile, albeit of lower-resolution, shows the fault-plane remains approximately planar to at least 7 km depth. Variably-tilted volcanic strata of known ages within the Surprise Valley basin demonstrate that ≥ 7 km of normal slip occurred along the Surprise Valley Fault within the past 8–4 Myr at a time-averaged slip rate of ~ 1–2 mm/yr, and that the fault rotated from an initial dip of as much as 60° to its present dip of 30°. The rotation of this isolated fault was not accomplished by domino-style rotation as seen elsewhere in the Basin-and-Range province, but rather by flexural rift-shoulder uplift aided by mid-crustal flow.

article i nfo Seismic shear-wave velocities are sensitive to the partial melts that should be present in the Himalayan orogen if low-viscosity channel flow is active at the present day. We analyzed regional earthquakes in the western Himalaya and Tibet recorded on 16 broadband seismometers deployed across the NW Indian Himalaya, from the Indian platform to the Karakoram Range. We used a multiple filter technique to calculate the group velocity dispersion of fundamental-mode Rayleigh waves, and then inverted the dispersion records to obtain separate one-dimensional shear-wave velocity models for five geologic provinces: the Tibetan plateau, Ladakh arc complex, Indus Tsangpo suture zone, Tethyan Himalaya, and Himalayan thrust belt. Our velocity models show a low-velocity layer (LVL) with 7-17% velocity reduction centered at ~30 km depth and apparently continuous from the Tethyan Himalaya to the Tibetan plateau. This LVL shows good spatial correspondence with observations of low resistivity from magnetotelluric studies along the same profile. Of the possible explanations for low velocity and low resistivity in the mid-crust, only the presence of melts or aqueous fluids (or both) satisfactorily explains both sets of observations. Elevated heat flow observed in the NW Himalaya implies that if aqueous fluids are present in the mid-crust, then the mid-crust is well above its solidus. Comparison of our results with laboratory measurements and theoretical models suggests 3-7% melt is present in a channel in the upper-middle crust of the NW Himalaya at the present day, and the physical conditions to enable active channel flow may be present.

A test of deep seismic reflection profiling across the central uplift or metamorphic belt of the Qiangtang (羌塘) terrane, Tibetan plateau, provides a first image of the crustal structure. Complex reflection patterns in the upper crust are interpreted as a series of folds and thrusts, and bivergent reflections in the lower crust may represent a convergence between the Indian and the Eurasian plates.