Iterative deconvolution of teleseismic P waves from the Thessaloniki (N. Greece) earthquake of June 20, 1978

Teleseismic long-period P waves from the June 20, 1978, Thessaloniki (N. Greece) earthquake (M _s=6.4) were modeled in an attempt to extract information about asperities or barriers on the fault plane. The analysis is based on the inversion method of complex P waves developed by Kikuchi and Kanamori (1982). A far-field source time function with a rise time of 2 sec and a process time of 5 sec is inferred, corresponding to a source dimension of about 10 km when a rupture velocity of 2 km/sec is assumed.The source depth of this shock, estimated by matching synthetic seismograms to observations, is found to be 8 km. The sum of the seismic moments of the individual subevents amounts to 3.3×10²⁵ dyn-cm.     

Probabilistic aftershock hazard assessment II: application of strong ground motion simulations

Probabilistic aftershock hazard assessment (PAHA) has been introduced by Wiemer (Geophys Res Lett 27:3405–3408, 2000). The method, in its original form, utilizes attenuation relations in evaluating peak ground velocity (PGV) exceedence probability. We substitute the attenuation relations together with their uncertainties by strong ground motion simulations for a set of scenarios. The main advantage of such an approach is that the simulations account for specific details of the aftershock source effects (faulting style, slip distribution, position of the nucleation point, etc.). Mean PGVs and their standard deviations are retrieved from the simulation results obtained by the new hybrid k-squared source model, and they are used for the PAHA analysis at a station under study. The model chosen for the testing purposes is inspired by the Izmit A25 aftershock (M _w  = 5.8) that occurred 26 days after the mainshock. The PAHA maps are compared with (1) those obtained by the use of attenuation relations and (2) the peak values of ten selected strong-motion recordings written by the aftershock at epicentral distances <50 km. We conclude that, although the overall hazard decay with increasing fault distance is similar, the PAHA maps obtained by the use of simulations exhibit remanent radiation pattern effect and prolongation in the strike direction due to the directivity effect pronounced for some of the scenarios. As regard the comparison with real data, we conclude that the PAHA maps agree with observed peak values due to appropriate attenuation model adopted in the analysis.     

Prediction of long-period ground motions from huge subduction earthquakes in Osaka,Japan

There is a high possibility of reoccurrence of the Tonankai and Nankai earthquakes along the Nankai Trough in Japan. It is very important to predict the long-period ground motions from the next Tonankai and Nankai earthquakes with moment magnitudes of 8.1 and 8.4, respectively, to mitigate their disastrous effects. In this study, long-period (>2.5 s) ground motions were predicted using an earthquake scenario proposed by the Headquarters for Earthquake Research Promotion in Japan. The calculations were performed using a fourth-order finite difference method with a variable spacing staggered-grid in the frequency range 0.05–0.4 Hz. The attenuation characteristics (Q) in the finite difference simulations were assumed to be proportional to frequency (f) and S-wave velocity (V _s) represented by Q = f · V _s / 2. Such optimum attenuation characteristic for the sedimentary layers in the Osaka basin was obtained empirically by comparing the observed motions during the actual M5.5 event with the modeling results. We used the velocity structure model of the Osaka basin consisting of three sedimentary layers on bedrock. The characteristics of the predicted long-period ground motions from the next Tonankai and Nankai earthquakes depend significantly on the complex thickness distribution of the sediments inside the basin. The duration of the predicted long-period ground motions in the city of Osaka is more than 4 min, and the largest peak ground velocities (PGVs) exceed 80 cm/s. The predominant period is 5 to 6 s. These preliminary results indicate the possibility of earthquake damage because of future subduction earthquakes in large-scale constructions such as tall buildings, long-span bridges, and oil storage tanks in the Osaka area.     

长周期地震动强度指标与SDOF体系残余变形相关性分析

