Ground-motion characterization for the probabilistic seismic hazard assessment in Turkey

This study describes the methodology implemented to establish the ground-motion logic-tree for national probabilistic seismic hazard map of Turkey for shallow active crustal regions. The presented procedure provides quantitative information to guide the hazard experts while establishing the logic tree to capture the epistemic uncertainty in ground-motion characterization. It uses non-data-driven and data-driven testing methods to identify and rank candidate ground-motion prediction equations (GMPEs) under a specific ground-motion database. The candidate GMPEs are subjected to visual inspection and are classified into center, body and range (CBR) spectral estimates for a proper consideration of epistemic uncertainty. The GMPEs classified into CBR are then used in a suite of seismic hazard sensitivity analysis to establish the most suitable GMPE logic-tree whose spectral estimates are not biased by any one of the GMPEs in the logic-tree structure. The sensitivity analysis considers normalized spectral ordinates and is not manipulated by the spectral amplitudes. The proposed procedure is inherited from the relevant studies of the Earthquake Model of the Middle East (EMME; www.efehr.org:8080/jetspeed/portal/emme.psml) regional seismic hazard project. This paper also highlights the similarities and differences in ground-motion characterization between EMME and our approach.     

The 2014 Earthquake Model of the Middle East: ground motion model and uncertainties

We summarize the main elements of a ground-motion model, as built in three-year effort within the Earthquake Model of the Middle East (EMME) project. Together with the earthquake source, the ground-motion models are used for a probabilistic seismic hazard assessment (PSHA) of a region covering eleven countries: Afghanistan, Armenia, Azerbaijan, Cyprus, Georgia, Iran, Jordan, Lebanon, Pakistan, Syria and Turkey. Given the wide variety of ground-motion predictive models, selecting the appropriate ones for modeling the intrinsic epistemic uncertainty can be challenging. In this respect, we provide a strategy for ground-motion model selection based on data-driven testing and sensitivity analysis. Our testing procedure highlights the models of good performance in terms of both data-driven and non-data-driven testing criteria. The former aims at measuring the match between the ground-motion data and the prediction of each model, whereas the latter aims at identification of discrepancies between the models. The selected set of ground models were directly used in the sensitivity analyses that eventually led to decisions on the final logic tree structure. The strategy described in great details hereafter was successfully applied to shallow active crustal regions, and the final logic tree consists of four models (Akkar and Ça?nan in Bull Seismol Soc Am 100:2978–2995, 2010; Akkar et al. in Bull Earthquake Eng 12(1):359–387, 2014; Chiou and Youngs in Earthq Spectra 24:173–215, 2008; Zhao et al. in Bull Seismol Soc Am 96:898–913, 2006). For other tectonic provinces in the considered region (i.e., subduction), we adopted the predictive models selected within the 2013 Euro-Mediterranean Seismic Hazard Model (Woessner et al. in Bull Earthq Eng 13(12):3553–3596, 2015). Finally, we believe that the framework of selecting and building a regional ground-motion model represents a step forward in ground-motion modeling, particularly for large-scale PSHA models.     

New Seismic Noise Models Obtained Using Very Broadband Stations

It has been two decades since the last comprehensive standard model of ambient earth noise was published Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The PETERSON model was updated by analyzing the absolute quietest conditions for stations within the GSN (Berger et al. in J Geophys Res 109, 2005; Mcnamara and Buland in Bull Seism Soc Am 94:1517–1527, 2004; Ringler et al. in Seismol Res Lett 81(4) doi:10.1785/gssrl.81.4.605, 2010). Unfortunately, both the original model and the updated models did not include any deployed station in North Africa and Middle East, which reflects the noise levels within the desert environment of those regions. In this study, a survey was conducted to create a new seismic noise model from very broadband stations which recently deployed in North Africa. For this purpose, 1 year of continuous recording of seismic noise data of the Egyptian National Seismic Network (ENSN) was analyzed in order to create a new noise model. Seasonal and diurnal variations in noise spectra were recorded at each station. Moreover, we constructed a new noise model for each individual station. Finally, we obtained a new cumulative noise model for all the stations. We compared the new high-noise model (EHNM) and new low-noise model (ELNM) with both the high-noise model (NHNM) and low-noise model (NLNM) of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The obtained noise levels are considerably lower than low-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993) at ultra long period band (ULP band), but they are still below the high-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The results of this study could be considered as a first step to create permanent seismic noise models for North Africa and Middle East regions.     

