Frequency characteristics of coda and determination of medium parameters in Lingwu region, China

In this paper, considering the influences of source spectrum, the scattering property of medium and instrument response on the dominant frequency of coda, a method of using the coda of local earthquake to determine the correlation length of medium andQ-value is given. We find the following formula as: {fx719-1} wheret* =t/Q, f is the dominant frequency of coda,u ₁ andu ₂ are the parameters depend on the correlation length and the corner frequency of the source spectrum respectively,I(f) is a function of instrument response. If the source parameter is given, we can obtain the correlation length andQ-value by means of the inversion of observed curves off-t of coda. We processed the data of coda wave of more than 40 earthquakes from 1982 to 1989 in Lingwu region, China, determined the correlation length andQ-value, and preliminarily studied the temporal change of correlation length before and after moderately strong earthquakes. The result suggests that there are indications that the correlation length of medium decreases before the moderate earthquake. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 62–70, 1992.     

The physical implication of coda amplitude ratio and its possible application

The physical implication of coda amplitude ratio is discussed in term of energy ratio. The digitized data recorded at the station of Beijing Telemetered Seismograph Network between 1989 and 1990 are used to calculate amplitude ratios of coda to direct S wave, and energy ratios. The spectral energy ratios are used to estimate the coda Q and mean free path l in the Beijing area, as well as the two quality factors Q _i and Q _S separately due to intrinsic absorption and scattering attenuation. The decay of seismic waves in their propagation seems mainly resulted from the intrinsic absorption in Beijing region. The temporal variations of amplitude ratio and energy ratio at Changli station during the above two years are inspected; some of them largely depart from their mean value. It may reflect the seismogenic process, but using the data lasting longer time with more case histories needs further study. This study is sponsored by the Key Project of State Science and Technology of China, No. 96-918.     

Study of the Change in Wave Velocity Ratio before and after Two Strong Earthquakes Using Airgun Source Data

Using the signals excited by the large-volume airgun source at the Binchuan transmitting seismic station from January to June, 2016, arrival-time data was acquired at four stations near the epicenter of the Eryuan M_S4.5 and M_S4.0 earthquakes on February 8, 2016, as well as the epicenter of the Yunlong M_S5.0 and Eryuan M_S4.6 earthquakes on May 18, 2016 through the waveform cross-correlation technique. The wave velocity ratio of the four stations was calculated using the single-station method. At the same time, the b-value and the focal mechanism consistency parameters of the study area were also calculated. The results show that:(1) the wave velocity ratio of each station experienced a process of decline-recovery-fast rise before the two strong earthquakes, and a significant quasi-synchronous rapid rise occurred within 3-12 days before the earthquake; (2) the timing of the rapid rise of the wave velocity ratio of the four stations before the Yunlong M_S5.0 and Eryuan M_S4.6 earthquakes were related to the epicentral distance. The station which observed the earliest increase in rapid rise is the farthest one from the epicenter, and the station where the rapid rise appeared in the latest is closest to the epicenter; (3) the form of change of the wave velocity ratio before the earthquake was different between stations located at different directions in the epicenter area; (4) the b-value and the focal mechanism consistency parameter which is commonly used to characterize the stress level both showed a downward trend before the two strong earthquakes, and were consistent with the change in the wave velocity ratio. According to the preliminary analysis, the wave velocity ratio obtained by using airgun source can better reflect the change in the stress state of the underground medium.     

中强地震前地震波参数异常的研究

本文研究了1988年5月26日库尔勒M_35.5地震和1991年2月25日柯坪M_s6.5地震前后的地震波初动符号、振幅比、尾波持续时间比和尾波衰减系数的变化特征,发现两次地震前上述参数均具有不同程度的异常。地震前(?)波初动一致性越强,持续时间越长,未来地震震级越大;振幅比突跳是即将发生地震的明显信号;首先出现尾波衰减系数异常的台站距未来地震震中较近。根据上述参数的异常特征,成功地预报了1991年6月6日和静M_s5.2地震。     

The <Emphasis Type="Italic">Q</Emphasis>-factor estimates for the crust and upper mantle in the vicinity of Sochi and Anapa (North Caucasus)

