Modeling Earthquake Sequences along the Manila Subduction Zone: Effects of Three-Dimensional Fault Geometry

Hongyu Yu,Yajing Liu,Hongfeng Yang,Jieyuan Ning

DOI: https://doi.org/10.1016/j.tecto.2018.01.025

IF: 3.66

2018-01-01

Tectonophysics

Abstract:To assess the potential of catastrophic megathrust earthquakes (M-w > 8) along the Manila Trench, the eastern boundary of the South China Sea, we incorporate a 3D non-planar fault geometry in the framework of rate-state friction to simulate earthquake rupture sequences along the fault segment between 15(center dot)N-19(center dot)N of northern Luzon Our simulation results demonstrate that the first-order fault geometry heterogeneity, the transitionalsegment (possibly related to the subducting Scarborough seamount chain) connecting the steeper south segment and the flatter north segment, controls earthquake rupture behaviors The strong along-strike curvature at the transitional-segment typically leads to partial ruptures of M-w-8 3 and similar to 8.similar to 7.8 al(o)similar to g the southern and northern segments respectively The entire fault occasionally ruptures in M-w similar to 8.8 ev(e)similar to ts when the cumulative stress in the transitional-segment is sufficiently high to overcome the geometrical inhibition. Fault shear stress evolution, represented by the S-ratio, is clearly modulated by the width of seismogemc zone niW). At a constant plate convergence rate, a larger W indicates on average lower mterseisinic stress loading rate and longer rupture recurrence period, and could slow down or sometimes stop ruptures that initiated from a narrower portion Moreover, the modeled mterseisinic slip rate before whole-fault rupture events is comparable with the coupling state that was inferred from the interplate seismicity distribution, suggesting the Manila trench could potentially rupture in a M8 + earthquake.