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.     

Deterministic Approach for the Seismic Microzonation of 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 which combines the modal summation technique, valid for laterally homogeneous anelastic media, with finite-differences technique, and optimizes the advantages of both methods. For recent earthquakes, it is possible to validate the modeling by comparing the synthetic seismograms with the records. We consider for our computations the frequency range from 0.05 to 1.0 Hz and control the synthetic signals against the accelerograms of the Magurele station, low-pass filtered with a cut-off frequency of 1.0 Hz of the 3 last major strong (M_w > 6) Vrancea earthquakes. Using the hybrid method with a double-couple seismic source approximation, scaled for the source dimensions and relatively simple regional (bedrock) and local structure models, we succeeded in reproducing the recorded ground motion in Bucharest at a satisfactory level for seismic engineering. Extending the modeling to the entire territory of the Bucharest area, we construct a new seismic microzonation map, where five different zones are identified by their characteristic response spectra.     

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.     

Dynamic properties of the Quaternary sedimentary rocks and their influence on seismic site effects. Case study in Bucharest City,Romania

Bucharest is one of the cities most affected by earthquakes in Europe. Situated at 150–170 km distance from Vrancea epicentral zone, Bucharest had suffered many damages due to high energy Vrancea intermediate-depth earthquakes. For example, the 4 March 1977 event produced the collapse of 32 buildings with 8–12 levels, while more than 150 old buildings with 6–9 levels were seriously damaged. The studies done after this earthquake had shown the importance of the surface geological structure upon ground motion parameters. New seismic measurements are performed in Bucharest area aiming at defining better elastic and dynamic properties of the shallow sedimentary rocks. Down-hole seismic measurements were performed in a number of 10 cased boreholes drilled in the Bucharest City area. Processing and interpretation of the data lead to the conclusion that shallow sedimentary rocks can be considered weak in the area, down to 150–200 m depth. Seismic wave velocity values and bulk density values presented in the paper associated with local geology are useful primary data in the seismic microzonation of Bucharest City. They are used as 1D models to derive transfer functions and response spectra for the stack of sedimentary rocks in several parts of Bucharest area, leading to a better knowledge of the local site amplification and associated frequency spectra. In a recent study the H/V spectral ratio using Nakamura’s method was applied on the seismic noise measurements in 22 sites in Bucharest City in order to derive the fundamental period associated with these sites. The values confirm the previous results, showing a dominant resonance in the period range of 1.25–1.75 s. The fundamental periods obtained with Nakamura’s method are in good agreement with those computed on the basis of geological and geotechnical data in boreholes, which show an increase of the fundamental period in the Bucharest area from south to north, in the same direction as the increase of the thickness of the Quaternary deposits above the Fratesti layer which is considered the bedrock in the area.     

Acceleration Time Histories for a Scenario Earthquake in Moscow at Sites with Different Soil Conditions

According to general seismic zoning maps, Moscow is in an area with a seismic intensity of 5, in which the maximum seismic effect is expected from remote deep-focal earthquakes in the Vrancea zone (Eastern Carpathians, Romania). In our previous studies, an earthquake with a hypocenter at a depth of 80–150 km in the Vrancea zone, a moment magnitude of M_w = 8.0, and a drop in stress of Δσ = 325 bar was used as a scenario earthquake for Moscow. A series of model acceleration time histories for ground vibrations was calculated for this earthquake for the reference local conditions of the Moskva seismic station (Moscow, Pyzhevskii per. 3). In this paper, these acceleration time histories are used to calculate the acceleration time histories and estimate the ground vibration parameters for an scenario earthquake at other sites on the territory of Moscow for which information on soil conditions is available. Since the epicentral distance is large (~1300 km), it can be assumed that changes in the shape and spectral content of the acceleration time histories at different sites in Moscow are only caused by different local conditions.     

A uniform approach to seismic site effect analysis in Bucharest, Romania

As a uniform approach to the assessment of ground motion variation within the Romanian capital Bucharest we analyze and compare strong motion records from analog recorders, weak motion data from a modern digital accelerometer network, and intensity observations of previous strong earthquakes. These different data sets allow to clearly characterize geographical trends in the distribution of ground shaking in the city for future earthquakes. Below 2 Hz the variability is small. Between 2 and 5 Hz, however, variations by a factor of 3–4 have to be expected. As the key source for the seismic hazard—the intermediate depth Vrancea earthquakes—remain at hypocentral distances in excess of 150 km from the city the ground motion variation must be predominantly attributed to site effects. This geometry of Vrancea sources to the site of Bucharest is ideal for the application of source-site separation techniques. However, despite this fact site effect amplification functions display a very large amount of aleatory uncertainty. In other words the standard source-site parameterization is too simple and we do not yet fully understand the cause and size of site effects.     

