Effects of Fault Roughness on Coseismic Slip and Earthquake Locations

Fault zone structure is well known to exert strong controls on earthquake properties including coseismic slip distribution, rupture propagation direction, and hypocenter location. It has also been well established that the principal slip surface, which accommodates the majority of earthquake displacement, exhibits roughness at all scales following self-affine fractal distributions. Here we explore the relationship between fault roughness and specific earthquake properties including coseismic slip distribution and hypocenter location based on long-term simulations of earthquake catalogs on fractally rough faults. We begin by using the von Karman autocorrelation function to procedurally generate single faults with different fractal roughness properties, which we place in a homogeneous elastic solid and apply pure right-lateral shear at a constant back slip rate with the earthquake simulator RSQSim. Running the simulations for 10,000 years each, we generate millions of earthquakes including thousands of events with Mw > 6.0, which rupture the surface. We show that the patterns of surface rupture in these large events follow self-affine fractal distributions with consistent fractal dimension related distinctly to the fractal dimension of the fault. In addition, the hypocenters of these large events occur in very specific predictable locations where the longest wavelength structure produces a stress asperity (i.e., restraining bend). The resulting patterns can explain many features observed on real fault systems, including clustered hypocenter locations, spatially variable coseismic slip distributions, and characteristic slip recurrence behavior. These results demonstrate a quantitative link between a directly measurable fault property-roughness-and the properties of future earthquakes.