Multiple scattering of seismic waves in fractured media: Velocity and effective attenuation of the coherent components ofP waves andS waves

We show that the multiple scattering by small fractures of seismic waves with wavelengths long compared to the fracture size and fracture spacing is indistinguishable from multiple-scattering effects produced by regular porosity, except for an orientation factor due to fracture alignment. The fractures reduce theP-wave andS-wave velocities and produce an effective attenuation of the coherent component of the seismic waves. The attenuation corresponds to 1000/Q of about unity for a Gaussian spectrum of fractures, and it varies with frequencyf asf ³. For a Kolmogorov spectrum of fractures of spectral index the attenuation is an order of magnitude or so larger and varies with frequency asf ^3-v The precise degree of attenuation depends upon the matrix properties, the fracture porosity, the degree of fracture anisotropy, the type of fluid filling the fractures, and the incidence angle of the wave.For fracture porosities less than about 15% theP-wave andS-wave velocities are decreased by the order of 5–10% with a lesser dependence on the type of fluid filling the fractures (gas, oil, or brine) and with a dependence on both the degree of anisotropy and the incident angle made by the wave. The tendency of fractures to occur perpendicularly to bedding suggests that the best way to measure seismically fractured rock behavior in situ is by using the travel-time delay and reflection amplitude. As both the offset and the azimuth of receivers vary from a shot, the travel-time delay and reflection amplitude should both show an elliptical pattern of behavior—the travel-time delay in response to the varying seismic speed, and the reflection amplitude in response to angular variations in the multiple scattering. Observations of attenuation at several frequencies should permit (a) determination of the spectrum of fractures (Gaussian versus Kolmogorovian) and (b) determination of the contribution of viscous damping to the effective attenuation.     

Local seismic observations at stations Nový Kostel and Skalná, and their iterpretation: Period 1986–93

Summary Systematic observations of natural seismic activity in the West-Bohemian earthquake-swarm region began at two autonomous seismological stations of the Geophysical Institute, the digital station Novy Kostel (NKC) and the analogue station Skalná (SKC), in May 1986 and December 1989, respectively. This paper presents the station data of NKC and SKC, the method of processing the records, and the database structure. It also includes the interpretation of observations made at these stations in the period 1986–1993. It was found that the seismic activity in the West Bohemian region, in the Vogtland (Saxony) and NE Bavaria was continuous. Between two strong earthquake swarms, the energy in this area was released in the form of a large number of micro-earthquakes of a markedly swarm-like nature, mostly concentrated in six focal regions.     

Seismic regime of the west Bohemian earthquake swarm region: Preliminary results

Summary The western part of the Bohemian Massif located between two tectonic units, the Moldanubian and the Saxo-Thuringian, is characterized by the re-occurrence of earthquake swarms. The focal region for these swarms includes the territory of West Bohemia and the adjacent territory of SE Saxony and NE Bavaria. During the most recent swarm in December 1985 – January 1986, more than 8000 small earthquakes were recorded; the two largest earthquakes with local magnitudes (M_L) of 4·6 and 4·1. This paper presents a summary of the seismic energy release in space and time for the western part of the Bohemian Massif, based on seismic observations of permanent seismic stations established in West Bohemia since 1986. It was found out that microearthquake activity, mostly of a swarm-like character, persisted between two macroseismically observed swarms. The foci of the microearthquakes predominantly cluster in six main epicentral zones, four of which are located in West Bohemia or in its immediate vicinity in Saxony. The remaining two are in Saxony and in Bavaria. The four epicentral zones in West Bohemia were studied in detail. It was found that the individual zones differ in size, in depth of hypocentres, in geometry, as well as in temporal activity. Moreover, it was found that the seismicity in the most active epicentral zone is closely related to the system of principal tectonic faults referred to as the Kruné Hory fault and the Mariánské Lázn fault.     

CrustalP-wave velocity and velocity-ratio study in northeast India by a microearthquake survey

