Comparison of average stress drop measures for ruptures with heterogeneous stress change and implications for earthquake physics

Abstract

Stress drop, a measure of static stress change in earthquakes, is the subject of numerous investigations. Stress drop in an earthquake is likely to be spatially varying over the fault, creating a stress drop distribution. Representing this spatial distribution by a single number, as commonly done, implies averaging in space. In this study, we investigate similarities and differences between three different averages of the stress drop distribution used in earthquake studies. The first one, Δσ¯¯¯¯¯M, is the commonly estimated stress drop based on the seismic moment and fault geometry/dimensions. It is known that Δσ¯¯¯¯¯M corresponds to averaging the stress drop distribution with the slip distribution due to uniform stress drop as the weighting function. The second one, Δσ¯¯¯¯¯A, is the simplest (unweighted) average of the stress drop distribution over the fault, equal to the difference between the average stress levels on the fault before and after an earthquake. The third one, Δσ¯¯¯¯¯E, enters discussions of energy partitioning and radiation efficiency; we show that it corresponds to averaging the stress drop distribution with the actual final slip at each point as the weighting function. The three averages, Δσ¯¯¯¯¯M, Δσ¯¯¯¯¯A, and Δσ¯¯¯¯¯E, are often used interchangeably in earthquake studies and simply called ‘stress drop’. Yet they are equal to each other only for ruptures with spatially uniform stress drop, which results in an elliptical slip distribution for a circular rupture. Indeed, we find that other relatively simple slip shapes—such as triangular, trapezoidal or sinusoidal—already result in stress drop distributions with notable differences between Δσ¯¯¯¯¯M, Δσ¯¯¯¯¯A, and Δσ¯¯¯¯¯E. Introduction of spatial slip heterogeneity results in further systematic differences between them, with Δσ¯¯¯¯¯E always being larger than Δσ¯¯¯¯¯M, a fact that we have proven theoretically, and Δσ¯¯¯¯¯A almost always being the smallest. In particular, the value of the energy-related Δσ¯¯¯¯¯E significantly increases in comparison to the moment-based Δσ¯¯¯¯¯M with increasing roughness of the slip distribution over the fault. Previous studies used Δσ¯¯¯¯¯M in place of Δσ¯¯¯¯¯E in computing the radiation ratio ηR that compares the radiated energy in earthquakes to a characteristic part of their strain energy change. Typical values of ηR for large earthquakes were found to be from 0.25 to 1. Our finding that Δσ¯¯¯¯¯E≥Δσ¯¯¯¯¯M allows us to interpret the values of ηR as the upper bound. We determine the restrictions placed by such estimates on the evolution of stress with slip at the earthquake source. We also find that Δσ¯¯¯¯¯E can be approximated by Δσ¯¯¯¯¯M if the latter is computed based on a reduced rupture area.