Web-based macroseismic survey in Italy: method validation and results

A new method of macroseismic surveys, based on voluntary collaboration through the Internet, has been running at the Istituto Nazionale di Geofisica e Vulcanologia (INGV) since July 2007. The macroseismic questionnaire is addressed to a single non-specialist; reported effects are statistically analysed to extrapolate a probabilistic estimate of Mercalli Cancani Sieberg and European Macroseismic Scale intensities for that observer. Maps of macroseismic intensity are displayed online in almost real time and are continuously updated when new data are made available. For densely inhabited zones, we have received reports of felt effects for even very small events (M = 2). Six earthquakes are presented here, showing the ability of the method to give fast and interesting results. The effects reported in questionnaires coming from three towns are carefully analysed and assigned intensities are compared with those derived from traditional macroseismic surveys, showing the reliability of our web-based method.     

On an anisotropic attenuation law of the macroseismic intensity

An anisotropic attenuation law of the macroseismic intensity has been deduced, congruent with an anisotropic modelling of macroseismic fields. The results, obtained by analysing a set of earthquakes in Eastern Sicily and Southern Calabria, show a greater adaptability to the observed data as compared with those deduced using isotropic attenuation laws modified to take the anisotropy into account.     

Atlas of macroseismic maps for French earthquakes with their principal characteristics

The SIRENE macroseismic database has been utilized to draw isoseismal maps for the 140 best-documented French earthquakes, characterized by epicentral intensities of at least V (MSK) and located in all parts of the country. A study of focal depths derived from available local intensity data using an intensity versus distance decay law (Sponheuer) shows that the focal depths of most of the events considered do not exceed about 10 km. Their distribution correlates fairly well with regional dynamic geology features. A relationship is then computed between magnitude, intensity and focal distance, based on 73 instrumenta]ly recorded earthquakes (M _L between 3.3 and 6.3) and on 217 mean radius values (from 2 to 380 km) for isoseismals of intensity VIII to III (MSK). This relationship is applied to historical earthquakes contained in the database SIRENE which are characterised by their intensity only. These results are used in the evaluation as well deterministic as probabilistic of the seismic hazard on the national territory.     

Distribution of effects in the urban area of Rome During the October 14, 1997 Umbria Marche earthquake

In the last years many studies were carried out in order to assess the correlation between nearsurface geology and seismic response of downtown Rome during local and regional earthquakes. Nevertheless, only historical macroseismic data and theoretical modelling could be used so far to estimate the ground motion in the historical part of the city. The occurrence of the October 14, 1997 Umbria Marche event (M_w = 5.7), located about 115 km far from Rome but largely felt in the city, allowed us to verify our knowledge of the seismic response on the basis of the present urban setting of Rome.A prompt macroseismic survey in the urban area within the ring-road (G.R.A.) was performed. 949 questionnaires were compiled by means of 1842 interviews: such a large amount of macroseismic information, for the first time available in Rome, was related to 669 observation points providing a picture of local intensity variations. The areal distribution of collected intensity points is quite homogeneous.The data analysis points out a preferential distribution of larger effects (65%) in correspondence of Holocene alluvium outcrops. Furthermore, a tight correspondence between remarkable intensities and alluvial sediments is pointed out in suburban areas characterized by the minor hydrographic network of the Tiber River. Macroseismic anomalies were outlined in sectors of the urban area where local hydrogeological conditions can be responsible for a significant deterioration of geomechanical properties.     

Relation between macroseismic intensity and instrumental parameters of strong motions — A statistical approach

One hundred and twenty-eight strong ground motion CALTECH (Earthquakes in the U.S.A., 1940–1971) records of five Californian earthquakes, recorded at ground level and in basements, for which the modified Mercalli macroseismic intensities are known, were statistically processed to calculate 165 wave parameters. Correlation of the peak values of particle acceleration, velocity and displacement, energy, impulse and root-mean-square amplitude of the vibrations, durations of vibrations for certain levels of the peak amplitude, and spectral energies of 10 frequency windows with macroseismic intensity are discussed from the viewpoint of classification of ground-motion records.Paper presented at the 21st General Assembly of the European Seismological Commission held in Sofia 1988.     

Macroseismic intensity attenuation for crustal sources on romanian territory and adjacent areas

The macroseismic intensity attenuation (m.i.a.) laws for the main seismic provinces with crustal foci of the Romanian territory and adjacent areas were established. The input data consist of MSK-64 intensities of 18 earthquakes with epicentral/maximum intensity in the range V to X degrees (MSK scale), whose isoseismal maps were available.The attenuation was analyzed as a function of distance and azimuth and from the three main attenuation formulas (logarithmic, exponential, and power-law) (the last one was preferred, as it best fits the observed data) and, consequently, were used for each seismogenic region. The bulk of derived equations is intended to be further used in the assessment of the seismic hazards of Romania.     

