Strong ground motion from the November 12, 2017, <Emphasis Type="Bold">M</Emphasis> 7.3 Kermanshah earthquake in western Iran

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 (D_b) was higher in the velocity domain than in the acceleration domain for the same modified Mercalli intensity (MMI) threshold. For example, D_b 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.     

Seismicity anomalies before the great earthquake of <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 in the Kunlun Pass and its significance to earthquake prediction

A great earthquake of M _S=8.1 took place in the west of Kunlun Pass on November 14, 2001. The epicenter is located at 36.2°N and 90.9°E. The analysis shows that some main precursory seismic patterns appear before the great earthquake, e.g., seismic gap, seismic band, increased activity, seismicity quiet and swarm activity. The evolution of the seismic patterns before the earthquake of M _S=8.1 exhibits a course very similar to that found for earthquake cases with M _S≥7. The difference is that anomalous seismicity before the earthquake of M _S=8.1 involves in the larger area coverage and higher seismic magnitude. This provides an evidence for recognizing precursor and forecasting of very large earthquake. Finally, we review the rough prediction of the great earthquake and discuss some problems related to the prediction of great earthquakes.     

Ground deformation effects from the ~M6 earthquakes (2014–2015) on Cephalonia–Ithaca Islands (Western Greece) deduced by GPS observations

The implications of the earthquakes that took place in the central Ionian Islands in 2014 (Cephalonia, M _w6.1, M _w5.9) and 2015 (Lefkas, M _w6.4) are described based on repeat measurements of the local GPS networks in Cephalonia and Ithaca, and the available continuous GPS stations in the broader area. The Lefkas earthquake occurred on a branch of the Cephalonia Transform Fault, affecting Cephalonia with SE displacements gradually decreasing from north (~100 mm) to south (~10 mm). This earthquake revealed a near N–S dislocation boundary separating Paliki Peninsula in western Cephalonia from the rest of the island, as well as another NW–SE trending fault that separates kinematically the northern and southern parts of Paliki. Strain field calculations during the interseismic period (2014–2015) indicate compression between Ithaca and Cephalonia, while extension appears during the following co-seismic period (2015–2016) including the 2015 Lefkas earthquake. Additional tectonically active zones with differential kinematic characteristics were also identified locally.     

Geodynamic and Seismotectonic Activity in Eastern Tibet in the 21st Century

At the beginning of the 21st century, a series of great earthquakes were recorded in northeastern Tibet, along the periphery of the Bayan Hara lithospheric block. An earthquake with M_S = 8.1 occurred within the East Kunlun fault zone in the Kunlun Mountains, which caused an extended surface rupture with left-lateral strike slip. An earthquake with M_S = 8 occurred in Wenchuan (China) on May 12, 2008, giving rise to an extended overthrust along the Lunmanshan fault zone. An earthquake with M_S = 7.1 occurred in Yushu (China) on April 14, 2010; its epicenter was on the Grazze–Yushu–Funchuoshan fault; a left-lateral strikeslip offset was observed on the surface. An earthquake with M_S = 7 occurred in the vicinity of Lushan on April 20, 2013; its epicenter was within the Lunmanshan fault zone, 103 km southwest of the zone of the catastrophic Wenchuan earthquake. An earthquake with M_S = 8.2 occurred in Nepal on April 25, 2015. Based on the CSN seismic catalog, the energy of all earthquakes in eastern Tibet at the end of the 20th and beginning of the 21st centuries was estimated. It was found that Tibet was seismically quiet from 1980 to 2000. The beginning of the 21st century has been marked by seismic activation with earthquake sources migrating southward to surround the Bayan Hara lithospheric block from every quarter. Therefore, this block can be regarded as one of the most seismically active regions of China.     

3D <Emphasis Type="Italic">v</Emphasis><Subscript>P</Subscript> and <Emphasis Type="Italic">v</Emphasis><Subscript>S</Subscript> models of southeastern margin of the Tibetan plateau from joint inversion of body-wave arrival times and surface-wave dispersion data

