Macroseismic Intensity Modelling of Earthquakes in Southern Spain: Attenuation and Tectonic Implications

A modelling of the observed macroseismic intensity of historical and instrumental earthquakes in southern Spain is proposed, with the aim of determining the macroseismic parameters for seismic hazard evaluation in a region in which the characterization of intensity distribution of seismic events shows different levels of difficulty referable to the complex faults system of the area in study. The adopted procedure allows an analytical determination of epicenters and principal attenuation directions of earthquakes with a double level of verification with reference to the maximum shaking area and structural lineaments of the region, respectively. The analyses, carried out on a suitable number of events, highlight, therefore, some elements for a preliminary characterization of a seismic zonation on the basis of the consistency between seismic intensity distribution of earthquakes and corresponding structural framework.     

Macroseismic field of the October 27, 2004 Vrancea (Romania) moderate subcrustal earthquake

On October 27, 2004, a moderate size earthquake occurred in the Vrancea seismogenic region (Romania). The Vrancea seismic zone is an area of concentrated seismicity at intermediate depths beneath the bending area of the southeastern Carpathians. The 2004 M _w?=?6 Vrancea subcrustal earthquake is the largest seismic event recorded in Romania since the 1990 earthquakes. With a maximum macroseismic intensity of VII Medvedev–Sponheuer–Kárník (MSK-64) scale, the seismic event was felt to a distance of 600 km from the epicentre. This earthquake caused no serious damage and human injuries. The main purpose of this paper is to present the macroseismic map of the earthquake based on the MSK-64 intensity scale. After the evaluation of the macroseismic effects of this earthquake, an intensity dataset has been obtained for 475 sites in the Romanian territory. Also, the maximum horizontal accelerations recorded in the area by the K2 network are compared to the intensity values.     

An Anisotropic Attenuation Law of Macroseismic Intensity Performed on Virtual Intensity Distribution of Seismogenic Zones

An anisotropic attenuation law, based on an anisotropic characterization of intensity distribution for seismogenic zones, is proposed. This approach, that distinguishes itself for its consistency to the observed data, initially reconfigured by filtering procedures, is particularly suitable for seismic hazard evaluation.     

Seismic hazard assessment in the volcanic region of Mt. Etna (Italy): a probabilistic approach based on macroseismic data applied to volcano-tectonic seismicity

In the framework of the UPStrat-MAFA project, a seismic hazard assessment has been undertaken in the volcanic region of Mt. Etna as a first step in studies aimed at evaluating the risk on an urban scale. The analysis has been carried out with the SASHA code which uses macroseismic data in order to calculate, starting from the site seismic history, the maximum intensity value expected in a given site with a probability of exceedance of 10 % (I_ref), for a fixed exposure time. Depending on the aims of the project, hazard is estimated for local volcano-tectonic seismicity and short exposure times (10 and 30 years), without taking into account the contribution of “regional” events characterized by much longer recurrence times. Results from tasks A, B and D of the project have produced an updated macroseismic dataset, better performing attenuation models and new tools for SASHA, respectively. The maps obtained indicate that the eastern flank of Etna, the most urbanized sector of the volcano, is characterized by a high level of hazard with I_ref values up to degree VIII EMS, and even IX EMS locally. The disaggregated data analysis allows recognizing the “design earthquake” and the seismogenic fault which most contribute to the hazard at a site-scale. The latter analysis is the starting point to select the scenario earthquake to be used in the analyses of tasks C and F of the project dealing with, respectively, synthetic ground motion simulations and the evaluation of the Disruption Index.     

Seismic hazard assessment in Central Asia: Outcomes from a site approach

