Universal precursor seismicity pattern before locked-segment rupture and evolutionary rule for landmark earthquakes

Despite extensive investigations, no precursor patterns for reliably predicting major earthquakes have thus far been identified. Seismogenic locked segments that can accumulate adequate strain energy to cause major earthquakes are highly heterogeneous and low brittle. The progressive cracking of the locked segments with these properties can produce an interesting seismic phenomenon: a landmark earthquake and a sequence of smaller subsequent earthquakes (pre-shocks) always arise prior to another landmark earthquake within a well-defined seismic zone and its current seismic period. Applying a mechanical model, magnitude constraint conditions, and case study data of 62 worldwide seismic zones, we show that two adjacent landmark earthquakes reliably occur at the volume-expansion point and peak-stress point (rupture) of a locked segment; thus, the former is an identified precursor for the latter. Such a precursor seismicity pattern before the locked-segment rupture has definite physical meanings, and it is universal regardless of the focal depth. Because the evolution of landmark earthquakes follows a deterministic rule described by the model, they are predictable. The results of this study lay a firm physical foundation for reliably predicting the occurrence of future landmark earthquakes in a seismic zone and can greatly improve our understanding of earthquake generation mechanism.     

Global strain rates in western to central Himalayas and their implications in seismic hazard assessment

The Himalayas has experienced varying rates of earthquake occurrence in the past in its seismo-tectonically distinguished segments which may be attributed to different physical processes of accumulation of stress and its release, and due diligence is required for its inclusion for working out the seismic hazard. The present paper intends to revisit the various earthquake occurrence models applied to Himalayas and examines it in the light of recent damaging earthquakes in Himalayan belt. Due to discordant seismicity of Himalayas, three types of regions have been considered to estimate larger return period events. The regions selected are (1) the North-West Himalayan Fold and Thrust Belt which is seismically very active, (2) the Garhwal Himalaya which has never experienced large earthquake although sufficient stress exists and (3) the Nepal region which is very seismically active region due to unlocked rupture and frequently experienced large earthquake events. The seismicity parameters have been revisited using two earthquake recurrence models namely constant seismicity and constant moment release. For constant moment release model, the strain rates have been derived from global strain rate model and are converted into seismic moment of earthquake events considering the geometry of the finite source and the rates being consumed fully by the contemporary seismicity. Probability of earthquake occurrence with time has been estimated for each region using both models and compared assuming Poissonian distribution. The results show that seismicity for North-West region is observed to be relatively less when estimated using constant seismicity model which implies that either the occupied accumulated stress is not being unconfined in the form of earthquakes or the compiled earthquake catalogue is insufficient. Similar trend has been observed for seismic gap area but with lesser difference reported from both methods. However, for the Nepal region, the estimated seismicity by the two methods has been found to be relatively less when estimated using constant moment release model which implies that in the Nepal region, accumulated strain is releasing in the form of large earthquake occurrence event. The partial release in second event of May 2015 of similar size shows that the physical process is trying to release the energy with large earthquake event. If it would have been in other regions like that of seismic gap region, the fault may not have released the energy and may be inviting even bigger event in future. It is, therefore, necessary to look into the seismicity from strain rates also for its due interpretation in terms of predicting the seismic hazard in various segments of Himalayas.     

Did the 7/9/1999 M5.9 Athens Earthquake Come with a Warning?