为分析两类长周期地震动作用下如何选取合理的强度指标作为结构抗震设计的输入,从2011年东日本9.0级地震、2003年十胜冲8.0级地震和2016年熊本7.3级地震数据库中选取90条可靠的远场长周期地震动。从1994年美国北岭6.7级地震和1999年台湾集集7.6级地震数据库中选取60条近场长周期地震动,以SDOF体系为研究对象,讨论了阻尼比、屈服刚度折减系数和强度折减系数对残余变形与18个地震动强度指标的相关系数的影响,对比分析了两类长周期地震作用下相关系数之间的异同。研究结果表明:在考虑近场长周期地震动作用时,建议应根据结构周期的大小来选择合适的强度指标作为地震动的输入。远场长周期地震动作用下,以PGA、PGV和PGD为代表的强度指标与残余变形的相关程度均较高,PGV稳定性略好于PGA和PGD,建议PGV作为地震动输入的控制指标。残余变形相关系数受阻尼比、屈服刚度折减系数以及强度折减系数影响不大。     

A seismological overview of long-period ground motion

Long-period ground motion has become an increasingly important consideration because of the recent rapid increase in the number of large-scale structures, such as high-rise buildings and oil storage tanks. Large subduction-zone earthquakes and moderate to large crustal earthquakes can generate far-source long-period ground motions in distant sedimentary basins with the help of path effects. Near-fault long-period ground motions are generated, for the most part, by the source effects of forward rupture directivity. Far-source long-period ground motions consist primarily of surface waves with longer durations than near-fault long-period ground motions. They were first recognized in the seismograms of the 1968 Tokachi-oki and 1966 Parkfield earthquakes, and their identification has been applied to the 1964 Niigata earthquake and earlier earthquakes. Even if there is no seismogram, we can identify far-source long-period ground motions through the investigation of tank damage by liquid sloshing.     

Constituents of Vertical-component Coda Waves at Long Periods

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⁻³ to (2.38±0.32) × 10⁻³ 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.     

Ray parameters of diffracted long period P and S waves and the velocities at the base of the mantle

The derivation of P and S velocities at the core-mantle boundary (CMB) from long-period diffracted waves by the use of the simple ray-theoretical formulav _CMB=r _c /p (v _CMB=velocity at the CMB;r _c =core radius;p=ray parameter) yields apparent velocity values which differ from the true velocities. Using a dominant period of about 20 sec for calculating theoretical seismograms, we found a linear relation between the apparent velocity and the average velocity in a transition zone at the base of the mantle with fixed velocity on top.The ray parameters determined from long-period earthquake data are found to be 4.540±0.035 and 8.427±0.072 sec/deg for P_diff and S_diff, respectively. These values yield apparent velocities of 13.378±0.103 for P and 7.207±0.062 km/sec for S waves. By means of the theoretical relation between apparent and average velocity and under the assumption of linear variation of velocity with depth, one can invert the apparent velocities into true CMB velocities of 13.736±0.170 and 7.320±0.124 km/sec. These results imply positive velocity gradients at the base of the mantle and hence no significant departures from adiabaticity and homogeneity.Contribution No. 211 of the Geophysical Institute, University of Karlsruhe.     

Broadband Ground Motion Reconstruction for the Kanto Basin during the 1923 Kanto Earthquake

