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: Crustal Structure Beneath the Wabash Valley Seismic Zone From the Joint Inversion of Receiver Functions and Surface‐Wave Dispersion: Implications for Continental Rifts and Intraplate Seismicity
Abstract Seismicity of several intraplate seismic zones in the North American midcontinent is believed to be related to reactivation of ancient faults in Precambrian continental rifts by the contemporary stress field. Existence of such a rift system beneath the Wabash Valley Seismic Zone (WVSZ) is not clear. Here we obtained a crustal structural image along a 300‐km‐long profile across WVSZ using a dense linear seismic array. We first calculated teleseismic receiver functions of stations and applied the Common‐Conversion‐Point stacking method to image crustal interfaces and the Moho. We then used ambient noise cross correlation to obtain phase and group velocities of Rayleigh and Love waves. Finally, we jointly inverted the receiver function and surface wave dispersion data to determine shear wave velocity structure along the profile. The results show a thick (50‐ to 60‐km) crust with a typical Proterozoic crustal layering: a 1‐ to 2‐km thick Phanerozoic sedimentary layer, an upper crust ∼15 km thick, and a 30‐ to 40‐km‐thick lower crust. The unprecedented high‐resolution image also reveals a 50‐km‐wide high‐velocity body above an uplifted Moho and several velocity anomalies in the upper and middle crust beneath the La Salle Deformation Belt. We interpreted them as features produced by magmatic intrusions in a failed, immature continental rift during the end of Precambrian. Current seismicity in WVSZ is likely due to reactivation of ancient faults of the rift system by a combination of stress fields from the far‐field plate motion and prominent crustal and upper mantle heterogeneities in the region.  more » « less
Award ID(s):
1661519
PAR ID:
10371843
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
124
Issue:
7
ISSN:
2169-9313
Page Range / eLocation ID:
p. 7028-7039
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. The New England Appalachians provide a fascinating window into a host of fundamental geological problems. These include the modification of crustal and mantle lithospheric structure via orogenesis, terrane accretion, and continental rifting, the evolution of individual terranes through processes such as channel flow and ductile extrusion, and the causes and consequences of the Northern Appalachian Anomaly (NAA), a prominent geophysical anomaly in the upper mantle. Recent and ongoing deployments of dense seismic arrays in New England are providing images of the crust and upper mantle in unprecedented detail, allowing us to address both new and longstanding science questions. These deployments include the Seismic Experiment for Imaging Structure beneath Connecticut (SEISConn, 2015-2019), the New England Seismic Transects (NEST, 2018-present), and the GEology of New England via Seismic Imaging Studies (GENESIS, 2022-present) arrays. Here we present results from these experiments that are shedding new light on the tectonic evolution of New England and the ways in which structures and processes in the upper mantle can affect the structure of the overlying lithosphere. These include detailed new images of crustal architecture beneath central and southern New England, including a sharp transition from thick (~48 km) crust Laurentia terranes to thin (~32 km) crust beneath Appalachian terranes. The character of this offset beneath the SEISConn and NEST arrays suggests an overlap of two Moho boundaries, forming an overthrust-type structure that may have resulted from reactivation of faults during the compression and shortening associated with the formation of the hypothesized Acadian Altiplano. Beneath SEISConn, there is evidence for multiple relict structures preserved in the lithosphere from past episodes of terrane accretion and suturing, as well as anisotropic layering that constrains the kinematics of past lithospheric deformation events. Beneath the NEST line in central New England, we infer a relatively shallow (~80 km) lithosphere-asthenosphere boundary above the NAA upper mantle geophysical anomaly, providing evidence for lithospheric thinning above a presumed asthenospheric upwelling. Finally, preliminary results suggest layered crustal anisotropy beneath the GENESIS array, perhaps corresponding to a past episode of channel flow in the mid-crust. 
    more » « less
  3. The New England Appalachians provide a fascinating window into a host of fundamental geological problems. These include the modification of crustal and mantle lithospheric structure via orogenesis, terrane accretion, and continental rifting, the evolution of individual terranes through processes such as channel flow and ductile extrusion, and the causes and consequences of the Northern Appalachian Anomaly (NAA), a prominent geophysical anomaly in the upper mantle. Recent and ongoing deployments of dense seismic arrays in New England are providing images of the crust and upper mantle in unprecedented detail, allowing us to address both new and longstanding science questions. These deployments include the Seismic Experiment for Imaging Structure beneath Connecticut (SEISConn, 2015-2019), the New England Seismic Transects (NEST, 2018-present), and the GEology of New England via Seismic Imaging Studies (GENESIS, 2022-present) arrays. Here we present results from these experiments that are shedding new light on the tectonic evolution of New England and the ways in which structures and processes in the upper mantle can affect the structure of the overlying lithosphere. These include detailed new images of crustal architecture beneath central and southern New England, including a sharp transition from thick (~48 km) crust Laurentia terranes to thin (~32 km) crust beneath Appalachian terranes. The character of this offset beneath the SEISConn and NEST arrays suggests an overlap of two Moho boundaries, forming an overthrust-type structure that may have resulted from reactivation of faults during the compression and shortening associated with the formation of the hypothesized Acadian Altiplano. Beneath SEISConn, there is evidence for multiple relict structures preserved in the lithosphere from past episodes of terrane accretion and suturing, as well as anisotropic layering that constrains the kinematics of past lithospheric deformation events. Beneath the NEST line in central New England, we infer a relatively shallow (~80 km) lithosphere-asthenosphere boundary above the NAA upper mantle geophysical anomaly, providing evidence for lithospheric thinning above a presumed asthenospheric upwelling. Finally, preliminary results suggest layered crustal anisotropy beneath the GENESIS array, perhaps corresponding to a past episode of channel flow in the mid-crust. 
    more » « less
  4. 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
  5. null (Ed.)
    Recent GPS studies show that the Indo-Burma subduction system is locked with the implication of a potential large-magnitude earthquake. To inform better seismic hazard models in the region, we need an improved understanding of the crustal structure and the dynamics of the Indo-Burma subduction system. The Bangladesh-India-Myanmar (BIMA) tripartite project deployed 60 broadband seismometers across the subduction system and have been continuously recording data for ~2 years. In this study, we computed receiver functions from 30 high-quality earthquakes (M≥5.9) with epicentral distances between 30º and 90º recorded by the array. The algorithm utilized ensures the uniqueness of the seismic model and provides an uncertainty estimate of every converted wave amplitude. We stacked all the receiver functions produced at each station along the entire transect to generate a cross-sectional model of the average crustal structure. The level of detail in the image is improved by computing higher frequency receiver functions up to 4 Hz. The results represent some of the strongest constraints on crustal structure across the subduction system. Beneath the Neogene accretionary prism's outer belt, we observe a primary conversion associated with the Ganges Brahmaputra Delta that ranges in depth from ~10 km near the deformation front up to ~12 km at the eastern boundary. From the eastern end of the Neogene accretionary prism to the Sagaing Fault, we image the Indian subducting slab and the Central Myanmar basin. The depth-extent of seismicity associated with the Wadati-Benioff zone is consistent with the locations of primary conversions from the subducting plate. We further verify the converted phases of the slab by analyzing azimuthal moveout variations. The Central Myanmar basin is roughly bowl-shaped in cross-section with a maximum thickness of ~15 km about halfway between the Kabaw and Sagaing faults. The average crustal thickness beneath the Ganges-Brahmaputra delta is ~20 km, most likely representing a transitional crust formed from thinning of the continental crust intruded and underplated by igneous rocks. In contrast, the average thickness of the continental crust beneath the Central Myanmar basin is ~40 km. Our results provide a baseline model for future geophysical investigations of the Indo-Burma subduction zone. 
    more » « less