东南亚地震活动与b 值的时空分布特征

东南亚地区是“21世纪海上丝绸之路”（以下简称“海洋丝路”）的重要组成部分,该区历史上曾发生十余次巨大地震,地震及其次生地质灾害是威胁东南亚地区经济社会发展和国际合作的主要自然灾害。系统梳理该区地震活动的时空分布特征及评估未来灾害风险格局,对于推进“一带一路”倡议实施及区域经济社会可持续发展具有重要意义。文章基于东南亚地区1900年以来M≥5地震的时空分布统计分析和地震b值计算,揭示出该区的地震活动在时间上表现出活跃期与平静期交替变化的特征;空间上表现出明显的聚集效应,成丛性强且主要集中在5个地震统计区内,其中印尼—马来多岛弧盆系地震区和菲律宾群岛地震区的地震活动最为活跃。总体而言,东南亚5个地震区的b值偏低,在0.42~0.91之间。该区内的地震b值也存在时空差异,受大地震事件、俯冲带年龄、活动断裂带和震源深度等众多因素影响,但主控因素在不同区域有所不同。地震b值时空变化特征对区域地震活动预测具有启示作用。上述认识为推进“海洋丝路”工程建设和“一带一路”防灾减灾对策提供了科学支撑。     

Crustal Thickness and Composition Variation along the Western Ghats of India through Teleseismic Receiver Function Analysis

Seismic-potential for the Andaman Arc System (AAS) is evaluated using a three-tier approach: (i) the seismic b-values derived using a revised and homogenized earthquake catalog for Mw ≥ 4.7, with uniform Mw, for the period 1964–2013 created from ISC Data, (ii) Free-air gravity-anomalies for AAS and their geologic interpretation and (iii) deep crustal structure from seismic reflection surveys and 3D seismic tomography results. Both long-term spatial b-value anomalies for the period 1964-2013 and temporal b-value anomalies for a shorter period (2002–2013) have been estimated. The b-value maps are interpreted in terms of the stress regime acting across AAS and stressed asperities inferred along the strike of AAS. Eight distinguishable seismic-zones are identified and their seismic potential is examined by temporal bvalue anomalies in producing moderate to large earthquakes. The latter demonstrates ‘low-high’ or ‘high-low’ couplet over years, and that a variation in b-value more than 20% compared to the previous year value is likely to produce an earthquake event with Mw ≥ 6.0. Some support to this interpretation comes from the regional Free-air gravity-anomalies and deep crustal structure interpreted from hypocentral distribution of earthquakes. The high b-values are seemingly associated with magma chambers or low velocity crust; creating asperity zones due to multiple batholithic intrusions at plate boundary. This geologic interpretation is evidenced by 3D P-wave seismic tomography and velocity heterogeneity study for AAS reported by us elsewhere.     

Probabilities of earthquake occurrences in Mainland Southeast Asia

The frequency–magnitude distributions of earthquakes are used in this study to estimate the earthquake hazard parameters for individual earthquake source zones within the Mainland Southeast Asia. For this purpose, 13 earthquake source zones are newly defined based on the most recent geological, tectonic, and seismicity data. A homogeneous and complete seismicity database covering the period from 1964 to 2010 is prepared for this region and then used for the estimation of the constants, a and b, of the frequency–magnitude distributions. These constants are then applied to evaluate the most probable largest magnitude, the mean return period, and the probability of earthquake of different magnitudes in different time spans. The results clearly show that zones A, B, and E have the high probability for the earthquake occurrence comparing with the other seismic zones. All seismic source zones have 100 % probability that the earthquake with magnitude ≤6.0 generates in the next 25 years. For the Sagaing Fault Zone (zones C), the next Mw 7.2–7.5 earthquake may generate in this zone within the next two decades and should be aware of the prospective Mw 8.0 earthquake. Meanwhile, in Sumatra-Andaman Interplate (zone A), an earthquake with a magnitude of Mw 9.0 can possibly occur in every 50 years. Since an earthquake of magnitude Mw 9.0 was recorded in this region in 2004, there is a possibility of another Mw 9.0 earthquake within the next 50 years.     

Seismicity, seismotectonics and seismic hazard of Italy