The regularities in the radiation and propagation of seismic waves in the regions of the North Caucasus are analyzed for estimating the ground motion parameters during the probable future strong earthquakes. Based on the records of the regional earthquakes with magnitudes M_W ~ 3.9–5.6 within epicentral distances up to ~300 km obtained during the period of digital measurements at the Sochi and Anapa seismic stations, the Q-factors in the vicinities of these sites are estimated at ~55 f^0.9 and ~90f^0.7, respectively. The estimates were obtained by the coda normalization method developed by Aki, Rautian, and other authors. This method is based on the phenomenon of suppression of the earthquake (source) effects and local (site) responses by coda waves in the S-wave spectra. The obtained Q-factor estimates can be used for forecasting the ground shaking parameters for the future probable strong earthquakes in the North Caucasus in the vicinities of Sochi and Anapa.     

Frequency dependent characteristics of coda wave quality factor in central and southcentral Alaska

Coda wave quality factor (Q _c ) was investigated by using digital data (100 sample sec^–1) recorded by a vertical component short-period station installed for this study. The station was located in the greater Fairbanks area in central Alaska. From several hundred earthquakes recorded by this station in about a year, 27 earthquakes were selected for the above study; 7 of these selected earthquakes were located along the Alaska Wadati-Benioff zone (Pacific plate). The other 20 earthquakes were located in the area of intraplate seismicity (North American plate). The data was filtered using 9 pass-bands with center frequency varying from 1.5 Hz to 16 Hz with octave bandwidth. The values ofQ _c obtained from the coda amplitude decay rates measured on the filtered data after corrections due to the recording instrument and source-receiver separation show appreciable frequency dependence. The value ofQ _c varies in the range of 253 and 1190 corresponding to the frequency interval from 1.5 Hz to 16 Hz for the study area. This variation is close to that reported by others for the Kanto region of Japan. Moreover, the characteristics ofQ _c obtained in the present case seem to be independent of epicentral distance and hypocentral depth.     

Changes in the early part of the seismic coda due to localized scatterers: The estimation ofQ in a stope environment

A single scattering model was used to analyse the temporary changes in the mean density of scattered waves in a discrete random medium. The model of the mean energy density, originally proposed bySato (1977) for spherical radiation and isotropic scattering, has been modified and applied to a medium in which the scatterers are confined to a specified volume. The time variation of the early part of the mean energy density function for the different source durations was investigated. The dominant effect on the theoretical mean energy density is caused by the specified volume containing scatterers. The duration of the source pulse influences the early part of the coda fort/t ₀<1.2, wheret is the lapse time measured from the source origin time, andt ₀is arrival time of the body wave.The analysis of the coda signal of micro-events occurring immediately in front of the face enables us to estimate the size of the fracture zone induced by the stope. The model of the mean energy density of coda for a medium containing scatterers close to the seismic source was used to analyse a large number of events recorded close to an advancing mine face in a deep level gold mine in South Africa. The coda decay rate has two trends: the first, with a steep decay of coda, is produced by a larger deviation of rock parameters and/or larger size of the scatterers; the second trend, which decays more slowly, has the corresponding mean-free path ranging from 20 m to 200 m. The analysis indicates that the rock mass about 15–20 m from the stope contains a large proportion of fractured and blocked rock, which is the source of scattering. The scattering of theS-wave was much stronger and more stable, with the mean-free path varying from 11 m to 45 m. This is due to the shorter wavelength of theS wave in comparison with theP wave. The quality factor for theP coda wave varies from 30 to 100 in the fracture zone of stope and outside this zone it has a value of 300. The quality factor of theS wave varies from 20 to 78 in the equivalent volume. For rock surrounding the stope the ratioQ _sp ^–1 /Q _ss ^–1 varied from 0.31 to 0.69. This suggests that the radii of scatterers are smaller than 3.5 m.     

Source Parameters for 1999 North Anatolian Fault Zone Aftershocks

We develop a data set of aftershock recordings of the 1999, M = 7.4 Izmit and M = 7.2 Duzce (Turkey) earthquakes to study their source parameters. We combined seismograms from 44 stations maintained by several sources (organizations) to obtain a unified data set of events (2.1 ≤ M_w ≤ 5.5). We calculate source parameters of these small earthquakes by two methods that use different techniques to address the difficulty in obtaining source spectra for small earthquakes subject to interference from site response. One method (program NetMoment (NM), Hutchings, 2004) uses spectra of direct S waves in a simultaneous inversion of local high-frequency network data to estimate seismic moment, source corner frequency (f_c), site attenuation (k) and whole-path Q. This approach takes advantage of the source commonality in all recordings for a particular earthquake by fitting a common Brune source spectrum to the data with a and individual k. The second approach (Mayeda et al., 2003) uses the coda method (CM) to obtain “nonmodel-based” source spectra and moment estimates from selected broadband recording sites. We found that both methods do well for events that allow the comparison with seismic moment estimates derived from waveform modeling. Also, source spectra obtained from the two methods are very closely matched for most of the events they have in common. We use an F test to examine the trade-off between k and f_c picks identified by the direct S-wave method. About half of the events could be constrained to have less than a 50% average uncertainty in f_c and k. We used these source spectra solutions to calculate energy and apparent stress and compare these to estimates from the selected “good quality” source spectra from CM. Both studies have values mutually consistent and show a similar increase in apparent stress with increasing moment. This result has added merit due to the independent approaches to calculate apparent stress. We conclude that both methods are at least partially validated by our study, and they both have usefulness for different circumstances of recording local small earthquakes. CM would work well in studies for which there is a broad magnitude range of events and NM works well for local events recorded by band-limited recorders.     

