Fine-scale Southern California Moho structure uncovered with distributed acoustic sensing

Moho topography yields insights into the evolution of the lithosphere and the strength of the lower crust. The Moho reflected phase (PmP) samples this key boundary and may be used in concert with the first arriving P phase to infer crustal thickness. The densely sampled station coverage of distributed acoustic sensing arrays allows for the observation of PmP at fine-scale intervals over many kilometers with individual events. We use PmP recorded by a 100-km-long fiber that traverses a path between Ridgecrest, CA and Barstow, CA to explore Moho variability in Southern California. With hundreds of well-recorded events, we verify that PmP is observable and develop a technique to identify and pick P-PmP differential times with high confidence. We use these observations to constrain Moho depth throughout Southern California, and we find that short-wavelength variability in crustal thickness is abundant, with sharp changes across the Garlock Fault and Coso Volcanic Field.

Crustal melting and continent uplift by mafic underplating at convergent boundaries

The thick crust of the southern Tibetan and central Andean plateaus includes high-conductivity, low-velocity zones ascribed to partial melt. The melt origin and effect on plateau uplift remain speculative, in particular if plateau uplift happens before continental collision. The East Anatolian Plateau (EAP) has experienced similar, more recent uplift but its structure is largely unknown. Here we present an 80 km deep geophysical model across EAP, constrained by seismic receiver functions integrated with interpretation of gravity data and seismic tomographic, magnetotelluric, geothermal, and geochemical models. The results indicate a 20 km thick lower crustal layer and a 10 km thick mid-crustal layer, which both contain pockets of partial melt. We explain plateau uplift by isostatic equilibration following magmatism associated with roll-back and break-off of the Neo-Tethys slab. Our results suggest that crustal thickening by felsic melt and mafic underplate are important for plateau uplift in the EAP, Andes and Tibet.

Seismic mapping of the central and southern segments of the Tanlu fault zone using P-wave receiver functions

The central and southern segments of the Tanlu fault zone, located in the collision boundary between the Yangtze Plate and the North China Craton, underwent complex tectonic deformation associated with the Pacific Plate subduction. The crustal thickness and Vp/Vs ratio are important parameters for comprehending tectonic evolution and geodynamic processes. By integrating a newly dense seismic array, our results based on P-wave receiver function analyses reveal the crust thickness varies significantly across the Tanlu fault zone. The Yangtze Plate is characterized by a thinner crust, contrasting a thicker crust imaged in the North China Craton, indicating the Tanlu fault zone is a prominent structural block boundary. The Vp/Vs ratio within the Tanlu fault zone is visibly higher than its surroundings, probably correlating with abundant fluids activity along the fault zone, as a result of the subduction of the Pacific Plate. Additionally, we observe a gradually deepening Moho discontinuity beneath the Sulu orogenic belt with a low-angle dip, presenting direct seismic evidence for supporting the ancient Yangtze Plate underthrusting beneath the Sulu orogenic belt.

Injection-induced basement seismicity beneath the Raton Basin: constraints from refined fault architectures and basin structure

Intraplate earthquakes induced by anthropogenic fluid injection present unexpected seismic risk to previously quiescent or low seismicity-rate regions. Despite many studies of induced seismicity, there are relatively few with detailed openly accessible constraints on the interaction between seismic sources and subsurface structures. In this study of the Raton Basin, we refine source observations from a dense nodal array and constrain basin structure using teleseismic receiver functions. The cross-correlation-based relocated hypocentres and a new set of focal mechanisms light up active fault segments and show clear spatiotemporal patterns. The geometric complexity of reactivated fault clusters appears greatest near higher rate injection wells. Simpler normal fault structure is found farther from injection wells and near abrupt structural transitions suggested by receiver functions. While less induced seismicity in the crystalline basement is expected when injection is >1 km from the top of the basement (like Raton), our receiver function analysis identified a basin thickness ~3 km beneath the nodal array and lateral variations in sedimentary structures. Our results explain potential fluid connectivity between the injection depths focused at ~1-1.5 km below the surface and basement fault activity that begins at ~3 km and reaches peak activity at ~4-8 km depths. This article is part of the theme issue 'Induced seismicity in coupled subsurface systems'.

