The reservoir-associated earthquakes of April 1983 in western Thailand: Source modeling and implications for induced seismicity

On 22 April 1983, a very large area of Thailand and part of Burma were strongly shaken by a rare earthquake (m _b=5.8,M _s=5.9). The epicenter was located at the Srinagarind reservoir about 190 km northwest of Bangkok, a relatively stable continental region that experienced little previous seismicity. The main shock was preceded by some foreshocks and followed by numerous aftershocks. The largest foreshock ofm _b=5.2 occurred 1 week before the main shock, and the largest aftershock ofm _b=5.3 took place about 3 hours after the main shock. Focal mechanisms of the three largest events in this earthquake sequence have been studied by other seismologists using firts-motion data. However, the solutions for the main shock and the largest aftershock showed significant inconsistency with known surface geology and regional tectonics. We reexamined the mechanisms of these three events by using teleseismicP-andS-waveforms and through careful readings ofP-wave first motions. The directions of theP axes in our study range from NNW-SSE to NNE-SSW, and nodal planes strike in the NW-SE to about E-W in agreement with regional tectonics and surface geology. The main shock mechanism strikes 255°, dips 48°, and slips 63.5°. The fault motions during the main shock and the foreshock are mainly thrust, while the largest aftershock has a large strike-slip component. The seismic moment and the stress drop of the mainshock are determined to be 3.86×10²⁴ dyne-cm and 180 bars, respectively. The occurrence of these thrust events appears to correlate with the unloading of the Srinagarind reservoir. The focal depths of the largest foreshock, the main shock, and the largest aftershock are determined to be 5.4 km, 8 km, and 22.7 km, respectively, from waveform modeling and relative location showing a downward migration of hypocenters of the three largest events during the earthquake sequence. Other characteristics of this reservoir-induced earthquake sequence are also discussed.     

汶川8.0级大震的余震最大强度用科里奥利力效应方法的判定

用科里奥利力效应预测强余震是一种震源物理的方法。回顾2008年汶川8.0级大震时用该方法判定余震最大强度的过程,半定性与综合判定为可能发生的最大余震强度为6.5级左右,实际发生了6.4级地震,与主震震级相差大于1级（M_主—M_余=1.6）。验证结果进一步说明该方法的科学性,给科里奥利力效应判定余震增加了一个可靠的判例。     

A catalog of aftershock sequences in Greece (1971–1997): Their spatial and temporal characteristics

A complete catalog of aftershock sequences is provided for main earthquakes with M_L 5.0, which occurred in the area of Greece and surrounding regions the last twenty-seven years. The Monthly Bulletins of the Institute of Geodynamics (National Observatory of Athens) have been used as data source. In order to get a homogeneous catalog, several selection criteria have been applied and hence a catalog of 44 aftershock sequences is compiled. The relations between the duration of the sequence, the number of aftershocks, the magnitude of the largest aftershock and its delay time from the main shock as well as the subsurface rupture length versus the magnitude of the main shock are calculated. The results show that linearity exists between the subsurface rupture length and the magnitude of the main shock independent of the slip type, as well as between the magnitude of the main shock (M) and its largest aftershock (Ma). The mean difference M–Ma is almost one unit. In the 40% of the analyzed sequences, the largest aftershock occurred within one day after the main shock.The fact that the aftershock sequences show the same behavior for earthquakes that occur in the same region supports the theory that the spatial and temporal characteristics are strongly related to the stress distribution of the fault area.     

The probability of occurrence of largest earthquakes in the European area

Summary TheGumbel's theory of largest values has been applied to the estimation of probability of occurrence and of return periods of largest earthquakes in the European area. For this study shallow shocks from the period 1901–1955 and from 15 earthquake zones were used. For each zone the largest magnitudes corresponding to one-year intervals were arranged in order of increasingM, grouped in classes and then the probabilitiesF(x _j) were calculated. The data plotted on the probability paper fit a straight line fairly well. The extrapolated lines yield the possibility of estimating large magnitudes which will be exceeded with a given probability, e.g. 1%. Such values were compared with largest magnitudes observed during the period 1901–1955. Their return periods indicate that in most regions the largest probable shock already occurred. Following the procedure ofEpstein-Lomnitz the coefficients and were calculated and compared with corresponding values ofa andb of the magnitude-frequency relation.     

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.     