In the process of updating existing PSHA maps in Central Asia, a first step is the evaluation of the seismic hazard in terms of macroseismic intensity by applying a data driven method. Following the Site Approach to Seismic Hazard Assessment (SASHA) [11], the evaluation of the probability of exceedance of any given intensity value over a fixed exposure time, is mainly based on the seismic histories available at different locations without requiring any a-priori assumption about seismic zonation. The effects of earthquakes not included in the seismic history can be accounted by propagating the epicentral information through a Intensity Prediction Equation developed for the analyzed area. In order to comply with existing building codes in the region that use macroseismic intensity instead of PGA, we evaluated the seismic hazard at 2911 localities using a macroseismic catalog composed by 5322 intensity data points relevant to 75 earthquakes in the magnitude range 4.6–8.3. The results show that for most of the investigated area the intensity having a probability of at least 10% to be exceeded in 50 years is VIII. The intensity rises to IX for some area struck by strong earthquakes in the past, like the Chou-Kemin-Chilik fault zone in northern Tien-Shan, between Kyrgyzstan and Kazakhstan, or in Gissar range between Tajikistan and Uzbekistan. These values are about one intensity unit less than those evaluated in the Global Seismic Hazard Assessment Program (GSHAP; Ulomov, The GSHAP Region 7 working group [29]). Moreover, hazard curves have been extracted for the main towns of Central Asia and the results compared with the estimates previously obtained. A good agreement has been found for Bishkek (Kyrgyzstan) and Dushanbe (Tajikistan), while a lower probability of occurrence of I=VIII has been obtained for Tashkent (Uzbekistan) and a larger one for I=IX in Almaty (Kazakhstan).     

根据烈度分布确定华北历史地震破裂区的经验准则及其应用

地震破裂区是地震时沿发震断裂带的同震错动面或破裂面在地表的垂直投影区域,指示了震源断层/破裂的位置与尺度。确定过去长期的强震/大地震破裂区是鉴别地震空区、研究与预测强震危险性的重要基础。对于现代强震,破裂区可运用多种现代技术方法确定,但对于历史强震,破裂区确定的方法需要探索与发展。以华北地区为例,研究利用烈度/等震线资料、结合地震构造与震区地表地质环境等信息确定历史强震破裂区的方法,并开展应用试验。结果表明:研究区现代地震破裂区延伸的烈度区间与极震区烈度、震区环境之间存在密切关系,基于这种关系建立了2条经验准则,可分别用于根据烈度分布确定华北2类震区环境(基岩区和厚层第四纪松散堆积覆盖区)历史强震破裂区的位置与延伸。文中还提出通过综合地震构造、现代小震/余震分布等信息,辅助确定历史强震破裂区横向宽度的思路与途径。作为应用试验,文中确定了5次历史地震的破裂区,结果表明本文发展的经验准则及相应方法适用于华北地区历史强震破裂区的确定。     

Sensitivity analysis on the input parameters in probabilistic seismic hazard assessment

A project has been implemented in recent years for assessing seismic hazard in the Italian territory on probabilistic bases, to be used as scientific background for the revision of the current seismic zonation. A consolidated approach was considered for the purpose; seismic hazard was estimated in terms of peak ground acceleration and macroseismic intensity. As the computer code employed allows the user to make specific choices on some input data, some rather unorthodox decisions were taken regarding earthquake catalogue completeness, seismicity rates, boundaries of the seismogenic zones, definition of the maximum magnitude, attenuation relation, etc. The overwhelming amount of geological and seismological data for Italy (just consider, for example, that the earthquake catalogue collects events which occurred over the last ten centuries) permits the operator to make different choices, more or less cautiously. It is quite interesting, then, to evaluate the influence of the specific choices on the final hazard results as a comparison to traditional possibilities. The tests performed clearly indicate the critical choices and quantify their contribution. In particular, we consider thorough comprehension of the space geometry of the earthquake source boundaries and the adequacy of the attenuation relation in modelling the radiation pattern very important.     

Probabilistic seismic hazard assessment in SE-Spain based on macroseismic site histories

