Earthquake hazard and risk assessment based on Unified Scaling Law for Earthquakes: Greater Caucasus and Crimea

We continue applying the general concept of seismic risk analysis in a number of seismic regions worldwide by constructing regional seismic hazard maps based on morphostructural analysis, pattern recognition, and the Unified Scaling Law for Earthquakes (USLE), which generalizes the Gutenberg-Richter relationship making use of naturally fractal distribution of earthquake sources of different size in a seismic region. The USLE stands for an empirical relationship log₁₀N(M, L)?=?A?+?B·(5 – M)?+?C·log₁₀L, where N(M, L) is the expected annual number of earthquakes of a certain magnitude M within a seismically prone area of linear dimension L. We use parameters A, B, and C of USLE to estimate, first, the expected maximum magnitude in a time interval at seismically prone nodes of the morphostructural scheme of the region under study, then map the corresponding expected ground shaking parameters (e.g., peak ground acceleration, PGA, or macro-seismic intensity). After a rigorous verification against the available seismic evidences in the past (usually, the observed instrumental PGA or the historically reported macro-seismic intensity), such a seismic hazard map is used to generate maps of specific earthquake risks for population, cities, and infrastructures (e.g., those based on census of population, buildings inventory). The methodology of seismic hazard and risk assessment is illustrated by application to the territory of Greater Caucasus and Crimea.     

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.     

Calibration of local magnitude scale for Hindukush continental subduction zone

Hindukush is an active subduction zone where at least one earthquake occurs on daily basis. For seismic hazard studies, it is important to develop a local magnitude scale using the data of local seismic network. We have computed local magnitude scale for Hindukush earthquakes using data from local network belonging to Center for Earthquake Studies (CES) for a period of three years, i.e. 2015–2017. A total of 26,365 seismic records pertaining to 2,683 earthquakes with magnitude 2.0 and greater, was used with hypocentral distance less than 600 km. Magnitude scale developed by using this data comes to be M_L = logA + 0.929logr + 0.00298r – 1.84. The magnitude determined through formulated relation was compared with that of standard relation for Southern California and relation developed by the same authors for local network for Northern Punjab. It was observed that Hindukush region has high attenuation as compared to that of Southern California and Northern Punjab which implies that Hindukush is tectonically more disturbed as compared to the said regions, hence, seismically more active as well. We have calculated station correction factors for our network. Station correction factors do not show any pattern which probably owes to the geological and tectonic complexity of this structure. Standard deviation and variance of magnitude residuals for CES network determined using Hutton and Boore scale and scale developed in this study were compared, it showed that a variance reduction of 44.1% was achieved. Average of magnitude residuals for different distance ranges was almost zero which showed that our magnitude scale was stable for all distances up to 600 km. Newly developed magnitude scale will help in homogenization of earthquake catalog. It has been observed that b-value of CES catalog decreases when magnitude is calculated by using newly developed magnitude scale.     

The 2014 seismic hazard model of the Middle East: overview and results

The Earthquake Model of Middle East (EMME) Project aimed to develop regional scale seismic hazard and risk models uniformly throughout a region extending from the Eastern Mediterranean in the west to the Himalayas in the east and from the Gulf of Oman in the south to the Greater Caucasus in the North; a region which has been continuously devastated by large earthquakes throughout the history. The 2014 Seismic Hazard Model of Middle East (EMME-SHM14) was developed with the contribution of several institutions from ten countries. The present paper summarizes the efforts towards building a homogeneous seismic hazard model of the region and highlights some of the main results of this model. An important aim of the project was to transparently communicate the data and methods used and to obtain reproducible results. By doing so, the use of the model and results will be accessible by a wide community, further support the mitigation of seismic risks in the region and facilitate future improvements to the seismic hazard model. To this end all data, results and methods used are made available through the web-portal of the European Facilities for Earthquake Hazard and Risk (www.efehr.org).     

Assessment of seismic hazard in the North Caucasus

The seismicity of the North Caucasus is the highest in the European part of Russia. The detection of potential seismic sources here and long-term prediction of earthquakes are extremely important for the assessment of seismic hazard and seismic risk in this densely populated and industrially developed region of the country. The seismogenic structures of the Iran-Caucasus-Anatolia and Central Asia regions, adjacent to European Russia, are the subjects of this study. These structures are responsible for the specific features of regional seismicity and for the geodynamic interaction with adjacent areas of the Scythian and Turan platforms. The most probable potential sources of earthquakes with magnitudes M = 7.0 ± 0.2 and 7.5 ± 0.2 in the North Caucasus are located. The possible macroseismic effect of one of them is assessed.     