Simulation of strong earthquake characteristics of a scenario earthquake (MS7.5) based on the enlightenment of 2022 MS6.9 earthquake in Menyuan

The Menyuan area is an important transportation hub in the Hexi Corridor. The Menyuan M_S6.9 earthquake that occurred on January 8, 2022 had a major impact on the local infrastructure and transportation of this region. Due to the high possibility of similar strong earthquakes occurring in this area in the future, preliminary assessment of the seismic intensity characteristics of destructive earthquakes in this region is essential for effective disaster control. This paper uses the empirical Green′s function (EGF) method as a numerical simulation tool to predict the ground motion intensity of Datong Autonomous County under the action of the scenario earthquake (M_S7.5). Seismic records of aftershocks of the 2016 Menyuan M_S6.4 earthquake were used as Green’s functions for this simulation. The uncertainties associated with various source parameters were considered, and 36 possible earthquake scenarios were simulated to obtain 72 sets of horizontal ground motions in Datong County. The obtained peak ground acceleration (PGA) vs. time histories of the horizontal ground motion were screened using the attenuation relationships provided by the fifth-edition of China's Seismic Ground Motion Parameter Zoning Map and the NGA-West2 dataset. Ultimately, 32 possible acceleration-time histories were selected for further analysis. The screened PGA values ranged from 78.8 to 153 cm/s². The uncertainty associated with the initial rupture point was found to greatly affect the results of the earthquake simulation. The average acceleration spectrum of the selected acceleration-time history exceeded the expected spectrum of a intermediate earthquake, which means that buildings in Datong County might sustain some damage should the scenario earthquake occur. This research can provide reliable ground motion input for urban earthquake damage simulation and seismic design in Datong County. Growing the dataset of small earthquakes recorded in this region will facilitate the large-scale simulation of ground motions under different earthquake scenarios.     

A preliminary empirical model for frequency-dependent attenuation of Fourier amplitude spectra In Serbia from the Vrancea earthquakes

We present the frequency-dependent attenuation model for Fourier amplitude spectra of strong earthquake ground motion in Serbia from intermediate depth earthquakes in the Vrancea source zone in Romania. The development of this type of scaling is the essential first step toward developing the corresponding attenuation and scaling equations for pseudo relative velocity spectra (PSV), which are necessary for seismic macro- and microzoning in the territory of Serbia. Such scaling is necessary because the Vrancea source zone produces large earthquakes with shaking that attenuates differently from the local earthquakes in Serbia. Development of such a scaling model is associated with several difficulties, the principal one being the lack of recorded strong motion accelerograms at epicentral distances exceeding 300 km. To reduce uncertainties with such scaling, we require our preliminary scaling equations to be consistent with independent estimates of seismic moment, stress drop, and radiated wave energy. In the future, when the recorded strong motion data from Vrancea earthquakes increases several-fold of what it is today, it will become possible to perform this analysis again, thus leading to more reliable and permanent scaling estimates.     

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.     

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.     

Simulation of ground motion in the Moscow region using the empirical Green’s function

Historically, the Moscow region regularly experienced rather weak but quite perceptible seismic vibrations produced by intermediate-depth earthquakes of the Vrancea zone (Romania), located at a distance of 1400 km from Moscow. The coincidence of a number of unique factors such as a slowly varying focal depth, predominant source mechanisms, weak attenuation of seismic radiation in the north-northeast direction provide favorable conditions for application of the empirical Green’s function method. Using the digital seismogram of the Vrancea Mw-5.8 earthquake recorded at the Moscow seismic station, we simulated synthetic seismograms of a scenario (expected maximum) earthquake with Mw = 8.0, by application of the empirical Green’s function method adjusted for the given conditions. The calculation procedure was verified using analog records of strong earthquakes available at the Moscow seismic station. Digital records of the Obninsk seismic station included in the Incorporated Research Institutions for Seismology (IRIS) system were used for additional control. Here, the scenario earthquake was modeled using the data on a much stronger earthquake of 1990 (MW = 6.9). It is shown that, despite a certain scatter (quite adequately assessed in the scope of the method), the ultimate estimates of expected seismic impacts are quite reliable and can be recommended for practical use.     

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.     

Investigation on regional attenuation of Vrancea (Romania) intermediate-depth earthquakes

