Mature strike‐slip faults are usually surrounded by a narrow zone of damaged rocks characterized by low seismic wave velocities. Observations of earthquakes along such faults indicate that seismicity is highly concentrated within this fault damage zone. However, the long‐term influence of the fault damage zone on complete earthquake cycles, that is, years to centuries, is not well understood. We simulate aseismic slip and dynamic earthquake rupture on a vertical strike‐slip fault surrounded by a fault damage zone for a thousand‐year timescale using fault zone material properties and geometries motivated by observations along major strike‐slip faults. The fault damage zone is approximated asan elastic layer with lower shear wave velocity than the surrounding rock. We find that dynamic wave reflections, whose characteristics are strongly dependent on the width and the rigidity contrast of the fault damage zone, have a prominent effect on the stressing history of the fault. The presence of elastic damage can partially explain the variability in the earthquake sizes and hypocenter locations along a single fault, which vary with fault damage zone depth, width and rigidity contrast from the host rock. The depth extent of the fault damage zone has a pronounced effect on the earthquake hypocenter locations, and shallower fault damage zones favor shallower hypocenters with a bimodal distribution of seismicity along depth. Our findings also suggest significant effects on the hypocenter distribution when the fault damage zone penetrates to the nucleation sites of earthquakes, likely being influenced by both lithological (material) and rheological (frictional) boundaries.
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- PAR ID:
- 10371671
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 127
- Issue:
- 7
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
null (Ed.)Earthquake prediction is the holy grail of seismology. Many previous studies have searched for robust precursory signals to inform us of imminent earthquakes, the most significant of which are seen in laboratory experiments as temporal changes in pressure and shear wave velocities during the seismic cycle. Similar changes are seen in natural faults and the surrounding structurally complex network of fractures with nested hierarchy of localized deformation, referred to as fault damage zone. However, little is known whether such temporal changes in material properties contains any precursory signals for imminent earthquakes.Conversely, the effect of precursory velocity changes on the seismic cycle is not well understood. By imposing shear wave velocity changes in fault damage zones, we investigate the effects of these precursors on multiple stages of the seismic cycle, including nucleation, coseismic, postseismic, and interseismic stages. We perform 2D fully dynamic earthquake cycle simulations with a fault-parallel damage zone for strike-slip fault systems with antiplane geometry. The fault is governed by rate-state-dependent friction laws, and the fault damage zone material is considered elastic. Our preliminary results show that the temporal onset of shear wave velocity drop causes a reduction in earthquake recurrence intervals over the seismic cycle. Furthermore, a dynamic earthquake rupture within the seismic cycle terminates much faster and abruptly in models with precursory velocity changes. We will also discuss how the precursory velocity changes affect the fault-slip behavior, including fast-slip, slow-slip, and aseismic creep, for different amplitudes of shear wave velocity changes at different compliance contrast of the fault damage zones. Our results highlight the importance of short and long-term monitoring of fault zone structures for better assessment of regional seismic hazard.more » « less
-
"Important physical observations in rupture dynamics such as static fault friction, short-slip, self-healing, and supershear phenomenon in cracks are studied. A continuum model of rupture dynamics is developed using the field dislocation mechanics (FDM) theory. The energy density function in our model encodes accepted and simple physical facts related to rocks and granular materials under compression. We work within a 2-dimensional ansatz of FDM where the rupture front is allowed to move only in a horizontal fault layer sandwiched between elastic blocks. Damage via the degradation of elastic modulus is allowed to occur only in the fault layer, characterized by the amount of plastic slip. The theory dictates the evolution equation of the plastic shear strain to be a Hamilton-Jacobi (H-J) equation, resulting in the representation of a propagating rupture front. A Central-Upwind scheme is used to solve the H-J equation. The rupture propagation is fully coupled to elastodynamics in the whole domain, and our simulations recover static friction laws as emergent features of our continuum model, without putting in by hand any such discontinuous criteria in our model. Estimates of material parameters of cohesion and friction angle are deduced. Short-slip and slip-weakening (crack-like) behaviors are also reproduced as a function of the degree of damage behind the rupture front. The long-time behavior of a moving rupture front is probed, and it is deduced that the equilibrium profiles under no shear stress are not traveling wave profiles under non-zero shear load in our model. However, it is shown that a traveling wave structure is likely attained in the limit of long times. Finally, a crack-like damage front is driven by an initial impact loading, and it is observed in our numerical simulations that an upper bound to the crack speed is the dilatational wave speed of the material unless the material is put under pre-stressed conditions, in which case supersonic motion can be obtained. Without pre-stress, intersonic (supershear) motion is recovered under appropriate conditions."more » « less
-
null (Ed.)The fault damage zone is a well-known structure of localized deformation around faults. Its material properties evolve over earthquake cycles due to coseismic damage accumulation and interseismic healing. We will present fully dynamic earthquake cycle simulations to show how the styles of earthquake nucleation and rupture propagation change as fault zone material properties vary temporally. First, we will focus on the influence of fault zone structural maturity quantified by near-fault seismic wave velocities in simulations. The simulations show that immature fault zones promote small and moderate subsurface earthquakes with irregular recurrence intervals, whereas mature fault zones host pulse-like earthquake rupture that can propagate to the surface, extend throughout the seismogenic zone, and occur at regular intervals. The interseismic healing in immature fault zones plays a key role in allowing the development of aseismic slip episodes including slow-slip events and creep, which can propagate into the seismogenic zone, and thus limit the sizes of subsequent earthquakes by releasing fault stress. In the second part, we will discuss how the precursory changes of seismic wave velocities of fault damage zones may affect earthquake nucleation process. Both laboratory experiments and seismic observations show that the abrupt earthquake failure can be preceded by accelerated fault deformation and the accompanying velocity reduction of near-fault rocks. We will use earthquake cycle simulations to systematically test the effects of timing and amplitudes of such precursory velocity changes. Our simulations will provide new insights into the interplay between fault zone structure and earthquake nucleation process, which can be used to guide future real-time monitoring of major fault zones.more » « less
-
more » « less
Abstract A devastating magnitude 7.2 earthquake struck Southern Haiti on 14 August 2021. The earthquake caused severe damage and over 2000 casualties. Resolving the earthquake rupture process can provide critical insights into hazard mitigation. Here we use integrated seismological analyses to obtain the rupture history of the 2021 earthquake. We find the earthquake first broke a blind thrust fault and then jumped to a disconnected strike‐slip fault. Neither of the fault configurations aligns with the left‐lateral tectonic boundary between the Caribbean and North American plates. The complex multi‐fault rupture may result from the oblique plate convergence in the region, so that the initial thrust rupture is due to the boundary‐normal compression and the following strike‐slip faulting originates from the Gonâve microplate block movement, orienting SW‐NE direction. The complex rupture development of the earthquake suggests that the regional deformation is accommodated by a network of segmented faults with diverse faulting conditions.
