Modeling of strong motion generation areas of the 2011 Tohoku,Japan earthquake using modified semi-empirical technique

Modification in the semi-empirical technique for the simulation of strong ground motion has been introduced to incorporate the strong motion generation areas (SMGA) in the modeled rupture plane. Strong motion generation areas identified within the rupture plane of the Tohoku earthquake of March 11, 2011 (M _w = 9.0), have been modeled using this modified technique. Two different source models having four and five SMGAs, respectively, are considered for modeling purpose. Strong motion records using modified semi-empirical technique have been simulated at two near-field stations located at epicentral distance of 137 and 140 km, respectively, using two different source models. Comparison of the observed and simulated acceleration waveforms is made in terms of root mean square error (RMSE) at both stations. Minimum root mean square error of the waveform comparison has been obtained at both the stations for source model having five SMGAs. Simulations from same rupture model have been made at other four stations lying at epicentral distance between 154 and 249 km. Comparison of observed and simulated records has been made in terms of RMSE in acceleration records, velocity records and response spectra at each six station. Simulations have been made at six other stations to obtain distribution of peak ground acceleration and peak ground velocity with hypocentral distance. Peak ground acceleration and velocity from simulated and observed records are compared at twelve stations surrounding the source of Tohoku earthquake. Comparison of waveforms and parameters extracted from observed and simulated strong motion records confirms the efficacy of the developed modified technique to model earthquake characterized by SMGAs.     

Near-field and far-field simulation of accelerograms of Sikkim earthquake of September 18, 2011 using modified semi-empirical approach

The semi-empirical approach for modeling of strong ground motion given by Midorikawa (Tectonophysics 218:287?C295, 1993) has been modified in the present paper for component wise simulation of strong ground motion. The modified approach uses seismic moment in place of attenuation relation for scaling of acceleration envelope. Various strong motion properties like directivity effect and dependence of peak ground acceleration with respect to surface projection of source model have been studied in detail in the present work. Recently, Sikkim earthquake of magnitude 6.9 (M _w ) that occurred on September 18, 2011 has been recorded at various near-field and far-field strong motion stations. The modified semi-empirical technique has been used to confirm the location and parameters of rupture responsible for this earthquake. Strong motion record obtained from the iterative modeling of the rupture plane has been compared with available strong motion records from near as well as far-field stations in terms of root mean square error between observed and simulated records. Several possibilities of nucleation point, rupture velocity, and dip of rupture plane have been considered in the present work and records have been simulated at near-field stations. Final selection of model parameters is based on root mean square error of waveform comparison. Final model confirms southward propagating rupture. Simulations at three near-field and twelve far-field stations have been made using final model. Comparison of simulated and observed record has been made in terms of peak ground acceleration and response spectra at 5?% damping. Comparison of simulated and observed record suggests that the method is capable of simulating record which bears realistic appearance in terms of shape and strong motion parameters. Present work shows that this technique gives records which matches in a wide frequency range for Sikkim earthquake and that too from simple and easily accessible parameters of the rupture plane.     

Modeling of strong motion generation area of the Uttarkashi earthquake using modified semiempirical approach

The semiempirical approach based on envelope summation method given by Midorikawa (Tectonophysics 218:287–295, 1993) has been modified in this paper for modeling of strong motion generation areas (SMGAs). Horizontal components of strong ground motion have been simulated using modifications in the semiempirical approach given by Joshi et al. (Nat Hazard 71:587–609, 2014). Various modifications in the technique account for finite rupture source, layering of earth, componentwise division of energy and frequency-dependent radiation pattern. In this paper, SMGAs of the Uttarkashi earthquake have been modeled. Two different isolated wave packets in the recorded accelerogram have been identified from recorded ground motion, which accounts for two different SMGAs in the entire rupture plane. The approximate locations of SMGAs within the rupture plane were estimated using spatio-temporal variation of 77 aftershocks. Source parameters of each SMGA were calculated from theoretical and observed source displacement spectra computed from two different wave packets in the record. The final model of rupture plane responsible for the Uttarkashi earthquake consists of two SMGAs, and the same has been used to simulate horizontal components of acceleration records at different station using modified semiempirical technique. Comparison of the observed and simulated acceleration records in terms of root mean square error confirms the suitability of the final source model for the Uttarkashi earthquake.     