This paper aims at investigating possible regional attenuation patterns in the case of Vrancea(Romania) intermediate-depth earthquakes.Almost 500 pairs of horizontal components recorded during 13 intermediate-depth Vrancea earthquakes are employed in order to evaluate the regional attenuation patterns.The recordings are grouped according to the azimuth with regard to the Vrancea seismic source and subsequently,Q models are computed for each azimuthal zone assuming similar geometrical spreading.Moreover,the local soil amplification which was disregarded in a previous analysis performed for Vrancea intermediate-depth earthquakes is now clearly evaluated.The results show minor differences between the four regions situated in front of the Carpathian Mountains and considerable differences in attenuation of seismic waves between the forearc and backarc regions(with regard to the Carpathian Mountains).Consequently,an average Q model of the type Q(f) = 115×f~(1.25) is obtained for the four forearc regions,while a separate Q model of the type Q(f) = 70×f~(0.90) is computed for the backarc region.These results highlight the need to evaluate the seismic hazard of Romania by using ground motion models which take into account the different attenuation between the forearc/backarc regions.     

Prediction of the seismic manifestations of Vrancea earthquakes in Moscow

According to the normative maps of the General Seismic Zoning in the Russian Federation, OSR-97, the Moscow metropolitan area is situated within the 5 point seismic zone. Of highest hazard priority for tall buildings in Moscow are the low-frequency vibrations proceeding from the deep sources of strong earthquakes that occur in the East Carpathians (the Vrancea zone, Romania) at a distance of approximately 1350 km from Moscow. Accelerations of the ground vibrations in Moscow are found from the analysis of seismic signals produced by Mw = 5.0 to Mw = 7.4 Vrancea earthquakes and recorded at the Moskva seismic station. Extrapolation of the parameters of the weak and moderate earthquakes towards stronger seismic events provides an estimate for the maximum expected horizontal accelerations of A_hor = 2.3 cm/s² in case of the Mw = 8.0 Vrancea earthquake. The synthetic accelerogram of the maximum possible effect on the benchmark soils of Moscow is calculated. The displacements of the ground are multidimensional and not necessarily oriented strictly towards the seismic source. These inferences suggest that the MSK-64 macroseismic scale be corrected and the Construction Norms and Regulations, SNIP II-7-81*, be updated with regard to the hazard assessment of low-frequency seismic effects of 5 point and weaker seismic events including those caused by distant earthquakes.     

基于条件均值反应谱的特大地震强震记录的选取及调整方法

在PEER地震动数据库(PGMD)的基础上, 结合近几年国内外特大地震的地面运动记录, 建立了地面运动数据库, 同时根据日本M_W9.0特大地震获得的141组记录进行统计回归建立加速度反应谱衰减关系, 并采用条件均值反应谱法, 即设定地震与结构概率需求结合的方法选择地面运动. 选波实例表明, 当设定地震为特大地震时, 基于条件均值反应谱法选取地面运动记录时, 扩展数据库中大震记录并建立符合大震记录加速度反应谱的衰减关系是十分必要与迫切的. 该思路为进一步研究结构动态时程分析中地面运动记录选取问题及所选记录提供了依据.      

Improved seismic risk estimation for Bucharest,based on multiple hazard scenarios and analytical methods

Bucharest, capital of Romania, is one of the most exposed big cities in Europe to seismic damage, due to the intermediate-depth earthquakes in the Vrancea region, to the vulnerable building stock and local soil conditions.This paper tries to answer very important questions related to the seismic risk at city scale that were not yet adequately answered. First, we analyze and highlight the bottlenecks of previous risk-related studies. Based on new researches in the hazard of Bucharest (recent microzonation map and ground-motion prediction equations, reprocessed real recorded data) and in vulnerability assessment (analytical methods, earthquake loss estimation software like SELENA and ELER, the recently implemented Near Real-Time System for Estimating the Seismic Damage in Romania) we provide an improved estimation of the number of buildings and population that could be affected, for different earthquake scenarios. A new method for enhancing the spatial resolution of the building stock data is used successfully.     

Simulation of Ground Motion from Strong Earthquakes of Kamchatka Region (1992–1993) at Rock and Soil Sites

