Deterministic earthquake damage and loss assessment for the city of Bucharest, Romania

On March 4, 1977, an earthquake with a moment magnitude M _w 7.4 at a hypocentral depth of 94 km hit the Vrancea region (Romania). In Bucharest alone, the earthquake caused severe damage to 33,000 buildings while 1,424 people were killed. Under the umbrella of the SAFER project, the city of Bucharest, being one of the larger European cities at risk, was chosen as a test bed for the estimation of damage and connected losses in case of a future large magnitude earthquake in the Vrancea area. For the conduct of these purely deterministic damage and loss computations, the open-source software SELENA is applied. In order to represent a large event in the Vrancea region, a set of deterministic scenarios were defined by combining ranges of focal parameters, i.e., magnitude, focal depth, and epicentral location. Ground motion values are computed by consideration of different ground motion prediction equations that are believed to represent earthquake attenuation effects in the region. Variations in damage and loss estimates are investigated through considering different sets of building vulnerability curves (provided by HAZUS-MH and various European authors) to characterize the damaging behavior of prevalent building typologies in the city of Bucharest.     

Microzonation of Bucharest: State-of-the-Art

— The 1940 (M_w=7.7) and 1977 (M_w=7.4) Vrancea earthquakes (Romania) inflicted heavy damage and casualties in Bucharest and the statistics indicate a recurrence interval of 25 years for M_w 7.0 events. Under these circumstances, the seismic microzonation represents important information for detailed urban planning that establishes an appropriate level of preparedness to the earthquake threat. This paper reviews the main studies concerning the seismicity of the Vrancea region, the site conditions of the city, the characterization of the building stock, and the codes of practice that regulate the antiseismic design. The first-order microzonation of Bucharest was performed starting from the existing database of structural and geotechnical parameters. New insights originating from direct instrumental observation and interpretation of the local effects as well as realistic numerical modeling that update and improve the input data necessary for a detailed microzoning map of the city are also discussed.     

Accounting for site effect in probabilistic assessment of seismic hazard for Romania and Bucharest: a case of deep seismicity in Vrancea zone

The paper presents recent achievements in evaluations of site-dependent seismic hazard in Romania and the capital city of Bucharest caused by the Vrancea focal zone (SE-Carpathians). The zone is characterized by a high rate of occurrence of large earthquakes in a narrow focal volume at depths 60–170 km. The database that was used for the hazard evaluation includes parameters of seismicity, ground-motion source scaling and attenuation models (Fourier amplitude spectra), and site-dependent spectral amplification functions. Ground-motion characteristics were evaluated on the basis of several hundred records from more than 120 small magnitude (M 3.5–5) earthquakes occurred in 1996–2001 and a few tens of acceleration records obtained during four large (M 7.4, 7.2, 6.9 and 6.3) earthquakes. The data provide a basis for probabilistic seismic hazard assessment in terms of peak ground acceleration, peak spectral acceleration and MSK intensity using Fourier amplitude spectra for various exceedance probabilities or average return periods. It has been shown that the influence of geological factors plays very important role in distribution of earthquake ground-motion parameters along the territory of Romania.     

A reconnaissance report for the Romanian earthquake of 4 March 1977

The engineering aspects of the 4 March 1977 Romanian earthquake are presented. They are based upon a field investigation conducted by the writers in Bucharest and in southern Romania in collaboration with members of the Building Research Institute of Romania, during the period 25–31 March 1977. This report covers general observations, data and evaluation on the character of the earthquake, structural damage inflicted by it, performance of different types of buildings during the earthquake and relief operations.     

Comparative performance of friction damped systems and base isolation systems for earthquake retrofit and aseismic design

