Abstract:The distribution of slip during an earthquake and how it propagates among faults in the subduction system play a major role in seismic and tsunami hazards, yet they are poorly understood because offshore observations are often lacking. Here we derive the slip distribution and rupture evolution during the 2016 Mw 7.9 Kaikōura (New Zealand) earthquake that reconcile the surface rupture, space geodetic measurements, seismological and tsunami waveform records. We use twelve fault segments, with eleven in the crust and one on the megathrust interface, to model the geodetic data and match the major features of the complex surface ruptures. Our modeling result indicates that a large portion of the moment is distributed on the subduction interface, making a significant contribution to the far field surface deformation and teleseismic body waves. The inclusion of local strong motion and teleseismic waveform data in the joint inversion reveals a unilateral rupture towards northeast with a relatively low averaged rupture speed of ∼1.5 km/s. The first 30 s of the rupture took place on the crustal faults with oblique slip motion and jumped between fault segments that have large differences in strike and dip. The peak moment release occurred at ∼65 s, corresponding to simultaneous rupture of both plate interface and the overlying splay faults with rake angle changes progressively from thrust to strike-slip. The slip on the Papatea fault produced more than 2 m of offshore uplift, making a major contribution to the tsunami at the Kaikōura station, while the northeastern end of the rupture can explain the main features at the Wellington station. Our inversions and simulations illuminate complex up-dip rupture behavior that should be taken into consideration in both seismic and tsunami hazard assessment. The extreme complex rupture behavior also brings new challenges to the earthquake dynamic simulations and understanding the physics of earthquakes.

Rupture Process of the 2011 Tohoku Earthquake from the Joint Inversion of Teleseismic and Gps Data

The 2016 Kaikōura Earthquake: Simultaneous Rupture of the Subduction Interface and Overlying Faults

Inversion of High-Rate (1 Sps) GPS Data for Rupture Process of the 11 March 2011 Tohoku Earthquake (Mw9.1)

Tomography of the Source Zone of the Great 2011 Tohoku Earthquake

Source Rupture Models for the Mw 9.0 2011 Tohoku Earthquake from Joint Inversions of High-Rate Geodetic and Seismic Data

Quick inversion of the rupture process of the 2009 March 19 Tonga earthquake

Fast Inversion of Rupture Process of the 14 April 2010 Yushu, Qinghai, Earthquake

Three sub-events composing the 2011 off the Pacific coast of Tohoku Earthquake ( M w 9.0) inferred from rupture imaging by back-projecting teleseismic P waves

Structural Heterogeneity Controlled Rupture Process of the 2021 Mw 7.1 Fukushima, Japan, Earthquake Revealed by Joint Inversion of Seismic and Geodetic Data

Fast inversion of the rupture process of 2009 July 15 off west coast of South Island of New Zealand earthquake

Inversion of Earthquake Rupture Process: Theory and Applications(*)

Rapidness and Robustness of Finite-Source Inversion of the 2011 Mw 9.0 Tohoku Earthquake by an Elliptical-Patches Method Using Continuous GPS and Acceleration Data

Fast in version of the rupture process of 2010 January 12 Haiti earthquake

Spatio-temporal Variation of the Source Mechanism of the 2008 Great Wenchuan Earthquake

Fast inversion of the rupture process of 27 February 2010 Chile earthquake

Automatic Inversions of Strong‐Motion Records for Finite‐Fault Models of Significant Earthquakes in and Around Japan

Fast inversion of the rupture process of 2009 September 29 Samoa Islands region earthquake

Fast inversion of rupture process based on strong motion data and the feasibility of its automation

Inversion of Near-Field Waveform Data for Earthquake Source Rupture Process ( II ): Inversion of Chi Chi Earthquake Source, 20 September, 1999, Taiwan

Fast inversion of the rupture process of 2009 October 7 Vanuatu earthquake

Estimating the rupture process of the 8 August 2017 Jiuzhaigou earthquake by inverting strong-motion data with IDS method