In many countries such as Spain earthquake databases still mainly comprise macroseismic data from felt effects. The full exploit of this information is of basic importance for seismic risk assessment and emergency planning, given the strict link between macroseismic intensity and damage. A probabilistic procedure specifically developed to handle macroseismic data, mostly relying on site information and seismogenic-source free, has been applied to evaluate seismic hazard in SE-Spain (Alicante-Murcia region). Present seismicity is moderate-low with largest magnitudes slightly over Mw5.0. The historical record includes very destructive earthquakes, maximum EMS98 intensities reaching IX–X and X in the nineteenth century (e.g., Torrevieja 1829 earthquake). Very recently, two events in the area on 11 May 2011 (Mw4.5, Mw5.2) killed nine people, injured 300, and produced important damage in the city of Lorca. Regional hazard maps for the area together with specific hazard curves at selected localities are obtained. Results are compared with the maximum observed intensities in the period 1300–2012, and with the values in the seismic hazard map from the Spanish Building Code in force. In general, the maximum felt intensity values are closer to the hazard values calculated for 2 % probability of exceedance in 50 years, using felt and expected intensity. The intensity-based probabilistic hazard maps obtained through the applied approach reduce the inherent smoothing of those based on standard probabilistic seismic hazard assessment approaches for the region, allowing identifying possible over- or sub-estimates of site hazard values, providing very valuable information for risk reduction strategies or for future updates of the building code hazard maps.     

面状发震构造在地震构造法中的应用：以大姚—姚安发震构造鉴定为例

核电厂地震安全性评价中的地震构造法,要求鉴定发震构造和划分地震构造区,在以往实践中,发震构造鉴定往往基于地表活动断裂构造,且表征为线状震源.当存在较强非随机分布的地震活动且难以找到清晰的地表活动断裂构造形迹时,地震构造法就难以合理地表现这些地震的危险性.本文以云南滇中大姚—姚安发震构造鉴定为例,探讨了在地表活动构造形迹不清,中强地震活动性较强的滇中大姚—姚安地区,采用面状发震构造来表征地震危险性的方法,讨论了在地震构造法中采用面状发震构造的必要性、鉴定思路和方法,并建议在今后的核工程地震危险性评价地震构造法中应充分考虑面状发震构造的应用.     

Assessment of macroseismic intensity in the Nile basin,Egypt

This work intends to assess deterministic seismic hazard and risk analysis in terms of the maximum expected intensity map of the Egyptian Nile basin sector. Seismic source zone model of Egypt was delineated based on updated compatible earthquake catalog in 2015, focal mechanisms, and the common tectonic elements. Four effective seismic source zones were identified along the Nile basin. The observed macroseismic intensity data along the basin was used to develop intensity prediction equation defined in terms of moment magnitude. Expected maximum intensity map was proven based on the developed intensity prediction equation, identified effective seismic source zones, and maximum expected magnitude for each zone along the basin. The earthquake hazard and risk analysis was discussed and analyzed in view of the maximum expected moment magnitude and the maximum expected intensity values for each effective source zone. Moderate expected magnitudes are expected to put high risk at Cairo and Aswan regions. The results of this study could be a recommendation for the planners in charge to mitigate the seismic risk at these strategic zones of Egypt.     

Seismic hazard assessment for the Itoiz dam site (Western Pyrenees,Spain)

This work summarises the seismic hazard analysis performed for the complete characterisation of strong ground-motion at the site of the Itoiz dam (Western Pyrenees, Spain). The hazard analysis includes the compilation of a composite catalogue from French and Spanish agencies, the definition of an original hybrid seismogenic source model (including zones and major faults) and the selection of ground motion prediction equations (GMPEs). Hazard results are provided as hazard curves and acceleration response spectra on rock for the 1000- and 5000-year return periods, which correspond respectively to the operating basis earthquake (OBE) and safety evaluation earthquake (SEE). The impact of truncating GMPEs at a number of standard deviations (epsilon) has been found not critical here for the return periods targeted. Subsequently, an analysis of the contribution of each source to total hazard and a hazard disaggregation analysis are performed in order to establish the earthquake-source parameters for both the OBE and SEE scenarios consistently with the seismotectonics of the region. The European Strong Motion database is then searched and a selection of records is proposed for each of the scenarios. Our results suggest that seismic hazard in the region is underestimated by the official Spanish seismic hazard map included in the current version of the code (NCSE-02), which is the reference document for the definition of seismic actions for dam projects in the whole Pyrenees.     

