Temporal evolution of the focal mechanism consistency of the 2021 Yangbi MS 6.4 earthquake sequence in Yunnan

Using the Cut And Paste (CAP) method, we invert the focal mechanism of 38 moderate earthquakes (M_S ≥ 3.0) recorded by Yunnan seismic network and analyze the corresponding focal mechanism consistency based on the minimum spatial rotation angle. Our results indicate that the M_S 6.4 mainshock is induced by a lateral strike slip fault (with a rake angle of ～ ?165°) and a little normal-faulting component event along a nearly vertical plane (dipping angle～ 79° and strike ～138°). Combining our results with high resolution catalog, we argue that the seismogenic fault of this earthquake sequence is a secondary fault western to the major Weixi-Qiaohou-Weishan fault. The focal mechanism evolution can be divided into three periods. During the first period, the foreshock sequence, the focal mechanism consistency is the highest (KA<36°); during the second period which is shortly after the mainshock, the focal mechanism shows strong variation with KA ranging from 8° to 110°; during the third period, the seismicity becomes weak and the focal mechanism of the earthquakes becomes more consistent than the second period (18°<KA<73°). We suggest that the KA, to some extent, represents the coherence between local tectonic stress regime and the stress state of each individual earthquake. Furthermore, high focal mechanism consistency and high linearity of seismic distribution may serve as indicators for the identification of foreshock sequence.     

Historical seismicity near Chagos: a complex deformation zone in the equatorial Indian Ocean

Over the last twenty years, Chagos Bank has a seismicity rate disproportionate to its supposed intraplate location. Earthquake relocation also shows a high seismicity rate in pre-WWSSN time (1912–1963), with seven events located off of the Central Indian Ridge, including large events in 1912 (M = 6.8) and 1944 (M = 7.2). This study uses the moment variance technique, a systematic search for the mechanism which best fits P, PP, SH, Love and Rayleigh amplitudes, to determine the focal mechanisms of two pre-WWSSN earthquakes. A test with a recent event of known mechanism demonstrates that accurate focal parameter determination is possible even when only a few good records are available. Moment variance analysis shows a thrust faulting mechanism for the 1944 event, northeast of Chagos Bank near the Chagos-Laccadive ridge, and a strike-slip focal mechanism for a smaller 1957 event west of Chagos Bank. The 1944 event, one of the largest oceanic “intraplate” earthquakes known (moment 1.4 × 10²⁷ dyne-cm), indicates that the Chagos seismicity reflects not an isolated occurrence of normal faulting as previously thought, but rather regional tectonic deformation extending northeast of Chagos Bank and including thrust, normal and strike-slip events. This seismicity and previously studied seismicity near the Ninetyeast Ridge and Central Indian Basin suggest a broad zone of deformation stretching across the equatorial Indian Ocean. This zone contains all known magnitude seven oceanic “intraplate” earthquakes not associated with subduction zones or continental margins, suggesting that elsewhere such extensive deformation occurs only along plate boundaries. This study proposes that a slow, diffuse plate boundary extends east from the Central Indian Ridge to the Ninetyeast Ridge and north to the Sumatra Trench. A recent plate motion study confirms this boundary and suggests that it separates the Australian plate from a single Indo-Arabian plate.     

The Gilbert Islands (Republic of Kiribati) earthquake swarm of 1981–1983

A major swarm of intraplate earthquakes at the southeastern end of the Gilbert Islands Chain (3.5°S, 177.5°E) commenced in December 1981 and lasted through March 1983. No seismicity had been reported in the historical record in this region prior to 1981, but during the swarm 217 events with m_b ? 4.0 were located by the NEIS, with 86 events having m_b ? 5.0. The source region is quite remote, and the uniform detection level for the NEIS is for m_b ? 4.8. A b-value of 1.35 is found for the swarm using the maximum likelihood method. Four events in the sequence were large enough (m_b = 5.6?5.9) to determine focal mechanisms teleseismically using body- and surface-wave analysis. These events are found to have a range of mechanisms, from predominantly thrust with a significant oblique component, to purely strike-slip. The compression axes are consistent for all four events, with horizontal orientation trending NNE-SSW. This orientation is perpendicular to the direction of plate motion. The events are located at depths between 15 and 20 km placing them deep in the oceanic crust or in the upper mantle. No obvious bathymetric feature can be related to the fault plane orientations, though there is an offset in the island chain near the epicenters. While some characteristics of the swarm suggest a magmatic origin, the nature of the focal mechanisms, the location of the swarm, and the large accumulated moment release of the sequence favor a tectonic origin.     

