Search

Sliding of natural and artificial slopes generally occurs during or following strong earthquakes. Such sliding is greatly affected by a combination of geological conditions and earthquake loading. Earthquake-induced landslides often cause more damage to infrastructure and human lives than the earthquake itself. Pseudo-static analysis is widely implemented as one of several design methods used in engineering practice to assess the seismic stability of natural and artificial slopes. However, the most important issue of pseudo-static analysis is to select the most appropriate method for measuring seismic coefficient. In order to investigate this, back analysis was conducted for surficial slides subjected to strong ground motion during the 2004 Mid-Niigata prefecture earthquake in Japan. This paper surveyed the stochastic properties of earthquake-induced surficial slides and clearly showed that the obtained results were applicable to back analysis of shear strength and seismic coefficient. In back analysis, soil properties such as soil strength and density and sliding depth were assumed as random variables owing to their uncertainties. Seismic coefficient is also assumed to be a random variable and varies with distance from the epicenter fault line. The analysis of 4504 recorded surficial slides clearly shows a unique relationship of landslide occurrence ratio with slope angle and distance from the epicenter fault line. This study’s results enhance the calculation of the shear strength of weathered soil covering slopes and the horizontal seismic coefficient through back analysis procedure. By considering possible stochastic properties of variables, some case studies were implemented in the back analysis. [ABSTRACT FROM AUTHOR]

Dense strong motion observation networks provided us with valuable data for studying strong motion generation from large earthquakes. From kinematic waveform inversion of seismic data, the slip distribution on the fault surface of large earthquakes is known to be spatially heterogeneous. Because heterogeneities in the slip and stress drop distributions control the generation of near-source ground motion, it is important to characterize these heterogeneities for past earthquakes in constructing a source model for reliable prediction of strong ground motion. The stress changes during large earthquakes on the faults recently occurring in Japan are estimated from the detailed slip models obtained by the kinematic waveform inversion. The stress drops on and off asperities are summarized on the basis of the stress change distributions obtained here. In this paper, we define the asperity to be a rectangular area whose slip is 1.5 or more times larger than the average slip over the fault according to the previous study for inland crustal earthquakes. The average static stress drops on the asperities of the earthquakes studied here are in the range 6-23 MPa, whereas those off the asperities are below 3 MPa. We compiled the stress drop on the asperities together with a data set from previous studies of other inland earthquakes in Japan and elsewhere. The static stress drop on the asperity depends on its depth, and we obtained an empirical relationship between the static stress drop and the asperity's depth. Moreover, surface-breaking asperities seemed to have smaller stress drops than buried asperities. Simple ground motion simulations using the characterized asperity source models reveal that deep asperities generate larger ground motion than shallow asperities, because of the different stress drops of the asperities. These characteristics can be used for advanced source modeling in strong ground motion prediction for inland crustal earthquakes. [ABSTRACT FROM AUTHOR]

A large number of debris flows occurred in the Wenchuan earthquake zone after the 12 May 2008 earthquake. The risks posed by these debris flows were rather high. An appropriate model is required to predict the possible runout distance and impacted area. This paper describes a study on the runout characteristics of the debris flows that occurred in the Wenchuan earthquake zone over the past four years. A total of 120 debris flows are analyzed. Separate multivariate regression models are established for the runout distances of hill-slope debris flows and channelized debris flows. The control variables include type of debris flow, debris flow volume, and elevation difference. Comparison of the debris flows occurring before and after the earthquake shows that the runout distance increased after the earthquake due to sufficient material supply and increased mobility of the source materials. In addition, the runout distances of annual debris flow events in 2008, 2010 and 2011 are analyzed and compared. There is a tendency that the runout distance decreases over time due to the decreasing source material volume and possible changes of debris flow type. Comparison between the debris flows in the earthquake zone and the debris flows in Swiss Alps, Canada, Austria, and Japan shows that the former have a smaller mobility. [ABSTRACT FROM AUTHOR]