Emergence of the semi-empirical technique of strong ground motion simulation: A review

High frequency ground motion simulation techniques are powerful tools for designing earthquake resistant structures in seismically active regimes. Simulation techniques also provide the synthetic strong ground motion in the regions where actual records are not available (Kumar et al. 2015).These techniques require several parameters of earthquake and other seismic information proceeding to the simulation. Practically estimation of parameters is a tough task, particularly in a region with limited information. This demands a simulation technique based on the easily estimated parameters for a new site. The purposes of this paper are to briefly review existing simulation techniques and to discuss in detail the new, simple and effective semi-empirical technique (Midorikawa 1993) of strong motion simulation.     

First-order rupture features of the 2011 MW 9.0 Tohoku (Japan) earthquake from surface waves

This study analyzed the rupture directivity of the 2011 Tohoku earthquake by using 100-s Rayleigh-wave travel-times, influenced by the finite source, to derive the fault parameters of the earthquake. The results demonstrated that the earthquake exhibited a slow rupture propagation with a rupture velocity of approximately 1.5–2.0 km/s and asymmetric bilateral faulting. The two rupture directions were N60°E and N127°E, with rupture lengths of approximately 276 km and 231 km, respectively. The rupture toward N60°E had a source duration of approximately 183 s, longer than that toward N127°E (approximately 156 s). Overall, the entire source duration of the earthquake faulting lasted approximately 183 s. Regarding historical seismicity in eastern Japan, the 2011 Tohoku earthquake not only ruptured a locked area in which large earthquakes have rarely occurred, but also ruptured the source regions of several historical earthquakes. With the exception of its slow rupture velocity and generation of a tsunami, the rupture features of the 2011 Tohoku earthquake were inconsistent with those of typical tsunami earthquakes.     

Seamount subduction and rupture characteristics of the March 11, 2011, Tohoku earthquake

We suggest that the spatial location of the 2011 Tohoku earthquake rupture and slip distribution on it was strongly influenced by the subduction of seamount chains. Subduction of seamounts across the Japan trench caused weak coupling on the plate interface which acted as barriers to the 2011 Tohoku earthquake rupture and thus delimited it.     

Destructive near-fault strong ground motion from the 2016 Kumamoto prefecture,Japan, M7.3 earthquake

The sequence of the 2016 Kumamoto, Japan, earthquake, which included an initial M6.5 foreshock on April 14, followed by a larger M7.3 mainshock on April 16, and subsequently occurred high aftershock activity, caused significant damage in Kumamoto and neighboring regions. The near-field strong motion record by strong motion network (K-NET and KiK-net) and the intensity meter network demonstrated clearly the characteristics of the strong ground motion developed by the shallow (H = 12 km), inland earthquake comprising short-time duration (<15–20 s) but large (>1G) ground accelerations. The velocity response spectra of the near-fault motion at Mashiki and Nishihara showed large levels (>300–550 cm/s) in the short-period range (T = 1–2 s), several times larger than that of the near-field record of the destructive 1995 Kobe earthquake (M7.3) and that of the 2004 Mid-Niigata earthquake (M6.8). This period corresponds to the collapse vulnerability of Japanese wooden-frame houses, and is the major cause of severe damage during the Kumamoto earthquake. The response spectra also showed extremely large levels (>240–340 cm/s) in the long-period (T > 3 s) band, which is potentially disastrous for high-rise buildings, large oil storage tanks, etc. to have longer resonant period. Such long-period motion was, for the most parts, developed by the static displacement of the fault movement rather than by the seismic waves radiating from the source fault. Thus, the extreme near-fault long-period motion was hazardous only close to the fault but it attenuated very rapidly away from the fault.     

March 11, 2011 M 9.0 Tohoku earthquake in Japan: Tectonic setting of source,macroseismic, seismological,and geodynamic manifestations