Moho Imaging with Fiber Borehole Strainmeters Based on Ambient Noise Autocorrelation

Moho tomography is important for studying the deep Earth structure and geodynamics, and fiber borehole strainmeters are broadband, low-noise, and attractive tools for seismic observation. Recently, many studies have shown that fiber optic seismic sensors can be used for subsurface structure imaging based on ambient noise cross-correlation, similar to conventional geophones. However, this array-dependent cross-correlation method is not suitable for fiber borehole strainmeters. Here, we developed a Moho imaging scheme for the characteristics of fiber borehole strainmeters based on ambient noise autocorrelation. S-wave reflection signals were extracted from the ambient noise through a series of processing steps, including phase autocorrelation (PAC), phase-weighted stacking (PWS), etc. Subsequently, the time-to-depth conversion crustal thickness beneath the station was calculated. We applied our scheme to continuous four-component recordings from four fiber borehole strainmeters in Lu'an, Anhui Province, China. The obtained Moho depth was consistent with the previous research results. Our work shows that this method is suitable for Moho imaging with fiber borehole strainmeters without relying on the number of stations.

Fault-scale crustal structure across the Dunhua-Mishan fault (Tanlu northern segment) constrained from teleseismic P-wave receiver functions

The Dunhua-Mishan fault, located in the northern segment of the Tanlu fault zone, experienced multiple tectonic processes associated with the effects of the Pacific Plate subduction and the Indo-Asia collision. The high-resolution fault-scale structure is critical for understanding the fault evolution and potential fault damage. However, the well-defined deep structure of the Dunhua-Mishan fault is still unclear due to the lack of the dense seismic array. In this study, we construct a high-resolution P-wave receiver function imaging based on linear dense seismic array across the fault. Our results reveal the strong Moho depth variation across the Dunhua-Mishan fault zone. The slightly higher Vp/Vs ratio values within the fault zone indicate the presence of a small amount of mafic crust composition. Interestingly, the significant double positive Ps converted phases are observed within the fault zone, which may represent double Moho discontinuities. The double Moho structure may be related to multiple significant tectonic activities in the Tanlu northern segment. These newly observed structures provide new seismic constraints on the formation and evolution of the Tanlu fault zone and probably reflect that the lithospheric structure of the Dunhua-Mishan fault has been modified by a series of tectonic processes.

Least squares collocation method in Moho depth determination in Iran using gravity gradient data

In this study, an approach using gravity observations utilizing the Least Squares Collocation (LSC) method is developed with the aim of mapping the depth spatial variability of the Mohorovičić discontinuity. This approach is based on a spherical two-layer isostatic model where the exterior gravity field only varies because of the shifting topographic masses and the related isostatic adjustment since it is believed that the Earth's core has a uniform density distribution. Assuming mass conservation between the Moho column of height δ R with respect to R m representing the mean Moho and following a Helmert condensation approach, the relationship between the surface layer density to the potential δ T can be obtained and δ R can be estimated via LSC from observed values of any functional derived from δ T . With such approach, the depth of Moho in the Iranian Plateau is estimated from T r r data generated by GOCO06S model reduced by topography, bathymetry and sediments effects by considering GEBCO2021 and CRUST1.0 models. The needed a-priori assumptions on R m and the density contrast Δ ρ are tuned so to obtain the best fit with seismic Moho depths reported by literature. 73 stations were matched with 3 km of standard deviation, which is coherent with the expected accuracy of the benchmark values. The remaining greater discrepancies showed to be clustered in defined areas like the Zagros chain and the reliefs along the Caspian coastline and the East borders.

Crustal structure and seismic anisotropy of rift basins in Somaliland

