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.     

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.     

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.     

Geomagnetic disturbances and earthquakes in the Vrancea zone

Based on the example of the Vrancea zone of concentrated seismicity, it is shown how the stress-strain state of the medium responds to a disturbance of the geomagnetic field. Geomagnetic conditions are examined in relation to earthquakes in the Vrancea zone in the period 1988–1996. It is established that the seismic energy release in the Vrancea zone is associated with differences (“gradients”) in the H component of the geomagnetic field. Such a gradient preceding earthquakes is shown to be the midnight polar substorm and the degree of its mid-latitude effect. The time interval from the maximum of the substorm development to a shock (τ, h) is directly related to the focal depth. The seismic characteristics K _en and h (km) are demonstrated to be related to morphological features of the substorm development, namely, its duration T (min), intensity, and background. Differences in the duration of polar substorms before crustal (shallow) and deep earthquakes are revealed. Morphological features of the spectrum of geomagnetic variations preceding the seismic energy release are established.     

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.     

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.     

Probabilistic seismic hazard assessment for Romania and sensitivity analysis: A case of joint consideration of intermediate-depth (Vrancea) and shallow (crustal) seismicity

The earthquake risk on Romania is one of the highest in Europe, and seismic hazard for almost half of the territory of Romania is determined by the Vrancea seismic region, which is situated beneath the southern Carpathian Arc. The region is characterized by a high rate of occurrence of large earthquakes in a narrow focal volume at depth from 70 to 160 km. Besides the Vrancea area, several zones of shallow seismicity located within and outside the Romanian territory are considered as seismically dangerous. We present the results of probabilistic seismic hazard analysis, which implemented the “logic tree” approach, and which considered both the intermediate-depth and the shallow seismicity. Various available models of seismicity and ground-motion attenuation were used as the alternative variants. Seismic hazard in terms of macroseismic intensities, peak ground acceleration, and response spectra was evaluated for various return periods. Sensitivity study was performed to analyze the impact of variation of input parameters on the hazard results. The uncertainty on hazard estimates may be reduced by better understanding of parameters of the Vrancea source zone and the zones of crustal seismicity. Reduction of uncertainty associated with the ground-motion models is also very important issue for Romania.     

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.     

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.     

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.     

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.     

Correlations of peak acceleration,velocity and displacement with earthquake magnitude,distance and site conditions

A brief review of proposed correlations between peak accelerations and earthquake magnitude and distance has been presented. It has been found that most investigators agree favourably on what should be the amplitude of peak accelerations for the distance range between about 20 and 200 km. For distances less than 20 km, there is significant disagreement in the predicted peak amplitudes, reflecting the lack of data there and the uncertainties associated with the extrapolation. Correlations of peak accelerations, peak velocities and peak displacements with earthquake magnitude, epicentral distance and the geologic conditions of the recording sites have been presented for 187 accelerograms recorded during 57 earthquakes. This data set describes strong earthquake ground motion in the Western United States during the period from 1933 to 1971. For large earthquakes, dependence of peak acceleration, velocity and displacement amplitudes on earthquake magnitude seems to be lost. This suggests that the amplitudes of strong ground motion close to a fault are scaled primarily by the maximum dislocation amplitudes and the stress drop, rather than the overall ‘size’ of an earthquake as measured by magnitude. The influence of geologic conditions at the recording station seems to be of minor importance for scaling peak accelerations, but it becomes noticeable for the peaks of velocity and even more apparent for the peaks of displacement.     

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.     

Increment of Macroseismic Intensity at Seismic Stations of the Kamchatka Region Relative to the Petropavlovsk Seismic Station

The paper presents the method and results of calculating the increment of macroseismic intensity at seismic stations of Kamchatka. Calculation is based on measurement of the relative level of maximum accelerations of intense earth vibrations in the phase of S-waves of comparatively strong regional earthquakes and the root-mean-square deviation of acceleration in the phase of P-waves of a strong distant earthquake. In the latter case, records of an earthquake with a magnitude of 9.1, which occurred in Japan on March 11, 2011, were used. The Petropavlovsk seismic station was used as the reference station. At the foundation of this station rests on rocky soil composed of siliceous shales. An estimate of the increment for the majority of digital stations is presented. Anomalously high intensity values were noted at a number of stations. The data obtained are used to assess the properties of soils in the investigated area. At several stations, the intensity of the horizontal component of soil vibrations above the intensity of the vertical component is much greater than the corresponding design value, which is probably due to the presence of resonant soil layers under these stations. The discrepancy in the incremental intensity estimates from records of intense oscillations from regional earthquakes and from records of a very strong remote earthquake obtained from sensors located in basements of heavy-frame concrete structures is revealed. To avoid distortion in recording ground vibrations, it is desirable to place seismic instruments far from such structures. The results obtained in the study can be used for seismic microzoning of construction sites in the investigated territory.     

Three-dimensional liquefaction potential analysis using geostatistical interpolation

