Spectral parameters of ground motion in different regions: comparison of empirical models

Fourier-amplitude spectrum is one of the most important parameters describing earthquake ground motion, and it is widely used for strong ground motion prediction and seismic hazard estimation. The relationships between Fourier-acceleration spectra, earthquake magnitude and distance were analysed for different seismic regions (the Caucasus and Taiwan island) on the basis of ground motion recordings of small to moderate (3.5≤M_L≤6.5) earthquakes. It has been found that the acceleration spectra of the most significant part of the records, starting from S-wave arrival, can be modelled accurately by the Brune's “ω-squared” point-source model. Parameters of the model are found to be region-dependent. Peak ground accelerations and response spectra for condition of rock sites were calculated using stochastic simulation technique and obtained models of source spectra. The modelled ground-motion parameters are compared with those predicted by recent empirical attenuation relationship for California.     

A note on computed magnitude ML from strong motion records

It is shown that a new definition of M_L^SM (Trifunac¹⁴), which is computed from recorded strong motion earthquake accelerograms, leads to stable and unbiased estimates of the local earthquake magnitude M_L (Richter⁸) for epicentral distance R < 100 km and for 3 M_nL < 6.5. Tables of the uniform estimates using this new M_L^SM are presented for all earthquakes contributing to the current data base of free field strong motion accelerograms in EQINFOS files (Lee and Trifunac⁶).     

2018年11月25日新疆博乐MS4.9地震余震检测与目录完整性研究

本文基于匹配滤波技术，通过SEPD（Seismic Events and Phase Detection）对2018年11月25日新疆博乐M_S4.9地震序列进行检测，检测出遗漏地震32条，84.4%地震为M_L0.0—1.0，9.4%地震小于M_L0.0，较地震目录中原有15条地震多213%，检测出的遗漏地震事件使地震目录更加完整。检测后的最小完整性震级由检测前的M_L1.6减至M_L0.8，地震目录最小完整性震级的减小有利于地震工作者对区域地震活动性作出更准确全面的结论，并使地震危险性分析更可靠。     

ML of the 1989 Loma Prieta earthquake

The use of regional attenuation in computing the local magnitude, M_L, from strong motion data gathered at distances less than 100 km may lead to systematic underestimates approaching 0·5 magnitude units (Trifunac & Herak, Soil Dynamics and Earthquake Engineering, 1992, 18, 229-41). The use of the attenuation law Att(Δ), for example, with synthetic estimates of Wood-Anderson seismometer response, during the Loma Preita earthquake, leads to estimates of M_L which agree with the surface wave and moment magnitudes, and which are essentially distance-independent.     

内蒙古中西部地区中小地震矩震级研究

基于S震相"S窗"内的波形信号识别、品质因子Q(f)和22个台站场地响应,利用2009~2016年3月内蒙古中西部地区地震的波形资料,反演了182次中小地震的震源波谱参数,得到这些小震的零频幅值及其拐角频率,据此计算了这些地震的地震矩M_0、矩震级M_W和应力降Δσ。利用回归分析方法得到了近震震级与矩震级、矩震级与应力降的关系式。分析表明,近震震级与矩震级、矩震级与应力降呈线性关系。可见,将矩震级纳入地震的快报与正式目录中,可以丰富地震观测报告内容,更好地为地震应急和地震科研服务。     

Experimental evidence for flexibility of a building foundation supported by concrete friction piles

This article explores the possibility to measure deformations of building foundations from measurements of ambient noise and strong motion recordings. The case under study is a seven-storey hotel building in Van Nuys, California. It has been instrumented by strong motion accelerographs, and has recorded several earthquakes, including the 1971 San Fernando (M_L=6.6, R=22 km), 1987 Whittier–Narrows (M_L=5.9, R=41 km), 1992 Landers (M_L=7.5, R=186 km), 1992 Big Bear (M_L=6.5, R=149 km), and 1994 Northridge (M_L=6.4, R=1.5 km) earthquake and its aftershocks (20 March: M_L=5.2, R=1.2 km; 6 December, 1994: M_L=4.3, R=11 km). It suffered minor structural damage in 1971 earthquake and extensive damage in 1994. Two detailed ambient vibration tests were performed following the Northridge earthquake, one before and the other one after the 20 March aftershock. These included measurements at a grid of points on the ground floor and in the parking lot surrounding the building, presented and analyzed in this article. The analysis shows that the foundation system, consisting of grade beams on friction piles, does not act as a “rigid body” but deforms during the passage of microtremor and therefore earthquake waves. For this geometrically and by design essentially symmetric building, the center of stiffness of the foundation system appears to have large eccentricity (this is seen both from the microtremor measurements and from the earthquake recordings). This eccentricity may have contributed to strong coupling of transverse and torsional responses, and to larger than expected torsional response, contributing to damage during the 1994 Northridge, earthquake.     

