Mapping finite-fault slip with spatial correlation between seismicity and point-source Coulomb failure stress change

Most earthquake energy release arises during fault slip many kilometers below the Earth's surface. Understanding earthquakes and their hazard requires mapping the geometry and amount of this slip. Such finite-fault maps are typically derived from measures of surface phenomena, such as seismic and geodetic ground motions. Here we introduce a procedure for mapping finite-fault slip directly from seismicity and aftershocks-phenomena occurring at depth around an earthquake rupture. For specified source and receiver faults, we map source-fault slip in 3D by correlation of point-source Coulomb failure stress change ({Δ}CFS) kernels across the distribution of seismicity around an earthquake. These seismicity-stress maps show relative, static fault slip compatible with the surrounding seismicity under the physics of {Δ}CFS; they can aid other slip inversions and aftershock forecasting, and be obtained for early instrumental earthquakes. We verify the maps recover synthetic fault slip and match independent estimates of slip for the 2004 Mw 6.0 Parkfield and 2021 Mw 6.0 Antelope Valley California earthquakes. For the 2018 Mw 7.1 Anchorage Alaska intra-slab earthquake, seismicity-stress maps, combined with multi-scale precise hypocenter relocation, resolve the enigma of which mainshock faulting plane ruptured (the gently east-dipping plane), and clarify slab structures activated in the energetic aftershock sequence.