Showing papers in "Journal of Seismology in 2018"

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

Abstract: 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 Mw 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.

Earthquake chemical precursors in groundwater: a review

Abstract: We review changes in groundwater chemistry as precursory signs for earthquakes. In particular, we discuss pH, total dissolved solids (TDS), electrical conductivity, and dissolved gases in relation to their significance for earthquake prediction or forecasting. These parameters are widely believed to vary in response to seismic and pre-seismic activity. However, the same parameters also vary in response to non-seismic processes. The inability to reliably distinguish between changes caused by seismic or pre-seismic activities from changes caused by non-seismic activities has impeded progress in earthquake science. Short-term earthquake prediction is unlikely to be achieved, however, by pH, TDS, electrical conductivity, and dissolved gas measurements alone. On the other hand, the production of free hydroxyl radicals (•OH), subsequent reactions such as formation of H2O2 and oxidation of As(III) to As(V) in groundwater, have distinctive precursory characteristics. This study deviates from the prevailing mechanical mantra. It addresses earthquake-related non-seismic mechanisms, but focused on the stress-induced electrification of rocks, the generation of positive hole charge carriers and their long-distance propagation through the rock column, plus on electrochemical processes at the rock-water interface.

Seismicity in the block mountains between Halle and Leipzig, Central Germany: centroid moment tensors, ground motion simulation, and felt intensities of two M ≈ 3 earthquakes in 2015 and 2017

Abstract: On April 29, 2017 at 0:56 UTC (2:56 local time), an MW = 2.8 earthquake struck the metropolitan area between Leipzig and Halle, Germany, near the small town of Markranstadt. The earthquake was felt within 50 km from the epicenter and reached a local intensity of I0 = IV. Already in 2015 and only 15 km northwest of the epicenter, a MW = 3.2 earthquake struck the area with a similar large felt radius and I0 = IV. More than 1.1 million people live in the region, and the unusual occurrence of the two earthquakes led to public attention, because the tectonic activity is unclear and induced earthquakes have occurred in neighboring regions. Historical earthquakes south of Leipzig had estimated magnitudes up to MW ≈ 5 and coincide with NW-SE striking crustal basement faults. We use different seismological methods to analyze the two recent earthquakes and discuss them in the context of the known tectonic structures and historical seismicity. Novel stochastic full waveform simulation and inversion approaches are adapted for the application to weak, local earthquakes, to analyze mechanisms and ground motions and their relation to observed intensities. We find NW-SE striking normal faulting mechanisms for both earthquakes and centroid depths of 26 and 29 km. The earthquakes are located where faults with large vertical offsets of several hundred meters and Hercynian strike have developed since the Mesozoic. We use a stochastic full waveform simulation to explain the local peak ground velocities and calibrate the method to simulate intensities. Since the area is densely populated and has sensitive infrastructure, we simulate scenarios assuming that a 12-km long fault segment between the two recent earthquakes is ruptured and study the impact of rupture parameters on ground motions and expected damage.

Characterisation of ground motion recording stations in the Groningen gas field

Abstract: The seismic hazard and risk analysis for the onshore Groningen gas field requires information about local soil properties, in particular shear-wave velocity (VS). A fieldwork campaign was conducted at 18 surface accelerograph stations of the monitoring network. The subsurface in the region consists of unconsolidated sediments and is heterogeneous in composition and properties. A range of different methods was applied to acquire in situ VS values to a target depth of at least 30 m. The techniques include seismic cone penetration tests (SCPT) with varying source offsets, multichannel analysis of surface waves (MASW) on Rayleigh waves with different processing approaches, microtremor array, cross-hole tomography and suspension P-S logging. The offset SCPT, cross-hole tomography and common midpoint cross-correlation (CMPcc) processing of MASW data all revealed lateral variations on length scales of several to tens of metres in this geological setting. SCPTs resulted in very detailed VS profiles with depth, but represent point measurements in a heterogeneous environment. The MASW results represent VS information on a larger spatial scale and smooth some of the heterogeneity encountered at the sites. The combination of MASW and SCPT proved to be a powerful and cost-effective approach in determining representative VS profiles at the accelerograph station sites. The measured VS profiles correspond well with the modelled profiles and they significantly enhance the ground motion model derivation. The similarity between the theoretical transfer function from the VS profile and the observed amplification from vertical array stations is also excellent.

GeoInt: the first macroseismic intensity database for the Republic of Georgia

