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.     

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.     

Linking ground motion measurements and macroseismic observations in France: a case study based on accelerometric and macroseismic databases

Comparison between accelerometric and macroseismic observations is made for three M _w?=?4.5 earthquakes, which occurred in north-eastern France and south-western Germany in 2003 and 2004. Scalar and spectral instrumental parameters are processed from the accelerometric data recorded by nine accelerometric stations located between 29 and 180 km from the epicentres. Macroseismic data are based on French Internet reports. In addition to the single questionnaire intensity, analysis of the internal correlation between the encoded answers highlights four predominant fields of questions bearing different physical meanings: (1) “vibratory motions of small objects”, (2) “displacement and fall of objects”, (3) “acoustic noise” and (4) “personal feelings”. Best correlations between macroseismic and instrumental observations are obtained when the macroseismic parameters are averaged over 10-km-radius circles around each station. Macroseismic intensities predicted by published peak ground velocity (PGV)–intensity relationships agree with our observed intensities, contrary to those based on peak ground acceleration (PGA). Correlation between the macroseismic and instrumental data for intensities between II and V (EMS-98) is better for PGV than for PGA. Correlation with the response spectra exhibits clear frequency dependence for all macroseismic parameters. Horizontal and vertical components are significantly correlated with the macroseismic parameters between 1 and 10 Hz, a range corresponding to both natural frequencies of most buildings and high energy content in the seismic ground motion. Between 10 and 25 Hz, a clear lack of correlation between macroseismic and instrumental observations exists. It could be due to a combination of the decrease in the energy signal above 10 Hz, a high level of anthropogenic noise and an increase in variability in soil conditions. Above 25 Hz, the correlation coefficients between the acceleration response spectra and the macroseismic parameters are close to the PGA correlation level.     

The 7 and 11 May 1984 earthquakes in Abruzzo-Latium (Central Italy): reappraisal of the existing macroseismic datasets according to the EMS98

The aim of this paper is to provide a complete and reliable macroseismic knowledge of the events that stroke a large area in Central Italy on 7 and 11 May 1984. Previous studies, together with original accounts integrated with new and unpublished information, have been gathered and examined in order to re-evaluate macroseismic intensities in terms of the European Macroseismic Scale (EMS98). New intensity maps have been compiled; the total number of localities with available information for both the shocks increases from 1254 of the previous study to 1576. On the basis of the new dataset, the macroseismic magnitude of the first shock is M_W 5.6 which is lower than the previous macroseismic computation (M_W 5.7). Moreover, the topic of assessing macroseismic intensity in the presence of multiple shocks has been also investigated, proposing an unconventional approach to presenting the macroseismic data: an overall picture of the cumulative effects produced by all the seismic sequence is given to support a partial but faithful reconstruction of the second shock. This approach is inspired by the common experience in interpreting historical seismic sequences and gives a picture of the impact of the 1984 events on the territory.     

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 foundations of seismic microzoning

The Russian seismic safety regulations are based on macroseismic evaluations in terms of seismic intensity values. According to the seismic scale, intensity value determines both the degree of destruction of buildings and engineering structures and such parameters of seismic impacts as peak ground accelerations and response spectra. Application of macroseismic approach in various aspects of seismic microzoning has led to achievements in earthquake engineering; these aspects are discussed. The advantages and disadvantages of the macroseismic approach are considered.     

On a homogeneous representation of seismic history in large regions

The possibility to obtain a more complete and unbiased long-term history of seismic shakings over large territories than is explicitly reported from inhabited localities is discussed in the paper. An approach proposed for this purpose consists in complementing the spatial distribution of the macroseismic effect of earthquakes by calculated intensities at localities where information on felt shakings is absent. The calculated intensity is obtained on the basis of data on the epicentral intensity and location of epicenters provided by earthquake catalogs. This approach is applied to the analysis of the history of seismic shaking in Spain. The calculated intensities are shown to be comparable in accuracy with the ordinary practice of intensity determinations at national seismological centers.     

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.     

Application of the ESI-2007 Scale for Estimating the Intensity of the Kultuk Earthquake,August 27, 2008 (South Baikal)

The paper deals with issues related to the testing of the ESI-2007 scale by using as an example the real seismic event that occurred on August 27, 2008 in South Baikal. The main objective of the paper is to carry out a comparative assessment of the earthquake’s intensity based on traditional macroseismic scales and environmental seismic intensity (ESI-2007) scales. The results of the macroseismic survey served as the initial data. Analysis has been made on the ESI-2007 scale in conformity with the requirements for seismic scales. Particular emphasis has been placed on the type (or rank) of the ESI-2007 scale. Such an investigation is one of the first cases of the application of a new scale by the example of a regional seismic event.     

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.     

Transfrontier macroseismic observations of the Ml = 5.4 earthquake of February 22, 2003 at Rambervillers,France