Surface fault traces,fault plane solution and spatial distribution of the aftershocks of the September 13, 1986 earthquake of Kalamata (Southern Greece)

A shallow earthquake ofM _S=6.2 occurred in the southern part of the Peloponnesus, 12 km north of the port of the city of Kalamata, which caused considerable damage. The fault plane solution of the main shock, geological data and field observations, as well as the distribution of foci of aftershocks, indicate that the seismic fault is a listric normal one trending NNE-SSW and dipping to WNW. The surface ruptures caused by the earthquake coincide with the trace of a neotectonic fault, which is located 2–3 km east of the city of Kalamata and which is related to the formation of Messiniakos gulf during the Pliocene-Quaternary tectonics. Field observations indicate that the earthquake is due to the reactivation of the same fault.A three-days aftershock study in the area, with portable seismographs, recorded many aftershocks of which 39 withM _S1.7 were very well located. The distribution of aftershocks forms two clusters, one near the epicenter of the main shock in the northern part of the seismogenic volume, and the other near the epicenter of the largest aftershock (M _S=5.4) in the southern part of this volume. The central part of the area lacks aftershocks, which probably indicates that this is the part of the fault which slipped smoothly during the earthquake.     

Space-time correlations ofb values with stress release

A total of 1503 events for a 2-month period associated with am _N 2.6 rockburst is investigated for possible space-time correlations between low magnitude (–1.1 to –0.4)b values and several estimates of stress (static stress drop, apparent stress, and dynamic stress). Spatial variations of decreasingb values were found to be well correlated with increasing stress release estimates for time intervals prior to the rockburst and following the aftershock sequence. The strongest correlation tob value was with the dynamic stress drop, having correlation coefficients of 0.87 and 0.79 for the two intervals, respectively. The rockburst was found to actually occur at the intersection of the spatial coordinates corresponding to the largest gradient inb value. Based on these correlations, we conclude that the low magnitude seismicity is an indicator of the stress state within the rock mass, and can be used to study and forecast stress patterns in the vicinity of an impending major event. Time variations, however, did not show the same clear correlations and these are discussed in terms of departure from steady state conditions. Regardless, our results favour the use ofb values in a spatial, context rather than in a time analysis approach, and we consider thatb values provide valuable information regarding the changing stress conditions within the seismogenic volume.     

Effect of an earthquake occurrence on the flow intensity of the subsequent events

Aftershocks or swarms indicate increase of the flow intensity in the vicinity of the initial earthquakes. By normalizing their number according to the dynamic range of the standard frequency magnitude distribution the increase or positive aftereffect property of the initial earthquakes can be compared for different magnitude intervals, periods of time or regions. After applying accurate formal algorithm of aftershock identification it is possible to study negative aftereffect of the main events (nonaftershocks) in the catalog.Negative aftereffect means decrease of the probability of successive events in a time-space vicinity of the main event, when the aftershocks are over. The negative effect is the most important part of the seismic cycle and seismic gaps approach. Global statistical test give high confidence level for the relative decrease in intensity of the flow of the events withM7 in the first 20–25 years after the events withM8 in their 1^o-vicinities in the total time period under study of approximately 60 years. The decrease approximates 32% of the undisturbed intensity of the flow ofM>7 events in the vicinities.Self-similar negative aftereffect was observed 3–7 years after 6M<7 events, it totals approximately 18% of the undisturbed intensity. Another type of self-similarity of seismic regime, with respect to the negative aftereffect, is the decrease of probabilities of aftershocks with large magnitudes in aftershock sequences. When we have adequate dynamic range in the catalog for the study of this property, for example, for main events withM7 in the catalog with low cut-off limitM=4, the statistical significance of the negative aftereffect is clear. However, the absolute value of the effect is also rather small, about 10%, which means that in 90% of the cases the aftershock sequences do not experience lack of energy due to the main shock energy release and follow a standard magnitude distribution for earthquakes in the entire catalog.The small values of the negative aftereffect apparently indicate partial stress relase by earthquakes and may explain short recurrence time intervals after major earthquakes observed periodically in different places.     

