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.     

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.     

Macroseismic field of the October 27, 2004 Vrancea (Romania) moderate subcrustal earthquake

On October 27, 2004, a moderate size earthquake occurred in the Vrancea seismogenic region (Romania). The Vrancea seismic zone is an area of concentrated seismicity at intermediate depths beneath the bending area of the southeastern Carpathians. The 2004 M _w?=?6 Vrancea subcrustal earthquake is the largest seismic event recorded in Romania since the 1990 earthquakes. With a maximum macroseismic intensity of VII Medvedev–Sponheuer–Kárník (MSK-64) scale, the seismic event was felt to a distance of 600 km from the epicentre. This earthquake caused no serious damage and human injuries. The main purpose of this paper is to present the macroseismic map of the earthquake based on the MSK-64 intensity scale. After the evaluation of the macroseismic effects of this earthquake, an intensity dataset has been obtained for 475 sites in the Romanian territory. Also, the maximum horizontal accelerations recorded in the area by the K2 network are compared to the intensity values.     

Criteria for Selecting and Adjusting Ground-Motion Models for Specific Target Regions: Application to Central Europe and Rock Sites

A vital component of any seismic hazard analysis is a model for predicting the expected distribution of ground motions at a site due to possible earthquake scenarios. The limited nature of the datasets from which such models are derived gives rise to epistemic uncertainty in both the median estimates and the associated aleatory variability of these predictive equations. In order to capture this epistemic uncertainty in a seismic hazard analysis, more than one ground-motion prediction equation must be used, and the tool that is currently employed to combine multiple models is the logic tree. Candidate ground-motion models for a logic tree should be selected in order to obtain the smallest possible suite of equations that can capture the expected range of possible ground motions in the target region. This is achieved by starting from a comprehensive list of available equations and then applying criteria for rejecting those considered inappropriate in terms of quality, derivation or applicability. Once the final list of candidate models is established, adjustments must be applied to achieve parameter compatibility. Additional adjustments can also be applied to remove the effect of systematic differences between host and target regions. These procedures are applied to select and adjust ground-motion models for the analysis of seismic hazard at rock sites in West Central Europe. This region is chosen for illustrative purposes particularly because it highlights the issue of using ground-motion models derived from small magnitude earthquakes in the analysis of hazard due to much larger events. Some of the pitfalls of extrapolating ground-motion models from small to large magnitude earthquakes in low seismicity regions are discussed for the selected target region.     

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.     

Modelling of Block-structure Dynamics: Parametric Study for Vrancea

—The dependence of synthetic earthquake catalog obtained by numerical modelling of block-structure dynamics for the Vrancea (Romania) seismoactive region on values of the model parameters is studied. The features of the synthetic seismicity, such as the spatial distribution of epicentres, the level of seismic activity, the relative activity of the different faults, are considered as functions of the directions of motions of the different blocks of the structure. The likelihood of the features of the synthetic and of the real seismicity could be used as criterion for the reconstruction of the direction of the tectonic motions.     

A new seismic hazard assessment in the region of Andalusia (Southern Spain)

A probabilistic seismic hazard assessment of Andalusia (Southern Spain) in terms of peak ground acceleration, PGA, and spectral accelerations, SA(T), is presented in this paper. In contrast to most of the previous studies in the region, which were performed for PGA, making use of Intensity-to-PGA relationships, hazard was here calculated in terms of magnitude, using published spectral ground-motion models. Moreover, we considered different ground-motion models for the Atlantic sources, since the attenuation of those motions seems to be slower, as evidenced in the case of the extensive macroseismic areas of earthquakes like those occurred in the years 1755, 1969 and 2007. A comprehensive review of the seismic catalogue and of the seismogenic models proposed for the region was carried out, including those for Northern Africa, which is part of the influence area. Hazard calculations were performed following the Probabilistic Seismic Hazard Assessment (PSHA) methodology using a logic tree, which accounts for six different seismic source zonings and five different ground-motion attenuation relationships. Hazard maps in terms of PGA and SA (0.2 s) and SA (1 s) and coefficient of variation (COV) maps, for the 475-year return period were first obtained in rock sites. A geotechnical classification and amplification factors were proposed and new hazard maps including local effects were represented, showing PGA values ranging from 24 to 370 cm/s² for the whole Andalusian territory, with the highest expected values (PGA > 300 cm/s²) in some parts of the Granada Province and in the town of Vélez Málaga. Lowest values (PGA < 50 cm/s²) correspond to some towns of the Huelva and Córdoba provinces. The inclusion of soil effects provides a more detailed picture of the actual hazard the region is subjected to.     

