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.     

Effect of frequency-dependent radiation pattern in the strong motion simulation of the 2011 Tohoku earthquake,Japan, using modified semi-empirical method

We perform a strong ground motion simulation using a modified semi-empirical technique (Midorikawa in Tectonophysics 218:287–295, 1993), with frequency-dependent radiation pattern model. Joshi et al. (Nat Hazards 71:587–609, 2014) have modified the semi-empirical technique to incorporate the modeling of strong motion generation areas (SMGAs). A frequency-dependent radiation pattern model is applied to simulate high-frequency ground motion more precisely. Identified SMGAs (Kurahashi and Irikura in Earth Planets Space 63:571–576, 2011) of the 2011 off the Pacific coast of Tohoku earthquake (M _w  = 9.0) were modeled using this modified technique. We analyzed the effect of changing seismic moment values of SMGAs on the simulated acceleration time series. Final selection of the moment values of SMGAs is based on the root-mean-square error (RMSE) of waveform comparison. Records are simulated for both frequency-dependent and constant radiation pattern function. Simulated records for both cases are compared with observed records in terms of peak ground acceleration, peak ground velocity and pseudo-acceleration response spectra at different stations. Comparison of simulated and observed records in terms of RMSE suggests that the method is capable of simulating record, which matches in a wide frequency range for this earthquake and bears realistic appearance in terms of shape and strong motion parameters. The results confirm the efficacy and suitability of rupture model defined by five SMGAs for the developed modified technique.     

Use of site amplification and anelastic attenuation for the determination of source parameters of the Sikkim earthquake of September 18, 2011, using far-field strong-motion data

An earthquake of magnitude 6.9 (M _w) occurred in the Sikkim region of India on September 18, 2011. This earthquake is recorded on strong-motion network in Uttarakhand Himalaya located about 900 km away from the epicenter of this earthquake. In this paper acceleration record from six far-field stations has been used to compute the source parameters of this earthquake. The acceleration spectra of ground motion at these far-field stations are strongly affected by both local site effects and near-site anelastic attenuation. In the present work the spectrum of S-phase recorded at these far-field stations has been corrected for anelastic attenuation at both source and site and the site amplification terms. Site amplifications at different stations and near-site shear wave attenuation factor have been computed by the technique of inversion of acceleration spectra given by Joshi et al. (Pure Appl Geophys 169:1821–1845, 2012a). For estimation of site amplification and shear wave quality factor [Q _β (f)] at the recording sites, ten local events recorded at various stations between July 2011 and December 2011 have been used. The obtained source spectrum from acceleration records is compared with the theoretical source spectrum defined by Brune (J Geophys Res 76:5002, 1970) at each station for both horizontal components of the records. Iterative forward modeling of theoretical source spectrum gives the average estimate of seismic moment (M _o), source radius (r _o) and stress drop (Δσ) as (3.2 ± 0.8) × 10²⁶ dyne cm, 13.3 ± 0.8 km and 59.2 ± 8.8 bars, respectively, for the Sikkim earthquake of September 18, 2011.     

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.     

Estimation of model parameter of Sumatra earthquake using empirical Green’s function technique and generation of hypothetical earthquake scenario for Andaman Island, India

Empirical Green??s function (EGF) technique is considered to be most effective technique for simulation of ground motions due to a finite earthquake source. In the present paper, this technique has been used to simulate ground motion due to a great earthquake. The coastal region of Sumatra Island has been visited by a great earthquake on December 26, 2004. This earthquake has been recorded at several broadband stations including a nearest broadband station PSI in Indonesia. The shear wave contributions in both horizontal components have been simulated at PSI station using EGF technique. The comparison of simulated and observed waveform has been made for various possibilities of rupture parameters in terms of root mean square error. The final rupture model supports rupture velocity of 3.0?km/s with nucleation point supporting northward propagating rupture that coincide with high-slip asperity defined by Sorensen et al. (Bull Seism Soc Am 97:S139?CS151, 2007). The final modeling parameters have been used to simulate record at MDRS station in coastal state of Tamilnadu, India. In an attempt to model a scenario of great earthquake in the Andaman Island, a hypothetical rupture plane is modeled in this region. The event occurred on August 10, 2008 of magnitude 6.2 (M _w ) recorded on strong motion array at Port Blair has been used as EGF to simulate records due to the hypothetical great earthquake. Possibilities of earthquake due to the oblique strike-slip and thrust mechanism have been modeled in the present paper. Several possibilities of nucleation point for both cases has been considered, and it is seen that variation of peak ground acceleration at Port Blair station for strike-slip and thrust mechanism is 126?C738 gals and 647?C2,571 gals, respectively, which indicate high seismic hazard potential of Andaman Island.     