The results of interpretation of seismological, geological, geophysical, geodetic, and macroseismic data on the source zone of the catastrophic Tohoku earthquake with M = 8.8–9.0 (from different estimates), which occurred March 11, 2011 off the eastern coast of Honshu Island, are reported. Consideration of the seismotectonic features of the Western Pacific; the distribution of epicenters and hypocenters of the main shock, fore- and aftershocks; the solution of focal mechanisms of the strongest shocks; and the data on directions of lateral and vertical displacement of the island surface makes it possible to contour the source region, reconstruct the structure of the source in the subsurface, and estimate the deformation of the lithosphere resulted from this great seismic event.     

Atypical soil hardening during the Tohoku earthquake of March 11, 2011 (Mw = 9.0)

Doklady Earth Sciences - Based on the records of KiK-net vertical arrays, models of soil behavior down to depths of ~100–200 m in the near-fault zones during the Tohoku earthquake are...     

Recovery of groundwater in the Sanriku region contaminated by the tsunami inundation from the 2011 Tohoku earthquake

Many wells in the Sanriku region used as sources for water supply systems were heavily contaminated by the tsunami of the 2011 great Tohoku earthquake on March 11 in 2011. To better understand the nature of the groundwater contamination by the tsunami inundation and to clarify the recovery process of contaminated groundwater at the study wells, groundwater monitoring has been conducted once or twice yearly since early summer in 2011. High and abnormal values of electric conductivity (EC), chloride ion concentration (CIC), Na⁺, Ca⁺, heavy metal ions, and heavier isotopes of the contaminated groundwater were also obtained in April and June 2011. The chemical elements have rapidly and exponentially decreased as a result of effective pumping of the contaminated groundwater from the study wells and because of abundant rainfall in 2011. In April 2015 (about 4 years after the tsunami inundation), the CIC and EC of the contaminated groundwater of two study wells in Minamisanriku town had reached pre-inundation values. The estimated residence times of groundwater of the two study wells were 105–118 days in the full-day pumping stage and 910–1000 days in the daytime-only pumping stage.     

Near-fault strong motion complexity of the 2000 Tottori earthquake (Japan) from a broadband source asperity model

The October 6/2000 Tottori earthquake that occurred in central Japan was an intermediate size strike-slip event that produced a very large number of near field strong motion recordings. The large amount of recorded data provides a unique opportunity for investigating a source asperity model of the Tottori earthquake that, combined with a hybrid strong motion simulation technique, is able to reproduce the observed broadband frequency near-fault ground motion.

We investigated the optimum source asperity parameters of the Tottori earthquake, by applying a Genetic Algorithm (GA) inversion scheme to optimise the fitting between simulated and observed response spectra and Peak Ground Acceleration (PGA) values. We constrained the initial model of our inversion by using the heterogeneous slip distribution obtained from a kinematic inversion of the source of previous studies. We used all the observed near-fault ground motions (−100 m) from the borehole strong motion network of Japan (KiK-Net), which are little affected by surficial geology (site effects).

The calculation of broadband frequency strong ground motion (0.1–10 Hz) is achieved by applying a hybrid technique that combines a deterministic simulation of the wave propagation for the low frequencies and a semi-stochastic modelling approach for the high frequencies. For the simulation of the high frequencies, we introduce a frequency-dependent radiation pattern model that efficiently removes the dependence of the pattern coefficient on the azimuth and take-off angle as the frequency increases. The good agreement between the observed and simulated broadband ground motions shows that our inversion procedure is successful in estimating the optimum asperity parameters of the Tottori earthquake and provides a good test for the strong ground motion simulation technique.

The ratio of background stress drop to average asperity stress drop from our inversion is nearly 50%, in agreement with the theoretical asperity model of Das and Kostrov [Das, S., Kostrov, B.V., 1986. Fracture of a single asperity on a finite fault: a model for weak earthquakes? Earthquake Source Mechanics, AGU, pp. 91–96.], and an empirical ratio of asperities to rupture area [Seismol. Res. Lett. 70 (1999) 59–80.].

The simulated radiation pattern is very complex for epicentral distances within half the fault length, but it approaches the radiation of a double-couple point source for larger distances.