Atténuation de l'intensité macrosismique pour la France métropolitaine : importance de l'intensité épicentrale

The main result of this work is to show that macroseismic intensity decay with distance strongly depends on the epicentral intensity. An attenuation law that takes this parameter into account is proposed for Metropolitan France, from the analysis of SISFRANCE macroseismic database. Such a model significantly reduces the difference between observed and theoretical intensities. A map of the attenuation variations is also set up for Metropolitan France. No major site effects are observed, but at a broad scale, young Alpine regions display a stronger attenuation than old Hercynian regions. To cite this article: P. Arroucau et al., C. R. Geoscience 338 (2006).     

Macroseismic Surveys in Theory and Practice

Macroseismology is the part of seismology that collects and evaluates non-instrumental data on earthquakes, i.e., effects on people, objects, buildings and nature. The methods that seismologists use for collecting and evaluating the macroseismic data are often based on long (trial-and-error) experience more than on some formal procedure. Until very recently manuals or guidelines on how to do a macroseismic survey were rare and often superficial. After an earthquake is felt in some region, the data are usually collected by means of questionnaires. Field survey is an obligatory tool that complements the questionnaires in the case of a damaging earthquake. An overview of the approaches to deriving the earthquake parameters (epicentre and barycentre, epicentral intensity, magnitude, depth, source parametres) from macroseismic data, as well as a review of some existing practices is given.     

Croation macroseismic database

Compilation of the digital Croatian Macroseismic Database is an ongoing project that started in 1995. Final version of the database will contain all the macroseismic data from the archive of the Department of Geophysics (including data for neighboring countries) without time or intensity threshold. The database has two subsets: intensity point database and the digital archive. A total of 426 earthquakes felt in Croatia are currently being analyzed. The job is carried out through several stages, each of them corresponding to digitalization of a different input data format (intensity maps, intensity point lists, questionnaires, etc.). 69 intensity maps of the earthquakes felt in Croatia, published in the Catalogue of Earthquakes in the Balkan Region were digitized first. The coordinates of intensity points were assigned to geographic locations (villages, cities, …) by consulting the database of latitudes and longitudes of all settlements in Croatia (database of settlements). 330 existing intensity lists (intensity and name of the place without coordinates) were entered into the database during the second step. The intensity maps which were not published so far will be scanned and digitized during the third stage of the project. The last stage will be input of primary data to the digital archive and reassessing of intensities according to the EMS macroseismic scale. The main goal of this paper is to describe the way of database compilation. Preliminary felt intensity map and number of intensity points map illustrate possibilities of macroseismic field analyze.     

A distribution-free analysis of magnitude-intensity relationships: an application to the Mediterranean region

A statistical analysis of the relations between macroseismic intensity and magnitude is presented. The examined data set contains earthquakes characterized by epicentral or maximum intensity ≥ VI which occurred in the Mediterranean region. As a first step, an empirical magnitude-intensity relationship has been determined by using the whole data set. Then, differences between experimental magnitude values and the ones expected on the basis of the empirical relationship have been correlated with some features related both to physical and data sources characteristics. On this basis, a distribution-free statistical approach has been developed to attempt a regionalization of the examined area, able to locally optimize the performances of magnitude-intensity relations. However, the results showed that data relative to larger events (intensity ≥ VII) are not sufficient to perform any reliable zonation of the area. Thus, well-constrained relationships determined for the whole Mediterranean region should be preferred to ill-defined local ones. Concerning smaller earthquakes (intensity VI), the analysis suggests that an efficient zonation could only be obtained if medium-scale variations (lower than 200 Km) are taken into account.     

Probabilistic macroseismic hazard assessment for Nicaragua (preliminary results)

A probabilistic macroseismic hazard assessment has been done for Nicaragua. For this, the most complete catalogue for Central America, compiled by NORSAR in Norway has been used. In this catalogue, empirical intensity attenuation relations were found. Using these empirical relations, magnitudes were changed to epicentral intensities expected in sites where no intensities had been reported. The calculated intensities from a polygon surrounding Nicaragua were used to assess the macroseismic hazard in the region. For the whole polygon, the cumulative intensity frequency was calculated resulting in a b-value of 0.60 for an intensity interval of V–IX. The time completeness was also studied indicating that, for strong events causing higher intensities (I ₀ VII), the catalogue is complete for events that have been recorded since 1840. The whole polygon was cut into independent seismotectonic regions where the statistical procedure (intensity frequency and time completeness) was done.     

