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.     

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.     

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.     

Attenuation of seismic waves from moderate earthquakes of the Vrancea zone

The paper presents results of studying the attenuation of ground accelerations from earthquakes of the Vrancea with magnitudes less that 6.0 at distance greater than 300 km in a narrow sector located northeast of the focal region, within the limits of which are located acting and planned nuclear power plants (NNPs). Attenuation of peak ground accelerations in individual sections along the Vrancea–IRIS OBN station path is approximated by individual functions. It is shown that for a rough estimation of the seismic effect of earthquakes from the Vrancea zone it is acceptable to use the empirical relationship obtained by F.F. Aptikaev. For a more accurate estimate, it needs to be modified by adding a constant, whose value depends on the specific conditions of the NPP sites. It is shown that the results of data analysis on moderate earthquakes can be extrapolated to an earthquake with the maximum possible magnitude for the Vrancea zone and used to estimate the maximum seismic effects at the sites of operating and planned NPPs.     

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.     

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.     

Regularities in the time variations of seismic parameters and their implications for prediction of strong earthquakes in Greece

Time variations in the parameters of seismic activity in two regions in Greece, which are known to have different geodynamical conditions, are analyzed using the FastBEE algorithm suggested in (Papadopoulos and Baskoutas, 2009). The study is based on the data on weak earthquakes that occurred in two local regions. One region pertains to the zone dominated by intensive compression stress field, while another is located in the region of a relatively lower intensity extension stress field. It is shown that in the zone of compression the seismic parameters exhibit anomalous temporal behavior before strong earthquakes with Ms ≥ 5.7, whereas in the zones of extension, similar anomalies precede earthquakes with lower magnitudes of up to Ms ≥ 4.9. The most informative parameters for the purposes of predicting strong seismic events are the released seismic energy in the form logE ^2/3 and the slope of the frequency-magnitude dependence, b-value. The seismic activity in the region, expressed in terms of the logarithmic number of earthquakes, per unit time in some cases does not exhibit any particular pattern of behavior before strong earthquakes. In the time series of the studied parameters, four stages in the seismic process are clearly distinguished before strong earthquakes. Typically, a strong earthquake has a low probability to occur within the first two stages. Instead, this probability arises at stage III and attains its maximum at the end of this stage coinciding with the occurrence of the strong earthquake. We suggest these features of the time series to be used for the assessment of seismic hazard and for the real-time prediction of strong earthquakes. The time variations in the b-value are found to be correlated with the time variations inlogE ^2/3. This correlation is closely approximated by the power-law function. The parameters of this function depend on the geodynamical features of the region and characterize the intensity and the type of the regional tectonic stresses. The results of our study show that the FastBEE algorithm can be successfully applied for monitoring seismic hazard and predicting strong earthquakes.     

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.     

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.     

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.     

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.     

Scenario predictions for potential near-field and far-field earthquakes affecting Hong Kong

Historical earthquakes noted in the written records of the South China region, including Hong Kong, are not well delineated along identified prominent fault sources. Despite the lack of any definitive, localised trend in the spatial distribution of seismic activity in the region, there does appear to be some major disparity in the seismic activity rates (especially for large magnitude earthquakes) between the near-field and the far-field regions of Hong Kong. Despite this observation, previous studies of the regional seismic activity and seismic ground motion hazard (the latter using a probabilistic seismic hazard assessment, PSHA) have considered very broad source zone regions, in which uniform levels of seismic activity have been assumed. The present paper further scrutinises this broad source zone (BSZ) approach by adopting a novel expanding circular disc (ECD) method to determine the rates of earthquake recurrence. Such a method is intended to counter-check previously developed models by determining earthquake scenario events in terms of magnitude–distance (M–R) pairs or combinations, having defined values of average return period. Unlike the BSZ approach, the ECD method specifically accounts for the supposed variations in the seismic activity rates between events in the near-field and the far-field of Hong Kong.The form of the developed method is particularly suited to the determination of design-level earthquake ground motions for bedrock sites, since it assumes a directionally-independent attenuation model as described in the companion paper. It is found that, whilst the BSZ approach may indicate the overall average levels of hazard that are representative of the South China region as a whole, it does not capture the large disparity in seismic activity rates between near-field and far-field events. This important feature is expected to have a significant impact upon engineering assessments of the seismic safety of structures in Hong Kong and elsewhere in the South China region. For example, it is found that for events with M≥6, the seismic activity rate (normalised by time and area) in the very far-field is around 3.5 times larger than in the near-field and medium-field of Hong Kong. The resulting design M-R combinations, covering a range of return periods from 70 to 2500 years, are limited, for very long return periods and for distant events, by the maximum credible earthquake (MCE) magnitude. Intensive research to determine this seismic hazard parameter is recommended, in order to refine further the results of the ECD analysis, which presently conservatively assumes the MCE to range between M=6 in the near-field of Hong Kong to M=8 in the very far-field, at distances greater than 280 km from Hong Kong.     