Seismogeodynamics of lineament structures in the mountainous regions bordering the Scythian-Turan plate

The Scythian-Turan platform, together with the Alpine Iran-Caucasus-Anatolia and Hercynian Central Tien Shan orogenic structures adjacent to it, represents a coherent seismogeodynamic system responsible for regional seismicity features in the territory under consideration. Investigations of the spatiotemporal and energy evolution of seismogeodynamic processes along the main lineament structures of the orogen reveal characteristic features directly related to the prediction of seismic hazard in this region, as well as in southern European Russia. These characteristics primarily include kinematic features in the sequences of seismic events of various magnitudes and an ordered migration of seismic activation, enabling the more or less reliable determination of the occurrence time intervals (years) and areas of forthcoming large earthquakes (magnitudes of 7.0 ± 0.2, 7.5 ± 0.2, and 8.0 ± 0.2).     

K-means cluster analysis and seismicity partitioning for Pakistan

Pakistan and the western Himalaya is a region of high seismic activity located at the triple junction between the Arabian, Eurasian and Indian plates. Four devastating earthquakes have resulted in significant numbers of fatalities in Pakistan and the surrounding region in the past century (Quetta, 1935; Makran, 1945; Pattan, 1974 and the recent 2005 Kashmir earthquake). It is therefore necessary to develop an understanding of the spatial distribution of seismicity and the potential seismogenic sources across the region. This forms an important basis for the calculation of seismic hazard; a crucial input in seismic design codes needed to begin to effectively mitigate the high earthquake risk in Pakistan. The development of seismogenic source zones for seismic hazard analysis is driven by both geological and seismotectonic inputs. Despite the many developments in seismic hazard in recent decades, the manner in which seismotectonic information feeds the definition of the seismic source can, in many parts of the world including Pakistan and the surrounding regions, remain a subjective process driven primarily by expert judgment. Whilst much research is ongoing to map and characterise active faults in Pakistan, knowledge of the seismogenic properties of the active faults is still incomplete in much of the region. Consequently, seismicity, both historical and instrumental, remains a primary guide to the seismogenic sources of Pakistan. This study utilises a cluster analysis approach for the purposes of identifying spatial differences in seismicity, which can be utilised to form a basis for delineating seismogenic source regions. An effort is made to examine seismicity partitioning for Pakistan with respect to earthquake database, seismic cluster analysis and seismic partitions in a seismic hazard context. A magnitude homogenous earthquake catalogue has been compiled using various available earthquake data. The earthquake catalogue covers a time span from 1930 to 2007 and an area from 23.00° to 39.00°N and 59.00° to 80.00°E. A threshold magnitude of 5.2 is considered for K-means cluster analysis. The current study uses the traditional metrics of cluster quality, in addition to a seismic hazard contextual metric to attempt to constrain the preferred number of clusters found in the data. The spatial distribution of earthquakes from the catalogue was used to define the seismic clusters for Pakistan, which can be used further in the process of defining seismogenic sources and corresponding earthquake recurrence models for estimates of seismic hazard and risk in Pakistan. Consideration of the different approaches to cluster validation in a seismic hazard context suggests that Pakistan may be divided into K?=?19 seismic clusters, including some portions of the neighbouring countries of Afghanistan, Tajikistan and India.     

A time probabilistic approach to seismic landslide hazard estimates in Iran

Understanding where seismically induced landslides are most likely to occur is crucial in land use planning and civil protection actions aimed at reducing property damage and loss of life in future earthquakes. For this purpose an approach proposed by Del Gaudio et al. [1] has been applied to the whole Iranian territory to provide the basis to assess location and temporal recurrence of conditions of seismic activation of slope failures, according to the Newmark's model [2]. Following this approach, occurrence probabilities for different levels of seismic shaking in a time interval of interest (50 years) were first obtained through a standard hazard estimate procedure. Then, empirical formulae in the form proposed by Jibson et al. [3] and calibrated for the main seismogenic Iranian regions were used to evaluate the slope critical acceleration (Ac)_x for which a prefixed probability exists that, under seismic shakings, Newmark's displacement D_N exceeds a threshold×corresponding to landslide triggering conditions. The obtained (Ac)_x values represent the minimum slope resistance required to limit the probability of landslide seismic triggering within the prefixed value. A map reporting the spatial distribution of these values gives comparative indications on regional different exposure of slopes to shaking capable of inducing failures and provides a reference for hazard estimate at local scale. The obtained results show that the exposure to landslide seismic induction is maximum in the Alborz Mountains region, where critical accelerations up to ∼0.1 g are required to limit the probability of seismic triggering of coherent type landslides within 10% in 50 years.     