A Probabilistic Assessment of Earthquake Hazard Parameters in NW Himalaya and the Adjoining Regions

The maximum likelihood estimation method is applied to study the geographical distribution of earthquake hazard parameters and seismicity in 28 seismogenic source zones of NW Himalaya and the adjoining regions. For this purpose, we have prepared a reliable, homogeneous and complete earthquake catalogue during the period 1500–2010. The technique used here allows the data to contain either historical or instrumental era or even a combination of the both. In this study, the earthquake hazard parameters, which include maximum regional magnitude (M _max), mean seismic activity rate (λ), the parameter b (or β?=?b/log e) of Gutenberg–Richter (G–R) frequency-magnitude relationship, the return periods of earthquakes with a certain threshold magnitude along with their probabilities of occurrences have been calculated using only instrumental earthquake data during the period 1900–2010. The uncertainties in magnitude have been also taken into consideration during the calculation of hazard parameters. The earthquake hazard in the whole NW Himalaya region has been calculated in 28 seismogenic source zones delineated on the basis of seismicity level, tectonics and focal mechanism. The annual probability of exceedance of earthquake (activity rate) of certain magnitude is also calculated for all seismogenic source zones. The obtained earthquake hazard parameters were geographically distributed in all 28 seismogenic source zones to analyze the spatial variation of localized seismicity parameters. It is observed that seismic hazard level is high in Quetta-Kirthar-Sulaiman region in Pakistan, Hindukush-Pamir Himalaya region and Uttarkashi-Chamoli region in Himalayan Frontal Thrust belt. The source zones that are expected to have maximum regional magnitude (M _max) of more than 8.0 are Quetta, southern Pamir, Caucasus and Kashmir-Himanchal Pradesh which have experienced such magnitude of earthquakes in the past. It is observed that seismic hazard level varies spatially from one zone to another which suggests that the examined regions have high crustal heterogeneity and seismotectonic complexity.     

雄安新区地震危险性评估

地震危险性评估对确定工程抗震设防等级、制定城市规划与减轻地震灾害等具有重要意义.由于传统分析方法存在诸多缺陷,用其对特定区域进行地震危险性评估可靠性差.本文基于孕震断层多锁固段脆性破裂理论,提出了利用岩石破裂自相似性预测研究区标志性预震的方法,并进而创立了一种以地震物理预测为基础的地震危险性评估新方法.雄安新区位于唐山地震区内,近邻运城地震区.基于唐山和运城地震区地震趋势分析结果及主要断裂展布与雄安新区位置关系,在唐山地震区内划分了保定研究区,在运城地震区内划分了行唐和十渡研究区.采用上述方法,预测了这三个区将发生的较大预震.根据某些学者提出的地震烈度经验公式,评估了唐山地震区发生标志性地震、预震及运城地震区发生预震导致的雄安新区地震烈度.结果表明,未来50年内,雄安新区抗震设防烈度从原7度调整为8度为宜.     

K-means cluster analysis and seismicity partitioning for Pakistan