A new 3D velocity model of the crust and upper mantle in the southeastern (SE) margin of the Tibetan plateau was obtained by joint inversion of body- and surface-wave data. For the body-wave data, we used 7190 events recorded by 102 stations in the SE margin of the Tibetan plateau. The surface-wave data consist of Rayleigh wave phase velocity dispersion curves obtained from ambient noise cross-correlation analysis recorded by a dense array in the SE margin of the Tibetan plateau. The joint inversion clearly improves the v _S model because it is constrained by both data types. The results show that at around 10 km depth there are two low-velocity anomalies embedded within three high-velocity bodies along the Longmenshan fault system. These high-velocity bodies correspond well with the Precambrian massifs, and the two located to the northeast of 2013 M _S 7.0 Lushan earthquake are associated with high fault slip areas during the 2008 Wenchuan earthquake. The aftershock gap between 2013 Lushan earthquake and 2008 Wenchuan earthquake is associated with low-velocity anomalies, which also acts as a barrier zone for ruptures of two earthquakes. Generally large earthquakes (M ≥ 5) in the region occurring from 2008 to 2015 are located around the high-velocity zones, indicating that they may act as asperities for these large earthquakes. Joint inversion results also clearly show that there exist low-velocity or weak zones in the mid-lower crust, which are not evenly distributed beneath the SE margin of Tibetan plateau.     

GPS time-series and its response to <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 Kunlunshan earthquake

In this paper, observation data in 25 GPS reference stations of China have been analyzed by calculating GPS position coordinate time-series with GIPSY. Result shows there is an obvious trend variation in such time-series. The trend variations of time series along the longitude and latitude coordinate reflect the motion of each position in the global-plate, in which the trend variation in the vertical direction reveals some large-scale construction information or reflects the local movement around the positions. The analysis also shows that such time-series have a variation cycle of nearly 1.02 a, but the reason still remains to be further studied. At the end of this paper, response of the time-series of M _S=8.1 Kunlunshan earthquake was analyzed, and the seismogenic process of M _S=8.1 Kunlunshan earthquake, according to the time proceeding and the feature of anomaly, was divided into 3 phases—changes in blocks with forces, strain accumulation, quick accumulation and slow release of energy. At the initial stage of seismogenic process of M _S=8.1 earthquake and at the imminent earthquake, coseismic process as well as during the post earthquake recovery, anomaly in vertical direction is always in a majority. The anomalous movement in vertical direction at the initial stage resulted in a blocking between faults, while at the middle stage of seismogenic process, the differential movement between blocks are in a majority, which is the major reason causing energy accumulating at the blocking stage of faults.     

Historical intermediate-depth earthquakes in the southern Aegean Sea Benioff zone: modeling their anomalous macroseismic patterns with stochastic ground-motion simulations

We model the macroseismic damage distribution of four important intermediate-depth earthquakes of the southern Aegean Sea subduction zone, namely the destructive 1926 M?=?7.7 Rhodes and 1935 M?=?6.9 Crete earthquakes, the unique 1956 M?=?6.9 Amorgos aftershock (recently proposed to be triggered by a shallow event), and the more recent 2002 M?=?5.9 Milos earthquake, which all exhibit spatially anomalous macroseismic patterns. Macroseismic data for these events are collected from published macroseismic databases and compared with the spatial distribution of seismic motions obtained from stochastic simulation, converted to macroseismic intensity (Modified Mercalli scale, I_MM). For this conversion, we present an updated correlation between macroseismic intensities and peak measures of seismic motions (PGA and PGV) for the intermediate-depth earthquakes of the southern Aegean Sea. Input model parameters for the simulations, such as fault dimensions, stress parameters, and attenuation parameters (e.g. back-arc/along anelastic attenuation) are adopted from previous work performed in the area. Site-effects on the observed seismic motions are approximated using generic transfer functions proposed for the broader Aegean Sea area on the basis of V_S30 values from topographic slope proxies. The results are in very good agreement with the observed anomalous damage patterns, for which the largest intensities are often observed at distances >?100 km from the earthquake epicenters. We also consider two additional “prediction” but realistic intermediate-depth earthquake scenarios, and model their macroseismic distributions, to assess their expected damage impact in the broader southern Aegean area. The results suggest that intermediate-depth events, especially north of central Crete, have a prominent effect on a wide area of the outer Hellenic arc, with a very important impact on modern urban centers along northern Crete coasts (e.g. city of Heraklion), in excellent agreement with the available historical information.     