Probabilistic seismic hazard assessment of Eastern Marmara Region

The objective of this study is to evaluate the seismic hazard in Eastern Marmara Region using an improved probabilistic seismic hazard assessment methodology. Two significant improvements over the previous seismic hazard assessment practices are accomplished in this study: advanced seismic source characterization models in terms of source geometry and recurrence relationships are developed, and improved global ground motion models (NGA-W1 models) are employed to represent the ground motion variability. Planar fault segments are defined and a composite magnitude distribution model is used for all seismic sources in the region to properly represent the characteristic behavior of the North Anatolian Fault without the need for an additional background zone. Multi-segment ruptures are considered using the rupture model proposed by the Working Group on California Earthquake Probabilities (2003). Events in the earthquake catalogue are attributed to the fault zones and scenario weights are determined by releasing the accumulated seismic energy. The uniform hazard spectra at 10 % probability of exceedance in 50 years hazard level for different soil conditions (soil and rock) are revealed for specific locations in the region (Adapazar?, Düzce, Gemlik, Izmit, Iznik and Sapanca). Hazard maps of the region for rock site conditions at the selected hazard levels are provided to allow the readers perform site-specific hazard assessment and develop site-specific design spectrum for local site conditions.     

Seismic vulnerability scenarios of Unreinforced Masonry churches in New Zealand

A vulnerability analysis of c.300 unreinforced Masonry churches in New Zealand is presented. The analysis uses a recently developed vulnerability index method (Cattari et al. in Proceedings of the New Zealand Society for Earthquake Engineering NZSEE 2015 conference, Rotorua, New Zealand, 2015a; b; SECED 2015 conference: earthquake risk and engineering towards a Resilient World, Cambridge; Goded et al. in Vulnerability analysis of unreinforced masonry churches (EQC 14/660)—final report, 2016; Lagomarsino et al. in Bull Earthq Eng, 2018), specifically designed for New Zealand churches, based on a widely tested approach for European historical buildings. It consists of a macroseismic approach where the seismic hazard is defined by the intensity and correlated to post seismic damage. The many differences in typologies of New Zealand and European churches, with very simple architectural designs and a majority of one nave churches in New Zealand, justified the need to develop a method specifically created for this country. A statistical analysis of the churches damaged during the 2010–2011 Canterbury earthquake sequence was previously carried out to develop the vulnerability index modifiers for New Zealand churches. This new method has been applied to generate seismic scenarios for each church, based on the most likely seismic event for 500 years return period, using the latest version of New Zealand’s National Seismic Hazard Model. Results show that highly vulnerable churches (e.g. stone churches and/or with a weak structural design) tend to produce higher expected damage even if the intensity level is lower than for less vulnerable churches in areas with slightly higher seismicity. The results of this paper provide a preliminary tool to identify buildings requiring in depth structural analyses. This paper is considered as a first step towards a vulnerability analysis of all the historical buildings in the country, in order to preserve New Zealand’s cultural and historical heritage.     

Seismic hazard map of the Middle East

The collaborative project Earthquake Model of the Middle East (EMME, 2010–2015) brought together scientists and engineers from the leading research institutions in the region and delivered state-of-the-art seismic hazard assessment covering Afghanistan, Armenia, Azerbaijan, Cyprus, Georgia, Iran, Iraq, Jordan, Lebanon, Palestine, Pakistan, Syria and Turkey. Their efforts have been materialized in the first homogenized seismic hazard model comprising earthquake catalogues, mapped active faults, strong motions databank, ground motion models and the estimated ground motion values for various intensity measure types and relevant return periods (e.g. 475–5000 years). The reference seismic hazard map of the Middle East, depicts the mean values of peak ground acceleration with a 10% chance of exceedance in 50 years, corresponding to a mean return period of 475 years. A full resolution poster is provided with this contribution.     

An Evaluation of Earthquake Hazard Potential for Different Regions in Western Anatolia Using the Historical and Instrumental Earthquake Data