A new seismic hazard analysis using FOSM algorithms

From recent lessons, it is evident that earthquake prediction is immature and impractical as of now. Under the circumstances, seismic hazard analysis is considered a more practical approach for earthquake hazard mitigation, by estimating the annual rate of earthquake ground motions (or seismic hazard) based on seismicity and other geological evidences. Like other earthquake studies for the high-seismicity region around Taiwan, this study aims to conduct a new seismic hazard assessment for the region using the well-established FOSM (first-order second-moment) algorithm, on the record of 55,000 earthquakes observed in the past 110 years. The new seismic hazard analysis from a different perspective shows that the annual rate for earthquake-induced PGA to exceed the current design value (i.e., 0.23g) in two major cities in Taiwan should be relatively low, with it no greater than 0.0006 per year. Besides, the FOSM estimates were found very close to those with Monte Carlo Simulation (MCS), mainly because the skewness of the three random variables (i.e., earthquake magnitude, location, and model error) considered in the probabilistic analysis is not very large.     

The correlation between the characteristics of seismic wave propagation in Western Caucasus and the geological–tectonic features of the region

The relationship between the characteristics of seismic waves in the Western Caucasus and the geological-tectonic structure of the region is studied for identifying the specificity of seismic propagation in the mountainous regions with a complicated geological structure and forecasting the characteristics of the propagation from the geological and tectonic data. The interpretation is presented for the estimates of the Q-factor of the medium (Q(f) ~ 55f^0.9 in the region of Sochi and Q(f) ~ 90f^0.7 in the region of Anapa), seismic wave enhancement in the upper crustal layers (A(f) ~ 1), and peak ground acceleration residuals, which were previously determined from the records of the local earthquakes and show the distributions of local variations in the parameters of seismic wave radiation and propagation. The obtained characteristics are interpreted in the context of the up-to-date information about the tectonic, geological, and deep structure of the epicentral zones in the Western Caucasus and neighboring territory of the Black Sea. The discrepancies revealed in the low-frequency behavior of the Q-factor in the vicinities of Sochi and Anapa is accounted for by the spatial scale and character of tectonic dislocations of the rocks in these regions. The local variations in the parameters of seismic radiation and propagation are probably related to the geological features of the region such as the fault structures, including the thrusts, shatter zones, oblique seismic boundaries, variations in the thickness and consolidation of the sedimentary cover, as well as the peculiarities in the structure and material composition of the basement.     

Some Probabilistic and Statistical Properties of the Seismic Regime of Zemmouri (Algeria) Seismoactive Zone

Statistical tests have been used to adjust the Zemmouri seismic data using a distribution function. The Pareto law has been used and the probabilities of various expected earthquakes were computed. A mathematical expression giving the quantiles was established. The extreme values limiting law confirmed the accuracy of the adjustment method. Using the moment magnitude scale, a probabilistic model was made to predict the occurrences of strong earthquakes. The seismic structure has been characterized by the slope of the recurrence plot γ, fractal dimension D, concentration parameter K_sr, Hurst exponents H_r and H_t. The values of D, γ, K_sr, H_r, and H_t diminished many months before the principal seismic shock (M = 6.9) of the studied seismoactive zone has occurred. Three stages of the deformation of the geophysical medium are manifested in the variation of the coefficient G% of the clustering of minor seismic events.     

The active fault system of SW Alps

Historical and active seismicity in the south-western Alps (France and Italy) shows the recurrence of relatively high-magnitude earthquakes (M ≥ 5.8), like the one that recently affected the Italian Apennine range (M = 6.3 on the 30th March 2009). However, up-to-date detailed mapping of the active fault network has been poorly established. The evaluation of seismological hazard in particular in the highly populated French and Italian coastal region cannot be done without this. Here, we present a detailed study of the main active fault system, based on geological observations along the south-western flank of the Alpine arc. This N140° right-lateral strike-slip active fault system runs along the edge of the Argentera-Mercantour range and can be followed down to the Mediterranean Sea. It is evidenced by (1) Holocene offsets of glacial geomorphology witnessing ongoing fault activity since 10 ka, (2) widespread recent (10–20 Ma) pseudotachylytes featuring long term activity of the faults, (3) active landslides along the main fault zone, (4) geothermal anomalies (hot springs) emerging in the active faults, (5) ongoing low-magnitude seismic activity and (6) localization of the main historical events. In the light of our investigations, we propose a new tectonic pattern for the active fault system in the south-western Alps.     

