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.     

Validation through HVSR measurements of a method for the quick detection of site amplification effects from intensity data: an application to a seismic area in Northern Italy

An instrumental validation is attempted of an innovative approach devoted to the quick individuation, from macroseismic data, of site amplification phenomena able to significantly modify seismic hazard levels expected on the basis of average propagation effects only. According to this methodology, two evaluations of hazard are performed at each investigated locality: the former, obtained by epicentral intensity data ‘reduced’ at the site through a probabilistic attenuation function and, the latter, computed by integrating such data with seismic effects actually observed at the site during past earthquakes. The comparison, for each locality, between these two hazard estimates allow to orientate the identification of those sites where local amplifications of earthquake ground motion could be significant. In order to check such methodology, indications obtained in this way from macroseismic data are compared with the estimates of transfer functions performed through the HVSR technique applied to microtremors. Results concerning municipalities located in a seismic area of Northern Italy indicate a good agreement between macroseismic and instrumental estimates.     

Damage to residential buildings in Hveragerði during the 2008 Ölfus Earthquake: simulated and surveyed results

This study analyses the performance of residential buildings in the town of Hveragerði in South Iceland during the 29 May 2008 Mw 6.3 Ölfus Earthquake. The earthquake occurred very close to the town, approximately 3–4 km from it. Ground shaking caused by the earthquake was recorded by a dense strong-motion array in the town. The array provided high-quality three-component ground acceleration data which is used to quantify a hazard scenario. In addition, surveys conducted in the town in the aftermath of the earthquake have provided information on macroseismic intensity at various locations in the town. Detailed information regarding the building stock in the town is collected, and their seismic vulnerability models are created by using building damage data obtained from the June 2000 South Iceland earthquakes. Damage to buildings are then simulated by using the scenario hazard and vulnerability models. Damage estimates were also obtained by conducting a survey. Simulated damage based on the scenario macroseismic intensity is found to be similar to damage estimated from survey data. The buildings performed very well during the earthquake—damage suffered was only 5 % of the insured value on the average. Correlation between actual damage and recorded ground-motion parameters is found to be statistically insignificant. No significant correlation of damage was observed, even with macroseismic intensity. Whereas significant correlation was observed between peak ground velocity and macroseismic intensity, neither of them appear to be good indicators of damage to buildings in the study area. This lack of correlation is partly due to good seismic capacity of buildings and partly due to the ordinal nature of macroseismic intensity scale. Consistent with experience from many past earthquakes, the survey results indicate that seismic risk in South Iceland is not so much due to collapse of buildings but rather due to damage to non-structural components and building contents.     

Probabilistic Microzonation of Urban Territories: A Case of Tashkent City

—?The problem of accounting for local soil effect on earthquake ground motion is especially urgent when assessing seismic hazard – recent needs of earthquake engineering require local site effects to be included into hazard maps. However, most recent works do not consider the variety of soil conditions or are performed for generalized site categories, such as “hard rock,”“soft soil” or “alluvium.” A technique of seismic hazard calculations on the basis of the Fourier Amplitude Spectra recently developed by the authors allows us to create hazard maps involving the influence of local soil conditions using soil/bedrock spectral ratios. Probabilistic microzoning maps may be constructed showing macroseismic intensity, peak ground acceleration, response and design spectra for various return periods (probability of exceedance), that allow optimization of engineering decisions. An application of this approach is presented which focused on the probabilistic microzoning of the Tashkent City.     

Application of SASHA to seismic hazard assessment for Portugal mainland

In the frame of the UPStrat-MAFA “Urban Disaster Prevention Strategies Using MAcroseismic Fields and FAult Sources” project, seismic hazard has been assessed in Portugal in terms of macroseismic intensity. Assessment has been performed by using a probabilistic approach based on the statistical analysis of local seismic histories (i.e., the record of seismic effects at each locality) performed by a new version of the SASHA code (D’Amico and Albarello in Res Lett 79(5):663–671, 2008) on purpose modified to account for this specific area of study. Local seismic histories are reconstructed by considering documented effects or indirect estimates deduced from epicentral information or numerical simulations. All these pieces of evidence are combined taking into account relevant uncertainty and statistical completeness of information locally available. Distribution of expected maximum intensity (i.e., the maximum intensity characterized by a fixed exceedance probability for a exposure time of 50 years) has been obtained and compared with the one deduced from alternative approaches.     

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.     