We applied the maximum likelihood method produced by Kijko and Sellevoll (Bull Seismol Soc Am 79:645–654, 1989; Bull Seismol Soc Am 82:120–134, 1992) to study the spatial distributions of seismicity and earthquake hazard parameters for the different regions in western Anatolia (WA). Since the historical earthquake data are very important for examining regional earthquake hazard parameters, a procedure that allows the use of either historical or instrumental data, or even a combination of the two has been applied in this study. By using this method, we estimated the earthquake hazard parameters, which include the maximum regional magnitude $ \hat{M}_{\max } , $ the activity rate of seismic events and the well-known $ \hat{b} $ value, which is the slope of the frequency-magnitude Gutenberg-Richter relationship. The whole examined area is divided into 15 different seismic regions based on their tectonic and seismotectonic regimes. The probabilities, return periods of earthquakes with a magnitude M?≥?m and the relative earthquake hazard level (defined as the index K) are also evaluated for each seismic region. Each of the computed earthquake hazard parameters is mapped on the different seismic regions to represent regional variation of these parameters. Furthermore, the investigated regions are classified into different seismic hazard level groups considering the K index. According to these maps and the classification of seismic hazard, the most seismically active regions in WA are 1, 8, 10 and 12 related to the Alia?a Fault and the Büyük Menderes Graben, Aegean Arc and Aegean Islands.     

A probabilistic seismic hazard assessment for the Turkish territory—part I: the area source model

The seismic zoning map of Turkey that is used in connection with the national seismic design code (versions issued both in 1997 and 2007) is based on a probabilistic seismic hazard assessment study conducted more than 20 years ago (Gülkan et al. in En son verilere göre haz?rlanan Türkiye deprem bölgeleri haritas?, Report No: METU/EERC 93-1, 1993). In line with the efforts for the update of the seismic design code, the need aroused for an updated seismic hazard map, incorporating recent data and state-of-the-art methodologies and providing ground motion parameters required for the construction of the design spectra stipulated by the new Turkish Earthquake Design Code. Supported by AFAD (Disaster and Emergency Management Authority of Turkey), a project has been conducted for the country scale assessment of the seismic hazard by probabilistic methods. The present paper describes the probabilistic seismic hazard assessment study conducted in connection with this project, incorporating in an area source model, all recently compiled data on seismicity and active faulting, and using a set of recently developed ground motion prediction equations, for both active shallow crustal and subduction regimes, evaluated as adequately representing the ground motion characteristics in the region. The area sources delineated in the model are fully parameterized in terms of maximum magnitude, depth distribution, predominant strike and dip angles and mechanism of possible ruptures. Resulting ground motion distributions are quantified and presented for PGA and 5 % damped spectral accelerations at T = 0.2 and 1.0 s, associated with return periods of 475 and 2475 years. The full set of seismic hazard curves was also made available for the hazard computation sites. The second part of the study, which is based on a fault source and smoothed seismicity model is covered in Demircioglu et al. in Bull Earthq Eng, (2016).     

Empirical ground-motion models for point- and extended-source crustal earthquake scenarios in Europe and the Middle East

This article presents the latest generation of ground-motion models for the prediction of elastic response (pseudo-) spectral accelerations, as well as peak ground acceleration and velocity, derived using pan-European databases. The models present a number of novelties with respect to previous generations of models (Ambraseys et al. in Earthq Eng Struct Dyn 25:371–400, 1996, Bull Earthq Eng 3:1–53, 2005; Bommer et al. in Bull Earthq Eng 1:171–203, 2003; Akkar and Bommer in Seismol Res Lett 81:195–206, 2010), namely: inclusion of a nonlinear site amplification function that is a function of $\text{ V }_\mathrm{S30}$ and reference peak ground acceleration on rock; extension of the magnitude range of applicability of the model down to $\text{ M }_\mathrm{w}$ 4; extension of the distance range of applicability out to 200 km; extension to shorter and longer periods (down to 0.01 s and up to 4 s); and consistent models for both point-source (epicentral, $\text{ R }_\mathrm{epi}$ , and hypocentral distance, $\text{ R }_\mathrm{hyp}$ ) and finite-fault (distance to the surface projection of the rupture, $\text{ R }_\mathrm{JB}$ ) distance metrics. In addition, data from more than 1.5 times as many earthquakes, compared to previous pan-European models, have been used, leading to regressions based on approximately twice as many records in total. The metadata of these records have been carefully compiled and reappraised in recent European projects. These improvements lead to more robust ground-motion prediction equations than have previously been published for shallow (focal depths less than 30 km) crustal earthquakes in Europe and the Middle East. We conclude with suggestions for the application of the equations to seismic hazard assessments in Europe and the Middle East within a logic-tree framework to capture epistemic uncertainty.     