Source Characteristics of the 22 January 2003 M w?=?7.5 Tecom��n, Mexico, Earthquake: New Insights

Aftershock locations, source parameters and slip distribution in the coupling zone between the overriding North American and subducted Rivera and Cocos plates were calculated for the 22 January 2003 Tecomán earthquake. Aftershock locations lie north of the El Gordo Graben with a northwest-southeast trend along the coast and superimposed on the rupture areas of the 1932 (M _w?=?8.2) and 1995 (M _w?=?8.0) earthquakes. The Tecomán earthquake ruptured the northwest sector of the Colima gap, however, half of the gap remains unbroken. The aftershock area has a rectangular shape of 42?±?2 by 56?±?2?km with a shallow dip of roughly 12° of the Wadati-Benioff zone. Fault geometry calculated with the Náb??lek (1984) inversion procedure is: (strike, dip, rake)?=?(277°, 27°, 78°). From the teleseimic body wave spectra and assuming a circular fault model, we estimated source duration of 20?±?2?s, a stress drop of 5.4?±?2.5?MPa and a seismic moment of 2.7?±?.7?×?10²⁰?Nm. The spatial slip distribution on the fault plane was estimated using new additional near field strong motion data (54?km from the epicenter). We confirm their main conclusions, however we found four zones of seismic moment release clearly separated. One of them, not well defined before, is located toward the coast down dip. This observation is the result of adding new data in the inversion. We calculated a maximum slip of 3.2?m, a source duration of 30?s and a seismic moment of 1.88?×?10²⁰?Nm.     

2013年山东莱州M4.6地震序列发震构造初探

采用双差定位法对山东莱州地震序列重新定位,通过CAP方法反演M4.6地震震源机制,在此基础上初步探讨莱州地震序列发震构造。结果显示:精确定位震中位置主要位于柞村—仙夼断裂的NW方向,深度剖面显示从SE方向到NW方向断层深度呈由浅逐渐变深的趋势,这均与柞村—仙夼断裂位置、走向、倾向特征较为吻合;M4.6地震震源机制解的节面Ⅰ与柞村—仙夼断裂走向、倾角较为接近。综合精确定位震中位置、剖面深度分布特征、M4.6地震震源机制解及宏观调查烈度分布等结果与柞村-仙夼断裂产状之间的关系,初步推测柞村—仙夼断裂可能为莱州地震序列的发震断层。     

The Kaliningrad, West Russia earthquakes on the 21st of September 2004—Surprise events in a very low-seismicity area

The earthquakes in Kaliningrad, West Russia on the 21st of September 2004 were unexpected in a seismically quiet area. The main shock of magnitude m_b = 5 was widely felt around the Baltic Sea. A comparison with some historic earthquakes in Northern Europe shows that its perceptibility area is smaller than that of the 1904 Oslo Graben earthquake of an estimated magnitude (M_L) 5.4 but larger than those of the 1759 Kattegat and 1819 Lurøy earthquakes. The latter are claimed to have had magnitudes (M_S) in the range of 5.7-6.0. An analysis of the Lg phase of the Kaliningrad earthquakes as recorded at a number of European stations accords only weakly with the macroseismic intensity pattern that shows fast attenuation towards west-northwest and southeast. The strike-slip focal mechanism of the main shock is discussed in the context of remnant glacial rebound stresses in generating present-day seismicity in N. Europe.     

2019年10月2日贵州沿河MS4.9地震震害特征及发震构造分析

贵州沿河M_S4.9地震发生在历史地震强度较低的上扬子地块凤冈SN向隔槽式褶皱变形区。通过地震地质背景分析、震害调查、震源机制解、断层调查和库区水位变化情况等,得到主要认识如下:由于震源深度浅、灾区老旧自建房抗震性能差,导致本次地震直接经济损失严重;本次地震主震的机制解为节面Ⅰ:走向61°/倾角35°/滑动角135°,节面Ⅱ:走向190°/倾角66°/滑动角63°,表现为走向NE、逆冲兼平移型运动方式;结合等震线走向及震中主要断层性质,判断NE向沿河断层为本次地震主震的发震构造,并进一步推测此次地震为水库诱发断层活化引起的地震。     

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

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

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

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