Estimation of seismic hazard with allowance for noncircular isoseismals of arbitrary orientation

A new mathematical model describing the field of macroseismic intensity has been elaborated. It is based on elliptic isoseismals. The orientation of the main axes of elliptic isoseismals depends on the direction of stretching of the main geological structures on the investigated territory.The new model of a macroseismic field was applied to the territory of Eastern Uzbekistan. Some results of macroseismic investigations of the effect of large regional earthquakes were used as initial data.A noncircular model of a macroseismic field was introduced into the integral of the seismic shakability of Riznichenko and, according to the model, a macroseismic shakability map for the territory of Eastern Uzbekistan was computed in isolines of the long-term mean return period of vibrations for the intensity I 8.Paper presented at the 21st General Assembly of the European Seismological Commission, Sofia, 1988.     

Unification of Macroseismic Data Collection Procedures: A Pilot Project for Border Earthquakes Assessment

The border area between Italy and Slovenia has been chosen for a proposed project on border earthquake's assessment (BEA). Attempts at the unification of macroseismic practice among the countries have been rare and not very successful, due to the different and often incompatible local traditions for the derivation of intensity data. With the publication of the new EMS-92 scale, the idea was born to try to unify the procedures in data collection on both sides of the border. The pilot project was proposed, consisting mainly of the preparation for eventual strong earthquakes in the future. A special civil engineering form for damage assessment was designed to be used in the damaged area. The path of this idea is presented.     

The determination of the epicentre by a vectorial modelling of macroseismic intensity distribution

A vectorial modelling of observed macroseismic intensity aimed at the analytical determination of the epicentre is proposed here. The methodology is based on the determination of a plane system of vectors which characterises the macroseismic intensity distribution. The epicentre of each seismic event considered is determined as the centre of this vector system by an analytical expression which is independent from all possible directions of seismic energy propagation. The analysis of the intensity distribution is carried out by a new model called a macroseismic plane, different from the one known as macroseismic field, formed by a set of small areas built around the observed intensity points; hence its name.With the proposed methodology, some earthquakes in southern Italy and eastern Sicily are analysed calculating their epicentres, also for distributions of observed intensity which are particularly complex.     

Site response of the Ganges basin inferred from re-evaluated macroseismic observations from the 1897 Shillong, 1905 Kangra,and 1934 Nepal earthquakes

We analyze previously published geodetic data and intensity values for the M _s = 8.1 Shillong (1897), M _s = 7.8 Kangra (1905), and M _s = 8.2 Nepal/Bihar (1934) earthquakes to investigate the rupture zones of these earthquakes as well as the amplification of ground motions throughout the Punjab, Ganges and Brahmaputra valleys. For each earthquake we subtract the observed MSK intensities from a synthetic intensity derived from an inferred planar rupture model of the earthquake, combined with an attenuation function derived from instrumentally recorded earthquakes. The resulting residuals are contoured to identify regions of anomalous intensity caused primarily by local site effects. Observations indicative of liquefaction are treated separately from other indications of shaking severity lest they inflate inferred residual shaking estimates. Despite this precaution we find that intensites are 1–3 units higher near the major rivers, as well as at the edges of the Ganges basin. We find evidence for a post-critical Moho reflection from the 1897 and 1905 earthquakes that raises intensities 1–2 units at distances of the order of 150 km from the rupture zone, and we find that the 1905 earthquake triggered a substantial subsequent earthquake at Dehra Dun, at a distance of approximately 150 km. Four or more M = 8 earthquakes are apparently overdue in the region based on seismic moment summation in the past 500 years. Results from the current study permit anticipated intensities in these future earthquakes to be refined to incorporate site effects derived from dense macroseismic data.     

Detecting Voronoi (area-of-influence) polygons

Voronoi, or area-of-influence, polygons are convex, space-filling polygons constructed around a set of points (Voronoi centers) such that each polygon contains all points closer to its Voronoi center than to the center of any other polygon. The relationship of Voronoi centers to edges of Voronoi polygons is used to test whether any convex tessellation consists of Voronoi polygons. This test amounts to finding Voronoi centers that best fit the given tessellation. Voronoi centers are found by solving two systems of linear equations. These equations represent (1) conditions on the slope of polygon edges relative to the slope of lines through Voronoi centers, and (2) conditions on the distance from edges to Voronoi centers. Least squares and constrained least-squares solutions are used to solve the two systems. Different methods of solution can provide insight as to how a tessellation varies from Voronoi polygons. A goodness-of-fit statistic is derived and examined by testing randomly generated convex tessellations. Some polygonal ice cracks provide an example of naturally occurring polygons that are approximated closely by Voronoi polygons.     