Ground motions of the 1923 Kanto Earthquake inside the Kanto Basin are numerically simulated in a wide frequency range (0?C10?Hz) based on new knowledge of the earthquake??s source processes, the sedimentary structure of the basin, and techniques for generating broadband source models of great earthquakes. The Kanto Earthquake remains one of the most important exemplars for ground motion prediction in Japan due to its size, faulting geometry, and location beneath the densely populated Kanto sedimentary basin. We reconstruct a broadband source model of the 1923 Kanto Earthquake from inversion results by introducing small-scale heterogeneities. The corresponding ground motions are simulated using a hybrid technique comprising the following four calculations: (1) low-frequency ground motion of the engineering basement, modeled using a finite difference method; (2) high-frequency ground motion of the engineering basement, modeled using a stochastic Green??s function method; (3) total ground motion of the engineering basement (i.e. 1?+?2); and (4) ground motion at the surface in response to the total basement ground motion. We employ a recently developed three-dimensional (3D) velocity structure model of the Kanto Basin that incorporates prospecting data, microtremor observations and measurements derived from strong ground motion records. Our calculations reveal peak ground velocities (PGV) exceeding 50?cm/s in the area above the fault plane: to the south, where the fault plane is shallowest, PGV reaches 150?C200?cm/s at the engineering basement and 200?C250?cm/s at the surface. Intensity 7, the maximum value in the Japan Meteorological Agency??s intensity scale, is calculated to have occurred widely in Sagami Bay, which corresponds well with observed house-collapse rates due to the 1923 event. The modeling reveals a pronounced forward directivity effect for the area lying above the southern, shallow part of the fault plane. The high PGV and intensity seen above the southeastern corner of the fault plane and further east are largely due to this effect. Waveforms above the fault plane contain both short- and long-period components, but the short-period components are not observed further afield. Away from the fault, long-period waves (>2?s) dominate the ground motion, and in areas where the base of the third layer is relatively deep, the predominant period is >5?s. Levels of long-period ground motion in the southern part of the study area, around Sagami Bay and the southern parts of Boso Peninsula and Tokyo Bay, exceed that recorded at Tomakomai during the 2003 Tokachi-oki earthquake, when large oil storage tanks collapsed in response to sloshing generated by strong long-period motions.     

Development of Long-period Ground Motions from the Nankai Trough, Japan, Earthquakes: Observations and Computer Simulation of the 1944 Tonankai (Mw 8.1) and the 2004 SE Off-Kii Peninsula (Mw 7.4) Earthquakes

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.     

Ground motion attenuation relation for small to moderate earthquakes in Fujian region, China

We collect 1974 broad-band velocity records of 94 earthquakes (ML=2.8~4.9, △=13~462 km) from seven stations of the Fujian Seismic Network from March 1999 to March 2007. Using real-time simulation, we obtain the corresponding acceleration and then adopt different models to analyze the seismic data. As a result, a new attenuation relationship between PGA and PGV of the small and moderate earthquakes on bedrock site in Fujian region is established. The Yongchun earthquake occurred recently verifies the attenuation relationship well. This paper provides a new approach for studying the ground motion attenuation relationship using velocity records.     

熊本MW7.0地震近场地表与井下地震动对比研究

选取日本熊本M_W7.0地震断层距小于200 km的82个近场KiK-net台站记录到的三分量记录数据进行基线校正后,获得近场地面运动水平向的峰值加速度PGA、峰值速度PGV及周期为0.2,1,2,3,5和10 s的加速度反应谱数据,并与美国NGA-West2的地震动预测模型相比较,研究熊本地震地表和井下地震动峰值及反应谱的衰减特征,通过比较KiK-net台站地表与井下记录结果,探讨浅层场地放大效应的影响。研究结果表明:① 对于井下观测结果,NGA-West2的地震动模型对PGA和短周期0.2 s的反应谱的预测值与井下观测值相比整体偏高,而PGV和较长周期地震动（如1,2和3 s的反应谱）的预测值与井下观测值较为吻合;② 地表观测记录的PGA,PGV和周期为0.2—3 s的反应谱残差整体上随v_S30对数值的增大呈线性减小的趋势,而周期为5 s和10 s的长周期部分,其场地效应的影响很小;③ 相对于井下记录,地表记录的地震动PGA,PGV和周期为0.2,1和2 s的反应谱有明显的放大,这种放大作用随浅层场地剪切波速的增大而减小;周期为3,5和10 s时长周期地震动的放大效应很小。      

Source Parameters of the Deadly M<Subscript>w</Subscript> 7.6 Kashmir Earthquake of 8 October, 2005