Compatible ground-motion prediction equations for damping scaling factors and vertical-to-horizontal spectral amplitude ratios for the broader Europe region

In a companion article Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4, 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b) present a new ground-motion prediction equation (GMPE) for estimating 5 %-damped horizontal pseudo-acceleration spectral (PSA) ordinates for shallow active crustal regions in Europe and the Middle East. This study provides a supplementary viscous damping model to modify 5 %-damped horizontal spectral ordinates of Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b) for damping ratios ranging from 1 to 50 %. The paper also presents another damping model for scaling 5 %-damped vertical spectral ordinates that can be estimated from the vertical-to-horizontal (V/H) spectral ratio GMPE that is also developed within the context of this study. For consistency in engineering applications, the horizontal and vertical damping models cover the same damping ratios as noted above. The article concludes by introducing period-dependent correlation coefficients to compute horizontal and vertical conditional mean spectra (Baker in J Struct Eng 137:322–331, 2011). The applicability range of the presented models is the same as of the horizontal GMPE proposed by Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b): as for spectral periods $0.01 \hbox { s}\le \,\hbox {T}\le \,4\hbox { s}$ as well as PGA and PGV for V/H model; and in terms of seismological estimator parameters $4\le \hbox {M}_\mathrm{w} \le 8, \hbox { R} \le 200 \hbox { km}, 150\hbox { m/s}\le \hbox { V}_\mathrm{S30}\le $ 1,200 m/s, for reverse, normal and strike-slip faults. The source-to-site distance measures that can be used in the computations are epicentral $(\hbox {R}_\mathrm{epi})$ , hypocentral $(\hbox {R}_\mathrm{hyp})$ and Joyner–Boore $(\hbox {R}_\mathrm{JB})$ distances. The implementation of the proposed GMPEs will facilitate site-specific adjustments of the spectral amplitudes predicted from probabilistic seismic hazard assessment in Europe and the Middle East region. They can also help expressing the site-specific design ground motion in several formats. The consistency of the proposed models together with the Akkar et al. (Bull Earthq Eng, doi:10.1007/s10518-013-9461-4 2013a; Bull Earthq Eng, doi:10.1007/s10518-013-9508-6, 2013b) GMPE may be advantageous for future modifications in the ground-motion definition in Eurocode 8 (CEN in Eurocode 8, Design of structures for earthquake resistance—part 1: general rules, seismic actions and rules for buildings. European Standard NF EN 1998-1, Brussels, 2004).     

Recent evolution and challenges in the Seismic Hazard Analysis of the Po Plain region, Northern Italy

SH evaluation in Po Plain region is revisited with a view to: ( $a)$ updating available seismic source and ground motion attenuation models with recent data and testing the sensitivity of the results, ( $b)$ exploiting available non-stationary models to check their predictive performance at particular sites and, ( $c)$ handling the ground motion hazard generated by composite fault systems. Update of area source models and GMPEs does not strongly affect SH evaluated at representative sites, while usefulness of predictions derived from non-stationary models needs to be assessed on a site-by-site basis. A hybrid (deterministic-probabilistic) approach seems to be a promising tool for assessing site-specific hazard contributed by composite sources.     

Empirical equations for the prediction of PGA and pseudo spectral accelerations using Iranian strong-motion data

A recently compiled, comprehensive, and good-quality strong-motion database of the Iranian earthquakes has been used to develop local empirical equations for the prediction of peak ground acceleration (PGA) and 5%-damped pseudo-spectral accelerations (PSA) up to 4.0 s. The equations account for style of faulting and four site classes and use the horizontal distance from the surface projection of the rupture plane as a distance measure. The model predicts the geometric mean of horizontal components and the vertical-to-horizontal ratio. A total of 1551 free-field acceleration time histories recorded at distances of up to 200 km from 200 shallow earthquakes (depth < 30 km) with moment magnitudes ranging from M_w 4.0 to 7.3 are used to perform regression analysis using the random effects algorithm of Abrahamson and Youngs (Bull Seism Soc Am 82:505–510, 1992), which considers between-events as well as within-events errors. Due to the limited data used in the development of previous Iranian ground motion prediction equations (GMPEs) and strong trade-offs between different terms of GMPEs, it is likely that the previously determined models might have less precision on their coefficients in comparison to the current study. The richer database of the current study allows improving on prior works by considering additional variables that could not previously be adequately constrained. Here, a functional form used by Boore and Atkinson (Earthquake Spect 24:99–138, 2008) and Bindi et al. (Bull Seism Soc Am 9:1899–1920, 2011) has been adopted that allows accounting for the saturation of ground motions at close distances. A regression has been also performed for the V/H in order to retrieve vertical components by scaling horizontal spectra. In order to take into account epistemic uncertainty, the new model can be used along with other appropriate GMPEs through a logic tree framework for seismic hazard assessment in Iran and Middle East region.     