The seismic performance of a steel framed structure equipped with (i) friction damping devices and (ii) base isolators is compared. A parametric study based on energy concepts is performed first using time-history dynamic analysis to determine the optimum properties of the two systems when excited by an earthquake whose energy is distributed over a relatively broad frequency band (1940 El Centro, N-S). Using these same properties, the responses of the two structural systems are then examined when excited by earthquakes whose power content essentially is concentrated at the low frequency end of the energy spectrum (1977 Romania, Bucharest, N-S and 1985 Mexico, SCT E-W). The results of the study show that, while both systems similarly reduce the response of conventional structures to the California earthquake, the friction damped structure exhibits a superior performance under the low frequency earthquakes. Very large shear forces and displacements are observed when the Romania and Mexico earthquakes are applied to the base isolated structure, indicating that the performance of a base isolated structure depends on the characteristics of the site earthquake. By comparison, friction damped structures are shown to behave favourably under the three earthquakes studied; this suggests that friction damping devices offer a more consistent way of protecting structures during earthquakes.     

Seismic vulnerability assessment and characterisation of the buildings on Faial Island, Azores

The earthquake of the 9th of July 1998 that hit in the central group of the Azores archipelago greatly affected the islands of Faial, Pico and S?o Jorge, reaching a magnitude of Mw 6.2 with the epicentre located about 15km northeast of the Faial Island. This earthquake allowed the collection of an unprecedented quantity of data concerning the characterisation of the building stock and the damage suffered by construction. This is the main purpose of this research, consisting essentially of three main aspects: (i) A detailed characterisation of the building stock, assigning a five category classification, from old traditional rubble stone masonry to reinforced concrete moment framed buildings; (ii) A detailed damage grade classification based on the different damage mechanisms observed; and, (iii) A seismic vulnerability assessment of the building stock. The results of the vulnerability assessment together with the building stock database and damage classification were integrated into a GIS tool, allowing the spatial visualation of damage scenarios, which is potentially useful for the planning of emergency response strategies and retrofitting priorities to mitigate and manage seismic risk.     

Vrancea Source Influence on Local Seismic Response in Bucharest

—?The mapping of the seismic ground motion in Bucharest, due to the strong Vrancea earthquakes, is carried out using a complex hybrid waveform modeling method that allows easy parametric tests. Starting from the actually available strong motion database, we can make realistic predictions for the possible ground motion. The basic information necessary for the modeling consists of: (a) The representative mechanisms for the strong subcrustal events, (b) the average regional structural model, and (c) the local structure for Bucharest. Two scenario earthquakes are considered and the source influence on the local response is analyzed in order to define generally valid ground motion parameters, to be used in the seismic hazard estimations. The source has its own (detectable) contribution on the ground motion and its effects on the local response in Bucharest are quite stable on the transversal component (T), while the radial (R) and vertical (V) components are sensitive to the scenario earthquake. Although the strongest local effects affect the T component, both observed and synthetic, a complete determination of the seismic input for the built environment requires the knowledge of all three components of motion (R, V, T). The damage observed in Bucharest for the March 4, 1977 Vrancea event, the strongest earthquake to strike the city in modern times, is in agreement with the synthetic signals and local response.     

Evaluation of the seismic hazard for 20 cities in Romania using Monte Carlo based simulations

This work focuses on the evaluation of the seismic hazard for Romania using earthquake catalogues generated by a Monte Carlo approach. The seismicity of Romania can be attributed to the Vrancea intermediate-depth seismic source and to 13 other crustal seismic sources. The recurrence times of large magnitude seismic events (both crustal and subcrustal), as well as the moment release rates are computed using simulated earthquake catalogues. The results show that the largest contribution to the overall moment release for the crustal seismic sources is from the seismic regions in Bulgaria, while the seismic regions in Romania contribute less than 5% of the overall moment release. In addition, the computations show that the moment release rate for the Vrancea subcrustal seismic source is about ten times larger than that of all the crustal seismic sources. Finally, the Monte Carlo approach is used to evaluate the seismic hazard for 20 cities in Romania with populations larger than 100,000 inhabitants. The results show some differences between the seismic hazard values obtained through Monte-Carlo simulation and those in the Romanian seismic design code P100-1/2013, notably for cities situated in the western part of Romania that are influenced by local crustal seismic sources.     