Rift margins provide insights into the processes governing the rupture of the continental lithosphere and the subsequence formation of sedimentary basins. The Proterozoic basement underlying Somaliland has been affected by multiple rifting; however, the crustal structure of these rifted basins remains unknown. This study utilized teleseismic receiver function analysis, Bayesian inversion, common conversion point imaging and 2D forward gravity modelling to examine the crust and upper mantle of Somaliland. The results indicate 36.8-38.2 km of crust in southern Somaliland, while the central and northern regions feature thinned crust (~ 21 km) with 5-6 km thick sediments. The joint analysis of radial and transverse components of receiver functions and shear wave splitting revealed fast axis directions trending to 50-56° in the upper mantle, indicating that azimuthal anisotropy is oriented in the regional Africa-Arabia plate motion. Such orientation may have resulted from lattice preferred orientation of olivine from the asthenospheric flow. Additionally, the fast polarization of the crust in central Somaliland is oriented at - 15°, indicating fossil deformation in the thinned crust related to the NW-SE trending Late Jurassic rift event. Further, the fast polarization for stations near the Gulf of Aden is oriented at 75-80°, suggesting crustal deformation associated with the Oligocene rift event. The crustal anisotropy at southern Somaliland revealed fast polarization oriented at - 85°, indicating a preserved far-field response of the WNW-ESE trending Late Cretaceous rift event. Overall, the study provides for the first-time insight into the rift-related extensional strain fabric in the crust and upper mantle anisotropy induced by asthenospheric flow in Somaliland.

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 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.

Receiver Function Imaging of the Crustal Structure Beneath Northern Taiwan Using Dense Linear Arrays

In order to realize the crustal structure in Taiwan, the receiver function method was used to analyze the teleseismic waveforms recorded by two orthogonal broadband linear arrays deployed in northern Taiwan in the east–west and south–north directions by the TAiwan Integrated GEodynamics Research (TAIGER) project from 2007 to 2009. By incorporation with Common Conversion Point (CCP) stacking, the Moho discontinuities beneath northern Taiwan were imaged. Based on the CCP stack of receiver functions in the east–west direction array, a collision boundary between the Philippine Sea Plate and the Eurasian Plate appears at the east of Taiwan. The Moho depth of the Eurasian Plate in west Taiwan is flat and 30 km; the Moho depth of the Philippine Sea Plate below the Central Mountain Range is about 55 km; in the east, the Moho depth of the Ryukyu arc is about 40 km. The south–north profile shows a brittle–ductile transition zone at depths of 15–20 km beneath central Taiwan from south to north. Moreover, the Moho depth of the Eurasian Plate is about 20–25 km in northern Taiwan. The Moho depth appears to deepen from north to south. The deepest Moho is located at the junction of the two profile lines, the Philippine Sea Plate, and has a depth of 60 km. According to these Moho depths, the crustal structure is thin and flat in the west part of northern Taiwan which is similar to the thin-skin model. However, the lithosphere is deformed and forms the mountain root in the east part which is similar to the lithospheric collision model.

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.

Modelling of crustal composition and Moho depths and their Implications toward seismogenesis in the Kumaon-Garhwal Himalaya

We image the lateral variations in the Moho depths and average crustal composition across the Kumaon–Garhwal (KG) Himalaya, through the H–K stacking of 1400 radial PRFs from 42 three-component broadband stations. The modelled Moho depth, average crustal Vp/Vs, and Poisson’s ratio estimates vary from 28.3 to 52.9 km, 1.59 to 2.13 and 0.17 to 0.36, respectively, in the KG Himalaya. We map three NS to NNE trending transverse zones of significant thinning of mafic crust, which are interspaced by zones of thickening of felsic crust. These mapped transverse zones bend toward the north to form a NE dipping zone of maximum changes in Moho depths, below the region between Munsiari and Vaikrita thrusts. The 1991 M_w6.6 Uttarakashi and 1999 M_w6.4 Chamoli earthquakes have occurred on the main Himalayan thrust (MHT), lying just above the mapped zone of maximum changes in Moho depths. Modelled large values of average crustal Vp/Vs (> 1.85) could be attributed to the high fluid (metamorphic fluids) pressure associated with the mid-crustal MHT. Additionally, the serpentinization of the lowermost crust resulted from the continent–continent Himalayan collision process could also contribute to the increase of the average crustal Vp/Vs ratio in the region.

Structural Characteristics of Moho Surface Based on Time Series Function of Natural Earthquakes

