In this study, we estimated source process of the 2016 Kumamoto earthquake from strong-motion data by using the multiple-time window linear kinematic waveform inversion method to discuss generation of strong motions and to explain crustal deformation pattern with a seismic source inversion model. A four-segment fault model was assumed based on the aftershock distribution, active fault traces, and interferometric synthetic aperture radar data. Three western segments were set to be northwest-dipping planes, and the most eastern segment under the Aso caldera was examined to be a southeast-dipping plane. The velocity structure models used in this study were estimated by using waveform modeling of moderate earthquakes that occurred in the source region. We applied a two-step approach of the inversions of 20 strong-motion datasets observed by K-NET and KiK-net by using band-pass-filtered strong-motion data at 0.05-0.5 Hz and then at 0.05-1.0 Hz. The rupture area of the fault plane was determined by applying the criterion of Somerville et al. (Seismol Res Lett 70:59-80, 1999) to the inverted slip distribution. From the first-step inversion, the fault length was trimmed from 52 to 44 km, whereas the fault width was kept at 18 km. The trimmed rupture area was not changed in the second-step inversion. The source model obtained from the two-step approach indicated 4.7 × 10 19 Nm of the total moment release and 1.8 m average slip of the entire fault with a rupture area of 792 km 2 . Large slip areas were estimated in the seismogenic zone and in the shallow part corresponding to the surface rupture that occurred during the Mj7.3 mainshock. The areas of the high peak moment rate correlated roughly with those of large slip; however, the moment rate functions near the Earth surface have low peak, bell shape, and long duration. These subfaults with long-duration moment release are expected to cause weak short-period ground motions. We confirmed that the southeast dipping of the most eastern segment is more plausible rather than northwest-dipping from the observed subsidence around the central cones of the Aso volcano.

No abstract

A two-stage scaling relationship of the source parameters for crustal earthquakes in Japan has previously been constructed, in which source parameters obtained from the results of waveform inversion of strong motion data are combined with parameters estimated based on geological and geomorphological surveys. A three-stage scaling relationship was subsequently developed to extend scaling to crustal earthquakes with magnitudes greater than M w 7.4. The effectiveness of these scaling relationships was then examined based on the results of waveform inversion of 18 recent crustal earthquakes (M w 5.4-6.9) that occurred in Japan since the 1995 Hyogo-ken Nanbu earthquake. The 2016 Kumamoto earthquake, with M w 7.0, was one of the largest earthquakes to occur since dense and accurate strong motion observation networks, such as K-NET and KiK-net, were deployed after the 1995 Hyogo-ken Nanbu earthquake. We examined the applicability of the scaling relationships of the source parameters of crustal earthquakes in Japan to the 2016 Kumamoto earthquake. The rupture area and asperity area were determined based on slip distributions obtained from waveform inversion of the 2016 Kumamoto earthquake observations. We found that the relationship between the rupture area and the seismic moment for the 2016 Kumamoto earthquake follows the second-stage scaling within one standard deviation (σ = 0.14). The ratio of the asperity area to the rupture area for the 2016 Kumamoto earthquake is nearly the same as ratios previously obtained for crustal earthquakes. Furthermore, we simulated the ground motions of this earthquake using a characterized source model consisting of strong motion generation areas (SMGAs) based on the empirical Green's function (EGF) method. The locations and areas of the SMGAs were determined through comparison between the synthetic ground motions and observed motions. The sizes of the SMGAs were nearly coincident with the asperities with large slip. The synthetic ground motions obtained using the EGF method agree well with the observed motions in terms of acceleration, velocity, and displacement within the frequency range of 0.3-10 Hz. These findings indicate that the 2016 Kumamoto earthquake is a standard event that follows the scaling relationship of crustal earthquakes in Japan.

No abstract

No abstract

