Development of an inversion method to extract information on fault geometry from teleseismic data

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 M-W 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 degrees to 240 degrees 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.