Remote sensing is a non-contact, long-distance detection technology. The reflection characteristics of a seismic wave can be used to detect remote and non-contact targets. Based on the reflection characteristics of a seismic wave, the underground structure in Tengchong Volcanic Area is explored. In order to further study the deep structure and magmatic activity of the crust in the volcanic area, we carried out a one-year mobile seismic observation. In this paper, nine broadband seismic stations were set up in the Tengchong Volcanic Area, and 3350 receiver function waveforms were collected. The crustal thickness, average wave velocity ratio, and Poisson’s ratio below these stations were calculated by the receiver function method, and the velocity structure near the Moho below these stations was evaluated. Combined with topographic data from SRTM3 (Shuttle Radar Topography Mission 3), this study reveals the dynamic relationship among crustal structure, crustal magmatism, and regional tectonic movement. Mantle upwelling plays an important role on the Moho uplift in the northern Tengchong Volcanic Area, and there are interconnected intracrustal magma chambers in the upper platform. The evaluation results of the Moho transition zone also indicate that the Dayingjiang fault is closely related to the tectonic activity of the Tengchong Volcanic fault.

Geophysical imaging of ophiolite structure in the United Arab Emirates

The Oman-United Arab Emirates ophiolite has been used extensively to document the geological processes that form oceanic crust. The geometry of the ophiolite, its extension into the Gulf of Oman, and the nature of the crust that underlies it are, however, unknown. Here, we show the ophiolite forms a high velocity, high density, >15 km thick east-dipping body that during emplacement flexed down a previously rifted continental margin thereby contributing to subsidence of flanking sedimentary basins. The western limit of the ophiolite is defined onshore by the Semail thrust while the eastern limit extends several km offshore, where it is defined seismically by a ~40-45°, east-dipping, normal fault. The fault is interpreted as the southwestern margin of an incipient suture zone that separates the Arabian plate from in situ Gulf of Oman oceanic crust and mantle presently subducting northwards beneath the Eurasian plate along the Makran trench.

Aseismic mid-crustal magma reservoir at Cleveland Volcano imaged through novel receiver function analyses

Processes related to eruptions at arc volcanoes are linked by structures that transect the entire crust. Imaging the mid- to lower-crustal portions (here, ~5–15 km and >15 km respectively) of these magmatic systems where intermediate storage may occur has been a longstanding challenge. Tomography, local seismic source studies, geodetic, and geochemical constraints, are typically most sensitive to shallow (<5 km) storage and/or have insufficient resolution at these depths. Geophysical methods are even further limited at frequently-erupting volcanoes where well-developed trans-crustal magmatic systems are likely to exist, due to a lack of deep seismicity. Here we show direct evidence for mid-crustal magma storage beneath the frequently erupting Cleveland volcano, Alaska, using a novel application of seismic receiver functions. We use P-s scattered waves from the Moho as virtual sources to investigate S-wave velocities between the Moho and the surface. Our forward modeling approach allows us to provide direct constraints on the geometry of low velocity regions beneath volcanoes despite having a comparatively sparse seismic network. Our results show clear evidence of mid-crustal magma storage beneath the depths of located volcanic seismicity. Future work using similar approaches will enable an unprecedented comparative examination of magmatic systems beneath sparsely instrumented volcanoes globally.

Water, Hydrous Melting, and Teleseismic Signature of the Mantle Transition Zone

Recent geophysical and petrological observations indicate the presence of water and hydrous melts in and around the mantle transition zone (MTZ), for example, prominent low-velocity zones detected by seismological methods. Experimental data and computational predictions describe the influence of water on elastic properties of mantle minerals. Using thermodynamic relationships and published databases, we calculated seismic velocities and densities of mantle rocks in and around the MTZ in the presence of water for a plausible range of mantle potential temperatures. We then computed synthetic receiver functions to explore the influence of different water distribution patterns on the teleseismic signature. The results may improve our understanding and interpretation of seismic observations of the MTZ.

Southern Africa crustal anisotropy reveals coupled crust-mantle evolution for over 2 billion years

