We develop and apply a technique for receiver functions that is analogous to “velocity analysis” in seismic reflection processing, in which a velocity model is found directly from the data. In the case of receiver functions, which represent waves converted from
This study presents an improved approach to common‐conversion point stacking of converted body waves that incorporates scattering kernels, accurate and efficient measurement of stack uncertainties, and an alternative method for estimating free surface seismic velocities. To better separate waveforms into the
- Award ID(s):
- 1829401
- NSF-PAR ID:
- 10381601
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 125
- Issue:
- 11
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract P type toS type (Ps ) orS type toP type (Sp ) rather than reflected waves, we jointly optimize the average correlation value for target features in common receiver gathers while perturbing shear wave velocities and Poisson's ratio in a process driven by simulated annealing. The technique simultaneously finds depths to major discontinuities and 1‐DS velocity profiles beneath each seismic station. Estimates of Poisson's ratio are used withS velocity models to estimateP velocities, and individual 1‐D profiles are interpolated to form a 3‐D velocity model. We stacked and depth migratedPs andSp receiver functions computed from data recorded by broadband stations in southeastern New Mexico and west Texas. Images reveal an abrupt increase in lithospheric thickness beneath the Great Plains craton with respect to the Rio Grande Rift.Ps results show that the Moho ranges from 36 to 60 km in the region, whileSp results show variations in the seismically determined lithosphere‐asthenosphere boundary of 75–112 km. Thickened crust and lithosphere overlie a high‐velocity anomaly determined previously via tomography. We speculate that the thickened crust and lithosphere maybe due to removal of the lower lithosphere by a process that is associated with east‐west extension and thus the northward propagation of the Rio Grande Rift. -
Abstract The Colorado Plateau and its surroundings serve as an archetypal case to investigate the interaction of mantle melting processes and lithospheric structure. It has been hypothesized that widespread Cenozoic volcanism indicates the encroachment of the convective upwelling of asthenosphere toward the Plateau center. In this study, we generate a Common Conversion Point (CCP) stack of S‐to‐p (Sp) receiver functions to image the locations of lithospheric discontinuities in the southwestern United States. Our results are broadly similar to prior work, showing a strong and continuous Negative Velocity Gradient (NVG) consistent with the Lithosphere‐Asthenosphere Boundary (LAB) over much of the study area. However, with several methodological improvements, we are able to obtain more reliable NVG depth picks below the Colorado Plateau where the LAB becomes weaker, deeper, and broader. We compare the inferred topography of NVGs with the locations of volcanoes, and find that the majority of recent volcanoes are co‐located with lithosphere that is ∼80 km thick. This appears to be the critical depth at which partial melt from upwelling asthenosphere pooling at the base of (or within) the lithosphere may percolate to the surface. We compare our CCP profiles with magma equilibration conditions determined from petrologic analysis and find good agreement between the depth of NVGs and depth of magma equilibration. This analysis provides insight into the progression of magmatism and lithospheric loss toward the center of the Colorado Plateau, and demonstrates how small‐scale processes like melting influence lithosphere‐asthenosphere interactions that persist over large temporal and spatial scales.
-
With the ongoing discussion of Earth structure under West Antarctica and how it relates to the extension and volcanism of the area, we explore the possibility of a hydrated or thermally perturbed mantle underneath the region. Using P-wave receiver functions, we focus on the Mantle Transition Zone (MTZ) and how its thickness fluctuates from the global average (240-260 km). Prior studies have explored the West Antarctic regions of Marie Byrd Land and the West Antarctic Rift, but we expand this to include ~3-5 years of recent, additional seismic data from the Amundsen Sea and Pine Island Bay regions. Several years of additional data from the Ronne-Fichtner Ice Shelf, Ellsworth Land, and Marie Byrd Land regions will help provide a more complete picture of the mantle transition zone. Data for this study was obtained from IRIS for earthquakes of a 5.5 magnitude or greater. We use an iterative, time domain deconvolution method, filtered with Gaussian widths of 0.5, 0.75, and 1.0. All events within their respective Gaussian filter have undergone quality check by removing waveforms that have lower than 85% fit and visually checking for clear outliers. We migrate the receiver functions to depth and stack, using both single station stacking and common conversion point (CCP) stacking. We migrate the CCP stacks assuming both 1D (AK-135) and 3D velocity models throughout the region. Preliminary results from single-station stacks beneath the Thurston Island and Amundsen Sea regions suggest that the MTZ thickness is similar to the global average and the depth to the transition zone appears to be depressed, with average transition zone boundaries appearing around 430 and 680 km. If the MTZ is thinner than the global average, it may be an indication for high temperature thermal anomalies or a plume under West Antarctica that may help explain the history of extension and uplift there. These results could be useful for glacial isostatic adjustment and/or geothermal heat flux models that attempt to understand ice sheet history and stability.more » « less
-
Abstract Constraints on chemical heterogeneities in the upper mantle may be derived from studying the seismically observable impedance contrasts that they produce. Away from subduction zones, several causal mechanisms are possible to explain the intermittently observed X‐discontinuity (X) at 230–350 km depth: the coesite‐stishovite phase transition, the enstatite to clinoenstatite phase transition, and/or carbonated silicate melting, all requiring a local enrichment of basalt. Africa hosts a broad range of terranes, from Precambrian cores to Cenozoic hotspots with or without lowermost mantle origins. With the absence of subduction below the margins of the African plate for >0.5 Ga, Africa presents an ideal study locale to explore the origins of the X. Traditional receiver function (RF) approaches used to map seismic discontinuities, such as common conversion‐point stacking, ignore slowness information crucial for discriminating converted upper mantle phases from surface multiples. By manually assessing depth and slowness stacks for 1° radius overlapping bins, normalized vote mapping of RF stacks is used to robustly assess the spatial distribution of converted upper mantle phases. The X is mapped beneath Africa at 233–340 km depth, revealing patches of heterogeneity proximal to mantle upwellings in Afar, Canaries, Cape Verde, East Africa, Hoggar, and Réunion with further observations beneath Cameroon, Madagascar, and Morocco. There is a lack of an X beneath southern Africa and strikingly, the magmatic eastern rift branch of the southern East African Rift. With no relationships existing between depth and amplitudes of observed X and estimated mantle temperatures, multiple causal mechanisms are required across a range of continental geodynamic settings.
-
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.