To estimate the parameters of ground motion in future strong earthquakes, characteristics of radiation and propagation of seismic waves in the Kamchatka region were studied. Regional parameters of radiation and propagation of seismic waves were estimated by comparing simulations of earthquake records with data recorded by stations of the Kamchatka Strong Motion Network. Acceleration time histories of strong earthquakes (M _w = 6.8–7.5, depths 45–55 km) that occurred near the eastern coast of Kamchatka in 1992–1993 were simulated at rock and soil stations located at epicentral distances of 67–195 km. In these calculations, the source spectra and the estimates of frequency-dependent attenuation and geometrical spreading obtained earlier for Kamchatka were used. The local seismic-wave amplification was estimated based on shallow geophysical site investigations and deep crustal seismic explorations, and parameters defining the shapes of the waveforms, the duration, etc. were selected, showing the best-fit to the observations. The estimated parameters of radiation and propagation of seismic waves describe all the studied earthquakes well. Based on the waveforms of the acceleration time histories, models of slip distribution over the fault planes were constructed for the studied earthquakes. Station PET can be considered as a reference rock station having the minimum site effects. The intensity of ground motion at the other studied stations was higher than at PET due to the soil response or other effects, primarily topographic ones. At soil stations INS, AER, and DCH the parameters of soil profiles (homogeneous pyroclastic deposits) were estimated, and nonlinear models of their behavior in the strong motion were constructed. The obtained parameters of radiation and propagation of seismic waves and models of soil behavior can be used for forecasting ground motion in future strong earthquakes in Kamchatka.     

Seismic Hazard and Ground Motion Characterization at the Itoiz Dam (Northern Spain)

This paper presents a new hazard-consistent ground motion characterization of the Itoiz dam site, located in Northern Spain. Firstly, we propose a methodology with different approximation levels to the expected ground motion at the dam site. Secondly, we apply this methodology taking into account the particular characteristics of the site and of the dam. Hazard calculations were performed following the Probabilistic Seismic Hazard Assessment method using a logic tree, which accounts for different seismic source zonings and different ground-motion attenuation relationships. The study was done in terms of peak ground acceleration and several spectral accelerations of periods coinciding with the fundamental vibration periods of the dam. In order to estimate these ground motions we consider two different dam conditions: when the dam is empty (T?=?0.1?s) and when it is filled with water to its maximum capacity (T?=?0.22?s). Additionally, seismic hazard analysis is done for two return periods: 975?years, related to the project earthquake, and 4,975?years, identified with an extreme event. Soil conditions were also taken into account at the site of the dam. Through the proposed methodology we deal with different forms of characterizing ground motion at the study site. In a first step, we obtain the uniform hazard response spectra for the two return periods. In a second step, a disaggregation analysis is done in order to obtain the controlling earthquakes that can affect the dam. Subsequently, we characterize the ground motion at the dam site in terms of specific response spectra for target motions defined by the expected values SA (T) of T?=?0.1 and 0.22?s for the return periods of 975 and 4,975?years, respectively. Finally, synthetic acceleration time histories for earthquake events matching the controlling parameters are generated using the discrete wave-number method and subsequently analyzed. Because of the short relative distances between the controlling earthquakes and the dam site we considered finite sources in these computations. We conclude that directivity effects should be taken into account as an important variable in this kind of studies for ground motion characteristics.     

Earthquake ground motion predictive equations for Garhwal Himalaya,India

Predictive equations based on the stochastic approach are developed for earthquake ground motions from Garhwal Himalayan earthquakes of 3.5≤M_w≤6.8 at a distance of 10≤R≤250 km. The predicted ground motion parameters are response spectral values at frequencies from 0.25 to 20 Hz, and peak ground acceleration (PGA). The ground motion prediction equations (GMPEs) are derived from an empirically based stochastic ground motion model. The GMPEs show a fair agreement with the empirically developed ground motion equations from Himalaya as well as the NGA equation. The proposed relations also reasonably predict the observed ground motion of two major Himalayan earthquakes from Garhwal Himalayan region. For high magnitudes, there is insufficient data to satisfactorily judge the relationship; however it reasonably predicts the 1991 Uttarkashi earthquake (M_w=6.8) and 1999 Chamoli earthquake (M_w=6.4) from Garhwal Himalaya region.     

Investigation on the stochastic simulation of strong ground motions for Bucharest area

This short article evaluates the stochastic method of ground motion simulation for Bucharest area using both the single-corner frequency model and recently introduced double-corner frequency models. A dedicated Q model is derived using ground motions obtained during the largest Vrancea earthquakes from the past 30 years. The simulated ground motions are tested against the observed data from the Vrancea earthquakes of August 1986 and May 1990. Moreover, the observed data are also compared against simulations obtained using the Q model derived by Oth et al. (2008). Finally, the results of the simulations show that the derived Q model is better suited for larger magnitude events, while the Q model of Oth et al. (2008) provides better results for smaller earthquakes.