Long-term earthquake triggering in the Southern and Northern Apennines

We argue that the study of long-range interaction between seismic sources in the peri-Adriatic regions may significantly contribute to estimating seismic hazard in Italy. This hypothesis is supported by the reconstruction of the geodynamic and tectonic settings in the Central Mediterranean region, the space–time distribution of major past earthquakes, and the quantification of post-seismic relaxation. The most significant evidence of long-distance interaction is recognized for the Southern Apennines, whose major earthquakes have almost regularly followed within a few years the largest events in the Montenegro-Albania zone since 1850. Statistical analyses of the post-1850 earthquake catalogues give a probability of about 10% that a major event in the Southern Apennines is not preceded by the occurrence of a strong event in the Southern Dinarides–Albanides within 3–5 years. Conversely, the probability of false alarms is relevant (50% within 3 years, 33% within 5 years). Northward, the tectonic setting and some patterns of regularity seen in major events suggest that the seismic activation of the main transtensional decoupling shear zones in the Central Apennines should influence the probability of major earthquakes in the Northern Apennines.     

Seismic hazard assessment in the megacity of Blida (Algeria) and its surrounding regions using parametric-historic procedure

The region of Blida is characterized by a relatively high seismic activity, pointed especially during the past two centuries. Indeed, it experienced a significant number of destructive earthquakes such as the earthquakes of March 2, 1825 and January 2, 1867, with intensity of X and IX, respectively. This study aims to investigate potential seismic hazard in Blida city and its surrounding regions. For this purpose, a typical seismic catalog was compiled using historical macroseismic events that occurred over a period of a few hundred years, and the recent instrumental seismicity dating back to 1900. The parametric-historic procedure introduced by Kijko and Graham (1998, 1999) was applied to assess seismic hazard in the study region. It is adapted to deal with incomplete catalogs and does not use any subjective delineation of active seismic zones. Because of the lack of recorded strong motion data, three ground prediction models have been considered, as they seem the most adapted to the seismicity of the study region. Results are presented as peak ground acceleration (PGA) seismic hazard maps, showing expected peak accelerations with 10% probability of exceedance in 50-year period. As the most significant result, hot spot regions with high PGA values are mapped. For example, a PGA of 0.44 g has been found in a small geographical area centered on Blida city.     

空间光滑地震活动性模型中光滑函数的比较研究

使用Frankel提出的基于空间光滑地震活动性模型的地震危险性分析方法,选择华南、华北、川滇3个地区的地震记录,比较分析了高斯、幂律和地震分形分布光滑函数3种光滑函数在不同地区的适用性．结果表明,使用交叉验证法可以为高斯光滑函数选取合适的相关距离c值,光滑得到的地震活动性模型能够真实反映研究区域的地震活动特征,根据活动性模型计算得出的峰值加速度(PGA)分布也符合人们对研究区域地震危险性的认识．幂律光滑函数适用于地震活动性较强的地区,且具有容易求取光滑参数的优点．光滑程度较低的幂律光滑函数不适用于地震活动性弱的地区,在该类地区应选择光滑程度较高的高斯光滑函数．地震分形分布光滑函数不适用于地震活动较强且地震活动强度差异较大的地区,其容易过分高估高震级地震对地震危险性的影响,而忽略了低震级地震对地震危险性的贡献．但对于地震活动较弱且地震活动强度差异较小的地区,可使用地震分形分布光滑函数,且同样具有容易求取光滑参数的优点．      

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.     

Continuous evaluation of seismic hazard induced by the deposit extraction in selected coal mines in Poland

A probabilistic relation between seismic activity and the volumeV of extracted deposits in mines is derived $$\Sigma E = C \cdot V^B ,$$ whereC andB are parameters characterizing mining works and the state of rock mass. Assuming that the measure of seismic hazard is the amount of seismic energy released in a given time interval, it is shown how the hazard can be evaluated continuously. The derived relations were tested in selected coal mines in Upper Silesia.     

Seismic response in archaeological areas: the case-histories of Rome