Strong ground motion modelling of causative fault for the 2002 Avaj earthquake, Iran

The suitability of a very fast method for obtaining synthesizing accelerograms has been demonstrated for a hybrid simulation technique of source wavelet and acceleration envelope waveform for the 2002 Avaj earthquake. This method is based on the amplitude modeled white noise and envelope waveform. The estimation of peak acceleration from a preliminary simulated record is based on using modeling parameters of rupture plane instead of empirical relations for peak acceleration. Based on comparison between observed and simulated strong ground motion data, a fair agreement is observed between simulated and observed records up to distances 40 km for peak acceleration and duration. The most important feature of the recorded strong motion is decay up to a distance of 40 km which is due to direct upgoing shear waves. At distance of 50 to 60 km peak acceleration increase, which is due to postcritical reflection from velocity gradient in the lower crust. A flat trend is observed for peak acceleration at distance of 60 to 100 km. The simulation indicates that the rupture is started at depth of 8 km and propagated from northwest to southeast. The causative fault for the 2002 Avaj earthquake shows similar mechanism to the 1962 Buin-Zahra earthquake.     

汶川Ms8.0地震强震动基线改正及其在位错反演中的初步应用

通过研究近场强震动记录, 发现汶川M_s8.0地震近场峰值加速度在空间上存在较明显的上盘效应和方向性效应, 与汶川引起的地质灾害空间分布具有较好的一致性.但在所有强震仪所记录的汶川M_s8.0地震同震加速度记录积分所得地壳同震速度中, 有的台站数据存在典型的线性偏移, 有的台站数据除线性偏移外还存在明显的非线性偏移.采用非线性基线改正方法处理汶川M_s8.0强震同震记录, 改正后所得同震位移明显要比线性基线改正更合乎实际情况.以强震动、GPS和InSAR同震位移处理结果做约束, 反演了汶川M_s8.0地震同震位错分布, 对于汶川M_s8.0地震主要同震破裂断裂(北川-映秀断裂), 强震动反演结果不仅较好地刻画了汶川M_s8.0地震同震主断裂上地表破裂空间分布详细变化特征, 同时也较好地反映北端破裂衰减情况, 该结果表明: 强震动资料可以为强震后的救援和灾害评估等工作提供具有参考价值的研究结果; 另一方面, 受数据数量的制约, 用强震动改正后位移反演所得位错分布中仅汉旺断裂南段存在较为明显位错, 强震仪布设时应更多地考虑是否相对均匀地分布在具有发震潜势的断裂周缘, 以期更好地在震后应急救灾中发挥更好的作用.      

汶川M_s 8.0级地震强震记录所揭示的地震断层特征分析

基于2008年5月12日汶川Ms8.0级地震的强震台站记录,对比发震断层两侧的峰值加速度与地表破裂带上同震位移的分布特点,探讨了地震动强度分布特征与地表破裂位移分布之间的相关关系。分析近断层典型台站的强震动记录时程特征,获得了强地震动记录中所包含的断层破裂过程和破裂习性信息,从强震观测记录的角度进一步证实了汶川地震主震的多次破裂特征。结果表明,汶川地震主震至少包含了4次地震破裂事件,最主要的前两次破裂事件分别对应映秀—北川断裂段和北川—南坝断裂段的破裂过程,后两次破裂事件释放的能量相对较小,应该是第二次破裂过程触发局部次级破裂所引起的。此外,垂直于断层的峰值加速度剖面揭示的发震断层的高倾角逆冲特性,与地震地质调查和小震精定位等确定的相应结果是一致的。     