Rediscovering Signal Complexity as a Teleseismic Discriminant

We re-examine the utility of teleseismic seismic complexity discriminants in a multivariate setting using United Kingdom array data. We measure a complexity discriminant taken on array beams by simply taking the logarithm of the ratio of the P-wave coda signal to that of the first arriving direct P wave (β_CF). The single station complexity discriminant shows marginal performance with shallow earthquakes having more complex signatures than those from explosions or deep earthquakes. Inclusion of secondary phases in the coda window can also degrade performance. However, performance improves markedly when two-station complexity discriminants are formed showing false alarm rates similar to those observed for network m_b − M_s. This suggests that multistation complexity discriminants may ameliorate some of the problems associated with m_b – M_s discrimination at lower magnitudes. Additionally, when complexity discriminants are combined with m_b – M_s there is a tendency for explosions, shallow earthquakes and deep earthquakes to form three distinct populations. Thus, complexity discriminants may follow a logic that is similar to m_b – M_s in terms of the separation of shallow earthquakes from nuclear explosions, although the underlying physics of the two discriminants is significantly different.     

Q c of three component seismograms of volcanic microearthquakes at Campi Flegrei volcanic area — Southern Italy

Digital recordings of three component microearthquake codas from shallow seismic events in the volcanic region of Campi Flegrei — Southern Italy — were used with an automatic technique to calculate the attenuation factorQ _c (codaQ) in the hypothesis of singleS toS backscattering.Results show the same value ofQ for each of the three components. This result is interpreted as due to isotropicS wave radiation pattern.A check of the coda method was performed using a single station method based on simple assumptions on the direct SH wave spectrum. Single stationQ was averaged over the stations and over the earthquakes. Results show that the two methods lead to comparable results.A frequency dependence quite different from that evaluated in active tectonic regions was found for coda attenuation, comparable to other volcanic areas throughout the world. This is interpreted as due to the presence of magma that affects anelasticity and scattering.     

单台Lg尾波Q_(Lg)~c的测量研究

将单台记录叠加频谱比法（简称SSR法）应用于震中距小于500km的区域地震Lg尾波为lHz下的值,n是与频率f有关的指数)的提取之中·为了检验方法的稳定性和可靠性,我们选用两组发生在不同地区的地震记录．一组是1989-10-18大同主震（ML＝5．9）的15个余震,北京遥测地震台网记录到的信噪比较高的Lg尾波序列。另一组是秦岭造山带数字流动地震观测台网记录到的1995年甘肃平凉两次地震信噪比较高的Lg尾波序列．每一组地震的震源参数基本相同,即它们到某一特定台站的路径几乎相同,得到的　可以认为是不受震源影响的地震一台站路径上的平均衰减效应·两区域　的测量结果表明：（1）尽管两个地震观测台网所用观测仪器不同,每个台站的观测条件不同,所记录的地震也不尽相同,但对Q₀和,的拟合结果线性程度都比较高,说明SSR法是稳定的;（2）同一震群的不同地震到某一特定台站测得的结果基本一致,即Q₀值及其平均值大都在其标准差范围之内,说明SSR法是可靠的。（3）对同一地震不同方位各台站的Lg尾波测量结果可以看出,穿过不同构造区的Q₀值存在明显差...     

Estimation of coda and shear wave attenuation in the volcanic area in SE Sabalan Mountain,NW Iran

