Rushan earthquake swarm in eastern China and its indications of fluid-triggered rupture

An extraordinary earthquake swarm occurred at Rushan on the Jiaodong Peninsula from October 1, 2013, onwards, and more than 12,000 aftershocks had been detected by December 31, 2015. All the activities of the whole swarm were recorded at the nearest station, RSH, which is located about 12 km from the epicenter. We examine the statistical characteristics of the Rushan swarm in this paper using RSH station data to assess the arrival time difference, \(t_{{{\text{S}} \,-\, {\text{P}}}}\), of Pg and Sg phases. A temporary network comprising 18 seismometers was set up on May 6, 2014, within the area of the epicenter; based on the data from this network and use of the double difference method, we determine precise hypocenter locations. As the distribution of relocated sources reveals migration of seismic activity, we applied the mean-shift cluster method to perform clustering analysis on relocated catalogs. The results of this study show that there were at least 16 clusters of seismic activities between May 6, 2014, and June 30, 2014, and that each was characterized by a hypocenter spreading process. We estimated the hydraulic diffusivity, D, of each cluster using envelope curve fitting; the results show that D values range between 1.2 and 3.5 m²/d and that approximate values for clusters on the edge of the source area are lower than those within the central area. We utilize an epidemic-type aftershock sequence (ETAS) model to separate external triggered events from self-excited aftershocks within the Rushan swarm. The estimated parameters for this model suggest that α = 1.156, equivalent to sequences induced by fluid-injection, and that the forcing rate (μ) implies just 0.15 events per day. These estimates indicate that around 3% of the events within the swarm were externally triggered. The fact that variation in μ is synchronous with swarm activity implies that pulses in fluid pressure likely drove this series of earthquakes.     

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.     

Location and source mechanism of the Karlsruhe earthquake of 24 September 2014

On 24 September 2014, a M_L 2.3 earthquake occurred southwest of the urban area of Karlsruhe, Germany, which was felt by a few people (maximum intensity I ₀?=?III). It was the first seismic event in this highly populated area since an I ₀?=?VII earthquake in 1948. Data of 35 permanent and temporary seismometers were analysed to localise the event and to determine the focal mechanism to compare it to previous seismicity. Restricting the data to P- and S-phases from 18 nearby stations and optimising the local earth model result in an epicentre in the southwest of the city at 48.986°N/8.302°E and in a hypocentral depth of 10 km. To calculate the focal mechanism, 22 P- and 5 S_H-polarities were determined that constrain a stable left lateral strike-slip focal mechanism with a minor thrusting component and nodal planes striking NE-SW and NW-SE. The epicentre lies in the vicinity of the I ₀?=?VII earthquake of 1948. Both events are part of the graben-parallel flower structure beneath the Upper Rhine Graben, parallel to the active Rastatt source zone, which runs 5 km further east and included the epicentre of the 1933 Rastatt I ₀?=?VII earthquake. The focal mechanisms of the 2014 and 1948 earthquakes show NE-SW striking nodal planes that dip to the southeast. However, for the 1948 event, a normal faulting mechanism was determined earlier. Taking the uncertainty of the epicentre and focal mechanism in 1948 and its fault dimensions into account, both events might have happened on the same fault plane.     

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.     

Seismicity study in Pournari reservoir area (W. Greece) 1981–2010

The spatial-temporal evolution of seismicity is examined, during the initial impoundment of Pournari reservoir located on Arachthos River (Western Greece), as well as for the next 30 years. The results show that, despite the relatively moderate-to-high seismicity from west to east, there is no remarkable earthquake in the vicinity before the first reservoir impoundment. Immediately after the impoundment (January 1981), and during the first 4 months, a considerable number of low-magnitude seismic events were recorded in the broader area of the dam. Moreover, two independent major events occurred on March 10, 1981 (M _L ?=?5.6) and April 10, 1981 (M _L ?=?4.7) with focal depths 13 and 10 km, respectively. The detailed analysis of the two corresponding aftershock sequences shows that they present different behaviors (e.g., larger b-value and lower magnitude of the main aftershock) than that of other aftershock sequences in Greece. This seismicity is probably due to triggering, via the water loading mechanism and the undrained response due to a flysch appearance on the reservoir basement. The activation of the thrust fault may be attributed to the bulging of evaporites that characterize the disordered structure of W. Greece, via possible water intake. The detailed processing of the recorded seismicity during the period 1982–2010, in comparison with the variations of Pournari Dam water level, shows an increase of shallow seismicity (h?≤?5 km) in the vicinity of the reservoir up to a 10-km distance—in contrast to the initial period, characterized by a number of deeper events due to the background response change from undrained to drained status.     