Seismic Hazard Estimates Using Ill-defined Macroseismic Data at Site

—?A new approach is proposed to the seismic hazard estimate based on documentary data concerning local history of seismic effects. The adopted methodology allows for the use of “poor” data, such as the macroseismic ones, within a formally coherent approach that permits overcoming a number of problems connected to the forcing of available information in the frame of “standard” methodologies calibrated on the use of instrumental data. The use of the proposed methodology allows full exploitation of all the available information (that for many towns in Italy covers several centuries) making possible a correct use of macroseismic data characterized by different levels of completeness and reliability. As an application of the proposed methodology, seismic hazard estimates are presented for two towns located in Northern Italy: Bologna and Carpi.     

Seismic hazard assessment in continental Ecuador

A seismic hazard assessment study of continental Ecuador is presented in this paper. The study begins with a revision of the available information on seismic events and the elaboration of a seismic catalog homogenized to magnitude Mw. Different seismic source definitions are revised and a new area-source model, based on geological and seismic data, is proposed. The available ground motion prediction equations for crustal and subduction sources are analyzed and selected for the tectonic environments observed in Ecuador. A probabilistic seismic hazard assessment approach is carried out to evaluate the exceedance probability of several levels of peak ground acceleration PGA and spectral accelerations SA (T) for periods (T) of 0.1, 0.2, 0.5, 1 and 2s. The resulting hazard maps for continental Ecuador are presented, together with the uniform hazard spectra of four province capital cities. Hazard disaggregation is carried out for target motions defined by the PGA values and SA (1s) expected for return periods of 475 and 2475 years, providing estimates for short-period and long-period controlling earthquakes.     

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).     

Comparison of conventional and Monte Carlo simulation-based probabilistic seismic hazard analyses for Shiraz city,southern Iran

Estimation of ground-motion amplitudes of different hazard levels is of paramount importance in planning of urban development of any metropolis. Such estimation can be computed through a probabilistic seismic hazard analysis (PSHA). This paper concentrates on the PSHA of an area located in Shiraz city, southern Iran. The area includes whole of Shiraz city (i.e., one of the largest and most populous cities of Iran) and its outskirts. Conventional and Monte Carlo simulation-based approaches are utilized to perform the PSHA of the studied area. Two areal seismic source models are delineated, and thence seismicity parameters of all zones associated with their corresponding uncertainties are computed. Uncertainties in ground-motion prediction are accounted for via three ground-motion prediction equations (GMPEs) within the logic tree framework. These GMPEs are applied to estimate bedrock ground shaking (V_s30?=?760 m/s) for several return periods (i.e., 75, 475, 975, and 2475 years). In general, the results of the two abovementioned PSHA approaches show relatively similar results. However, the Monte Carlo simulation-based approach overpredicts bedrock spectral accelerations at periods of 0.4–2.5 s compared to the conventional PSHA approach for return periods of 475, 975, and 2475 years.     

Do French macroseismic intensity observations agree with expectations from the European Seismic Hazard Model 2013?

Probabilistic seismic hazard assessments are the basis of modern seismic design codes. To test fully a seismic hazard curve at the return periods of interest for engineering would require many thousands of years’ worth of ground-motion recordings. Because strong-motion networks are often only a few decades old (e.g. in mainland France the first accelerometric network dates from the mid-1990s), data from such sensors can be used to test hazard estimates only at very short return periods. In this article, several hundreds of years of macroseismic intensity observations for mainland France are interpolated using a robust kriging-with-a-trend technique to establish the earthquake history of every French mainland municipality. At 24 selected cities representative of the French seismic context, the number of exceedances of intensities IV, V and VI is determined over time windows considered complete. After converting these intensities to peak ground accelerations using the global conversion equation of Caprio et al. (Ground motion to intensity conversion equations (GMICEs): a global relationship and evaluation of regional dependency, Bulletin of the Seismological Society of America 105:1476–1490, 2015), these exceedances are compared with those predicted by the European Seismic Hazard Model 2013 (ESHM13). In half of the cities, the number of observed exceedances for low intensities (IV and V) is within the range of predictions of ESHM13. In the other half of the cities, the number of observed exceedances is higher than the predictions of ESHM13. For intensity VI, the match is closer, but the comparison is less meaningful due to a scarcity of data. According to this study, the ESHM13 underestimates hazard in roughly half of France, even when taking into account the uncertainty in the conversion from intensity to acceleration. However, these results are valid only for the acceleration range tested in this study (0.01 to 0.09 g).     