Prior to the 7/9/1999 M_S = 5.9 Athens earthquake, regional seismicity has exhibited a power-law increase, of the form = K+A(t_c - t)ⁿ, where is estimated using an expression log = cM + d and t_c is the time of the culminating event. Such changes appeared after the 17/8/1999 M7.4 Izmit event. We quantified the performance of the power law vs. the null hypothesis of constant seismic release rates, by defining the curvature C as the ratio of the power law fit RMS/linear fit RMS, so that the smaller C is, the better the power law behaviour. By mapping C, we have established a critical radius of 110 km and observed that the region of correlated accelerating seismic release extended from the N. Aegean, through Euboea and Attica to the SW Peloponnese. A few days prior to the Athens event, min(C) was centred at the epicentral area and numerical simulation yielded t_c = 1999.676 and predicted M_S = 5.77. Seismicity rates returned to normal (quasi-constant) after the Athens event. We interpret this effect as critical point behaviour, following remote excitation of a broad area by stress redistribution due to the Izmit event which, at Athens, has triggered `premature' failure of a fault nearing its load bearing capacity. If this is correct, we have documented a case of remote earthquake triggering by another earthquake, as well as insight into the mechanisms producing it. As a corollary, we note that a large event may beget another large event in its broader region of interaction, which may be preceded by characteristic precursory seismicity changes.     

Use of tsunami waveforms for earthquake source study

Tsunami waveforms recorded on tide gauges, like seismic waves recorded on seismograms, can be used to study earthquake source processes. The tsunami propagation can be accurately evaluated, since bathymetry is much better known than seismic velocity structure in the Earth. Using waveform inversion techniques, we can estimate the spatial distribution of coseismic slip on the fault plane from tsunami waveforms. This method has been applied to several earthquakes around Japan. Two recent earthquakes, the 1968 Tokachi-oki and 1983 Japan Sea earthquakes, are examined for calibration purposes. Both events show nonuniform slip distributions very similar to those obtained from seismic wave analyses. The use of tsunami waveforms is more useful for the study of unusual or old earthquakes. The 1984 Torishima earthquake caused unusually large tsunamis for its earthquake size. Waveform modeling of this event shows that part of the abnormal size of this tsunami is due to the propagation effect along the shallow ridge system. For old earthquakes, many tide gauge records exist with quality comparable to modern records, while there are only a few good quality seismic records. The 1944 Tonankai and 1946 Nankaido earthquakes are examined as examples of old events, and slip distributions are obtained. Such estimates are possible only using tsunami records. Since tide-gauge records are available as far back as the 1850s, use of them will provide unique and important information on long-term global seismicity.     

地震迁移的类型、特征及机制讨论

地震活动的迁移是指地震沿着某一方向循序地发生,是地震活动总体无序中局部出现的有序结构。通过系统分析区域上典型的地震迁移现象可以发现,迁移可划分为沿断裂走向的纵向迁移、垂直断裂走向的横向迁移与岩石圈尺度的深源迁移三大类。结合具体的实例分析,可初步归纳出纵向迁移（包括单向、双向、反复和跳跃式迁移等常见形式）、横向迁移和深源地震迁移的主要特征,并初步估算出了不同类型迁移的速度值,其中沿全球板块边界纵向迁移平均速度约为V=569 km/a,沿亚板块边界的平均速度约为V=120 km/a,沿大陆内部断裂带平均速度约为V=50 km/a。横向迁移相对比较复杂,初步发现在东亚存在两种速度分别为约20 km/a、50 km/a的迁移现象。而深源地震迁移速度的全球平均值大约为360 km/a。地震的纵向与横向迁移都存在不同层次和级别,也存在多种不同频率、能量与速度的形变波与迁移现象,这很可能是区域上大地震丛集活动过程中断层相互作用、地震应力触发和岩石圈尺度的形变波传播等因素共同作用的结果,而这种大空间尺度上的地震迁移现象的存在及其所具有的规律性特征,显然可为开展区域地震危险性分析提供新的思路和方法参考。      

Time-dependent analysis of aftershock events and structural impacts on intraplate crustal seismicity of the Van earthquake (Mw 7.1, 23 October 2011), E-Anatolia

The Van earthquake (M _W 7.1, 23 October 2011) in E-Anatolia is typical representative of intraplate earthquakes. Its thrust focal character and aftershock seismicity pattern indicate the most prominent type of compound earthquakes due to its multifractal dynamic complexity and uneven compressional nature, ever seen all over Turkey. Seismicity pattern of aftershocks appears to be invariably complex in its overall characteristics of aligned clustering events. The population and distribution of the aftershock events clearly exhibit spatial variability, clustering-declustering and intermittency, consistent with multifractal scaling. The sequential growth of events during time scale shows multifractal behavior of seismicity in the focal zone. The results indicate that the extensive heterogeneity and time-dependent strength are considered to generate distinct aftershock events. These factors have structural impacts on intraplate seismicity, suggesting multifractal and unstable nature of the Van event. Multifractal seismicity is controlled by complex evolution of crustal-scale faulting, mechanical heterogeneity and seismic deformation anisotropy. Overall seismicity pattern of aftershocks provides the mechanism for strain softening process to explain the principal thrusting event in the Van earthquake. Strain localization with fault weakening controls the seismic characterization of Van earthquake and contributes to explain the anomalous occurrence of aftershocks and intraplate nature of the Van earthquake.     

Tectonic stress,seismicity, and seismic hazard in the southeastern Carpathians

Intermediate-depth earthquakes in the Vrancea region occur in response to stress generation due to descending lithosphere beneath the southeastern Carpathians. In this article, tectonic stress and seismicity are analyzed in the region on the basis of a vast body of observations. We show a correlation between the location of intermediate-depth earthquakes and the predicted localization of maximum shear stress in the lithosphere. A probabilistic seismic hazard assessment (PSHA) for the region is presented in terms of various ground motion parameters on the utilization of Fourier amplitude spectra used in engineering practice and risk assessment (peak ground acceleration, response spectra amplitude, and seismic intensity). We review the PSHA carried out in the region, and present new PSHA results for the eastern and southern parts of Romania. Our seismic hazard assessment is based on the information about the features of earthquake ground motion excitation, seismic wave propagation (attenuation), and site effect in the region. Spectral models and characteristics of site-response on earthquake ground motions are obtained from the regional ground motion data including several hundred records of small and large earthquakes. Results of the probabilistic seismic hazard assessment are consistent with the features of observed earthquake effects in the southeastern Carpathians and show that geological factors play an important part in the distribution of the earthquake ground motion parameters.     

Earthquake activity related to seismic cycles in a model for a heterogeneous strike-slip fault

We investigate the evolution of seismicity within large earthquake cycles in a model of a discrete strike-slip fault in elastic solid. The model dynamics is governed by realistic boundary conditions consisting of constant velocity motion of regions around the fault, static/kinetic friction and dislocation creep along the fault, and 3D elastic stress transfer. The fault consists of brittle parts which fail during earthquakes and undergo small creep deformation between events, and aseismic creep cells which are characterized by high ongoing creep motion. This mixture of brittle and creep cells is found to generate realistic aftershock sequences which follow the modified Omori law and scale with the mainshock size. Furthermore, we find that the distribution of interevent times of the simulated earthquakes is in good agreement with observations. The temporal occurrence, however, is magnitude-dependent; in particular, the small events are clustered in time, whereas the largest earthquakes occur quasiperiodically. Averaging the seismicity before several large earthquakes, we observe an increase of activity and a broadening scaling range of magnitudes when the time of the next mainshock is approached. These results are characteristics of a critical point behavior. The presence of critical point dynamics is further supported by the evolution of the stress field in the model, which is compatible with the observation of accelerating moment release in natural fault systems.     

强震的孕育规律与孕震模式

针对强震能否预测以及如何预测的科学难题,建立孕震断层多锁固段脆性破裂理论,发现强震孕育过程的指数规律：sf（k）=1.48ksc,其中sf（k）和sc分别为第k个锁固体断裂点与第一个锁固段膨胀起点对应的累加Be-nioff应变,可以利用锁固段在其变形膨胀点处开始发生的震群事件（加速性地震活动前兆）预测未来大震,并给出了强震四要素相关预测方法。通过对诸多历史强震（如邢台地震、海城地震、汶川地震、玉树地震等）的回溯性检验分析表明：强震可以预测,且其孕震过程都遵循着上述简单的共性力学规律。在此基础上,归纳出4种典型强震的孕震模式,即大震震级呈＂大—小—大＂型,大震震级呈连续上升型,锁固段快速连续破裂型与标准型。此外,根据相关强震预测理论方法,对有关抗震救灾未来研究的方向提出如下建议：建议加强活动断裂位置精确定位、性质判定的地震地质研究,并开展孕震区锁固段（闭锁区域）判识的地质与地球物理研究等。     

Seismic reversal pattern for the 1999 Chi-Chi, Taiwan, MW 7.6 earthquake

We in this study have calculated the standard normal deviate Z-value to investigate the variations in seismicity patterns in the Taiwan region before and after the Chi-Chi earthquake. We have found that the areas with relatively high seismicity in the eastern Taiwan became abnormally quiet before the Chi-Chi earthquake while the area in the central Taiwan with relatively low seismicity showed unusually active. Such a spatially changing pattern in seismicity strikingly demonstrates the phenomenon of “seismic reversal,” and we here thus present a complete, representative cycle of “seismic reversal” embedding in the changes of seismicity patterns before and after the Chi-Chi earthquake.     

Temporal behavior of the background seismicity rate in central Japan, 1998 to mid-2003

We studied the temporal behavior of the background shallow seismicity rate in 700 circular areas across inland Japan. To search for and test the significance of the possible rate changes in background seismicity, we developed an efficient computational method that applies the space–time ETAS model proposed by Ogata in 1998 to the areas. Also, we conducted Monte Carlo tests using a simulated catalog to validate the model we applied. Our first finding was that the activation anomalies were found so frequently that the constant background seismicity hypothesis may not be appropriate and/or the triggered event model with constraints on the parameters may not adequately describe the observed seismicity. However, quiescence occasionally occurs merely by chance. Another outcome of our study was that we could automatically find several anomalous background seismicity rate changes associated with the occurrence of large earthquakes. Very significant seismic activation was found before the M6.1 Mt. Iwate earthquake of 1998. Also, possible seismic quiescence was found in an area 150 km southwest of the focal region of the M7.3 Western Tottori earthquake of 2000. The seismicity rate in the area recovered after the mainshock.     

Seismically induced landslide displacements: a predictive model

Newmark's model for predicting earthquake-induced landslide displacements provides a simple way to predict the coseismic displacements affecting a sliding mass subject to earthquake loading. In this model, seismic slope stability is measured in terms of critical acceleration, which depends on the mechanical soil properties, pore-pressure distribution, and slope geometry. The triggering seismic forces are investigated in terms of energy radiation from the source, propagation, and site effects, based on 190 accelerometric recordings from 17 Italian earthquakes with magnitudes between 4.5 and 6.8. The method is based on the calibration of relations having the general form of an attenuation law that relates the energy of the seismic forces to the dynamic shear resistances of the sliding mass to propagate the expected landslide displacements as an inverse function of the distance from the fault rupture; the amount of displacement computed through these relations provides a criterion to predict the occurrence of slope failures. Finally, maps showing, in a deterministic and a probabilistic way, the potential of seismically induced landslide displacements are displayed as a tool to provide seismic landslide scenarios and earthquake-induced landslide hazard maps, respectively.     

Seismic precursors to a 2017 Nuugaatsiaq,Greenland, earthquake–landslide–tsunami event

High-frequency (5–20 Hz) seismic signals precursory to and embedded within the June 17, 2017 M_L?=?4 earthquake–landslide event are analyzed. This event in western Greenland generated a tsunami in Karrat fjord inundating Nuugaatsiaq village 32 km distant. Spectrogram and wavelet analyses of seismic data from the Greenland Ice Sheet Monitoring Network (GLISN) corroborate observations of seismic precursors at Nuugaatsiaq reported by Poli (Geophys Res Lett 44:8832–8836, 2017) and Caplan-Auerbach (in: AGU fall meeting abstracts, 2017) and reveal additional high-frequency arrivals being generated after the apparent initiation of fault rupture. New observations of seismic precursors 181 km from the Event at Upernavik, Greenland are correlated with those seen at Nuugaatsiaq. Wavelet analysis presents?>?100 significant energy peaks accelerating up to and into the earthquake–landslide event. The precursor events show a distinct, power law distribution, characterized by b values of ~?2.4. Results are compared and contrasted with small precursors observed in the studies of a natural chalk cliff landslide at Mesnil-Val, Haute Normandie, France. The earthquake–landslide appears to have been initiated by seismic precursors located at the fault scarp, leading to a small seismic foreshock and small landslide initiation, followed by a larger earthquake at the fault scarp, precipitating the primary landslide into the Karrat Fjord, which caused the subsequent tsunami.

     

A database for worldwide seismicity quantification

The first step in a seismicity analysis usually consists of defining the seismogenic units, seismic zones or individual faults. The worldwide delimitation of these zones involves an enormous effort and is often rather subjective. Also, a complete recording of faults will not be available for a long time yet. The seismicity model presented in this paper therefore is not based on individually defined seismic zones but rather on the assumption that each point in a global 1/2° grid of coordinates represents a potential earthquake source. The corresponding seismogenic parameters are allocated to each of these points. The earthquake occurrence frequency, one of the most important parameters, is determined purely statistically by appropriately spreading out the positions of past occurrences. All the other significant seismicity characteristics, such as magnitude-frequency relations, maximum possible magnitude and attenuation laws including the dependence on focal depth are determined in a global 1/2° grid of co-ordinates. This method of interpreting seismicity data allows us to establish a transparent, sufficiently precise representation of seismic hazard which is ideally suited for computer-aided risk analyses.     

震源区能量积累和释放过程的熵模型基本特征

在地震活动区的局部地壳地震活动性很大程度上是随机的,但在某些情况下,小的地方震的震级时间序列却具有确定性的分量,此分量很可能与一个大地震的成核有关。当小地震事件中最大的事件变小,最小的事件变大,并且它们的差别不断地减小,这个分量在地震记录上就表现为由震级的两个反向实时趋势产生的所谓能量楔。在一个大的成核事件的震源区,利用相图法,笔者依据非线性动力学已经解释了地震过程的演化和小震的大小分布。模拟地震过程的这种新的处理方法和数学模型已经被应用于来自世界各地区的大批地震目录数据,特别是中国的地震数据。     

The environmental effects of the 1743 Salento earthquake (Apulia,southern Italy): a contribution to seismic hazard assessment of the Salento Peninsula

The aim of this study was to provide a contribution to seismic hazard assessment of the Salento Peninsula (Apulia, southern Italy). It is well known that this area was struck by the February 20, 1743, earthquake (I ₀ = IX and M _w = 7.1), the strongest seismic event of Salento, that caused the most severe damage in the towns of Nardò (Lecce) and Francavilla Fontana (Brindisi), in the Ionian Islands (Greece) and in the western coast of Albania. It was also widely felt in the western coast of Greece, in Malta Islands, in southern Italy and in some localities of central and northern Italy. Moreover, the area of the Salento Peninsula has also been hit by several low-energy and a few high-energy earthquakes over the last centuries; the instrumental recent seismicity is mainly concentrated in the western sector of the peninsula and in the Otranto Channel. The Salento area has also experienced destructive seismicity of neighboring regions in Italy (the Gargano Promontory in northern Apulia, the Southern Apennines chain, the Calabrian Arc) and in the Balkan Peninsula (Greece and Albania). Accordingly, a critical analysis of several documentary and historical sources, as well as of the geologic–geomorphologic ground effects due to the strong 1743 Salento earthquake, has been carried out by the authors in this paper; the final purpose has been to re-evaluate the 1743 MCS macroseismic intensities and to provide a list of newly classified localities according to the ESI-07 scale on the base of recognized Earthquake Environmental Effects. The result is a quite different damage scenario due to this earthquake that could raise the seismic potential currently recognized for the Salento area, and consequently upgrade the seismic hazard classification of the Salento. Indeed it is important to remind that currently, despite the intense earthquake activity recorded not only in the Otranto Channel, but especially in Greece and Albania, this area is classified in the least dangerous category of the Seismic Classification of the Italian territory (IV category).     

Cautious neotectonic hypotheses for assessing the seismic hazard in northeastern Italy

A simplified tectonic scheme for hazard purposes was recently adopted for northeastern Italy, introducing large generalized seismogenic areas containing systems of complex geometry faults. This scheme considers only major faults with documented seismic activity. In the present analysis, a different tectonic scheme, with linear elements as seismogenic sources, is presented. The assessment of the regional seismic hazard is done with the fault rupture model, its most important advantage being the recognition that the length of fault rupture during an earthquake is an important consideration in probabilistic calculations of seismic hazard. Moreover, some structures with no associated seismicity but with notable neotectonic activity are considered, and their contribution to the results investigated. Important uncertainties such as those in the maximum possible magnitude of future earthquakes, in the location of the fault, in the focal depth, and in the attenuation law are accounted for in the calculations and their influence studied. The results identify a seismic belt running from Lake Garda to Friuli and along the Yugoslav coast and are very similar to those already known for Friuli, with the largest values corresponding to the zone around Gemona. Some slight differences in the shape of the areas of equal acceleration are probably due to the delineation of the seismic sources of the proposed model. For a cautious elaboration, some neotectonic lines without present seismicity were added into the fault model. Their contribution is negligible in the areas of highest acceleration, but increases remarkably in the areas where acceleration is not expected to exceed the medium values.     

储集层中高压流体引爆强地震的机理——以5.12汶川地震为例

目前产生地震的机制仍以弹性回跳说为主:地震是因为断层错断使岩层的弹性能释放而引发.但越来越多的学者开始质疑,仅断层错断后的弹性能,是否真能达到实际地震所释放的巨大能量.因此,有必要探讨地震初动后破坏性强震的性质及其真正的能量来源.文章根据沉积地层中的储集层及其压力的特点分析得出,储集层内含有大量的高压流体,其压力在一定...     

岩石I型裂纹定向扩展规律试验研究

为精准获得岩石I型裂纹扩展演化全过程,采用一种简易裂纹定向扩展装置开展了不同岩性试样裂纹扩展试验研究,借助声发射及数字散斑技术对裂纹扩展全过程进行了监测,建立了裂纹定向扩展力学模型,分析了裂纹扩展过程中声发射及变形场的演化规律,提出了评价岩石I型裂纹扩展难易程度的能量指标CE,探讨了I型裂纹定向起裂扩展机制。结果表明：该简易裂纹定向扩展装置能够有效实现I型裂纹沿预定方向稳定扩展,其起裂角均小于10º,同时通过简化力学模型计算得到白砂岩、灰砂岩的裂纹扩展峰值强度与巴西劈裂抗拉强度相比偏差分别为22.76%、7.53%;根据变形场演化规律,可将裂纹扩展分为微裂隙发育（散斑变形场分区不明显）、主控裂纹孕育（散斑变形场出现分区现象）和主控裂纹扩展3个阶段;声发射演化过程可分为平静期、缓增期、急增期和降低期4个阶段,由于灰砂岩相较于白砂岩质地更致密、更坚硬,导致其声发射平静期长,而后3个阶段持续时间短;将载荷−位移曲线峰前与峰后的面积之比定义为评价岩石I型裂纹扩展难易程度的能量指标CE,计算得到灰砂岩、白砂岩的CE分别为13～16、1～2,表明CE可有效评价岩石I型裂纹扩展难易程度;岩石I型裂纹起裂扩展机制可概况为：在加载峰值前裂隙尖端受最大拉应力作用,存储的弹性能快速增加、耗散能缓慢增加,但在加载峰值后裂隙尖端存储弹性能超过其储能极限迅速释放,此时输入能大部分转化为耗散能促进主控裂纹快速扩展。后续将对裂纹定向扩展试验装置进一步优化改进,以期为裂纹扩展机制、岩石破坏前兆信息、裂纹止裂原理等研究提供一种新方法,同时为工程现场煤岩层定向爆破、压裂、止裂等相关技术优化提供理论指导。