During the last six years, National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5–8 broadband seismographs and 10–20 accelerographs in the Kachchh seismic zone, Gujarat with a prime objective to monitor the continued aftershock activity of the 2001 M_w 7.7 Bhuj mainshock. The reliable and accurate broadband data for the 8 October M_w 7.6 2005 Kashmir earthquake and its aftershocks from this network as well as Hyderabad Geoscope station enabled us to estimate the group velocity dispersion characteristics and one-dimensional regional shear velocity structure of the Peninsular India. Firstly, we measure Rayleigh-and Love-wave group velocity dispersion curves in the period range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of Peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8 km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

Lessons from the 2003 Tokachi-oki, Japan, earthquake for prediction of long-period strong ground motions and sloshing damage to oil storage tanks

The 2003 Tokachi-oki earthquake (M _w 8.0) in northern Japan generated large-amplitude long-period (4–8 s) ground motions in the Yufutsu sedimentary basin, causing severe damage to seven large oil storage tanks with floating roof structures because of severe sloshing of oil. The 30,000–40,000-m³ tanks having suffered the severe damage such as fires and sinking of floating roofs experienced the sloshing with large amplitudes exceeding 3 m in which the fundamental mode was predominant. The second mode of sloshing was also excited in the 110,000-m³ tanks in which their floating roofs sank into oil, indicating that the higher modes of sloshing as well as the fundamental mode should be considered in damage prediction. The strong ground motion recordings demonstrated the earthquake dependency of predominant periods and the substantial spatial variation of the long-period shaking observed within the Yufutsu basin, meaning the necessity of source- and site-specific prediction of long-period strong ground motions. The two-dimensional numerical modeling suggested the importance of detailed structures of soft near-surface sediments as well as deep basin structure for accurate prediction of long-period strong ground motions in deep sedimentary basins.     

用体波合成地震图方法确定渤海、永善两大地震的震源参数.

用体波合成地震图与观测地震图拟合的方法确定了渤海地震（1969.7.18. M=7.4）和云南永善地震（1974.5.11. M=7.1）的震源参数.针对震中距为30至90的远震 P+pP+sP 波,采用 WKBJ 近似和实路径积分,做出了20个 WWSSN 台站的渤海地震和7个WWSSN 台站的永善地震的全波理论合成地震图.与实际记录到的长周期垂直分量地震图进行比较,以试错法确定两大地震的震源参数如下:渤海地震:h=25km,=205,=85,=-145,M0=2.41019Nm,R=25,L=36km,vR=6km/s;永善地震由两个同时发生的 h=15km 的子破裂组成,其一为=45,=65,=0,M0=1.91019Nm,R=225,L=10km,vR=1km/s;其二为=200,=75,=170,M0=1.31019Nm,R=30,L=23km,vR=2km/s.这里 h 为震源深度,为走向,为倾角,为滑动角,M0 为总地震矩,R 为断裂扩展方向,L 为断层长度,vR 为破裂传播速度.      

南北地震带南段水平向地震动衰减特征

地震动衰减关系作为抗震救灾的重要依据,一直都是现代地震学研究的重点之一.本研究使用南北地震带南段区域2009—2016年共217个地震事件获得随震中距变化的水平向地震动峰值速度(PGV)和地震动峰值加速度(PGA)经验衰减关系,并计算场地响应.研究结果显示PGV衰减关系的拟合效果较PGA更好,两者的距离衰减系数会随事件矩震级增大呈线性减小;相较于大矩震级事件,中矩震级事件在近场可能产生较衰减关系理论值更大的PGV和PGA,同时衰减关系的拟合标准差会随事件矩震级的增大而减小.进行场地响应校正后的PGV和PGA更加符合经验衰减关系,PGA的场地响应影响较PGV更强但两者的趋势一致,并与该区域前人计算得到的地壳Qs值分布对应,表明地壳介质放大或压制地震波振幅和其传递地震波能量的能力是相关联的.本文结果一定程度上揭示了南北地震带南段的地震动强度衰减特征,为未来中国西南部的抗震减灾工作提供了重要的参考.     

Finite-element simulations of long-period ground motions: Japanese subduction-zone earthquakes and the 1906 San Francisco earthquake

