Journal of Geophysical Research - Solid Earth

Similar documents
Shattering a plate boundary: the 2016 Mw 7.8 Kaikōura earthquake

Buried Strike Slip Faults: The 1994 and 2004 Al Hoceima, Morocco Earthquakes.

The Seismic Hazardscape of New Zealand

Journal of Geophysical Research Letters Supporting Information for

The March 11, 2011, Tohoku-oki earthquake (Japan): surface displacement and source modelling

Case Study 1: 2014 Chiang Rai Sequence

Calculation of Focal mechanism for Composite Microseismic Events

Journal of Geophysical Research (Solid Earth) Supporting Information for

Case Study of Japan: Crustal deformation monitoring with GNSS and InSAR

Basics of the modelling of the ground deformations produced by an earthquake. EO Summer School 2014 Frascati August 13 Pierre Briole

to: Interseismic strain accumulation and the earthquake potential on the southern San

Haiti Earthquake (12-Jan-2010) co-seismic motion using ALOS PALSAR

Two Contrasting InSAR Studies of Recent Earthquakes in Tibet

Tsunami Simulation of 2009 Dusky Sound Earthquake in New Zealand

Learning goals - January 16, Describe the geometry of a fault (1) orientation of the plane (strike and dip) (2) slip vector

Inversion of Earthquake Rupture Process:Theory and Applications

Joint inversion of InSAR and broadband teleseismic waveform data with ABIC: application to the 1997 Manyi, Tibet earthquake

3-Dimension Deformation Mapping from InSAR & Multiaperture. Hyung-Sup Jung The Univ. of Seoul, Korea Zhong Lu U.S. Geological Survey, U.S.A.

Rupture Characteristics of Major and Great (M w 7.0) Megathrust Earthquakes from : 1. Source Parameter Scaling Relationships

Coseismic fault slip of the 2008 M w 7.9 Wenchuan earthquake estimated from InSAR and GPS measurements

GEO-HAZARDS STUDIES BY PALSAR INTERFEROMETRY PI No: 390 Xiaoli Ding 1, Lei Zhang 1, Guangcai Feng 1, Mi Jiang 1, and Jun Hu 2, 1 The Hong Kong Polytec

SUPPLEMENTARY INFORMATION

Ground surface deformation of L Aquila. earthquake revealed by InSAR time series

Preliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code.

This is a repository copy of Complex multi-fault rupture during the 2016 Mw 7.8 Kaikōura earthquake, New Zealand.

Supplementary Material

3D temporal evolution of displacements recorded on Mt. Etna from the 2007 to 2010 through the SISTEM method

INVESTIGATION OF EARTHQUAKE-CYCLE DEFORMATION IN TIBET FROM ALOS PALSAR DATA PI 168 Roland Bürgmann 1, Mong-Han Huang 1, Isabelle Ryder 2, and Eric Fi

THE VERTICAL COSEISMIC DEFORMATION FIELD OF THE WENCHUAN EARTHQUAKE BASED ON THE COMBINATION OF GPS AND INSAR

Supporting information for the main manuscript: Strain budget of the Ecuador-Colombia subduction zone: a stochastic view

crustal structure experiment beneath Wairarapa - Wellington area: results from SAHKE

Urban land and infrastructure deformation monitoring by satellite radar interferometry

SOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION

Was the February 2008 Bukavu seismic sequence associated with magma intrusion?

14 S. 11/12/96 Mw S. 6/23/01 Mw S 20 S 22 S. Peru. 7/30/95 Mw S. Chile. Argentina. 26 S 10 cm 76 W 74 W 72 W 70 W 68 W

Short Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood, California, Earthquake

DETECTION OF CRUSTAL MOVEMENTS IN THE 2011 TOHOKU, JAPAN EARTHQUAKE USING HIGH-RESOLUTION SAR DATA

Auxiliary Material for The 2010 Maule, Chile earthquake: Downdip rupture limit revealed by space geodesy

Burst overlapping of ALOS-2 PALSAR-2 ScanSAR-ScanSAR interferometry

