Afterslip Moment Scaling and Variability From a Global Compilation of Estimates

Aseismic afterslip is postseismic fault sliding that may significantly redistribute crustal stresses and drive aftershock sequences. Afterslip is typically modeled through geodetic observations of surface deformation on a case-by-case basis, thus questions of how and why the afterslip moment varies between earthquakes remain largely unaddressed. We compile 148 afterslip studies following 53 M(w)6.0-9.1 earthquakes, and formally analyze a subset of 88 well-constrained kinematic models. Afterslip and coseismic moments scale near-linearly, with a median Spearman's rank correlation coefficient (CC) of 0.91 after bootstrapping (95% range: 0.89-0.93). We infer that afterslip area and average slip scale with coseismic moment as Mo2/3 ${M}_{o}<^>{2/3}$ and Mo1/3 ${M}_{o}<^>{1/3}$, respectively. The ratio of afterslip to coseismic moment (M-rel) varies from 300% (interquartile range: 9%-32%). M-rel weakly correlates with M-o (CC: -0.21, attributed to a publication bias), rupture aspect ratio (CC: -0.31), and fault slip rate (CC: 0.26, treated as a proxy for fault maturity), indicating that these factors affect afterslip. M-rel does not correlate with mainshock dip, rake, or depth. Given the power-law decay of afterslip, we expected studies that started earlier and spanned longer timescales to capture more afterslip, but M-rel does not correlate with observation start time or duration. Because M-rel estimates for a single earthquake can vary by an order of magnitude, we propose that modeling uncertainty currently presents a challenge for systematic afterslip analysis. Standardizing modeling practices may improve model comparability, and eventually allow for predictive afterslip models that account for mainshock and fault zone factors to be incorporated into aftershock hazard models.