Abstract: Our work is intended to present the new macroseismic intensity database for the Republic of Georgia—hereby named GeoInt—which includes earthquakes from the historical (from 1250 B.C. onwards) to the instrumental era. Such database is composed of 111 selected earthquakes and related 3944 intensity data points (IDPs) for 1509 different localities, reported in the Medvedev-Sponheuer-Karnik scale (MSK). Regarding the earthquakes, the MS is in the 3.3–7 range and the depth is in the 2–36 km range. The entire set of IDPs is characterized by intensities ranging from 2–3 to 9–10 and covers an area spanning from 39.508° N to 45.043° N in a N-S direction and from 37.324° E to 48.500° E in an E-W direction, with some of the IDPs located outside the Georgian border, in the (i) Republic of Armenia, (ii) Russian Federation, (iii) Republic of Turkey, and (iv) Republic of Azerbaijan. We have revised each single IDP and have reevaluated and homogenized intensity values to the MSK scale. In particular, regarding the whole set of 3944 IDPs, 348 belong to the Historical era (pre-1900) and 3596 belong to the instrumental era (post-1900). With particular regard to the 3596 IDPs, 105 are brand new (3%), whereas the intensity values for 804 IDPs have been reevaluated (22%); for 2687 IDPs (75%), intensities have been confirmed from previous interpretations. We introduce this database as a key input for further improvements in seismic hazard modeling and seismic risk calculation for this region, based on macroseismic intensity; we report all the 111 earthquakes with available macroseismic information. The GeoInt database is also accessible online at and will be kept updated in the future.

Influence of wind turbines on seismic stations in the upper rhine graben, SW Germany

Abstract: By analysing long- and short-term seismological measurements at wind farms close to the town of Landau, SW Germany, we present new insights into ground motion signals from wind turbines (WTs) at local seismic stations. Because of their need to be located in similar regions with sparsely anthropogenic activities, wind turbines impact seismic stations and their recordings in a way that is not yet fully understood by researchers. To ensure the undisturbed recording tasks of a regional seismic array or a single station by a protected area around those endangered stations, it is very important to investigate the behavior of WTs as a seismic source. For that reason, we calculate averaged one-hour long spectra of the power spectral density (PSD) before and after the installation of a new wind farm within the investigated area. These PSD are ordered according to the rotation speed. We observe a clear increase of the PSD level after the WT installation in a frequency range of 0.5 to 10 Hz up to a distance of 5.5 km away from the WT. By analysing seismic borehole data, we also observe a decrease of the PSD of wind dependent signals with depth. The impact of wind-dependent signals is found to be much more pronounced for the shallower station (150 m depth) than for the deeper one (305 m depth). Using short-term profile measurements, we fit a power-law decay proportional to 1/r b to the main WT-induced PSD peaks and differentiate between near-field and far-field effects of ground motions. For low frequencies in the range from 1 to 4 Hz, we determine a b value of 0.78 to 0.85 for the far field, which is consistent with surface waves. The b value increases (up to 1.59) with increasing frequencies (up to 5.5 Hz), which is obviously due to attenuating effects like scattering or anelasticity. These results give a better understanding of the seismic wavefield interactions between wind turbines (or wind farms) with nearby seismic stations, including borehole installations, in a sedimentary setting.

Earthquake induced liquefaction hazard, probability and risk assessment in the city of Kolkata, India: its historical perspective and deterministic scenario