The upper crustal (20 km)P-wave velocity beneath the Shillong Plateau and Nowgong area has been studied by the time-distance plot method. TheP-arrival data of the shallow (20 km) microearthquakes from three temporary networks are used, and the average velocity is found to be 5.55 km/s. The velocity ratio (V _p /V _s ) for the upper crust (0–20 km) as well as for the lower crust (21–40 km) are determined by the Wadati-plot method and station-by-station method. The average value obtained by the two methods is compatible; theV _p /V _s ranges between 1.74 to 1.76. A generalized seismic velocity model of the area is suggested by this study, which has been very useful for microearthquake location.     

Intrinsic Qp seismic attenuation from the rise time of microearthquakes: a local scale application at Rio‐Antirrio,Western Greece

We determine the attenuation structure of a three‐dimensional medium based on first pulse‐width measurements from microearthquake data. Ninety‐five microearthquakes from a seventy stations local network were considered in this study. Measurements of the first half cycle of the wave, the so‐called rise time τ were carried out on high quality velocity seismograms and inverted to estimate the P‐waves intrinsic quality factor Q_p. The results of this investigation indicate that first pulse width data from a local microearthquake network permit retrieval with sufficient accuracy of the heterogeneous Q_p structure. The inferred attenuation variability corresponds to the known geological formations in the region.     

水口水库诱发地震的微震前兆信息与趋势预测

福建水口水库地震是我国当前正在发生的一个大型水库诱发地震实例。针对水库地震序列中普遍含有大量微震的特点,探讨微震活动可能包含的前兆信息。研究发现,“微震相对频度”在水口水库地震趋势预测中具有明显的前兆性质,文中还介绍了“微震相对频度”作为一种前兆指标在该水库地震趋势预测实践中的应用效果,对“微震相对频度”的物理意义以及深入开展“微地震研究”的前景进行了讨论。     

Site term from single-station sigma analysis of S-waves in western Turkey

The main aim of this study is to obtain site terms from single-station sigma analysis and to compare them with the site functions resulting from different techniques. The dataset consists of 1764 records from 322 micro- and moderate-size local earthquakes recorded by 29 stations in western Turkey. Median models were derived from S-wave Fourier amplitude spectra for selected 22 frequencies, by utilizing the MLR procedure which performs the maximum likelihood (ML) estimation of mixed models where the fixed effects are treated as random (R) effects with infinite variance. At this stage, b (geometrical spreading coefficient) and Q (quality factor) values were decomposed, simultaneously. The residuals of the median models were examined by utilizing the single-station sigma analysis to obtain the site terms of 29 stations. Sigma for the median models is about 0.422 log10 units and decreases to about 0.308, when the site terms from the single-station sigma analysis were considered (27% reduction). The event-corrected within-event standard deviations for each frequency are rather stable, in the range 0.19–0.23 log10 units with an average value of 0.20 (±?0.01). The site terms from single-station sigma analysis were compared with the site function estimates from the horizontal-to-vertical-spectral-ratio (HVSR) and generalized inversion (INV) techniques by Akyol et al. (2013) and Kurtulmu? and Akyol (2015), respectively. Consistency was observed between the single-station sigma site terms and the INV site transfer functions. The results imply that the single-station sigma analysis could separate the site terms with respect to the median models.     

A note on the stress-dilatancy relation for simulated fault gouge

Theoretical constraints on the stress-dilation relation for a deforming Coulomb material requirev ifC=0 andv sin^-1( _m / _m ) always, wherev is the dilation angle, is the friction angle,C is cohesion, _m is the maximum shear stress, and _m is the mean effective stress. Recent laboratory measurements of friction and dilatancy of simulated fault gouge show that small amplitude shear-load cycling causes compaction and consolidation. Comparison of the data with theory indicates that such load cycling produces: (1) increased coefficient of friction (or friction angle), (2) increased cohesion, and (3) increased dilatancy rate (or dilation angle). Under certain conditions of load cycling without significant plastic shear strain accumulation ( ^p <0.005) we find thatv exceeds both and, in contrast to theory, sin^-1( _m / _m ). This result is interpreted in terms of enhanced cohesion and overconsolidation, which lead to residual stresses within the gouge. An analogy is drawn between these special loading conditions and those extant on natural faults. In particular, our results imply that jostling and minor stress variations associated with microearthquakes may produce strengthening of fault gouge and changes in the fault zone's stress-dilatancy relation. Hence, compaction associated with microseismicity may lead to subsequent dilation of fault gouge, even for faults with large displacement rates and large net offsets (e.g., San Andreas). In regions where such dilation persists over sufficient displacements (on the order of the critical slip distance for seismic faulting) it may tend to inhibit unstable slip.     

