An official website of the United States government
Abstract The impact of past tectonic events on the formation and modification of continental lithosphere remains as an open question of fundamental importance. Eastern North America provides a complete record of supercontinent assembly and breakup over the past 1.3 Ga, serving as a natural laboratory for our understanding of continental crust and mantle lithosphere and for integrating geologic and geophysical observations. In this study, we used teleseismic Ps receiver functions to image the detailed distribution of crustal thickness beneath eastern North America. The radial‐component receiver functions were calculated from seismic waveforms recorded by a total of 659 broadband stations during 2010–2017, yielding a high‐resolution image of Moho depth distribution. The depths of the Moho and intracrustal layers vary within and across the major tectonic units. Specifically, there are distinct differences in crustal thickness between the northern and southern Grenville Province. A dipping intracrustal feature can be seen within the central Grenville Province, with the depth increasing eastward from 5 to 27 km. The Moho depth decreases southeastward across the Grenville‐Appalachian boundary, with a sharp Moho offset of up to 12–15 km in the central segment and a more gradual variation to the north and south. The thickness difference between the southern and northern Grenville‐aged crusts suggests different tectonic and/or exhumation histories during and after the Grenville Orogeny. The low‐angle eastward dipping crustal feature is interpreted to be a Grenville‐aged collisional structure. Differences in the steepness of the Moho offset along the strike of Appalachians probably reflect variation of the steepness of the subsurface boundary between Laurentia and accreted terranes with different intensities of postorogenic modification. The observed spatial relation between the geologically defined tectonic boundaries and crustal thickness variations provides new constraints on the depth extent of the tectonic units within the crust.
more »
« less
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
- 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
-
-
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
-
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
-
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 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
