Characterized source model of the M7.3 2016 Kumamoto earthquake by the 3D reciprocity GFs inversion with special reference to the velocity pulse at KMMH16

• Published in: Earth, planets, and space Vol. 75; no. 1; pp. 16 - 18
• Main Authors: Petukhin, Anatoly, Kawase, Hiroshi, Nagashima, Fumiaki, Ito, Eri
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer Berlin Heidelberg 01.12.2023; Springer; Springer Nature B.V; SpringerOpen
• Subjects: 2016 AD; 2016 Kumamoto earthquake; 4. Seismology; Analysis; Characterized source model; Earth and Environmental Science; Earth Sciences; Earthquake damage; Earthquakes; Effects of Surface Geology on Seismic Motion (ESG): General State-of-Research; Environmental aspects; Epidermal growth factor; Fault lines; Faults (Geology); Geology; Geophysics/Geodesy; Green's function; Green's functions; Ground motion; Japan; Mathematical models; Modelling; Parameters; Recipe for strong-motion prediction; Reciprocity; Reciprocity method; Residential areas; Search methods; Seismic activity; Seismology; Slip velocity; Source inversion; Strong-motion generation area; Velocity; Waveforms; Japan; Source inversion; Characterized source model; Source time function; 2016 Kumamoto earthquake; Reciprocity method; Recipe for strong-motion prediction; Strong-motion generation area
• ISSN: 1880-5981; 1343-8832; 1880-5981
• DOI: 10.1186/s40623-023-01768-w

The 2016 Kumamoto earthquakes caused severe damage centering on the Mashiki residential area. The velocity waveforms at station KMMH16 in Mashiki, during the M 7.3 mainshock, show large pulses. We found that severe damage in Mashiki may be the result of the strong westward velocity pulse. The question raised is how the near-fault ground motions with strong velocity pulse at KMMH16 were generated during the mainshock. We focus on the characterized source model with Strong Motion Generation Areas (SMGA). Empirical Green’s function (EGF) method is widely used for source modeling in this case. However, in case that the target site is located just near the fault in nodal plane of source mechanism (like KMMH16), mechanism of the EGF event should perfectly fit mechanism of the mainshock, which is a rare case. Therefore, instead of using EGFs, we used theoretical 3D Green’s functions. Our approach is a nonlinear source inversion. This method requires calculation of waveforms and comparison with observations for many source models. To accelerate these calculations, we use pre-calculated GFs by the reciprocity method in the JIVSM velocity structure model. By comparison with aftershock records, we validated this structure for periods as short as 1.5 s. Target sites are limited to sites close to the fault: KMM005, KMM006, KMMH14, and KMMH16. First, we look for an initial SMGA source model by the grid search method applied to relatively long-period (> 3 s) waveforms and coarse grid of source parameters. Then, we tune that source model to fit observed short-period waveforms with the simplex search method. Necessary physical constraints for the range of the source parameters are applied here. The important point in our inversion scheme is to describe the Kostrov-like slip velocity functions inside each SMGAs by using two triangles. The resulting source model agrees well with other inversion results. We found that the observed westward pulse at KMMH16 is the result of the constructive interference of two pulses from SMGA1 and SMGA2, located in Hinagu fault and southwestern segment of Futagawa fault. Graphical Abstract