Calibration of magnitude scales for earthquakes of the Mediterranean

In order to provide the tools for uniform size determination for Mediterranean earthquakes over the last 50-year period of instrumental seismology, we have regressed the magnitude determinations for 220 earthquakes of the European-Mediterranean region over the 1977–1991 period, reported by three international centres, 11 national and regional networks and 101 individual stations and observatories, using seismic moments from the Harvard CMTs. We calibrate M(M₀) regression curves for the magnitude scales commonly used for Mediterranean earthquakes (M_L, M_WA, m_b, M_S, M_LH, M_LV, M_D, M); we also calibrate static corrections or specific regressions for individual observatories and we verify the reliability of the reports of different organizations and observatories. Our analysis shows that the teleseismic magnitudes (m_b, M_S) computed by international centers (ISC, NEIC) provide good measures of earthquake size, with low standard deviations (0.17–0.23), allowing one to regress stable regional calibrations with respect to the seismic moment and to correct systematic biases such as the hypocentral depth for M_S and the radiation pattern for m_b; while m_b is commonly reputed to be an inadequate measure of earthquake size, we find that the ISC m_b is still today the most precise measure to use to regress M_W and M₀ for earthquakes of the European-Mediterranean region; few individual observatories report teleseismic magnitudes requiring specific dynamic calibrations (BJI, MOS). Regional surface-wave magnitudes (M_LV, M_LH) reported in Eastern Europe generally provide reliable measures of earthquake size, with standard deviations often in the 0.25–0.35 range; the introduction of a small (±0.1–0.2) static station correction is sometimes required. While the Richter magnitude M_L is the measure of earthquake size most commonly reported in the press whenever an earthquake strikes, we find that M_L has not been computed in the European-Mediterranean in the last 15 years; the reported local magnitudes M_WA and M_L do not conform to the Richter formula and are of poor quality and little use, with few exceptions requiring ad hoc calibrations similar to the M_S regression (EMSC, ATH). The duration magnitude M_D used by most seismic networks confirms that its use requires accurate station calibrations and should be restricted only to events with low seismic moments.     

岷县漳县6.6级地震及余震强震动记录特征分析

对2013年7月22日甘肃岷县—漳县地震获取的273条主余震加速度记录进行格式转换、基线校正和滤波等常规处理,分析该地震主震(M_S6.6)和余震(M_S5.6)两次地震记录的幅值、持时以及反应谱特征,发现M_S6.6主震记录的PGA范围在0.728～177.5 gal间,M_S5.6余震记录的PGA范围在0.732～69.3 gal间;将观测数据与霍俊荣和第五代《中国地震动参数区划图》地震动衰减关系进行对比,发现霍俊荣衰减关系更吻合于此次地震的主余震加速度衰减;绘制主余震5%～95%重要持时分布图,并针对主震62MXT反应谱和本地设计谱以及近些年国内主要强震震中反应谱开展比较分析;最后研究土层台、基岩台、相同台站各个震级反应谱的特征。     

Determination of source duration and moment of earthquakes in the Indian Peninsula using love waveforms