Mechanism of the M<Subscript><Emphasis Type="Italic">L</Emphasis></Subscript>4.0 25 April 2016 earthquake in southwest of France in the vicinity of the Lacq gas field

The source mechanism of the M_L 4.0 25 April 2016 Lacq earthquake (Aquitaine Basin, South-West France) is analyzed from the available public data and discussed with respect to the geometry of the nearby Lacq gas field. It is one of the biggest earthquakes in the area in the past few decades of gas extraction and the biggest after the end of gas exploitation in 2013. The routinely obtained location shows its hypocenter position inside the gas reservoir. We first analyze its focal mechanism through regional broad-band seismograms recorded in a radius of about 50 km epicentral distances and obtain EW running normal faulting above the reservoir. While the solution is stable using regional data only, we observe a large discrepancy between the recorded data on nearby station URDF and the forward modeling up to 1 Hz. We then look for the best epicenter position through performing wave propagation simulations and constraining the potential source area by the peak ground velocity (PGV). The resulting epicentral position is a few to several km away to the north or south direction with respect to station URDF such that the simulated particle motions are consistent with the observation. The initial motion of the seismograms shows that the epicenter position in the north from URDF is preferable, indicating the north-east of the Lacq reservoir. This study is an application of full waveform simulations and characterization of near-field ground motion in terms of an engineering factor such as PGV. The finally obtained solution gives a moment magnitude of M_w 3.9 and the best focal depth of 4 km, which corresponds to the crust above the reservoir rather than its interior. This position is consistent with the tendency of Coulomb stress change due to a compaction at 5 km depth in the crust. Therefore, this earthquake can be interpreted as a relaxation of the shallow crust due to a deeper gas reservoir compaction so that the occurrence of similar events cannot be excluded in the near future. It would be necessary to continue monitoring such local induced seismicity in order to better understand the reservoir/overburden behavior and better assess the local seismic hazard even after the end of gas exploitation.     

Gravity anomaly and crustal density structure in Jilantai rift zone and its adjacent region

This paper deals with the interpretation of Bouguer gravity anomalies measured along a 250 km long Suhaitu-Etuokeqi gravity profile located at the transitional zone of the Alxa and Ordos blocks where geophysical characteristics are very complex. The analysis is carried out in terms of the ratio of elevation and Bouguer gravity anomaly, the normalized full gradient of a section of the Bouguer gravity anomaly (G ^h ) and the crustal density structure reveal that (1) the ratio of highs and lows of elevation and Bouguer gravity anomaly is large between Zhengyiguan fault (F4) and Helandonglu fault (F6), which can be explained due to crustal inhomogeneities related to the uplift of the Qinghai-Tibet block in the northeast; (2) the main active faults correspond to the G ^h contour strip or cut the local region, and generally show strong deformation characteristics, for example the Bayanwulashan mountain front fault (F1) or the southeast boundary of Alxa block is in accord with the western change belt of G ^h , a belt about 10 km wide that extends to about 30 km; (3) Yinchuan-Pingluo fault (F8) is the seismogenic structure of the Pingluo M earthquake, and its focal depth is about 15 km; (4) the Moho depth trend and Bouguer gravity anomaly variation indicates that the regional gravity field is strongly correlated with the Moho discontinuity.     

The 2014–2015 earthquake series in the northern Upper Rhine Graben,Central Europe