Near-source strong motion database catalog for Iran

This paper discusses a newly developed high-quality integrated dataset of shallow earthquake ground motions that occurred in Iran, from 1976 to 2013. A total of 860 three-component strong motion records are processed from 183 earthquake events, moment magnitudes 5.0?≤?M _w ?≤?7.4, and rupture distances of R _RUP  ≤ 120 km. Strong motion data from Iran having special tectonic features and shallow earthquakes with depths less than 35 km are included. This paper presents a thorough procedure used to collect and to generate a database following the Next-Generation Attenuation-West research projects. This database can be used in the development and ranking of ground motion models and for seismological and engineering hazard and risk analyses. Unprocessed strong motion records are obtained from the Iranian Strong Motion Network (ISMN). The time series collected were thoroughly examined through several rounds of quality reviews. The newly generated database includes the peak ground acceleration, peak ground velocity, and pseudo-spectral acceleration for the 5% damped with periods ranging from 0.01 to 10 s. The database also includes ground motion information and source characterization and parameters. This study is the near-source compiled ground motion database that can be used for Iran, and it is consistent with standard worldwide databases.     

Hybrid attenuation model for estimation of peak ground accelerations in the Kutch region, India

The Kutch region of Gujarat in India is the locale of one of the most devastating earthquake of magnitude (M _w) 7.7, which occurred on January 26, 2001. Though, the region is considered as seismically active region, very few strong motion records are available in this region. First part of this paper uses available data of strong motion earthquakes recorded in this region between 2006 and 2008 years to prepare attenuation relation. The developed attenuation relation is further used to prepare synthetic strong motion records of large magnitude earthquakes using semiempirical simulation technique. Semiempirical simulation technique uses attenuation relation to simulate strong ground motion records of any target earthquake. The database of peak ground acceleration obtained from simulated records is used together with database of peak ground acceleration obtained from observed record to develop following hybrid attenuation model of wide applicability in the Kutch region: $$ \begin{aligned} \ln \left( {\text{PGA}} \right) & = - 2.56 + 1.17 \, M_{\text{w}} - \, 0.015R - 0.0001\ln \left( {E + 15} \right) \\ &\quad 3.0 \le M_{\text{w}} \le 8.2;\quad 12 \le R \le 120;\quad {\text{std}} . {\text{ dev}}.(\sigma ): \pm 0.5 \\ \end{aligned} $$ ln ( PGA ) = ? 2.56 + 1.17 M w ? 0.015 R ? 0.0001 ln ( E + 15 ) 3.0 ≤ M w ≤ 8.2 ; 12 ≤ R ≤ 120 ; std . dev . ( σ ) : ± 0.5 In the above equation, PGA is maximum horizontal ground acceleration in gal, M _w is moment magnitude of earthquake, R is hypocentral distance, and E is epicentral distance in km. The standard deviation of residual of error in this relation is 0.5. This relation is compared with other available relations in this region, and it is seen that developed relation gives minimum root mean square error in comparison with observed and calculated peak ground acceleration from same data set. The applicability of developed relation is further checked by testing it with the observed peak ground acceleration from earthquakes of magnitude (M _w), 3.6, 4.0, 4.4, and 7.7, respectively, which are not included in the database used for regression analysis. The comparison demonstrates the efficacy of developed hybrid attenuation model for calculating peak ground acceleration values in the Kutch region.     

基于包络的汶川大地震高频地震波辐射区域反演及近场加速度合成

合理地估计汶川破坏区域的地震动有助于地震灾害的研究.通过利用芦山地震记录建立的加速度包络衰减关系和汶川地震近场30个台站的加速度包络,基于线源模型,采用差分进化方法反演了汶川地震断层面上高频 (＞1 Hz) 辐射区域分布.结果表明:断层面上高频辐射分布很不均匀,辐射较强的区域主要位于:(1) 产生较大地表破裂的映秀、北川和南坝区域;(2) 映秀和北川等凹凸体的周边区域,包括震中东北侧60~90 km区域、北川和南坝东北侧30 km处;(3) 断层破裂停止的东北端约30 km长的区域.其中,破裂贯穿到地表的映秀、北川和南坝是低频和高频辐射都很强的区域.对于无观测记录场点,选择其临近且场地条件类似的台站加速度提取平稳随机过程,结合高频辐射分布和衰减关系得到的包络,合成了加速度时程,可为汶川地震结构震害分析提供地震动输入.      

