skip to main content

Title: Multicycle Dynamics of the Aksay Bend Along the Altyn Tagh Fault in Northwest China: 2. The Realistically Complex Fault Geometry

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
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Page Range / eLocation ID:
p. 1120-1137
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    I perform multicycle dynamic simulations of a restraining double bend that is simplified from the Aksay bend along the Altyn Tagh Fault in northwest China. The earthquake cycle includes a coseismic dynamic rupture phase and an interseismic stress loading and relaxation phase. The double bend fault system has two fault strands and each strand comprises a stem segment that strikes parallel to the maximum shear loading direction and a bend segment that deviates from the loading direction. I find that the restraining double bend is an effective barrier to dynamically propagating ruptures over multiple earthquake cycles, with a probability of rupture jumping from one strand to the other varying from 0% to 10%, depending on parameter values of the rate dependence of friction during the dynamic phase and the viscosity for stress relaxation during the interseismic phase. The critical condition for the primary rupture on the first strand to jump onto the second strand is that the preshear stress is close enough to the shear strength over a certain size of the fault patch on the second strand. These results on the simplified double bend provide important insights into rupture behavior of geometrically complex faults over multiple earthquake cycles in general and a reference to explore rupture behavior of the Aksay bend in particular.

    more » « less
  2. 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

    We develop a finite element dynamic earthquake simulator, EQsimu, to model multicycle dynamics of 3-D geometrically complex faults. The fault is governed by rate- and state-dependent friction (RSF). EQsimu integrates an existing finite element code EQdyna for the coseismic dynamic rupture phase and a newly developed finite element code EQquasi for the quasi-static phases of an earthquake cycle, including nucleation, post-seismic and interseismic processes. Both finite element codes are parallelized through Message Passing Interface to improve computational efficiency and capability. EQdyna and EQquasi are coupled through on-fault physical quantities of shear and normal stresses, slip-rates and state variables in RSF. The two-code scheme shows advantages in reconciling the computational challenges from different phases of an earthquake cycle, which include (1) handling time-steps ranging from hundredths of a second to a fraction of a year based on a variable time-stepping scheme, (2) using element size small enough to resolve the cohesive zone at rupture fronts of dynamic ruptures and (3) solving the system of equations built up by millions of hexahedral elements.

    EQsimu is used to model multicycle dynamics of a 3-D strike-slip fault with a bend. Complex earthquake event patterns spontaneously emerge in the simulation, and the fault demonstrates two phases in its evolution. In the first phase, there are three types of dynamic ruptures: ruptures breaking the whole fault from left to right, ruptures being halted by the bend, and ruptures breaking the whole fault from right to left. As the fault bend experiences more ruptures, the zone of stress heterogeneity near the bend widens and the earthquake sequence enters the second phase showing only repeated ruptures that break the whole fault from left to right. The two-phase behaviours of this bent fault system suggest that a 10° bend may conditionally stop dynamic ruptures at the early stage of a fault system evolution and will eventually not be able to stop ruptures as the fault system evolves. The nucleation patches are close to the velocity strengthening region. Their sizes on the two fault segments are different due to different levels of the normal stress.

    more » « less
  4. Abstract

    Understanding mechanical conditions that lead to complexity in earthquakes is important to seismic hazard analysis. In this study, we simulate physics‐based multicycle dynamic models of the San Andreas fault (Carrizo through San Bernardino sections) and the San Jacinto fault (Claremont and Clark strands). We focus on a complex fault geometry based on the Southern California Earthquake Center Community Fault Model and its effect over multiple earthquake cycles. Using geodetically derived strain rates, we validate the models against geologic slip rates and recurrence intervals at various paleoseismic sites. We find that the interactions among fault geometry, dynamic rupture and interseismic stress accumulation produce stress heterogeneities, leading to rupture segmentation and variability in earthquake recurrence. Our models produce earthquakes with rupture extents similar to a recent comprehensive paleoseismic catalog. The “earthquake gates” of the Big Bend and the Cajon Pass occasionally impede dynamic ruptures. The angle of compression, which is the subtraction of the maximum shear strain rate direction from the local fault strike, can better determine the likelihood of the impedance of restraining bends to dynamic ruptures. Because the Big Bend has an angle of compression of ∼20°, ruptures that traverse the Big Bend, like the 1857 Fort Tejon earthquake, are more frequent than expected based on empirical relations which predict the ∼40° restraining bend to terminate most ruptures. Our models indicate that large ruptures tend to initiate north of the Big Bend and propagate southwards, similar to the 1857 earthquake, providing critical information for ground shaking assessment in the region.

    more » « less
  5. Abstract

    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