Toward a ground-motion logic tree for probabilistic seismic hazard assessment in Europe

The Seismic Hazard Harmonization in Europe (SHARE) project, which began in June 2009, aims at establishing new standards for probabilistic seismic hazard assessment in the Euro-Mediterranean region. In this context, a logic tree for ground-motion prediction in Europe has been constructed. Ground-motion prediction equations (GMPEs) and weights have been determined so that the logic tree captures epistemic uncertainty in ground-motion prediction for six different tectonic regimes in Europe. Here we present the strategy that we adopted to build such a logic tree. This strategy has the particularity of combining two complementary and independent approaches: expert judgment and data testing. A set of six experts was asked to weight pre-selected GMPEs while the ability of these GMPEs to predict available data was evaluated with the method of Scherbaum et al. (Bull Seismol Soc Am 99:3234?C3247, 2009). Results of both approaches were taken into account to commonly select the smallest set of GMPEs to capture the uncertainty in ground-motion prediction in Europe. For stable continental regions, two models, both from eastern North America, have been selected for shields, and three GMPEs from active shallow crustal regions have been added for continental crust. For subduction zones, four models, all non-European, have been chosen. Finally, for active shallow crustal regions, we selected four models, each of them from a different host region but only two of them were kept for long periods. In most cases, a common agreement has been also reached for the weights. In case of divergence, a sensitivity analysis of the weights on the seismic hazard has been conducted, showing that once the GMPEs have been selected, the associated set of weights has a smaller influence on the hazard.     

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.     

Seismic hazard assessment of the Province of Murcia (SE Spain): analysis of source contribution to hazard

A probabilistic seismic hazard assessment of the Province of Murcia in terms of peak ground acceleration (PGA) and spectral accelerations [SA(T)] is presented in this paper. In contrast to most of the previous studies in the region, which were performed for PGA making use of intensity-to-PGA relationships, hazard is here calculated in terms of magnitude and using European spectral ground-motion models. Moreover, we have considered the most important faults in the region as specific seismic sources, and also comprehensively reviewed the earthquake catalogue. Hazard calculations are performed following the Probabilistic Seismic Hazard Assessment (PSHA) methodology using a logic tree, which accounts for three different seismic source zonings and three different ground-motion models. Hazard maps in terms of PGA and SA(0.1, 0.2, 0.5, 1.0 and 2.0 s) and coefficient of variation (COV) for the 475-year return period are shown. Subsequent analysis is focused on three sites of the province, namely, the cities of Murcia, Lorca and Cartagena, which are important industrial and tourism centres. Results at these sites have been analysed to evaluate the influence of the different input options. The most important factor affecting the results is the choice of the attenuation relationship, whereas the influence of the selected seismic source zonings appears strongly site dependant. Finally, we have performed an analysis of source contribution to hazard at each of these cities to provide preliminary guidance in devising specific risk scenarios. We have found that local source zones control the hazard for PGA and SA(T ≤ 1.0 s), although contribution from specific fault sources and long-distance north Algerian sources becomes significant from SA(0.5 s) onwards.     

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.     

Seismic Hazard Estimate at the Iberian Peninsula

—?Seismic hazard at the Iberian Peninsula has been evaluated by using a methodology which combines both zonified and non-zonified probabilistic methods. Seismic sources are used when considering zones where certain calculation parameters may be considered homogeneous, as in zonified methods, while, on the other hand, earthquakes are considered wherever it has taken place, as in non-zonified methods. The methodology which is applied in this paper has been originally used to calculate the seismic hazard maps in the United States. In our case, it has been necessary to adapt the method to the specific features of the seismicity in the Iberian Peninsula and its geographical surroundings, not only with respect to its distribution and characteristics, but also with respect to the properties of the seismic catalog used.¶Geographically, the main feature of the result is the fact that it reflects both historical seismicity and current seismic clusters of the region. Despite the smoothing, maps show marked differences between several seismic zones; these differences becoming more noticeable as exposure time increases. Maximum seismic hazard is found to be in the southwestern region of the Peninsula, especially in the area of the Cape St. Vicent, and around Lisbon. The uncertainty of the results, without considering that due to the attenuation laws, as deduced from the other evaluation parameters, is quite stable, being more sensitive to the parameters b and m _max of the Gutenberg-Richter relation.     

