Spectral determination of source parameters in the Marmara Region

Ever since the 1999 Kocaeli earthquake, in which the Kandilli Observatory and Earthquake Research Institute (KOERI) was not able to correctly reflect the magnitude size in its preliminary report because of the saturation effect, a rapid and accurate determination of the earthquake becomes a very important issue. Therefore, in the framework of this study, an automatic determination of the moment magnitude was performed by using the displacement spectra of selected earthquakes in the Marmara Region. For this purpose, 39 three-component broadband stations from KOERI seismic network which recorded 174 earthquakes with magnitudes 2.5?≤?M?≤?5.0 in between 2006–2009 were used. Due to the importance of quality factor in determination of the moment magnitude with spectral analysis method, the quality factor was calculated for the whole region in the beginning. Source spectrum which was obtained by converting the velocity records to displacement spectra and moment magnitudes of earthquakes were determined by fitting this spectrum to classical Brune model. For this aim, an automatic procedure was utilized which based on minimizing the differences between observed and synthetic source spectra identified by the S waves. Besides moment magnitude and location parameters, some source parameters such as seismic moment, spectral level, corner frequency and stress drop were also calculated. Application of the method proves that determining the seismic moment from the source spectra is applicable not only for earthquakes with small magnitude but also moderate earthquakes as well.     

内蒙古中西部地区震源参数研究

选取内蒙古测震台网2009—2016年内蒙古中西部M_L≥2.8中小地震波形记录,采用多台联合反演方法,计算该区中小地震震源谱参数,获得地震矩M₀、矩震级M_W、应力降Δσ和震源尺度r,利用线性回归,分析近震震级M_L、地震矩M₀、矩震级M_W、应力降Δσ和震源尺度r之间的关系。结果表明,各参数之间存在一定线性关系,内蒙古中西部地区应力降模型为增加应力降模型（ISD）。     

Source parameters of some significant earthquakes near Koyna dam,India

Seismic moment, stress drop, source radius, and fault dislocation have been determined for nineteen significant earthquakes in the Konya area, India using displacement spectra of shear waves computed from strong-motion accelerograph records. Though the stress-drop shows a definite increasing trend with the seismic moment, its correlations with source radius and corner frequency indicate that a constant stress drop of 170 bars represents a good mean value to describe the source mechanism of the Koyna dam earthquakes. An empirical relationship also has been established between magnitude and seismic moment.     

浙江珊溪水库地震震源参数特征研究

利用2014年1月-2015年12月浙江省区域数字地震台网记录的1 360次珊溪水库M_L0.5-4.4地震资料,采用Brune模式,将台站获取的速度记录进行几何扩散校正、介质衰减校正及仪器校正后得到速度谱,对速度谱在频率域进行傅里叶积分得到震源位移谱,之后使用遗传算法计算得到拐角频率和震源位移谱零频极限值,进而计算出该区的地震矩为1010-1014 N·m,震源破裂半径为33-550 m,拐角频率为2.4-39.7 Hz,地震应力降为0.04-6.74 MPa,视应力为0-2.75 MPa.在此基础上,对各参数特征及其关系进行系统性分析,结果表明:各震源参数之间表现出一定的对数线性或半对数线性关系;空间上,应力降和拐角频率在库区西北段较高,东南段较低;应力降与地震矩在空间上呈反向关系,拐角频率与破裂半径在空间上呈负相关.      

内蒙古中西部地区中小地震矩震级研究

基于S震相"S窗"内的波形信号识别、品质因子Q(f)和22个台站场地响应,利用2009~2016年3月内蒙古中西部地区地震的波形资料,反演了182次中小地震的震源波谱参数,得到这些小震的零频幅值及其拐角频率,据此计算了这些地震的地震矩M_0、矩震级M_W和应力降Δσ。利用回归分析方法得到了近震震级与矩震级、矩震级与应力降的关系式。分析表明,近震震级与矩震级、矩震级与应力降呈线性关系。可见,将矩震级纳入地震的快报与正式目录中,可以丰富地震观测报告内容,更好地为地震应急和地震科研服务。     

