We use 3D dynamic rupture simulations to discover a previously un‐described effect of asymmetric topography on the earthquake process. With the Cajon Pass along the San Andreas Fault as an example, we find that asymmetric topography generates an alternating normal stress perturbation around the rupture front, near the free surface. When topography lies to the right of the propagating right‐lateral front, the normal stress perturbation is clamping ahead of the rupture front and unclamping behind. When topography alternates to the left, the perturbation reverses sign. The process is analogous to the normal stress variations on dip‐slip faults. While this effect does not strongly affect rupture propagation and slip in our current parametrization, it requires explanation and exploration. An understanding of the normal stress perturbation due to asymmetric topography will help prevent its misattribution to other sources and lead to a better understanding of the interplay of multiple processes during earthquakes.
A Synthesis of Fracture, Friction and Damage Processes in Earthquake Rupture Zones
Abstract We review properties and processes of earthquake rupture zones based on field studies, laboratory observations, theoretical models and simulations, with the goal of assessing the possible dominance of different processes in different parts of the rupture and validity of commonly used models. Rupture zones may be divided into front , intermediate , and tail regions that interact to different extents. The rupture front is dominated by fracturing and granulation processes and strong dilatation, producing faulting products that are reworked by subsequent sliding behind. The intermediate region sustains primarily frictional sliding with relatively high slip rates that produce appreciable stress transfer to the propagating front. The tail region further behind is characterized by low slip rates that effectively do not influence the propagating front, although it (and the intermediate region) can spawn small offspring rupture fronts. Wave-mediated stress transfer can also trigger failures ahead of the rupture front. Earthquake ruptures are often spatially discontinuous and intermittent with a hierarchy of asperity and segment sizes that radiate waves with different tensorial compositions and frequency bands. While different deformation processes dominating parts of the rupture zones can be treated effectively with existing constitutive relations, a more appropriate analysis of earthquake processes would require a model that combines aspects of fracture, damage-breakage, and frictional frameworks.
more »
« less
- Award ID(s):
- 2122168
- NSF-PAR ID:
- 10421629
- Date Published:
- Journal Name:
- Pure and Applied Geophysics
- Volume:
- 179
- Issue:
- 12
- ISSN:
- 0033-4553
- Page Range / eLocation ID:
- 4323 to 4339
- 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
-
Abstract. Substantial insight into earthquake source processes has resulted from considering frictional ruptures analogous to cohesive-zone shear cracks from fracture mechanics. This analogy holds for slip-weakening representations of fault friction that encapsulate the resistance to rupture propagation in the form of breakdown energy, analogous to fracture energy, prescribed in advance as if it were a material property of the fault interface. Here, we use numerical models of earthquake sequences with enhanced weakening due to thermal pressurization of pore fluids to show how accounting for thermo-hydro-mechanical processes during dynamic shear ruptures makes breakdown energy rupture-dependent. We find that local breakdown energy is neither a constant material property nor uniquely defined by the amount of slip attained during rupture, but depends on how that slip is achieved through the history of slip rate and dynamic stress changes during the rupture process. As a consequence, the frictional breakdown energy of the same location along the fault can vary significantly in different earthquake ruptures that pass through. These results suggest the need to reexamine the assumption of predetermined frictional breakdown energy common in dynamic rupture modeling and to better understand the factors that control rupture dynamics in the presence of thermo-hydro-mechanical processes.more » « less
-
Fluids are known to trigger a broad range of slip events, from slow, creeping transients to dynamic earthquake ruptures. Yet, the detailed mechanics underlying these processes and the conditions leading to different rupture behaviors are not well understood. Here, we use a laboratory earthquake setup, capable of injecting pressurized fluids, to compare the rupture behavior for different rates of fluid injection, slow (megapascals per hour) versus fast (megapascals per second). We find that for the fast injection rates, dynamic ruptures are triggered at lower pressure levels and over spatial scales much smaller than the quasistatic theoretical estimates of nucleation sizes, suggesting that such fast injection rates constitute dynamic loading. In contrast, the relatively slow injection rates result in gradual nucleation processes, with the fluid spreading along the interface and causing stress changes consistent with gradually accelerating slow slip. The resulting dynamic ruptures propagating over wetted interfaces exhibit dynamic stress drops almost twice as large as those over the dry interfaces. These results suggest the need to take into account the rate of the pore-pressure increase when considering nucleation processes and motivate further investigation on how friction properties depend on the presence of fluids.more » « less
-
more » « less
Abstract Localized frictional sliding on faults in the continental crust transitions at depth to distributed deformation in viscous shear zones. This brittle‐ductile transition (BDT), and/or the transition from velocity‐weakening (VW) to velocity‐strengthening (VS) friction, are controlled by the lithospheric thermal structure and composition. Here, we investigate these transitions, and their effect on the depth extent of earthquakes, using 2D antiplane shear simulations of a strike‐slip fault with rate‐and‐state friction. The off‐fault material is viscoelastic, with temperature‐dependent dislocation creep. We solve the heat equation for temperature, accounting for frictional and viscous shear heating that creates a thermal anomaly relative to the ambient geotherm which reduces viscosity and facilitates viscous flow. We explore several geotherms and effective normal stress distributions (by changing pore pressure), quantifying the thermal anomaly, seismic and aseismic slip, and the transition from frictional sliding to viscous flow. The thermal anomaly can reach several hundred degrees below the seismogenic zone in models with hydrostatic pressure but is smaller for higher pressure (and these high‐pressure models are most consistent with San Andreas Fault heat flow constraints). Shear heating raises the BDT, sometimes to where it limits rupture depth rather than the frictional VW‐to‐VS transition. Our thermomechanical modeling framework can be used to evaluate lithospheric rheology and thermal models through predictions of earthquake ruptures, postseismic and interseismic crustal deformation, heat flow, and the geological structures that reflect the complex deformation beneath faults.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s00024-022-03168-9