Scenario of the January 1997 West Bohemia Earthquake Swarm

In order to learn more about the nature of the dynamic processes taking place in the West Bohemia/Vogtland earthquake swarm region, we investigated the temporal and spatial variations of the source mechanisms of the January 1997 swarm beneath Nový Kostel (NKC). Visual analyses of WEBNET seismograms of over 800 events revealed that a specific feature of this swarm was the occurrence of eight classes of multiplet events. The result of single-source, absolute moment tensor inversion of the P and SH peak amplitudes of a subset of 70 events representing all multiplet classes indicated that eight statistically significant types of mechanisms occurred during the swarm. Two of them, types A and B in our denotation, comprised all M _L 1.3 events and predominated in the swarm. Type A were pure strike-slip mechanisms or strike-slip mechanisms containing a small normal component, with a nearly pure double-couple source. For class B events, oblique-thrust faulting and non-double-couple components significant at a fairly high confidence level were typical. Type A events predominated in the southern subcluster of the swarm, whereas most of type B events occurred in the subcluster northwards from NKC. This indicates that two major seismogenic planes were active during the swarm. The swarm essentially developed in four phases: in the first, type A events prevailed and the southern plane was active; during the second, characterised by the occurrence of both type A and B events (the former in the southern, the latter predominantly in the northern subcluster), the activity of the swarm culminated; in the third and fourth, the occurrence of type B events in the northern plane predominated, and only weak single events occurred southwards from NKC. Mechanisms of types AB , C , D , E , F and G , which were typical for M _L 1.2 events, occurred randomly throughout the swarm. Type AB events were identified in both the southern and northern clusters, type C , E , F and G mechanisms only southwards from NKC. Type D events exhibited a large scatter of hypocentres which fell in neither the southern nor the northern cluster. Focal mechanisms like those reported in this study and with analogous temporal and spatial variations were observed by other authors already fifteen years ago in the 1985/86 earthquake swarm and may, therefore, be typical for the region under study.     

井下记录的PP和SS震相

本文在系统分析静海台井下地震记录图的基础上,发现并解决了具有一定实际和理论意义的两个问题:其一是证明了最初误认为是干扰的“双脉冲”图形实际上均是微地震。从而指出研究静海周围地区的地震活动性应考虑微震活动问题。其二确认静海周围地区微震波形的4个主要震相是P、PP、S和SS,并对地面反射波PP和SS能够观测到的条件作了初步分析,其结果进一步说明井下地震记录较之地面记录更为“逼真”,复杂和丰富。     

Amplitude Response of a LCR Gravimeter with Feedback at Periods of Microseisms and Earthquake Waves

The amplitude response of a LCR gravimeter with the SRW-E feedback was determined on a vertical vibrating platform. The ink-pen recorder was connected parallel with the digital voltmeter input to obtain an analog response of the gravimeter to the harmonic motion of the base with a peak-to-peak amplitude of 10 m and periods excited in an interval of 4 - 10 s as by the ground motion of meteorological microseisms, and in the interval 10 -100 s as by surface waves of distant strong earthquakes. In the first interval, an unexpected maximum of the amplitude response was observed with the double amplitude of apparent 6,5 Gal (6.5 × 10 _–8 m/s ² ) at a period of 4.8 s, and a baseline shift with the amplitude of–64Gal was observed at the same period. The value of this direct component cannot be separated from the effect of the Earth's gravity field. In the second interval, the amplitude response of the gravimeter displayed one expected maximum at a period of 40 s with the double amplitude of 8.1Gal. At the usual level of microseisms with a peak-to-peak amplitude of up to 2 m the fluctuation of the gravimeter on the direct baseline shifted by –4.9 Gal was estimated at 1 Gal. With typical Rayleigh surface waves with periods of 20 s and double amplitudes of up to 100 m, the fluctuation reached 67 Gal.     

