3D finite-difference dynamic-rupture modeling along nonplanar faults

Proper understanding of seismic emissions associated with the growth of complexly shaped microearthquake networks and larger-scale nonplanar fault ruptures, both in arbitrarily heterogeneous media, requires accurate modeling of the underlying dynamic processes. We present a new 3D dynamic-rupture, finite-difference model called the finite-difference, fault-element (FDFE) method; it simulates the dynamic rupture of nonplanar faults subjected to regional loads in complex media. FDFE is based on a 3D methodology for applying dynamic-rupture boundary conditions along the fault surface. The fault is discretized by a set of parallelepiped fault elements in which specific boundary conditions are applied. These conditions are applied to the stress tensor, once transformed into a local fault reference frame. Numerically determined weight functions multiplying particle velocities around each element allow accurate estimates of fault kinematic parameters (i.e., slip and slip rate) independent of faulting mechanism. Assuming a Coulomb-like slip-weakening friction law, a parametric study suggests that the FDFE method converges toward a unique solution, provided that the cohesive zone behind the rupture front is well resolved (i.e., four or more elements inside this zone). Solutions are free of relevant numerical artifacts for grid sizes smaller than approximately 70 In. Results yielded by the FDFE approach are in good quantitative agreement with those obtained by a semianalytical boundary integral method along planar and nonplanar parabola-shaped faults. The FDFE method thus provides quantitative, accurate results for spontaneous-rupture simulations on intricate fault geometries.