PUBLICATIONS Journal of Geophysical Research: Solid Earth RESEARCH ARTICLE Transpressional Rupture Cascade of the 2016 Mw 10.1002/2017JB015168 7.8 Kaikoura Earthquake, New Zealand Key Points: Wenbin Xu1 , Guangcai Feng2, Lingsen Meng3 , Ailin Zhang3, Jean Paul Ampuero4 , • Complex coseismic ground 5 6 deformation can be explained by slip Roland Bürgmann , and Lihua Fang on six crustal fault segments 1 2 • Rupture process across multiple faults Department of Land Surveying and Geo-informatics, Hong Kong Polytechnic University, Hong Kong, China, School of 3 likely resulted from a triggering Geosciences and Info-Physics, Central South University, Changsha, China, Department of Earth Planetary and Space cascade between crustal faults Sciences, University of California, Los Angeles, CA, USA, 4Seismological Laboratory, California Institute of Technology, • Rupture speed was overall slow, but Pasadena, CA, USA, 5Department of Earth and Planetary Science, University of California, Berkeley, CA, USA, 6Institute of locally faster along individual fault segments Geophysics, China Earthquake Administration, Beijing, China Supporting Information: Abstract Large earthquakes often do not occur on a simple planar fault but involve rupture of multiple • Supporting Information S1 • Data Set S1 geometrically complex faults. The 2016 Mw 7.8 Kaikoura earthquake, New Zealand, involved the rupture of • Data Set S2 at least 21 faults, propagating from southwest to northeast for about 180 km. Here we combine space • Data Set S3 geodesy and seismology techniques to study subsurface fault geometry, slip distribution, and the kinematics of the rupture. Our finite-fault slip model indicates that the fault motion changes from predominantly Correspondence to: W. Xu, G. Feng, and L. Meng, right-lateral slip near the epicenter to transpressional slip in the northeast with a maximum coseismic surface [email protected]; displacement of about 10 m near the intersection between the Kekerengu and Papatea faults. Teleseismic [email protected]; back projection imaging shows that rupture speed was overall slow (1.4 km/s) but faster on individual fault [email protected] segments (approximately 2 km/s) and that the conjugate, oblique-reverse, north striking faults released the largest high-frequency energy. We show that the linking Conway-Charwell faults aided in propagation of Citation: rupture across the step over from the Humps fault zone to the Hope fault. Fault slip cascaded along the Xu, W., Feng, G., Meng, L., Zhang, A., Ampuero, J. P., Bürgmann, R., & Fang, L. Jordan Thrust, Kekerengu, and Needles faults, causing stress perturbations that activated two major (2018). Transpressional rupture cascade conjugate faults, the Hundalee and Papatea faults. Our results shed important light on the study of of the 2016 Mw 7.8 Kaikoura earthquake, earthquakes and seismic hazard evaluation in geometrically complex fault systems. New Zealand. Journal of Geophysical Research: Solid Earth, 123, 2396–2409. https://doi.org/10.1002/2017JB015168 1. Introduction Received 29 OCT 2017 Accepted 5 FEB 2018 On 14 November 2016, an earthquake with moment magnitude (Mw) 7.8 took place in the Kaikoura region in Accepted article online 12 FEB 2018 the Marlborough fault system (MFS), South Island, New Zealand (Figure 1) (Hamling et al., 2017; Kaiser et al., Published online 10 MAR 2018 2017). The region straddles a wide zone of active crustal deformation associated with the oblique collision between the Pacific and Australian tectonic plates (Van Dissen & Yeats, 1991). The seismic moment tensor solutions for the Kaikoura earthquake indicate an oblique-reverse focal mechanism. The hypocenter of the earthquake was located at about 15 km depth near the town of Culverden (Kaiser et al., 2017). Long-period seismological analysis of the earthquake indicates that the event initiated as a small strike-slip rupture during the first 60 s followed by a major 20 s burst of moment release, resulting in a long rupture duration (Duputel & Rivera, 2017). More than 2,000 aftershocks occurred within three days, four of which had Mw > 6 (Figure 1). A large number of shallow aftershocks occurred across a broad area north of Canterbury (Kaiser et al., 2017). Field investigations carried out by the Institute of Geological and Nuclear Sciences suggested that at least 14 major faults and seven minor faults including 11 previously unmapped faults ruptured from the Humps fault zone in the south to Cape Campbell ending offshore on the Needles fault (Figure 1) (Litchfield et al., 2016). The total length of the ruptured faults exceeds 150 km making it the largest historic event in the region. Field observations found that the largest surface fault offsets occurred along the Kekerengu fault with right-lateral slip reaching a maximum of about 11 m (Litchfield et al., 2016). Vertical displacements (6–8m) were observed on the previously unmapped Papatea fault with a conjugate geometry to the Kekerengu fault (Litchfield et al., 2016). Continuous and campaign Global Positioning System (GPS) stations captured the coseismic ground deformation (Hamling et al., 2017). Stations near the epicenter show that the ground displacements were dominated by horizontal motions. ©2018. American Geophysical Union. Multiple studies have been carried out to understand the complex faulting along the 2016 Kaikoura rupture All Rights Reserved. using different data sets and methods (e.g., field observations, geodesy, seismology, satellite optical and XU ET AL. 2396 Journal of Geophysical Research: Solid Earth 10.1002/2017JB015168 −41˚30' Australian Plate WF HST NF AWF CF AF AF Pacific Plate 2016/11/13 −42˚00' 11:52:45 KF JT PF 2016/11/13 13:31:25 HF PKF −42˚30' SJF HDF CCF 2016/11/13 HFZ 13:31:25 2016/11/14 HST 00:34:22 2016/11/13 11:02:56 −43˚00' 172˚30' 173˚00' 173˚30' 174˚00' 174˚30' Figure 1. The regional setting and location of the 2016 Kaikoura earthquake in the northern South Island, New Zealand. The red lines show the location of the modeled faults. The dark green lines show the fault ruptures mapped. The mapped major faults are shown in black. AF, Alpine fault; WF, Wairau fault; AWF, Awatere fault; CF, Clarence fault; NF, Needles fault; KF, Kekerengu fault; PF, Papatea fault; JT, Jordan Thrust; HF, Hope fault; HDF, Hundalee fault; SJF, Stone Jug fault; CCF, Conway-Charwell fault zone; HFZ, Humps fault zone; PKF, Point Kean Fault; HST, Hikurangi subduction thrust. The red focal mechanism represents the epicenter of 2016 Kaikoura earthquake determined by the U.S. Geological Survey. The green focal mechanisms represent four Mw > 6 aftershocks. The cyan dots represent 2383 relocated aftershocks from the first 5.5 days (Kaiser et al., 2017). The dashed purple box outlines the region shown in Figure 2. The inset shows the study area (red dashed rectangle). The grey and cyan rectangles show the coverage of ALOS-2 and Sentinel-1, respectively. The white lines indicate the plate boundary faults. The yellow line represents the geometry of the HST from Williams et al. (2013). radar imagery, and tsunami recordings). Kääb et al. (2017) showed the localized horizontal coseismic displacement field along the Kekerengu fault using images collected by the PlanetScope optical cubesat constellation. Kääb et al. (2017) observed ~10 m horizontal surface offsets on the Kekerengu fault and ~7 m on the Papatea fault. Clark et al. (2017) found highly variable coseismic coastal deformation ranging from À2.5 to 6.5 m from airborne light detection and ranging differencing. Hollingsworth et al. (2017) correlated Landsat-8 images from before and after the earthquake to obtain a complete horizontal coseismic displacement field covering the whole rupture. Hollingsworth et al. (2017) also analyzed seismic waveforms to determine fault slip on the Kekerengu fault and on a deeper offshore fault. Hamling et al. (2017) proposed a model involving slip on over 20 crustal fault segments and the Hikurangi subduction thrust to explain the observed coseismic displacements from both GPS and radar interferometry (interferometric synthetic aperture radar [InSAR]). Using a multiarray P wave back projection (BP) method, Zhang et al. (2017) showed the earthquake ruptured unilaterally along multiple faults to the northeast and found that the peak ground acceleration is well correlated with inferred release of short-period energy. Modeling teleseismic body waves and regional tsunami recordings, Bai et al. (2017) proposed simultaneous rupture on the subduction thrust and the upper crustal faults. This is also supported by joint XU ET AL. 2397 Journal of Geophysical Research: Solid Earth 10.1002/2017JB015168 inversion of teleseismic waveform and space geodetic data (Wang et al., 2018). Heidarzadeh and Satake (2017), however, suggested that the tsunami was likely triggered by submarine landslides and fault slip on the subduction thrust. Cesca et al. (2017) examined regional and teleseismic data to reconstruct the ruptured fault geometry. Cesca et al. (2017) found that the upper crustal faults are connected with a low-angle splay thrust fault rather than the Hikurangi subduction thrust. Holden et al. (2017) studied local strong motion and high-rate GPS data and suggested that the Hikurangi subduction thrust most likely did not slip during the earthquake. Using only 21 local strong motion data, Zheng et al. (2018) estimated the source rupture process on multiple fault segments from the kinematic waveform inversion. In this study, we combine seismic and geodetic data to study ground deformation of the 2016 Kaikoura earth- quake and its relation to the geometrical complexity of faults and the rupture process at depth. We use InSAR data and radar offset measurements from multiple platforms to generate a complete three-dimensional coseismic displacement field. We analyze geodetic data to identify the ruptured faults, to constrain the source parameters and fault geometry. We examine the resolving power of our geodetic data to answer whether onshore geodetic data are sensitive to fault slip on the Hikurangi subduction thrust and the Point Kean fault.