Study on the Moment Magnitude of Small and Moderate Earthquakes Located in the Inner Mongolia Region

Based on digital seismic waveform data from Inner Mongolia Digital Seismic Network, the source spectrum parameters of 182 small and moderate earthquakes from January, 2009 to September, 2016 are derived, and the seismic moment M₀ and moment magnitude M_W of the earthquakes are calculated. The M_L-M_W relationship and the relationship between stress drop and magnitude are obtained using the linear regression method. It is clear that incorporating the moment magnitude into the seismic quick report catalog and the official earthquake catalog can enrich earthquake observation report content, thus providing better service for earthquake emergency and earthquake scientific research.     

云南地区中小地震静力学和动力学参数定标关系

Source spectra for moderate and small earthquakes are obtained after removing the path effect, site effect, and instrument response, etc. in the observed S-wave spectra. Based on the Brune source model...     

Use of long-period surface waves for rapid determination of earthquake-source parameters

A method for rapid retrieval of earthquake-source parameters from long-period surface waves is developed. With this method, the fault geometry and seismic moment can be determined immediately after the surface wave records have been retrieved. Hence, it may be utilized for warning of tsunamis in real time. The surface wave spectra are inverted to produce either a seismic moment tensor (linear) or a fault model (nonlinear). The method has been tested by using the IDA (International Deployment of Accelerographs) records. With these records the method works well for the events larger than M_s = 6, and is useful for investigating the nature of slow earthquakes.For events deeper than 30 km, all of the five moment tensor elements can be determined. For very shallow events (d ? 30 km) the inversion becomes ill-conditioned and two of the five source moment tensor elements become unresolvable. This difficulty is circumvented by a two-step inversion. In the first step, the unresolvable elements are constrained to be zero to yield a first approximation. In the second step, additional geological and geophysical data are incorporated to improve the first approximation. The effect of the source finiteness is also included.     

Toroidal free oscillations of the Earth observed by a ring laser system: a comparative study

In this study, we explore the potential of measuring systematically the Earth's free oscillations using ring laser gyro (RLG) vertical axis rotational records. The RLG that we use is the vertical axis G-ring laser system of the Geodetic Observatory Wettzell (Germany). In 2009, its signal-to-noise ratio was considerably improved over the broadband frequency range of seismic measurements. Since then, three large magnitude earthquakes have occurred (Samoa Islands 2009; Maule, Chile, 2010; and Tohoku, Japan, 2011), leading to the first direct observations of rotational ground motions induced by toroidal free oscillations of the Earth. We compare these G-ring laser observations with synthetic seismograms computed by summing normal modes. We also analyse amplitude spectra of real and synthetic data to aid in the interpretation of the observations. We show that several toroidal modes are detected by the G-ring laser for earthquakes with a moment magnitude M _W????8.0 and that our observations are in reasonable agreement with the synthetic spectra. We also report evidence for mode coupling in both translation and rotation spectra.     

On the dynamics and origin of the Pugachevo group of gas-water-lithoclast (“mud”) volcanoes on Sakhalin Island: Visual observations and orohydrography

Visual observations of the activities of three gas-water-lithoclast (“mud”) volcanoes in the Pugachevo group in the Makarov District of Sakhalin Island were used to analyze the dynamics of these volcanoes. Our analysis confirmed the previous inference that a multirank periodicity is present in the eruptions of the volcanoes, with the periodicity being comparatively frequent and low in magnitude for the eruptions (every 1–2 years) at the Glavnyi or Central volcano and a very long but violent (every 65–70 years) at all the three volcanoes. Comparison of this periodicity with a similar periodicity observed in natural seismicity in the adjacent, highly seismic Uglegorsk District shows that the two phenomena are relatively independent. The natural seismicity in the form of large earthquakes is superimposed upon the gas-water-lithoclast volcanism mostly to serve as a trigger. The high-rank periodicity of 65–70 years noted above may be related to 11–22-year cycles of solar activity. Analysis of the orohydrographic setting in the area shows that the location of the present low-lying basin with volcanoes inside was formerly occupied by a dome uplift with a major gas field that controlled the generation of these volcanoes; it continues to drive their activity, causing a continued gradual subsidence or collapse of the deep basin with volcanoes, thus providing evidence of the high gas potential in the area of study.     