Aki (1969) first modeled coda waves of a local earthquake as a superposition of scattered surface waves. This paper attempts to clarify the constituents of surface-wave coda at long periods at very long lapse times. For a large earthquake of magnitude 7 or larger, vertical component oscillation in periods from 90 to 180 s persists for more than 20 hours from the earthquake origin time. Although the early portion of the coda envelope is successfully modeled by assuming incoherent scattered Rayleigh waves by heterogeneities distributed all over the Earth, the later potion of the observed coda envelope (roughly later than 35,000 s) has systematically larger amplitude than theoretical prediction. To clarify the cause of this discrepancy, we studied the constituents of vertical-component seismograms of three large earthquakes recorded by the F-net in Japan using the f- k power spectral analysis. We found that the direct and scattered fundamental-mode Rayleigh waves of velocity about 3.7 km/s are dominant in the earlier part of each envelope. It justifies the use of a scattering model of the fundamental Rayleigh waves for synthesizing the envelope. At lapse times later than 20,000 s–35,000 s, higher modes with phase velocities around 20 km/s become dominant. The transition time to the dominance of higher modes is found to become earlier for a deeper focus earthquake. The small coda attenuation factor from (1.90±0.23) × 10−3 to (2.38±0.32) × 10−3 estimated from later coda envelopes recorded at IRIS stations distributed worldwide also agrees with the attenuation factor of spheroidal modes according to PREM. We may interpret that higher mode waves are uniformly distributed at large lapse time due to large velocity dispersion and/or scattering and they dominate over the fundamental mode waves because of smaller attenuation in the lower mantle. The coda attenuation measurement proposed by Aki is found to be useful even for long periods and at very large lapse times. [ABSTRACT FROM AUTHOR]

We propose a data mining method to extract physical phenomena from the ground vibration velocity fluctuation (GVF) which is a daily observed weak signal and has no specific waveform. Our aim is to detect any traces of the return of ground motion to its normal state after earthquakes, and the evolution of ground motion towards earthquakes in the GVF. We focus on the GVF recorded near the epicenter of the magnitude 9 Great East Japan Earthquake (GEJE). Although the GVF can be discussed in the framework of the master equation which represents the thermodynamics of a stochastically fluctuating non-equilibrium system, the relation between the thermodynamic parameters and seismology-related phenomena is hardly known. So, we introduce cellular automata (CA) to analogically interpret the physical phenomena of interest as thermodynamic parameters, and show the thermodynamic equivalence between GVF and CA. These CA-derived thermodynamic parameters are then expected to deliver physical phenomena similar to those of the target when evaluated on the GVF data. In the demonstrated example of data mining, the stress relaxation signal immediately after GEJE is discovered from the recorded GVF data by deriving the thermodynamic parameters representing the stress relaxation of CA and evaluating these thermodynamic parameters for the recorded GVF data. Furthermore, the elastic rigidity and viscosity of the underground structures are estimated from the stress relaxation signal. The novelty of this study lies in drawing attention to the GVF, and combining standard techniques of CA with GVF through the use of the master equation. [ABSTRACT FROM AUTHOR]

By interpreting the 2D/3D seismic survey data acquired in the surrounding ocean areas of the Northeast (NE) Japan Arc, we clarified the detailed geological structure and demonstrated that the basic structure in the hanging-wall plate of the subduction system consists of many structural blocks (segments) separated by NW–SE trending large transcurrent faults (strike-slip faults). This structural configuration showed a close relationship with the distribution of foreshocks, mainshock, and aftershocks, coseismic slip models of the 2011 M9.0 Tohoku-Oki megathrust earthquake, coseismic slip area of M-7 class earthquakes, quasi-static slip rates, back slip rate, and seismic tomography images. In addition, the coseismic slip models revealed that the trenchward forearc of the structural blocks between the Offshore Hidaka tectonic line and the Honjo-Sendai tectonic line fitted well with the coseismic slip area of the 2011 Tohoku-Oki earthquake. These findings suggest that the structural blocks bounded by these two tectonic lines slipped rapidly trenchward when the mainshock occurred. The M7 earthquakes were also concentrated along these two tectonic lines, thereby suggesting a close relationship between seismic activity and the inherited geological structure of the overriding plate in the NE Japan forearc. [ABSTRACT FROM AUTHOR]

An approach to organization of a system for automatic earthquake prediction is developed. This approach is not bound to any particular type of data, and all data on seismotectonic processes are displayed using the spatial and space-time fields. The prediction field is learned from the feature fields and the sample of predicted earthquakes. The learning method proposed here is the least alarm method. The results of modeling of earthquake prediction for the regions of Japan and Mediterranean are presented. Modeling implied the analysis of prediction fields reflecting stationary and dynamic properties of the seismic process. The results of modeling demonstrate the efficiency of application of the least alarm method to earthquake prediction. [ABSTRACT FROM AUTHOR]