A first generation of probabilistic seismic hazard maps of the Italian country are presented. They are based on seismogenic zoning deriving from a kinematic model of the structural tectonic units and on an earthquake catalogue with the foreshock and aftershock events filtered out. The following ground motion parameters have been investigated and mapped using attenuation equations based on strong-motion recordings of Italian earthquakes: peak ground acceleration and velocity; Arias intensity; strong motion duration; and the pseudovelocity and pseudoacceleration spectral values at 14 fixed frequencies both for the vertical and the largest horizontal component. A Poissonian model of earthquake occurrence is assumed as a default and the hazard maps are presented in terms of ground motion values expected to be exceeded at a 10% probability level in 50 years (return period 475 years) according to the requirement of Eurocode 8 for the seismic classification of national territories, as well as in terms of exceedance probabilities of selected ground motion values. Finally, as a tentative study, the use of hybrid methods (implementing both seismogenic zones and structures), renewal processes (including earthquake forecasting) and the influence of site effects (as the basis for the planning of earthquake scenarios) were explored.     

Seismic coupling and uncoupling at subduction zones

Seismic coupling has been used as a qualitative measure of the “interaction” between the two plates at subduction zones. Kanamori (1971) introduced seismic coupling after noting that the characteristic size of earthquakes varies systematically for the northern Pacific subduction zones. A quantitative global comparison of many subduction zones reveals a strong correlation of earthquake size with two other variables: age of the subducting lithosphere and convergence rate. The largest earthquakes occur in zones with young lithosphere and fast convergence rates, while zones with old lithosphere and slow rates are relatively aseismic for large earthquakes. Results from a study of the rupture process of three great earthquakes indicate that maximum earthquake size is directly related to the asperity distribution on the fault plane (asperities are strong regions that resist the motion between the two plates). The zones with the largest earthquakes have very large asperities, while the zones with smaller earthquakes have small scattered asperities. This observation can be translated into a simple model of seismic coupling, where the horizontal compressive stress between the two plates is proportional to the ratio of the summed asperity area to the total area of the contact surface. While the variation in asperity size is used to establish a connection between earthquake size and tectonic stress, it also implies that plate age and rate affect the asperity distribution. Plate age and rate can control asperity distribution directly by use of the horizontal compressive stress associated with the “preferred trajectory” (i.e. the vertical and horizontal velocities of subducting slabs are determined by the plate age and convergence velocity). Indirect influences are many, including oceanic plate topography and the amount of subducted sediments.All subduction zones are apparently uncoupled below a depth of about 40 km, and we propose that the basalt to eclogite phase change in the down-going oceanic crust may be largely responsible. This phase change should start at a depth of 30–35 km, and could at least partially uncouple the plates by superplastic deformation throughout the oceanic crust during the phase change.     

成矿过程奇异性与矿床多重分形分布

本文将成矿过程与地震、火山、滑坡、洪水、暴雨、森林火灾等一系列非线性地球系统过程称之为奇异性过程（Singular Processes）。奇异性过程的基本特征是在相对很小的时间间隔或者空间范围中产生巨大能量释放或者物质超常富集或堆积。探讨了奇异性过程的“资源”与“灾害”双重效应;从地幔对流、板块构造等自组织临界性系统探讨了成矿过程的奇异性以及成矿结果服从多重分形分布规律,如矿床规模与个数、矿床品位与个数等;介绍了可以刻画成矿域空间结构的不均匀性、尺度独立性和白组织临界性的多重分形模型,如度量元素富集强度的局部奇异性模型,度量成矿域结构各向异性的广义自相似性模型,以及度量分形结构序次的分形谱系模型。这些非线性模型为矿产资源定量预测提供了新的思路和工具。     

青藏高原东北缘地震断层的研究

青藏高原东北缘新构造运动强烈,地震活动频繁。调查资料表明,伴随每次大地震的发生,都将产生规模不等的地震断层。本文研究了该区地震断层的分布规律、类型以及与先存断裂的关系等,同时对研究区不同地段的地壳稳定性进行了初步评价。      

Probabilistic seismic hazard analysis for Sumatra, Indonesia and across the Southern Malaysian Peninsula

