skip to main content

Title: Dynamic Rupture Models, Fault Interaction and Ground Motion Simulations for the Segmented Húsavík‐Flatey Fault Zone, Northern Iceland

The Húsavík‐Flatey Fault Zone (HFFZ) is the largest strike‐slip fault in Iceland and poses a high seismic risk to coastal communities. To investigate physics‐based constraints on earthquake hazards, we construct three fault system models of varying geometric complexity and model 79 3‐D multi‐fault dynamic rupture scenarios in the HFFZ. By assuming a simple regional prestress and varying hypocenter locations, we analyze the rupture dynamics, fault interactions, and the associated ground motions up to 2.5 Hz. All models account for regional seismotectonics, topo‐bathymetry, 3‐D subsurface velocity, viscoelastic attenuation, and off‐fault plasticity, and we explore the effect of fault roughness. The rupture scenarios obey earthquake scaling relations and predict magnitudes comparable to those of historical events. We show how fault system geometry and segmentation, hypocenter location, and prestress can affect the potential for rupture cascading, leading to varying slip distributions across different portions of the fault system. Our earthquake scenarios yield spatially heterogeneous near‐field ground motions modulated by geometric complexities, topography, and rupture directivity, particularly in the near‐field. The average ground motion attenuation characteristics of dynamic rupture scenarios of comparable magnitudes and mean stress drop are independent of variations in source complexity, magnitude‐consistent and in good agreement with the latest regional empirical ground motion models. However, physics‐based ground motion variability changes considerably with fault‐distance and increases for unilateral compared to bilateral ruptures. Systematic variations in physics‐based near‐fault ground motions provide important insights into the mechanics and potential earthquake hazard of large strike‐slip fault systems, such as the HFFZ.