南北地震带两次典型强震周边自然电场变化及机理探讨

根据场源不同可将地电场E分为大地电场E_T和自然电场E_SP。空间S_q电流系和固体潮一般被认为是大地电场E_T的起源;自然电场E_SP源于地下介质的物理化学作用,其局部性变化相对稳定。2013年南北地震带相继发生了芦山M_S7.0地震（30.3°N,103.0°E）和岷县、漳县M_S6.6地震（34.5°N,104.2°E）。两次地震周边约400 km范围内的地电场台站（成都、汉王）数据表明震前自然电场E_SP出现明显的小幅度突跳或大幅度跃变,同时在时间上具有准同步性,其他地电场台站也有类似变化。基于大地电场岩体裂隙水（电荷）渗流（移动）模型,对两次地震周边自然电场E_SP的变异机理进行探讨,尝试解析其物理过程。     

Activation of seismicity in Central and South Asia after the Makran earthquakes: Possible acceleration of preparation of large seismic events in the Tien Shan region

Two zones of seismicity (ten events with M _w = 7.0–7.7) stretching from Makran and the Eastern Himalaya to the Central and EasternTien Shan, respectively, formed over 11 years after the great Makran earthquake of 1945 (M _w = 8.1). Two large earthquakes (M _w = 7.7) hit theMakran area in 2013. In addition, two zones of seismicity (M ≥ 5.0) occurred 1–2 years after theMakran earthquake in September 24, 2013, stretching in the north-northeastern and north-northwestern directions. Two large Nepal earthquakes struck the southern extremity of the “eastern” zone (April 25, 2015, M _w = 7.8 and May 12, 2015, M _w = 7.3), and the Pamir earthquake (December 7, 2015, M _w = 7.2) occurred near Sarez Lake eastw of the “western” zone. The available data indicate an increase in subhorizontal stresses in the region under study, which should accelerate the possible preparation of a series of large earthquakes, primarily in the area of the Central Tien Shan, between 70° and 79° E, where no large earthquakes (M _w ≥ 7.0) have occurred since 1992.     

Gravity variation before Kunlun mountain pass western <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 earthquake

The relation between the gravity variation features and M _S=8.1 earthquake in Qinghai-Xizang monitoring area is analyzed preliminarily, by using spatial dynamic variation results of regional gravity field from absolute gravity and relative gravity observation in 1998 and 2000. The results show that: 1) M _S=8.1 earthquake in Kulun mountain pass western occurred in the gravity variation high gradient near gravity’s high negative variation; 2) The main tectonic deformation and energy accumulation before M _S=8.1 earthquake are distributed at south side of the epicenter; 3) The range of gravity’s high negative variation at east of the M _S=8.1 earthquake epicenter relatively coincides with that rupture region according to field geology investigation; 4) Gravity variation distribution in high negative value region is just consistent with the second shear strain’s high value region of strain field obtained from GPS observation.     

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.     

Effect of zonal E × B plasma drift on electron density in the low-latitude ionospheric <Emphasis Type="Italic">F</Emphasis> region at a solar activity maximum near vernal equinox