Abstract: Liquefaction-induced ground failure is one amongst the leading causes of infrastructure damage due to the impact of large earthquakes in unconsolidated, non-cohesive, water saturated alluvial terrains. The city of Kolkata is located on the potentially liquefiable alluvial fan deposits of Ganga-Bramhaputra-Meghna Delta system with subsurface litho-stratigraphic sequence comprising of varying percentages of clay, cohesionless silt, sand, and gravel interbedded with decomposed wood and peat. Additionally, the region has moderately shallow groundwater condition especially in the post-monsoon seasons. In view of burgeoning population, there had been unplanned expansion of settlements in the hazardous geological, geomorphological, and hydrological conditions exposing the city to severe liquefaction hazard. The 1897 Shillong and 1934 Bihar-Nepal earthquakes both of M w 8.1 reportedly induced Modified Mercalli Intensity of IV–V and VI–VII respectively in the city reportedly triggering widespread to sporadic liquefaction condition with surface manifestation of sand boils, lateral spreading, ground subsidence, etc., thus posing a strong case for liquefaction potential analysis in the terrain. With the motivation of assessing seismic hazard, vulnerability, and risk of the city of Kolkata through a consorted federal funding stipulated for all the metros and upstart urban centers in India located in BIS seismic zones III, IV, and V with population more than one million, an attempt has been made here to understand the liquefaction susceptibility condition of Kolkata under the impact of earthquake loading employing modern multivariate techniques and also to predict deterministic liquefaction scenario of the city in the event of a probabilistic seismic hazard condition with 10% probability of exceedance in 50 years and a return period of 475 years. We conducted in-depth geophysical and geotechnical investigations in the city encompassing 435 km2 area. The stochastically synthesized bedrock ground motion for both the 1897 and 1934 earthquakes on non-linear analysis of local site conditions through DEEPSOIL Geotechnical analysis package present surface level peak ground acceleration of the order of 0.05–0.14 g for the 1934 Bihar-Nepal earthquake while for the 1897 Shillong earthquake it is found to be in the range of 0.03–0.11 g. The factor of safety (FOS) against liquefaction, the probability of liquefaction (P L), the liquefaction potential index (LPI), and the liquefaction risk index are estimated under the influence of these two earthquakes wherein the city is classified into severe (LPI > 15), high (5 < LPI ≤ 15), moderate (0 < LPI ≤ 5), and non-liquefiable (LPI = 0) susceptibility zones. While the 1934 Bihar-Nepal earthquake induced moderate to severe liquefaction hazard condition in the city in mostly the deltaic plain and interdistributary marsh geomorphologic units with 13.5% sites exhibiting moderate hazard with a median LPI of 1.8, 8.5% sites depicting high with a median LPI of 9.1 and 4% sites with a median LPI of 18.9 exhibiting severe hazard condition, 1897 Shillong earthquake induced mostly non-liquefaction condition with very few sites depicting moderate and high liquefaction hazard. A conservative liquefaction hazard scenario of the city on the other hand estimated through deterministic approach for 10% probability of exceedance in 50 years predicts a high hazard zone in the 3.5–19 m depth region with FOS < 1 and P L > 65% comprising of coarse-grained sediments of sand, silty sand, and clayey silty sand in mostly the deltaic plain geomorphologic unit with 39.1% sites depicting severe liquefaction hazard with a median LPI of 28.3. A non-linear regression analysis on both the historical and deterministic liquefaction scenarios in P L versus LPI domain with ± 1 standard deviation confidence bound generated a cubic polynomial relationship between the two liquefaction hazard proxies. This study considered a bench mark for other cities in the country and elsewhere forms an integral part of the mega-seismic microzonation endeavors undertaken in all the earthquake-prone counties in the world.

Strong ground motion from the November 12, 2017, M 7.3 Kermanshah earthquake in western Iran

Abstract: In this paper, we analyzed the strong ground motion from the November 12, 2017, Kermanshah earthquake in western Iran with moment magnitude (M) of 7.3. Nonlinear and linear amplification of ground motion amplitudes were observed at stations with soft soil condition at hypocentral distances below and above 100 km, respectively. Observation of large ground motion amplitudes dominated with long-period pulses on the strike-normal component of the velocity time series suggests a right-lateral component of movement and propagation of rupture towards southeast. Comparison of the horizontal peak ground acceleration (PGA) from the M 7.3 earthquake with global PGA values showed a similar decay in ground motion amplitudes, although it seems that PGA from the M 7.3 Kermanshah earthquake is higher than global values for NEHRP site class B. We also found that the bracketed duration (Db) was higher in the velocity domain than in the acceleration domain for the same modified Mercalli intensity (MMI) threshold. For example, Db reached ~ 30 s at the maximum PGA while it was ~ 50 s at the maximum peak ground velocity above the threshold of MMI = 5. Although the standard design spectrum from Iranian Code of Practice for Seismic Resistant Design of Buildings (standard No. 2800) seems to include appropriate values for the design of structures with fundamental period of 1 s and higher, it is underestimated for near-field ground motions at lower periods.

Updated earthquake catalogue for seismic hazard analysis in Pakistan

Abstract: A reliable and homogenized earthquake catalogue is essential for seismic hazard assessment in any area. This article describes the compilation and processing of an updated earthquake catalogue for Pakistan. The earthquake catalogue compiled in this study for the region (quadrangle bounded by the geographical limits 40–83° N and 20–40° E) includes 36,563 earthquake events, which are reported as 4.0–8.3 moment magnitude (MW) and span from 25 AD to 2016. Relationships are developed between the moment magnitude and body, and surface wave magnitude scales to unify the catalogue in terms of magnitude MW. The catalogue includes earthquakes from Pakistan and neighbouring countries to minimize the effects of geopolitical boundaries in seismic hazard assessment studies. Earthquakes reported by local and international agencies as well as individual catalogues are included. The proposed catalogue is further used to obtain magnitude of completeness after removal of dependent events by using four different algorithms. Finally, seismicity parameters of the seismic sources are reported, and recommendations are made for seismic hazard assessment studies in Pakistan.

The relation between M E , M L and M w in theory and numerical simulations for small to moderate earthquakes