The ground motion hazard for Sumatra and the Malaysian peninsula is calculated in a probabilistic framework, using procedures developed for the US National Seismic Hazard Maps. We constructed regional earthquake source models and used standard published and modified attenuation equations to calculate peak ground acceleration at 2% and 10% probability of exceedance in 50 years for rock site conditions. We developed or modified earthquake catalogs and declustered these catalogs to include only independent earthquakes. The resulting catalogs were used to define four source zones that characterize earthquakes in four tectonic environments: subduction zone interface earthquakes, subduction zone deep intraslab earthquakes, strike-slip transform earthquakes, and intraplate earthquakes. The recurrence rates and sizes of historical earthquakes on known faults and across zones were also determined from this modified catalog. In addition to the source zones, our seismic source model considers two major faults that are known historically to generate large earthquakes: the Sumatran subduction zone and the Sumatran transform fault. Several published studies were used to describe earthquakes along these faults during historical and pre-historical time, as well as to identify segmentation models of faults. Peak horizontal ground accelerations were calculated using ground motion prediction relations that were developed from seismic data obtained from the crustal interplate environment, crustal intraplate environment, along the subduction zone interface, and from deep intraslab earthquakes. Most of these relations, however, have not been developed for large distances that are needed for calculating the hazard across the Malaysian peninsula, and none were developed for earthquake ground motions generated in an interplate tectonic environment that are propagated into an intraplate tectonic environment. For the interplate and intraplate crustal earthquakes, we have applied ground-motion prediction relations that are consistent with California (interplate) and India (intraplate) strong motion data that we collected for distances beyond 200 km. For the subduction zone equations, we recognized that the published relationships at large distances were not consistent with global earthquake data that we collected and modified the relations to be compatible with the global subduction zone ground motions. In this analysis, we have used alternative source and attenuation models and weighted them to account for our uncertainty in which model is most appropriate for Sumatra or for the Malaysian peninsula. The resulting peak horizontal ground accelerations for 2% probability of exceedance in 50 years range from over 100% g to about 10% g across Sumatra and generally less than 20% g across most of the Malaysian peninsula. The ground motions at 10% probability of exceedance in 50 years are typically about 60% of the ground motions derived for a hazard level at 2% probability of exceedance in 50 years. The largest contributors to hazard are from the Sumatran faults.     

辽宁海城及其邻区地震b值空间分布特征

震级-频度分布（FMD）是地震学研究中最重要的经验公式之一,相关系数b是构造学和地震危险性评估的重要因子,具有表征前震和余震的特性。辽宁省地震多发生在金州断裂附近,自1975年海城7.3级地震发生后,与金州断裂交汇的海城河—大洋河断裂开启活跃模式,其东南端岫岩附近在1999年又发生5.6级地震。近年来盖州附近地震活动也在增强。因此,本文利用b值空间分布特征对海城及其邻区的应力分布特点进行研究。震源定位准确与否直接影响b值计算,双差定位后的数据与常规目录相比具有更高的精度,但是完整性有一定下降。本文收集了中国地震台网1981—2005年的辽宁省地震目录,并进行双差定位,比较分析了常规目录数据和双差数据的b值分布差异,认为在地震密集区,双差定位后的数据可以被用来获得更准确的b值。对主要研究区进行网格划分,使用双差数据,得到b值的水平和垂直分布特征。结果表明：b值为0.6~1.8,随深度增加而降低;岫岩和盖州震区具有较低的b值,意味着具有较高的地震危险性;浑河震区与海城河—大洋河断裂东南方向具有较高的b值,说明该区域未来发生大地震的概率很低;与金州断裂交汇区域的b值在1.0附近,说明该地区应力暂时处于稳定状态,未来具有较低的地震危险性。     

太阳活动周与全球大震

为研究太阳活动与全球大震的关系, 引入一个无量纲的"地震能量函数√G", 并分析研究了1681—2011故年间全球M≥7.0大震的能量释放的时间序列.由此发现全球大震在太阳活动周4个阶段的分布和活动度, 随震级的强度而异.提出地壳对太阳风暴加卸载响应模式, 用于解释此现象: 通过考察最近331 a, 得出全球共发生了10个M≥9.0超级巨震的时空分布特征, 特别是太阳活动峰年期间没有发生过超级巨震.该研究结果可为判断全球大震提供参考.      

地震成因综述

本文从地质、地球物理、地球化学和能量等方面分析了地震的成因。源于地核地幔的流体携带大量热能,为岩浆起源、地震形成和地热田提供了充足的能量,然而岩石聚集的应变能不足以产生中等以上的地震。大地震(M≥6.0)绝大部分分布在海沟、火山岛弧和大陆裂谷带等拉张性构造带,如环太平洋海沟、东印度洋海沟、大洋中脊、非洲裂谷、地中海-黑海-里海-波斯湾、欧亚大陆中部的伊塞克湖-阿拉湖-乌布苏湖-库苏古尔湖-贝加尔湖裂谷。流体在地球深部物质运动、地壳运动、地震和火山活动中扮演着重要作用。全球到处发育的隐爆角砾岩表明隐爆作用的普遍性。深部流体向上运移、向地表逃逸的过程中发生爆炸,在地球内部产生了不同震级和震源深度的地震。因此,隐爆应该是产生地震的主要机制。地震成因的隐爆模型不仅能够更好地解释不连续、各向异性的非弹性介质中发生的各类地震,譬如中深源震、震群、慢地震和非双力偶性地震等,而且能够更好地诠释全球地震、火山和地热带在空间上的吻合以及隐爆角砾岩型矿藏的形成。     