The quality factors of coda and shear waves have been estimated for the SE Sabalan Mountain, geothermal region in northwestern Iran. We have analyzed 65 local earthquakes with magnitude of 2.8 to 6.1 and 2.8 to 5 for shear and coda wave quality factor estimation, respectively. These events were recorded on five stations installed by Building and Housing Research Center Network. Coda normalization and Spectral decay methods have been used to estimate the frequency dependence attenuation relation for shear wave, and single back-scattering method for coda waves. We have observed that the coda normalization method has supplied significantly higher Q _S values as compared to the spectral method. The results show that, in general, Q values are significantly smaller for the entire frequency range as compared to tectonically active areas and are close to the values for volcanic areas.     

Characteristics of radiation and propagation of seismic waves in the Baikal Rift Zone estimated by simulations of acceleration time histories of the recorded earthquakes

The regularities in the radiation and propagation of seismic waves within the Baikal Rift Zone in Buryatia are studied to estimate the ground motion parameters from the probable future strong earthquakes. The regional parameters of seismic radiation and propagation are estimated by the stochastic simulation (which provides the closest agreement between the calculations and observations) of the acceleration time histories of the earthquakes recorded by the Ulan-Ude seismic station. The acceleration time histories of the strongest earthquakes (M _W ~ 3.4–4.8) that occurred in 2006–2011 at the epicentral distances of ~96–125 km and had source depths of ~8–12 km have been modeled. The calculations are conducted with estimates of the Q-factor which were previously obtained for the region. The frequency-dependent attenuation and geometrical spreading are estimated from the data on the deep structure of the crust and upper mantle (velocity sections) in the Ulan-Ude region, and the parameters determining the wave forms and duration of acceleration time histories are found by fitting. These parameters fairly well describe all the considered earthquakes. The Ulan-Ude station can be considered as the reference bedrock station with minimum local effects. The obtained estimates for the parameters of seismic radiation and propagation can be used for forecasting the ground motion from the future strong earthquakes and for constructing the seismic zoning maps for Buryatia.     

Coda Q c Attenuation and Source Parameter Analysis in Friuli (NE Italy) and its Vicinity

—?The digital data acquired by 16 short-period seismic stations of the Friuli-Venezia-Giulia seismic network for 56 earthquakes of magnitude 2.3–4.7 which occurred in and near NE Italy have been used to estimate the coda attenuation Q _c and seismic source parameters. The entire area under study has been divided into five smaller regions, following a criterion of homogeneity in the geological characteristics and the constrains imposed by the distribution of available events. Standard IASPEI routines for coda Q _c determination have been used for the analysis of attenuation in the different regions showing a marked anomaly in the values measured across the NE border between Friuli and Austria for Q ₀ value. A large variation exists in the coda attenuation Q _c for different regions, indicating the presence of great heterogeneities in the crust and upper mantle of the region. The mean value of Q _c (f) increases from 154–203 at 1.5?Hz to 1947–2907 at 48?Hz frequency band with large standard deviation estimates.¶Using the same earthquake data, the seismic-moment, M ₀, source radius, r and stress-drop, Δσ for 54 earthquakes have been estimated from P- and S-wave spectra using the Brune's seismic source model. The earthquakes with higher stress-drop (greater than 1?Kbar) occur at depths ranging from 8 to 14?km.     

The study of the radiation characteristics and propagation of seismic waves in the North Caucasus by modeling the accelerograms of the recorded earthquakes

For evaluating the parameters of the vibrations of the Earth’s surface in the case of strong earthquakes, which are possible in the future, the regular patterns of the emission and propagation of seismic waves in the North Caucasus regions are investigated. The regional parameters of emission and propagation of seismic waves are evaluated by solution of the inverse problems of stochastic modeling of the accelerograms of the earthquakes, recorded by the seismic station in Sochi. The horizontal components of the strongest earthquakes (M _w ～ 3.9?5.6), that occurred in 2002–2006 within a radius of ～300 km from the seismic station, with source depths up to 60 km are modeled. For calculations of accelerograms, estimates of the quality are used, obtained earlier for this region in the form: Q(f) ～ 80 ～ f ^0.9. The parameter settings are carried out, which determine the shapes of the source spectra, the amplification of the seismic waves in the Earth’s crust, the weakening of the waves at high frequencies (κ), the parameters that determine the shape and duration of accelerograms, etc. Sufficiently good agreement of the calculated and recorded accelerograms is obtained, the regional characteristics of emission and propagation of seismic waves, which can be used for prediction of the parameters of strong motions in the North Caucasus, are evaluated; however, in the future these characteristics should be studied in more detail.     

Source Scaling of Inland Crustal Earthquake Sequences in Japan Using the S-Wave Coda Spectral Ratio Method

