A new <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> estimation parameter for use in earthquake early warning systems

We propose a method that employs the squared displacement integral (ID2) to estimate earthquake magnitudes in real time for use in earthquake early warning (EEW) systems. Moreover, using τ _c and P _d for comparison, we establish formulas for estimating the moment magnitudes of these three parameters based on the selected aftershocks (4.0 ≤ M _s  ≤ 6.5) of the 2008 Wenchuan earthquake. In this comparison, the proposed ID2 method displays the highest accuracy. Furthermore, we investigate the applicability of the initial parameters to large earthquakes by estimating the magnitude of the Wenchuan M _s 8.0 mainshock using a 3-s time window. Although these three parameters all display problems with saturation, the proposed ID2 parameter is relatively accurate. The evolutionary estimation of ID2 as a function of the time window shows that the estimation equation established with ID2 ^Ref determined from the first 8-s of P wave data can be directly applicable to predicate the magnitudes of 8.0. Therefore, the proposed ID2 parameter provides a robust estimator of earthquake moment magnitudes and can be used for EEW purposes.     

Body wave attenuation characteristics in the crust of Alborz region and North Central Iran

Attenuation of P and S waves has been investigated in Alborz and north central part of Iran using the data recorded by two permanent and one temporary networks during October 20, 2009, to December 22, 2010. The dataset consists of 14,000 waveforms from 380 local earthquakes (2 < M _L < 5.6). The extended coda normalization method (CNM) was used to estimate quality factor of P (Q _P) and S waves (Q _S) at seven frequency bands (0.375, 0.75, 1.5, 3, 6, 12, 24 Hz). The Q _P and Q _S values have been estimated at lapse times from 40 to 100 s. It has been observed that the estimated values of Q _P and Q _S are time independent; therefore, the mean values of Q _P and Q _S at different lapse times have been considered. The frequency dependence of quality factor was determined by using a power-law relationship. The frequency-dependent relationship for Q _P was estimated in the form of (62 ± 7)f ^{(1.03 ± 0.07)} and (48 ± 5)f ^{(0.95 ± 0.07)} in Alborz region and North Central Iran, respectively. These relations for Q _S for Alborz region and North Central Iran have estimated as (83 ± 8)f ^{(0.99 ± 0.07)} and (68 ± 5)f ^{(0.96 ± 0.05)}, respectively. The observed low Q values could be the results of thermoelastic effects and/or existing fracture. The estimated frequency-dependent relationships are comparable with tectonically active regions.     

Frictional and viscous effects on the nucleation phase of an earthquake

We study the frictional and viscous effects on earthquake nucleation, especially for the nucleation phase, based on a one-degree-of-freedom spring-slider model with friction and viscosity. The frictional and viscous effects are specified by the characteristic displacement, U _c, and viscosity coefficient, η, respectively. Simulation results show that friction and viscosity can both lengthen the natural period of the system and viscosity increases the duration time of motion of the slider. Higher viscosity causes a smaller amplitude of lower velocity motion than lower viscosity. A change of either U _c (under large η) or η (under large U _c) from a large value (U _ch for U _c and η _h for η) to a small one (U _cl for U _c and η _l for η) in two stages during sliding can result in a clear nucleation phase prior to the P-wave. The differences δU _c = U _ch ? U _cl and δη = η _h ? η _l are two important factors in producing a nucleation phase. The difference between the nucleation phase and the P-wave increases with either δU _c or δη. Like seismic observations, the peak amplitude of P-wave, which is associated with the earthquake magnitude, is independent upon the duration time of nucleation phase. A mechanism specified with a change of either η or U _c from a larger value to a smaller one due to temporal variations in pore fluid pressure and temperature in the fault zone based on radiation efficiency is proposed to explain the simulation results and observations.     

A study of Guptkashi,Uttarakhand earthquake of 6 February 2017 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>5.3) in the Himalayan arc and implications for ground motion estimation

The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M₀?=?1.12×10¹⁷ Nm, M _w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ _s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω^?2 source model) by attenuation parameters Q(f)?=?500f^0.9, κ?=?0.04 s, and f_max?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M _w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.     

Source parameters of the 27th of June 2015 Gulf of Aqaba earthquake

