Please use this identifier to cite or link to this item:
Title: An MCMC multiple point sources inversion scheme and its application to the 2016 Kumamoto Mw 6.2 earthquake
Authors: Shi, Qibin
Wei, Shengji
Chen, Meng
Keywords: Time-series Analysis
Inverse Theory
Science::Geology::Volcanoes and earthquakes
Issue Date: 2018
Source: Shi, Q., Wei, S., & Chen, M. (2018). An MCMC multiple point sources inversion scheme and its application to the 2016 Kumamoto Mw 6.2 earthquake. Geophysical Journal International, 215(2), 737-752. doi:10.1093/gji/ggy302
Series/Report no.: Geophysical Journal International
Abstract: While single point source is an oversimplified representation of medium to large earthquakes, finite fault models in many cases over parameterize the inversion due to the lack of sufficient near field data. Multiple point source solutions can fill in the gap in-between these two representations. Here, we propose a Markov-Chain-Monte-Carlo multiple point source inversion scheme, in combination with the advantage of the cut-and-paste technique, which cuts the seismogram into Pnl and surface portions and allows different time-shifts for each segment to align data with the synthetics. We apply the approach to the Mw 6.2 foreshock in the 2016 Kumamoto earthquake sequence by using the strong-motion observations within 100 km. We are able to perform the inversion at relatively high-frequency ranges (0.02–0.28 Hz for Pnl and 0.02–0.20 Hz for surface waves) with confidence on the velocity model built on the Mw 5.4 path calibration event. Our results show that the rupture was mainly composed of three subevents, with total duration of about 12 s and total Mw of 6.2. The strikes of three subevents agree well with surface fault mapping where the Futagawa Fault intersects the Hinagu Fault with ∼30° difference in strike. Our solution shows that the first subevent dips to the southeast, while the second and the third subevents, located ∼3 km the north and ∼4 km to the southwest of the first subevent, dip to the northwest. The focal mechanism of the first subevent shows remarkable agreement with the first-motion solution. The fault geometry also shows well consistence with the relocated aftershocks, which delineate a SE-dipping fault around the first subevent and two NW-dipping faults to the north and south, respectively, corresponding to the second and the third subevent. The sum of moment tensor of subevents that have not only different geometry but also rake angles shows strong Compensated Linear Vector Dipole (CLVD) component (48–50 per cent). With a local 1D crustal model, a full-moment-tensor inversion using regional long-period waveform data also detects strong CLVD component of this earthquake. In contrast, using the PREM model results in an almost pure double-couple solution. In short, we have precisely resolved the rupture process, the intricate fault geometry, and the strong CLVD component with strong-motion data. This highlights the importance of extracting the relatively high-frequency information from the waveform data and with accurate velocity model in seismic source analyses of large earthquakes.
ISSN: 0956-540X
DOI: 10.1093/gji/ggy302
Rights: © 2018 The Author(s). Published by Oxford University Press on behalf of The Royal Astronomical Society. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.
Fulltext Permission: open
Fulltext Availability: With Fulltext
Appears in Collections:ASE Journal Articles
EOS Journal Articles

Files in This Item:
File Description SizeFormat 
An MCMC multiple point sources inversion scheme ...pdf10.44 MBAdobe PDFThumbnail

Citations 20

checked on Sep 3, 2020


checked on Oct 17, 2020

Page view(s)

checked on Oct 24, 2020


checked on Oct 24, 2020

Google ScholarTM




Items in DR-NTU are protected by copyright, with all rights reserved, unless otherwise indicated.