Weng H. H., Ampuero Jean-Paul. (2019). The dynamics of elongated earthquake ruptures. Journal of Geophysical Research : Solid Earth, 124 (8), p. 8584-8610. ISSN 2169-9313.
Source
Journal of Geophysical Research : Solid Earth, 2019,
124 (8), p. 8584-8610 ISSN 2169-9313
The largest earthquakes propagate laterally after saturating the fault's seismogenic width and reach large length-to-width ratios L/W. Smaller earthquakes can also develop elongated ruptures due to confinement by heterogeneities of initial stresses or material properties. The energetics of such elongated ruptures is radically different from that of conventional circular crack models: they feature width-limited rather than length-dependent energy release rate. However, a synoptic understanding of their dynamics is still missing. Here we combine computational and analytical modeling of long ruptures in three dimension (3D) and 2.5D (width-averaged) to develop a theoretical relation between the evolution of rupture speed and the along-strike distribution of fault stress, fracture energy, and rupture width. We find that the evolution of elongated ruptures in our simulations is well described by the following rupture-tip-equation-of-motion: G(c) = G(0) (1 - (v)over dot(r)W/v(s)(2) gamma/A alpha(P)(s)) where G(c) is the fracture energy, G(0) is the steady state energy release rate, v(s) is the S wave speed, v(r) is the rupture speed, (v)overdot(r) = dv(r)/dt is the rupture acceleration, and gamma/A alpha(P)(s) is a known function of rupture speed. The steady energy release rate is limited by rupture width as G(0) = gamma Delta tau W-2/mu, where gamma is a geometric factor, Delta tau is the stress drop (spatially smoothed over a length scale smaller than W), and mu is the shear modulus. If G(c) is a constant and exactly balanced by G(0), the rupture can in principle propagate steadily at any speed. If G(c) increases with rupture speed, steady ruptures have a well-defined speed and are stable. When G(c) not equal G(0), the rupture acquires an inertial effect: the rupture-tip-equation-of-motion depends explicitly on rupture acceleration. This inertial effect does not exist in the classical theory of dynamic rupture in 2-D unbounded media and in unbounded faults in 3D, but emerges in 2-D bounded media or, as shown here, as a consequence of the finite rupture width in 3D. These findings highlight the essential role of the seismogenic width on rupture dynamics. Based on the rupture-tip-equation-of-motion we define the rupture potential, a function that determines the size of next earthquake, and we propose a conceptual model that helps rationalize one type of "supercycles" observed on segmented faults. More generally, the theory developed here can yield relations between earthquake source properties (final magnitude, moment rate function, radiated energy) and the heterogeneities of stress and strength along the fault, which can then be used to extract statistical information on fault heterogeneity from source time functions of past earthquakes or as physics-based constraints on finite-fault source inversion and on seismic hazard assessment.
