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.
Fault‐damage zones comprise multiscale fracture networks that may slip dynamically and interact with the main fault during earthquake rupture. Using 3D dynamic rupture simulations and scale‐dependent fracture energy, we examine dynamic interactions of more than 800 intersecting multiscale fractures surrounding a listric fault, emulating a major listric fault and its damage zone. We investigate 10 distinct orientations of maximum horizontal stress, probing the conditions necessary for sustained slip within the fracture network or activating the main fault. Additionally, we assess the feasibility of nucleating dynamic rupture earthquake cascades from a distant fracture and investigate the sensitivity of fracture network cascading rupture to the effective normal stress level. We model either pure cascades or main fault rupture with limited off‐fault slip. We find that cascading ruptures within the fracture network are dynamically feasible under certain conditions, including: (a) the fracture energy scales with fracture and fault size, (b) favorable relative pre‐stress of fractures within the ambient stress field, and (c) close proximity of fractures. We find that cascading rupture within the fracture network discourages rupture on the main fault. Our simulations suggest that fractures with favorable relative pre‐stress, embedded within a fault damage zone, may lead to cascading earthquake rupture that shadows main fault slip. We find that such off‐fault events may reach moment magnitudes up to
- NSF-PAR ID:
- 10496075
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 3
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
null (Ed.)SUMMARY Earthquake ruptures are generally considered to be cracks that propagate as fracture or frictional slip on pre-existing faults. Crack models have been used to describe the spatial distribution of fault offset and the associated static stress changes along a fault, and have implications for friction evolution and the underlying physics of rupture processes. However, field measurements that could help refine idealized crack models are rare. Here, we describe large-scale laboratory earthquake experiments, where all rupture processes were contained within a 3-m long saw-cut granite fault, and we propose an analytical crack model that fits our measurements. Similar to natural earthquakes, laboratory measurements show coseismic slip that gradually tapers near the rupture tips. Measured stress changes show roughly constant stress drop in the centre of the ruptured region, a maximum stress increase near the rupture tips and a smooth transition in between, in a region we describe as the earthquake arrest zone. The proposed model generalizes the widely used elliptical crack model by adding gradually tapered slip at the ends of the rupture. Different from the cohesive zone described by fracture mechanics, we propose that the transition in stress changes and the corresponding linear taper observed in the earthquake arrest zone are the result of rupture termination conditions primarily controlled by the initial stress distribution. It is the heterogeneous initial stress distribution that controls the arrest of laboratory earthquakes, and the features of static stress changes. We also performed dynamic rupture simulations that confirm how arrest conditions can affect slip taper and static stress changes. If applicable to larger natural earthquakes, this distinction between an earthquake arrest zone (that depends on stress conditions) and a cohesive zone (that depends primarily on strength evolution) has important implications for how seismic observations of earthquake fracture energy should be interpreted.more » « less
-
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
-
more » « less
SUMMARY Earthquakes come in clusters formed of mostly aftershock sequences, swarms and occasional foreshock sequences. This clustering is thought to result either from stress transfer among faults, a process referred to as cascading, or from transient loading by aseismic slip (pre-slip, afterslip or slow slip events). The ETAS statistical model is often used to quantify the fraction of clustering due to stress transfer and to assess the eventual need for aseismic slip to explain foreshocks or swarms. Another popular model of clustering relies on the earthquake nucleation model derived from experimental rate-and-state friction. According to this model, earthquakes cluster because they are time-advanced by the stress change imparted by the mainshock. This model ignores stress interactions among aftershocks and cannot explain foreshocks or swarms in the absence of transient loading. Here, we analyse foreshock, swarm and aftershock sequences resulting from cascades in a Discrete Fault Network model governed by rate-and-state friction. We show that the model produces realistic swarms, foreshocks and aftershocks. The Omori law, characterizing the temporal decay of aftershocks, emerges in all simulations independently of the assumed initial condition. In our simulations, the Omori law results from the earthquake nucleation process due to rate and state friction and from the heterogeneous stress changes due to the coseismic stress transfers. By contrast, the inverse Omori law, which characterizes the accelerating rate of foreshocks, emerges only in the simulations with a dense enough fault system. A high-density complex fault zone favours fault interactions and the emergence of an accelerating sequence of foreshocks. Seismicity catalogues generated with our discrete fault network model can generally be fitted with the ETAS model but with some material differences. In the discrete fault network simulations, fault interactions are weaker in aftershock sequences because they occur in a broader zone of lower fault density and because of the depletion of critically stressed faults. The productivity of the cascading process is, therefore, significantly higher in foreshocks than in aftershocks if fault zone complexity is high. This effect is not captured by the ETAS model of fault interactions. It follows that a foreshock acceleration stronger than expected from ETAS statistics does not necessarily require aseismic slip preceding the mainshock (pre-slip). It can be a manifestation of a cascading process enhanced by the topological properties of the fault network. Similarly, earthquake swarms might not always imply transient loading by aseismic slip, as they can emerge from stress interactions.
-
more » « less
Abstract We use high‐resolution earthquake locations to characterize the three‐dimensional structure of active faults in California and how it evolves with fault structural maturity. We investigate the distribution of aftershocks of several recent large earthquakes that occurred on continental strike slip faults of various structural maturity (i.e. various cumulative fault displacement, length, initiation age and slip rate). Aftershocks define a tabular zone of shear deformation surrounding the mainshock rupture plane. Comparing this to geological observations, we conclude that this results from the re‐activation of secondary faults. We observe a rapid fall off of the number of aftershocks at a distance range of 0.06‐0.22 km from the main fault surface of mature faults, and 0.6‐1.0 km from the fault surface of immature faults. The total width of the active shear deformation zone surrounding the main fault plane reaches 1.0‐2.5 km and 6‐9 km for mature and immature faults, respectively. We find that the width of the shear deformation zone decreases as a power law with cumulative fault displacement. Comparing with a dynamic rough fault model, we infer that the narrowing of the shear deformation zone agrees quantitatively with earlier estimates of the smoothing of faults with displacement, both of which are aspects of fault wear. We find that earthquake stress drop decreases with fault displacement and hence with increased smoothness and/or slip rate. This may result from fault healing or the effect of roughness on friction.