Since March 2014, an unusually large amount of earthquakes occur southeast of the city of Darmstadt in the northern Upper Rhine Graben. During the period, until April 2015, we have recorded 356 earthquakes with magnitudes ranging from M_L?=??0.6 to 4.2. We identified two source clusters separated laterally by about 5 km. The hypocentres within these clusters are aligned vertically extending over a depth range from 1 to 8 km with a lateral extent of about 1 to 2 km. Focal mechanisms show left-lateral strike-slip movements; b values are changing with time between b?=?0.6 and b?=?0.9. This is the first time in almost 150 years that such high earthquake rates have been observed in the region. Historical accounts dating back to the nineteenth century report of over 2000 felt earthquakes over a time span from 1869 to 1871. From these, maximum intensities of VII have been estimated. Other seismic activities in the region were reported in the 1970s. The observations of the 2014–2015 earthquake series do not completely match a typical main shock–aftershock sequence or a typical earthquake swarm. Especially the activity at the beginning of the earthquake series may be considered as a mixture of a main shock–aftershock sequence and a short-lasting swarm event. Whether or not the time gap between the current seismic activity, which actually takes place at the same locations as parts of the seismic swarm in 1869–1871, and the seismic activity in the nineteenth century or the seismic activity in the 1970s can be interpreted as a seismic cycle remains unclear.     

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.     

An Evaluation of Coulomb Stress Changes from Earthquake Productivity Variations in the Western Gulf of Corinth,Greece

Spatial and temporal evolution of the stress field in the seismically active and well-monitored area of the western Gulf of Corinth, Greece, is investigated. The highly accurate and vast regional catalogues were used for inverting seismicity rate changes into stress variation using a rate/state-dependent friction model. After explicitly determining the physical quantities incorporated in the model (characteristic relaxation time, fault constitutive parameters, and reference seismicity rates), we looked for stress changes across space and over time and their possible association with earthquake clustering and fault interactions. We focused our attention on the Efpalio doublet of January 2010 (M = 5.5 and M = 5.4), with a high aftershock productivity, and attempted to reproduce and interpret stress changes prior to and after the initiation of this seismicity burst. The spatial distribution of stress changes was evaluated after smoothing the seismological data by means of a probability density function (PDF). The inverted stress calculations were compared with the calculations derived from an independent approach (elastic dislocation model) and this comparison was quantified. The results of the two methods are in good agreement (up to 80 %) in the far field, with the inversion technique providing more robust results in the near field, where they are more sensitive to the uncertainties of coseismic slip distribution. It is worth mentioning that the stress inversion model proved to be a very sensitive stress meter, able to detect even small stress changes correlated with spatio–temporal earthquake clustering. Data analysis was attempted from 1975 onwards to simulate the stress changes associated with stronger earthquakes over a longer time span. This approach revealed that only M > 5.5 events induce considerable stress variations, although in some cases there was no evidence for such stress changes even after an M > 5.5 earthquake.     

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.     

Source parameters of seismic events potentially associated with damage in block 33/34 of the Kiirunavaara mine (Sweden)

Forty-six mining-induced seismic events with moment magnitude between ?1.2 and 2.1 that possibly caused damage were studied. The events occurred between 2008 and 2013 at mining level 850–1350 m in the Kiirunavaara Mine (Sweden). Hypocenter locations were refined using from 6 to 130 sensors at distances of up to 1400 m. The source parameters of the events were re-estimated using spectral analysis with a standard Brune model (slope ?2). The radiated energy for the studied events varied from 4.7 × 10^?1 to 3.8 × 10⁷ J, the source radii from 4 to 110 m, the apparent stress from 6.2 × 10² to 1.1 × 10⁶ Pa, energy ratio (E _s/E _p) from 1.2 to 126, and apparent volume from 1.8 × 10³ to 1.1 × 10⁷ m³. 90% of the events were located in the footwall, close to the ore contact. The events were classified as shear/fault slip (FS) or non-shear (NS) based on the E _s/E _p ratio (>10 or <10). Out of 46 events 15 events were classified as NS located almost in the whole range between 840 and 1360 m, including many events below the production. The rest 31 FS events were concentrated mostly around the production levels and slightly below them. The relationships between some source parameters and seismic moment/moment magnitude showed dependence on the type of the source mechanism. The energy and the apparent stress were found to be three times larger for FS events than for NS events.     