The long-term stability of Precambrian continental lithosphere depends on the rheology of the lithospheric mantle as well as the coupling between crust and mantle lithosphere, which may be inferred by seismic anisotropy. Anisotropy has never been detected in cratonic crust. Anisotropy in southern Africa, detected by the seismological SKS-splitting method, usually is attributed to the mantle due to asthenospheric flow or frozen-in features of the lithosphere. However, SKS-splitting cannot distinguish between anisotropy in the crust and the mantle. We observe strong seismic anisotropy in the crust of southern African cratons by Receiver Function analysis. Fast axes are uniform within tectonic units and parallel to SKS axes, orogenic strike in the Limpopo and Cape fold belts, and the strike of major dyke swarms. Parallel fast axes in the crust and mantle indicate coupled crust-mantle evolution for more than 2 billion years with implications for strong rheology of the lithosphere.

The long-term stability of Precambrian continental lithosphere depends on the rheology of the lithospheric mantle as well as the coupling between crust and mantle lithosphere. Here, the authors study crustal seismic anisotropy to reveal that the crust and mantle lithosphere of southern Africa have been coupled since cratonisation, over 2 billion years ago.

Combined Gravimetric-Seismic Moho Model of Tibet

Substantial progress has been achieved over the last four decades to better understand a deep structure in the Himalayas and Tibet. Nevertheless, the remoteness of this part of the world still considerably limits the use of seismic data. A possible way to overcome this practical restriction partially is to use products from the Earth’s satellite observation systems. Global topographic data are provided by the Shuttle Radar Topography Mission (SRTM). Global gravitational models have been derived from observables delivered by the gravity-dedicated satellite missions, such as the Gravity Recovery and Climate Experiment (GRACE) and the Gravity field and steady-state Ocean Circulation Explorer (GOCE). Optimally, the topographic and gravity data should be combined with available results from tomographic surveys to interpret the lithospheric structure, including also a Moho relief. In this study, we use seismic, gravity, and topographic data to estimate the Moho depth under orogenic structures of the Himalayas and Tibet. The combined Moho model is computed based on solving the Vening Meinesz–Moritz (VMM) inverse problem of isostasy, while incorporating seismic data to constrain the gravimetric solution. The result of the combined gravimetric-seismic data analysis exhibits an anticipated more detailed structure of the Moho geometry when compared to the solution obtained merely from seismic data. This is especially evident over regions with sparse seismic data coverage. The newly-determined combined Moho model of Tibet shows a typical contrast between a thick crustal structure of orogenic formations compared to a thinner crust of continental basins. The Moho depth under most of the Himalayas and the Tibetan Plateau is typically within 60–70 km. The maximum Moho deepening of ~76 km occurs to the south of the Bangong-Nujiang suture under the Lhasa terrane. Local maxima of the Moho depth to ~74 km are also found beneath Taksha at the Karakoram fault. This Moho pattern generally agrees with the findings from existing gravimetric and seismic studies, but some inconsistencies are also identified and discussed in this study.

Probing layered arc crust in the Lesser Antilles using receiver functions

Oceanic arcs can provide insight into the processes of crustal growth and crustal structure. In this work, changes in crustal thickness and composition along the Lesser Antilles Arc (LAA) are analysed at 10 islands using receiver function (RF) inversions that combine seismological data with vP/vS ratios estimated based on crustal lithology. We collected seismic data from various regional networks to ensure station coverage for every major island in the LAA from Saba in the north to Grenada in the south. RFs show the subsurface response of an incoming signal assuming horizontal layering, where phase conversions highlight discontinuities beneath a station. In most regions of the Earth, the Mohorovičić discontinuity (Moho) is seismically stronger than other crustal discontinuities. However, in the LAA we observe an unusually strong along-arc variation in depth of the strongest discontinuity, which is difficult to explain by variations in crustal thickness. Instead, these results suggest that in layered crust, especially where other discontinuities have a stronger seismic contrast than the Moho, H-k stacking results can be easily misinterpreted. To circumvent this problem, an inversion modelling approach is introduced to investigate the crustal structure in more detail by building a one-dimensional velocity-depth profile for each island. Using this method, it is possible to identify any mid-crustal discontinuity in addition to the Moho. Our results show a mid-crustal discontinuity at about 10-25 km depth along the arc, with slightly deeper values in the north (Montserrat to Saba). In general, the depth of the Moho shows the same pattern with values of around 25 km (Grenada) to 35 km in the north. The results suggest differences in magmatic H2O content and differentiation history of each island.