Large earthquakes at shallow depths commonly excite long-period ground motions in distant sedimentary basins, thereby inflicting damage upon large-scale structures. For example, the 2003 Tokachi-oki earthquake, Japan, damaged oil tanks in the Yufutsu Basin, located 250 km from the epicenter. Similar long-range effects were also observed during the 2004 earthquake off Kii Peninsula, Japan. In this study, we present the results of simulations of these earthquakes undertaken using the finite element method (FEM) with a voxel mesh. In addition, to examine whether the 1906 San Francisco earthquake excited long-period ground motions in the Los Angeles-area basins, we performed long-period ground motion simulations of most of the California region. The FEM simulations confirmed the importance of path effects for the development of long-period ground motions.     

Source Parameters of the 8 October, 2005 M<Subscript>w</Subscript>7.6 Kashmir Earthquake

During the last six years, the National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5 broadband seismographs and 10 accelerographs in the Kachchh seismic zone, Gujarat, with the prime objective to monitor the continued aftershock activity of the 2001 M_w7.7 Bhuj mainshock. The reliable and accurate broadband data for the M_w 7.6 (8 Oct., 2005) Kashmir earthquake and its aftershocks from this network, as well as from the Hyderabad Geoscope station, enabled us to estimate the group velocity dispersion characteristics and the one-dimensional regional shear-velocity structure of peninsular India. Firstly, we measure Rayleigh- and Love-wave group velocity dispersion curves in the range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8-km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

Three-dimensional simulations of strong ground motion in the Shidian basin based upon the spectral-element method

The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical profiles within the basin are used to analyze the generation and propagation of surface waves and their relation to the basin structures in both the time and frequency domain.The amplification effects are analyzed by the distribution of peak ground velocity(PGV)and cumulative kinetic energy(Ek) in the basin.The results show that in the 3D basin,the excitation of the fundamental and higher surface wave modes are similar to that of the 2D model.Small bowls in the basin have great influence on the amplification and distribution of strong ground motion,due to their lateral resonances when the wavelengths of the lateral surface waves are comparable to the size of the bowls.Obvious basin edge effects can be seen at the basin edge closer to the source for constructive interference between direct body waves and the basin-induced surface waves.The Ek distribution maps show very large values in small bowls and some corners in the basin due to the interference of waves propagating in different directions.A high impedance contrast model can excite more surface wave modes,resulting in longer shaking durations as well as more complex seismograms and PGV and Ek distributions.     

2020年10月22日四川北川MS4.7地震强震动特征分析

2020年10月22日11时03分37秒四川省绵阳市北川县发生MS4.7地震,四川强震动台网与预警烈度速报台网在震区建成较密集的台站,获取了532组三分量加速度记录,有助于开展区域地震动衰减和地震动特征研究.本文对强震记录进行常规处理后计算出强震动记录的相关参数,利用克里金插值方法得到地震动峰值加速度PGA和峰值速度PGV的空间分布图,长轴呈北西—南东方向.分析强震动记录PGA、PGV随距离的衰减规律,与常用衰减关系预测值进行对比,此次地震PGA的衰减特性与俞言祥和汪素云(2006)提出的中国西部地区水平向基岩加速度衰减关系有较好的一致性.北川MS4.7地震获得的密集强震动记录为建立区域衰减关系,以及开展基于经验格林函数方法(EGFM)再现大震强地震动场展布等研究提供了重要的数据支撑.     

Ground-motion modelling at regional distances for earthquakes in a continental interior,I. Theory and observations

To understand the ground-motion contribution by multiple-mode surface-wave arrivals, the surface-wave theory required for predicting ground-motion time histories is discussed. The adequacy of the theory in accounting for observed earthquake ground motion is tested by comparing theoretically predicted long-period seismograms with real seismograms for two earthquakes in the central United States. The agreement between the predicted and observed maximum vertical component L_g ground velocities and accelerations in the 2- to 3-sec period range is excellent over a distance range of 100 to 2,000 km.