Deterministic Earthquake Scenarios for the City of Sofia

— The city of Sofia is exposed to a high seismic risk. Macroseismic intensities in the range of VIII – X (MSK) can be expected in the city. The earthquakes that can influence the hazard in Sofia originate either beneath the city or are caused by seismic sources located within a radius of 40 km. The city of Sofia is also prone to the remote Vrancea seismic zone in Romania, and particularly vulnerable are the long-period elements of the built environment. The high seismic risk and the lack of instrumental recordings of the regional seismicity make the use of appropriate credible earthquake scenarios and ground-motion modelling approaches for defining the seismic input for the city of Sofia necessary. Complete synthetic seismic signals, due to several earthquake scenarios, were computed along chosen geological profiles crossing the city, applying a hybrid technique, which combines the modal summation technique and finite differences. The modelling takes into account simultaneously the geotechnical properties of the site, the position and geometry of the seismic source and the mechanical properties of the propagation medium. Acceleration, velocity and displacement time histories and related quantities of earthquake engineering interest (e.g., response spectra, ground-motion amplification along the profiles) have been supplied. The approach applied in this study allows us to obtain the definition of the seismic input at low cost, exploiting large quantities of existing data (e.g. geotechnical, geological, seismological). It may be efficiently used to estimate the ground motion for the purposes of microzonation, urban planning, retrofitting or insurance of the built environment, etc.     

Numerical simulation of the earthquake generation process

A numerical algorithm is proposed for the simulation of the earthquake process during a seismic cycle. The algorithm is based on a heterogeneous discrete model of the fault plane and assumes there are two kinds of seismicity: background crack-like earthquakes and asperity-like events. An active zone of the fault contains an asperity distribution with a characteristic elementary area. The background seismicity randomly develops shear stress-free surfaces which tend to surround the asperities as in a 2D percolation process. The model parameters are taken from observations on the Vrancea (Romania) intermediate depth seismic region. The results emphasize the significant role of the geometry in the mechanism of the seismic failure. The algorithm predicts the nonlinear behavior in the frequency-magnitude distribution, the decrease of theb-slope associated with the asperity-like events, the magnitude range of major earthquakes, and their recurrence times.     

The application of seismic data with different precision in the determination of seismicity parameters

ＴｈｅａｐｐｌｉｃａｔｉｏｎｏｆｓｅｉｓｍｉｃｄａｔａｗｉｔｈｄｉｆｆｅｒｅｎｔｐｒｅｃｉｓｉｏｎｉｎｔｈｅｄｅｔｅｒｍｉｎａｔｉｏｎｏｆｓｅｉｓｍｉｃｉｔｙｐａｒａｍｅｔｅｒｓＸｕｅ－ＳｈｅｎＪＩＮ（金学申）；Ｙｉｎｇ－ＨｕａＤＡＩ（戴英华）；Ｊｕ...     

Current empirical ground-motion prediction equations for Europe and their application to Eurocode 8

The first ground-motion prediction equation derived from European and Middle Eastern strong-motion data was published more than 30 years ago; since then strong-motion networks and the resulting databank of accelerograms in the region have expanded significantly. Many equations for the prediction of peak ground-motion parameters and response spectral ordinates have been published in recent years both for the entire Euro-Mediterranean and Middle Eastern region as well as for individual countries within this region. Comparisons among empirical ground-motion models for these parameters, developed using large regional datasets, do not support the hypothesis of there being significant differences in earthquake ground-motions from one area of crustal seismicity to another. However, there are certain regions within Europe—affected by different tectonic regimes—for which the existing pan-European equations may not be applicable. The most recent European equations make it possible to now implement overdue modifications to the presentation of seismic design actions in Eurocode 8 that allow an improved approximation to the target uniform hazard spectrum (UHS). Using these recent equations, this study outlines a new approach via which an approximation to the UHS may be constructed using hazard maps calculated for peak ground velocity and the corner period T _D in addition to the maps for peak ground acceleration that underpin the current stipulations of Eurocode 8.     