We estimate corner frequencies and stress drops for 298 events ranging from M _w 3.2–7.0 in 17 inland crustal earthquake sequences in Japan to investigate the source scaling and variation in stress drops. We obtain the source spectral ratio from observed records by the S-wave coda spectral ratio method. The advantage of using the S-wave coda is in obtaining much more stable source spectral ratios than using direct S-waves. We carefully examine the common shape of the decay of coda envelopes between event pair records. The corner frequency and stress drop are estimated by modeling the observed source spectral ratio with the omega-square source spectral model. We investigate the dependences of stress drops on some tectonic effects such as regionality, focal mechanism, and source depth. The principal findings are as follows: (1) a break in self-similar source scaling is found in our dataset. Events larger than M _w 4.5 show larger stress drops than those of smaller events. (2) Stress drops of aftershocks are mostly smaller than those of mainshocks in each sequence. (3) There are no systematic differences between stress drops of events occurring inside and outside the Niigata-Kobe Tectonic Zone in Japan. (4) Clear dependence of the faulting type on stress drops cannot be seen. (5) Stress drops of aftershocks depend on their source depth. (6) The crack size obtained from the corner frequency corresponds to the total rupture area of heterogeneous slip models for large events.     

天山中段2次6级地震前钻孔应变高频异常分析

2014年以来,天山中段分量钻孔应变仪空前增多,这些高采样率的应变观测资料蕴含着丰富的构造信息。如何从高采样率观测资料中提取有效的前兆异常信息,是分析研究人员亟待解决的问题。文章通过对天山中段分量钻孔应变观测数据进行S变换和超限率分析发现,在天山中段2次6级地震前有5套应变资料出现高频信息异常。这些异常均在震前出现,随后达到峰值,临震前或地震后衰减,其中短周期异常信号主要集中在10～720 min频段,且S变换与超限率分析结果具有很好的同步性。结合精河地震震源区及附近的GPS分析结果,发现高频异常信息的分布与该地区地壳运动场具有很好的一致性,进一步验证了高频信息异常的可信度。     

Scattering and attenuation of high-frequency body waves (1–25 Hz) in the lithosphere

A strong frequency dependence of apparent Q_β^?1 of shear waves was found for the frequency range from 1 to 25 Hz using band-pass filtered records of about 900 earthquakes occurring in the central Japan area with focal depths from 0 to 150 km. The method used for estimating Q_β is a single-station method based on elimination of the source effect from S waves by the use of coda spectra. The validity of the method is confirmed by agreement between the results obtained independently using data from two stations in the same area.     

基于CWI方法的门源MS6.9地震余震定位研究

门源发生M_S6.9地震后,余震目录中含有大量的单台事件,为了解决单台事件的定位不准确性,利用尾波干涉(CWI)的震源定位方法将门源台记录到震后8小时内的101次单台事件进行重定位,该方法主要通过计算事件对尾波部分的间距不确定性矩阵,从而估计聚类中具有相似震源机制的地震事件对之间的间距,通过优化这些事件对的间距来解决聚类中一组事件的相对位置。最后,对定位结果进行评估,得到相对可靠的定位结果。     

Using Airgun Source Signals to Study Regional Wave Velocity Changes before and after the Yunlong MS5.0 and Yangbi MS5.1 Earthquakes

On the basis of the airgun source signals recorded by the stations from January, 2016 to June, 2017, we use cross-correlation detection technology to obtain the characteristics of the stable phase travel time change of each station. We used the Yunlong M_S5.0 and Yangbi M_S5.1 earthquakes as samples. According to regional characteristics, 13 stations with high signal-to-noise ratios and complete data were selected (including 3 fixed stations and 10 active source stations). They are divided into four regions, and on the basis of the GNSS baseline data, the characteristics of regional wave velocity changes before and after the earthquake are analyzed. The results show that the station phase travel time change and the regional stress characteristics represented by the GNSS baseline data have good correlation in the short-term. Due to different degrees of regional stress, there are differences in the travel time changes of different stations in the four regions. Before the Yunlong M_S5.0 and Yangbi M_S5.1 earthquakes, with regional stress adjustment, there is an upward trend in the travel time changes of related stations in the adjacent areas of up to 0.02s. The difference is that there are differences in the time nodes and duration of the travel time anomalies, and there is a reverse descent process after the Yangbi M_S5.1 earthquake. There are different degrees of travel time fluctuations in the relevant stations before and after the two earthquakes, but the fluctuation range before and after the earthquake was small. Compared with the water level change of the reservoir, the adjustment of the regional stress is more likely to have a substantial impact on the travel time changes of the relevant stations.