Re-appraisal of seismicity and seismotectonics in the north-eastern Algeria Part I: Review of historical seismicity

As large destructive seismic events are not frequent in Algeria, anexhaustive knowledge of the historical seismicity is required to have arealistic view of seismic hazard in this part of the world. This research workpresents a critical reappraisal of seismicity in the north-eastern Algeria forseismotectonic and seismic hazard purposes. This part of work focuses onthe seismicity of pre-1900 period for the area under consideration[33^°N-38^°N, 4^°E-9.5^°E]. By going back tothe available documentary sources and evaluating and analysing the eventsin geographical, cultural and historical context, it has been possible toidentify 111 events, from 1850–1899, which are not reported in therecent Algerian catalogue. Several spurious events, reported in standardlistings, have been deleted and nine unknown events have been discovered.It is quite clear that macroseismic information derived from press reportsand published documents in Algeria, under certain conditions, is veryincomplete, even for destructive earthquakes, located in the countrysideaway from communication centres. One of the reasons for this iscensorship, noticeable during the colonisation period. Critical analysis ofnewly collected information has allowed the determination and/or theimprovement of the macroseismic parameters of each event, such aslocation, maximum epicentral intensity and magnitude to produce anearthquake catalogue as homogeneous and complete as the available data,for the zone under study. The criteria used in this research are explainedand eight historical earthquakes have been the subject of retrospectivemacroseismic field construction.The investigation of historical earthquakes is one of the most important taskin studying seismotectonic for seismic hazard evaluation purposes.     

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

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.     

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.     

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.     

Seismic Hazard Assessment of the East Thuringian Region/Germany – Case Study

East Thuringia/Germany, especially the region Gera-Ronneburg, is part of the large Kyffhäuser-Jachymov-Fault-Zone and displays moderate seismicity. However, its seismic hazard is significantly higher than that of the surrounding area including the Vogtland/Northern Bohemian region. The earthquake catalogue of Germany contains for this region besides the well-investigated Central German Earthquake (March 1872, I ₀ =VII-VIII) entries of up to I ₀ =VIII (14th century). Epicentral intensities and coordinates of these historical earthquakes are considered as uncertain. In seismic hazard analysis historical events which are uncertain are often neglected. But, especially in regions of moderate seismicity and infrequent larger earthquakes, the time window considered should be extended as far as possible. Apart from the necessity to study the historical sources of the strongest 14th century earthquakes, we investigate the influence of these events on the seismic hazard, taking into account the uncertainties of their size and location. Generally, the investigations clearly reveal the importance of defining source regions on the one hand and the significance of the local relevant attenuation function on the other hand. A further important point in seismic hazard assessment is the strong influence of the geological site conditions on seismic hazard (amplification or damping phenomena). For both points the well-known Central German Earthquake (1872) supplies important information.     

Response spectral attenuation relationships for Sumatran-subduction earthquakes and the seismic hazard implications to Singapore and Kuala Lumpur

Singapore and Kuala Lumpur, the capital of Malaysia, may well represent the classic examples of area with low seismic hazard but with high consequence. Both cities are located in a low-seismicity region of Southeast Asia, where active seismic sources are located more than 300 km away. Seismic designs have not been implemented in this seemingly low-hazard region though distant earthquakes in Sumatra had frequently shaken high-rise structures in the two cities. Several studies have been conducted to systematically assess the seismic hazards of Singapore and the Malay Peninsula. The present research particularly addresses issues in deriving a new set of attenuation relationships of peak ground acceleration (PGA), peak ground velocity (PGV) and response spectral acceleration (RSA) for the Sumatran-subduction earthquakes. To be relevant for the seismic hazard assessment of the remote metropolises, the derived attenuation relationships cover a long distance range from 150 to 1500 km. The attenuation relationships are derived using synthetic seismograms that account for source and path effects. The uncertainties in rupture parameters, such as stress drop, strike, dip and rake angles, have been defined according to the regional geological and tectonic settings as well as the ruptures of previous earthquakes. The seismic potential of the Sumatran subduction zone are high in the region from 2°N to 5°S as there has been no recurrence of great thrust events since 1861. A large event with M_w greater than 7.8 in this particular subduction zone may be capable of generating destructive ground motions in Singapore and Kuala Lumpur, even at a distance of 700 km.