We constructed the rupture process model for the 2016 Kumamoto, Japan, earthquake from broadband teleseismic body waveforms (P-waves) by using a novel waveform inversion method that takes into account the uncertainty of Green's function. The estimated source parameters are: seismic moment = 5.1 × 1019 Nm (Mw = 7.1), fault length = 40 km, and fault width = 15 km. The mainshock rupture mainly propagated northeastward from the epicenter, for about 30 km, along an active strike-slip fault. The rupture propagation of the mainshock decelerated and terminated near the southwest side of the Aso volcano; the aftershock activity was low around the northeastern edge of the major slip area. Our results suggest that the rupture process of the mainshock and the distribution of aftershocks were influenced by the high-temperature area around the magma chamber of Mt. Aso.

International audienc

SUMMARY Teleseismic waveforms contain information on fault slip evolution during an earthquake, as well as on the fault geometry. A linear finite-fault inversion method is a tool for solving the slip-rate function distribution under an assumption of fault geometry as a single or multiple-fault-plane model. An inappropriate assumption of fault geometry would tend to distort the solution due to Green’s function modelling errors. We developed a new inversion method to extract information on fault geometry along with the slip-rate function from observed teleseismic waveforms. In this method, as in most previous studies, we assumed a flat fault plane, but we allowed arbitrary directions of slip not necessarily parallel to the assumed fault plane. More precisely, the method represents fault slip on the assumed fault by the superposition of five basis components of potency-density tensor, which can express arbitrary fault slip that occurs underground. We tested the developed method by applying it to real teleseismic P waveforms of the MW 7.7 2013 Balochistan, Pakistan, earthquake, which is thought to have occurred along a curved fault system. The obtained spatiotemporal distribution of potency-density tensors showed that the focal mechanism at each source knot was dominated by a strike-slip component with successive strike angle rotation from 205° to 240° as the rupture propagated unilaterally towards the south-west from the epicentre. This result is consistent with Earth’s surface deformation observed in optical satellite images. The success of the developed method is attributable to the fact that teleseismic body waves are not very sensitive to the spatial location of fault slip, whereas they are very sensitive to the direction of fault slip. The method may be a powerful tool to extract information on fault geometry along with the slip-rate function without requiring detailed assumptions about fault geometry.

The rupture process of the 2014 Iquique, Chile earthquake is inverted from teleseismic P wave data applying a novel formulation that takes into account the uncertainty of Green's function, which has been a major error source in waveform inversion. The estimated seismic moment is 1.5 × 10 21 Nm (Mw = 8.1), associated with a 140 km long and 140 km wide fault rupture along the plate interface. The source process is characterized by unilateral rupture propagation. During the first 20 s, the dynamic rupture front propagated from the hypocenter to the large asperity located about 50 km southward, crossing a remarkably active foreshock area at high velocity (of about 3.0 km/s), but small and irregular seismic moment release rate. Our result may suggest that the 20 s long initial phase was influenced by the stress drop due to the foreshock activity near the main shock hypocenter. Moreover, the 2 week long swarm-like foreshock activity migrating roughly at 5 km/day toward the main shock hypocenter, and possibly associated slow slip, contributed to the stress accumulation prior to the Mw 8.1 megaquake. The main shock initial rupture phase might have triggered the rupture of the large asperity, which had large fracture energy.

We compared spatiotemporal slip‐rate and high‐frequency (around 1 Hz) radiation distributions from teleseismic P wave data to infer the seismic rupture process of the 2015 Gorkha, Nepal, earthquake. For these estimates, we applied a novel waveform inversion formulation that mitigates the effect of Green's functions uncertainty and a hybrid backprojection method that mitigates contamination by depth phases. Our model showed that the dynamic rupture front propagated eastward from the hypocenter at 3.0 km/s and triggered a large‐slip event centered about 50 km to the east. It also showed that the large‐slip event included a rapid rupture acceleration event and an irregular deceleration of rupture propagation before the rupture termination. Heterogeneity of the stress drop or fracture energy in the eastern part of the rupture area, where aftershock activity was high, inhibited rupture growth. High‐frequency radiation sources tended to be in the deeper part of the large‐slip area, which suggests that heterogeneity of the stress drop or fracture energy there may have contributed to the damage in and around Kathmandu.

Rupture propagation of an earthquake strongly influences potentially destructive ground shaking. Variable rupture behaviour is often caused by complex fault geometries, masking information on fundamental frictional properties. Geometrically smoother ocean transform fault (OTF) plate boundaries offer a favourable environment to study fault zone dynamics because strain is accommodated along a single, wide zone (up to 20 km width) offsetting homogeneous geology comprising altered mafic or ultramafic rocks. However, fault friction during OTF ruptures is unknown: no large (Mw>7.0) ruptures had been captured and imaged in detail. In 2016, we recorded an Mw 7.1 earthquake on the Romanche OTF in the equatorial Atlantic on nearby seafloor seismometers. We show that this rupture had two phases: (1) up and eastwards propagation towards the weaker ridge-transform intersection (RTI), then (2) unusually, back-propagation westwards at super-shear speed toward the fault's centre. Deep slip into weak fault segments facilitated larger moment release on shallow locked zones, highlighting that even ruptures along a single distinct fault zone can be highly dynamic. The possibility of reversing ruptures is absent in rupture simulations and unaccounted for in hazard assessments.

Teleseismic waveforms contain information on fault slip evolution during an earthquake, as well as on the fault geometry. A linear finite-fault inversion method is a tool for solving the slip-rate function distribution under an assumption of fault geometry as a single or multiple-fault-plane model. An inappropriate assumption of fault geometry would tend to distort the solution due to Green's function modelling errors. We developed a new inversion method to extract information on fault geometry along with the slip-rate function from observed teleseismic waveforms. In this method, as in most previous studies, we assumed a flat fault plane, but we allowed arbitrary directions of slip not necessarily parallel to the assumed fault plane. More precisely, the method represents fault slip on the assumed fault by the superposition of five basis components of potency-density tensor, which can express arbitrary fault slip that occurs underground. We tested the developed method by applying it to real teleseismic P waveforms of the Mw 7.7 2013 Balochistan, Pakistan, earthquake, which is thought to have occurred along a curved fault system. The obtained spatiotemporal distribution of potency-density tensors showed that the focal mechanism at each source knot was dominated by a strike-slip component with successive strike angle rotation from 205° to 240° as the rupture propagated unilaterally towards the south-west from the epicentre. This result is consistent with Earth's surface deformation observed in optical satellite images. The success of the developed method is attributable to the fact that teleseismic body waves are not very sensitive to the spatial location of fault slip, whereas they are very sensitive to the direction of fault slip. The method may be a powerful tool to extract information on fault geometry along with the slip-rate function without requiring detailed assumptions about fault geometry.

No abstract