The Charlevoix Seismic Zone (CSZ) is located along the early Paleozoic St. Lawrence rift zone in southeastern Quebec at the location of a major Devonian impact structure. The impact structure superimposed major, steeply dipping basement faults trending approximately N35 • E. Approximately 250 earthquakes are recorded each year and are concentrated within and beneath the impact structure. Most M4+ earthquakes associated with the rift faults occurred outside the impact structure. Apart from the unique distribution of earthquakes, stress inversion of focal mechanisms shows stress rotations within the CSZ, and in the CSZ relative to the stress orientation determined from borehole breakouts. The primary goal of this research is to investigate the combined effects of the preexisting structures and regional stresses on earthquake activity and stress rotations in the CSZ. We approach this using PyLith, a finite-element code for simulations of crustal deformation. Adopting the results from recent hypocenter relocation and 3-D tomography studies, we modify the locations and dips of the rift faults and assess the effect of the new fault geometries on stress distributions. We also discuss the effects of resolved velocity anomalies. We find that the observed stress rotation is due to the combined effect of the rift faults and the impact structure. One-dimensional velocity models of the CSZ with an embedded impact structure and a combination of 65 • -40 • -40 • and constant 70 • fault dip models with a very low friction coefficient of 0.3 and cohesion of 0 MPa can explain the observed seismicity and more than 50% of the stress rotations. Key Points:• The CSZ seismicity and stress orientations are the combined effects of the rift faults and the impact structure • Planar faults we considered can explain the observed seismicity but only 50% of the stress rotation • An impact structure 4 times less elastically stiff than the surrounding crust can explain the seismicity in the CSZ Supporting Information:• Supporting Information S1 Figure 2. (a) Model domain with the impact structure (gray) and crust (red). The three rift faults (black lines), dimensions, boundary conditions and the total amount of displacement are annotated. The outer 10-km-thick layer of crust (green) is added to contain fault edges within the domain, a requirement by PyLith. (b) The orientation of the model domain (red dotted box) and the associated coordinate axes (red arrows) relative to the geographic reference frame. The red dotted box is not the actual model domain but just a box that shows its "orientation." Cyan circles on the geographic map show the inner and outer impact structure boundaries and black solid lines trace rift faults (see Figure 1).for the region. We then try to correlate modeled stress distributions with the recently-relocated hypocenters (Powell & Lamontagne, 2017) and discuss implications for the nonuniform stress orientation around the impact structure that was previously recognized.
Accurately modeling time-dependent coseismic crustal deformation as observed on high-rate Global Navigation Satellite System (HR-GNSS) lends insight into earthquake source processes and improves local earthquake and tsunami early warning algorithms. Currently, time-dependent crustal deformation modeling relies most frequently on simplified 1D radially symmetric Earth models. However, for shallow subduction zone earthquakes, even low-frequency shaking is likely affected by the many strongly heterogeneous structures such as the subducting slab, mantle wedge, and the overlying crustal structure. We demonstrate that including 3D structure improves the estimation of key features of coseismic HR-GNSS time series, such as the peak ground displacement (PGD), the time to PGD (tPGD), static displacements (SD), and waveform cross-correlation values. We computed 1D and 3D synthetic, 0.25 Hz and 0.5 Hz waveforms at HR-GNSS stations for four M7.3+ earthquakes in Japan using MudPy and SW4, respectively. From these synthetics, we computed intensity-measure residuals between the synthetic and observed GNSS waveforms. Comparing 1D and 3D residuals, we observed that the 3D simulations show better fits to the PGD and SD in the observed waveforms than the 1D simulations for both 0.25 Hz and 0.5 Hz simulations. We find that the reduction in PGD residuals in the 3D simulations is a combined effect of both shallow and deep 3D structures; hence incorporating only the upper 30 km 3D structure will still improve the fit to the observed PGD values. Our results demonstrate that 3D simulations significantly improve models of GNSS waveform characteristics and will not only help understand the underlying processes, but also improve local tsunami warning.
A refined model of fault structure in the active New Madrid Seismic Zone (NMSZ) is presented based on relocated hypocenters and application of a statistical clustering method to determine fault planes. Over 200 earthquakes are recorded in the NMSZ every year, but the three‐dimensional (3‐D) fault structure is difficult to determine because the zone is covered by thick, Mississippi Embayment sediment. The distribution of earthquakes in the NMSZ indicates four major arms of seismicity, suggesting the presence of a northeast‐southwest trending strike‐slip fault system with a major northwest trending, contractional stepover fault. Relocation of 4,131 earthquakes using HypoDD resulted in major improvement in the depiction of fault structure in the NMSZ including three‐dimensional structural variations along the along the Reelfoot fault (RF) and a well‐defined intersection of the Axial and RFs. Optimal Anisotropic Dynamic Clustering analysis of the relocated hypocenters produced a fault model consisting of 12, well resolved planes. The RF is continuous along strike from the northern end to the Ridgely fault, located south of the intersection with the Axial fault (AF). The crosscutting Ridgely fault may serve as a barrier to rupture propagation along the entire fault. The strike‐slip arms of seismicity are well resolved and correspond to near vertical planes. Three planes are resolved in the southern part of the AF and are associated with the Osceola igneous complex.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.