Experimental vulnerability curves for the residential buildings of Iran

Iran is one of the most seismically active countries of the world located on the Alpine-Himalayan earthquake belt. More than 180,000 people were killed due to earthquakes in Iran during the last five decades. Considering the fact that most Iranians live in masonry and non-engineered houses, having a comprehensive program for decreasing the vulnerability of society holds considerable importance. For this reason, loss estimation should be done before an earthquake strikes to prepare proper information for designing and selection of emergency plans and the retrofitting strategies prior to occurrence of earthquake. The loss estimation process consists of two principal steps of hazard analysis and vulnerability assessment. After identifying the earthquake hazard, the first step is to evaluate the vulnerability of residential buildings and lifelines and also the social and economic impacts of the earthquake scenarios. Among these, residential buildings have specific importance, because their destruction will disturb the daily life and result in casualties. Consequently, the vulnerability assessment of the buildings in Iran is important to identify the weak points in the built environment structure. The aim of this research is to prepare vulnerability curves for the residential buildings of Iran to provide a proper base for estimating probable damage features by future earthquakes. The estimation may contribute fundamentally for better seismic performance of Iranian societies. After a brief review of the vulnerability assessment methods in Iran and other countries, through the use of the European Macroseismic method, a model for evaluating the vulnerability of the Iranian buildings is proposed. This method allows the vulnerability assessment for numerous sets of buildings by defining the vulnerability curves for each building type based on the damage observations of previous earthquakes. For defining the vulnerability curves, a building typology classification is presented in this article, which is representative of Iranian building characteristics. The hazard is described in terms of the macroseismic intensity and the EMS-98 damage grades have been considered for classifying the physical damage to the buildings. The calculated vulnerability indexes and vulnerability curves show that for engineered houses there is not any notable difference between the vulnerability of Iranian and Risk-UE building types. For the non-engineered houses, the vulnerability index of brick and steel structures is less than the corresponding values of the other unreinforced masonry buildings of Iran. The vulnerability index of unreinforced and masonry buildings of Iran are larger than the values of the similar types in Risk-UE and so the Iranian buildings are more vulnerable in this regard.     

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.     

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.     

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.     

The record of historic earthquakes in lake sediments of Central Switzerland

Deformation structures in lake sediments in Central Switzerland can be attributed to strong historic earthquakes. The type and spatial distribution of the deformation structures reflect the historically documented macroseismic intensities thus providing a useful calibration tool for paleoseismic investigations in prehistoric lake sediments.The Swiss historical earthquake catalogue shows four moderate to strong earthquakes with moment magnitudes of M_w=5.7 to M_w=6.9 and epicentral intensities of I₀=VII to I₀=IX that affected the area of Central Switzerland during the last 1000 years. These are the 1964 Alpnach, 1774 Altdorf, 1601 Unterwalden, and 1356 Basel earthquakes. In order to understand the effect of these earthquakes on lacustrine sediments, four lakes in Central Switzerland (Sarner See, Lungerer See, Baldegger See, and Seelisberg Seeli) were investigated using high-resolution seismic data and sediment cores. The sediments consist of organic- and carbonate-rich clayey to sandy silts that display fine bedding on the centimeter to millimeter scale. The sediments are dated by historic climate and environmental records, ¹³⁷Cs activity, and radiocarbon ages. Deformation structures occur within distinct zones and include large-scale slumps and rockfalls, as well as small-scale features like disturbed and contorted lamination and liquefaction structures. These deformations are attributed to three of the abovementioned earthquakes. The spatial distribution of deformation structures in the different lakes clearly reflects the historical macroseismic dataset: Lake sediments are only affected if they are situated within an area that underwent groundshaking not smaller than intensity VI to VII. We estimate earthquake size by relating the epicentral distance of the farthest liquefaction structure to earthquake magnitude. This relationship is in agreement with earthquake size estimations based on the historical dataset.     

Test of Source-Parameter Inversion of Intensity Data

We demonstrate that the approximate source kinematics of the San Fernando, 1971 earthquake can be back-predicted by analysing its macroseismic intensity data set (felt reports) objectively and quantitatively. This is done by inverting either the data set of the intensity values observed in all sites, or the intensities tessellated with the Voronoi polygons technique. It is shown that the kinematic characteristics found following our method (epicentral coordinates, source depth, seismic moment, rupture length, Mach number, fault plane solution) match those determined by other authors, via instrumental measurements, rather well. The prerequisite for obtaining these results is that local amplification must not affect groups of neighboring sites. It was possible to invert the U.S.G.S. ``felt reports' for the source because this data set is sufficiently uncontaminated by local site responses, and retains relevant regional traces of source effects. Isoseismal maps cannot be safely used for this task, because qualitative drawing criteria give subjective results. Isoseismals, based on incomplete space frequency samplings, give rise to spurious effects, whereas the Voronoi polygons produce easy-to-grasp, quantitative and objective, representations of macroseismic intensity data. The tests performed, up to now on a series of earthquakes, suggest that the combined use of tessellation and of our KF model is promising mostly for inverting intensities of preinstrumental earthquakes.     

