Homogenization of Earthquake Catalog for Northeast India and Adjoining Region

A catalog for northeast India and the adjoining region for the period 1897–2009 with 4,497 earthquakes events is compiled for homogenization to moment magnitude M _w,GCMT in the magnitude range 3–8.7. Relations for conversion of m _b and M _s magnitudes to M _w,GCMT are derived using three different methods, namely, linear standard regression, inverted standard regression (ISR) and orthogonal standard regression (OSR), for different magnitude ranges based on events data for the catalog period 1976–2006. The OSR relations for M _s to M _w,GCMT conversion derived in this paper have significantly lower errors in regression parameters compared to the relations reported in other studies. Since the number of events with magnitude ≥7 for this region is scanty, we, therefore, considered whole India region to obtain the regression relationships between M _w,GCMT and M _s,ISC. A relationship between M _w,GCMT and M _w,NEIC is also obtained based on 17 events for the range 5.2 ≤ magnitude ≤ 6.6. A unified homogeneous catalog prepared using the conversion relations derived in this paper can serve as a reference catalog for seismic hazard assessment studies in northeast India and the adjoining region.     

The Conversion Relationship between the Local Earthquake Magnitude and Surface Wave Magnitude in the Inner Mongolia Digital Seismic Network

Using 116 earthquakes over M_L3.8 in the Inner Mongolia region from 2008 to 2015, the local earthquake magnitude M_L and surface wave magnitude M_S are remeasured. Based on norm linear regression(SR1 and SR2) and norm(OR) orthogonal regression method, we established the conversion relationship between M_L and M_S. The results were tested with Gaussian disturbance. The result shows that the orthogonal regression method(OR) result has the best fitting curve, and the conversion relation is M_S=0.96 M_L-0.10. The difference between our result and Guo Lücan's(M_S=1.13 M_L-1.08) may be caused by regional tectonic characteristics. M_(S Inner Mongolia) value is significantly higher than the M_(S empirical) value, with an average difference of 0.23, the difference distribution of empirical relation and the rectified relation is in the range of 0.2-0.3.     

The empirical relationships between M s, m b and M L for China and vicinity

Data from 753 earthquakes are used to determine a relationship between surface-wave magnitude (M _s) and bodywave magnitude (m _b), and from 541 earthquakes to determine a relationship between surface-wave magnitude (M _s) and local magnitude (M _L) for China and vicinity: M _s=0.9883 m _b-0.0420, M _s=0.9919 M _L-0.1773. The relationship of M _s versus m _b is obtained for 292 events occurred in the Chinese mainland in the time period from 1964 to 1996, 291 events occurred in Taiwan in the time period from 1964 to 1995 and 170 events occurred in the surrounding area. Standard deviation of the fitting is 0.445. Relationship of M _s versus M _L is obtained for 36 events occurred in the Chinese mainland, 293 events occurred in Taiwan, China and 212 events occurred in the surrounding area. The total amount is 541 events. Standard deviation of the fitting is 0.4673. The uncertainties of the converted M _s in different magnitude intervals can be estimated using complementary cumulative distribution function (CCDF). In the relationship of M _s versus m _b, taking ±0.25 as a range of uncertainties, in magnitude interval m _b 4.0–4.9, the probabilities for the converted M _s taken value less than (M _s-0.25) and more than (M _s+0.25) are 17% and 27% respectively. Similarly, we have probabilities for m _b 5.0–5.9 are 34% and 20% and that for m _b 6.0–6.9 are 11% and 47%. In the relationship of M _s versus M _L, if the range of uncertainties is still taken as ±0.25, the corresponding probabilities for magnitude interval M _L 4.0–4.9 are 22% and 38%, for M _L 5.0–5.9 are 20% and 15% and for magnitude interval M _L 6.0–6.9, are 15% and 29%, respectively. The relationships developed in this paper can be used for the conversion of one magnitude scale into another magnitude scales conveniently. The estimation of uncertainties described in this paper is more accurate and more objective than the usual estimation expressed by deviation. The estimations described in this paper indicate various dispersions in different magnitude intervals of original data. The estimations of uncertainties described by probabilities can be well connected with the total estimations of uncertainties in seismic hazard assessment.     

