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: Variations of Crustal Thickness in Alaska and Northwestern Canada: Implications for Crustal Modification and Accretion of Tectonic Units
Abstract The northwestern part of North America has recorded multiple tectonic events, such as terrane accretion, strike‐slip motion, and subduction of the Pacific and Yakutat plates, providing an iconic setting to investigate the tectonic evolution of the continental crust. In this study we analyze the receiver functions at seismic stations deployed during 1999–2022 to estimate the crustal thickness, as well as possible slab signature, in Alaska and northwestern Canada. The Moho signal can be clearly detected within the continental region. Specifically, in northwestern Canada, the thickest crust is observed beneath the Cordilleran Deformation Front, which marks the structural boundary between the North American Craton and the North American Margin. We observe a few distinct offsets in the Moho depth located both within the tectonic units and approximately across the major faults between the tectonic units. We provide a first‐order estimate of the depth gradient of the Moho offsets based on the horizontal distance of the two closest seismic stations across the offsets. We propose that the Moho offsets reflect the cumulative impact of the accretionary orogenies and post‐orogenic tectonic events on crustal modification. The continental Moho signal is weak or obscure in Aleutian and southcentral Alaska, and the oceanic Moho within the subducting plates is likely detected. This study provides new seismic insights into understanding the impacts of the tectonic events on continental formation and evolution.  more » « less
Award ID(s):
1751974
PAR ID:
10467663
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
128
Issue:
10
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. 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
  4. Abstract The chemical composition of the deep continental crust is key to understanding the formation and evolution of the continental crust. Constraining the chemical composition of present‐day deep continental crust is, however, limited by indirect accessibility. This paper presents a modeling method for constraining deep crustal chemical structures from observed crustal seismic structures. We compiled a set of published composition models for the continental crust to construct functional relationships between seismic wave speed and major oxide content in the crust. Phase equilibria and compressional wave speeds (VP) for each composition model were calculated over a range of depths and temperatures of the deep crust. For conditions within the alpha(α)‐quartz stability field, robust functional relationships were obtained betweenVPand major oxide contents of the crust. Based on these relationships, observedVPof the deep crust can be inverted to chemical compositions for regions with given geotherms. We provide a MATLAB code for this process (CalcCrustComp). We apply this method to constrain compositions from deep crustalVPof global typical tectonic settings and the North China Craton (NCC). Our modeling results suggest that the lower crust in subduction‐related and rifting‐related tectonic settings may be more mafic than platforms/shields and orogens. The lowVPsignature in the deep crust of the NCC can be explained by intermediate crustal compositions, higher water contents, and/or higher temperatures. The chemical structure obtained by this method can serve as a reference model to further identify deep crustal features. 
    more » « less
  5. Deep‐crustal magma plumbing at arc volcanoes controls the volume, frequency, and composition of magma being transported to and stored in the upper crust. However, the mid‐to‐lower crust remains a challenging region to image. We explore the mid‐to‐lower crustal velocity structure beneath the Christiana‐Santorini‐Kolumbo Volcanic Field (CSKVF) to better understand how an established stratovolcano and flanking volcano (Santorini and Kolumbo) are fed through the mid‐to‐lower crust. We use active‐source seismic data to obtain a P‐wave velocity model of the crust below the CSKVF. We invert direct and reflected P phases to cover the entire depth extent of the crust and solve for the Moho interface depth. Our model requires a curved Moho interface representative of crustal thickening via underplating. Results show a highVpanomaly in the lower crust under Santorini and a mid‐crustal lowVpanomaly offset from both Santorini and Kolumbo. We find that accumulation of magma is located under the local extensional basin in the upper mid‐crust (<10 km) but is offset at deeper depths. We find evidence for melt storage at 11–13 km depth feeding volcanism at the Kolumbo volcanic chain. This melt is also a plausible source for the 2025 seismic swarm and dike intrusion. Resolution is limited in the mid‐crust below the Santorini caldera, leaving Santorini's mid‐crustal magma plumbing unconstrained. We think it likely that Santorini and Kolumbo have entirely separate crustal plumbing systems and mantle sources, but allow the possibility of a connection in the mid or lower crust. 
    more » « less