华北地区中等强度地震前近震中区微震活化-寂静图象

本文研究了1970——1979年华北地区11个中等强度地震前区域微震时-空分布图象.发现大部分地震前1——2年,在其周围70——120公里范围内(简称震兆区),微震活动出现活化——寂静图象.同时,震前微震活动似有一个逐渐扩散到日趋收缩直至发震的过程.利用泊松分布检验活化、寂静阶段地震次数发生的概率,其结果是活化阶段绝大多数地震的概率小于5%.而寂静阶段只有近半数地震的概率小于5%.另外,发现大多数地震前震兆区的累加频度曲线的速率有一个加速——减小的过程,而外围区的速率基本保持稳定.最后,给出典型的异常图象模式并予以初步解释.      

Classification of seismic events using linear discriminant function (LDF) in the Sakarya region,Turkey

The Sakarya prefecture is an interesting area with various seismicity types. This activity comes from earthquakes occurring at the North Anatolian Fault Zone and from a few quarry blast areas in the region. These quarry blast recordings produce errors in the determination of active faults and mapping of the microearthquake activity. Therefore, to recognize the tectonic activity in the region, we need to be able to discriminate between earthquakes and quarry blasts in the catalogues. In this study, a statistical analysis method (linear discriminant function) has been applied to classify seismic events occurring in the Sakarya region. We used 110 seismic events that were recorded by Sakarya University Seismic Station between 2012 and 2014. Time and frequency variant parameters, maximum S wave and maximum P wave amplitude ratio (S/P), the spectral ratio (Sr), maximum frequency (f_max), and total signal duration of the waveform were used for discrimination analyses. The maximum frequency (f_max) versus time duration of the seismic signal gives a higher classification percentage (94%) than the other discriminants. At the end of this study, 41 out of 110 events (44%) are determined as quarry blasts, and 62 (56%) are considered as earthquakes.     

Spatio-temporal variations of seismicity in the southern Peru and northern Chile seismic gaps

The spatio-temporal variation of seismicity in the southern Peru and northern Chile seismic gaps is analyzed with teleseismic data (m _b 5.5) between 1965 and 1991, to investigate whether these gaps present the precursory combination of compressional outer-rise and tensional downdip events observed in other subduction zones. In the outer-rise and the inner-trench (0 to 100 km distance from the trench) region, lower magnitude (5.0m _b <5.5) events were also studied. The results obtained show that the gaps in southern Peru and northern Chile do not present compressional outer-rise events. However, both gaps show a continuous, tensional downdip seismicity. For both regions, the change from compressional to tensional regime along the slab occurs at a distance of about 160 km from the trench, apparently associated with the coupled-uncoupled transition of the interplate contact zone. In southern Peru, an increase of compressional seismicity near the interplate zone and of tensional events (5.0m _b 6.3) in the outer-rise and inner-trench regions is observed between 1987 and 1991. A similar distribution of seismicity in the outer-rise and inner-trench regions is observed with earthquakes (m _b <5.5). In northern Chile there is a relative absence of compressional activity (m _b 5.5) near the interplate contact since the sequence of December 21, 1967. After that, only a cluster of low-magnitude compressional events has been located in the area 50 to 100 km from the trench. The compressional activity occurring near the interplate zone in both seismic gaps represents that a seismic preslip is occurring in and near the plate contact. Therefore, if this seismic preslip is associated with the maturity of the gap, the fact that it is larger in southern Peru than in northern Chile may reflect that the former gap is more mature than the latter. However, the more intense downdip tensional activity and the absence of compressional seismicity near the contact zone observed in northern Chile, may also be interpreted as evidence that northern Chile is seismically more mature than southern Peru. Therefore, the observed differences in the distribution of stresses and seismicity analyzed under simple models of stress accumulation and transfer in coupled subduction zones are not sufficient to assess the degree of maturity of a seismic gap.     