The instability of the <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> and <Emphasis Type="Italic">M</Emphasis><Subscript>L</Subscript> comparison for earthquakes in Greece for the period 1969 to 2007

We use 576 earthquakes of magnitude, M _w, 3.3 to 6.8 that occurred within the region 33° N–42.5° N, 19° E–30° E in the time period 1969 to 2007 to investigate the stability of the relation between moment magnitude, M _w, and local magnitude, M _L, for earthquakes in Greece and the surrounding regions. We compare M _w to M _L as reported in the monthly bulletins of the National Observatory of Athens (NOA) and to M _L as reported in the bulletins of the Seismological Station of the Aristotle University of Thessaloniki. All earthquakes have been analyzed through regional or teleseismic waveform inversion, to obtain M _w, and have measured maximum trace amplitudes on the Wood–Anderson seismograph in Athens, which has been in operation since 1964. We show that the Athens Wood–Anderson seismograph performance has changed through time, affecting the computed by NOA M _L by at least 0.1 magnitude units. Specifically, since the beginning of 1996, its east–west component has been recording systematically much larger amplitudes compared to the north–south component. From the comparison between M _w and M _L reported by Thessaloniki, we also show that the performance of the sensors has changed several times through time, affecting the calculated M _L’s. We propose scaling relations to convert the M _L values reported from the two centers to M _w. The procedures followed here can be applied to other regions as well to examine the stability of magnitude calculations through time.     

Empirical model for estimating Fourier amplitude spectra of ground acceleration in Taiwan region

A collection of ground‐motion recordings (1070 acceleration records) of moderate (5.1⩽M_L⩽6.5) earthquakes obtained during the execution of the Taiwan Strong Motion Instrumentation Program (TSMIP) since 1991 was used to study source scaling model and attenuation relations for a wide range of earthquake magnitudes and distances and to verify the models developed recently for the Taiwan region. The results of the analysis reveal that the acceleration spectra of the most significant part of the records, starting from S‐wave arrival, can be modelled accurately using the Brune's ω‐squared source model with magnitude‐dependent stress parameter Δσ, that should be determined using the recently proposed regional relationships between magnitude (M_L) and seismic moment (M₀) and between M₀ and Δσ. The anelastic attenuation Q of spectral amplitudes with distance may be described as Q=225 ƒ^1.1 both for deep (depth more than 35 km) and shallow earthquakes. The source scaling and attenuation models allow a satisfactory prediction of the peak ground acceleration for magnitudes 5.1⩽M⩽6.5 and distances up to about 200 km in the Taiwan region, and may be useful for seismic hazard assessment. Copyright © 2000 John Wiley & Sons, Ltd.     

Seismic source-parameter relations for earthquakes in Greece

A data set of nineteen, mainly shallow, moderate to large earthquakes, which occurred in the Aegean and the surrounding area, has been used to derive empirical relations for kinematic fault parameters. Thus the relations between seismic momentM ₀ and magnitudeM _s andm _b and betweenM ₀ andM _s and fault dimensionsS andL have been determined. From these relations and theoretical ones it was deduced that earthquakes in the Aegean and the surrounding events, chiefly interplate, are characterized by low average stress drop values. Values of ranging from 1 to 30 bar are consistent with the data. It was also found that, in general terms, most of the data obey the geometrical similarity conditionL=2w, whereL is the fault length measured along the strike andw is its width measured along the dip. For strike-slip faults, however, the conditionL=4w seems to hold.     

Empirical Global Relations Converting M S and m b to Moment Magnitude

