Our method significantly reduces the number of costly simulations required, allowing us to precisely quantify uncertainties and relationships between how faults slip and how surrounding rocks deform. We find a strong link between fault slip properties and surrounding rock deformation: areas with less frictional weakening on the fault coincide with areas where surrounding rocks experience more extensive damage. These findings align with geological evidence that suggests the fault is less mature and more fragmented in some regions.

We also identify a shallow slip deficit, where less fault slip occurs near the Earth's surface than deeper down, consistent with previous studies. Our results highlight the effectiveness of combining realistic earthquake simulations with observational data, providing a more accurate and comprehensive understanding of earthquake behavior and hazard potential.