Real-time risk assessment in seismic early warning and rapid response: a feasibility study in Bishkek (Kyrgyzstan)

Earthquake early warning systems (EEWS) are considered to be an effective, pragmatic, and viable tool for seismic risk reduction in cities. While standard EEWS approaches focus on the real-time estimation of an earthquake’s location and magnitude, innovative developments in EEWS include the capacity for the rapid assessment of damage. Clearly, for all public authorities that are engaged in coordinating emergency activities during and soon after earthquakes, real-time information about the potential damage distribution within a city is invaluable. In this work, we present a first attempt to design an early warning and rapid response procedure for real-time risk assessment. In particular, the procedure uses typical real-time information (i.e., P-wave arrival times and early waveforms) derived from a regional seismic network for locating and evaluating the size of an earthquake, information which in turn is exploited for extracting a risk map representing the potential distribution of damage from a dataset of predicted scenarios compiled for the target city. A feasibility study of the procedure is presented for the city of Bishkek, the capital of Kyrgyzstan, which is surrounded by the Kyrgyz seismic network by mimicking the ground motion associated with two historical events that occurred close to Bishkek, namely the 1911 Kemin (M?=?8.2; ±0.2) and the 1885 Belovodsk (M?=?6.9; ±0.5) earthquakes. Various methodologies from previous studies were considered when planning the implementation of the early warning and rapid response procedure for real-time risk assessment: the Satriano et al. (Bull Seismol Soc Am 98(3):1482–1494, 2008) approach to real-time earthquake location; the Caprio et al. (Geophys Res Lett 38:L02301, 2011) approach for estimating moment magnitude in real time; the EXSIM method for ground motion simulation (Motazedian and Atkinson, Bull Seismol Soc Am 95:995–1010, 2005); the Sokolov (Earthquake Spectra 161: 679–694, 2002) approach for estimating intensity from Fourier amplitude spectra; and the Tyagunov et al. (Nat Hazard Earth Syst Sci 6:573–586, 2006) approach for risk computation. Innovatively, all these methods are jointly applied to assess in real time the seismic risk of a particular target site, namely the city of Bishkek. Finally, the site amplification and vulnerability datasets considered in the proposed methodology are taken from previous studies, i.e., Parolai et al. (Bull Seismol Soc Am, 2010) and Bindi et al. (Soil Dyn Earthq Eng, 2011), respectively.     

History of Modern Earthquake Hazard Mapping and Assessment in California Using a Deterministic or Scenario Approach