more » « less
Award ID(s):
2121666 2121568
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Mature faults with large cumulative slip often separate rocks with dissimilar elastic properties and show asymmetric damage distribution. Elastic contrast across such bimaterial faults can significantly modify various aspects of earthquake rupture dynamics, including normal stress variations, rupture propagation direction, distribution of ground motions, and evolution of off‐fault damage. Thus, analyzing elastic contrasts of bimaterial faults is important for understanding earthquake physics and related hazard potential. The effect of elastic contrast between isotropic materials on rupture dynamics is relatively well studied. However, most fault rocks are elastically anisotropic, and little is known about how the anisotropy affects rupture dynamics. We examine microstructures of the Sandhill Corner shear zone, which separates quartzofeldspathic rock and micaceous schist with wider and narrower damage zones, respectively. This shear zone is part of the Norumbega fault system, a Paleozoic, large‐displacement, seismogenic, strike‐slip fault system exhumed from middle crustal depths. We calculate elastic properties and seismic wave speeds of elastically anisotropic rocks from each unit having different proportions of mica grains aligned sub‐parallel to the fault. Our findings show that the horizontally polarized shear wave propagating parallel to the bimaterial fault (with fault‐normal particle motion) is the slowest owing to the fault‐normal compliance and therefore may be important in determining the elastic contrast that affects rupture dynamics in anisotropic media. Following results from subshear rupture propagation models in isotropic media, our results are consistent with ruptures preferentially propagated in the slip direction of the schist, which has the slower horizontal shear wave and larger fault‐normal compliance.

    more » « less
  2. Abstract

    Spontaneous rupture simulations along geometrically rough faults have been shown to produce realistic far‐field spectra and comparable fits with ground motion metrics such as spectral accelerations and peak motions from Ground Motion Prediction Equations (GMPEs), but they are too computationally demanding for use with physics‐based probabilistic seismic hazard analysis efforts. Here, we present our implementation of a kinematic rupture generator that matches the characteristics of, at least in a statistical sense, rough‐fault spontaneous rupture models. To this end, we analyze ~100 dynamic rupture simulations on strike‐slip faults withMwranging from 6.4 to 7.2. We find that our dynamic simulations follow empirical scaling relationships for strike‐slip events and provide source spectra comparable to a source model withω−2decay. To define our kinematic source model, we use a regularized Yoffe function parameterized in terms of slip, peak‐time, rise‐time, and rupture initiation time. These parameters are defined through empirical relationships with random fields whose one‐ and two‐point statistics are derived from the dynamic rupture simulations. Our rupture generator reproduces Next Generation Attenuation (NGA) West2 GMPE medians and intraevent standard deviations of spectral accelerations with periods as short as 0.2 s for ensembles of ground motion simulations. Our rupture generator produces kinematic source models forM6.4–7.2 strike‐slip scenarios that can be used in broadband physics‐based probabilistic seismic hazard efforts or to supplement data in areas of limited observations for the development of future GMPEs.

    more » « less
  3. ABSTRACT We explore the response of ground motions to topography during large crustal fault earthquakes by simulating several magnitude 6.5–7.0 rupture scenarios on the Seattle fault, Washington State. Kinematic simulations are run using a 3D spectral element code and a detailed seismic velocity model for the Puget Sound region. This model includes realistic surface topography and a near-surface low-velocity layer; a mesh spacing of ∼30 m at the surface allows modeling of ground motions up to 3 Hz. We simulate 20 earthquake scenarios using different slip distributions and hypocenter locations on a planar fault surface. Results indicate that average ground motions in simulations with and without topography are similar. However, shaking amplification is common at topographic highs, and more than a quarter of all sites experience short-period (≤2 s) ground-motion amplification greater than 25%–35%, compared with models without topography. Comparisons of peak ground velocity at the top and bottom of topographic features demonstrate that amplification is sensitive to period, with the greatest amplifications typically manifesting near a topographic feature’s estimated resonance frequency and along azimuths perpendicular to its primary axis of elongation. However, interevent variability in topographic response can be significant, particularly at shorter periods (<1 s). We do not observe a clear relationship between source centroid-to-site azimuths and the strength of topographic amplification. Overall, our results suggest that although topographic resonance does influence the average ground motions, other processes (e.g., localized focusing and scattering) also play a significant role in determining topographic response. However, the amount of consistent, significant amplification due to topography suggests that topographic effects should likely be considered in some capacity during seismic hazard studies. 
    more » « less
  4. ABSTRACT The July 2019 Ridgecrest, California, earthquake sequence involved two large events—the M 6.4 foreshock and the M 7.1 mainshock that ruptured a system of intersecting strike-slip faults. We present analysis of space geodetic observations including Synthetic Aperture Radar and Global Navigation Satellite System data, geological field mapping, and seismicity to constrain the subsurface rupture geometry and slip distribution. The data render a complex pattern of faulting with a number of subparallel as well as cross-cutting fault strands that exhibit variations in both strike and dip angles, including a “flower structure” formed by shallow splay faults. Slip inversions are performed using both homogeneous and layered elastic half-space models informed by the local seismic tomography data. The inferred slip distribution suggests a moderate amount of the shallow coseismic slip deficit. The peak moment release occurred in the depth interval of 3–4 km, consistent with results from previous studies of major strike-slip earthquakes, and the depth distribution of seismicity in California. We use the derived slip models to investigate stress transfer and possible triggering relationships between the M 7.1 mainshock and the M 6.4 foreshock, as well as other moderate events that occurred in the vicinity of the M 7.1 hypocenter. Triggering is discouraged for the average strike of the M 7.1 rupture (320°) but encouraged for the initial orientation of the mainshock rupture suggested by the first-motion data (340°). This lends support to a scenario according to which the earthquake rupture nucleated on a small fault that was more optimally oriented with respect to the regional stress and subsequently propagated along the less-favorably oriented pre-existing faults, possibly facilitated by dynamic weakening. The nucleation site of the mainshock experienced positive dynamic Coulomb stress changes that are much larger than the static stress changes, yet the former failed to initiate rupture. 
    more » « less
  5. Abstract

    We study rupture behavior of the Aksay bend along the Altyn Tagh fault in northwest China over multiple earthquake cycles. A finite element method is used to numerically simulate spontaneous rupture during a coseismic process, and a viscoelastic model is used to analytically compute the fault stresses during an interseismic process. We find that the Aksay bend is an effective barrier to halt dynamically propagating ruptures from either side of the bend within a range of model parameters, with statistically only about 10% of ruptures jumping across the bend and propagating through almost the entire local fault system. Secondary complexities in fault geometry within the bend, in particular those portions that align relatively well with the regional strike of the fault system, play a critical role in these occasionally jumping ruptures. Well‐developed fault patches with shear stress close to shear strength allow dynamically propagating ruptures to penetrate into the bend and are more susceptible to the dynamic triggering that enables rupture to jump across the bend onto the other strand. We identify additionally nine large rupture scenarios with different occurrences, and most of them rupture one strand outside the bend with triggered slip on some portions of the same or the other strand within the bend. Slip rate distributions from the models show significantly reduced fault slip within the bend and a permanently locked portion on the south strand near the peak of the Altun Mountains. These findings have important implications for seismic hazard assessments of complex fault systems worldwide.

    more » « less