We applied three-dimensional geostatistical interpolation to evaluate the extent of liquefiable materials at two sites that liquefied during the 1994 Northridge Earthquake. The sites were the Balboa Blvd site and the Wynne Ave. site located in the alluvial San Fernando Valley. The estimated peak ground accelerations at the sites are 0.84 g (Balboa Blvd) and 0.51 g (Wynne Ave.). These sites were chosen because surface effects due to liquefaction were not predicted using available techniques based on thickness and depth of liquefiable layers (Ishihara [Ishihara K. Stability of natural deposits during earthquakes. Proceedings of the 11th international conference on soil mechanics and foundation engineering, vol. 1. Rotterdam, The Netherlands: A.A. Balkema; 1985. p. 321–76.]) and the Liquefaction Potential Index (Iwasaki et al. [Iwasaki T, Tatsuoka F, Tokida K, Yasuda S. A practical method for assessing soil liquefaction potential based on case studies at various sites in Japan. In: Proceedings of the second international conference on microzonation, San Francisco; 1978. p. 885–96.]). During the earthquake, both sites experienced surface effects including ground cracking and extension as a result of liquefaction. Foundations and buried utilities were damaged at both sites. The sites were investigated after the event by researchers with the United States Geologic Survey using standard penetration tests (SPT) and cone penetration tests. In this paper, liquefaction potential was estimated for each soil sample using results from SPTs according to the updated Seed and Idriss simplified procedure. The probability of liquefaction was estimated by applying an indicator transform to the results of the liquefaction potential calculation. We compared our results to detailed geologic mapping of the sites performed by other researchers. Using geostatistical interpolation to estimate the probability of liquefaction is a useful supplement to geologic evaluation of liquefaction potential. The geostatistical analysis provides an estimate of the continuous volume of liquefiable soil along with an assessment of confidence in an interpolation. The probability of liquefaction volumes compare well with those predicted using geologic interpretations.     

Seismicity and Seismic Hazard in Alexandria (Egypt) and its Surroundings

— Alexandria City has suffered great damage due to earthquakes from near and distant sources, both in historical and recent times. Sometimes the source of such damages is not well known. Seismogenic zones such as the Red Sea, Gulf of Aqaba-Dead Sea Hellenic Arc, Suez-Cairo-Alexandria, Eastern-Mediterranean-Cairo-Faiyoum and the Egyptian costal area are located in the vicinity of this city. The Egyptian coastal zone has the lowest seismicity, and therefore, its tectonic setting is not well known. The 1998 Egyptian costal zone earthquake is a moderate complex source. It is composed of two subevents separated by 4 sec. The first subevent initiated at a depth of 28 km and caused a rupture of strike (347°), dip (29°) and slip (125°). The second subevent occurred at a shallower depth (24 km) and has a relatively different focal parameter (strike 334°, dip 60° and slip 60°). The available focal mechanisms strongly support the manifestation of a complex stress regime from the Hellenic Arc into the Alexandria offshore area. In the present study a numerical modeling technique is applied to estimate quantitative seismic hazard in Alexandria. In terms of seismic hazard, both local and remote earthquakes have a tremendous affect on this city. A local earthquake with magnitude M_s = 6.7 at the offshore area gives peak ground acceleration up to 300 cm/sec². The total duration of shaking expected from such an earthquake is about three seconds. The Fourier amplitude spectra of the ground acceleration reveals that the maximum energy is carried by the low frequency (1–3 Hz), part of the seismic waves. The largest response spectra at Alexandria city is within this frequency band. The computed ground accelerations due to strong earthquakes in the Hellenic Arc, Red Sea and Gulf of Aqaba are very small (less than 10 cm/sec²) although with long duration (up to 3 minutes).     

Low-frequency transfer of seismic energy by superficial soil deposits and soft rocks

Recordings of recent strong earthquakes obtained on alluvial sites show that the maximum horizontal accelerations tend towards a limit of about 0.45 to 0.50g, associated with large displacements. by contrast, vertical accelerations do not appear to be subject to such a limit (1.7g for the 1979 Imperial Valley earthquake). Theoretical linear elasticity models, when applied to superficial layers of low strength, seem to be inadequate for the prediction of near-field ground motions in alluvial deposits. A good approximation for the horizontal component of certain Imperial Valley records was, however, obtained through a non-linear approach, using local soil properties together with a reasonable hypothesis for motion at the base of the superficial layers in question that included large values of acceleration for high-frequency shear waves.     

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.     

A spatial cross‐correlation model of spectral accelerations at multiple periods

Many seismic loss problems (such as disruption of distributed infrastructure and losses to portfolios of structures) are dependent upon the regional distribution of ground‐motion intensity, rather than intensity at only a single site. Quantifying ground‐motion over a spatially‐distributed region therefore requires information on the correlation between the ground‐motion intensities at different sites during a single event. The focus of the present study is to assess the spatial correlation between ground‐motion spectral accelerations at different periods. Ground motions from eight well‐recorded earthquakes were used to study the spatial correlations. On the basis of obtained empirical correlation estimates, we propose a geostatistics‐based method to formulate a predictive model that is suitable for simulation of spectral accelerations at multiple sites and multiple periods, in the case of crustal earthquakes in active seismic regions. While the calibration of this model and investigation of its implications were somewhat complex, the model itself is very simple to use for making correlation predictions. A user only needs to evaluate a simple equation relying on three sets of coefficients provided here to compute a correlation coefficient for spectral values at two periods and at a specified separation distance. These results may then be used in evaluating the seismic risk of portfolios of structures with differing fundamental periods. Copyright © 2012 John Wiley & Sons, Ltd.     

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.