Advanced real-time acquisition of the Vrancea earthquake early warning system

The Vrancea seismogenic zone in Romania represents a peculiar source of seismic hazard, which is a major concern in Europe, especially to neighboring regions of Bulgaria, Serbia and Republic of Moldavia. Earthquakes in the Carpathian–Pannonian region are confined to the crust, except the Vrancea zone, where earthquakes with focal depth down to 200 km occur. One of the cities most affected by earthquakes in Europe is Bucharest. Situated at 140–170 km distance from Vrancea epicenter zone, Bucharest encountered many damages due to high energy Vrancea intermediate-depth earthquakes; the March 4, 1977 event (M_w=7.2) produced the collapse of 36 buildings with 8–12 levels, while more than 150 old buildings were seriously damaged. A dedicated set of applications and a method to rapidly estimate magnitude in 4–5 s from detection of P wave in the epicenter were developed. They were tested on all recorded data. The magnitude error for 77.9% of total considered events is in the interval [−0.3, +0.3] magnitude units. This is acceptable taking into account that the magnitude is computed from only 3 stations in a 5 s time interval (1 s delay is caused by data packing). The ability to rapidly estimate the earthquake magnitude combined with powerful real-time software, as parts of an early warning system, allows us to send earthquake warning to Bucharest in real time, in about 5 s after detection in the epicenter. This allows 20–27 s warning time to automatically issue preventive actions at the warned facility.     

Rapid Earthquake Information for Bucharest

Advances in sensor techniques, digital acquisition, communication and computing allowed establishment of the field of real-time seismology (Kanamori et al., 1997) and the design of earthquake information systems that provide rapid information on ground motion and potential damage in an area for which a ground motion estimation data base is available and where the inventory of buildings and infrastructure and its vulnerability to shaking is known. We discuss a system for the Romanian capital of Bucharest, where an early warning system with about 30 seconds lead time, shake maps and a sophisticated damage estimation tool for Bucharest have been developed in the frame of the Collaborative Research Center on Strong Earthquakes ().     

Using the damage from 2010 Haiti earthquake for calibrating vulnerability models of typical structures in Port-au-Prince (Haiti)

After the 2010 Haiti earthquake, that hits the city of Port-au-Prince, capital city of Haiti, a multidisciplinary working group of specialists (seismologist, geologists, engineers and architects) from different Spanish Universities and also from Haiti, joined effort under the SISMO-HAITI project (financed by the Universidad Politecnica de Madrid), with an objective: Evaluation of seismic hazard and risk in Haiti and its application to the seismic design, urban planning, emergency and resource management. In this paper, as a first step for a structural damage estimation of future earthquakes in the country, a calibration of damage functions has been carried out by means of a two-stage procedure. After compiling a database with observed damage in the city after the earthquake, the exposure model (building stock) has been classified and through an iteratively two-step calibration process, a specific set of damage functions for the country has been proposed. Additionally, Next Generation Attenuation Models (NGA) and \(\hbox {Vs}^{30}\) models have been analysed to choose the most appropriate for the seismic risk estimation in the city. Finally in a next paper, these functions will be used to estimate a seismic risk scenario for a future earthquake.     

Ground-motion prediction equations for the intermediate depth Vrancea (Romania) earthquakes

We present the regional ground-motion prediction equations for peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration (PSA), and seismic intensity (MSK scale) for the Vrancea intermediate depth earthquakes (SE-Carpathians) and territory of Romania. The prediction equations were constructed using the stochastic technique on the basis of the regional Fourier amplitude spectrum (FAS) source scaling and attenuation models and the generalised site amplification functions. Values of considered ground motion parameters are given as the functions of earthquake magnitude, depth and epicentral distance. The developed ground-motion models were tested and calibrated using the available data from the large Vrancea earthquakes. We suggest to use the presented equations for the rapid estimation of seismic effect after strong earthquakes (Shakemap generation) and seismic hazard assessment, both deterministic and probabilistic approaches.     