应用近场速度记录快速估算川滇地区中强地震强度

采用Wu等(2004)提出的地震强度测定方法,应用近场速度波形记录,测定地震参数M_(ew),实现了川滇地区中强地震(M_W≥5.7)强度近实时快速测定。M_(ew)是一个与矩震级M_W对接的地震参数,与传统的应用振幅测定震级的方法相比,不会存在震级饱和问题,与全球矩心矩张量(GCMT)的矩震级M_W具有较好的一致性。该方法解决了应用近场波形资料和矩张量反演方法求取大地震矩震级比较困难的问题,可为区域中强地震速报提供近实时矩震级提供参考,从而提高地震速报的时效性和准确性,同时为强震加速度数据的应用奠定基础。     

Moment Magnitude and Its Calculation

In this paper, we give a brief introduction to the proposal and development history of the earthquake magnitude concept. Moment magnitude M_W is the best physical quantity for measuring earthquakes. Compared with other magnitude scales used traditionally, moment magnitude is not saturated for all earthquakes, regardless of big and small earthquakes, deep and shallow earthquakes, far field and near field seismic data, geodetic and geological data, moment magnitude can be measured, and can be connected with well-known magnitude scales such as surface wave magnitude M_S. Moment magnitude is a uniform magnitude scale, which is suitable for statistics with wide magnitude range. Moment magnitude is the preferred magnitude selected by the International Seismological community, and it is preferred by the departments responsible for publishing seismic information to the public.Moment magnitude is a uniform magnitude scale, which is suitable for statistics with wide magnitude range. Moment magnitude is a preferred magnitude for international seismology, it is preferred by the agency responsible for providing information about earthquakes to the public. We provide all formulas used in the calculation of moment magnitude, and the calculation steps in detail. We also analyzed some problems and rules to solve these problems by using different formulas and numerical value calculation steps.     

Volcanic activity and large earthquakes

Several investigators have noted an increase in volcanic activity prior to several large interplate earthquakes in the circum-Pacific area. From these data, I have examined the temporal relationship between the beginning of this increased volcanic activity and the occurrence of large earthquakes. I find that the time (year) of the beginning of enhanced volcanic activity depends on the magnitude of the subsequent large earthquake, for earthquakes with magnitude > 7.5. The volcanoes which exhibit this enhanced activity are situated on the landward side of the rupture zone of the large earthquake, and this study therefore suggests a casual physical relationship between the two phenomena. I also find that this increase in volcanic activity occurs a few years before an increase in the occurrence of smaller earthquakes, noted by several invesigators, in the rupture zones of large earthquakes.     

河北唐山地区中小地震震源参数的测定

采用谱分析的方法，利用首都圈数字遥测地震台网的观测资料，计算了唐山地区2002年ML2．0以上的58个中小地震的波谱，进而估算了地震矩、拐角频率、应力降、震源破裂半径，并讨论了它们之间的相互关系以及拐角频率和应力降的时间分布特征。分析表明，地震矩与震级、拐角频率与震级的关系密切，证实结果可信、稳定。     

Site response and source spectra of S waves in the Zagros region, Iran

S wave amplitude spectra from shallow earthquakes with magnitudes ranging between 4.2 and 6.2 in the Zagros region of Iran that occurred between 1998 and 2008 are used to examine source parameters and site response of S waves. A generalized inversion scheme has been used to separate the source, propagation path, and local site effects from S wave spectra. For removing the trade-off between source and site terms and propagation effects (including geometric and anelastic attenuation), the spectral amplitudes of the records used were corrected for attenuation and geometrical spreading function using a path model proposed by Zafarani and Soghrat (Bull Seism Soc Am 102:2031–2045, 2012) for the region. We assume a Brune’s point source model to retrieve source parameters like corner frequency, moment magnitude, and high-frequency fall off coefficient, for each event. When the source spectra are interpreted in terms of Brune’s model, the average stress drops obtained are about 7.1 and 5.9 MPa (71 and 59 bars), respectively for the eastern and western Zagros regions. Stress drops range from 1.4 to 35.0 MPa (14 to 350 bars), with no clear dependence on magnitude. The results in terms of stress drop and S wave seismic energy indicate that the Zagros events are more similar to interplate earthquakes of western North America than to intraplate events of eastern North America. The method also provides us with site responses for all 40 stations individually and is an interesting alternative to other methods, such as the H/V method. A new empirical relationship between body-wave magnitudes and moment magnitude has been proposed for the Iranian plateau using derived seismic moment from the inversion.     