Correlation between ground motion parameters and lining damage indices for mountain tunnels

In order to study the interdependencies between the overall lining damage indices for mountain tunnels and ground motion parameters, a total of 89 ground motion records were selected for analysis. These records were divided into three groups: (1) near-field ground motion with velocity pulses, (2) near-field ground motion without velocity pulses, and (3) far-field ground motion. Calculations of Pearson’s linear correlation coefficient have been used to clarify the grade of interrelationships between ground motion parameters. Further, nonlinear dynamic analyses were carried out on mountain tunnels to determine overall lining damage indices under seismic actions. Based on numerical results, linear correlation coefficients between ground motion parameters and overall lining damage indices for mountain tunnels were then calculated and analyzed. Overall lining damage indices were found to be highly correlated to the velocity-related seismic parameters but not so well correlated to spectra-related parameters.     

Flood risk assessment in a coastal lagoon under present and future scenarios: Ria de Aveiro case study

Iranian strong motion records as well as detailed conditions of their instrument sites and the characteristics of their causative seismic sources are compiled and processed. The dataset consists of 2286 three-component records from 461 Iranian earthquakes with at least two high-quality records having moment magnitude from 3.9 to 7.3. These records are about 20% of the Iranian database and are suitable for seismic hazard analysis and engineering applications. Perhaps for the first time in the literature, the distance to the surface projection of the fault is reported for a great number of records corresponding to earthquakes with M > 6.0. The raw accelerations are processed using the wavelet de-noising method. Having corrected and filtered these raw data, the pseudospectral accelerations are calculated for each of the three components of time series, separately. In addition to the ground motion parameters, a large and comprehensive list of metadata characterizing the recording conditions of each record is also developed. Moreover, careful revision of the characteristics of the earthquakes such as location, magnitude, style of faulting and fault rupture plane geometry, if available, is carried out using the best available information in a scientifically sound manner. Finally, we also focus on special ground motion records including records with peak ground acceleration (PGA) >300 cm/s² and distances less than 30 km. These are “exceptional” records in the Iranian dataset and include less than 2% of the selected dataset.     

Strong ground motion attenuation in Aswan area, Egypt

Recent and paleo seismicity indicate that moderate seismic activity is relatively large for Aswan area. This is a warning on the possibility of occurrence of earthquakes in the future too. No strong motion records are available in Aswan area for engineers to rely upon. Consequently, the seismological modeling is an alternative approach till sufficient instrumental records around Aswan become available. In the present study, we have developed new ground motion attenuation relationship for events spanning 4.0?≤? M _w?≤?7.0 and distance to the surface projection of the fault up to 100 km for Aswan based on a statistically simulated seismological model. We generated suites of ground motion time histories using stochastic technique. The ground motion attenuation relation describes the dependence of the strength of the ground motions on the earthquake magnitude and distance from the earthquake. The proposed equation for peak ground acceleration (PGA) for the bed rock is in the form of: $ {\mathbf{log}}{\text{ }}\left( {{\mathbf{PGA}}/{\mathbf{gal}}} \right){\text{ }} = {\mathbf{1}}.{\mathbf{24}} + {\mathbf{0}}.{\mathbf{358}}{M_{\mathbf{w}}} - {\text{ }}{\mathbf{log}}\left( {\mathbf{R}} \right){\text{ }}-{\text{ }}{\mathbf{0}}.{\mathbf{008}}{\text{ }}{\mathbf{R}}{\text{ }} + {\text{ }}{\mathbf{0}}.{\mathbf{22}}{\text{ }}{\mathbf{P}} $ . Where PGA is the peak ground acceleration in gal (cm/s²); M_w, its moment magnitude; R is the closest distance between the rupture projection and the site of interest; and the factor P is a dummy variable. It is observed that attenuation of strong motion in Aswan is correlated with those used before in Egypt.     

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.     