Detailed assessment of structural characteristics of Turkish RC building stock for loss assessment models

Assessment of the seismic vulnerability of the building stock in the earthquake-prone Marmara region of Turkey is of growing importance since such information is needed for reliable estimation of the losses that possible future earthquakes are likely to induce. The outcome of such loss assessment exercises can be used in planning of urban/regional-scale earthquake protection strategies; this is a priority in Turkey, particularly following the destructive earthquakes of 1999. Considering the size of the building inventory, Istanbul and its surrounding area is a case for which it is not easy to determine the structural properties and characteristics of the building stock. In this paper, geometrical, functional and material properties of the building stock in the northern Marmara Region, particularly around Istanbul, have been investigated and evaluated for use in loss estimation models and other types of statistic- or probability-based studies. In order to do that, the existing reinforced concrete (RC) stock has been classified as ‘compliant’ or ‘non-compliant’ buildings, dual (frame-wall) or frame structures and emergent or embedded-beam systems. In addition to the statistical parameters such as mean values, standard deviations, etc., probability density functions and their goodness-of-fit have also been investigated for all types of parameters. Functionalities such as purpose of use and floor area properties have been defined. Concrete properties of existing and recently constructed buildings and also characteristics of 220 and 420 MPa types of steel have been documented. Finally, the financial effects of retrofitting operations and damage repair have been investigated.     

Seismic response estimation method for earthquake‐damaged RC buildings

This study proposes a procedure for identifying spectral response curves for earthquake‐damaged areas in developing countries without seismic records. An earthquake‐damaged reinforced concrete building located in Padang, Indonesia was selected to illustrate the identification of the maximum seismic response during the 2009 West Sumatra earthquake. This paper summarizes the damage incurred by the building; the majority of the damage was observed in the third story in the span direction. The damage was quantitatively evaluated using the damage index R according to the Japanese guidelines for post‐earthquake damage evaluation. The damage index was also applied to the proposed spectral response identification method. The seismic performance of the building was evaluated by a nonlinear static analysis. The analytical results reproduced a drift concentration in the third story. The R‐index decreased with an increase in the story drift, which provided an estimation of the maximum response of the building during the earthquake. The estimation was verified via an earthquake response analysis of the building using ground acceleration data, which were simulated based on acceleration records of engineering bedrock that considered site amplification. The maximum response estimated by the R‐index was consistent with the maximum response obtained from the earthquake response analysis. Therefore, the proposed method enables the construction of spectral response curves by integrating the identification results for the maximum responses in a number of earthquake‐damaged buildings despite a lack of seismic records. Copyright © 2016 The Authors. Earthquake Engineering & Structural Dynamics published by John Wiley & Sons Ltd.     

Geological, geophysical and seismological criteria for local response evaluation in Bucharest urban area

Vrancea major intermediate-depth earthquakes produced extreme damage in Bucharest city, located at about 165 km epicenter distance. Our purpose is to investigate the influence of local geological conditions upon the seismic motion in Bucharest in case of large (M>7) Vrancea earthquakes. Two input data sets are used: (a) geological, geotechnical and geophysical information, including in situ measurements, and (b) acceleration recordings of Vrancea earthquakes. Local response evaluation based on first dataset is confirmed by the spectral analysis of the earthquake records. Two main features are outlined: non-stationarity of ground motion dynamic amplification from one event to other and inadequacy of limiting the investigation depth to uppermost 30 m to evaluate ground dynamic characteristics. Consequently (1) we cannot extrapolate the ground motion response determined for moderate and small earthquakes to anticipate the effects of the large Vrancea shocks and (2) the local response is controlled by the entire package of Quaternary deposits which are significantly deeper than 30 m depth beneath Bucharest Area.     

The role of the urban system dysfunction in the assessment of seismic risk in the Mt. Etna area (Italy)