Historical seismogram reproductions for the source parameters determination of the 1902, Atushi (Kashgar) earthquake

The majority of original seismograms recorded at the very beginning of instrumental seismology (the early 1900s) did not survive till present. However, a number of books, bulletins, and catalogs were published including the seismogram reproductions of some, particularly interesting earthquakes. In case these reproductions contain the time and amplitude scales, they can be successfully analyzed the same way as the original records. Information about the Atushi (Kashgar) earthquake, which occurred on August 22, 1902, is very limited. We could not find any original seismograms for this earthquake, but 12 seismograms from 6 seismic stations were printed as example records in different books. These data in combination with macroseismic observations and different bulletins information published for this earthquake were used to determine the source parameters of the earthquake. The earthquake epicenter was relocated at 39.87° N and 76.42° E with the hypocenter depth of about 18 km. We could further determine magnitudes m _B = 7.7 ± 0.3, M _S = 7.8 ± 0.4, M _W = 7.7 ± 0.3 and the focal mechanism of the earthquake with strike/dip/rake ? 260°± 20/30°± 10/90°± 10. This study confirms that the earthquake likely had a smaller magnitude than previously reported (M8.3). The focal mechanism indicates dominant thrust faulting, which is in a good agreement with presumably responsible Tuotegongbaizi-Aerpaleike northward dipping thrust fault kinematic, described in previous studies.     

Activation of seismicity in Central and South Asia after the Makran earthquakes: Possible acceleration of preparation of large seismic events in the Tien Shan region

Two zones of seismicity (ten events with M _w = 7.0–7.7) stretching from Makran and the Eastern Himalaya to the Central and EasternTien Shan, respectively, formed over 11 years after the great Makran earthquake of 1945 (M _w = 8.1). Two large earthquakes (M _w = 7.7) hit theMakran area in 2013. In addition, two zones of seismicity (M ≥ 5.0) occurred 1–2 years after theMakran earthquake in September 24, 2013, stretching in the north-northeastern and north-northwestern directions. Two large Nepal earthquakes struck the southern extremity of the “eastern” zone (April 25, 2015, M _w = 7.8 and May 12, 2015, M _w = 7.3), and the Pamir earthquake (December 7, 2015, M _w = 7.2) occurred near Sarez Lake eastw of the “western” zone. The available data indicate an increase in subhorizontal stresses in the region under study, which should accelerate the possible preparation of a series of large earthquakes, primarily in the area of the Central Tien Shan, between 70° and 79° E, where no large earthquakes (M _w ≥ 7.0) have occurred since 1992.     

The Campania-Lucania (southern Italy) earthquake of 23 November 1980

We present the main seismological results of our study of the Campania-Lucania earthquake of 23 November 1980. A complete set of far field and local data has been analysed. From long-period body waves data we determine the fault plane solution (φ₁ = 140°,δ₁ = 60°,φ₂ = 75°,δ₂ = 54°), a depth of 15 km and calculate a seismic moment of 6 × 10²⁵ dyne cm and a source duration of 6 s. From data of a local network deployed immediately after the event we determine aftershock locations: they are aligned in a direction NW-SE that fit extremely well with the focal solution determined above. We can choose as fault plane the plane striking 140° and dipping at 60° and the event is a normal event with a large component of left-lateral strike slip. The source area evaluated from this aftershock distribution 14 km × 40 km is quite suitable for an earthquake of a seismic moment of 6 × 10²⁵ dyne cm.     

Implications of 2007’s Earthquake Activity in Eğirdir Lake (SW Anatolia) Based on Moment Tensor Solutions and Inversion of Stress State

We analyzed the waveforms of the small- to moderate-sized earthquakes that took place in the northern part of the inner Isparta Angle (IA) to retrieve their source parameters and combine these results with the focal mechanism solutions of the larger events that occurred in 2007 in E?irdir Lake at the apex of IA. In total, source mechanisms of 20 earthquakes within the magnitude range 3.5 < M < 5.0 were calculated using a regional moment tensor inversion technique. The inversion of the focal mechanisms yields an extensional regime with a NNE–SSW (N38°E) trending σ ₃ axis. Inversion results are related to a mainly WNW–ESE oriented normal fault beneath E?irdir Lake. The R value of a NNE–SSW extensional regime is 0.562 showing a triaxial stress state in the region. The current stress regime results from complex subduction processes such as slab pull, slab break-off, roll-back and/or retreating mechanism along the Hellenic and Cyprus arcs and the southwestward extrusion of the Anatolian block since the early Pliocene.     