Abstract: This article addresses the question of how moment magnitude, Mw, local magnitude, ML, and energy magnitude, ME, of small to moderate earthquakes would scale with each other in an ideal elastic medium and why they scale differently in reality. It turns out that ML, in the way it is commonly determined, is a poor and inconsistent measure of earthquake size. For moderate to large events, the situation could be improved with a new magnitude, MA, that is not biased by the Wood-Anderson response. Mw is robust and can in principle be determined in a consistent way over the entire imaginable magnitude range, but it accounts only for the static component of the earthquake source. Could ME be an alternative? A simple conceptual analysis reveals the relation between radiated seismic energy, seismic moment, static stress drop, apparent stress, radiation efficiency and rupture velocity. Numerical simulations based on an extended-source model show that, to estimate the radiated seismic energy (and thus also ME), one needs to take into account the effects of rupture directivity and radiation pattern. Anelastic attenuation and the bandwidth of the recording instrument determine the magnitude limit below which pulse widths and corner frequencies remain nearly constant. Below this limit, all information about the dynamics of the source is lost and attempts to correct for this in the presence of noise and aleatory signal variability, even with well-calibrated attenuation models, are probably futile. We thus have to accept the fact that ME (just as ML) is a different measure of earthquake size for small earthquakes than for moderate to large earthquakes.

Time-dependent seismic hazard analysis for the Greater Tehran and surrounding areas

Abstract: This study presents a time-dependent approach for seismic hazard in Tehran and surrounding areas. Hazard is evaluated by combining background seismic activity, and larger earthquakes may emanate from fault segments. Using available historical and paleoseismological data or empirical relation, the recurrence time and maximum magnitude of characteristic earthquakes for the major faults have been explored. The Brownian passage time (BPT) distribution has been used to calculate equivalent fictitious seismicity rate for major faults in the region. To include ground motion uncertainty, a logic tree and five ground motion prediction equations have been selected based on their applicability in the region. Finally, hazard maps have been presented.

Site classification of Indian strong motion network using response spectra ratios

Abstract: In the present study, we tried to classify the Indian strong motion sites spread all over Himalaya and adjoining region, located on varied geological formations, based on response spectral ratio. A total of 90 sites were classified based on 395 strong motion records from 94 earthquakes recorded at these sites. The magnitude of these earthquakes are between 2.3 and 7.7 and the hypocentral distance for most of the cases is less than 50 km. The predominant period obtained from response spectral ratios is used to classify these sites. It was found that the shape and predominant peaks of the spectra at these sites match with those in Japan, Italy, Iran, and at some of the sites in Europe and the same classification scheme can be applied to Indian strong motion network. We found that the earlier schemes based on description of near-surface geology, geomorphology, and topography were not able to capture the effect of sediment thickness. The sites are classified into seven classes (CL-I to CL-VII) with varying predominant periods and ranges as proposed by Alessandro et al. (Bull Seismol Soc Am 102:680-695 2012). The effect of magnitudes and hypocentral distances on the shape and predominant peaks were also studied and found to be very small. The classification scheme is robust and cost-effective and can be used in region-specific attenuation relationships for accounting local site effect.

Probabilistic seismic hazard assessment of southern part of Ghana

Abstract: This paper presents a seismic hazard map for the southern part of Ghana prepared using the probabilistic approach, and seismic hazard assessment results for six cities. The seismic hazard map was prepared for 10% probability of exceedance for peak ground acceleration in 50 years. The input parameters used for the computations of hazard were obtained using data from a catalogue that was compiled and homogenised to moment magnitude (Mw). The catalogue covered a period of over a century (1615–2009). The hazard assessment is based on the Poisson model for earthquake occurrence, and hence, dependent events were identified and removed from the catalogue. The following attenuation relations were adopted and used in this study—Allen (for south and eastern Australia), Silva et al. (for Central and eastern North America), Campbell and Bozorgnia (for worldwide active-shallow-crust regions) and Chiou and Youngs (for worldwide active-shallow-crust regions). Logic-tree formalism was used to account for possible uncertainties associated with the attenuation relationships. OpenQuake software package was used for the hazard calculation. The highest level of seismic hazard is found in the Accra and Tema seismic zones, with estimated peak ground acceleration close to 0.2 g. The level of the seismic hazard in the southern part of Ghana diminishes with distance away from the Accra/Tema region to a value of 0.05 g at a distance of about 140 km.

On the possible correlation between the Gutenberg-Richter parameters of the frequency-magnitude relationship

Abstract: The study of the Gutenberg-Richter (GR) parameters a and b has been very important to describe and characterize the seismicity over the different seismic provinces around the world. As far as we know, the possible correlation between the GR parameters a and b has not received enough attention. Bayrak et al. reported the a and b values for 27 active seismic regions around the boundaries of the main tectonic plates of the world. From these data, we found that there exists a positive correlation between the a and b parameters (R = 0.85, R2 = 0.72). On the other hand, we made around 150 computer runs of a spring-block model proposed by Olami et al. (Phys Rev Lett 68(8):1244–1247, 1992). This model roughly emulates the interaction between two fault planes and it reaches a self-organized critical state. With these simulations, we also found that the a and b parameters are positively correlated. Motivated by these results, we propose an analytical demonstration that indeed a and b are positively correlated. In addition, we discuss on other possible applications of the spring-block model to actual seismicity and to frictional experiments made with sandpapers.