Modern earthquake ground motion hazard mapping in California began following the 1971 San Fernando earthquake in the Los Angeles metropolitan area of southern California. Earthquake hazard assessment followed a traditional approach, later called Deterministic Seismic Hazard Analysis (DSHA) in order to distinguish it from the newer Probabilistic Seismic Hazard Analysis (PSHA). In DSHA, seismic hazard in the event of the Maximum Credible Earthquake (MCE) magnitude from each of the known seismogenic faults within and near the state are assessed. The likely occurrence of the MCE has been assumed qualitatively by using late Quaternary and younger faults that are presumed to be seismogenic, but not when or within what time intervals MCE may occur. MCE is the largest or upper-bound potential earthquake in moment magnitude, and it supersedes and automatically considers all other possible earthquakes on that fault. That moment magnitude is used for estimating ground motions by applying it to empirical attenuation relationships, and for calculating ground motions as in neo-DSHA (Zuccolo et al., 2008). The first deterministic California earthquake hazard map was published in 1974 by the California Division of Mines and Geology (CDMG) which has been called the California Geological Survey (CGS) since 2002, using the best available fault information and ground motion attenuation relationships at that time. The California Department of Transportation (Caltrans) later assumed responsibility for printing the refined and updated peak acceleration contour maps which were heavily utilized by geologists, seismologists, and engineers for many years. Some engineers involved in the siting process of large important projects, for example, dams and nuclear power plants, continued to challenge the map(s). The second edition map was completed in 1985 incorporating more faults, improving MCE??s estimation method, and using new ground motion attenuation relationships from the latest published results at that time. CDMG eventually published the second edition map in 1992 following the Governor??s Board of Inquiry on the 1989 Loma Prieta earthquake and at the demand of Caltrans. The third edition map was published by Caltrans in 1996 utilizing GIS technology to manage data that includes a simplified three-dimension geometry of faults and to facilitate efficient corrections and revisions of data and the map. The spatial relationship of fault hazards with highways, bridges or any other attribute can be efficiently managed and analyzed now in GIS at Caltrans. There has been great confidence in using DSHA in bridge engineering and other applications in California, and it can be confidently applied in any other earthquake-prone region. Earthquake hazards defined by DSHA are: (1) transparent and stable with robust MCE moment magnitudes; (2) flexible in their application to design considerations; (3) can easily incorporate advances in ground motion simulations; and (4) economical. DSHA and neo-DSHA have the same approach and applicability. The accuracy of DSHA has proven to be quite reasonable for practical applications within engineering design and always done with professional judgment. In the final analysis, DSHA is a reality-check for public safety and PSHA results. Although PSHA has been acclaimed as a better approach for seismic hazard assessment, it is DSHA, not PSHA, that has actually been used in seismic hazard assessment for building and bridge engineering, particularly in California.     

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.     

A Crustal Model for the Eastern Alps Region and a New Moho Map in Southeastern Europe

In the last two decades, south-central Europe and the Eastern Alps have been widely explored by many seismic refraction experiments (e.g., CELEBRATION 2000, ALP 2002, SUDETES 2003). Although quite detailed images are available along linear profiles, a comprehensive, three-dimensional crustal model of the region is still missing. This limitation makes this region a weak spot in continental-wide comprehensive representations of crustal structure. To improve on this situation, we select and collect 37 published active-source seismic lines in this region. After geo-referencing each line, we sample them along vertical profiles—every 50?km or less along the line—and derive P-wave velocities in a stack of homogeneous layers (separated by discontinuities: depth of crystalline basement, top of lower crust, and Moho). We finally merge the information using geostatistical methods, and infer S-wave velocity and density using empirical scaling relations. We present here the resulting crustal model for a region encompassing the Eastern Alps, Dinarides, Pannonian basin, Western Carpathians and Bohemian Massif, covering the region within $45^{\circ}\text{--}51^{\circ}\hbox{N}$ and $11^{\circ} \text{--} 22^{\circ}\hbox{E}$ with a resolution of $0.2^{\circ} \times 0.2^{\circ}.$ We are also able to extend and update the map of Moho depth in a wider region within $35^{\circ}\text{--}51^{\circ}\hbox{N}$ and $12^{\circ}\text{--}45^{\circ}\hbox{E},$ gathering Moho values from the collected seismic lines, other published dataset and using the European plate reference EPcrust as a background. All the digitized profiles and the resulting model are available online.     

Attenuation Tomography Based on Strong Motion Data: Case Study of Central Honshu Region, Japan

Three-dimensional frequency dependent S-wave quality factor (Q_β(f)) value for the central Honshu region of Japan has been determined in this paper using an algorithm based on inversion of strong motion data. The method of inversion for determination of three-dimensional attenuation coefficients is proposed by Hashida and Shimazaki (J Phys Earth. 32, 299–316, 1984) and has been used and modified by Joshi (Curr Sci. 90, 581–585, 2006; Nat Hazards. 43, 129–146, 2007) and Joshi et al. (J. Seismol. 14, 247–272, 2010). Twenty-one earthquakes digitally recorded on strong motion stations of Kik-net network have been used in this work. The magnitude of these earthquake ranges from 3.1 to 4.2 and depth ranging from 5 to 20 km, respectively. The borehole data having high signal to noise ratio and minimum site effect is used in the present work. The attenuation structure is determined by dividing the entire area into twenty-five three-dimensional blocks of uniform thickness having different frequency-dependent shear wave quality factor. Shear wave quality factor values have been determined at frequencies of 2.5, 7.0 and 10 Hz from record in a rectangular grid defined by 35.4°N to 36.4°N and 137.2°E to 138.2°E. The obtained attenuation structure is compared with the available geological features in the region and comparison shows that the obtained structure is capable of resolving important tectonic features present in the area. The proposed attenuation structure is compared with the probabilistic seismic hazard map of the region and shows that it bears some remarkable similarity in the patterns seen in seismic hazard map.     