Pakistan and the western Himalaya is a region of high seismic activity located at the triple junction between the Arabian, Eurasian and Indian plates. Four devastating earthquakes have resulted in significant numbers of fatalities in Pakistan and the surrounding region in the past century (Quetta, 1935; Makran, 1945; Pattan, 1974 and the recent 2005 Kashmir earthquake). It is therefore necessary to develop an understanding of the spatial distribution of seismicity and the potential seismogenic sources across the region. This forms an important basis for the calculation of seismic hazard; a crucial input in seismic design codes needed to begin to effectively mitigate the high earthquake risk in Pakistan. The development of seismogenic source zones for seismic hazard analysis is driven by both geological and seismotectonic inputs. Despite the many developments in seismic hazard in recent decades, the manner in which seismotectonic information feeds the definition of the seismic source can, in many parts of the world including Pakistan and the surrounding regions, remain a subjective process driven primarily by expert judgment. Whilst much research is ongoing to map and characterise active faults in Pakistan, knowledge of the seismogenic properties of the active faults is still incomplete in much of the region. Consequently, seismicity, both historical and instrumental, remains a primary guide to the seismogenic sources of Pakistan. This study utilises a cluster analysis approach for the purposes of identifying spatial differences in seismicity, which can be utilised to form a basis for delineating seismogenic source regions. An effort is made to examine seismicity partitioning for Pakistan with respect to earthquake database, seismic cluster analysis and seismic partitions in a seismic hazard context. A magnitude homogenous earthquake catalogue has been compiled using various available earthquake data. The earthquake catalogue covers a time span from 1930 to 2007 and an area from 23.00° to 39.00°N and 59.00° to 80.00°E. A threshold magnitude of 5.2 is considered for K-means cluster analysis. The current study uses the traditional metrics of cluster quality, in addition to a seismic hazard contextual metric to attempt to constrain the preferred number of clusters found in the data. The spatial distribution of earthquakes from the catalogue was used to define the seismic clusters for Pakistan, which can be used further in the process of defining seismogenic sources and corresponding earthquake recurrence models for estimates of seismic hazard and risk in Pakistan. Consideration of the different approaches to cluster validation in a seismic hazard context suggests that Pakistan may be divided into K?=?19 seismic clusters, including some portions of the neighbouring countries of Afghanistan, Tajikistan and India.     

Intensity attenuation in the Campania region,Southern Italy

Ground motion prediction equations (GMPE) in terms of macroseismic intensity are a prerequisite for intensity-based shake maps and seismic hazard assessment and have the advantage of direct relation to earthquake damage and good data availability also for historical events. In this study, we derive GMPE for macroseismic intensity for the Campania region in southern Italy. This region is highly exposed to the seismic hazard related to the high seismicity with moderate- to large-magnitude earthquakes in the Appenninic belt. The relations are based on physical considerations and are easy to implement for the user. The uncertainties in earthquake source parameters are accounted for through a Monte Carlo approach and results are compared to those obtained through a standard regression scheme. One relation takes into account the finite dimensions of the fault plane and describes the site intensity as a function of Joyner–Boore distance. Additionally, a relation describing the intensity as a function of epicentral distance is derived for implementation in cases where the dimensions of the fault plane are unknown. The relations are based on an extensive dataset of macroseismic intensities for large earthquakes in the Campania region and are valid in the magnitude range M _w = 6.3–7.0 for shallow crustal earthquakes. Results indicate that the uncertainties in earthquake source parameters are negligible in comparison to the spread in the intensity data. The GMPE provide a good overall fit to historical earthquakes in the region and can provide the intensities for a future earthquake within 1 intensity unit.     

Seismic hazard assessment for Iceland in terms of macroseismic intensity using a site approach

We present the results of probabilistic seismic hazard assessment for Iceland in the framework of the EU project UPStrat-MAFA using the so-called site approach implemented in the SASHA computational code. This approach estimates seismic hazard in terms of macroseismic intensity by basically relying on local information about documented effects of past seismic events in the framework of a formally coherent and complete treatment of intensity data. In the case of Iceland, due to the lack of observed intensities for past earthquakes, local seismic histories were built using indirect macroseismic estimates deduced from epicentral information through an empirical attenuation relationship in probabilistic form. Seismic hazard was computed for four exceedance probabilities for an exposure time of 50 years, equivalent to average return periods of 50, 200, 475 and 975 years. For some localities, further return periods were examined and deaggregation analysis was performed. Results appear significantly different from previous seismic hazard maps, though just a semi-qualitative comparison is possible because of the different shaking measure considered (peak ground acceleration versus intensity), and the different computational methodology and input data used in these studies.     

