skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Incorporating H‐ κ Stacking With Monte Carlo Joint Inversion of Multiple Seismic Observables: A Case Study for the Northwestern US
Abstract Accurately determining the seismic structure of the continental deep crust is crucial for understanding its geological evolution and continental dynamics in general. However, traditional tools such as surface waves often face challenges in solving the trade‐offs between elastic parameters and discontinuities. In this work, we present a new approach that combines two established inversion techniques, receiver function H‐κstacking and joint inversion of surface wave dispersion and receiver function waveforms, within a Bayesian Monte Carlo (MC) framework to address these challenges. Demonstrated by synthetic tests, the new method greatly reduces trade‐offs between critical parameters, such as the deep crustal Vs, Moho depth, and crustal Vp/Vs ratio. This eliminates the need for assumptions regarding crustal Vp/Vs ratios in joint inversion, leading to a more accurate outcome. Furthermore, it improves the precision of the upper mantle velocity structure by reducing its trade‐off with Moho depth. Additional notes on the sources of bias in the results are also included. Application of the new approach to USArray stations in the Northwestern US reveals consistency with previous studies and identifies new features. Notably, we find elevated Vp/Vs ratios in the crystalline crust of regions such as coastal Oregon, suggesting potential mafic composition or fluid presence. Shallower Moho depth in the Basin and Range indicates reduced crustal support to the elevation. The uppermost mantle Vs, averaging 5 km below Moho, aligns well with the Pn‐derived Moho temperature variations, offering the potential of using Vs as an additional constraint to Moho temperature and crustal thermal properties.  more » « less
Award ID(s):
2322632 2145410 1945856 2317868
PAR ID:
10519781
Author(s) / Creator(s):
; ;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
129
Issue:
7
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    SUMMARY Interfaces are important part of Earth’s layering structure. Here, we developed a new model parametrization and iterative linearized inversion method that determines 1-D crustal velocity structure using surface wave dispersion, teleseismic P-wave receiver functions and Ps and PmP traveltimes. Unlike previous joint inversion methods, the new model parametrization includes interface depths and layer Vp/Vs ratios so that smoothness constraint can be conveniently applied to velocities of individual layers without affecting the velocity discontinuity across the interfaces. It also allows adding interface-related observation such as traveltimes of Ps and PmP in the joint inversion to eliminate the trade-off between interface depth and Vp/Vs ratio and therefore to reduce the uncertainties of results. Numerical tests show that the method is computationally efficient and the inversion results are robust and independent of the initial model. Application of the method to a dense linear array across the Wabash Valley Seismic Zone (WVSZ) produced a high-resolution crustal image in this seismically active region. The results show a 51–55-km-thick crust with a mid-crustal interface at 14–17 km. The crustal Vp/Vs ratio varies from 1.69 to 1.90. There are three pillow-like, ∼100 km apart high-velocity bodies sitting at the base of the crust and directly above each of them are a low-velocity anomaly in the middle crust and a high-velocity anomaly in the upper crust. They are interpreted to be produced by mantle magmatic intrusions and remelting during rifting events in the end of the Precambrian. The current diffuse seismicity in the WVSZ might be rooted in this ancient distributed rifting structure. 
    more » « less
  2. Abstract This study provides a global assessment of the abundance of the major oxides in the deep continental crust. The combination of geochemistry and seismology better constrains the composition of the middle and lower continental crust better than either discipline can achieve alone. The inaccessible nature of the deep crust (typically >15 km) forces reliance on analog samples and modeling results to interpret its bulk composition, evolution, and physical properties. A common practice relates major oxide compositions of small‐ to medium‐scale samples (e.g., medium to high metamorphic grade terrains and xenoliths) to large scale measurements of seismic velocities (Vp, Vs, Vp/Vs) to determine the composition of the deep crust. We provide a framework for building crustal models with multidisciplinary constraints on composition. We present a global deep crustal model that documents compositional changes with depth and accounts for uncertainties in Moho depth, temperature, and physical and chemical properties. Our 3D compositional model of the deep crust uses the USGS Global Seismic Structure Catalog (Mooney, 2015) and a compilation of geochemical analyses on amphibolite and granulite facies lithologies (Sammon & McDonough, 2021,https://doi.org/10.1029/2021JB022791). We find a SiO2gradient from 61.2 ± 7.3 to 53.3 ± 4.8 wt.% from the middle to the base of the crust, with the equivalent lithological gradient ranging from quartz monzonite to gabbronorite. In addition, we calculate trace element abundances as a function of depth from their correlations with major oxides. From here, other lithospheric properties, such as Moho heat flux ( mW/m2), are derived. 
    more » « less
  3. Abstract In this study, we perform a 2‐frequency sequential receiver function stacking investigation in Southern California. The resulting Moho depths exhibit similar patterns to previous studies while the crystalline crustal Vp/Vs values show more regional variations. Most Vp/Vs variations can be explained by compositional differences. We observe a dichotomy in Moho depth, Vp/Vs, and crustal strain rates between the Peninsular Ranges and Southern San Andreas Fault system. Comparisons between strain rates, Vp/Vs, and temperature suggest that crustal compositional variations may have played a more critical role in influencing the crustal strain rate variations in the Peninsular Ranges and Southern San Andreas than temperature. The structural and compositional variations provide a new insight into the causes of the migration of the Southern San Andreas Fault system and the formation of the “Big Bend.” 
    more » « less
  4. Santorini volcano in the South Aegean Volcanic Arc has a detailed history of ongoing volcanic and seismic activity, making it a prime location for studying magma storage and transport at arc volcanoes. The shallow magmatic system (<5 km depth) is well constrained by geophysical studies, but the deeper crustal structure is not. Located 15 km NE of Santorini, the Kolumbo seamount is also an active edifice, with consistently more seismicity and hydrothermal venting than Santorini. Geochemical studies indicate that Santorini and Kolumbo are fed by separate mantle and crustal magma sources, but prior seismic studies suggest otherwise (Dimitriadis et al, 2010; McVey et al, 2020). This study addresses the nature of lower-crustal magma structure beneath arc volcanoes and whether evolved volcanoes and nearby vents are connected through their plumbing. Tomographic inversion of P-wave Moho reflection (PmP) and turning P-wave (Pg) traveltimes is used to create 3-D models of Moho depth and P-wave velocity (Vp) down to depths of ~25 km. The PROTEUS experiment provides an exceptionally dense and large aperture traveltime dataset from an amphibious array of ~150 seismometers and ~14,000 active marine sources. The data are ~33,000 manually picked PmP arrivals and ~256,000 Pg arrivals from existing studies. Results show a low Vp anomaly extending from the Moho to the surface. This anomaly starts at the base of the crust under the NW Santorini caldera and extends up to the east. It is most pronounced at 10-15 km depth, where it is offset from both Santorini and Kolumbo. Limited resolution prevents imaging of a connection between this mid-crustal anomaly and the known shallow magma storage region under the Santorini caldera. A high-velocity core beneath Santorini is not found, a feature interpreted at other volcanoes as a cooled intrusive complex. Because no additional low Vp anomalies are found in the lower crust, we infer that a common mantle source and mid-crustal plumbing system is actively feeding both Santorini and Kolumbo. The spatial offset and elongated nature of magma storage implies a complex relationship between evolving magmatic structures and tectonics. 
    more » « less
  5. Abstract Earthquakes are known to occur beneath southern Tibet at depths up to ∼95 km. Whether these earthquakes occur within the lower crust thickened in the Himalayan collision or in the mantle is a matter of current debate. Here we compare vertical travel paths expressed as delay times between S and P arrivals for local events to delay times of P-to-S conversions from the Moho in receiver functions. The method removes most of the uncertainty introduced in standard analysis from using velocity models for depth location and migration. We show that deep seismicity in southern Tibet is unequivocally located beneath the Moho in the mantle. Deep seismicity in continental lithosphere occurs under normally ductile conditions and has therefore garnered interest in whether its occurrence is due to particularly cold temperatures or whether other factors are causing embrittlement of ductile material. Eclogitization in the subducting Indian crust has been proposed as a cause for the deep seismicity in this area. Our observation of seismicity in the mantle, falling below rather than within the crustal layer with proposed eclogitization, requires revisiting this concept and favors other embrittlement mechanisms that operate within mantle material. 
    more » « less