An official website of the United States government
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
Constraining the crustal composition of the continental U.S. using seismic observables
The composition of the crust is one of the most uncertain and controversial components of continental estimates due to (1) limited direct access and (2) inconsistent indirect assessments. Here we show that by combining high-resolution shear velocity (Vs) models with newly measured with newly measured ratio of compressional wave velocity (Vp) and Vs, or Vp/Vs ratio, for the crystalline crust, a 3-D composition (SiO2 wt%) model of the continental crust can be derived with uncertainty estimates. Comparing the model with local xenolith data shows consistency at mid and lower crustal depths. The spatial patterns in bulk SiO2 content correlate with major geological provinces, including the footprints of Cenozoic and Mesozoic mafic volcanism in the western U.S., and offer new insight into the composition and evolution of the continental U.S.
more »
« less
- PAR ID:
- 10410449
- Date Published:
- Journal Name:
- GSA Bulletin
- ISSN:
- 0016-7606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Near‐surface seismic velocity structure plays a critical role in ground motion amplification during large earthquakes. In particular, the local Vp/Vs ratio strongly influences the amplitude of Rayleigh waves. Previous studies have separately imaged 3D seismic velocity and Vp/Vs ratio at seismogenic depth, but lack regional coverage and/or fail to constrain the shallowest structure. Here, we combine three datasets with complementary sensitivity in a Bayesian joint inversion for shallow crustal shear velocity and near‐surface Vp/Vs ratio across Southern California. Receiver functions–including with an apparent delayed initial peak in sedimentary basins, and long considered a nuisance in receiver function imaging studies–highly correlate with short‐period Rayleigh wave ellipticity measurements and require the inclusion of a Vp/Vs parameter. The updated model includes near‐surface low shear velocity more in line with geotechnical layer estimates, and generally lower than expected Vp/Vs outside the basins suggesting widespread shallow fracturing and/or groundwater undersaturation.more » « less
-
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
-
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
-
Abstract The composition of the lower continental crust, as well as its formation, growth, and evolution, remains a fundamental subject to be understood. In this study, we carry out a comparative and integrative analysis of seismic tomographic models, teleseismic receiver function results, and Airy isostasy in order to investigate the properties of the lower continental crust in eastern North America. We extract the depths for Vs of 4.0 km/s, 4.2 km/s, and 4.5 km/s from three selected tomographic models and calculate the differences between the Vs depth contours and the Moho depth defined by receiver functions. We then calculate the Airy isostatic Moho depth and its misfit with the receiver‐function‐defined Moho. Our analysis reveals three key features: (a) the deepening of the Vs depth contours and the strong negative Airy misfit within the U.S. Grenville Province; (b) a seismically faster‐than‐average and compositionally denser‐than‐average lowermost crust in the eastern North American Craton and the Grenville Province; and (c) the thickest, seismically fastest, and densest lowermost crust beneath the southern Grenville Front, the southern Grenville‐Appalachian boundary, and the U.S.‐Canada national border. We suggest that the lower crust of the craton and the Grenville Province has densified through garnet‐forming metamorphic reactions during and after orogenesis, contributing to the widely distributed fast‐velocity layer. The lower crust beneath the tectonic boundaries could have experienced more extensive garnet growth during orogenesis and emplacement of mafic magma. This study provides new constraints on the seismic and compositional properties of the lower crust in eastern North America.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1130/B36229.1