The existence of several magnitude scales used by seismological centers all over the world and the compilation of earthquake catalogs by many authors have rendered globally valid relations connecting magnitude scales a necessity. This would allow the creation of a homogeneous global earthquake catalog, a useful tool for earthquake research. Of special interest is the definition of global relations converting different magnitude scales to the most reliable and useful scale of magnitude, the moment magnitude, M _W. In order to accomplish this, a very large sample of data from international seismological sources (ISC, NEIC, HRVD, etc.) has been collected and processed. The magnitude scales tested against M _W are the surface wave magnitude, M _S, the body wave magnitude, m _b, and the local magnitude, M _L. The moment magnitudes adopted have been taken from the CMT solutions of HRVD and USGS. The data set used in this study contains 20,407 earthquakes, which occurred all over the world during the time period 1.1.1976–31.5.2003, for which moment magnitudes are available. It is shown that well-defined relations hold between M _W and m _b and M _S and that these relations can be reliably used for compiling homogeneous, with respect to magnitude, earthquake catalogs.     

Scaling relations of moment magnitude,local magnitude,and duration magnitude for earthquakes originated in northeast India

In this study, we aim to improve the scaling between the moment magnitude (M _W), local magnitude (M _L), and the duration magnitude (M _D) for 162 earthquakes in Shillong-Mikir plateau and its adjoining region of northeast India by extending the M _W estimates to lower magnitude earthquakes using spectral analysis of P-waves from vertical component seismograms. The M _W-M _L and M _W-M _D relationships are determined by linear regression analysis. It is found that, M _W values can be considered consistent with M _L and M _D, within 0.1 and 0.2 magnitude units respectively, in 90 % of the cases. The scaling relationships investigated comply well with similar relationships in other regions in the world and in other seismogenic areas in the northeast India region.     

Comparison between different earthquake magnitudes determined by China Seismograph Network