Enhancing seismic P phase arrival picking based on wavelet denoising and kurtosis picker

Abstract: P phase arrival picking of weak signals is still challenging in seismology. A wavelet denoising is proposed to enhance seismic P phase arrival picking, and the kurtosis picker is applied on the wavelet-denoised signal to identify P phase arrival. It has been called the WD-K picker. The WD-K picker, which is different from those traditional wavelet-based pickers on the basis of a single wavelet component or certain main wavelet components, takes full advantage of the reconstruction of main detail wavelet components and the approximate wavelet component. The proposed WD-K picker considers more wavelet components and presents a better P phase arrival feature. The WD-K picker has been evaluated on 500 micro-seismic signals recorded in the Chinese Yongshaba mine. The comparison between the WD-K pickings and manual pickings shows the good picking accuracy of the WD-K picker. Furthermore, the WD-K picking performance has been compared with the main detail wavelet component combining-based kurtosis (WDC-K) picker, the single wavelet component-based kurtosis (SW-K) picker, and certain main wavelet component-based maximum kurtosis (MMW-K) picker. The comparison has demonstrated that the WD-K picker has better picking accuracy than the other three-wavelet and kurtosis-based pickers, thus showing the enhanced ability of wavelet denoising.

The instrumental seismicity of the Barents and Kara sea region: relocated event catalog from early twentieth century to 1989

Abstract: We have relocated seismic events registered within the Barents and Kara sea region from early twentieth century to 1989 with a view to creating a relocated catalog. For the relocation, we collected all available seismic bulletins from the global network using data from the ISC Bulletin (International Seismological Centre), ISC-GEM project (International Seismological Centre–Global Earthquake Model), EuroSeismos project, and by Soviet seismic stations from Geophysical Survey of the Russian Academy of Sciences. The location was performed by applying a modified method of generalized beamforming. We have considered several travel time models and selected one with the best location accuracy for ground truth events. Verification of the modified method and selection of the travel time model were performed using data on four nuclear explosions that occurred in the area of the Novaya Zemlya Archipelago and in the north of the European part of Russia. The modified method and the Barents travel time model provide sufficient accuracy for event location in the region. The relocation procedure was applied to 31 of 36 seismic events registered within the Barents and Kara sea region.

Risk-targeted maps for Romania

Abstract: Romania has one of the highest seismic hazard levels in Europe. The seismic hazard is due to a combination of local crustal seismic sources, situated mainly in the western part of the country and the Vrancea intermediate-depth seismic source, which can be found at the bend of the Carpathian Mountains. Recent seismic hazard studies have shown that there are consistent differences between the slopes of the seismic hazard curves for sites situated in the fore-arc and back-arc of the Carpathian Mountains. Consequently, in this study we extend this finding to the evaluation of the probability of collapse of buildings and finally to the development of uniform risk-targeted maps. The main advantage of uniform risk approach is that the target probability of collapse will be uniform throughout the country. Finally, the results obtained are discussed in the light of a recent study with the same focus performed at European level using the hazard data from SHARE project. The analyses performed in this study have pointed out to a dominant influence of the quantile of peak ground acceleration used for anchoring the fragility function. This parameter basically alters the shape of the risk-targeted maps shifting the areas which have higher collapse probabilities from eastern Romania to western Romania, as its exceedance probability increases. Consequently, a uniform procedure for deriving risk-targeted maps appears as more than necessary.

Tutorial review of seismic surface waves’ phenomenology

Abstract: In recent years, surface wave seismology has become one of the leading directions in seismological investigations of the Earth’s structure and seismic sources. Various applications cover a wide spectrum of goals, dealing with differences in sources of seismic excitation, penetration depths, frequency ranges, and interpretation techniques. Observed seismic data demonstrates the great variability of phenomenology which can produce difficulties in interpretation for beginners. This tutorial review is based on the many years’ experience of authors in processing and interpretation of seismic surface wave observations and the lectures of one of the authors (ALL) at Workshops on Seismic Wave Excitation, Propagation and Interpretation held at the Abdus Salam International Center for Theoretical Physics (Trieste, Italy) in 1990–2012. We present some typical examples of wave patterns which could be encountered in different applications and which can serve as a guide to analysis of observed seismograms.

Spatial variations of b -value and crustal stress in the Pyrenees