Forecasting the Locations of Future Large Earthquakes: An Analysis and Verification

The objective of this paper is to quantify the use of past seismicity to forecast the locations of future large earthquakes and introduce optimization methods for the model parameters. To achieve this the binary forecast approach is used where the surface of the Earth is divided into l° × l° cells. The cumulative Benioff strain of m ≥ m _c earthquakes that occurred during the training period, ΔT _tr, is used to retrospectively forecast the locations of large target earthquakes with magnitudes ≥m _T during the forecast period, ΔT _for. The success of a forecast is measured in terms of hit rates (fraction of earthquakes forecast) and false alarm rates (fraction of alarms that do not forecast earthquakes). This binary forecast approach is quantified using a receiver operating characteristic diagram and an error diagram. An optimal forecast can be obtained by taking the maximum value of Pierce’s skill score.     

Universality of the Seismic Moment-frequency Relation

—We analyze the seismic moment-frequency relation in various depth ranges and for different seismic regions, using Flinn-Engdahl's regionalization of global seismicity. Three earthquake lists of centroid-moment tensor data have been used the Harvard catalog, the USGS catalog, and the Huang et al. (1997) catalog of deep earthquakes. The results confirm the universality of the β-values and the maximum moment for shallow earthquakes in continental regions, as well as at and near continental boundaries. Moreover, we show that although fluctuations in earthquake size distribution increase with depth, the β-values for earthquakes in the depth range of 0–500 km exhibit no statistically significant regional variations. The regional variations are significant only for deep events near the 660 km boundary. For declustered shallow earthquake catalogs and deeper events, we show that the worldwide β-values have the same value of 0.60 ± 0.02. This finding suggests that the β-value is a universal constant. We investigate the statistical correlations between the numbers of seismic events in different depth ranges and the correlation of the tectonic deformation rate and seismic activity (the number of earthquakes above a certain threshold level per year). The high level of these correlations suggests that seismic activity indicates tectonic deformation rate in subduction zones. Combined with the universality of the β-value, this finding implies little if any variation in maximum earthquake seismic moment among various subduction zones. If we assume that earthquakes of maximum size are similar in different depth ranges and the seismic efficiency coefficient, χ, is close to 100% for shallow seismicity, then we can estimate χ for deeper earthquakes for intermediate earthquakes χ≈ 5%, and χ≈ 1% for deep events. These results may lead to new theoretical understanding of the earthquake process and better estimates of seismic hazard.     

近震全波形反演2017年九寨沟M7.0地震序列震源机制解

搜集了四川地震台网的波形资料,采用全波形反演2017年8月8日九寨沟M7.0地震序列震源机制解.反演结果显示,九寨沟主震矩震级为M_W6.36,震源深度为22 km,节面I走向为150°,倾角为80°,滑动角为-20°;节面Ⅱ走向为244°,倾角为70°,滑动角为-169°.余震主要分布在14~22 km深度范围内,震源机制以走滑型为主,其中正断型地震2个,逆冲型地震2个,走滑型地震24个,混合型地震8个.断层面优势方向为SSE向,与塔藏断裂和虎牙断裂走向基本一致,但与塔藏断裂最南段存在明显差异.倾角变化集中在60°~80°,滑动角主要分布在0°附近,表明九寨沟地震序列主要受SSE走向、近似直立的左旋走滑断层控制.P轴优势方位为SEE向,仰角主要分布在30°以内,与区域应力场基本一致.震源区的机制类型和应力状态均存在空间分段差异.本文推测此次九寨沟M7.0地震序列可能发生在虎牙断裂向北延伸的隐伏断裂上,但不排除地震引起了塔藏断裂南段和虎牙断裂以北隐伏断裂同时破裂的可能.     

Asperity Distribution of the 1952 Great Kamchatka Earthquake and its Relation to Future Earthquake Potential in Kamchatka