The 2003 Ml = 5.4 Rambervillers earthquake, north-east of France, is the largest seismic event recorded north of the Alps since the 1992 Ms = 5.3, I₀ = VII, Roermond earthquake, Netherlands. With a maximum macroseismic intensity of VI-VII EMS-98, the 2003 event was broadly felt to a distance of 300 km from the epicentre. It provides a unique opportunity to test and compare the different procedures used in France, Germany and Switzerland when evaluating macroseismic intensities. The main purpose of this paper is to present a common transfrontier macroseismic map based on the EMS-98 intensity scale. Maximum horizontal accelerations recorded in the area are compared to the intensity values, and we propose to use a differential technique to re-estimate the magnitude of the 1682 Remiremont, I₀ = VIII, earthquake, which occurred 40 km south of Rambervillers.     

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 seismic risk assessment based on macroseismic intensity and spectrum approaches using ‘SeisVARA’

A comparative study of risk assessment methodologies based on macroseismic intensity and response spectrum approaches is presented. To facilitate the comparative study, a spreadsheet-based software tool ‘SeisVARA’ is developed for the estimation of earthquake risk, in which the seismic hazard can be specified either in terms of macroseismic intensity, or peak ground acceleration in combination with the spectral shape and soil amplification model of various earthquake building codes, or in terms of inelastic response spectra using the ‘next generation attenuation relationships’. A comparison of these different approaches is conducted for a typical city in northern India. In addition, the effect of different parameters, e.g., level of PGA, spectral shape, source-site parameters, and soil amplification models, is studied. It is observed that not only the different approaches result in widely varying damage and loss estimates, but also the variation of parameters within a given approach can result in considerable differences.     

An ordered probit model for seismic intensity data

Seismic intensity, measured through the Mercalli–Cancani–Sieberg (MCS) scale, provides an assessment of ground shaking level deduced from building damages, any natural environment changes and from any observed effects or feelings. Generally, moving away from the earthquake epicentre, the effects are lower but intensities may vary in space, as there could be areas that amplify or reduce the shaking depending on the earthquake source geometry, geological features and local factors. Currently, the Istituto Nazionale di Geofisica e Vulcanologia analyzes, for each seismic event, intensity data collected through the online macroseismic questionnaire available at the web-page www.haisentitoilterremoto.it. Questionnaire responses are aggregated at the municipality level and analyzed to obtain an intensity defined on an ordinal categorical scale. The main aim of this work is to model macroseismic attenuation and obtain an intensity prediction equation which describes the decay of macroseismic intensity as a function of the magnitude and distance from the hypocentre. To do this we employ an ordered probit model, assuming that the intensity response variable is related through the link probit function to some predictors. Differently from what it is commonly done in the macroseismic literature, this approach takes properly into account the qualitative and ordinal nature of the macroseismic intensity as defined on the MCS scale. Using Markov chain Monte Carlo methods, we estimate the posterior probability of the intensity at each site. Moreover, by comparing observed and estimated intensities we are able to detect anomalous areas in terms of residuals. This kind of information can be useful for a better assessment of seismic risk and for promoting effective policies to reduce major damages.     

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.     

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.     

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.     

中国大陆强震灾害范围的统计研究

强震的灾害范围是非常重要的地震学参数,它在地震灾害快速评估,地震保险及防震减灾等方面有着广泛的应用。本文全面,系统地收集了中国大陆宏观等震线资料,对没有地理坐标的等震线进行了数字化。精确计算了从１３０３年至１９９４年共１８３次强震的灾害范围,并对误差因素进行了分析。     

Field survey around strong motion stations and its implications on the seismic intensity in the Lushan earthquake on April 20, 2013

The M _s7.0 Lushan earthquake on April 20, 2013 is another destructive event in China since the M _s8.0 Wenchuan earthquake in 2008 and M _s7.1 Yushu earthquake in 2010. A large number of strong motion recordings were accumulated by the National Strong Motion Observation Network System of China. The maximum peak ground acceleration (PGA) at Station 51BXD in Baoxing Country is recorded as ?1,005.3 cm/s², which is even larger than the maximum one in the Wenchuan earthquake. A field survey around three typical strong motion stations confirms that the earthquake damage is consistent with the issued map of macroseismic intensity. For the oscillation period 0.3–1.0 s which is the common natural period range of the Chinese civil building, a comparison shows that the observed response spectrums are considerably smaller than the designed values in the Chinese code and this could be one of the reasons that the macroseismic intensity is lower than what we expected despite the high amplitude of PGAs. The Housner spectral intensities from 16 stations are also basically correlated with their macroseismic intensities, and the empirical distribution of spectral intensities from Lushan and Wenchuan Earthquakes under the Chinese scale is almost identical with those under the European scale.