扰动强度递增斜坡动力学演变规律的振动台试验

“5.12”汶川大地震和“4.20”芦山地震均触发了大量的崩塌、滑坡。实震资料显示,不同地震烈度区地震触发崩塌滑坡规模的整体分布规律会发生变化。这一统计层面的认知亟待得到物理试验的验证。在自组织临界状态理论的概念框架下,开展了振动台砂堆模型试验。试验表明：输入地震波峰值加速度（PGA）为0.075g～0.125g时,落砂量与发生频率的关系可用幂律描述;PGA增加到0.15g～0.25g时,该关系服从对数正态分布;PGA增加到0.35g～0.45g时,该关系具有正态分布特征。元胞自动机模拟试验结果表明,随扰动强度增加,砂堆模型的动力学特性也经历了幂律-幂律弱化-正态分布的演变过程。按照物理学中的普适性原理,汶川、芦山地震Ⅸ度区崩塌滑坡规模与出现频率之间所呈现负幂律分布的现象,以及汶川地震Ⅺ度区所呈现的对数正态分布,可能是具有普适性意义的规律。这些认识可望为不同烈度区地震触发崩塌滑坡灾势预测提供科学依据。     

Spatial-temporal variability of seismic hazard in peninsular India

This paper examines the variability of seismic activity observed in the case of different geological zones of peninsular India (10°N–26°N; 68°E–90°E) based on earthquake catalog between the period 1842 and 2002 and estimates earthquake hazard for the region. With compilation of earthquake catalog in terms of moment magnitude and establishing broad completeness criteria, we derive the seismicity parameters for each geologic zone of peninsular India using maximum likelihood procedure. The estimated parameters provide the basis for understanding the historical seismicity associated with different geological zones of peninsular India and also provide important inputs for future seismic hazard estimation studies in the region. Based on present investigation, it is clear that earthquake recurrence activity in various geologic zones of peninsular India is distinct and varies considerably between its cratonic and rifting zones. The study identifies the likely hazards due to the possibility of moderate to large earthquakes in peninsular India and also presents the influence of spatial rate variation in the seismic activity of this region. This paper presents the influence of source zone characterization and recurrence rate variation pattern on the maximum earthquake magnitude estimation. The results presented in the paper provide a useful basis for probabilistic seismic hazard studies and microzonation studies in peninsular India.     

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.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

Multifractality in seismic sequences of NW Himalaya

Multifractal behaviour of interevent time sequences is investigated for the earthquake events in the NW Himalaya, which is one of the most seismically active zones of India and experienced moderate to large damaging earthquakes in the past. In the present study, the multifractal detrended fluctuation analysis (MF-DFA) is used to understand the multifractal behaviour of the earthquake data. For this purpose, a complete and homogeneous earthquake catalogue of the period 1965–2013 with a magnitude of completeness M _w 4.3 is used. The analysis revealed the presence of multifractal behaviour and sharp changes near the occurrence of three earthquakes of magnitude (M _w ) greater than 6.6 including the October 2005, Muzaffarabad–Kashmir earthquake. The multifractal spectrum and related parameters are explored to understand the time dynamics and clustering of the events.

     

Seismic hazard map of Coimbatore using subsurface fault rupture

