Attenuation of modified mercalli intensity in New Zealand earthquakes

This paper describes a comprehensive study of the attenuation of Modified Mercalli intensity in New Zealand earthquakes, which has resulted in significantly different and new findings when compared with those of earlier work. The current study used recently revised magnitudes for the 30 events in the carefully selected data set. Magnitudes ranged from M_L = 5.0 to M_s = 7.8. Special effort was also put into establishing the depths of the events, which ranged from very shallow to 65 km. An expression of a form also used in peak ground acceleration attenuation studies was adopted, I = a + bM + cr + dlog₁₀r with r being taken as the mean distance from the centre of the fault rupture surface to each isoseismal. A two-step stratified regression analysis was used because it results in more realistic estimates of standard errors, although its mean curves were virtually the same as those derived from the standard one-step method. It was found possible to model accurately the attenuation of the deepest events from regressions of the data of the shallowest events, but the inverse was not true. In a study of the influence of source mechanisms it was found that the attenuation was the same for events with normal and strike-slip faulting. In contrast, events with reverse fault mechanisms were found to have higher intensities (for the same M and r) than normal and strike-slip events, by a factor corresponding to that found by Campbell in a study of peak ground accelerations, although the statistical test was narrowly short of the conventional significance level.     

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.     

The 2008 West Bohemia earthquake swarm in the light of the WEBNET network

A swarm of earthquakes of magnitudes up to M _L = 3.8 stroke the region of West Bohemia/Vogtland (border area between Czechia and Germany) in October 2008. It occurred in the Novy Kostel focal zone, where also all recent earthquake swarms (1985/1986, 1997, and 2000) took place, and was striking by a fast sequence of macroseismically observed earthquakes. We present the basic characteristics of this swarm based on the observations of a local network WEBNET (West Bohemia seismic network), which has been operated in the epicentral area, on the Czech territory. The swarm was recorded by 13 to 23 permanent and mobile WEBNET stations surrounding the swarm epicenters. In addition, a part of the swarm was also recorded by strong-motion accelerometers, which represent the first true accelerograms of the swarm earthquakes in the region. The peak ground acceleration reached 0.65 m/s². A comparison with previous earthquake swarms indicates that the total seismic moments released during the 1985/1986 and 2008 swarms are similar, of about 4E16 Nm, and that they represent the two largest swarms that occurred in the West Bohemia/ Vogtland region since the M _L = 5.0 swarm of 1908. Characteristic features of the 2008 swarm are its short duration (4 weeks) and rapidity and, consequently, the fastest seismic moment release compared to previous swarms. Up to 25,000 events in the magnitude range of 0.5 < M _L < 3.8 were detected using an automatic picker. A total of nine swarm phases can be distinguished in the swarm, five of them exceeding the magnitude level of 2.5. The magnitude–frequency distribution of the complete 2008 swarm activity shows a b value close to 1. The swarm hypocenters fall precisely on the same fault portion of the Novy Kostel focal zone that was activated by the 2000 swarm (M _L ≤ 3.2) in a depth interval from 6 to 11 km and also by the 1985/1986 swarm (M _L ≤ 4.6). The steeply dipping fault planes of the 2000 and 2008 swarms seem to be identical considering the location error of about 100 m. Furthermore, focal mechanisms of the 2008 swarm are identical with those of the 2000 swarm, both matching an average strike of 170° and dip of 80° of the activated fault segment. An overall upward migration of activity is observed with first events at the bottom and last events at the top of the of the activated fault patch. Similarities in the activated fault area and in the seismic moments released during the three largest recent swarms enable to estimate the seismic potential of the focal zone. If the whole segment of the fault plane was activated simultaneously, it would represent an earthquake of M _L ~5. This is in good agreement with the estimates of the maximum magnitudes of earthquakes that occurred in the West Bohemia/Vogtland region in the past.     

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.     

喀什基准台宽频带井下和地表地震计记录数据对比研究

对比分析喀什基准台麦盖提井下和相邻地表宽频带地震计记录的背景噪声、麦盖提ML3.1近震和棉兰老岛ML6.6远震波形记录以及功率谱,并进一步研究了相关的震级差和台站下方结构.研究结果表明,井下和地表地震计记录到的波形特征和功率谱既存在相似性,又具有差异性;井下地震计记录的近、远震波形更清晰、简洁,且功率谱也相对较低.麦盖...     

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.     

Revision of magnitudes of large shallow earthquakes, 1897–1912