Joint inversion of GNSS and teleseismic data for the rupture process of the 2017 <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript>6.5 Jiuzhaigou,China, earthquake

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 × 10¹⁸ N·m (M_w?≈?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 M_w7.9 Wenchuan, China, earthquake. Therefore, the Wenchuan earthquake promoted the occurrence of the 8 August 2017 Jiuzhaigou earthquake.     

Fault slip and identification of the second fault plane in the Varzeghan earthquake doublet

An intraplate earthquake doublet, with 11-min delay between the events, devastated the city of Varzeghan in northwestern Iran on August 11, 2012. The first Mw 6.5 strike-slip earthquake, which occurred after more than 200 years of low seismicity, was followed by an Mw 6.4 oblique thrust event at an epicentral separation of about 6 km. While the first event can be associated with a distinct surface rupture, the absence of a surface fault trace and no clear aftershock signature makes it challenging to identify the fault plane of the second event. We use teleseismic body wave inversion to deduce the slip distribution in the first event. Using both P and SH waves stabilize the inversion and we further constrain the result with the surface rupture extent and the aftershock distribution. The obtained slip pattern shows two distinct slip patches with dissimilar slip directions where aftershocks avoid high-slip areas. Using the estimated slip for the first event, we calculate the induced Coulomb stress change on the nodal planes of the second event and find a preference for higher Coulomb stress on the N-S nodal plane. Assuming a simple slip model for the second event, we estimate the combined Coulomb stress changes from the two events on the focal planes of the largest aftershocks. We find that 90% of the aftershocks show increased Coulomb stress on one of their nodal planes when the N-S plane of the second event is assumed to be the correct fault plane.     

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.     

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.     

Simulation of spatial and temporal properties of aftershocks by means of the fiber bundle model

The rupture processes of any heterogeneous material constitute a complex physical problem. Earthquake aftershocks show temporal and spatial behaviors which are consequence of the heterogeneous stress distribution and multiple rupturing following the main shock. This process is difficult to model deterministically due to the number of parameters and physical conditions, which are largely unknown. In order to shed light on the minimum requirements for the generation of aftershock clusters, in this study, we perform a simulation of the main features of such a complex process by means of a fiber bundle (FB) type model. The FB model has been widely used to analyze the fracture process in heterogeneous materials. It is a simple but powerful tool that allows modeling the main characteristics of a medium such as the brittle shallow crust of the earth. In this work, we incorporate spatial properties, such as the Coulomb stress change pattern, which help simulate observed characteristics of aftershock sequences. In particular, we introduce a parameter (P) that controls the probability of spatial distribution of initial loads. Also, we use a “conservation” parameter (π), which accounts for the load dissipation of the system, and demonstrate its influence on the simulated spatio-temporal patterns. Based on numerical results, we find that P has to be in the range 0.06 < P < 0.30, whilst π needs to be limited by a very narrow range (0.60 < π < 0.66) in order to reproduce aftershocks pattern characteristics which resemble those of observed sequences. This means that the system requires a small difference in the spatial distribution of initial stress, and a very particular fraction of load transfer in order to generate realistic aftershocks.     

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.     

Upper mantle structure beneath central Western Europe from data on phase velocities of surface waves

The deep structure of the upper mantle is determined from data on phase velocities of Love and Rayleigh waves measured by a differential method on traces between two stations in central Western Europe. One-dimensional velocity structures are first constructed from data of each pair of stations, after which two-dimensional distributions of SH and SV velocities are calculated by the method of two-dimensional tomography from S wave velocities at fixed depths. The results are presented in the form of 2-D vertical structures of the average S wave velocity (S = (SV + SH)/2) constructed along profiles crossing the region in directions of the best resolution. The main structural features are a higher velocity zone at depths of 60–80 km in the area (48°–50°N, 9°–11°E) and a lower velocity zone in the western part of the region at depths of 100–150 km, probably extending farther beyond the studied area.     

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.