Rome is affected by earthquakes associated to three different seismogenic districts: the Central Apennines area, the Colli Albani volcanic area and the Roman area. The major effects were exclusively due to Apennine seismicity and reached in some cases felt intensities up to VII–VIII degree (MCS scale). The predominant role in the damage distribution seems to be played by the local geological conditions. The historical centre of the city is characterized by the presence of two geomorphologic domains: the alluvial plain of Tiber river and the topographic relieves of Roman Hills, where tradition indicates the first site of the city foundation. In particular, the right river side is characterized by the outcropping of the regional bedrock along the Monte Mario–Gianicolo ridge, while the eastern relieves are the remnants of the Sabatini and Albani volcanic plateau, deeply eroded by the Tiber river and its tributaries during the last glacial low-stand (Würm). These domains are characterized by a large difference in seismic response, due to the high impedance contrast between Holocene coarse deposits filling the Tiber Valley and sedimentary and volcanic Plio–Pleistocene units. Seismic damage observed in 150 monuments of downtown Rome was indicating a significant concentration on alluvial recent deposits. This result was confirmed by the geographical distribution of conservation and retrofitting activities subsequent to main earthquakes, mostly related to local geological conditions. The cases of Marcus Aurelius' Column and Colosseum confirmed the influence of the Holocene alluvial network in local seismic response. During 2500 years of history, the monuments of Rome have `memorized' the seismic effects of historical earthquakes. In some cases, the integration of historical and geological research and macroseismic observations may provide original and useful indications to seismologists to define the seismic response of the city. Local site effects represent a serious threat for historical buildings in Rome and in other historical towns with similar cultural heritage and geological characteristics, as in the Mediterranean region, even in areas that are not affected by a local seismic activity.     

Deep-seated faults and lineaments, and the location of earthquake epicenters in the Russian northeast on land

This paper summarizes the available geological and geophysical material for faults as regards their role in the seismic process. The entirety of the geological and geophysical evidence is used to reveal hidden faults, which are important in influencing the spatial distribution of earthquakes, and to produce a map of the major earthquake-generating faults and lineaments in the Russian northeast. As well as the occurrence of earthquakes at known faults that have surface expression, we find that seismicity tends to occur at the hidden faults and lineaments we have identified, as well as at intersections of faults. We made a quantitative assessment of the relationship of seismicity to tectonic fragmentation of the crust, correlating the density and discordance measure for faults to indicators of seismic activity (rate and energy release of earthquakes per unit area) for the southeast flank of the Okhotsk-Lena seismic region. The results obtained in this study revealed some features in the spatial distribution of earthquakes occurring on land in the Okhotsk-Lena seismic region: the maximum level of seismic activity occurs in areas with moderate values of the discordance measure for faults (12 < ‖D‖ ≤ 18) as identified from gravity data and in zones of increased horizontal gradients of the lines of equal discordance. At these locations, the greatest probability of earthquake occurrence for events of energy class K ≥ 12 corresponds to moderate values of the density of faults visible at the surface (0.12 < τ ≤ 0.16 km^?1).     

Probabilistic seismic hazard assessment for the two layer fault system of Antalya (SW Turkey) area

Southwest Turkey, along Mediterranean coast, is prone to large earthquakes resulting from subduction of the African plate under the Eurasian plate and shallow crustal faults. Maximum observed magnitude of subduction earthquakes is Mw = 6.5 whereas that of crustal earthquakes is Mw = 6.6. Crustal earthquakes are sourced from faults which are related with Isparta Angle and Cyprus Arc tectonic structures. The primary goal of this study is to assess seismic hazard for Antalya area (SW Turkey) using a probabilistic approach. A new earthquake catalog for Antalya area, with unified moment magnitude scale, was prepared in the scope of the study. Seismicity of the area has been evaluated by the Gutenberg-Richter recurrence relationship. For hazard computation, CRISIS2007 software was used following the standard Cornell-McGuire methodology. Attenuation model developed by Youngs et al. Seismol Res Lett 68(1):58–73, (1997) was used for deep subduction earthquakes and Chiou and Youngs Earthq Spectra 24(1):173–215, (2008) model was used for shallow crustal earthquakes. A seismic hazard map was developed for peak ground acceleration and for rock ground with a hazard level of a 10% probability of exceedance in 50 years. Results of the study show that peak ground acceleration values on bedrock change between 0.215 and 0.23 g in the center of Antalya.