The variations in the density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift perpendicularly to the magnetic (B) and electric (E) fields in the direction geomagnetic west-geomagnetic east have been studied using the three-dimensional nonstationary theoretical model of electron and ion densities (N _e and N _i ) and temperatures (T _e and T _i ) in the low-latitude and midlatitude ionospheric F region and plasmasphere. The method of numerical calculations of N _e , N _i , T _e , and T _i , including the advantages of the Lagrangian and Eulerian methods, is used in the model. A dipole approximation of the geomagnetic field (B), taking into account the non-coincidence of the geographic and geomagnetic poles and differences between the positions of the Earth’s and geomagnetic dipole centers, is accepted in the calculations. The calculated NmF2 and altitudes of the F2 layer maximum (hmF2) have been compared with these quantities measured at 16 low-latitude ionospheric sounding stations during the geomagnetically quiet period October 11–12, 1958. This comparison made it possible to correct the input model parameters: the NRLMSISE-00 model [O], the meridional component of the neutral wind velocity according to the HWW90 model, and the meridional component of the equatorial plasma drift due to the electric field specified by the empirical model. It has been indicated that the effect of the zonal E × B plasma drift on NmF2 can be neglected under daytime conditions and changes in NmF2 and hmF2 under the action of this drift are insignificant under nighttime conditions north of 25° and south of ?26° geomagnetic latitude. The effect of the zonal E × B plasma drift on NmF2 and hmF2 is most substantial in the nightside ionosphere approximately from ?20° to 20° geomagnetic latitude, and the neglect of this drift results in an up to 2.4-fold underestimation of NmF2. The found dependence of the effect of the zonal E × B plasma drift on NmF2 and hmF2 on geomagnetic latitude is related to the longitudinal asymmetry of B, asymmetry of the neutral wind about the geomagnetic equator, and changes in the meridional E × B plasma drift at a change in geomagnetic longitude.     

Moment tensor inversion and source rupture process of the September 27, 2003 <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=7.9 earthquake occurred in the border area of China,Russia and Mongolia

We conducted moment tensor inversion and studied source rupture process for M _S=7.9 earthquake occurred in the border area of China, Russia and Mongolia on September 27 2003, by using digital teleseismic P-wave seismograms recorded by long-period seismograph stations of the global seismic network. Considering the aftershock distribution and the tectonic settings around the epicentral area, we propose that the M _S=7.9 earthquake occurred on a fault plane with the strike of 127°, the dip of 79° and the rake of 171°. The rupture process inversion result of M _S=7.9 earthquake shows that the total rupture duration is about 37 s, the scalar moment tensor is M ₀=0.97×10²⁰ N·m. Rupture mainly occurred on the shallow area with 110 km long and 30 km wide, the location in which the rupture initiated is not where the main rupture took place, and the area with slip greater than 0.5 m basically lies within 35 km deep middle-crust under the earth surface. The maximum static slip is 3.6 m. There are two distinct areas with slip larger than 2.0 m. We noticed that when the rupture propagated towards northwest and closed to the area around the M _S=7.3 hypocenter, the slip decreased rapidly, which may indicate that the rupture process was stopped by barriers. The consistence of spatial distribution of slip on the fault plane with the distribution of aftershocks also supports that the rupture is a heterogeneous process owing to the presence of barriers.     

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.     

Estimation of maximum mass velocities in the focal zones of strong earthquakes based on measured displacements of rock blocks: A case study of Kyrgyzstan

An important task in seismic hazard assessment is estimation of the intensity and frequency of extremely strong earthquake effects, in particular, peak ground velocities (PGV). Earlier, a method was proposed to evaluate PGV values based on the magnitude of displacements of rock blocks (Rodkin et al., 2012). In this study, this method is used to analyze field data on the source zones of the August 19, 1992, M_S = 7.3 Susamyr earthquake and the January 3, 1911, M_w = 7.9 Kemin earthquake, and estimate maximum ground shaking at the upper construction site of the Upper Naryn series of hydropower plants, Kyrgyz Republic. It is shown that the resulting estimates are consistent with data obtained through other techniques. Therefore, the new approach can be recommended to estimate earthquake effects.     