Sendai Earthquake NE Japan March 11, Some explanatory slides Bob Stern, Dave Scholl, others updated March

The 2003, M W 7.2 Fiordland Earthquake, and its nearsource aftershock strong motion data

Modelling Strong Ground Motions for Subduction Events in the Wellington Region, New Zealand

Co-seismic slip from the July 30, 1995, M w 8.1 Antofagasta, Chile, earthquake as constrained by InSAR and GPS observations

Tomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima

Characterization of stress changes in subduction zones from space- and ground-based geodetic observations

Dear editors and reviewer(s), thank for your comments and suggestions. Replies as follows:

GEOPHYSICAL RESEARCH LETTERS, VOL. 37, L02304, doi: /2009gl041835, 2010

Data Repository of Paper: The role of subducted sediments in plate interface dynamics as constrained by Andean forearc (paleo)topography

Geodesy (InSAR, GPS, Gravity) and Big Earthquakes

Geometry of co-seismic surface ruptures and tectonic meaning of the 23 October 2011 Mw 7.1 Van earthquake (East Anatolian Region, Turkey)

Susitna-Watana Hydroelectric Project (FERC No ) Site-Specific Seismic Hazard Study Study Plan Section 16.6

The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets

A preliminary report on the 2018 Lombok region Indonesia earthquakes. August 9. Athanassios Ganas, Varvara Tsironi and Sotirios Valkaniotis

Radar Remote Sensing: Monitoring Ground Deformations and Geohazards from Space

4600 Rickenbacker Cswy, Miami 33149, FL. U.S.A., (2)

Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source

JCR (2 ), JGR- (1 ) (4 ) 11, EPSL GRL BSSA

Vertical to Horizontal (V/H) Ratios for Large Megathrust Subduction Zone Earthquakes

Received 17 June 2011; revised 20 March 2013; accepted 7 June 2013; published 27 August 2013.

Jocelyn Karen Campbell

SUPPLEMENTARY INFORMATION

RELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE

Case Study 2: 2014 Iquique Sequence

Tomographic imaging of P wave velocity structure beneath the region around Beijing

27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies

Coseismic slip model of the 2007 August Pisco earthquake (Peru) as constrained by Wide Swath radar observations

LECTURES 10 and 11 - Seismic Sources Hrvoje Tkalčić

Introduction to Displacement Modeling

27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies

The April 6 th 2009, L Aquila (Italy) earthquake: DInSAR analysis and seismic source model inversion

Rupture Process of the Great 2004 Sumatra-Andaman Earthquake

EXTENDING INTERFEROMETRIC SYNTHETIC APERTURE RADAR MEASUREMENTS FROM ONE TO TWO DIMENSIONS

Shallow intraplate earthquakes in Western Australia observed by Interferometric Synthetic Aperture Radar

Earthquakes and Seismotectonics Chapter 5

Scientific Research on the Cascadia Subduction Zone that Will Help Improve Seismic Hazard Maps, Building Codes, and Other Risk-Mitigation Measures

Contents of this file

The Mw 6.2 Leonidio, southern Greece earthquake of January 6, 2008: Preliminary identification of the fault plane.

Modeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation

Widespread Ground Motion Distribution Caused by Rupture Directivity during the 2015 Gorkha, Nepal Earthquake

NUMERICAL SIMULATIONS FOR TSUNAMI FORECASTING AT PADANG CITY USING OFFSHORE TSUNAMI SENSORS

Available online at ScienceDirect. Procedia Technology 16 (2014 )

Volcanoes Session. Chairs: Y. Fukushima & G. Puglisi

2 Model and Methods We employ a fully-dynamic rupture simulation using the Fast-Domain-Partitioning Boundary Integral equation method (FDP-BIEM)

Study on the feature of surface rupture zone of the west of Kunlunshan pass earthquake ( M S 811) with high spatial resolution satellite images

9th Workshop on Three-Dimensional Modelling of Seismic Waves Generation, Propagation and their Inversion

Crustal deformation by the Southeast-off Kii Peninsula Earthquake