Declustering of Iran earthquake catalog (1983–2017) using the epidemic-type aftershock sequence (ETAS) model

The main goal of this article is to decluster Iranian plateau seismic catalog by the epidemic-type aftershock sequence (ETAS) model and compare the results with some older methods. For this purpose, Iranian plateau bounded in 24°–42°N and 43°–66°E is subdivided into three major tectonic zones: (1) North of Iran (2) Zagros (3) East of Iran. The extracted earthquake catalog had a total of 6034 earthquakes (Mw?>?4) in the time span 1983–2017. The ETAS model is an accepted stochastic approach for seismic evaluation and declustering earthquake catalogs. However, this model has not yet been used to decluster the seismic catalog of Iran. Until now, traditional methods like the Gardner and Knopoff space–time window method and the Reasenberg link-based method have been used in most studies for declustering Iran earthquake catalog. Finally, the results of declustering by the ETAS model are compared with result of Gardner and Knopoff (Bull Seismol Soc Am 64(5):1363–1367, 1974), Uhrhammer (Earthq Notes 57(1):21, 1986), Gruenthal (pers. comm.) and Reasenberg (Geophys Res 90:5479–5495, 1985) declustering methods. The overall conclusion is difficult, but the results confirm the high ability of the ETAS model for declustering Iranian earthquake catalog. Use of the ETAS model is still in its early steps in Iranian seismological researches, and more parametric studies are needed.     

Seismic risk and damage prediction: case of the buildings in Constantine city (Algeria)

Located at the North-Eastern part of Algeria (Tellian Atlas), Constantine has crucial administrative, economic, scientific and cultural importance. It has continuously experienced significant urban evolutions during the different periods of its history. The city is located in an active seismic region within Algeria and has been struck in the past by several moderate and strong earthquakes. The strongest earthquake recorded since the beginning of instrumental seismology took place on October 27, 1985 with a magnitude M \(_\mathrm{S}=\) 5.9. Constantine presents a high seismic risk, because of its dense housing and high population density (2,374 inhabitants/km \(^{2})\) . This requires a risk assessment in order to take preventive measures and reduce the losses in case of potential major earthquake. For this purpose, a scenario based approach is considered. The building damage assessment methodology adopted for the Algerian context is adapted from HAZUS approach. In the present case, the effective Algerian seismic code response spectrum (RPA 99/2003) is considered as a seismic hazard model. The prediction of the expected damages is performed for a set of almost 29,000 buildings.     

Slope instabilities triggered by the 11th May 2011 Lorca earthquake (Murcia,Spain): comparison to previous hazard assessments and proposition of a new hazard map and probability of failure equation

The Lorca Basin has been the object of recent research aimed at studying the phenomena of earthquake-induced landslides and its assessment in the frame of different seismic scenarios. However, it has not been until the 11th May 2011 Lorca earthquakes when it has been possible to conduct a systematic approach to the problem. In this paper we present an inventory of slope instabilities triggered by the Lorca earthquakes which comprises more than 100 cases, mainly rock and soil falls of small size (1–100  \(\hbox {m}^{3}\) ). The distribution of these instabilities is here compared to two different earthquake-triggered landslide hazard maps: one considering the occurrence of the most probable earthquake for a 475-years return period in the Lorca Basin \((\hbox {M}_{\mathrm{w}}=5.0)\) based on both low- and high-resolution digital elevation model (DEM); and a second one matching the occurrence of the \(\hbox {M}_{\mathrm{w}}=5.2\) 2011 Lorca earthquake, which was performed using the higher resolution DEM. The most frequent Newmark displacements related to the slope failures triggered by the 2011 Lorca earthquakes are lower than 2 cm in both the hazard scenarios considered. Additionally, the predicted Newmark displacements were correlated to the inventory of slope instabilities to develop a probability of failure equation. The fit seems to be very good since most of the mapped slope failures are located on the higher probability areas. The probability of slope failure in the Lorca Basin for a seismic event similar to the \(\hbox {M}_{\mathrm{w}}\) 5.2 2011 Lorca earthquake can be considered as very low (0–4 %).