Neville Palmer

Complex multifault rupture during the 2016 Mw 7.8 Kaikōura earthquake, New Zealand

An earthquake with a dozen faults The 2016 moment magnitude (Mw) 7.8 Kaikōura earthquake was one of the largest ever to hit New Zealand. Hamling et al. show with a new slip model that it was an incredibly complex event. Unlike most earthquakes, multiple faults ruptured to generate the ground shaking. A remarkable 12 faults ruptured overall, with the rupture jumping between faults located up to 15 km away from each other. The earthquake should motivate rethinking of certain seismic hazard models, which do not presently allow for this unusual complex rupture pattern. Science, this issue p. eaam7194 At least 12 faults spaced up to 15 kilometers apart ruptured during the magnitude 7.8 Kaikōura earthquake. INTRODUCTION On 14 November 2016 (local time), northeastern South Island of New Zealand was struck by a major moment magnitude (Mw) 7.8 earthquake. The Kaikōura earthquake was the most powerful experienced in the region in more than 150 years. The whole of New Zealand reported shaking, with widespread damage across much of northern South Island and in the capital city, Wellington. The earthquake straddled two distinct seismotectonic domains, breaking multiple faults in the contractional North Canterbury fault zone and the dominantly strike-slip Marlborough fault system. RATIONALE Earthquakes are conceptually thought to occur along a single fault. Although this is often the case, the need to account for multiple segment ruptures challenges seismic hazard assessments and potential maximum earthquake magnitudes. Field observations from many past earthquakes and numerical models suggest that a rupture will halt if it has to step over a distance as small as 5 km to continue on a different fault. The Kaikōura earthquake’s complexity defies many conventional assumptions about the degree to which earthquake ruptures are controlled by fault segmentation and provides additional motivation to rethink these issues in seismic hazard models. RESULTS Field observations, in conjunction with interferometric synthetic aperture radar (InSAR), Global Positioning System (GPS), and seismology data, reveal the Kaikōura earthquake to be one of the most complex earthquakes ever recorded with modern instrumental techniques. The rupture propagated northward for more than 170 km along both mapped and unmapped faults before continuing offshore at the island’s northeastern extent. A tsunami of up to 3 m in height was detected at Kaikōura and at three other tide gauges along the east coast of both the North and South Islands. Geodetic and geological field observations reveal surface ruptures along at least 12 major crustal faults and extensive uplift along much of the coastline. Surface displacements measured by GPS and satellite radar data show horizontal offsets of ~6 m. In addition, a fault-bounded block (the Papatea block) was uplifted by up to 8 m and translated south by 4 to 5 m. Modeling suggests that some of the faults slipped by more than 20 m, at depths of 10 to 15 km, with surface slip of ~10 m consistent with field observations of offset roads and fences. Although we can explain most of the deformation by crustal faulting alone, global moment tensors show a larger thrust component, indicating that the earthquake also involved some slip along the southern end of the Hikurangi subduction interface, which lies ~20 km beneath Kaikōura. Including this as a fault source in the inversion suggests that up to 4 m of predominantly reverse slip may have occurred on the subduction zone beneath the crustal faults, contributing ~10 to 30% of the total moment. CONCLUSION Although the unusual multifault rupture observed in the Kaikōura earthquake may be partly related to the geometrically complex nature of the faults in this region, this event emphasizes the importance of reevaluating how rupture scenarios are defined for seismic hazard models in plate boundary zones worldwide. Observed ground deformation from the 2016 Kaikōura, New Zealand, earthquake. (A and B) Photos showing the coastal uplift of 2 to 3 m associated with the Papatea block [labeled in (C)]. The inset in (A) shows an aerial view of New Zealand. Red lines denote the location of known active faults. The black box indicates the Marlborough fault system

New Zealand GPS velocity field: 1995–2013

ABSTRACT We collate nearly two decades of campaign GPS data gathered at over 900 sites throughout New Zealand to release a New Zealand nationwide GPS velocity field. The data span the entire North and South islands of New Zealand with a typical spacing of 10–20 km and a denser network (c. 2–8 km spacing) in the Wellington region, central Taupo Volcanic Zone and parts of the Arthurs Pass area. The dataset provides the most comprehensive-to-date view of crustal deformation within the Australia–Pacific plate boundary zone in the New Zealand region. We discuss the data acquisition, processing and derivation of the velocities and uncertainties. We also undertake corrections for earthquake displacements to obtain a velocity field that is largely representative of interseismic deformation between 1995 and 2013.

Crustal deformation and stress transfer during a propagating earthquake sequence: The 2013 Cook Strait sequence, central New Zealand

The 2013 Cook Strait earthquake sequence began on 18 July 2013 with two foreshocks of Mw 5.7 and Mw 5.8 and culminated in the Mw 6.6 Cook Strait and Lake Grassmere events on 21 July and 16 August, respectively. Located ∼50 km south of New Zealands capital, Wellington, the earthquakes generated the most significant ground shaking in the Wellington and Marlborough regions in recent decades. During the first event, located under Cook Strait, continuously recording GPS instruments across central New Zealand recorded up to 5 cm of horizontal displacement. Modeling suggests that the rupture was 25 km long with up to 90 cm of dextral strike slip. The second event, located 20 km to the southwest, caused displacements of up to 25 cm at GPS sites located around the Clifford bay area. In addition, two interferograms from RADARSAT-2 and TerraSAR-X showed up to 30 cm of line-of-sight displacement in the vicinity of Lake Grassmere. Modeling indicates predominantly dextral strike slip of up to 2.1 m. Coulomb Stress changes induced by the earlier foreshocks suggest that the Cook Strait event was triggered by the preceding events and that the Lake Grassmere event was subsequently triggered by the Cook Strait earthquake.

Off-axis magmatism along a subaerial back-arc rift: Observations from the Taupo Volcanic Zone, New Zealand.

A study of the growth of a large off-axis magma body along the Taupo Volcanic Zone. Continental rifting and seafloor spreading play a fundamental role in the generation of new crust. However, the distribution of magma and its relationship with tectonics and volcanism remain poorly understood, particularly in back-arc settings. We show evidence for a large, long-lived, off-axis magmatic intrusion located on the margin of the Taupo Volcanic Zone, New Zealand. Geodetic data acquired since the 1950s show evidence for uplift outside of the region of active extension, consistent with the inflation of a magmatic body at a depth of ~9.5 km. Satellite radar interferometry and Global Positioning System data suggest that there was an increase in the inflation rate from 2003 to 2011, which correlates with intense earthquake activity in the region. Our results suggest that the continued growth of a large magmatic body may represent the birth of a new magma chamber on the margins of a back-arc rift system.

First geodetic observations using new VLBI stations ASKAP-29 and WARK12M

We report the results of a successful 7-hour 1.4GHz Very Long Baseline Interferometry (VLBI) experiment using two new stations, ASKAP-29 located in Western Australia and WARK12M located on the North Island of New Zealand. This was the first geodetic VLBI observing session with the participation of these new stations. We have determined the positions of ASKAP-29 and WARK12M. Random errors on position estimates are 150-200mm for the vertical component and 40-50mm for the horizontal component. Systematicerrorscausedbytheunmodeledionospherepathdelaymayreach1.3mfortheverticalcomponent.