The M 7.3 Kumamoto earthquake that occurred at 1:25 JST on April 16, 2016, not only triggered aftershocks in the vicinity of the epicenter, but also triggered earthquakes that were 50-100 km away from the epicenter of the main shock. The active seismicity can be divided into three regions: (1) the vicinity of the main faults, (2) the northern region of Aso volcano (50 km northeast of the mainshock epicenter), and (3) the regions around three volcanoes, Yufu, Tsurumi, and Garan (100 km northeast of the mainshock epicenter). Notably, the zones between these regions are distinctively seismically inactive. The electric resistivity structure estimated from one-dimensional analysis of the 247 broadband (0.005-3000 s) magnetotelluric and telluric observation sites clearly shows that the earthquakes occurred in resistive regions adjacent to conductive zones or resistive-conductive transition zones. In contrast, seismicity is quite low in electrically conductive zones, which are interpreted as regions of connected fluids. We suggest that the series of the earthquakes was induced by a local accumulated stress and/or fluid supply from conductive zones. Because the relationship between the earthquakes and the resistivity structure is consistent with previous studies, seismic hazard assessment generally can be improved by taking into account the resistivity structure. Following on from the 2016 Kumamoto earthquake series, we suggest that there are two zones that have a relatively high potential of earthquake generation along the western extension of the MTL. [ABSTRACT FROM AUTHOR]

We investigated the development of a distinct later phase observed at stations near the Japan Trench associated with shallow, outer-rise earthquakes off the coast of Sanriku, northern Japan based on the analysis of three-component broadband seismograms and FDM simulations of seismic wave propagation using a heterogeneous structural model of the Japan Trench subduction zone. Snapshots of seismic wave propagation obtained through these simulations clearly demonstrate the complicated seismic wavefield constructed by a coupling of the ocean acoustic waves and the Rayleigh waves propagating within seawater and below the sea bottom by multiple reflections associated with shallow subduction zone earthquakes. We demonstrated that the conversion to the Rayleigh wave from the coupled ocean acoustic waves and the Rayleigh wave as they propagate upward along the slope of seafloor near the coast is the primary cause of the arrival of the distinct later phase at the station near the coast. Through a sequence of simulations using different structural models of the Japan Trench subduction zone, we determined that the thick layer of seawater along the trench and the suddenly rising sea bottom onshore of the Japanese island are the major causes of the distinct later phase. The results of the present study indicate that for realistic modeling of seismic wave propagation from the subduction zone earthquakes, a high-resolution bathymetry model is very crucial, although most current simulations do not include a water column in their simulation models. [ABSTRACT FROM AUTHOR]

The tectonic and seismological aspects of the Great Japan Earthquake, which occurred on March 11, 2011 ( M = 9.0), at the Pacific margin of the northeastern part of Honshu Island, are discussed. The structure and seismotectonic data, seismicity, and the reccurence rate of the great ( M ≥ 7.6) earthquakes throughout history and in modern times are represented. It is shown that the reccurence rate of the great events is about 40 years, and that of megaearthquakes is 1000 years or more. A seismic gap of about 800 km in length is found in the region under study, located to the south of latitude 39° N and full of aftershocks to the megaearthquake of March 11, 2011. This event is probably connected with the deep thrust along the Benioff zone and its structural front (Oyashio nappe at the middle Pacific continental slope). The aftershock sequences of this megaearthquake and the Sumatra-Andaman (2004) megaearthquake are compared. It is found that several of their key characteristics (the number of aftershocks, the magnitude of the strongest aftershock, and the time of its occurrence) for 25 days are comparable for both cases with a significant difference in the energies of aftershock processes. A probable scenario for the origination of a repeated shock with M ∼ 8.0 in the Japan Trench is discussed. [ABSTRACT FROM AUTHOR]

