Search

Slow earthquakes are episodic slow fault slips. They form a fundamental component of interplate deformation processes, along with fast, regular earthquakes. Recent seismological and geodetic observations have revealed detailed slow earthquake activity along the Japan Trench—the subduction zone where the March 11, 2011, moment magnitude (Mw) 9.0 Tohoku-Oki earthquake occurred. In this paper, we review observational, experimental, and simulation studies on slow earthquakes along the Japan Trench and their research history. By compiling the observations of slow earthquakes (e.g., tectonic tremors, very-low-frequency earthquakes, and slow slip events) and related fault slip phenomena (e.g., small repeating earthquakes, earthquake swarms, and foreshocks of large interplate earthquakes), we present an integrated slow earthquake distribution along the Japan Trench. Slow and megathrust earthquakes are spatially complementary in distribution, and slow earthquakes sometimes trigger fast earthquakes in their vicinities. An approximately 200-km-long along-strike gap of seismic slow earthquakes (i.e., tectonic tremors and very-low-frequency earthquakes) corresponds with the huge interplate locked zone of the central Japan Trench. The Mw 9.0 Tohoku-Oki earthquake ruptured this locked zone, but the rupture terminated without propagating deep into the slow-earthquake-genic regions in the northern and southern Japan Trench. Slow earthquakes are involved in both the rupture initiation and termination processes of megathrust earthquakes in the Japan Trench. We then compared the integrated slow earthquake distribution with the crustal structure of the Japan Trench (e.g., interplate sedimentary units, subducting seamounts, petit-spot volcanoes, horst and graben structures, residual gravity, seismic velocity structure, and plate boundary reflection intensity) and described the geological environment of the slow-earthquake-genic regions (e.g., water sources, pressure–temperature conditions, and metamorphism). The integrated slow earthquake distribution enabled us to comprehensively discuss the role of slow earthquakes in the occurrence process of the Tohoku-Oki earthquake. The correspondences of the slow earthquake distribution with the crustal structure and geological environment provide insights into the slow-earthquake-genesis in the Japan Trench and imply that highly overpressured fluids are key to understanding the complex slow earthquake distribution. Furthermore, we propose that detailed monitoring of slow earthquake activity can improve the forecasts of interplate seismicity along the Japan Trench. [ABSTRACT FROM AUTHOR]

A modified k-nearest-neighbors method is introduced that provides an efficient nonlinear estimate of the intensity of a point process (field) based on locations of events. The method is applied to perform a detailed statistical analysis of the spatial structure of the seismic intensity field (in the context of this paper, the term "intensity" refers to the intensity of the seismic field, and should not be mistaken for the intensity of ground motion caused by an earthquake). The proposed method requires neither a preliminarily delineated area nor a normalization procedure for the estimates. In contrast to many interpolation methods, estimates based on the proposed method are statistically justified. The "uncertainty relation" between the spatial smoothing effect and random errors is established for the method in an explicit form. A procedure for choosing the number of nearest neighbors k* that controls the effective radius of smoothing is suggested. The proposed method is applied to analyze the seismic intensity field in two seismogenic regions surrounding the Kuril Islands and Japan during the period from 1904 to 2011. Spots of increased seismic activity in these regions are localized by the new method, and some quantitative statistical characteristics of these spots are determined and discussed. [ABSTRACT FROM AUTHOR]

It is an acquired result that, in order to enable realistic earthquake hazard assessment and reliable estimation of the ground motion response to an earthquake, three-dimensional velocity models have to be considered. In this paper we present a new analytical procedure for determining the seismic wavefield in a 3D anelastic model based on the combination of the ray theory with the modal summation method. The proposed procedure has been validated by considering the three-dimensional model of the Kanto basin (Japan) available in the literature. Three-dimensional simulations were performed for the 1990 Odawara earthquake. The results obtained were compared both with recorded signals and with simulations available in the literature for several stations located within the Kanto basin. Besides the advantage of being a useful tool for assessment of seismic hazard and seismic risk reduction, it is characterized by high efficiency, in fact, once the study region is identified and the 3D model is constructed, the computation, at each station, of the three components of the synthetic signal (displacement, velocity, and acceleration) takes less than 3 h on a 2 GHz CPU. [ABSTRACT FROM AUTHOR]