Love waves recorded by a long-period seismograph at New Delhi (NDI) from seven earthquakes of magnitude 4.3 to 5.2 in Koyna and Bhatsa on the western coast and one earthquake in Ongole on the eastern coast of the Indian Peninsula have been used to determine the seismic moment for each of the earthquakes by waveform modeling. Transverse component of the synthetic seismogram shows that the maximum amplitude of waveform decreases with an increase of source duration. Thus for an evaluation of the seismic moment by equating the amplitude level of the observed and synthetic waveforms, we must know the source duration. The synthetic seismogram also indicates that a short source duration gives rise to a small but sharp pulse and this pulse is interpreted as anLg wave. Comparison of the observed and synthetic waveforms has been used for a simultaneous evaluation of the source duration and seismic moment. The source durations are found to vary between 2.2 and 4.4 s; for earthquakes with a magnitude range between 4.3 and 5.2 these durations are slightly higher than normal. We obtain moment (M ₀) of Ongole earthquake (M _L =5.1)as 1.7×10²⁴ dyne-cm; moments of Koyna and Bhatsa earthquakes (4.3M _L 5.2) on the western coast lie between 0.7×10²³ and 3.6×10²³ dyne-cm. Moment (M ₀)-magnitude (M _L ) relation logM ₀=1.5M _L +16.0 for the western United States region agrees as well, in general, with the results for the earthquakes in the Indian Peninsula.     

青海祁连MS5.2地震微震检测与目录完备性研究

利用模板匹配方法对2015年11月23日青海省祁连县M_S5.2地震进行遗漏地震检测研究,由于主震后短时间内目录中遗漏事件较多,故对主震后1天的连续波形进行检测。主震后1天内青海测震台网记录到的余震个数(包括单台)共62个,选取主震后M_L1.0以上余震30个作为模板事件,通过匹配滤波的方式扫描出遗漏地震31个,约为台网目录给出的0.5倍。基于包络差峰值振幅与震级的线性关系估测检测事件的震级参数,最后将检测后的余震目录与台网余震目录在主震后1天内的最小完备震级进行对比分析,结果发现检测后最小完备震级从M_L1.2降到了M_L0.7,得到青海测震台网在祁连地区最小完整性震级为M_L0.7。     

Modelling of Strong Ground Motions from 1991 Uttarkashi,India, Earthquake Using a Hybrid Technique

We present a simple and efficient hybrid technique for simulating earthquake strong ground motion. This procedure is the combination of the techniques of envelope function (Midorikawa et al. Tectonophysics 218:287–295, 1993) and composite source model (Zeng et al. Geophys Res Lett 21:725–728, 1994). The first step of the technique is based on the construction of the envelope function of the large earthquake by superposition of envelope functions for smaller earthquakes. The smaller earthquakes (sub-events) of varying sizes are distributed randomly, instead of uniform distribution of same size sub-events, on the fault plane. The accelerogram of large event is then obtained by combining the envelope function with a band-limited white noise. The low-cut frequency of the band-limited white noise is chosen to correspond to the corner frequency for the target earthquake magnitude and the high-cut to the Boore’s f _max or a desired frequency for the simulation. Below the low-cut frequency, the fall-off slope is 2 in accordance with the ω² earthquake source model. The technique requires the parameters such as fault area, orientation of the fault, hypocenter, size of the sub-events, stress drop, rupture velocity, duration, source–site distance and attenuation parameter. The fidelity of the technique has been demonstrated by successful modeling of the 1991 Uttarkashi, Himalaya earthquake (Ms 7). The acceptable locations of the sub-events on the fault plane have been determined using a genetic algorithm. The main characteristics of the simulated accelerograms, comprised of the duration of strong ground shaking, peak ground acceleration and Fourier and response spectra, are, in general, in good agreement with those observed at most of the sites. At some of the sites the simulated accelerograms differ from observed ones by a factor of 2–3. The local site geology and topography may cause such a difference, as these effects have not been considered in the present technique. The advantage of the technique lies in the fact that detailed parameters such as velocity-Q structures and empirical Green’s functions are not required or the records of the actual time history from the past earthquakes are not available. This method may find its application in preparing a wide range of scenarios based on simulation. This provides information that is complementary to the information available in probabilistic hazard maps.     

Moment Magnitude (M W) and Local Magnitude (M L) Relationship for Earthquakes in Northeast India

An attempt has been made to examine an empirical relationship between moment magnitude (M _W) and local magnitude (M _L) for the earthquakes in the northeast Indian region. Some 364 earthquakes that were recorded during 1950–2009 are used in this study. Focal mechanism solutions of these earthquakes include 189 Harvard-CMT solutions (M _W?≥?4.0) for the period 1976–2009, 61 published solutions and 114 solutions obtained for the local earthquakes (2.0?≤?M _L?≤?5.0) recorded by a 27-station permanent broadband network during 2001–2009 in the region. The M _W–M _L relationships in seven selected zones of the region are determined by linear regression analysis. A significant variation in the M _W–M _L relationship and its zone specific dependence are reported here. It is found that M _W is equivalent to M _L with an average uncertainty of about 0.13 magnitude units. A single relationship is, however, not adequate to scale the entire northeast Indian region because of heterogeneous geologic and geotectonic environments where earthquakes occur due to collisions, subduction and complex intra-plate tectonics.     