Abstract: The seismicity in the Pyrenees is continuous and well surveyed since more than 20 years. We use the catalogue of seismicity between 1997 and 2013 to explore the spatial variations of the b-value, which corresponds to the slope of the frequency-magnitude distribution of the earthquakes. Especially, variations of the b-value characterise the state of stress of the crust, possibly highlighting a deficit of large earthquake occurrence. We estimate the differential crustal stress from the b-value using a relationship published by Scholz (Geophys Res Lett 42:1399–1402, 2015). We also estimate the stress drop variations by determining a power law which links the magnitude to the seismic source radius in the Pyrenees. We focus on the depth variations and we analyse vertical profiles of b-value, differential stress and stress drop, first in the whole Pyrenean belt, then in 10 subregions. The b-values are generally smaller than 1, except in the uppermost 3–5 km where the obtained high values could be linked to the presence of fluids. Downward, the b-values decrease slowly or remain constant until a depth of increase, which could correspond to the brittle-ductile limit of the crust. We propose that this depth and the regional and vertical variations of the b-values are related to the regional tectonic context and possibly to the density heterogeneities. We also suggest that stress drop and differential stress are linearly correlated and that the stress drop is at least 1.8‰ of the differential stress.

Joint inversion of GNSS and teleseismic data for the rupture process of the 2017 M w 6.5 Jiuzhaigou, China, earthquake

Abstract: The spatio-temporal slip distribution of the earthquake that occurred on 8 August 2017 in Jiuzhaigou, China, was estimated from the teleseismic body wave and near-field Global Navigation Satellite System (GNSS) data (coseismic displacements and high-rate GPS data) based on a finite fault model. Compared with the inversion results from the teleseismic body waves, the near-field GNSS data can better restrain the rupture area, the maximum slip, the source time function, and the surface rupture. The results show that the maximum slip of the earthquake approaches 1.4 m, the scalar seismic moment is ~ 8.0 × 1018 N·m (Mw ≈ 6.5), and the centroid depth is ~ 15 km. The slip is mainly driven by the left-lateral strike-slip and it is initially inferred that the seismogenic fault occurs in the south branch of the Tazang fault or an undetectable fault, a NW-trending left-lateral strike-slip fault, and belongs to one of the tail structures at the easternmost end of the eastern Kunlun fault zone. The earthquake rupture is mainly concentrated at depths of 5–15 km, which results in the complete rupture of the seismic gap left by the previous four earthquakes with magnitudes > 6.0 in 1973 and 1976. Therefore, the possibility of a strong aftershock on the Huya fault is low. The source duration is ~ 30 s and there are two major ruptures. The main rupture occurs in the first 10 s, 4 s after the earthquake; the second rupture peak arrives in ~ 17 s. In addition, the Coulomb stress study shows that the epicenter of the earthquake is located in the area where the static Coulomb stress change increased because of the 12 May 2017 Mw7.9 Wenchuan, China, earthquake. Therefore, the Wenchuan earthquake promoted the occurrence of the 8 August 2017 Jiuzhaigou earthquake.

Ground motions around a deep semielliptic canyon with a horizontal edge subjected to incident plane SH waves

Abstract: We study scattering of antiplane shear waves induced by a deep semielliptic canyon with a horizontal edge We employ the region-point-matching technique to cope with the problem considered Through an auxiliary boundary, a part of the circumference of a semiellipse, the whole analyzed region is divided into two subregions We express the displacement fields in terms of Mathieu functions We unify two distinct elliptic coordinates via a simple coordinate transformation relation Integration of the coordinate transformation relation into the region-point-matching technique simplifies the procedure for constructing simultaneous equations Imposing the continuity conditions and traction-free ones, we obtain the expansion coefficients Frequency-domain results demonstrate ground motion variability based on several key factors Ground surface responses under seismic shaking are also simulated in the time domain

Revised seismotectonic model of NE Italy and W Slovenia based on focal mechanism inversion

Abstract: The study is focusing on the stress and strain inversions from focal mechanisms in a revised seismotectonic zonation of northeastern Italy and western Slovenia. The recent increase of monitoring capability of the local seismic network, the updated geological-structural model of the area, and the novelties emerged from studies on the spatial organization of the seismicity allowed a redefinition of the seismotectonic zones. The stress and strain tensors inversion is inferred from 203 focal mechanisms, corresponding to earthquakes occurred between 1984 and 2016 with coda-duration magnitude range from 2.0 to 5.6. The inverted stress domains reveal an articulated picture of the interaction of the Adria microplate with the Eurasian plate. A dominant strike-slip stress field characterizes the eastern part of the area, while the seismotectonic zones of the central part are undergoing to thrusting regime. The stress pattern inferred in the western part of the study area outlines a complex picture with prevailing strike-slip regime and dominant compression only in a seismotectonic zone. The comparison of stress and strain tensor orientations evidences a relative uniformity of the crustal strength in the eastern and northwestern zones of the study area. The central and western zones appear to be characterized by planes of mechanical weakness not favorably oriented for failure with respect to the stress tensor.