近源地震动峰值加速度衰减关系影响因素分析

本文收集了丰富的强震资料, 以峰值加速度为例, 采用简单且体现近场峰值加速度PGA震级饱和和距离饱和特性的衰减模型, 研究了表征震级与距离饱和效应的R0(M)的性质。R0与震级相关, 同时与震源性质、地震波频谱有关。在单个地震的R0(M)的求取中, 由于R0与系数d几乎呈线性关系, 所以要求单个地震R0的值, 必须先根据理论约束确定d的大小。在检验衰减方程的预测效果时, 不仅要判断衰减曲线是否反映了实测资料的平均变化趋势, 而且要判断实测资料是否绝大多数落在84%及16%概率水平的预测曲线之内(之间).     

Rupture line and damage effects in Aqaba area, Egypt

Strong-motion parameters, peak ground acceleration, peak ground velocity, and peak ground displacement depend on several factors, such as the source of earthquake, distance between the source and site, and the characteristics of that site. Five seismographs and two accelographs were installed by the Egyptian Geological Survey team along the western side of the Gulf of Aqaba from Taba to Sharm El-Sheikh to record seismic events during the period from September 1995 to June 1996. During this period, two events were recorded by the accelographs. The two accelographs were located on the surface of the basement rocks, the first near the epicenter and the second at a farther distance. However, the farthest accelograph recorded higher values compared to the nearer one. Fault mechanics are an important factor in determining the values of strong motion parameters, where the direction of the rupture line plays an important role in detecting the values of strong motion parameters, the strong motion parameters and damage effects seem to be attenuated very fast in the direction perpendicular to the rupture line. This can be interpreted by the fact that the farthest accelograph lies at the extension of the fault rupture (azimuth = 30°), while the nearer one was perpendicular to the strike of the fault rupture.     

Region-specific deterministic and probabilistic seismic hazard analysis of Kanpur city

A seismic hazard map of Kanpur city has been developed considering the region-specific seismotectonic parameters within a 500-km radius by deterministic and probabilistic approaches. The maximum probable earthquake magnitude (M _max) for each seismic source has been estimated by considering the regional rupture characteristics method and has been compared with the maximum magnitude observed \(\left ({M_{\max }^{\text {obs}}}\right )\), \(M_{\max }^{\text {obs}} +0.5\) and Kijko method. The best suitable ground motion prediction equations (GMPE) were selected from 27 applicable GMPEs based on the ‘efficacy test’. Furthermore, different weight factors were assigned to different M _max values and the selected GMPE to calculate the final hazard value. Peak ground acceleration and spectral acceleration at 0.2 and 1 s were estimated and mapped for worst-case scenario and 2 and 10% probability of exceedance for 50 years. Peak ground acceleration (PGA) showed a variation from 0.04 to 0.36 g for DSHA, from 0.02 to 0.32 g and 0.092 to 0.1525 g for 2 and 10% probability in 50 years, respectively. A normalised site-specific design spectrum has been developed considering three vulnerable sources based on deaggregation at the city center and the results are compared with the recent 2011 Sikkim and 2015 Nepal earthquakes, and the Indian seismic code IS 1893.     

Modeling and synthesis of strong ground motion

Success of earthquake resistant design practices critically depends on how accurately the future ground motion can be determined at a desired site. But very limited recorded data are available about ground motion in India for engineers to rely upon. To identify the needs of engineers, under such circumstances, in estimating ground motion time histories, this article presents a detailed review of literature on modeling and synthesis of strong ground motion data. In particular, modeling of seismic sources and earth medium, analytical and empirical Green’s functions approaches for ground motion simulation, stochastic models for strong motion and ground motion relations are covered. These models can be used to generate realistic near-field and far-field ground motion in regions lacking strong motion data. Numerical examples are shown for illustration by taking Kutch earthquake-2001 as a case study.     

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.