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.

Bathymetric survey images structure off Sumatra

Limits of the seismogenic zone in the epicentral region of the 26 December 2004 great Sumatra-Andaman earthquake: Results from seismic refraction and wide-angle reflection surveys and thermal modeling

Multi-fault rupture and successive triggering during the 2012 Mw 8.6 Sumatra offshore earthquake

Seismic imaging of forearc backthrusts at northern Sumatra subduction zone

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

Rupture of the 2004 Sumatra-Andaman Earthquake Inferred from Direct P-wave Imaging

Rheological Structure Beneath the Southern Sumatra Constrained from Postseismic Deformation of the 2007 Mw8.4 Sumatra Earthquake

Teleseismic traveltime tomography of northern Sumatra

Tsunami scenario triggered by a submarine landslide offshore of northern Sumatra Island and its hazard assessment

Teleseismic Double‐difference Relocation of Earthquakes along the Sumatra‐Andaman Subduction Zone Using a 3‐D Model

Imaging of 3-D seismic velocity structure of Southern Sumatra region using double difference tomographic method

Mantle Wedge Heterogeneities and Afterslip Distribution Following the 2004 Mw 9.2 and 2005 Mw 8.6 Sumatra Earthquakes

Similarities in the rupture process and cascading asperities between neighboring fault patches and seismic implications: The 2002–2009 Sumbawa (Indonesia) earthquakes with moment magnitudes of 6.2–6.6

Spatiotemporal characteristics of the energy radiation sources of the three great earthquakes near Sumatra Island in September 2007

Source Parameters of the Great Sumatran Megathrust Earthquakes of 1797 and 1833 Inferred from Coral Microatolls

Rupture and Frequency-Dependent Seismic Radiation of the 2012mw8.6 Sumatra Strike-Slip Earthquake

Slab Morphology Beneath Northern Sumatra Revealed by Regional and Teleseismic Traveltime Tomography

Rupture process of the 2010<i>M<sub>w</sub></i>7.8 Mentawai tsunami earthquake from joint inversion of near-field hr-GPS and teleseismic body wave recordings constrained by tsunami observations

Extent, duration and speed of the 2004 Sumatra-Andaman earthquake imaged by the Hi-Net array

A note on the seismicity of Sumatra, western Sunda Arc, Indonesia, in relation to the potential for back-arc thrusting