On 14th November 2016, the northeastern South Island of New Zealand was struck by a major Mw 7.8 earthquake. Field observations, in conjunction with InSAR, GPS, and seismology reveal this to be one of the most complex earthquakes ever recorded. The rupture propagated northward for more than 170 km along both mapped and unmapped faults, before continuing offshore at its northeastern extent. Geodetic and field observations reveal surface ruptures along at least 12 major faults, including possible slip along the southern Hikurangi subduction interface, extensive uplift along much of the coastline and widespread anelastic deformation including the ~8 m uplift of a fault-bounded block. This complex earthquake defies many conventional assumptions about the degree to which earthquake ruptures are controlled by fault segmentation, and should motivate re-thinking of these issues in seismic hazard models.One Sentence Summary: Major earthquake rips through evolving fault zone, defying conventional wisdom regarding the degree of fault segmentation during earthquake ruptures.

Supershear earthquakes with rupture velocity exceeding shear-wave speeds, previously observed in laboratory experiments and large strike-slip events, often have an initial sub-shear stage before they transition to supershear. In this study, integrated geophysical observations of the 2018 Mw 7.5 Palu, Indonesia earthquake, provide robust evidence of an early and persistent supershear rupture speed. Slowness-enhanced back-projection (SEBP) of teleseismic data provides a sharp image of the rupture process, consistently across multiple arrays. The inferred rupture path agrees with the surface rupture trace inferred from the net surface displacement field derived by sub-pixel InSAR image correlation. The SEBP results indicate a sustained rupture velocity of 4.1 km/s from the rupture initiation to the end, despite large fault bends. The persistent supershear speed is further validated by evidence of far-field Rayleigh Mach waves in regional seismograms. The short or absent supershear transition distance can be caused by high initial shear stress or short critical slip-weakening distance, and promoted by fault roughness near the hypocenter. Steady rupture propagation at a supershear speed considered to be unstable, lower than the Eshelby speed, could result from the presence of a damaged fault zone.

The 2015 Mw8.3 Illapel, Chile earthquake is the latest megathrust event on the central segment of that subduction zone. It generated strong ground motions and a large (up to 11 m runup) tsunami which prompted the evacuation of more than 1 million people in the first hours following the event. Observations during recent earthquakes suggest that these phenomena can be associated with rupture on different parts of the megathrust. The deep portion generates strong shaking while slow, large slip on the shallow fault is responsible for the tsunami. It is unclear whether all megathrusts can have shallow slip during coseismic rupture and what physical properties regulate this. Here we show that the Illapel event ruptured both deep and shallow segments with substantial slip. We resolve a kinematic slip model using regional geophysical observations and analyze it jointly with teleseismic backprojection. We find that the shallow and deep portions of the megathrust are segmented and have fundamentally different behavior. We forward calculate local tsunami propagation from the resolved slip and find good agreement with field measurements, independently validating the slip model. These results show that the central portion of the Chilean subduction zone has accumulated a significant shallow slip deficit and indicates that, given enough time, shallow slip might be possible everywhere along the subduction zone.

The September 2018, M w 7.5 Sulawesi earthquake occurring on the Palu-Koro strike-slip fault system was followed by an unexpected localized tsunami. We show that direct earthquakeinduced uplift and subsidence could have sourced the observed tsunami within Palu Bay. To this end, we use a physics-based, coupled earthquake-tsunami modeling framework tightly constrained by observations. The model combines rupture dynamics, seismic wave propagation, tsunami propagation and inundation. The earthquake scenario, featuring sustained supershear rupture propagation, matches key observed earthquake characteristics, including the moment magnitude, rupture duration, fault plane solution, teleseismic waveforms and inferred horizontal ground displacements. The remote stress regime reflecting regional transtension applied in the model produces a combination of up to 6 m left-lateral slip and up to 2 m normal slip on the straight fault segment dipping 65 East beneath Palu Bay. The time-dependent, 3D seafloor displacements are translated into bathymetry perturbations with a mean vertical offset of 1.5 m across the submarine fault segment. This sources a tsunami with wave amplitudes and periods that match those measured at the Pantoloan wave gauge and inundation that reproduces observations from field surveys. We conclude that a source related to earthquake displacements is probable and that landsliding may not have been the primary source of the tsunami. These results have important implications for submarine strike-slip fault systems worldwide. Physics-based modeling offers rapid response specifically in tectonic settings that are currently underrepresented in operational tsunami hazard assessment.

The 25 April 2015 M w 7.8 Gorkha earthquake caused more than 8000 fatalities and widespread building damage in central Nepal. The Italian Space Agency's COSMO-SkyMed Synthetic Aperture Radar (SAR) satellite acquired data over Kathmandu area four days after the earthquake and the Japan Aerospace Exploration Agency's Advanced Land Observing Satellite-2 SAR satellite for larger area nine days after the mainshock. We used these radar observations and rapidly produced damage proxy maps (DPMs) derived from temporal changes in Interferometric SAR coherence. Our DPMs were qualitatively validated through comparison with independent damage analyses by the National Geospatial-Intelligence Agency and the United Nations Institute for Training and Research's United Nations Operational Satellite Applications Programme, and based on our own visual inspection of DigitalGlobe's WorldView optical pre-versus postevent imagery. Our maps were quickly released to responding agencies and the public, and used for damage assessment, determining inspection/imaging priorities, and reconnaissance fieldwork.

No abstract

The 2016 Kumamoto earthquake sequence occurred on the Futagawa–Hinagu fault zone near the Aso volcano on Kyushu island. The sequence was initiated with two major (Mw ≥ 6.0) foreshocks, and the mainshock (Mw = 7.0) occurred 25 h after the second major foreshock. We combine GPS, strong motion, synthetic aperture radar images, and surface offset data in a joint inversion to resolve the kinematic rupture process of the mainshock and coseismic displacement of the foreshocks. The joint inversion results reveal a unilateral rupture process for the mainshock involving sequential rupture of four major asperities. The slip area of the foreshocks and mainshock and the aftershock loci form a detailed complementary pattern. The mainshock rupture terminates near the rim of the caldera, leaving a ~10 km long gap of aftershocks. This area is characterized by high temperature and low shear wave velocity, density, and resistivity, which may be related to the partially melted geothermal condition. Ductile material property near the volcano may act as a “material barrier” to the dynamic rupture. Topographic weight of the caldera increases compressional normal stress on the fault plane, which may behave as a “stress barrier.” Long‐term seismic hazard and deformation behaviors related to these two types of barriers are discussed in terms of the associated frictional mechanism. Significant postseismic creeps observed near the volcano area indicates a velocity strengthening frictional behavior near the rupture termination, which confirms that the “material barrier” mechanism is likely the dominant rupture termination mechanism.

Rapid shortening in convergent mountain belts is often accommodated by slip on faults at multiple levels in upper crust, but no geodetic observation of slip at multiple levels within hours of a moderate earthquake has been shown before. Here we show clear evidence of fault slip within a shallower thrust at 5–10 km depth in SW Taiwan triggered by the 2016 Mw 6.4 MeiNong earthquake at 15–20 km depth. We constrain the primary coseismic fault slip with kinematic modeling of seismic and geodetic measurements and constrain the triggered slip and fault geometry using synthetic aperture radar interferometry. The shallower thrust coincides with a proposed duplex located in a region of high fluid pressure and high interseismic uplift rate, and may be sensitive to stress perturbations. Our results imply that under tectonic conditions such as high‐background stress level and high fluid pressure, a moderate lower crustal earthquake can trigger faults at shallower depth.