Seismotectonics of south-west Dominican Republic using recent data

Abstract: The Dominican Republic has a high level of seismic activity, and a new seismic network has been installed to improve the detection of earthquakes. The network has been operated by Instituto Politecnico Loyola since 2012. It uses six new stations of its own, as well as 17 other stations publically available. In this study, we investigate in particular the south-west of the country where no seismic stations have been operating before. We find an area, SE of the Enriquillo Plantain Garden fault (responsible for the 2010 Haiti M7.0 earthquake), with a high seismic activity (M = 2–4) compared to the surrounding area. This shallow seismicity (except two events around 40 km depth) is not seen in any global catalogs, and it does not seem to be associated with any known faults. However, the region has been under rapid recent uplift since margins between hills and valleys are filled with massive alluvial fans. We made 24 new fault plane solutions using P-polarities and S/P amplitude ratios. The solutions show mainly reverse mechanisms and the P-axis directions are mainly NNE oriented, which is in agreement with a published strain direction from GPS measurements. We conclude that the main cause of the seismicity in our study area is the push of the Beata Ridge against Enriquillo basin and the rest of the Gonave microplate, reflected in the mainly reverse focal mechanisms.

New data on earthquake focal mechanisms in the Laptev Sea region of the Arctic-Asian seismic belt

Abstract: We consider 16 earthquakes with Mw = 4.2–5.2 that occurred in the south-eastern part of the Laptev Sea shelf, Lena River Delta, and North Verkhoyanye (Russia) in 1990–2014. Focal mechanisms, scalar seismic moments, moment magnitudes, and hypocentral depths of the seismic events have been calculated from the data on amplitude spectra of surface waves and P wave first-motion polarities. The obtained results sufficiently implement the existing dataset on reliable earthquake source parameters for the study region and prove the change of the stress-strain state of the crust from extension on the Laptev Sea shelf to compression on the continent providing finer spatial details of the deformation field in the transition zones such as Buor-Khaya Bay and the Lena River Delta.

Evaluation of local site effect from microtremor measurements in Babol City, Iran

Abstract: Every year, numerous casualties and a large deal of financial losses are incurred due to earthquake events. The losses incurred by an earthquake vary depending on local site effect. Therefore, in order to conquer drastic effects of an earthquake, one should evaluate urban districts in terms of the local site effect. One of the methods for evaluating the local site effect is microtremor measurement and analysis. Aiming at evaluation of local site effect across the city of Babol, the study area was gridded and microtremor measurements were performed with an appropriate distribution. The acquired data was analyzed through the horizontal-to-vertical noise ratio (HVNR) method, and fundamental frequency and associated amplitude of the H/V peak were obtained. The results indicate that fundamental frequency of the study area is generally lower than 1.25 Hz, which is acceptably in agreement with the findings of previous studies. Also, in order to constrain and validate the seismostratigraphic model obtained with this method, the results were compared with geotechnical, geological, and seismic data. Comparing the results of different methods, it was observed that the presented geophysical method can successfully determine the values of fundamental frequency across the study area as well as local site effect. Using the data obtained from the analysis of microtremor, a microzonation map of fundamental frequency across the city of Babol was prepared. This map has numerous applications in designing high-rise building and urban development plans.

Do French macroseismic intensity observations agree with expectations from the European Seismic Hazard Model 2013

Abstract: Probabilistic seismic hazard assessments are the basis of modern seismic design codes. To test fully a seismic hazard curve at the return periods of interest for engineering would require many thousands of years’ worth of ground-motion recordings. Because strong-motion networks are often only a few decades old (e.g. in mainland France the first accelerometric network dates from the mid-1990s), data from such sensors can be used to test hazard estimates only at very short return periods. In this article, several hundreds of years of macroseismic intensity observations for mainland France are interpolated using a robust kriging-with-a-trend technique to establish the earthquake history of every French mainland municipality. At 24 selected cities representative of the French seismic context, the number of exceedances of intensities IV, V and VI is determined over time windows considered complete. After converting these intensities to peak ground accelerations using the global conversion equation of Caprio et al. (Ground motion to intensity conversion equations (GMICEs): a global relationship and evaluation of regional dependency, Bulletin of the Seismological Society of America 105:1476–1490, 2015), these exceedances are compared with those predicted by the European Seismic Hazard Model 2013 (ESHM13). In half of the cities, the number of observed exceedances for low intensities (IV and V) is within the range of predictions of ESHM13. In the other half of the cities, the number of observed exceedances is higher than the predictions of ESHM13. For intensity VI, the match is closer, but the comparison is less meaningful due to a scarcity of data. According to this study, the ESHM13 underestimates hazard in roughly half of France, even when taking into account the uncertainty in the conversion from intensity to acceleration. However, these results are valid only for the acceleration range tested in this study (0.01 to 0.09 g).