青藏高原东北缘地震活动性广义帕累托模型的全域敏感性分析

基于广义帕累托分布构建地震活动性模型,因其输入参数取值难以避免不确定性,导致依据该模型所得的地震危险性估计结果具有不确定性。鉴于此,本文选取青藏高原东北缘为研究区,提出了基于全域敏感性分析的地震危险性估计的不确定性分析流程和方法。首先,利用地震活动性广义帕累托模型,进行研究区地震危险性估计;然后,选取地震记录的起始时间和震级阈值作为地震活动性模型的输入参数,采用具有全域敏感性分析功能的E-FAST方法,对上述两个参数的不确定性以及两参数之间的相互作用对地震危险性估计不确定性的影响进行定量分析。结果表明：地震危险性估计结果（不同重现期的震级重现水平、震级上限及相应的置信区间）对两个输入参数中的震级阈值更为敏感;不同重现期的地震危险性估计结果对震级阈值的敏感程度不同;对不同的重现期而言,在影响地震危险性估计结果的不确定性上,两个输入参数之间存在非线性效应,且非线性效应程度不同。本文提出的不确定性分析流程和方法,可以推广应用于基于其它类型地震活动性模型的地震危险性估计不确定性分析。     

应用地震层析成象技术研究深部速度结构

罗马尼亚Ｖａｎｃｅａ地震区是大陆上发生与板块磁撞和削减有关的中深部地震活动的地区之一。本文介绍了应用地震层析成象方法研究该地区深部速度结构成果。在研究中使用了地方和区域地震所记录的４３３个浅源和中部地震的到时资料反演求解深至２００ｋｍ的三维速度结构,在走时和射线路径的计算中利用了有效的三维射线跟踪技术,在反演中采用ＬＳＱＲ算法,高分辨率的地震层析图象揭示了速度结构的广泛不均匀性,结果表明,地震层析     

Comparison of conventional and Monte Carlo simulation-based probabilistic seismic hazard analyses for Shiraz city,southern Iran

Estimation of ground-motion amplitudes of different hazard levels is of paramount importance in planning of urban development of any metropolis. Such estimation can be computed through a probabilistic seismic hazard analysis (PSHA). This paper concentrates on the PSHA of an area located in Shiraz city, southern Iran. The area includes whole of Shiraz city (i.e., one of the largest and most populous cities of Iran) and its outskirts. Conventional and Monte Carlo simulation-based approaches are utilized to perform the PSHA of the studied area. Two areal seismic source models are delineated, and thence seismicity parameters of all zones associated with their corresponding uncertainties are computed. Uncertainties in ground-motion prediction are accounted for via three ground-motion prediction equations (GMPEs) within the logic tree framework. These GMPEs are applied to estimate bedrock ground shaking (V_s30?=?760 m/s) for several return periods (i.e., 75, 475, 975, and 2475 years). In general, the results of the two abovementioned PSHA approaches show relatively similar results. However, the Monte Carlo simulation-based approach overpredicts bedrock spectral accelerations at periods of 0.4–2.5 s compared to the conventional PSHA approach for return periods of 475, 975, and 2475 years.     

A macroseismic intensity prediction equation for intermediate depth earthquakes in the Vrancea region,Romania

The use of shake maps in terms of macroseismic intensity in earthquake early warning systems as well as intensity based seismic hazard assessments provides a valuable supplement to typical studies based on recorded ground motion parameters. A requirement for such applications is ground motion prediction equations (GMPE) in terms of macroseismic intensity, which have the advantages of good data availability and the direct relation of intensity to earthquake damage. In the current study, we derive intensity prediction equations for the Vrancea region in Romania, which is characterized by the frequent occurrence of large intermediate depth earthquakes giving rise to a peculiar anisotropic ground shaking distribution. The GMPE have a physical basis and take the anisotropic intensity distribution into account through an empirical regional correction function. Furthermore, the relations are easy to implement for the user. Relations are derived in terms of epicentral, rupture and Joyner–Boore distance and the obtained relations all provide a new intensity estimate with an uncertainty of ca. 0.6 intensity units.