The Japan Meteorological Agency (JMA) unified seismic catalog has been widely used for research and disaster prevention purposes for more than 20 years. Since the introduction in April 2016 of an improved method of automatic hypocenter determinations (PF method), the number of detected earthquakes has almost doubled due to a decrease in the completeness magnitude around the Tohoku region, where seismicity has been very active in the aftermath of the 2011 Tohoku earthquake. Automatically processed hypocenters of small events, accepted without manual modification, now make up approximately 70% of new events in the JMA unified catalog. In this paper, we show that the introduction of automated processing did not systematically bias the quality of the JMA unified catalog. Approximately 90% of automatically processed hypocenters were less than 1 km from their manually reviewed locations in inland and shallow areas. We also considered the use of automated event characterization in real-time monitoring of earthquake sequences using the example of the April 2016 Kumamoto earthquake sequence, when the PF method could have supplied the catalog with about 70,000 events in real time over the course of 2 months. We show that the PF method is capable of monitoring the migration or expansion of the hypocentral distribution and can support statistical analyses such as variations of the b-value distribution. Further improvements in automatic hypocenter determination will contribute to a better understanding of seismicity as well as rapid risk assessment, especially in cases of swarms and aftershocks. [ABSTRACT FROM AUTHOR]

On 6 September 2018, a Mj 6.7 (Mw 6.6) earthquake with a focal depth of approximately 37 km and maximum seismic intensity of 7.0 on the Japan Meteorological Agency (JMA) scale (corresponding to approximately X on the Modified Mercalli Intensity (MMI) scale) struck the central and eastern Iburi regions of Hokkaido, Northern Japan, 1 day after the Typhoon Jebi passed through the region. Thousands of landslides were triggered and significant losses resulted from the earthquake sequence. A detailed landslide inventory map, including 5625 coseismic landslides, was delineated on the basis of 3307 published landslide sites. Most of the coseismic landslides are translational landslides of small to medium scale with high mobility and long run-out distance. Thirty-six people were killed by the landslides despite the afflicted area being sparsely populated. It is found that all the 5625 landslides spread in an elliptic area extending NNW/SSE, running approximately parallel to the strike of (active) faults in this region. The preferred aspect of the landslide-affected area is southerly, running nearly perpendicular to the NNW/SSE striking (active) faults. Most coseismic landslides are distributed in regions with seismic intensity of 7.0 to 8.0 (MMI Scale), with peak ground acceleration (PGA) of 0.4 to 0.7 g. Most of the coseismic landslides occurred at the elevation between 100 and 250 m, and the slope angle between 15° and 35°. Miocene sedimentary rock is the predominant bedrock type identified in the landslide area. Slope failures were triggered in stratified pyroclastic fall deposits, in the combination of strong seismic ground motion and intense antecedent precipitation. Sliding zone liquefaction phenomena were confirmed in the field investigation. The relationship between the old landslides (slope failures occurred prior to the Iburi earthquake) and the coseismic landslides is also discussed in this study. [ABSTRACT FROM AUTHOR]

A probabilistic seismic hazard assessment was performed for the Japanese islands and surrounding areas. Seismic hazard parameters characteristic of the seismic history of the regions were obtained. The probability of occurrence of a large M ≥ 7 earthquake within a 10- and 50-year period was also calculated. Regions of very high levels of hazard occur where the Pacific, Phillipine and Eurasian Plates meet. High probabilities of occurrence of a large M ≥ 7 earthquake within a 10- and 50-year period occur within the region where the Pacific Plate subducts with the Eurasian Plate. [ABSTRACT FROM AUTHOR]

A seismic waveform inversion algorithm is proposed for the estimation of elastic soil properties using low amplitude, downhole array recordings. Based on a global optimization scheme in the wavelet domain, complemented by a local least-square's fit operator in the frequency domain, the hybrid scheme can efficiently identify the optimal solution vicinity in the stochastic search space, whereas the best-fit model detection is substantially accelerated through the local deterministic inversion. Results presented for selected aftershocks of the M w 7.0 Sanriku-Minami earthquake in Japan, recorded by the Kik-Net Strong Motion Network, illustrate robustness of the impedance structure estimation. By contrast, the attenuation structure is shown to be sensitive to the frequency content of seismic input data, attributed to the deterministic description of the continuum in the forward model that cannot simulate late arrivals of multiple-scattered energy. Sensitivity analyses illustrate that for the same forward model, results can be substantially different based on the definition of the objective function. It is concluded that even for engineering purposes, inversion should aim to decouple intrinsic and scattering attenuation mechanisms. [ABSTRACT FROM AUTHOR]