NOTES AND CORRESPONDENCE Segmented Faulting Process of Chelungpu Thrust: Implication of SAR Interferograms

Line of Sight Displacement from ALOS-2 Interferometry: Mw 7.8 Gorkha Earthquake and Mw 7.3 Aftershock

Apparent Slow Oceanic Transform Earthquakes Due to Source Mechanism Bias

Regional Geodesy. Shimon Wdowinski. MARGINS-RCL Workshop Lithospheric Rupture in the Gulf of California Salton Trough Region. University of Miami

7 Ground Motion Models

Inquiry: Sumatran earthquakes with GPS Earth Science Education

FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT THE GEYSERS GEOTHERMAL FIELD

SOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION

DETECTION OF CRUSTAL DEFORMATION OF THE NORTHERN PAKISTAN EARTHQUAKE BY SATELLITE DATA. Submitted by Japan **

3D MODELING OF EARTHQUAKE CYCLES OF THE XIANSHUIHE FAULT, SOUTHWESTERN CHINA

Case study of Japan: Reference Frames in Practice

Geophysical Research Letters. Supporting Information for

Incorporating simulated Hikurangi subduction interface spectra into probabilistic hazard calculations for Wellington

Representative ground-motion ensembles for several major earthquake scenarios in New Zealand

Transcription:

Journal of Geophysical Research - Solid Earth Supporting Information for Transpressional Rupture Cascade of the 2016 M w 7.8 Kaikoura Earthquake, New Zealand Wenbin Xu 1*, Guangcai Feng 2*, Lingsen Meng 3*, Ailin Zhang 3, Jean Paul Ampuero 4, Roland Bürgmann 5, Lihua Fang 6 1Department of Land Surveying and Geo-informatics, Hong Kong Polytechnic University, China 2School of Geosciences and Info-Physics, Central South University, Changsha 410083, Hunan, China 3Department of Earth Planetary and Space Sciences, University of California, Los Angeles, CA 90025, USA 4Seismological Laboratory, California Institute of Technology, 1200 E. California Blvd., Pasadena, California 91125, USA 5Department of Earth and Planetary Science, University of California, Berkeley, California, USA 6Institute of Geophysics, China Earthquake Administration, Beijing 100081, China *e-mail: wenbin.xu@polyu.edu.hk; fredgps@csu.edu.cn; meng@epss.ucla.edu 1

Contents of this file Figures S1 to S13 Tables S1 to S3 Introduction This document contains supplementary figures and tables. Figures S1 shows a subset of China Array stations. Figure S2 shows back projections of aftershocks at different depths. Figure S3 shows synthetic rupture scenarios of varying focal mechanisms. Figure S4-S7 show the coseismic InSAR data. Figure S8 shows the three dimensional surface displacement field. Figure S9 shows the difference between the coastal uplift lidar data (Clark et al., 2017) and the vertical displacement. Figure S10- S11 show the resolvability of InSAR data on the upper crustal faults and the Hikurangi subduction thrust, respectively. Figure S12 shows the complete GPS data and modeling results based on model I and model II. Figure S13 shows the checkerboard test for slip on the Point Kean fault of model III. Table S1 contains the SAR data used in the study. Table S2 summarizes the data used in modeling. Table S3 shows the comparison between estimated dip and rake angles and the fault database. 2

Figure S1 A subset of China Array stations involved in the back-projection analysis. 3

Figure S2 Back projections of aftershocks assuming different depths. Stars indicate the catalogue locations of the events. Colors denote events and symbols represent focal depths assumed in the BP. We back-project aftershock data at 10 km, 15 km, and 20 km and compare the obtained horizontal locations. We find that the location variation due to depth uncertainty is a minor factor compared with the spatial biases due to path effect. 1

Figure S3 Synthetic rupture scenarios of varying focal mechanisms. BP radiators (circles) are sized by power and color-coded by time. Beach balls are also color-coded by time. Each beach ball indicates the corresponding mechanisms of egf for each particular input subsource. 2