Probabilistic seismic hazard assessment for Romania and sensitivity analysis: A case of joint consideration of intermediate-depth (Vrancea) and shallow (crustal) seismicity

The earthquake risk on Romania is one of the highest in Europe, and seismic hazard for almost half of the territory of Romania is determined by the Vrancea seismic region, which is situated beneath the southern Carpathian Arc. The region is characterized by a high rate of occurrence of large earthquakes in a narrow focal volume at depth from 70 to 160 km. Besides the Vrancea area, several zones of shallow seismicity located within and outside the Romanian territory are considered as seismically dangerous. We present the results of probabilistic seismic hazard analysis, which implemented the “logic tree” approach, and which considered both the intermediate-depth and the shallow seismicity. Various available models of seismicity and ground-motion attenuation were used as the alternative variants. Seismic hazard in terms of macroseismic intensities, peak ground acceleration, and response spectra was evaluated for various return periods. Sensitivity study was performed to analyze the impact of variation of input parameters on the hazard results. The uncertainty on hazard estimates may be reduced by better understanding of parameters of the Vrancea source zone and the zones of crustal seismicity. Reduction of uncertainty associated with the ground-motion models is also very important issue for Romania.     

Geodetic and geological evidence of active tectonics in south-western Sicily (Italy)

Integrated geological, geodetic and marine geophysical data provide evidence of active deformation in south-western Sicily, in an area spatially coincident with the macroseismic zone of the destructive 1968 Belice earthquake sequence. Even though the sequence represents the strongest seismic event recorded in Western Sicily in historical times, focal solutions provided by different authors are inconclusive on possible faulting mechanism, which ranges from thrusting to transpression, and the seismogenic source is still undefined. Interferometric (DInSAR) observations reveal a differential ground motion on a SW–NE alignment between Campobello di Mazara and Castelvetrano (CCA), located just west of the maximum macroseismic sector. In addition, new GPS campaign-mode data acquired across the CCA alignment documents NW–SE contractional strain accumulation. Morphostructural analysis allowed to associate the alignment detected through geodetic measurements with a topographic offset of Pleistocene marine sediments. The on-land data were complemented by new high-resolution marine geophysical surveys, which indicate recent contraction on the offshore extension of the CCA alignment. The discovery of archaeological remains displaced by a thrust fault associated with the alignment provided the first likely surface evidence of coseismic and/or aseismic deformation related to a seismogenic source in the area. Results of the integrated study supports the contention that oblique thrusting and folding in response to NW–SE oriented contraction is still active. Although we are not able to associate the CCA alignment to the 1968 seismic sequence or to the historical earthquakes that destroyed the ancient Greek city of Selinunte, located on the nearby coastline, our result must be incorporated in the seismic hazard evaluation of this densely populated area of Sicily.     

Seismic Zonation of the Delhi Region for Bedrock Ground Motion

An endeavor is made to compute peak ground horizontal accelerations at bedrock level in the Delhi region due to the seismogenic sources present around Delhi. The entire area is divided into six seismogenic sources for which seismic hazard analysis is carried out using the complete and extreme part of the seismicity data. Maximum likelihood estimates of hazard parameters viz., seismic activity rate , b value and maximum probable earthquake M _max are made for each zone. The return periods and the probabilities of occurrence of various magnitudes for return periods of 50, 100 and 1000 years are also computed for each zone. The peak ground acceleration (PGA) values for 20% exceedance in 50 years are then computed for the Delhi region from each zone. The maximum PGA value considering all the zones is 0.34 g, which is due to the Mathura fault zone. The seismogenic zones V and VI, i.e., Mathura fault zone and the Sohna fault zone are observed to be contributing maximum PGA values in the Delhi region governing the isoacceleration contours computed for the region. The seismic zonation map for the PGA values at the bedrock level is obtained for the Delhi region. This can be used directly as input for the microzonation of ground motion at the surface by incorporating the local site conditions.