On macroseismic attenuation parameters in the comparison of theoretical seismic hazard maps for Southern Calabria and Eastern Sicily (Southern Italy)

For two different zonations of the Calabro-Sicilian region the macroseismic intensity attenuation parameters are computed using the Grandori model.Some modifications to this relation are proposed here in order to allow its applicability also to those cases in which the data set available for each source zone does not present three successive epicentral intensity levels (I ₀) as required for its application.Maps of theoretical distribution of maximum expected intensities were plotted for both zonations adopted, in order to test the reliability, also using other attenuation models.The results of the analysis show how the proposed modifications to the Grandori relation allow the macroseismic intensity attenuation modelling even in cases in which only one intensity level is available.Finally, the comparison of theoretical maximum expected intensity distribution maps, computed for different attenuation models and seismogenic zonations, shows that the maps plotted using the Grandori model and the zonation adopted by G.N.D.T. (Gruppo Nazionale Difesa dai Terremoti - C.N.R. Italy) are more reliable.     

Systematization of active faults for the assessment of the seismic hazard

A systematization of active faults has been developed based on the progress of scientists from the leading countries in the world in the study of seismotectonics and seismic hazard problems. It is underlain by the concept of the fault-block structure of the geological-geophysical environment governed by the interaction of differently oriented active faults, which are divided into two groups—seismogenic and nonseismogenic faults. In seismogenic fault zones, the tectonic stress accumulated is relieved by means of strong earthquakes. Nonseismogenic fault zones are characterized by creep displacement or short-term, oscillatory, and reciprocal movements, which are referred to local superintense deformations of the Earth’s crust (according to the terminology used by Yu.O. Kuz’min). For a situation when a strong earthquake happens, a subgroup of seismodistributing faults has been identified that surround the seismic source and affect the distribution of the seismic waves and, as a consequence, the pattern of the propagation of the coseismic deformations in the fault-block environment. Seismodistributing faults are divided into transit and sealing faults. Along transit faults, secondary coseismic effects (landfalls, landslides, ground fractures, liquefaction, etc) are intensified during earthquakes. In the case of sealing faults, enhancement of the coseismic effects can be observed on the disjunctive limb nearest to the epicenter, whereas, on the opposite limb, the intensity of such effects appreciably decreases. Seismogenic faults or their systems are associated with zones of earthquake source origination (ESO), which include concentrated seismicity regions. In such zones, each earthquake source is related to the evolution of a fault system. ESO zones also contain individual seismogenic sources being focuses of strong earthquakes with M of ≥5.5 in the form of ruptures, which can be graphically represented in 2D or 3D as a surface projection of the source. Depending on the type of data based on which they are identified, individual seismogenic sources are divided into geological-geophysical and macroseismic sources. The systematization presented is the theoretical basis for and the concept of the relational database that is being developed by the authors as an information system for the generation of seismotectonic GIS projects required for the subsequent analysis of the seismic hazard and the assessment of the probability of the origination of macroseismic earthquake effects in a predetermined location.     

Recognition of earthquake-related damage in archaeological sites: Examples from the Dead Sea fault zone

Archaeological structures that exhibit seismogenic damage expand our knowledge of temporal and spatial distribution of earthquakes, afford independent examination of historical accounts, provide information on local earthquake intensities and enable the delineation of macroseismic zones. They also illustrate what might happen in future earthquakes. In order to recover this information, we should be able to distinguish earthquake damage from anthropogenic damage and from other natural processes of wear and tear. The present paper reviews several types of damage that can be attributed with high certainty to earthquakes and discusses associated caveats. In the rare cases, where faults intersect with archaeological sites, offset structures enable precise determination of sense and size of slip, and constrain its time. Among the characteristic off-fault damage types, I consider horizontal shifting of large building blocks, downward sliding of one or several blocks from masonry arches, collapse of heavy, stably-built walls, chipping of corners of building blocks, and aligned falling of walls and columns. Other damage features are less conclusive and require additional evidence, e.g., fractures that cut across several structures, leaning walls and columns, warps and bulges in walls. Circumstantial evidence for catastrophic earthquake-related destruction includes contemporaneous damage in many sites in the same area, absence of weapons or other anthropogenic damage, stratigraphic data on collapse of walls and ceilings onto floors and other living horizons and burial of valuable artifacts, as well as associated geological palaeoseismic phenomena such as liquefaction, land- and rock-slides, and fault ruptures. Additional support may be found in reliable historical accounts. Special care must be taken in order to avoid circular reasoning by maintaining the independence of data acquisition methods.