Evidence for partial melt in the crust beneath Mt. Paektu (Changbaishan), Democratic People's Republic of Korea and China

Mt. Paektu (also known as Changbaishan) is an enigmatic volcano on the border between the Democratic People's Republic of Korea (DPRK) and China. Despite being responsible for one of the largest eruptions in history, comparatively little is known about its magmatic evolution, geochronology, or underlying structure. We present receiver function results from an unprecedented seismic deployment in the DPRK. These are the first estimates of the crustal structure on the DPRK side of the volcano and, indeed, for anywhere beneath the DPRK. The crust 60 km from the volcano has a thickness of 35 km and a bulk V P/V S of 1.76, similar to that of the Sino-Korean craton. The V P/V S ratio increases ~20 km from the volcano, rising to >1.87 directly beneath the volcano. This shows that a large region of the crust has been modified by magmatism associated with the volcanism. Such high values of V P/V S suggest that partial melt is present in the crust beneath Mt. Paektu. This region of melt represents a potential source for magmas erupted in the last few thousand years and may be associated with an episode of volcanic unrest observed between 2002 and 2005.

Seismic Evidence for a Geosuture between the Yangtze and Cathaysia Blocks, South China

South China, composed of the Yangtze and Cathaysia Blocks and the intervening Jiangnan orogenic belt, has been central to the debate on the tectonic evolution of East Asia. Here we investigate the crustal structure and composition of South China from seismic data employing the H-k stacking technique. Our results show that the composition and seismic structure of the crust in the Jiangnan orogenic belt are identical to those of the Cathaysia Block. Our data reveal a distinct contrast in the crustal structure and composition between the two flanks of the Jiujiang-Shitai buried fault. We propose that the Jiujiang-Shitai buried fault defines a geosuture between the Yangtze and Cathaysia Blocks, and that the felsic lower crust of the Cathaysia Block and the Jiangnan orogenic belt may represent fragments derived from the Gondwana supercontinent.

The role of crustal quartz in controlling Cordilleran deformation

Large-scale deformation of continents remains poorly understood more than 40 years after the plate tectonic revolution. Rock flow strength and mass density variations both contribute to stress, so both are certain to be important, but these depend (somewhat nebulously) on rock type, temperature and whether or not unbound water is present. Hence, it is unclear precisely how Earth material properties translate to continental deformation zones ranging from tens to thousands of kilometres in width, why deforming zones are sometimes interspersed with non-deforming blocks and why large earthquakes occasionally rupture in otherwise stable continental interiors. An important clue comes from observations that mountain belts and rift zones cyclically form at the same locations despite separation across vast gulfs of time (dubbed the Wilson tectonic cycle), accompanied by inversion of extensional basins and reactivation of faults and other structures formed in previous deformation events. Here we show that the abundance of crustal quartz, the weakest mineral in continental rocks, may strongly condition continental temperature and deformation. We use EarthScope seismic receiver functions, gravity and surface heat flow measurements to estimate thickness and seismic velocity ratio, v(P)/v(S), of continental crust in the western United States. The ratio v(P)/v(S) is relatively insensitive to temperature but very sensitive to quartz abundance. Our results demonstrate a surprising correlation of low crustal v(P)/v(S) with both higher lithospheric temperature and deformation of the Cordillera, the mountainous region of the western United States. The most plausible explanation for the relationship to temperature is a robust dynamical feedback, in which ductile strain first localizes in relatively weak, quartz-rich crust, and then initiates processes that promote advective warming, hydration and further weakening. The feedback mechanism proposed here would not only explain stationarity and spatial distributions of deformation, but also lend insight into the timing and distribution of thermal uplift and observations of deep-derived fluids in springs.

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.