Mapping seismic risk: the current crisis

The seismic risks to which populations are exposed should be estimated reliably for mitigation and preparation of response to disastrous earthquakes. Three parameters need to be known: Population numbers, properties of the built environment, and the seismic hazard. If we focus on large cities, we can say that at least one of these is known satisfactorily, namely the population, but not the other two. In the developing world, the numbers of buildings in a city are known only approximately, their distribution into building types (resistance to shaking) has to be assumed, and the distribution of types throughout the city is unknown. Recent verification of the world seismic hazard map has shown that it is grossly misleading: Instrumental measurements of accelerations due to six earthquakes were about three times larger, on average, than the maximum likely accelerations shown on the map; the macroseismic intensities reported for the last 60 earthquakes with M ≥ 7.5 were all significantly larger than expected, based on the hazard map (by 2.3 intensity units for the 12 deadliest earthquakes); and calculations of losses of life based on the hazard map underestimate the losses sustained in the 12 recent earthquakes with more than 1,000 fatalities by two to three orders of magnitude. This means that the seismic risk in most of the approximately 1,000 large cities at risk in the developing world is unknown. To remedy this intolerable situation, models for the built environment in cities need to be constructed, using cost-effective analyses of satellite images, and worst case scenario estimates of the losses in case of the nearest maximum credible earthquake.     

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.     

Problems in assessing macroseismic intensity from historical earthquake records

This paper tackles the problems connected with the use of macroseismic scales, with particular attention to MSK. Intensity estimates are in principle affected by uncertainties which are seldom quantified in practice. In the case of two localities (Scisciano and Orciano) damaged by two Italian earthquakes (1741, Fabriano and 1846, Orciano), which have a good historical documentation, the difficulties of obtaining the data necessary for MSK evaluation have been analysed. In order to assess the maximum range of uncertainty the data have been analysed in two extreme scenarios. Finally a comparison between intensity assigned using the MSK and MCS scales is analysed.     

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.     

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.     

Seismogeological effects and macroseismic fields

Several seismogeological effects of earthquakes on the environment have been collected by means of the ING macroseismic procedure. The fast and easy collection of information about ground and hydrological effects constitutes a complementary instrument in assessing the effects on human structures and the environment, and, ultimately, establishing intensity values. In this paper the distribution of seismogeological effects within intensity values has been verified, using data from three moderate Italian earthquakes; the density of the same effects are then discussed in comparison with isoseismal maps.     

Hazard mapping based on macroseismic data considering the influence of geological conditions

The object of this study is to consider directly the influence of regional geological conditions on the assessment of seismic hazard. It is assumed that macroseismic data at individual locations contain, in an average way, the influence of geological conditions.A Data Base referring to 199 historical (5) and instrumental (194, in the 1947–1993 period) events with macroseismic information in 1195 locations of Portugal was built. For any given seismic event, whenever macroseismic information was available at a location (town, village, etc.), an EMS-92 intensity value was estimated. To each one of those locations a geological unit, representing the most common type of soil, was assigned, based on the Geological Portuguese Map at a scale 1:500 000; the geological units were grouped into three categories: soft, intermediate and hard soils.The Data Base was used to determine the attenuation laws in terms of macroseismic intensity for the three different geological site conditions, using multiple linear regression analysis. The reasonability of the laws was tested by (i) checking residual distributions and (ii) comparing the map of isoseismals of important earthquakes with the isoseismals generated by the attenuation curves derived for each one of the three different soil classes, taking into consideration the soil class of each site. The main results of attenuation modeling are: high dispersion on macroseismic intensity data; all the models predict intensity values, for short hypocentral distances, lower than the ones observed; and for some important analyzed earthquakes and for the observed range of distances, the models confirm the expectancy that macroseismic intensity increases from hard to soft soil.The approach to obtain the hazard assessment at each location consisted in the use of the attenuation law specifically derived for the class of soil of that particular location. This method, which considers the influence of the regional geology, was illustrated with the mapping of hazard for the country for several return periods. Comparison with previous maps not taking into consideration the regional geological conditions emphasizes the importance of this new parameter. It can be concluded that (i) soil segmentation is clearly the cause for hazard increase in the region to the north of Lisbon, especially at sites with soft and intermediate soils as the ones in lower Tagus valley; the maximum increase on hazard is, in any case, less than one degree; (ii) when geological conditions are disregarded in the attenuation regression analysis, hazard pattern is similar to the one obtained for the case of hard soil everywhere.     