峰值速度和加速度对环境剪应力的依赖性

从导出的地震定标律和地震破裂过程的断裂力学模式出发,得到了震源平均位移、速度和加速度谱的表达式,进而又推导出震源的峰值位移d_m,速度v_m和加速度a_m的表达式:d_m=k_dM_0~2/~3τ_0~(2/3),v=k_vM_0~(1/3)τ_0~(4/3),a_m=k_aτ_0~2式中M_0是地震矩,τ_0是环境剪应力值,k_d,k_v,k_a为适当的常数.我们选用了66个地震的观测资料,这些地震的矩震级范围包括了从1级左右的极微震,3-5级的小震,直到6-7级的大震;地震矩从10~9-10~(20)Nm,跨越了10个数量级,并用这些地震检验了上述公式. 令所有地震的平均应力为5MPa,定出常数k_d,k_v,k_a,进而由速度和加速度观测资料求得66个地震的环境剪应力τ_0值,这些数值相当稳定.多数极微震的τ_0值在2-4MPa之间;小震的τ_0值多数在4-8MPa左右;大震的τ_0值为10MPa左右。τ_0值对震源深度和断层类型有明显的依赖性.一般深度很浅的小震和极微震,τ_0值很低;正断层地震的τ_0值相对较低;逆断层地震的τ_0值较高;走滑断层地震的τ_0值则居中.     

An investigation of seismicity for the Central Anatolia region,Turkey

The aim of this study is to investigate the seismicity of Central Anatolia, within the area restricted to coordinates 30–35° longitude and 38–41° latitude, by determining the “a” and “b” parameters in a Gutenberg–Richter magnitude–frequency relationship using data from earthquakes of moment magnitude (Mw)?≥?4.0 that occurred between 1900 and 2010. Based on these parameters and a Poisson model, we aim to predict the probability of other earthquakes of different magnitudes and return periods (recurrence intervals). To achieve this, the study area is divided into six seismogenic zones, using spatial distributions of earthquakes greater than Mw?≥?4.0 with active faults. For each seismogenic zone, the a and b parameters in the Gutenberg–Richter magnitude–frequency relationship were calculated by the least squares method. The probability of occurrence and return periods of various magnitude earthquakes were calculated from these statistics using the Poisson method.     

基于广义S变换的地震与爆炸识别

为了进一步增强区域台网针对天然地震与人工爆炸事件的识别能力,本文利用广义S变换方法,围绕河北省三河采石场爆破以及周边发生的天然地震波形记录展开研究,总结了地震与爆炸在时频谱中的差异性。研究表明,受震源机制影响,天然地震时频谱图的频带范围更宽广,能量团分布也较人工爆炸更为复杂。为了消除震中距的影响,本文将求解的时频谱保存为规格相等的灰度图像,通过滑动窗口计算图像的灰度一致性得到新的识别判据?谱图二阶矩。对三河地区已知地震和爆炸事件的测试表明,谱图二阶矩对单台波形记录的识别率可达91%,多台求取平均值之后的识别率超过98%,因此该判据具有良好的应用前景。      

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.     

甘肃地区震源参数的相关性研究

通过对甘肃数字地震台网2008年1月-2009年6月记录的ML 2.0-ML 6.0中小地震,S波观测谱中噪音、仪器、路径、场地等影响的逐一消除,利用遗传算法反演得到甘肃地区273条地震的震源谱,进而根据Brune的圆盘震源模型,计算地震矩、应力降、拐角频率、震源半径等震源参数,并研究甘肃地区震源参数的特征及相似关系....