On the 27 June 2015, at 15:34:03 UTC, a moderate-sized earthquake of M _w 5.0 occurred in the Gulf of Aqaba. Using teleseismic P waves, the focal mechanism of the mainshock was investigated by two techniques. The first technique used the polarities of the first P wave onsets, and the second technique was based on the normalized waveform modeling technique. The results showed that the extension stress has a NE orientation with a shallow southward plunge while the compression stress has a NW trend with a nearly shallow westward plunge, obtaining a strike-slip mechanism. This result agrees well with the typical consequence of crustal deformation resulting from the ongoing extensional to shear stress regime in the Gulf of Aqaba (NE-SW extension and NW-SE compression). The grid search method over a range of focal depths indicates an optimum solution at 15 ± 1 km. To identify the causative fault plane, the aftershock hypocenters were relocated using the local waveform data and the double-difference technique. Considering the fault trends, the spatial distribution of relocated aftershocks demarcated a NS-oriented causative fault, in consistence with one of the nodal planes of the focal mechanism solution, emphasizing the dominant stress regime in the region. Following the Brune model, the estimates of source parameters exhibited fault lengths of 0.29 ≤ L ≤ 2.48 km, moment magnitudes of 3.0 ≤ M _w ≤ 5.0, and stress drops of 0.14 ≤ Δσ < 1.14 MPa, indicating a source scaling similar to the tectonic earthquakes related to plate boundaries.     

Attenuation coefficients of Rayleigh and <Emphasis Type="Italic">Lg</Emphasis> waves

Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q _e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s^???1 (for Rayleigh waves) and 8 ×10^???3 s^???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ _i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.     

Detection of Antarctic oscillation signals in Earth’s oblateness variations

Variations of Earth’s oblateness (J ₂) reflect a large scale mass redistribution within the Earth system. The climate effect causing J ₂ interannual variations is still not clear, though previous studies indicated it may be related to EI Niño–Southern Oscillation (ENSO) and Pacific Decadal Oscillation (PDO). However, we have a new discovery of the significant Antarctic oscillation (AAO) signals in J ₂ interannual variations, especially on 4–6 year scales based on cross wavelet and wavelet coherence analysis with 95% confidence test during 1979–2012. The results additionally indicate that the close phase relationship between J ₂ and AAO (AAO leading J ₂ variations by 3 ± 2 months in phase) is far superior to that between J ₂ and ENSO/PDO on 4–6 year scales. In this work, we discuss, for the first time, a possible geophysical mechanism of AAO effecting J ₂ variations. The investigations are based on the definition of AAO and its spatial–temporal behavior influencing the large-scale mass movement. Finally, an approximate quantitative estimate of the AAO imprint on J ₂ with an emphasis on the atmospheric contribution is made.     

Near-source attenuation of high-frequency body waves beneath the New Madrid Seismic Zone

Attenuation characteristics in the New Madrid Seismic Zone (NMSZ) are estimated from 157 local seismograph recordings out of 46 earthquakes of 2.6?≤?M?≤?4.1 with hypocentral distances up to 60 km and focal depths down to 25 km. Digital waveform seismograms were obtained from local earthquakes in the NMSZ recorded by the Center for Earthquake Research and Information (CERI) at the University of Memphis. Using the coda normalization method, we tried to determine Q values and geometrical spreading exponents at 13 center frequencies. The scatter of the data and trade-off between the geometrical spreading and the quality factor did not allow us to simultaneously derive both these parameters from inversion. Assuming 1/R ^1.0 as the geometrical spreading function in the NMSZ, the Q _P and Q _S estimates increase with increasing frequency from 354 and 426 at 4 Hz to 729 and 1091 at 24 Hz, respectively. Fitting a power law equation to the Q estimates, we found the attenuation models for the P waves and S waves in the frequency range of 4 to 24 Hz as Q _P?=?(115.80?±?1.36) f ^{(0.495?±?0.129)} and Q _S?=?(161.34?±?1.73) f ^{(0.613?±?0.067)}, respectively. We did not consider Q estimates from the coda normalization method for frequencies less than 4 Hz in the regression analysis since the decay of coda amplitude was not observed at most bandpass filtered seismograms for these frequencies. Q _S/Q _P?>?1, for 4?≤?f?≤?24 Hz as well as strong intrinsic attenuation, suggest that the crust beneath the NMSZ is partially fluid-saturated. Further, high scattering attenuation indicates the presence of a high level of small-scale heterogeneities inside the crust in this region.     