2020年和林格尔ML4.5地震微震匹配定位及发震构造探讨

内蒙古和林格尔地处鄂尔多斯块体北缘阴山地震带内,历史上6级以上强震频发.2020年3月30日和林格尔发生ML4.5地震,打破了自2005年以来阴山地震带ML4.0以上地震的长期平静.研究此次地震序列的发震构造对区域应力状态和地震危险性分析有重要作用,然而内蒙古地震台网台站较为稀疏,相对于华北其他地区地震监测能力较低,对...     

Evaluation of models for Fourier amplitude spectra for the Taiwan region

In the Taiwan region, the empirical spectral models for estimating ground-motion parameters were obtained recently on the basis of recordings of small to moderate (5.0≤M_L≤6.5) earthquakes. A large collection of acceleration records from the M_L=7.3 Chi-Chi earthquake (21 September, 1999) makes it possible to test the applicability of the established relationships in the case of larger events. The comparison of ground-motion parameters (Fourier amplitude spectra, peak accelerations and response spectra), which were calculated using the models, and the observed data demonstrates that the models could provide an accurate prediction for the case of the Chi-Chi earthquake and the largest aftershocks. However, there are some peculiarities in the ground-motion frequency content and attenuation that, most probably, are caused by the features of the rupture process of the large shallow earthquake source.     

Critical evaluation of strong motion in Kocaeli and Düzce (Turkey) earthquakes

Turkey was struck by two major events on August 17th and November 12th, 1999. Named Kocaeli (M_w=7.4) and Düzce (M_w=7.2) earthquakes, respectively, the two earthquakes provided the most extensive strong ground motion data set ever recorded in Turkey. The strong motion stations operated by the General Directorate of Disaster Affairs, the Kandilli Observatory and Earthquake Research Institute of Bogazici University and Istanbul Technical University have produced at least 27 strong motion records for the Kocaeli earthquake within 200 km of the fault. Kocaeli earthquake has generated six motions within 20 km of the fault adding significantly to the near-field database of ground motions for M_w>=7.0 strike–slip earthquakes. The paper discusses available strong motion data, studies their attenuation characteristics, analyses time domain, as well as spectral properties such as spectral accelerations with special emphasis on fault normal and fault parallel components and the elastic attenuation parameter, kappa. A simulation of the Kocaeli earthquake using code FINSIM is also presented.     

Estimation of strong ground motions from hypothetical earthquakes on the Cascadia subduction zone,Pacific Northwest

Strong ground motions are estimated for the Pacific Northwest assuming that large shallow earthquakes, similar to those experienced in southern Chile, southwestern Japan, and Colombia, may also occur on the Cascadia subduction zone. Fifty-six strong motion recordings for twenty-five subduction earthquakes ofM _s7.0 are used to estimate the response spectra that may result from earthquakesM _w<81/4. Large variations in observed ground motion levels are noted for a given site distance and earthquake magnitude. When compared with motions that have been observed in the western United States, large subduction zone earthquakes produce relatively large ground motions at surprisingly large distances. An earthquake similar to the 22 May 1960 Chilean earthquake (M _w 9.5) is the largest event that is considered to be plausible for the Cascadia subduction zone. This event has a moment which is two orders of magnitude larger than the largest earthquake for which we have strong motion records. The empirical Green's function technique is used to synthesize strong ground motions for such giant earthquakes. Observed teleseismicP-waveforms from giant earthquakes are also modeled using the empirical Green's function technique in order to constrain model parameters. The teleseismic modeling in the period range of 1.0 to 50 sec strongly suggests that fewer Green's functions should be randomly summed than is required to match the long-period moments of giant earthquakes. It appears that a large portion of the moment associated with giant earthquakes occurs at very long periods that are outside the frequency band of interest for strong ground motions. Nevertheless, the occurrence of a giant earthquake in the Pacific Northwest may produce quite strong shaking over a very large region.     