2019年12月26日湖北应城M4.9地震震源参数及余震活动性研究

2019年12月26日湖北应城发生M4.9有感地震,其震感波及武汉大部分地区。为了分析该地震的发震构造及余震活动性,本文利用波形拟合方法测定了不同速度模型下该地震的震源机制解和矩心深度,并用Bootstrapping抽样反演技术评价反演结果;此外,利用模板匹配技术匹配主震和目录余震波形,获取了更为完整的余震目录。结果显示,应城地震以走滑为主,矩心深度7.5km左右,矩震级M_W4.67;应城地震有1个前震和17个余震,余震序列缺少M2~4事件,表明应城地震为孤立型地震,M2以下地震的b值为0.8。     

基于Gutenberg-Richter定律快速估算最大余震震级:以2017年九寨沟MS7.0地震为例

针对九寨沟M_S7.0地震之后不同时间段的余震序列目录,利用推定最大余震震级,给出了实际最大余震震级的估计值。结果表明,推定最大余震震级随主震后时间尺度的延长而趋于稳定,且该值与实际发生的最大余震的震级一致。需要强调的是,就九寨沟地震序列而言,当余震数据较为完备时,采用主震后较短时间段内（1～2天）的余震目录就可以较准确地估算出主震区域内可能发生的最大余震震级。实际上,主震后12h（0.5天）的余震数据已完全可以给出最大余震震级的有效下限。此外,计算中我们采用了里氏震级M_L和面波震级M_S的余震目录,结果显示,2种震级类型目录的估算结果完全一致,表明利用推定最大余震震级估算实际最大余震震级的方法不受震级类型的影响。据此,该最大余震震级快速评估方法可进一步推广应用于我国大陆地区中强震后强余震灾害分析评估中。目前的拟合技术也显示出随着测震技术的不断进步以及余震识别能力的提高,快速评估方法可以在主震后短时间（<1天）内准确地预测可能发生的最大余震震级。     

An application of the time- and magnitude-predictable model to long-term earthquake prediction in eastern Anatolia

In order to estimate the recurrence intervals for large earthquakes occurring in eastern Anatolia, this region enclosed within the coordinates of 36^∘–42^∘N, 35^∘–45^∘E has been separated into nine seismogenic sources on the basis of certain seismological and geomorphological criteria, and a regional time- and magnitude-predictable model has been applied for these sources. This model implies that the magnitude of the preceding main shock which is the largest earthquake during a seismic excitation in a seismogenic source governs the time of occurrence and the magnitude of the expected main shock in this source. The data belonging to both the instrumental period (M_S≥ 5.5) until 2003 and the historical period (I₀≥ 9.0 corresponding to M_S≥ 7.0) before 1900 have been used in the analysis. The interevent time between successive main shocks with magnitude equal to or larger than a certain minimum magnitude threshold were considered in each of the nine source regions within the study area. These interevent times as well as the magnitudes of the main shocks have been used to determine the following relations:

Properties of the Aftershock Sequences of the 2003 Bingöl, M D = 6.4, (Turkey) Earthquake

Aftershock sequences of the magnitude M _{W =}6.4 Bingöl earthquake of 1 May, 2003 (Turkey) are studied to analyze the spatial and temporal variability of seismicity parameters of the b value of the frequency-magnitude distribution and the p value describing the temporal decay rate of aftershocks. The catalog taken from the KOERI contains 516 events and one month’s time interval. The b value is found as 1.49 ± 0.07 with Mc =3.2. Considering the error limits, b value is very close to the maximum b value stated in the literature. This larger value may be caused by the paucity of the larger aftershocks with magnitude M _D ≥ 5.0. Also, the aftershock area is divided into four parts in order to detect the differences in b value and the changes illustrate the heterogeneity of the aftershock region. The p value is calculated as 0.86 ± 0.11, relatively small. This small p value may be a result of the slow decay rate of the aftershock activity and the small number of aftershocks. For the fitting of a suitable model and estimation of correct values of decay parameters, the sequence is also modeled as a background seismicty rate model. Constant background activity does not appear to be important during the first month of the Bingöl aftershock sequences and this result is coherent with an average estimation of pre-existing seismicity. The results show that usage of simple modified Omori law is reasonable for the analysis. The spatial variability in b value is between 1.2 and 1.8 and p value varies from 0.6 to 1.2. Although the physical interpretation of the spatial variability of these seismicity parameters is not straightforward, the variation of b and p values can be related to the stress and slip distribution after the mainshock, respectively. The lower b values are observed in the high stress regions and to a certain extent, the largest b values are related to Holocene alluvium. The larger p values are found in some part of the aftershock area although no slip occurred after the main shock and it is interpreted that this situation may be caused by the alluvium structure of the region. These results indicate that the spatial distribution in b and p values are generally related to the rupture mechanism and material properties of an aftershock area.     