Inter-station cross-correlation functions estimated using continuous ambient noise or microtremor records were used to extract the seismic wave propagation characteristics of the Osaka sedimentary basin, Japan. Temporary continuous observations were conducted at 15 sites in the Osaka basin between 2011 and 2013. The data were analyzed using seismic interferometry. The target period range was 2-8 s. Cross-correlations between all of the possible station pairs were calculated and stacked to produce a year-long data set, and Rayleigh wave signals in the vertical and radial components and Love wave signals in the transverse component were identified from the results. Simulation of interstation Green's functions using the finite difference method was conducted to check the performance of the current three-dimensional velocity structure model. The measured time lag between the observed and theoretical Green's functions was less than 2 s for most station pairs, which is less than the wave period of interest in the target frequency range. Group velocity tomography was applied to group delay times estimated by means of multiple filter analysis. The estimated group velocities for longer periods of 5-8 s exhibited spatial variation within the basin, which is consistent with the bedrock depth distribution; however, the group velocities for shorter periods of 2-3 s were almost constant over the studied area. The waveform and group velocity information obtained by seismic interferometry analysis can be useful for future reconstruction of a three-dimensional velocity structure model in the Osaka basin.

The Hakodate Plain in the southern part of the Oshima Peninsula, Hokkaido, Japan, is a sedimentary basin surrounded by mountains. The vertical displacements caused by active faults along the western margin of the Hakodate Plain increased the basin depth in the western part of the Hakodate Plain. Small-to-large-sized microtremor array surveys were conducted at five sites in the Hakodate Plain to estimate the S-wave velocity structure down to the seismic bedrock for each site to develop a detailed velocity structure model. A new three-dimensional velocity structure model of the Hakodate Plain was developed by integrating the results of microtremor array surveys and other existing geophysical explorations data. This three-dimensional velocity model was modeled as a stack of homogeneous isotropic layers to facilitate its incorporation into the present nation-wide three-dimensional velocity model for ground motion prediction. The bottom depth of the Quaternary sediments is deep along the western margin of the Hakodate Plain. The total thickness of the Quaternary and Neogene sedimentary layers reaches 2.9 km in the western Hakodate Plain. The proposed velocity model was validated by gravity anomaly modeling and ground motion simulation of a moderate-sized inland earthquake. The location of the low-gravity anomaly around the coastline of the Hakodate Bay was improved using the new model. The numerical ground motion simulation using FDM also demonstrated that the amplification and long duration observed in the western part of the Hakodate Plain were reproduced effectively using this velocity model. The spatial variation in long-period ground motion amplifications (period > 1 s) is discussed based on numerical simulations utilizing our three-dimensional velocity model. The highest amplifications at periods of 4 and 5 s were expected in the southwestern Hakodate Plain. The amplification at a period of 3 s was relatively high near the western margin of the plain. Conversely, the spatial characteristics below 2 s were quite complex due to interference of the seismic wavefield inside the basin structure. Variation due to the source location was also relatively high in the shorter period range. Graphical abstract

Ground motions near the source area of the mainshock of the 2013 Lushan earthquake (Mw 6.6) in the Sichuan Province in China were reproduced using the characterized source model and the empirical Green’s function method (EGFM). The best-fit characterized source model consists of one strong motion generation area (SMGA) and a background area. The synthesized ground motions of the characterized source model were in fairly good agreement with the observed ones in the frequency range from 0.5 to 30.0 Hz at eleven strong-motion stations. The relationships between the SMGA and the seismic moment, and those between the flat amplitude of acceleration source spectrum in the short-period and the seismic moment, for the 2013 Lushan earthquake (Mw 6.6), almost followed the empirical scaling relationships of inner fault parameters so far developed for the crustal earthquakes. The reason for the largest PGA (>1g) in the strong-motion observation history of China recorded at the 051BXD strong-motion station was investigated from the source and site effects. Through comparing the synthesized ground motions at the station for three different rupture-starting-point and comparing the earthquake H/V spectral ratios between the mainshock and the EGF event, neither the significant directivity effect nor the strong nonlinear effect was observed on the ground motions at the 051BXD. The predominant period around 0.1 s of the Fourier spectra of both the mainshock and the EGF event suggested that the site effect should be responsible for the large record.