Intensity distribution of M 4.9 Haryana–Delhi border earthquake

In this study, we have prepared an intensity map based on macroseismic survey and all the available information from print and electronic media of damage and other effects due to March 05, 2012, M 4.9 Bahadurgarh (Haryana–Delhi border) earthquake and interpreted them to obtain modified Mercalli intensities (MMI) at over 62 locations surrounding the Haryana and Delhi. We have cross-checked the damage information from print and electronic media in the field at 25 sites within 110 km surrounding the epicenter for validation. Based on the questionnaire which is used in macroseismic survey and personal judgment, intensities were assigned accordingly as per physical survey at 25 sites and for rest based on media reporting. A maximum intensity of VI was assigned to this seismic event. Isoseismals of V and VI have been fully covered in the field observations. Beside this, some of the points have also been covered for isoseismal IV and isoseismal III and rest are based on media report only. The intensity map reveals several interesting features. Elliptically elongated shape of intensity map shows that most of the slightly damaged areas are concentrated toward the northwestern side of the epicenter having intensity V which may be due to directivity or site effects. A regression relation has also been derived between intensity and epicentral distance. The derived attenuation relation will be useful for assessing damage of a potential future earthquake (earthquake scenario–based planning purposes) for the Delhi NCR region.     

Verification of macroseismic methods on five ML > 5 instrumental earthquakes in France

Macroseismic data available for five of the most recent M_L > 5 earthquakes that occurred in the Pyrenees and in the Alps, were analyzed using the Sponheuer and the Levret relationship to estimate depth and magnitude respectively. The aim of this paper is to verify if simple and robust macroseismic methods used on recent instrumental earthquakes may provide a good tool to calibrate historical events in France. The excellent agreement found between macroseismic and instrumental estimates shows that macroseismic data of historical events may provide the means to lengthen the instrumental catalogue and better constrain the recurrence rates of earthquakes in moderate seismic rate regions.     

Historical intensity VIII earthquakes along the Rhone valley (Valais, Switzerland): primary and secondary effects

In recent years the upper Rhone Valley has been one of the most intensively investigated regions by the Swiss Seismological Service. The high seismicity in the region encourages research in the seismological field and one main focus has been historical seismology. This report presents the state of the art of our historical investigations by giving an overview of the effects of four damaging earthquakes with intensity larger than VII, for which a fairly large number of documents could be found and analyzed. The overview includes the events of 1584 (Aigle, epicentral intensity VIII), 1755 (Brig, epicentral intensity VIII), 1855 (Visp, epicentral intensity VIII), and 1946 (Sierre, epicentral intensity VIII for the main shock and intensity VII for the largest aftershock). The paper focuses mainly on primary and secondary effects in the epicentral region, providing the key data and a general characterization of the event. Generally, primary effects such as the reaction of the population and impact on buildings took more focus in the past. Thus building damage is more frequently described in historic documents. However, we also found a number of sources describing secondary effects such as landslides, snow avalanches, and liquefaction. Since the sources may be useful, we include citations of these documents. The 1584 Aigle event, for example, produced exceptional movements in Lake Geneva, which can be explained by an expanded sub aquatic slide with resultant tsunami and seiche. The strongest of the aftershocks of the 1584 event triggered a destructive landslide covering the villages Corbeyrier and Yvorne, VD. All macroseismic data on the discussed events are accessible through the webpage of the Swiss Seismological Service (http://www.seismo.ethz.ch).     

Seismicity and seismotectonic implications in the southern Baltic Sea area

An updated list of earthquakes and earthquake parameters (location, homogenized magnitude, macroseismic data) for the southern Baltic Sea area reveals activity north of the sea, whereas there are very few epicentres in the sea itself and in the region south of it. This is the first study to combine seismological data for the whole region to cover also the sea. Macroseismic data for the 1930 earthquake were reinvestigated leading to an intensity of V–VI (MM or MSK scale), a radius of perceptibility of 135 km and an unusually big focal depth of about 40 km. It is difficult to correlate individual earthquakes with specific faults, but some seismotectonic relations are suggested, e.g. for the Tornquist zone, the predominant structure of the region. Only few reliable focal-mechanism solutions exist. Possible seismogenic processes (ridge push, isostasy, etc) are discussed.     

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