A procedure for seismic risk assessment is applied to the Mt. Etna area (eastern Sicily, Italy) through assessment of urban system dysfunction following the occurrence of an earthquake. The tool used is based on the Disruption Index as a concept implemented in Simulator QuakeIST, which defines urban disruption following a natural disaster. The first element of the procedure is the definition of the seismic input, which is based on information about historical seismicity and seismogenic faults. The second element is computation of seismic impact on the building stock and infrastructure in the area considered. Information on urban-scale vulnerability was collected and a geographic information system was used to organise the data relating to buildings and network systems (e.g., building stock, schools, strategic structures, lifelines). The central idea underlying the definition of the Disruption Index is identification and evaluation of the impact on a target community through the physical elements that most contribute to severe disruption. The procedure applied in this study (i.e., software and data) constitutes a very useful operational tool to drive the development of strategies to minimise risks from earthquakes.     

A PRELIMINARY STUDY ON THE METHOD OF SEISMIC INTENSITY ASSESSMENT BASED ON RESIDENTIAL BUILDING DATA AND HIGH RESOLUTION REMOTE SENSING IMAGES

After destructive earthquakes, the assessment result of seismic intensity is an important decision-making basis for emergency rescue, recovery and reconstruction. This job requires higher timeliness by government and society. Because remote sensing technology is not affected by the terrible traffic conditions on the ground after the earthquake, large-scale seismic damage information in the earthquake area can be collected in a short time by the remote sensing image. The remote sensing technique plays a more and more important role in rapid acquisition of seismic damage information, emergency rescue decision-making, seismic intensity assessment and other work. On the basis of previous studies, this paper proposes a new method to assess seismic intensity by using remote sensing image, i.e. to interpret the building collapse rate of a residential quarter after an earthquake by high-resolution remote sensing images. If there already are detailed building data and building structure vulnerability matrix data of a residential area, we can calculate the building collapse rate under any intensity values in this residential area by using the theory of earthquake damage prediction. Assuming that the building collapse rate interpreted by remote sensing is equal to the building collapse rate predicted by using the existing data, it will be easy to calculate the actual seismic intensity of the residential area in this earthquake event. Based on this idea, according to the relevant standard specifications issued by China Earthquake Administration, this paper puts forward some functional models, such as the calculation model of building collapse rate based on remote sensing, the data matrix model of residential building structure, the prediction function matrix model of residential building collapse rate and the prediction model of residential building collapse rate. A formula for calculating seismic intensity by using remote sensing interpretation of collapse rate is also proposed. To test and verify the proposed method, this paper takes two neighboring blocks of Jiegu Town after the Yushu M7.1 earthquake in Qinghai Province as an example. The building structure matrix of the study block was constructed by using pre-earthquake 0.6m resolution satellite remote sensing image(QuickBird, acquired on November 6, 2004), post-earthquake 0.2m aerial remote sensing image(acquired by National Bureau of Surveying and Mapping, April 15, 2010) and some field investigation data. The building collapse rate in the two blocks was calculated by using the interpretation results of seismic damage from the Remote Sensing Technology Coordinating Group of China Seismological Bureau. The seismic damage matrix of building structures in Yushu area is constructed by using the abundant scientific data of the scientific investigation team of the project “Comprehensive Scientific Investigation of the Yushu M7.1 Earthquake in Qinghai Province” of China Seismological Bureau. On this basis, the collapse rate prediction function of different structures in Yushu area is constructed. According to the prediction function of collapse rate and the building structure matrix of the two blocks, the building collapse rate under different intensity values is predicted, and the curve of intensity-collapse rate function is drawn. By comparing the building collapse rate interpreted by remote sensing and the intensity-collapse rate function curve of this two blocks, the seismic intensity of both blocks are calculated to be the same value: Ⅸ degree, which is consistent with the results of the field scientific investigation of the earthquake. The validation shows that the method proposed in this paper can effectively avoid the influence caused by the difference of seismic performance of buildings and accurately evaluate seismic intensity when using remote sensing technique. The method has certain application value for earthquake emergency work.     