The rupture velocity and rise time have a significant influence on the Peak Ground Velocity (PGV) distribution around the fault. An increase in rupture velocity produces a similar effect on the ground motion as a reduction in rise time.     

Explanation of liquefaction in after shock of the 2011 great east Japan earthquake using numerical analysis

During the 2011 Great East Japan Earthquake, severe liquefaction occurred in reclaimed ground in Urayasu city, Chiba prefecture. This liquefaction provided important lessons for us to re-recognize the liquefaction mechanism. A distinct feature of the liquefaction in this earthquake is that severe liquefaction happened not only in the main shock but also in an aftershock with a maximum acceleration of 25 gal. In some areas, liquefaction happened in the aftershock is even more serious than that happened in the main shock. In this paper, focus is placed on the characteristic features in the occurrence of liquefaction and consequent ground settlement. Based on the observed data, a series of dynamic–static analyses, considering not only the earthquake loading but also static loading during the consolidation after the earthquake shocks, are conducted in a sequential way just the same as the scenario in the earthquake. The calculation is conducted with 3D soil–water coupling finite element–finite difference analyses based on a cyclic elasto-plastic constitutive model. From the results of analyses, it is recognized that small sequential earthquakes, which cannot cause liquefaction of a ground in an independent earthquake vibration, cannot be neglected when the ground has already experienced liquefaction after a major vibration. In addition, the aftershock has great influence on the long-term settlement of low permeability soil layer. The observed and predicted liquefaction and settlements are compared and discussed carefully. It is confirmed that the numerical method used in this study can describe the ground behavior correctly under repeated earthquake shocks.     

The large aftershocks triggered by the 2011 Mw 9.0 Tohoku-Oki earthquake,Japan

The Mw 9.0 Tohoku-Oki earthquake that occurred off the Pacific coast of Japan on March 11, 2011, was followed by thousands of aftershocks, both near the plate interface and in the crust of inland eastern Japan. In this paper, we report on two large, shallow crustal earthquakes that occurred near the Ibaraki-Fukushima prefecture border, where the background seismicity was low prior to the 2011 Tohoku-Oki earthquake. Using densely spaced geodetic observations (GPS and InSAR datasets), we found that two large aftershocks in the Iwaki and Kita-Ibarake regions (hereafter referred to as the Iwaki earthquake and the Kita-Ibarake earthquake) produced 2.1 m and 0.44 m of motion in the line-of-sight (LOS), respectively. The azimuth-offset method was used to obtain the preliminary location of the fault traces. The InSAR-based maximum offset and trace of the faults that produced the Iwaki earthquake are consistent with field observations. The fault location and geometry of these two earthquakes are constrained by a rectangular dislocation model in a multilayered elastic half-space, which indicates that the maximum slips for the two earthquakes are 3.28 m and 0.98 m, respectively. The Coulomb stress changes were calculated for the faults following the 2011 Mw 9.0 Tohoku-Oki earthquake based on the modeled slip along the fault planes. The resulting Coulomb stress changes indicate that the stresses on the faults increased by up to 1.1 MPa and 0.7 MPa in the Iwaki and Kita-Ibarake regions, respectively, suggesting that the Tohoku-Oki earthquake triggered the two aftershocks, supporting the results of seismic tomography.     

Comparison of high-rate GPS,strong-motion records and their joint use for earthquake monitoring: a case study of the 2011 Mw 9.0 Tohoku earthquake

We explore the use of high-rate GPS, strong-motion records and their joint use for earthquake monitoring, the data collected during 2011 Mw 9.0 Tohoku earthquake was studied. We compared the recorded co-seismic movement, analyzed the displacements in both time-frequency domain. Meanwhile, the comparison of P wave detection was shown and the dynamic velocity waves were discussed. The results suggest that the GPS-only solution is good for low-frequency signal, and the strong-motion-only solution is good for high-frequency signal, thus, the integration of two datasets best complement the advantages of each, more details of co-seismic motions and broader frequency band of seismic signals. This is crucial for earthquake monitoring and early warning.