In previous research, trace amplitudes of surface wave maxima recorded by undamped Milne seismographs were used to determine the surface-wave magnitudes M_s of large shallow earthquakes which occurred prior to 1912. For this purpose, the effective gain of these instruments was calibrated by using the surface-wave magnitudes M_s(GR) which were calculated from the unpublished worksheets for Seismicity of the Earth of Gutenberg and Richter. In this paper, the real quality of M_s(GR) is critically re-evaluated by using independent sets of data. It is found that M_s(GR) for the period 1904–1909 is considerably overestimated. The average excess from the real M_s is 0.5 units for 1904–1906, 0.4 for 1907, 0.3 for 1908–1909 and 0.0 for 1910–1912. This overestimation is so systematic and large that the previous results are all redetermined. The average effective gain of Milne instruments is revised to be 21.9; previously, the gain depended on M_s. This revision results in systematic reduction in the previously assigned magnitudes. The revised values of M_s for 264 shallow earthquakes, with M_s=6.8 and over in the period 1897–1912 inclusive, are listed. The present revision is large enough to preclude the possibility of the high activity of large shallow earthquakes around the turn of the century. The present results have a direct effect on all the magnitude catalogues of shallow earthquakes which occurred prior to 1909.     

Rediscovering Signal Complexity as a Teleseismic Discriminant

We re-examine the utility of teleseismic seismic complexity discriminants in a multivariate setting using United Kingdom array data. We measure a complexity discriminant taken on array beams by simply taking the logarithm of the ratio of the P-wave coda signal to that of the first arriving direct P wave (β_CF). The single station complexity discriminant shows marginal performance with shallow earthquakes having more complex signatures than those from explosions or deep earthquakes. Inclusion of secondary phases in the coda window can also degrade performance. However, performance improves markedly when two-station complexity discriminants are formed showing false alarm rates similar to those observed for network m_b − M_s. This suggests that multistation complexity discriminants may ameliorate some of the problems associated with m_b – M_s discrimination at lower magnitudes. Additionally, when complexity discriminants are combined with m_b – M_s there is a tendency for explosions, shallow earthquakes and deep earthquakes to form three distinct populations. Thus, complexity discriminants may follow a logic that is similar to m_b – M_s in terms of the separation of shallow earthquakes from nuclear explosions, although the underlying physics of the two discriminants is significantly different.     

Characteristics of strong ground motion for typical Australian intra-plate earthquakes and their relationship with the recommended response spectra

Global epicentre maps show that the majority of earthquakes are inter-plate, although moderate to large earthquakes do occur intra-plate, i.e. within the plates. The seismicity of the Australian continent is typical of intra-plate environments and a magnitude M_L 6 earthquake has an average return period of about 5 years. Recordings of Australian intra-plate earthquakes are investigated here to characterise their frequency content, peak acceleration and duration.Due to lack of quality strong motion records of large intra-plate earthquakes at short distances, synthetic seismograms are commonly used for testing structural behaviour. An empirical Green's Function method (Geophys. Res. Lett., 5 (1978), 1–4; Proceedings of the Third International Microzonation Conference, Seattle, USA, vol. 1, (1982), pp. 447–458.) is chosen to simulate a large earthquake by summation in time of a number of smaller earthquakes or sub-events, each given a slightly different origin time to represent more realistically the propagation of a rupture along an assumed fault plane. In the first instance, recordings on rock of the magnitude M_L 2.3 aftershock of the 29 December 1989 Newcastle earthquake were used as sub-events to simulate the main shock of magnitude M_L 5.6. Validation studies for events recorded elsewhere in Australia are also considered.The response spectra of such synthetic events will be compared with the recommended spectra developed empirically from a statistical analysis of strong motion data for magnitude 5.4–6.5 intra-plate earthquakes recorded in other parts of the world and normalised to a peak ground velocity of 50 mm/s which is typical for a return period of 500 years in Australia (Australasian Structural Engineering Conference, Auckland, New Zealand, (1998), pp. 439–444.). Preliminary results from this comparison with the response spectra recommended for the Building Code of Australia show that the synthetic waveforms produced by this method are realistic and can be used to represent ground motion during typical Australian intra-plate earthquakes.     

Local magnitude,surface wave magnitude and seismic energy

Summary The local magnitude M_L at the seismological station Pruhonice (PRU) was converted into surface wave magnitude M_S using the formula M_S=–3.2+1.45 M_L and the seismic wave energy was estimated using the relation log E (Joule)=1.2+2.0 M_L. It was proposed to apply the same conversion formulae at seismological stations Kaperské Hory (KHC) and Berggiesshübel (BRG) where the calibrating functions for local magnitudes were determined for the same set of earthquakes with common reference magnitudes as in the case of the PRU station.     

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.     

Calibration of Local Magnitude M<Subscript>L</Subscript> in the Azores Archipelago Based on Recent Digital Recordings