Figure S4 Coseismic ALOS-2 data of the 2016 Kaikoura earthquake and the modelling result. First row: Observed line-of-sight (LOS) displacement map. Middle row: Model I prediction. Last row: Residual. The black lines represent the fault geometry used in the modelling. Color scale is the same for all subpanels. 3

Figure S5 Coseismic descending ALOS-2 offsets of the 2016 Kaikoura earthquake and the modelling result. Range offsets are shown in the left column and azimuth offsets shown in right column. Negative range-offset values indicate the ground motion away from the satellite. Positive azimuth-offset values indicate the ground motion along the flight direction of the satellite. The same color scale is used for all subpanels. 4

Figure S6 Same as Fig. S5 but for the descending Sentinel-1 offsets. 5

Figure S7 Same as Fig. S5 but for the ascending Sentinel-1 offsets. 6

Figure S8 Three dimensional surface displacement field derived from image offsets of the Sentinel-1 and ALOS-2 data (first row). Model I predictions (second row). Residuals (last row). 7

Figure S9 The comparison between the downsampled coastal uplift lidar data (Clark et al., 2017) and our vertical displacement. The systematic difference of ~0.4 cm between these two types of data is probably related to their spatial resolution difference and accuracy. 8

Figure S10 (a). The checkerboard test of model I. The result demonstrates that InSAR data have a very good resolvability for shallow fault slip, but poorer resolvability for deep slip. (b). Slip error estimated from 100 perturbed datasets generated using the Monte Carlo method. 9

Figure S11 (a). The checkerboard test for slip on Hikurangi subduction thrust of model II. The results (b), (c) demonstrate that InSAR data have a poor resolution for deep slip on the subduction thrust. (d). Slip error estimated from 100 perturbed datasets generated using the Monte Carlo method. 10

Figure S12 The complete coseismic GPS data (black) of the 2016 Kaikoura earthquake and the modelling result (red). The first and second row showing the horizontal and vertical GPS displacements, respectively. 11

Figure S13 Same as Fig. S10 but for the checkerboard test for slip on the Point Kean fault of model III. The results (b), (c) demonstrate that InSAR data have a poor resolution for slip on the offshore Point Kean fault. 12

Table S1 SAR images used in this study Sensor Path Master date (yyyy/mm/dd) Slave date (yyyy/mm/dd) Before (days) After (days) Perp. Baseline (m) Orbit ALOS-2 195 2016/10/23 2016/12/04 21 20-150 Desc. ALOS-2 102 2016/08/11 2016/12/01 91 17 130 Asc. Sentinel-1 73 2016/09/15 2016/11/16 58 2-10 Desc. Sentinel-1 52 2016/11/03 2016/11/15 11 1-10 Asc. 13

Table S2 Dataset used in generating three dimensional surface displacement and modeling Data set Reduced data points Uncertainty (m 2 ) Des. ALOS-2 4572 0.04 Asc. ALOS-2 5135 0.04 Des. ALOS-2 Azo 1108 0.43 Des. ALOS-2 Rng 691 0.28 Des. S-1 Azo 428 1.21 Des. S-1 Rng 927 0.81 Asc. S-1 Azo 420 0.85 Asc. S-1 Rng 1137 0.37 14

Table S3 Comparison with dip and rake angles from the fault database (Litchfield et al., 2014) and this study Fault Database This study, model I Dip [min, max], Best Rake [min, max], Best Dip [min, max], Best Rake [min, max], Best Hope Conway [60, 80], 70 [169, 176], 173 70 [150, 172], 162 Jordan [28, 48], 37 [90, 110], 100 [30, 50] [145, 161], 153 Kekerengu 1 [50, 70], 60 [125, 145], 135 50 [141, 160], 151 Kekerengu 2 [70, 90], 80 [143, 163], 153 50 [150, 172], 160 Needles [70, 90], 80 [147, 167], 157 [30, 50] [150, 172], 160 Hope-Jordan Thrust- Kekerengu- Needles [28, 90] [90, 176], 144 [30, 50] [130, 174], 158 Hundalee [40, 70], 55 [90, 94], 90 50 [64, 90], 75 15