Temporal and spatial distributions of precursory seismicity rate changes in the Thailand-Laos-Myanmar border region: implication for upcoming hazardous earthquakes

To study the prospective areas of upcoming strong-to-major earthquakes, i.e., M _w  ≥ 6.0, a catalog of seismicity in the vicinity of the Thailand-Laos-Myanmar border region was generated and then investigated statistically. Based on the successful investigations of previous works, the seismicity rate change (Z value) technique was applied in this study. According to the completeness earthquake dataset, eight available case studies of strong-to-major earthquakes were investigated retrospectively. After iterative tests of the characteristic parameters concerning the number of earthquakes (N) and time window (T _w ), the values of 50 and 1.2 years, respectively, were found to reveal an anomalous high Z-value peak (seismic quiescence) prior to the occurrence of six out of the eight major earthquake events studied. In addition, the location of the Z-value anomalies conformed fairly well to the epicenters of those earthquakes. Based on the investigation of correlation coefficient and the stochastic test of the Z values, the parameters used here (N = 50 events and T _w  = 1.2 years) were suitable to determine the precursory Z value and not random phenomena. The Z values of this study and the frequency-magnitude distribution b values of a previous work both highlighted the same prospective areas that might generate an upcoming major earthquake: (i) some areas in the northern part of Laos and (ii) the eastern part of Myanmar.     

Interaction of geometry and mechanical property of trapezoidal sedimentary basins with incident <Emphasis Type="Italic">SH</Emphasis> waves

This paper investigated the effects of basin geometry and material property on the response of 2D trapezoidal sediment-filled basin to incident plane SH waves. Ten basin configurations with different geometries were developed, and then their seismic responses to both Ricker wavelets and seismic records were simulated by using an explicit finite difference scheme. The definition of deep/shallow basin, the precondition for the observation of prominent surface waves and the influential area of edge effects of the shallow basin were discussed quantitatively in this study. The followings were concluded: in the common velocity contrast range (v _s1/v _s2 < 10), the fundamental frequency a basin with W/H > 3.0 can be estimated approximately by 1D theory. The complexity of peak ground acceleration distribution pattern, the width of the most affected section as well as the amplitude of ground motion in the Edge Region increase with incident frequency. Prominent surface waves can only be observed when the incident wavelength is shorter than the critical wavelength λ _c . The interaction between incident wave and basin dynamic property plays a dominant role on the peak ground acceleration amplitude while the interaction between incident wave and geometry plays a more significant role on the peak ground acceleration distribution. For very shallow basin, different areas along the basin width are affected to different extents. Only a limited area close to the basin edge is influenced significantly. It is more feasible to propose spectral aggravation factor for different surface zones respectively than a uniform constant as a tool to calibrate the 1D-based design spectrum so as to take the basin effects into account.     

Three-dimensional admittance analysis of lithospheric elastic thickness over the Louisville Ridge

Using bathymetry and altimetric gravity anomalies, a 1° 9 1° lithospheric effective elastic thickness(Te) model over the Louisville Ridge and its adjacent regions is calculated using the moving window admittance technique. For comparison, three bathymetry models are used: general bathymetric charts of the oceans, SIO V15.1,and BAT_VGG. The results show that BAT_VGG is more suitable for calculating T e than the other two models. T e along the Louisville Ridge was re-evaluated. The southeast of the ridge has a medium Te of 10–20 km, while Te increases dramatically seaward of the Tonga-Kermadec trench as a result of the collision of the Pacific and IndoAustralian plates.     

Thermoremanent and chemical magnetization of exsolution products of nanosized titanomagnetites