The dense recordings of the K-NET and KiK-net nationwide strong motion network of 1,189 accelerometers show clearly the radiation and propagation properties of the strong ground motions associated with the 2011 off-the-Pacific Coast-of-Tohoku, Japan (Mw = 9.0) earthquake. The snapshots of seismic wave propagation reveal strong ground motions from this earthquake that originate from three large slips; the first two slips occurred over the plate interface of off-Miyagi at the southwest and the east of the hypocenter, and the third one just beneath the northern end of Ibaraki over the plate interface or in the crust. Such multiple shocks of this event caused large accelerations (maximum 1-2 G) and prolonged ground shaking lasting several minutes with dominant high-frequency ( T < 1 s) signals over the entire area of northern Japan. On the other hand, ground motions of relatively longer-period band ( T = 1-2 s), which caused significant damage to wooden-frame houses, were about 1/2-1/3 of those observed near the source area of the destructive 1995 Kobe, Japan (M = 7.3) earthquake. Also, the long-period ( T = 6-8 s) ground motion in the Kanto (Tokyo) sedimentary basin was at an almost comparable level of those observed during the recent Mw = 7 inland earthquakes, but not as large as that from the former M = 8 earthquakes. Therefore, the impact of the strong ground motion from the present M = 9 earthquake was not as large as expected from the previously M = 7-8 earthquakes and caused strong motion damage only to short-scale construction and according to instruments inside the buildings, both have a shorter ( T < 1 s) natural period. [ABSTRACT FROM AUTHOR]

The article discusses a study of co-seismic energy release and the relevant region of the magnitude 9 earthquake in Tohoku, Japan on March 11, 2011. The rebound of a continental plate against a subducting oceanic plate at the Japan Trench is the cause of the earthquake. The co-seismic displacements were fitted using the Bi-Harmonic Interpolation Spline method. The adjacent elastic unloading body reportedly provides the energy dissipation associated with earthquake.

If aseismic slip occurs on a fault or its deeper extension, both seismicity and crustal deformation around the source would be affected. Anomalous phenomena of this kind are revealed from earthquake occurrence data and geodetic records during a period of 10 years leading up to the March 2005 M 7.0 earthquake west of Fukuoka that occurred off the northern coast of Kyushu, Japan. Seismicity rate anomalies (quiescence and activation) took place relative to the rates expected by the ETAS model in a number of seismic zones in and around the Kyushu District. The seismic zone of the relative quiescence and activation consistently corresponds to the zone of the negative and positive ΔCFS (Coulomb failure stress change), respectively, assuming the precursory aseismic slips on the M 7.0 source fault. In addition, we consider the time series of geodetic baseline distances between permanent GPS stations in the Kyushu District for the same period, which also supports the possible precursory slips rather than the known slow slips beneath the Bungo Straight, off the eastern coast of Kyushu. [ABSTRACT FROM AUTHOR]

We construct a single hazard function from multiple predictive parameters independently developed for moderate earthquakes in Kanto, Japan, during a learning period from 1990 to 1999, and applied to a testing period from 2000 to 2005. Here, we consider as predictive parameters the a and b values of the Gutenberg–Richter relation, the ν value (change in b value), and the Every Earthquake a Precursor According to Scale (EEPAS) model rate. To study the correlations among the parameters, we prepare two groups of space–time coordinate sets for assessment, namely the background and conditional groups selected from the learning period. The background group contains ten thousand sets of coordinates randomly selected from the space–time volume of our study. The conditional group contains 33 sets of space–time coordinates corresponding to the epicenters of the target earthquakes ( M ≥ 5.0) just before their times of occurrence. Each parameter for the background group is transformed so that its distribution conforms to the standard Normal function. The mean and variance of the conditional distribution is then estimated after applying the same transformation to the conditional group. Using the means and variances of b values, ν values and EEPAS rates and the correlation matrices in the background and conditional distributions, we construct a combined hazard function following the procedure developed for normally distributed parameters. The information gain per event ( IGpe) of the new hazard function is 0.26 and 0.3 units larger than that of the EEPAS rate for the learning and testing period, respectively. The R-test confirms the statistical significance of the difference in the IGpe value for the testing period. [ABSTRACT FROM AUTHOR]

