The Mackenzie Mountains (MMs) in the Yukon and Northwest Territories, Canada, are an enigmatic mountain range. They are currently uplifting (Leonard et al., 2008,
Seismic tomography of the crust is an essential tool for studying the three‐dimensional structure of magmatic plumbing systems feeding active volcanoes, but it is often limited in resolution by the absence of deep local seismicity. Teleseismic receiver functions can be used to illuminate local structural variations, but typically do not account for the effects of three‐dimensional velocity heterogeneities. Here we harness the complementary strengths of both techniques by processing
- NSF-PAR ID:
- 10362191
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 47
- Issue:
- 22
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract https//doi.org/10.1029/2007JB005456 ), yet are about 700 km from the nearest plate boundary. Their arcuate shape is distinct and extends over 100 km eastward from the general trend of the Northern Canadian Cordillera. To better assess the cause and conditions of the current uplift, we processed ambient seismic noise data from a linear array of broadband seismographs crossing the mountains, along with other regional seismic stations, to estimate Rayleigh wave phase velocities between 6 and 40 s periods. From this, we estimated phase velocity dispersion and performed a tomographic inversion to estimateV S . Tomography reveals a low‐velocity structure that extends upward from the base of the ∼50–66 km thick lithosphere to the upper crust, and we hypothesize that inferred low density and low rigidity associated with theV S anomaly localizes the ongoing uplift and thrust‐dominated seismicity of the MMs. Additionally, we find relatively low crustal velocities that extend to the west of the MMs, suggesting that strain transfer from the Gulf of Alaska plate boundary plays a driving role as the crust translates to the northeast and buckles up against the craton consistent with the orogenic float hypothesis of Mazzotti and Hyndman (2002,https//doi.org/10.1130/0091-7613(2002)030〈0495:YCASTA〉2.0.CO;2 ). Finally, we observe lithospheric azimuthal anisotropy with an NW‐SE fast direction. This is nearly orthogonal to teleseismic shear wave splitting measurements in the central MMs, and suggests that asthenosphere flow and lithospheric strain are not aligned in this region. -
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
-
more » « less
Abstract We conduct a joint inversion of teleseismic receiver functions and Rayleigh wave phase velocity dispersion from both ambient noise and earthquakes using data from 79 seismic stations in southern Africa, which is home to some of the world's oldest cratons and orogenic belts. The area has experienced two of the largest igneous activities in the world (the Okavango dyke swarm and Bushveld mafic intrusion) and thus is an ideal locale for investigating continental formation and evolution. The resulting 3‐D shear wave velocities for the depth range of 0–100 km and crustal thickness measurements show a clear spatial correspondence with known geological features observed on the surface. Higher than normal mantle velocities found beneath the southern part of the Kaapvaal craton are consistent with the basalt removal model for the formation of cratonic lithosphere. In contrast, the Bushveld complex situated within the northern part of the craton is characterized by a thicker crust and higher crustal
V p /V s but lower mantle velocities, which are indicative of crustal underplating of mafic materials and lithospheric refertilization by the world's largest layered mafic igneous intrusion. The thickened crust and relatively low elevation observed in the Limpopo belt, which is a late Archean collisional zone between the Kaapvaal and Zimbabwe cratons, can be explained by eclogitization of the basaltic lower crust. The study also finds evidence for the presence of a stalled segment of oceanic lithosphere beneath the southern margin of the Proterozoic Namaqua‐Natal mobile belt. -
more » « less
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
Abstract Passive seismic methods for imaging the discontinuity structure of Earth have primarily focused on differences in vertically and radially polarized energy in the coda of earthquake‐generated body waves (e.g., receiver functions). To convert the timing of scattered wave arrivals to depth, three parameters must be known or inferred: depth or layer thickness (
H ),P ‐wave velocity (V P), andS ‐wave velocity (V S). A common way to solve for these parameters is throughH ‐κ stacking analysis, in which layer thickness and the ratio betweenV PandV S(κ ) is calculated while holding one of the velocity parameters constant. However, this assumption biases estimates of layer properties and leads to uncertainties that are not appropriately quantified. As these results are commonly used as starting models for more complex seismic or geodynamic analyses, these assumptions can propagate much further than the initial study. In this study, we introduce independent observations from body‐wave autocorrelations that can help constrain this underdetermined problem.P ‐wave autocorrelation allows for the recovery of the Moho‐reflectedP ‐wave phase from teleseismic earthquakes, which is removed during deconvolution in the calculation of receiver functions. As the Moho‐reflectedP ‐wave is independent ofV S, this constraint allows us to create a system of equations that better quantifies the thickness,V P, andV Sof a layer and produces a more appropriate estimation of associated uncertainties. We apply this to 88 seismic stations that are spatially distributed throughout the United States to obtain a model of crustal variability that is unbiased by a priori assumptions of velocity structure.