Fourier spectrum of the seiches of an ideal rectangular basin with sinusoidally varying level

Summary If the level of an ideal rectangular water basin varies sinusoidally with amplitude 2h ₁ and angular frequency ₀ the modes of the free oscillations of the basin split into multiplets with spacing ₀: the shape and the height of the peaks are given as functions of the mean levelh ₀ and of the lengthl of the basin and of the amplitude 2h ₁.Istituto di Geodesia e Geofisica, Universita' di Trieste, Italy.- This research has been supported by E.R.O.-U.S. Army, under contract DA-91-591-3538.     

Microearthquake observations around Thien Dam in northwest Himalayas

A microearthquake survey has been conducted around Thien Dam (32°26.5′N, 75°43.5′E) during March–April 1983 by operating seven high gain stations. The data recorded at these stations have been supplemented by other seismological stations in the region under the Beas and Salal Projects. Surficial velocity structure has been obtained from waves generated from blasts in the dam area and hypocentres were determined using the program HYPO 71.Microearthquake distribution shows four main zones of active seismicity. Deeper microearthquakes with focal depths greater than 10 km below the dam were attributed to the Satlitta fault. Composite focal mechanism solution shows that the Satlitta fault at this depth is of normal type. Shallow microearthquakes below the dam are due to the Suruin-Mastgarh anticline which is of the thrust type. Microearthquakes located northeast of the dam are apparently associated with a ‘back’ thrust.     

库尔勒5.8级地震前地震波特性的某些短期变化

本文对新疆1978年4月22日库尔勒5.8级地震前地震波的动力学特性和运动学特性的某些短期变化进行了初步研究。重点讨论了5.8级地震前通过孕震区的小地震的直达波振幅比、振动持续时间比以及波速比的短期异常变化,並且,结合震前短期地震活动异常图象作了对比分析。实际观测结果表明,较大地震发生前,通过孕震区的小震的地震波特性可能出现一些较明显的短期异常变化,因而,对此进行研究无疑对于探讨地震孕育过程和地震的预报途径都是有意义的。     

AbnormalSq behaviour

It is well-known that the amplitude and phase of theSq(H) variation show considerable variability from day to day. In this paper we consider one aspect of the phase variability—that associated with AQDs. AQDs (or abnormal quiet days) are defined as magnetically quiet days where the maximum excursion ofH at a mid-latitude station on the poleward side of the focus occurs outside the normal time range 0830–1330 LST. Such days exhibit properties, many of which appear quite distinct from the properties of the normalSq(H) variation. The properties of AQDs, and the proposals that have been made to explain them, are considered in detail. The consequences of these proposals and some problems which need to be addressed in order to obtain a fuller understanding of the dynamics of the ionosphere on AQDs are also discussed.     

Seismic source-parameter relations for earthquakes in Greece

A data set of nineteen, mainly shallow, moderate to large earthquakes, which occurred in the Aegean and the surrounding area, has been used to derive empirical relations for kinematic fault parameters. Thus the relations between seismic momentM ₀ and magnitudeM _s andm _b and betweenM ₀ andM _s and fault dimensionsS andL have been determined. From these relations and theoretical ones it was deduced that earthquakes in the Aegean and the surrounding events, chiefly interplate, are characterized by low average stress drop values. Values of ranging from 1 to 30 bar are consistent with the data. It was also found that, in general terms, most of the data obey the geometrical similarity conditionL=2w, whereL is the fault length measured along the strike andw is its width measured along the dip. For strike-slip faults, however, the conditionL=4w seems to hold.     

On the propagation of finite Love Waves

Summary The problem of the propagation of finite Love Waves in a heterogeneous elastic half space lying over a homogeneous elastic half space, using the quasilinear stress-strain relation due toS. Ferhst [4] is considered in detail. The variations of the parameter in the layer assumed to be of the form 1= 0e ^z, 0e ^z where is a constant andz is distance measured from the surface into the layer.