Regionalization and dependence of coda Q on frequency and lapse time in the seismically active Peloritani region (northeastern Sicily, Italy)

Abstract: The Peloritani region is one of the most seismically active regions in Italy and, consequently, the quantification of attenuation of the medium plays an important role for seismic risk evaluation. Moreover, it is necessary for the prediction of earth ground motion and future seismic source studies. An in depth analysis has been made here to understand the frequency and lapse time dependence of attenuation characteristics of the region by using the coda of local earthquakes. A regionalization is likewise performed in order to investigate the spatial variation of coda Q across the whole region. Finally, our results are jointly interpreted with those obtained from recently published 3D velocity tomographies for further insights.

Introduction of conditional mean spectrum and conditional spectrum in the practice of seismic safety evaluation in China

Abstract: This paper aims at implementing and introducing the use of conditional mean spectrum (CMS) and conditional spectrum (CS) as the main input parameters in the practice of seismic safety evaluation (SSE) in China, instead of the currently used uniform hazard spectrum (UHS). For this purpose, a procedure for M-R-epsilon seismic hazard deaggregation in China was first developed. For illustration purposes, two different typical sites in China, with one to two dominant seismic zones, were considered as examples to carry out seismic hazard deaggregation and illustrate the construction of CMS/CS. Two types of correlation coefficients were used to generate CMS and the results were compared over a vibration period range of interest. Ground motion records were selected from the NSMONS (2007–2015) and PEER NGA-West2 databases to correspond to the target CMS and CS. Hazard consistency of the spectral accelerations of the selected ground motion records was evaluated and validated by computing the annual exceedance probability rate of the response spectra and comparing the results to the hazard curve corresponding to each site of concern at different periods. The tools developed in this work and their illustrative application to specific case studies in China are a first step towards the adoption of CMS and CS into the practice of seismic safety evaluation in this country.

Determination of GMPE functional form for an active region with limited strong motion data: application to the Himalayan region

Abstract: Advancement in the seismic networks results in formulation of different functional forms for developing any new ground motion prediction equation (GMPE) for a region. Till date, various guidelines and tools are available for selecting a suitable GMPE for any seismic study area. However, these methods are efficient in quantifying the GMPE but not for determining a proper functional form and capturing the epistemic uncertainty associated with selection of GMPE. In this study, the compatibility of the recent available functional forms for the active region is tested for distance and magnitude scaling. Analysis is carried out by determining the residuals using the recorded and the predicted spectral acceleration values at different periods. Mixed effect regressions are performed on the calculated residuals for determining the intra- and interevent residuals. Additionally, spatial correlation is used in mixed effect regression by changing its likelihood function. Distance scaling and magnitude scaling are respectively examined by studying the trends of intraevent residuals with distance and the trend of the event term with magnitude. Further, these trends are statistically studied for a respective functional form of a ground motion. Additionally, genetic algorithm and Monte Carlo method are used respectively for calculating the hinge point and standard error for magnitude and distance scaling for a newly determined functional form. The whole procedure is applied and tested for the available strong motion data for the Himalayan region. The functional form used for testing are five Himalayan GMPEs, five GMPEs developed under NGA-West 2 project, two from Pan-European, and one from Japan region. It is observed that bilinear functional form with magnitude and distance hinged at 6.5 M w and 300 km respectively is suitable for the Himalayan region. Finally, a new regression coefficient for peak ground acceleration for a suitable functional form that governs the attenuation characteristic of the Himalayan region is derived.

Characterizing the active tectonics in the Oran region (Algeria) and recasting the 1790 earthquake

Abstract: In this work, we reappraise the seismogenic potential of the geologic structures in the western Tell Atlas of Algeria, considered active host to moderate to low magnitude earthquakes. The direct identification of active faults is generally a difficult task in northern Algeria. The active tectonics in the Oran Plio-Quaternary age basin (Northwestern Algeria) is analyzed and characterized through a morpho-structural study combining topographic, geomorphologic, geological, and neotectonic data. Folds and fault scarps affecting Quaternary deposits show that the region is affected by compressional deformation still active nowadays, as shown by the recorded seismic activity. Our new observations enable a better understanding of the present seismotectonic context of the Oran region, particularly with regard to the magnitude and source of the 1790 Oran damaging event. The obtained result helps to shed some light on the elusive active tectonics characterizing this coastal area, and to assess regional seismic hazard, particularly in coastal zones where large seismogenic areas straddle the onshore–offshore zones.