Statistical features of aftershock distribution size for moderate and large earthquakes in Chinese mainland

Based on data of earthquake sequences with MS≥5.0 in Chinese mainland from 1970 to 2004,for different se-quence types and different rupture modes of the main shock,the relationship between aftershock distribution size Rand the magnitude of the main shock M0 has been studied statistically.Considering the rupture mode of the mainshock,we give the quantitative statistical relationships between R and M0 under 95%confidence level for differentsequence types.Qualitatively,lgR,the logarithm of the aftershock distribution size,is positively correlative to theM0,but the data distribution is dispersed.Viewing from different sequence types,the correlation between R and M0is very weak for isolated earthquake type(IET)sequence,R distributes in the range from 5 to 60 km;For main-shock-aftershock type(MAT),lgR is positively correlative to M0;For multiple main shock type(MMT),the core-lation between lgR and M0 is not very obvious when M0≤6.2 and R distributes in the range from 5 to 70 km,whileit shows a linear correlation when M0≥6.3.The statistical results also show that the occupational ratios of differentsequence types for strike-slip and oblique slip are almost the same.But for dip-slip(mostly are thrust mechanisms),the ratio of MAT is higher than that of IET and MMT.Comparing with previous results,it indicates that,when M0is large enough,R is mainly determined by M0 and there is almost no relationship with the rupture mode of themain shock.     

Relocation of the MS5.7 Earthquake Sequence in Songyuan, Jilin Province, 2018 and the Analysis of Its Seismogenic Structure

The 2018,Songyuan,Jilin M_S5. 7 earthquake occurred at the intersection of the FuyuZhaodong fault and the Second Songhua River fault. The moment magnitude of this earthquake is M_W5. 3,the centroid depth by the waveform fitting is 12 km,and it is a strike-slip type event. In this paper,with the seismic phase data provided by the China Earthquake Network, the double-difference location method is used to relocate the earthquake sequence,finally the relocation results of 60 earthquakes are obtained. The results show that the aftershock zone is about 4. 3km long and 3. 1km wide,which is distributed in the NE direction. The depth distribution of the seismic sequence is 9km-10 km. 1-2 days after the main shock,the aftershocks were scattered throughout the aftershock zone,and the largest aftershock occurred in the northeastern part of the aftershock zone. After 3-8 days,the aftershocks mainly occurred in the southwestern part of the aftershock zone. The profile distribution of the earthquake sequence shows that the fault plane dips to the southeast with the dip angle of about 75°. Combined with the regional tectonic setting,focal mechanism solution and intensity distribution,we conclude that the concealed fault of the Fuyu-Zhaodong fault is the seismogenic fault of the Songyuan M_S5. 7 earthquake. This paper also relocates the earthquake sequence of the previous magnitude 5. 0 earthquake in 2017. Combined with the results of the focal mechanism solution,we believe that the two earthquakes have the same seismogenic structure,and the earthquake sequence generally develops to the southwest. The historical seismic activity since 2009 shows that after the magnitude 5. 0 earthquake in 2017,the frequency and intensity of earthquakes in the earthquake zone are obviously enhanced,and attention should be paid to the development of seismic activity in the southwest direction of the earthquake zone.     

Statistical Analysis of Triggered Seismicity in the Kresna Region of SW Bulgaria (1904) and the Umbria-Marche Region of Central Italy (1997)