—The 1952 Kamchatka earthquake is among the largest earthquakes of this century, with an estimated magnitude of M _w = 9.0. We inverted tide gauge records from Japan, North America, the Aleutians, and Hawaii for the asperity distribution. The results show two areas of high slip. The average slip is over 3 m, giving a seismic moment estimate of 155×10²⁰Nm, or M _w = 8.8. The 20th century seismicity of the 1952 rupture zone shows a strong correlation to the asperity distribution, which suggests that the large earthquakes (M > 7) are controlled by the locations of the asperities and that future large earthquakes will also recur in the asperity regions.     

褶皱构造中的地震——2017年三峡库区巴东M4.3地震序列成因讨论

2017年巴东M4.3地震序列为一发生在褶皱构造翼部的地震活动.地震精定位结果显示,整个地震序列呈NE向展布,震源深度较浅,M4.3和M4.1地震的震源深度分别为3.8 km和2.4 km.地震序列表现出NW浅、SE深的特征.震源机制解反演结果表明,两次较大地震的发震节面走向NEE,倾向SE,运动性质为走滑.震中区断裂构造分析发现,没有与发震节面吻合的断裂构造.结合地震定位结果和地质剖面分析,M4.3地震发生在三叠系岩溶及裂隙发育的厚层灰岩中,而M4.1地震及大多数小震活动则主要发生在巴东组红层（软弱滑脱层）中.地震发生时,三峡水库正处于低水位腾库容期间,库水卸荷使得库区原来由于荷载压实作用产生的压应力出现了局部回弹,进而引起部分裂隙的扩张,造成部分裂隙发育的岩体发生失稳滑动.此外,裂隙的扩张也为流体的渗透扩散提供了有利通道.流体渗透扩散一方面使得孔隙压力增大,有效应力降低;另一方面则使得软弱地层发生软化泥化,使得岩体发生顺层或切层的失稳滑动,造成了地震的发生.     

Body-wave waveform constraints on the source parameters of the Yangjiang,China, earthquake of July 25, 1969: A devastating earthquake in a stable continental region

On July 25, 1969, anM _s 5.9 earthquake shook the Kwangtung Province near Yangjiang, China. Casualties and extensive damage were reported in the epicentral area. The Yangjiang earthquake occurred within the passive margin located along China's southeast coast. This stable continental setting is seismically one of the most quiet regions of China; historic records indicate this earthquake to have been the first devastating one in the area. A remarkable feature of this earthquake sequence is that its foreshocks and aftershocks are relatively small in terms of number and magnitude despite the relatively large main shock.Waveforms of long-and short-period teleseismicP and long-period teleseismicSH waves have been modeled to estimate the source parameters for this earthquake. The focal mechanism determined is predominately strike-slip with a small normal component (strike=263°, dip=79°, rake=194°) and is in good agreement with observations in the epicentral area regarding dislocation direction of fissures, distribution of aftershocks, and isoseismals. The average seismic moment is 5.15×10²⁴ dyne-cm and the focal depth is estimated to be 9 km. The earthquake is characterized by a rather short source-time function and a high stress drop of approximately 380 bars. The observed pattern of the foreshock-main shock-aftershock sequence is interpreted in terms of the large stress drop associated with the main shock and the material properties at and near the source. TheP-axis orientation lies roughly NW-SE and is consistent with the maximum compressive stress observed along the coast of southeast China. This stress orientation is consistent with the NW subducting Philippine Sea plate and the SE ward push from central China due to the convergence between Indian and Eurasian subcontinents.     

Paleomagnetic results from the Skaergaard intrusion,East Greenland

The natural remanent magnetization of 22 out of a total of 31 oriented cores from the layered series of the Skaergaard gabbroic intrusion (age: 55 m.y.) in East Greenland shows good stability in thermal and AF testing. The average direction of 22 AF and 9 thermally treated specimens isD = 170°,I = ?59°,α₉₅ = 4.2 before correction for tilt. The mean directions after rotation around strike to horizontal and after rotation to original attitudes suggested by others yields poorer population statistics. It is therefore concluded that flexuring took place between solidification and acquisition of remanent magnetization, a range in temperature of about 500°C which may represent an interval of somewhat less than 250,000 years. No evidence for secular variation is observed which may also suggest slow cooling through the blocking temperature range. The polarity is reversed and the pole position without “tilt correction” is 165°E, 61°N,d_m = 6.2,d_p = 4.6, which is similar to pole positions reported by others for the overlying slightly older basalt.