高频GNSS实时地震学与地震预警研究现状

为实现从注重灾后救助向注重灾前预防转变,如何提高地震灾害监测预警和风险防范能力成为我们关注的重点.本文给出了国际上GNSS位移记录、强震动加速度记录、测震速度记录在地震预警中的应用现状,并总结了各自的特点,归纳出围绕高频GNSS地震学在震级与破裂过程实时反演中的几个需要进一步研究的关键问题：（1）引入北斗系统,基于高频GNSS（GPS/BDS）双系统的实时位移解算方法来提高实时单站位移解算精度,使实时解算精度达到厘米级;（2）开展强震仪加速度记录基线偏移校正研究,弥补地震近场GNSS站密度不足问题;（3）强震仪加速度记录与GNSS位移记录特点不同,开展强震仪加速度数据与GNSS位移数据实时融合处理研究,快速获得包含丰富地震形变和速率的波形数据;（4）测震学方法可快速估算震级,但是在强震发生时会出现震级饱和现象,造成震级估算偏低.需要开展基于GNSS位移时间序列的多种方法相结合的实时震级估算方法研究,通过与地震学方法比较和结合,来得到精度高、计算快的震级估值算法;（5）基于高频GNSS、断层初始模型快速选取、断层尺度、参数自适应调整是快速判断断层破裂方向的基础,在断层破裂过程自适应准实时反演算法方面需要进一步加强.通过国内外研究现状调研、分析,表明基于高频GNSS地震学的震级快速确定、震源破裂过程准实时反演算法的发展将对我国地震预警系统从"二网融合"到"三网融合"提供坚实的技术支撑.     

Bi-normalized response spectra and seismic intensity in Bucharest for 1986 and 1990 Vrancea seismic events

The Vrancea subcrustal earthquakes of August 30,1986 and May 30,1990 are the two most recent seismic events that have occurred in Romania with moment magnitudes M W ≥ 6.9.The spectral analysis of the strong ground motions recorded in Bucharest reveals that despite small differences in magnitude between the 1986 and 1990 earthquakes,their frequency contents are very different,sometimes even opposing.The main focus of this study is to conduct a comparative analysis of the response spectra in terms of the bi-normalized response spectra(BNRS) proposed by Xu and Xie(2004 and 2007) for strong ground motions recorded in Bucharest during these two seismic events.The mean absolute acceleration and relative velocity response spectra for the two earthquakes are discussed and compared.Furthermore,the mean bi-normalized absolute acceleration and normalized relative velocity response spectra with respect to the control period T C are computed for the ground motions recorded in Bucharest in 1986 and 1990.The predominant period T P is also used in this study for the normalization of the spectral period axis.Subsequently,the methodology proposed by Martinez-Perreira and Bommer(1998) is applied in order to estimate the seismic intensity of the two events.The results are discussed and several conclusions regarding the possibility of using the bi-normalized response spectra(BNRS) are given.     

基于广义帕累托分布的地震震级分布尾部特征分析

极值理论在地震危险性分析中有着重要应用, 发震震级超过某一阈值的超出量分布可以近似为广义帕累托分布. 基于广义帕累托分布给出了若干地震活动性参数的估计公式, 包括强震震级分布、 地震复发周期和重现水平、 期望重现震级、 地震危险性概率和潜在震级上限等; 以云南地区震级资料为基础数据, 讨论了阈值选取、 模型拟合诊断和参数估计; 在此基础上计算了该地区的地震活动性参数. 结果表明, 广义帕累托分布较好地刻画了强震震级分布, 通过超阈值(POT)模型计算的复发周期与实际复发间隔统计基本一致, 高分位数估计在一定阈值范围内表现稳定, 为工程抗震中潜在震级上限的确定提供了一种途径.