A way of estimating the characteristic slip displacement

During the ruptures of an earthquake,the strain energy.△E,.will be transferred into,at least,three parts,i.e..the seismic radiation energy(E_s),fracture energy(E_g),and frictional energy(E_f),that is,△E = E_s + E_g + E_f.Friction,which is represented by a velocity- and state-dependent friction law by some researchers,controls the three parts.One of the main parameters of the law is the characteristic slip displacement.D_c.It is significant and necessary to evaluate the reliable value of D_c from observed and inverted seismic data.Since D_c controls the radiation efficiency.η_R = E_s/(E_s+ E_g),the value of η_r is a good constraint of estimating D_c.Integrating observed data and inverted results of source parameters from recorded seismograms.the values of E_s and E_g of an earthquake can be measured,thus leading to the value of η_R.The constraint used to estimate the reliable value of D_c will be described in this work.An example of estimates of D_c.based on the observed and inverted values of source parameters of the September 20,1999 M_S 7.6 Chi-Chi(Ji-Ji).Taiwan region,earthquake will be presented.     

The Tolud Burst of Seismicity and the Earthquake of November 30, 2012 (<Emphasis Type="Italic">M</Emphasis><Subscript>C</Subscript> = 5.4, <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> = 4.8) that Accompanied the Start of the 2012‒2013 Tolbachik Eruption

This paper reports a study of the Tolud earthquake sequence; the sequence was a burst of shallow seismicity between November 28 and December 7, 2012; it accompanied the initial phase in the Tolbachik Fissure Eruption of 2012?2013. The largest earthquake (the Tolud earthquake of November 30, 2012, to be referred to as the Tolud Earthquake in what follows, with K_S = 11.3, M_L = 4.9, M_C = 5.4, and M_W = 4.8) is one of the five larger seismic events that have been recorded at depths shallower than 10 km beneath the entire Klyuchevskoi Volcanic Cluster in 1961?2015. It was found that the Tolud earthquake sequence was the foreshock–aftershock process of the Tolud Earthquake. This is one of the larger seismicity episodes ever to have occurred in the volcanic areas of Kamchatka. Data of the Kamchatka seismic stations were used to compute some parameters for the Tolud Earthquake and its largest (M_L = 4.3) aftershock; the parameters include the source parameters and mechanisms, and the moment magnitudes, since no information on these is available at the world seismological data centers. The focal mechanisms for the Tolud Earthquake and for its aftershock are consistent with seismic ruptures at a tension fault in the rift zone. Instrumental data were used to estimate the intensity of shaking due to the Tolud Earthquake. We discuss the sequence of events that was a signature of the time-dependent seismic and volcanic activity that took place in the Tolbachik zone in late November 2012 and terminated in the Tolud burst of seismicity. Based on the current ideas of the tectonics and magma sources for the Tolbachik volcanic zone, we discuss possible causes of these earthquakes.     

Mechanism of formation of <Emphasis Type="Italic">Q</Emphasis>-disturbances in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> region of the equatorial ionosphere

Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F ₂ region of the equatorial ionosphere under quiet geomagnetic conditions (Q-disturbances) is interpreted. It is shown that the observed variations in the velocity of the vertical E × B plasma drift, related to the zonal E _y component of the electric field, are responsible for the formation of Q-disturbances. The plasma recombination at altitudes of the lower part of the F ₂ region and the dependence of the rate of this process on heliogeophysical conditions compose the mechanism of Q-disturbance formation at night. The daytime positive Q-disturbances are caused exclusively by a decrease in the upward E × B drift, and this type of disturbances could be related to the known phenomenon of counter electrojet. Possible causes of formation of the daytime negative Q-disturbances are discussed.     

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.     

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.