—The development of the digital seismic network in the Azores Archipelago during recent years made it possible to obtain the amplitudes (waveform) of recorded motion in a large set of stations. With this new data, maximum amplitudes of the Wood Anderson seismograph are computed, for each station/component, which, together with epicentral distances, allows for the estimation of local magnitude M _L . We used data recorded in 8 digital permanent three-component stations, with inter-stations distances up to 300 km, in the period June 1998 – June 2000, corresponding to a set of 1315 events with magnitude (M _L or M _D ) 2<M<5.8 and epicenters located in the Azores region, to estimate the coefficients of the equation to compute M _L , as well as to determine the corrections to be applied to each station. The new set of parameters, formed by attenuation coefficients and station corrections, were introduced in the calculations of the M _L , leading to smaller dispersions in the analyzed dataset. We also conclude that the attenuation in the first 150 km is similar to the California values, although higher for longer distances.     

A regional surface wave magnitude scale for the earthquakes of Russia’s Far East

The modified scale M _s(20R) is developed for the magnitude classification of the earthquakes of Russia’s Far East based on the surface wave amplitudes at regional distances. It extends the applicability of the classical Gutenberg scale M _s(20) towards small epicentral distances (0.7°–20°). The magnitude is determined from the amplitude of the signal that is preliminarily bandpassed to extract the components with periods close to 20 s. The amplitude is measured either for the surface waves or, at fairly short distances of 0.7°–3°, for the inseparable wave group of the surface and shear waves. The main difference of the M _s(20R) scale with the traditional M _s(BB) Soloviev–Vanek scale is its firm spectral anchoring. This approach practically eliminated the problem of the significant (up to–0.5) regional and station anomalies characteristic of the M _s(BB) scale in the conditions of the Far East. The absence of significant station and regional anomalies, as well as the strict spectral anchoring, make the M _s(20R) scale advantageous when used for prompt decision making in tsunami warnings for the coasts of Russia’s Far East.     

Synthesis of magnitude and focal mechanism computations for the M ≥ 4.5 earthquakes in France for the period 1995–2000

This paper proposes a synthesis of thestudies made in terms of source parametersevaluation for the last earthquakes oflocal magnitude greater than 4.5 whichoccurred in or nearby France during thelast five years. Focal mechanisms andseismic moments have been computed for thethree most important events, largely feltby the population: St Paul de Fenouillet(February 18^th 1996, M_L 5.6),Annecy-Epagny (July 15^th 1996, M_L5.2) and St-Béat (October 4^th1999, M_L 4.8). These focal mechanismshave been obtained either by regionalmoment tensor inversion or from firstmotion polarities and are compared withcomplementary studies made on theseearthquakes. In addition, for the otherearthquakes of local magnitude greater than4.5 which occurred nearby French borderssince the beginning of the recording ofbroadband data by the RéNaSS(Réseau National de SurveillanceSismique, French Seismological Survey) inmid–1995, several magnitude calculationsconcerning the following earthquakes arepresented: Pamplona (25/2/1996, M_L4.7), Aoste (31/03/1996, M_L 4.6),Imperia (24/2/1997, M_L 4.5),Barcelonette (31/10/1997, M_L 4.8),Pamplona (27/10/1998, M_L 4.9), andBonifacio (26/4/2000, M_L 4.5). Localmagnitudes are usually higher thanthe M_b magnitudes reported by the PDE(Preliminary Determination of Epicenters),while the extension of the Msz scale toregional magnitudes and the Mw magnitudesderived from seismic moments give smallervalues. The relative importance of thevarious earthquakes in terms of surfacewave magnitude or seismic moment does notalways agree with that implied by localmagnitudes.     

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.     

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

Simulating the Taipei basin response by numerical modeling of wave propagation

Taipei, the capital of Taiwan, suffered from destructive earthquakes four times during the 20th century (M _L = 7.3 on April 15, 1909; M _L = 6.8 on November 15, 1986; the Chi–Chi M _L = 7.3 earthquake on September 21, 1999; and M _L = 6.8 on March 31, 2002). Analysis of recorded data shows a strong dependence of spectral amplification in the Taipei Basin on earthquake depth and azimuth. At low frequencies (f < 3 Hz) significant larger amplifications are observed for shallow earthquakes as compared to intermediate depth events. The former ones also display strong azimuthal dependence. As structures with large response periods such as bridges and tall buildings are sensitive to these low frequencies the understanding of the associated wave effects within the basin and their role for site effect amplification is critical. The tool we employ is 3D finite-difference modeling of wave propagation of incident wave fronts. The available detailed model of the basin allows studying the wave effects. Modeling clearly reveals that basin edge effects as observed in data are related to surface wave generation at the basin edges with a high degree of azimuthal dependency. The reproduced site amplification effects are in qualitative agreement with the observations from strong motion data.     

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，地震目录最小完整性震级的减小有利于地震工作者对区域地震活动性作出更准确全面的结论，并使地震危险性分析更可靠。     

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

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