A version of the restricted trigger model is used to analyse the temporal behaviour of some aftershock sequences. The conditional intensity function of the model is similar to that of the Epidemic Type Aftershock-Sequence (ETAS) model with the restriction that only the aftershocks of magnitude bigger than or equal to some threshold M_tr can trigger secondary events. For this reason we have named the model Restricted Epidemic Type Aftershock-Sequence (RETAS) model. Varying the triggering threshold we examine the variants of the RETAS model which range from the Modified Omori Formula (MOF) to the ETAS model, including such models as limit cases. In this way we have a quite large set of models in which to seek the model that fits best an aftershock sequence bringing out the specific features of the seismotectonic region struck by the crisis. We have applied the RETAS model to the analysis of two aftershock sequences: The first is formed by the events which followed the strong earthquake of M=7.8 which occurred in Kresna, SW Bulgaria, in 1904. The second includes three main shocks and a large swarm of minor shocks following the quake of 26 September 1997 in the Umbria-Marche region, central Italy. The MOF provides the best fit to the sequence in Kresna; that leads to the thought that just the stress field changes due to the very strong main shock generate the whole sequence. On the contrary, the complex behaviour of the seismic sequence in Umbria-Marche appears when we make the threshold magnitude vary. Setting the cut-off magnitude M₀=2.9 the best fit is provided by the ETAS model, while if we raise the threshold magnitude M₀=3.6 and set M_tr=5.0, the RETAS model turns out to be the best model. In fact, observing the time distribution of this reduced data set, it appears more evident that especially the strong secondary events are followed by a cluster of aftershocks.     

Evidence for transform faulting in the Ionian sea: The Cephalonia island earthquake sequence of 1983

Accurate locations of aftershocks of the January 17, 1983 (M _s=7.0) main shock in the Ionian islands have been determined, as well as fault plane solutions for this main shock and its largest aftershock, which are interpreted as a right-lateral, strike-slip motion with a thrust component, on a fault striking in about a NE-SW direction.This is considered as a transform fault in the northwesternmost part of the Hellenic arc.     

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

基于指数函数的川滇地区余震序列衰减规律研究

较大的余震可能造成额外损失并有二次触发建筑物受灾的风险。为研究余震序列衰减规律,文章尝试采用指数衰减模型拟合分析5个不同地区余震序列,并借助修正赤池信息准则、贝叶斯信息准则与调整后R²,分析其与传统余震衰减模型的性能。结果表明,指数模型描述余震序列衰减规律的能力与修正的大森余震模型、修正的拉伸指数模型接近。尤其对于四川长宁M_S6.0余震序列和云南彝良M_S5.7余震序列,指数模型表现优于其他两种模型。指数模型参数具有明确的物理意义：参数A与r之和能够准确代表强震后的实际初始余震数,5个余震序列初始余震数偏差均小于1.70%;参数k可作为反映余震序列衰减快慢的特征值,k值越大则余震序列衰减越慢,其值与主震震级呈反比例关系。     

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.     

新源、和静交界ML6.8级地震序列应力降和震源机制研究

Based on the digital waveform data recorded by Xinjiang Digital Seismic Network for the Xinyuan-Hejing M_L6.8 earthquake sequences of June 30,2012,this paper analyzes the stress drops of earthquake sequences and the correlation coefficients of focal mechanisms significant for strong aftershocks.Firstly,the source parameters of the Xinyuan-Hejing M_L6.8 earthquake sequences are obtained by applying the spectrum analysis and the Brunes source model.Then,the correlation coefficients of spectral amplitudes are calculated using the low-frequency spectral amplitude recorded by the same station for the different events.Finally,based on the results of the correlation coefficients of spectral amplitudes,the events with similar focal mechanisms are grouped using the clustering method.The results show that:(1)The stress drop values show a steady trend in the aftershock sequence calm period and the stress drop values show a rise-fall in strong aftershocks.(2)The moving average correlation coefficient of amplitude spectrum begins to spread after the main shock.It shows that the correlation decreases between the main shock and the aftershocks in mechanisms.(3)The results of focal mechanism groups show that the earthquake sequences are mainly strike slips.The stress distribution of the main pressure axis is nearly NS,which is the same as the structural stress field.(4)The magnitude and mechanism show that there is an agreement before the strong aftershock,which shows that the regional stress field is enhanced.