The Umbria–Marche case: some suggestions for the Italian seismic zonation

Probabilistic seismic hazard maps, in terms of spectral acceleration and uniform hazard response spectra at given sites, considering local soil conditions, represent a much more complete estimate of the seismic hazard than the traditional maps in terms of peak ground acceleration or macroseismic intensity. This is particularly true when the requests of urban planners and engineers have to be met. The present analysis shows how some hazard parameters, such as the effective peak acceleration and the spectrum intensity, can well synthesise the overall information available from traditional probabilistic studies, but also suggests that soil condition is a first-order ingredient for effective seismic hazard mapping at national level. Three Italian towns, damaged by the 1997 Umbria–Marche earthquake sequence, are considered as example to demonstrate that: (1) soil condition dependent uniform hazard spectra well approximate actual spectra recorded during some events of the seismic sequence; (2) for these localities, the design spectrum of the present Italian seismic code does not seem adequate.

These considerations have induced the Italian scientific community to propose an updating of the national seismic zonation on the basis of several hazard parameters, that are described in this paper.     

Probabilistic seismic hazard assessment for South-Eastern France

The accurate evaluation and appropriate treatment of uncertainties is of primary importance in modern probabilistic seismic hazard assessment (PSHA). One of the objectives of the SIGMA project was to establish a framework to improve knowledge and data on two target regions characterized by low-to-moderate seismic activity. In this paper, for South-Eastern France, we present the final PSHA performed within the SIGMA project. A new earthquake catalogue for France covering instrumental and historical periods was used for the calculation of the magnitude-frequency distributions. The hazard model incorporates area sources, smoothed seismicity and a 3D faults model. A set of recently developed ground motion prediction equations (GMPEs) from global and regional data, evaluated as adequately representing the ground motion characteristics in the region, was used to calculate the hazard. The magnitude-frequency distributions, maximum magnitude, faults slip rate and style-of-faulting are considered as additional source of epistemic uncertainties. The hazard results for generic rock condition (Vs30 = 800 m/s) are displayed for 20 sites in terms of uniform hazard spectra at two return periods (475 years and 10,000 years). The contributions of the epistemic uncertainties in the ground motion characterizations and in the seismic source characterization to the total hazard uncertainties are analyzed. Finally, we compare the results with existing models developed at national scale in the framework of the first generation of models supporting the Eurocode 8 enforcement, (MEDD 2002 and AFPS06) and at the European scale (within the SHARE project), highlighting significant discrepancies at short return periods.     

A probabilistic seismic hazard assessment for the Turkish territory: part II—fault source and background seismicity model

Over the years, several local and regional seismic hazard studies have been conducted for the estimation of the seismic hazard in Turkey using different statistical processing tools for instrumental and historical earthquake data and modeling the geologic and tectonic characteristics of the region. Recently developed techniques, increased knowledge and improved databases brought the necessity to review the national active fault database and the compiled earthquake catalogue for the development of a national earthquake hazard map. A national earthquake strategy and action plan were conceived and accordingly with the collaboration of the several institutions and expert researchers, the Revision of Turkish Seismic Hazard Map Project (UDAP-Ç-13-06) was initiated, and finalized at the end of 2014. The scope of the project was confined to the revision of current national seismic hazard map, using the state of the art technologies and knowledge of the active fault, earthquake database, and ground motion prediction equations. The following two seismic source zonation models are developed for the probabilistic earthquake hazard analysis: (1) Area source model, (2) Fault and spatial smoothing seismic source model (FSBCK). In this study, we focus on the development and the characterization of the Fault Source model, the background spatially smoothed seismicity model and intrinsic uncertainty on the earthquake occurrence-rates-estimation. Finally, PSHA results obtained from the fault and spatial smoothed seismic source model are presented for 43, 72, 475 and 2475 years return periods (corresponding to 69, 50, 10, and 2% probability of exceedance in 50 years) for PGA and 5% damped spectral accelerations at 0.2 and 1.0 s.     

Estimates of site seismicity rates using ill-defined macroseismic data