This study presents the future seismic hazard map of Coimbatore city, India, by considering rupture phenomenon. Seismotectonic map for Coimbatore has been generated using past earthquakes and seismic sources within 300 km radius around the city. The region experienced a largest earthquake of moment magnitude 6.3 in 1900. Available earthquakes are divided into two categories: one includes events having moment magnitude of 5.0 and above, i.e., damaging earthquakes in the region and the other includes the remaining, i.e., minor earthquakes. Subsurface rupture character of the region has been established by considering the damaging earthquakes and total length of seismic source. Magnitudes of each source are estimated by assuming the subsurface rupture length in terms of percentage of total length of sources and matched with reported earthquake. Estimated magnitudes match well with the reported earthquakes for a RLD of 5.2% of the total length of source. Zone of influence circles is also marked in the seismotectonic map by considering subsurface rupture length of fault associated with these earthquakes. As earthquakes relive strain energy that builds up on faults, it is assumed that all the earthquakes close to damaging earthquake have released the entire strain energy and it would take some time for the rebuilding of strain energy to cause a similar earthquake in the same location/fault. Area free from influence circles has potential for future earthquake, if there is seismogenic source and minor earthquake in the last 20 years. Based on this rupture phenomenon, eight probable locations have been identified and these locations might have the potential for the future earthquakes. Characteristic earthquake moment magnitude (M _w ) of 6.4 is estimated for the seismic study area considering seismic sources close to probable zones and 15% increased regional rupture character. The city is divided into several grid points at spacing of 0.01° and the peak ground acceleration (PGA) due to each probable earthquake is calculated at every grid point in city by using the regional attenuation model. The maximum of all these eight PGAs is taken for each grid point and the final PGA map is arrived. This map is compared to the PGA map developed based on the conventional deterministic seismic hazard analysis (DSHA) approach. The probable future rupture earthquakes gave less PGA than that of DSHA approach. The occurrence of any earthquake may be expected in near future in these eight zones, as these eight places have been experiencing minor earthquakes and are located in well-defined seismogenic sources.     

Statistical aspects of Parkfield earthquake sequence and Parkfield prediction experiment

This note discusses three interconnected statistical problems concerning the Parkfield sequence of moderate earthquakes and the Parkfield prediction experiment: (a) Is it possible that the quasi-periodic Parkfield sequence of characteristic earthquakes is no uncommon, specific phenomenon (the research hypothesis), but can be explained by a preferential selection from available earthquake catalogs? To this end we formulate the null hypothesis (earthquakes occur according to the Poisson process in time and their size follows the Gutenberg-Richter relation). We test whether the null hypothesis can be rejected as an explanation for the Parkfield sequence. (b) If the null hypothesis cannot be refuted, what is the probability of magnitude m ≥ 6 earthquake occurrence in the Parkfield region? (c) The direct goal of the Parkfield experiment is the registration of precursory phenomena prior to a m6 earthquake. However, in the absence of the characteristic earthquake, can the experiment resolve which of the two competing hypotheses is true in a reasonable time? Statistical analysis is hindered by an insufficiently rigorous definition of the research model and inadequate or ambiguous data. However, we show that the null hypothesis cannot be decisively rejected. The quasi-periodic pattern of intermediate size earthquakes in the Parkfield area is a statistical event likely to occur by chance if it has been preferentially selected from available earthquake catalogs. The observed magnitude-frequency curves for small and intermediate earthquakes in the Parkfield area agree with the theoretical distribution computed on the basis of a modified Gutenberg-Richter law (gamma distribution), using deformation rates for the San Andreas fault. We show that the size distribution of the Parkfield characteristic earthquakes can also be attributed to selection bias. According to the null hypothesis, the yearly probability of a m ≥ 6 earthquake originating in the Parkfield area is less than 1%, signifying that several more decades of observation may be needed before the expected event occurs. By its design, the Parkfield experiment cannot be expected to yield statistically significant conclusions on the validity of the research hypothesis for many decades.     

Seismicity assessment in and around Syria based on instrumental data: application of Gumbel distributions and Gutenberg-Richter relationship

An instrumental earthquake catalog covering the time span between 1903 and 2007 and for the area bounded by 32°N–38°N and 35°E–43°E has been compiled in this research. The catalog has a magnitude of completeness (M _c ) with 3.5. Least squares and statistical probability Gumbel’s techniques with different approaches have been applied on the instrumental events in order to assess the average recurrence time periods for different earthquake magnitudes. The constants a and b of Gutenberg-Richter and the average recurrence times have been computed firstly for the study area and secondly for the central and northern parts of Dead Sea fault system. The different statistical computations using Knopoff and Kagan formalism are generally in agreement and suggest an average recurrence time of 203 years for an earthquake of magnitude 7 for the region. The occurrence of large well-documented historical earthquakes in Lebanon and western Syria, the existence of active fault segments, the absence of large earthquakes during the study period, the increasing number of the low-magnitude earthquakes, and the continued accumulation of the strain since 1900 indicate therefore the probability of an earthquake occurrence of a large magnitude. This should be permanently taken into consideration in seismic hazard assessment on the local and regional scales.