Based on the theory of two-phase interacting nanoparticles, the formation of thermoremanent and chemical remanent magnetization in nanosized titanomagnetites is modeled. It is shown that the value of thermoremanent magnetization barely depends on the degree of titanomagnetite exsolution whereas, chemical remanent magnetization which emerges during the exsolution increases up to at most the value of thermoremanent magnetization. The values of the ratio of thermoremanent to ideal magnetization, R _t , fall within the limits 0.8 ≤ R _t ≤ 1. The analogous ratio of chemical remanent magnetization to the ideal R _c are below R _t at all stages of the exsolution. Besides, the magnetic interaction between the nanoparticles reduces the values of thermoremanent and chemical magnetization but barely affects the ratio.     

<Emphasis Type="Italic">Q</Emphasis><Subscript>c</Subscript> and <Emphasis Type="Italic">Q</Emphasis><Subscript>S</Subscript> wave attenuation of South African earthquakes

Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M _L ?≤?4.4. Up to five seismograph stations were utilised to determine Q _c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q ₀ and α were determined for the attenuation relation Q _c(f)?=?Q ₀ f ^α . The result was Q ₀?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t _s???t ₀) (time from origin time t ₀ to the start of coda analysis window is 1.9 times the S-travel time, t _s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q _c was not calculated. The derived Q _c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q _S was Q ₀?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q _c result.     

Gradual late stage deepening of Gega ice-dammed lake,Tsangpo gorge,southeastern Tibet,indicated by preliminary sedimentary rock magnetic properties

Gega lake, in southeastern Tibet, was formed by the blocking of Yarlung Tsangpo gorge by a glacier and is a well-known example of geomorphological damming. However, the evolution of the damming process at the site is still not understood in detail. Here, we use measurements of multiple magnetic parameters of the sediments from the Yusong (YS) 3 section, which is well-dated by optically stimulated luminescence, to provide a detailed history of the late stage of Gega dammed lake since 17.0 ka. Low-frequency field magnetic susceptibility (χ_lf) increases upwards gradually from 25 to 79?×?10^?8 m³/kg above 5.5 m, but other magnetic properties, such as frequency-dependent magnetic susceptibility (χ_fd % and χ_fd), susceptibility of anhysteretic remanent magnetization (χ_ARM), and saturation isothermal remanent magnetization (SIRM) did not show a similar degree of enhancement. The magnetic grain size indicators of χ_ARM/χ_lf, χ_ARM/SIRM, and χ_lf/SIRM all indicate a trend of increasing magnetic grain size from the 5.5 m to the top of section; however, the bulk sediment grain size decreases gradually within the same interval. The total organic carbon is very low (0.2–0.7%), and thus, it is unlikely that the sedimentary environment is sufficiently strongly reducing to lead to the dissolution of magnetic minerals. Therefore, we infer that the coarsening-upwards of the magnetic grain size, and the increasing magnetic susceptibility from 5.5 m, reflect the gradual preferential preservation of magnetic minerals, caused by the deepening of the lake. Thus, we interpret the magnetic record of the section as reflecting the gradual deepening of the late stage of Gega ice-dammed lake. Last, the abrupt disappearance of the lake may have resulted in an outburst flood.     

Seasonal variation of carbon exchange of typical forest ecosystems along the eastern forest transect in China

The long-term and continuous carbon fluxes of Changbaishan temperate mixed forest (CBS), Qianyanzhou subtropical evergreen coniferous forest (QYZ), Dinghushan subtropical evergreen mixed forest (DHS) and Xishuangbana tropical rainforest (XSBN) have been measured with eddy covariance techniques. In 2003, different responses of carbon exchange to the environment appeared across the four ecosystems. At CBS, the carbon exchange was mainly determined by radiation and temperature. 0°C and 10°C were two important temperature thresholds; the former determined the length of the growing season and the latter affected the magnitude of carbon exchange. The maximum net ecosystem exchange (N _EE) of CBS occurred in early summer because maximum ecosystem photosynthesis (G _PP) occurred earlier than maximum ecosystem respiration (R _e). During summer, QYZ experienced severe drought and N _EE decreased significantly mainly as a result of the depression of G _PP. At DHS and XSBN, N _EE was higher in the drought season than the wet season, especially the conversion between carbon sink and source occurring during the transition season at XSBN. During the wet season, increased fog and humid weather resulted from the plentiful rainfall, the ecosystem G _PP was dispressed. The Q ₁₀ and annual respiration of XSBN were the highest among the four ecosystems, while the average daily respiration of CBS during the growing season was the highest. Annual N _EE of CBS, QYZ, DHS and XSBN were 181.5, 360.9, 536.2 and ?320.0 g·C·m^?2·a^?1, respectively. From CBS to DHS, the temperature and precipitation increased with the decrease in latitude. The ratio of N _EE/R _e increased with latitude, while R _e/G _PP, ecosystem light use efficiency (L _UE), precipitation use efficiency and average daily G _PP decreased gradually. However, XSBN usually escaped such latitude trend probably because of the influence of the south-west monsoon climate which does not affect the other ecosystems. Long-term measurement and more research were necessary to understand the adaptation of forest ecosystems to climate change and to evaluate the ecosystem carbon balance due to the complexity of structure and function of forest ecosystems.     

An evaluation of earthquake hazard parameters in the Iranian Plateau based on the Gumbel III distribution

The Gumbel’s third asymptotic distribution (GIII) of the extreme value method is employed to evaluate the earthquake hazard parameters in the Iranian Plateau. This research quantifies spatial mapping of earthquake hazard parameters like annual and 100-year mode beside their 90 % probability of not being exceeded (NBE) in the Iranian Plateau. Therefore, we used a homogeneous and complete earthquake catalogue during the period 1900–2013 with magnitude M _w ? ≥?4.0, and the Iranian Plateau is separated into equal area mesh of 1° late?×?1° long. The estimated result of annual mode with 90 % probability of NBE is expected to exceed the values of M _w 6.0 in the Eastern part of Makran, most parts of Central and East Iran, Kopeh Dagh, Alborz, Azerbaijan, and SE Zagros. The 100-year mode with 90 % probability of NBE is expected to overpass the value of M _w 7.0 in the Eastern part of Makran, Central and East Iran, Alborz, Kopeh Dagh, and Azerbaijan. The spatial distribution of 100-year mode with 90 % probability of NBE uncovers the high values of earthquake hazard parameters which are frequently connected with the main tectonic regimes of the studied area. It appears that there is a close communication among the seismicity and the tectonics of the region.     

Evaluation of the damping modification factor for structures subjected to near-fault ground motions

The damping modification factor (DMF) has been extensively used in earthquake engineering to describe the variation of structural responses due to varied damping ratios. It is known that DMFs are dependent not only on structural dynamic properties but also on characteristics of ground motions. DMFs regulated in current seismic codes are generally developed based on far-fault ground motions and are inappropriately used in structural design where pulse-like near-fault ground motions are involved. In this paper, statistical investigation of the DMF is performed based on 50 carefully selected pulse-like near-fault ground motions. It is observed that DMFs for pulse-like ground motions exhibit significant dependence on the pulse period T _p in a specific period range. If the period of the structure in response is close to the pulse period, the DMF attains the same level as that derived from far-fault ground motions; as the period of the structure is considerably larger or smaller than the pulse period T _p , the response reduction effect by the increased damping ratio is generally small, except for large earthquakes with long pulse periods, which exhibit significant reduction of response for structures with periods smaller than T _p . Based on the statistical results of DMFs, the empirical formulas for estimating DMFs for displacement, velocity and acceleration spectra are proposed, the effect of structural period, pulse period and damping ratio are considered in the formulas, and the formulas are designed to satisfy the specific reliability requirement in the period range of 0.1 < T/T _p  < 1, which is of engineering interest.     

Probabilistic aftershock hazard analysis,two case studies in West and Northwest Iran

Aftershock hazard maps contain the essential information for search and rescue process, and re-occupation after a main-shock. Accordingly, the main purposes of this article are to study the aftershock decay parameters and to estimate the expected high-frequency ground motions (i.e., Peak Ground Acceleration (PGA)) for recent large earthquakes in the Iranian plateau. For this aim, the Ahar-Varzaghan doublet earthquake (August 11, 2012; M _N =6.5, M _N =6.3), and the Ilam (Murmuri) earthquake (August 18, 2014 ; M _N =6.2) have been selected. The earthquake catalogue has been collected based on the Gardner and Knopoff (Bull Seismol Soc Am 64(5), 1363-1367, 1974) temporal and spatial windowing technique. The magnitude of completeness and the seismicity parameters (a,??b) and the modified Omori law parameters (P,??K,??C) have been determined for these two earthquakes in the 14, 30, and 60 days after the mainshocks. Also, the temporal changes of parameters (a,??b,??P,??K,??C) have been studied. The aftershock hazard maps for the probability of exceedance (33%) have been computed in the time periods of 14, 30, and 60 days after the Ahar-Varzaghan and Ilam (Murmuri) earthquakes. For calculating the expected PGA of aftershocks, the regional and global ground motion prediction equations have been utilized. Amplification factor based on the site classes has also been implied in the calculation of PGA. These aftershock hazard maps show an agreement between the PGAs of large aftershocks and the forecasted PGAs. Also, the significant role of b parameter in the Ilam (Murmuri) probabilistic aftershock hazard maps has been investigated.     

Region-specific <Emphasis Type="Italic">S</Emphasis>-wave attenuation for earthquakes in northwestern Iran

In this study, continuous wavelet transform is applied to estimate the frequency-dependent quality factor of shear waves, Q _S , in northwestern Iran. The dataset used in this study includes velocigrams of more than 50 events with magnitudes between 4.0 and 6.5, which have occurred in the study area. The CWT-based method shows a high-resolution technique for the estimation of S-wave frequency-dependent attenuation. The quality factor values are determined in the form of a power law as Q _S (f)?=?(147?±?16)f ^0.71?±?0.02 and (126?±?12)f ^0.73?±?0.02 for vertical and horizontal components, respectively, where f is between 0.9 and 12 Hz. Furthermore, in order to verify the reliability of the suggested Q _S estimator method, an additional test is performed by using accelerograms of Ahar-Varzaghan dual earthquakes on August 11, 2012, of moment magnitudes 6.4 and 6.3 and their aftershocks. Results indicate that the estimated Q _S values from CWT-based method are not very sensitive to the numbers and types of waveforms used (velocity or acceleration).     

The ratio between corner frequencies of source spectra of <Emphasis Type="Italic">P</Emphasis>- and <Emphasis Type="Italic">S</Emphasis>-waves—a new discriminant between earthquakes and quarry blasts

This study analyzes and compares the P- and S-wave displacement spectra from local earthquakes and explosions of similar magnitudes. We propose a new approach to discrimination between low-magnitude shallow earthquakes and explosions by using ratios of P- to S-wave corner frequencies as a criterion. We have explored 2430 digital records of the Israeli Seismic Network (ISN) from 456 local events (226 earthquakes, 230 quarry blasts, and a few underwater explosions) of magnitudes Md?=?1.4–3.4, which occurred at distances up to 250 km during 2001–2013 years. P-wave and S-wave displacement spectra were computed for all events following Brune’s source model of earthquakes (1970, 1971) and applying the distance correction coefficients (Shapira and Hofstetter, Teconophysics 217:217–226, 1993; Ataeva G, Shapira A, Hofstetter A, J Seismol 19:389-401, 2015), The corner frequencies and moment magnitudes were determined using multiple stations for each event, and then the comparative analysis was performed.The analysis showed that both P-wave and especially S-wave displacement spectra of quarry blasts demonstrate the corner frequencies lower than those obtained from earthquakes of similar magnitudes. A clear separation between earthquake and explosion populations was obtained for ratios of P- to S-wave corner frequency f ₀(P)/f ₀(S). The ratios were computed for each event with corner frequencies f ₀ of P- and S-wave, which were obtained from the measured f ₀ ^I at individual stations, then corrected for distance and finally averaged. We obtained empirically the average estimation of f ₀(P)/f ₀(S)?=?1.23 for all used earthquakes, and 1.86 for all explosions. We found that the difference in the ratios can be an effective discrimination parameter which does not depend on estimated moment magnitude M _w .The new multi-station Corner Frequency Discriminant (CFD) for earthquakes and explosions in Israel was developed based on ratios P- to S-wave corner frequencies f ₀(P)/f ₀(S), with the empirical threshold value of the ratio for Israel as 1.48.