By linear regression and orthogonal regression methods, comparisons are made between different magnitudes (lo-cal magnitude ML, surface wave magnitudes MS and MS7, long-period body wave magnitude mB and short-period body wave magnitude mb) determined by Institute of Geophysics, China Earthquake Administration, on the basis of observation data collected by China Seismograph Network between 1983 and 2004. Empirical relations between different magnitudes have been obtained. The result shows that: 1 As different magnitude scales reflect radiated energy by seismic waves within different periods, earthquake magnitudes can be described more objectively by using different scales for earthquakes of different magnitudes. When the epicentral distance is less than 1 000 km, local magnitude ML can be a preferable scale; In case M<4.5, there is little difference between the magnitude scales; In case 4.5MS, i.e., MS underestimates magnitudes of such events, therefore, mB can be a better choice; In case M>6.0, MS>mB>mb, both mB and mb underestimate the magnitudes, so MS is a preferable scale for deter-mining magnitudes of such events (6.08.5, a saturation phenomenon appears in MS, which cannot give an accurate reflection of the magnitudes of such large events; 2 In China, when the epicentral distance is less than 1 000 km, there is almost no difference between ML and MS, and thus there is no need to convert be-tween the two magnitudes in practice; 3 Although MS and MS7 are both surface wave magnitudes, MS is in general greater than MS7 by 0.2~0.3 magnitude, because different instruments and calculation formulae are used; 4 mB is almost equal to mb for earthquakes around mB4.0, but mB is larger than mb for those of mB≥4.5, because the periods of seismic waves used for measuring mB and mb are different though the calculation formulae are the same.     

Retardations in fault creep rates before local moderate earthquakes along the San Andreas fault system,central California

Records of shallow aseismic slip (fault creep) obtained along parts of the San Andreas and Calaveras faults in central California demonstrate that significant changes in creep rates often have been associated with local moderate earthquakes. An immediate postearthquake increase followed by gradual, long-term decay back to a previous background rate is generally the most obvious earthquake effect on fault creep. This phenomenon, identified as aseismic afterslip, usually is characterized by above-average creep rates for several months to a few years. In several cases, minor step-like movements, called coseismic slip events, have occurred at or near the times of mainshocks. One extreme case of coseismic slip, recorded at Cienega Winery on the San Andreas fault 17.5 km southeast of San Juan Bautista, consisted of 11 mm of sudden displacement coincident with earthquakes ofM _L =5.3 andM _L =5.2 that occurred 2.5 minutes apart on 9 April 1961. At least one of these shocks originated on the main fault beneath the winery. Creep activity subsequently stopped at the winery for 19 months, then gradually returned to a nearly steady rate slightly below the previous long-term average.The phenomena mentioned above can be explained in terms of simple models consisting of relatively weak material along shallow reaches of the fault responding to changes in load imposed by sudden slip within the underlying seismogenic zone. In addition to coseismic slip and afterslip phenomena, however, pre-earthquakeretardations in creep rates also have been observed. Onsets of significant, persistent decreases in creep rates have occurred at several sites 12 months or more before the times of moderate earthquakes. A 44-month retardation before the 1979M _L =5.9 Coyote Lake earthquake on the Calaveras fault was recorded at the Shore Road creepmeter site 10 km northwest of Hollister. Creep retardation on the San Andreas fault near San Juan Bautista has been evident in records from one creepmeter site for the past 5 years. Retardations with durations of 21 and 19 months also occurred at Shore Road before the 1974 and 1984 earthquakes ofM _L =5.2 andM _L =6.2, respectively.Although creep retardation remains poorly understood, several possible explanations have been discussed previously. (1) Certain onsets of apparent creep retardation may be explained as abrupt terminations of afterslip generated from previous moderate-mainshock sequences. (2) Retardations may be related to significant decreases in the rate of seismic and/or aseismic slip occurring within or beneath the underlying seismogenic zone. Such decreases may be caused by changes in local conditions related to growth of asperities, strain hardening, or dilatancy, or perhaps by passage of stress-waves or other fluctuations in driving stresses. (3) Finally, creep rates may be lowered (or increased) by stresses imposed on the fault by seismic or aseismic slip on neighboring faults. In addition to causing creep-rate increases or retardations, such fault interactions occasionally may trigger earthquakes.Regardless of the actual mechanisms involved and the current lack of understanding of creep retardation, it appears that shallow fault creep is sensitive to local and regional effects that promote or accompany intermediate-term preparation stages leading to moderate earthquakes. A strategy for more complete monitoring of fault creep, wherever it is known to occur, therefore should be assigned a higher priority in our continuing efforts to test various hypotheses concerning the mechanical relations between seismic and aseismic slip.     

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

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

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

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

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.     

The ground motion excited by the Olyutorskii earthquake of April 20, 2006 and by its aftershocks based on digital recordings

We studied broadband digital records of the M _W = 7.6 Olyutorskii earthquake of April 20, 2006 and its aftershocks at local and regional distances. We have made a detailed analysis of data on peak ground motion velocities and accelerations due to aftershocks based on records of two digital seismic stations, Tilichiki (TLC) and Kamenskoe (KAM). The first step in this analysis was to find the station correction for soil effects at TLC station using coda spectra. The correction was applied to the data to convert them to the reference bedrock beneath the Kamenskoe station. The second step involved multiple linear regression to derive average relationshis of peak amplitude to local magnitude M_L and distance R for the Koryak Upland conditions. The data scatter about the average relationshis is comparatively low (0.22–0.25 log units). The acceleration amplitudes for M _L = 5, R = 25 km are lower by factors of 2–3 compared with those for eastern Kamchatka, the western US, or Japan. A likely cause of this anomaly could be lower stress drops for the aftershocks.     

Quantification relations between the sourc eparameters

The various useful source-parameter relations between seismic moment and common use magnitude lg(M ₀) andM _s,M _L,m _b; between magnitudesMs andM _L,M _s andm _b,M _L andm _b; and between magnitudeM _s and lg(L) (fault length), lg (W) (fault width), lg(S) (fault area), lg(D) (average dislocation);M _L and lg(f _c) (corner frequency) have been derived from the scaling law which is based on an “average” two-dimensional faulting model of a rectangular fault. A set of source-parameters can be estimated from only one magnitude by using these relations. The average rupture velocity of the faultV _r=2.65 km/s, the total time of ruptureT(s)=0.35L (km) and the average dislocation slip rateD=11.4 m/s are also obtained. There are four strong points to measure earthquake size with the seismic moment magnitudeM _w.

1. The seismic moment magnitude shows the strain and rupture size. It is the best scale for the measurement of earthquake size.
2. It is a quantity of absolute mechanics, and has clear physical meaning. Any size of earthquake can be measured. There is no saturation. It can be used to quantify both shallow and deep earthquakes on the basis of the waves radiated.
3. It can link up the previous magnitude scales.
4. It is a uniform scale of measurement of earthquake size. It is suitable for statistics covering a broad range of magnitudes. So the seismic moment magnitude is a promising magnitude and worth popularization.

     

2017年3月渤海地震序列微震检测与发震构造分析

地震序列发震构造研究是区域地震活动性和地震危险性分析的重要基础.2017年3月渤海海域发生地震序列活动,该序列发生在郯城-庐江断裂带与张家口-渤海地震带的交汇部位,区域构造较为复杂.然而在渤海海域,连续运行的固定地震监测仪器难以布设,导致地震监测能力相对较弱.本文首先采用模板匹配方法对序列遗漏地震进行检测,再使用波形互...     

Precise determination of focal parameters for July 20, 1995 M_L=4.1 earthquake sequence in the Huailai basin

IntroductionSince the late 1970s, the quickly developed global digital seismograph network has been providing high quality recordings of large earthquakes in global scale, based on which digital seismology has made great progress. Compared with large earthquakes, moderate and small sized shocks have more frequent occurrence, and comprise clues to geological tectonics and tectonic stress field in a region. Preceding and following a large earthquake, usually occur numbers of small events that im…     

青海祁连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。     

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.     

Coal Mining Induced Seismicity in the Ruhr Area, Germany

Over the last 25 years mining-induced seismicity in the Ruhr area has continuously been monitored by the Ruhr-University Bochum. About 1,000 seismic events with local magnitudes between 0.7 ≤ M _L ≤ 3.3 are located every year. For example, 1,336 events were located in 2006. General characteristics of induced seismicity in the entire Ruhr area are spatial and temporal correlation with mining activity and a nearly constant energy release per unit time. This suggests that induced stresses are released rapidly by many small events. The magnitude–frequency distribution follows a Gutenberg–Richter relation which is a result from combining distributions of single longwalls that themselves show large variability. A high b-value of about 2 was found indicating a lack of large magnitude events. Local analyses of single longwalls indicate that various factors such as local geology and mine layout lead to significant differences in seismicity. Stress redistribution acts very locally since differences on a small scale of some hundreds of meters are observed. A regional relation between seismic moment M ₀ and local magnitude M _L was derived. The magnitude–frequency distribution of a single longwall in Hamm was studied in detail and shows a maximum at M _L = 1.4 corresponding to an estimated characteristic source area of about 2,200 m². Sandstone layers in the hanging or foot wall of the active longwall might fail in these characteristic events. Source mechanisms can mostly be explained by shear failure of two different types above and below the longwall. Fault plane solutions of typical events are consistent with steeply dipping fracture planes parallel to the longwall face and nearly vertical dislocation in direction towards the goaf. We also derive an empirical relation for the decay of ground velocity with epicenter distance and compare maximum observed ground velocity to local magnitude. This is of considerable public interest because about 30 events larger than M _L ≥ 1.2 are felt each month by people living in the mining regions. Our relations, for example, indicate that an event in Hamm with a peak ground velocity of 6 mm/s which corresponds to a local magnitude M _L between 1.7 and 2.3 is likely to be felt within about 2.3 km radius from the event.