A new approach to the problem of site seismic hazard analysis is proposed, based on intensity data affected by uncertainties. This approach takes into account the ordinal and discrete character of intensities, trying to avoid misleading results due to the assumption that intensity can be treated as a real number (continuous distribution estimators, attenuation relationships, etc.). The proposed formulation is based on the use of a distribution function describing, for each earthquake, the probability that site seismic effects can be described by each possible intensity value. In order to obtain site hazard estimates where local data are lacking, the dependence of this distribution function with the distance from the macroseismic epicenter and with epicentral intensity is examined. A methodology has been developed for the purpose of combining such probabilities and estimating site seismicity rates which takes into account the effect of uncertainties involved in this kind of analysis. An application of this approach is described and discussed.     

A framework for the evaluation of ground motion selection and modification procedures

This study develops a framework to evaluate ground motion selection and modification (GMSM) procedures. The context is probabilistic seismic demand analysis, where response history analyses of a given structure, using ground motions determined by a GMSM procedure, are performed in order to estimate the seismic demand hazard curve (SDHC) for the structure at a given site. Currently, a GMSM procedure is evaluated in this context by comparing several resulting estimates of the SDHC, each derived from a different definition of the conditioning intensity measure (IM). Using a simple case study, we demonstrate that conclusions from such an approach are not always definitive; therefore, an alternative approach is desirable. In the alternative proposed herein, all estimates of the SDHC from GMSM procedures are compared against a benchmark SDHC, under a common set of ground motion information. This benchmark SDHC is determined by incorporating a prediction model for the seismic demand into the probabilistic seismic hazard analysis calculations. To develop an understanding of why one GMSM procedure may provide more accurate estimates of the SDHC than another procedure, we identify the role of ‘IM sufficiency’ in the relationship between (i) bias in the SDHC estimate and (ii) ‘hazard consistency’ of the corresponding ground motions obtained from a GMSM procedure. Finally, we provide examples of how misleading conclusions may potentially be obtained from erroneous implementations of the proposed framework. Copyright © 2014 John Wiley & Sons, Ltd.     

Probabilistic earthquake hazard analysis for Cairo,Egypt

Cairo is the capital of Egypt and the largest city in the Arab world and Africa, and the sixteenth largest metropolitan area in the world. It was founded in the tenth century (969 ad) and is 1046 years old. It has long been a center of the region’s political and cultural life. Therefore, the earthquake risk assessment for Cairo has a great importance. The present work aims to analysis the earthquake hazard of Cairo as a key input’s element for the risk assessment. The regional seismotectonics setting shows that Cairo could be affected by both far- and near-field seismic sources. The seismic hazard of Cairo has been estimated using the probabilistic seismic hazard approach. The logic tree frame work was used during the calculations. Epistemic uncertainties were considered into account by using alternative seismotectonics models and alternative ground motion prediction equations. Seismic hazard values have been estimated within a grid of 0.1°?×?0.1 ° spacing for all of Cairo’s districts at different spectral periods and four return periods (224, 615, 1230, and 4745 years). Moreover, the uniform hazard spectra have been calculated at the same return periods. The pattern of the contour maps show that the highest values of the peak ground acceleration is concentrated in the eastern zone’s districts (e.g., El Nozha) and the lowest values at the northern and western zone’s districts (e.g., El Sharabiya and El Khalifa).     

Macroseismic Parameters of Seismogenic Zones of Calabria and Sicily for Seismic Hazard Evaluation

For the seismic hazard evaluation of the region including the southern Calabro-Peloritanian Arc and southeastern Sicily, the determination of the macroseismic virtual intensity distributions has been carried out, characteristic of the seismogenic zones that fall within the area in study, starting with the structural framework of the region and from the analysis of the observed intensity effected through suitable filters. The macroseismic parameters, derived from such virtual distributions and used for seismic hazard evaluation, are not only a reference for eventual subsequently deeper knowledge referable to the need for a better characterization of the reference modelling, but distinguish themselves as an essential instrument for the definition of seismic hazard scenarios correlated to seismic events that take place in single seismogenic zones.     

The Grandori model: Anisotropic formulation and its application in the plotting of seismic hazard maps

The seismic hazard for the Calabro-Sicilian area is evaluated using an anisotropic formulation of the Grandori attenuation law. For each macroseismic field two main directions are identified: minimum and maximum attenuation of the macroseismic intensity. The results of the investigation show that the anisotropic formulation improves the compatibility level of the model (with respect to the isotropic one) with the intensities observed and produces probabilistic expected intensities which compare favourably with the values of seismic history in the investigated area when the zonation defined by the Messina University research group was used.