Strong ground motions recorded in central Tokyo during the 1944 Tonankai Mw8.1 earthquake occurring in the Nankai Trough demonstrate significant developments of very large (>10 cm) and prolonged (>10 min) shaking of long-period (T > 10–12 s) ground motions in the basin of Tokyo located over 400 km from the epicenter. In order to understand the process by which such long-period ground motions developed in central Tokyo and to mitigate possible future disasters arising from large earthquakes in the Nankai Trough, we analyzed waveform data from a dense nation wide strong-motion network (K-NET and KiK-net) deployed across Japan for the recent SE Off-Kii Peninsula (Mw 7.4) earthquake of 5 September 2004 that occurred in the Nankai Trough. The observational data and a corresponding computer simulation for the earthquake clearly demonstrate that such long-period ground motion is primarily developed as the wave propagating along the Nankai Trough due to the amplification and directional guidance of long-period surface waves within a thick sedimentary layer overlaid upon the shallowly descending Philippine Sea Plate below the Japanese Island. Then the significant resonance of the seismic waves within the thick cover of sedimentary rocks of the Kanto Basin developed large and prolonged long-period motions in the center of Tokyo. The simulation results and observed seismograms are in good agreement in terms of the main features of the long-period ground motions. Accordingly, we consider that the simulation model is capable of predicting the long-period ground motions that are expected to occur during future Nankai Trough M 8 earthquakes. [ABSTRACT FROM AUTHOR]

In this study, the focal mechanism solutions and spatial distributions of 100 intermediate-depth earthquakes are examined and discussed. These earthquakes occurred from 1985 to 1993 mostly in the upper and lower planes of the double seismic zone beneath the southeastern part of Hokkaido. The focal-mechanism solutions along with the spatial distribution of their epicenters as a function of depth and the minimum- and maximum-compressive stress directions confirm the following essential features; A double-planed deep seismic zone between the depths from 60 to 120 km has been clearly revealed in this region and the number of earthquakes located in the upper plane are less than those for the lower one. Low-angle thrust fault mechanism is dominant to a depth of about 55–60 km in the zone. At depths deeper than about 60 km focal mechanism solutions are mostly characterized by down-dip extension type or down-dip compression type in both lower and upper planes of the zone. Most of the earthquakes deeper than 60 km have dip-slip components and some of the solutions have significant strike-slip components, as well. Although the down-dip extensional-type events are dominant, seven events that are located on the lower plane are constrained as down-dip compression mechanism. The T axes of earthquake focal mechanisms on the lower plane show extension at about N40°E and N35°W, and the dip-angle between 25° and 70°, almost coincide with the dip of the subducting Pacific Plate in this region. We also investigate the statistical features of the stress directions of the earthquakes occurring both on the lower and upper plane after converting all the events onto the plane having the same strike and dip angle as the subducting lithosphere in the region. Seventy-two percent of the minimum compressive stress directions are oriented mainly in the NS or NNE- SSW directions (within 45° of the dip direction) at the plate surface of the lower plane and are almost parallel to the dip of the local slab in this region. The number of earthquakes located in the upper plane is less than that for the lower one. Those which occurred on the upper plane show various types of mechanism solutions. Weak systematic orientation is found in the direction distribution of stress in the upper plane. [ABSTRACT FROM AUTHOR]

The Pacific plate converges with northeastern Eurasia at a rate of 8-9 m per century along the Kamchatka, Kuril and Japan trenches[SUP1]. Along the southern Kuril trench, which faces the Japanese island of Hokkaido, this fast subduction has recurrently generated earthquakes with magnitudes of up to ∼8 over the past two centuries[SUP2-6]. These historical events, on rupture segments 100-200 km long, have been considered characteristic of Hokkai-do's plate-boundary earthquakes[SUP7,8]. But here we use deposits of prehistoric tsunamis to infer the infrequent occurrence of larger earthquakes generated from longer ruptures. Many of these tsunami deposits form sheets of sand that extend kilometres inland from the deposits of historical tsunamis. Stratigraphic series of extensive sand sheets, intercalated with dated volcanic-ash layers, show that such unusually large tsunamis occurred about every 500 years on average over the past 2,000-7,000 years, most recently ∼350 years ago. Numerical simulations of these tsunamis are best explained by earthquakes that individually rupture multiple segments along the southern Kuril trench. We infer that such multi-segment earthquakes persistently recur among a larger number of single-segment events. [ABSTRACT FROM AUTHOR]

Non-volcanic tremor is a weak, extended duration seismic signal observed episodically on some major faults, often in conjunction with slow slip events. Such tremor may hold the key to understanding fundamental processes at the deep roots of faults, and could signal times of accelerated slip and hence increased seismic hazard. The mechanism underlying the generation of tremor and its relationship to aseismic slip are, however, as yet unresolved. Here we demonstrate that tremor beneath Shikoku, Japan, can be explained as a swarm of small, low-frequency earthquakes, each of which occurs as shear faulting on the subduction-zone plate interface. This suggests that tremor and slow slip are different manifestations of a single process. [ABSTRACT FROM AUTHOR]