2014年金寨ML3.9震群序列遗漏地震检测及发震构造分析

针对2014年8月-2015年1月安徽金寨发生的ML3.9震群,利用匹配滤波技术补充台网目录遗漏的地震事件,再利用波形互相关震相检测技术标定P波和S波到时,进一步采用双差定位方法对震群进行重定位,结合震源机制解等分析此次震群活动可能的发震构造.计算结果显示,通过互相关扫描检测到1376个地震台网常规分析遗漏的地震,数量...     

冀鲁豫交界地区重力场变化与区内小震关系探讨

冀鲁豫交界地区历史上曾发生过多次中强地震,是一直受关注的地震危险区。利用该区1995—2010年共20期的流动重力观测资料,在统一起算基准、消除系统误差的基础上进行拟稳平差,计算重力场累积变化量;通过对重力场变化的分析,研究测区内重力场动态变化与测区内发生的一系列4级左右地震的关系。研究结果认为:冀鲁豫地区发生M_L4.0以上地震前后震区的重力场变化具有明显的"震前持续上升-震后反向恢复"的特征。     

Self-organized Fractal Seismicity of Reservoir Triggered Earthquakes in the Koyna-Warna Seismic Zone,Western India

— Analysis of the Koyna-Warna earthquake catalog (1968–1996) shows that on an average there is a positive correlation between the b value (decrease) and fractal dimensions (decrease in both D₂^s and D₂^t) of earthquake epicenters 0.5 and 2.5 years prior to 1973 (M5.2) and 1980 (M5.5) events, respectively, except a negative correlation for about five years (1988–1993) prior to the 1993/1994 sequence (M5.4). This positive correlation indicates a weaker clustering, or that the epicenters tend to fill the two-dimensional plane. While the origin of the negative correlation seems to be that during periods of large events (low b value), there is strong clustering around the main shock epicenter (high fractal dimension). Interestingly, during the last year (1995–1996) of the studied period both the b value and correlation dimensions rose significantly, suggesting that stress release occurs through increased levels of low magnitude and increasingly scattered seismicity, suggesting an increased risk of larger magnitude events. Incidentally, during 2000 three earthquakes of magnitude M 5.0, one earthquake of M 4.0, 45 earthquakes of magnitude M 3.0–3.9, and several thousand earthquakes of M < 3 have occurred in the region. Thus it can be inferred that at local scales the relationship yields both positive and negative correlation that appears to be controlled by different modes of failure within the active fault complex.Acknowledgement. The authors are grateful to Dr. B.K. Rastogi of NGRI for providing the catalog of Koyna earthquakes and for useful scientific discussions. The comments of Dr. I. G. Main have improved the quality of paper for which we extend to him our sincere thanks. One of the authors (AOM) thanks the Third World Academy of Science and the Council of Scientific and Industrial Research, India for the Postdoctoral Fellowship award under which this work was carried out.     

Correlating Precursory Declines in Groundwater Radon with Earthquake Magnitude

Both studies at the Antung hot spring in eastern Taiwan and at the Paihe spring in southern Taiwan confirm that groundwater radon can be a consistent tracer for strain changes in the crust preceding an earthquake when observed in a low‐porosity fractured aquifer surrounded by a ductile formation. Recurrent anomalous declines in groundwater radon were observed at the Antung D1 monitoring well in eastern Taiwan prior to the five earthquakes of magnitude (M_w): 6.8, 6.1, 5.9, 5.4, and 5.0 that occurred on December 10, 2003; April 1, 2006; April 15, 2006; February 17, 2008; and July 12, 2011, respectively. For earthquakes occurring on the longitudinal valley fault in eastern Taiwan, the observed radon minima decrease as the earthquake magnitude increases. The above correlation has been proven to be useful for early warning local large earthquakes. In southern Taiwan, radon anomalous declines prior to the 2010 M_w 6.3 Jiasian, 2012 M_w 5.9 Wutai, and 2012 M_L 5.4 Kaohsiung earthquakes were also recorded at the Paihe spring. For earthquakes occurring on different faults in southern Taiwan, the correlation between the observed radon minima and the earthquake magnitude is not yet possible.