1984年5月21日南黄海地震的震源参数——用IDA记录测定浅震震源参数

本文讨论了利用国际部署加速度仪台网(IDA)的超长周期瑞雷波记录,反演震源参数的有关理论和方法.对于震源深度大于30km的地震,地震矩张量的五个元素都可以测定;对于浅源(h30km)地震,其中两个元素变为无解[1,2,3].本文以1984年5月21日南黄海浅源地震为例,研究了在浅震难题的情况下,如何更多地提取震源信息的问题,测定了主震的断层机制为:走向北35西,倾角52,滑动方向54,地震矩为1.21018Nm,压应力主轴方向为北74.4东,是一个走向北北西的左旋逆断层.本文还阐述了分布源模型问题。结果表明,采用分布源与采用点源模型对南黄海这种中强地震反演结果的影响不大.分析了震源有限性问题,提出了进行修正的一种明确表达式.      

Moment Tensor Inversion of Regional Phases: Application to a Mine Collapse

—?The aim of our study consists of analyzing potentially non-double-couple seismic events recorded at regional distances. In order to define the nature of the seismic source, a moment tensor inversion is carried out as this method is general enough not to initially constrain the source mechanism. In this paper we present an application to a seismic event induced by a mine collapse which occurred near the town of Halle in Germany. Because of its induced nature, many parameters such as the location and geometry of this seismic source are known. This information allows us to test the influence of inadequate propagation modeling on the moment tensor obtained from the inversion. Green's functions have been computed with the reflectivity method in a flat layered medium, using the European model EurID (Du et? al., 1998; Dufumier et al., 1997). From the inversion of P-wave seismograms recorded by the German Regional Seismic Network will, we obtained a source time function which can be decomposed into two subevents. The first one has a large isotropic part and a deviatoric mechanism with near vertical nodal planes. No volume change is observed for the second subevent, but a deviatoric component opposite of the first one. The addition of S-waves does not change the results of the inversion which are stable. Surface waves were not used because of their poor dispersion curves. Based on the moment tensor obtained from these inversions, the physical process at the source is compatible with a large cavity collapse.     

The 2011 Lorca seismic series: Temporal evolution,faulting parameters and hypocentral relocation

The Lorca 2011 seismic series was recorded by an unprecedented set of high quality on scale broadband seismograms and strong motion accelerograms. The waveforms from permanent and temporary broadband seismic networks deployed in the region by different institutions allowed to invert regional moment tensor for the fore, main and largest aftershock of the complete seismic sequence. Using double-difference algorithm we have performed a precise relocation of the seismic series, where body wave travel times from strong ground motion accelerograms were included in the data set. Regional moment tensor inversion for the three main events show similar oblique-reverse faulting regime with a northeast-southwest fault orientation. The scalar seismic moment, moment magnitude and focal depth retrieved from the inversion yield the following values for each event: \(\hbox {Mo}=6.5\times 10^{16}\) Nm (Mw = 5.2) for the mainshock, \(\hbox {Mo}= 9.6 \times 10^{15}\) Nm (Mw = 4.6) for the foreshock and \(\hbox {Mo}=7.3\times 10^{14}\) Nm (Mw = 3.9) for the large aftershock. The centroid depths range between 4 and 6 km. The double-difference relocation of the seismic series shows significant epicentral differences with the preliminary routine location. The epicentral solutions given by this relocation show a seismic sequence distributed following a NE–SW strike, subparallel to the Alhama de Murcia fault and compatible with the faulting parameters inverted from the moment tensor analysis. The hypocenters of the series generate a subvertical trend in depth distribution, being concentrated between 2 and 6 km. The depth distribution of the main events, which range from 4.6 to 5.5 km, is in good relationship with the faulting and depth parameters deduced from the moment tensor inversion technique. The regional moment tensor solutions for the three largest earthquakes, the epicentral distribution and the focal depths show good relationship with the surface geometry and tectonic regime of the Alhama de Murcia fault. The stress drop deduced for the mainshock gives a value ranging between 58 and 85 bars, which does not support the idea of a high stress drop release as a main factor contributing to the high ground acceleration recorded at Lorca. The PGA values observed at Lorca, which contributed to the high damage independently of structural deficiencies, could be generated mainly by shallowness and proximity to the seismic source together with a directivity effect in the seismic radiation.     

The determination of the best fitting and the extreme source models by a global inversion method

In this paper, a global inversion method is developed for seismic moment tensor inversion by using the body wave forms. The algorithm depends on neither the selection of starting model nor the forms of objective function and constraints. When the error function, measure of the difference between the observed and synthetic waveforms, is chosen as the objective function, the best fitting source model is found; when a certain combination in seismic moment tensor elements is selected as the objective function and the values of error function are constrained in a suitable bound, the extreme source models can be produced by minimizing or maximizing this combination. By changing the form of the combination of moment tensor elements, a variety of different source characteristics can be considered. Therefore the extreme solution provides an estimation of the uncertainty in the best fitting source model. The seismic waveform data was used to evaluate the effectiveness of this algorithm. This research was supported by the National Natural Science Foundation of China.     

Sensitivity of the strong ground motion time histories to a finite source model: A case study for the January 12, 2010 Haiti earthquake (Mw=7.0)

We analyze the waveforms generated by the January 12, 2010 Haiti earthquake (M_w=7.0) for its source characteristics. A 60 to 25 km source model is retrieved by the Kikuchi and Kanamori finite source inversion technique that uses broadband teleseismic body wave records. The derived rupture model points out unilateral rupture propagation commenced at the eastern side of the fault plane where the major seismic moment release occurred. The rupture front propagated westward and terminated at a site where the largest aftershocks occurred. Our estimates yield a seismic moment of M_o=8.17×10¹⁹ N m released on a 60 km-long fault plane. A patch at the eastern side of the ruptured fault plane inferred as a region of maximum moment release.     

2013年前郭MS5.8震群矩张量研究

2013年吉林前郭M_S5.8震群为爆发性震群,目前余震活动仍然在持续.基于吉林、辽宁、黑龙江、内蒙古四省地震台网记录的前郭震群波形资料,利用波形信噪比、震源类型、台站及速度模型组合的指标选择最佳的反演方案,应用矩张量的三种反演模式,对序列中5个M_S≥5.0地震进行矩张量反演研究,获得了全矩张量、偏矩张量和纯双力偶的矩张量.使用F-test对地震的三种模式的矩张量反演结果进行显著性检验来确定最佳反演模式.结果显示,5个地震的最优矩张量解均为全矩张量模式反演获得的结果,其双力偶分量仅有20%~65%,矩心深度位于地下3~4 km处,地震在Hudson震源类型图上的投影远离双力偶震源类型区域.这些结果表明,震源类型并非典型的构造地震,推断前郭地震可能是与人类活动有关的诱发地震.     

非对称地震矩张量时间域反演:理论与方法

本文在对称地震矩张量反演的基础上,进一步研究了非对称地震矩张量时间域反演的理论与方法,结果表明:非对称地震矩张量反演与对称地震矩张量反演类似,只需将对称地震矩张量反演方法略加改动,即增加3个待解参数,便可实现非对称地震矩张量反演.为了判断非对称地震矩张量反演相对于对称地震矩张量反演是否存在过度拟合,运用了AIC准则 (赤池信息准则).为了定量地描述地震矩张量之间的差异,引入了地震矩张量的矢量表示法.通过分析格林函数与地震矩张量各分量之间的关系,得出:在非对称地震矩张量反演时,若仅用垂直向地动位移数据,将无法区分M_xy与M_yx这两个分量, 需要同时运用垂直向与水平向地动位移数据进行联合反演才能区分M_xy与M_yx; 若采用不同的速度结构模型或不同的格林函数计算方法,则需重新评估地震矩张量各分量的分辨度问题.为检验非对称地震矩张量反演方法的可行性, 利用合成地震图进行了一系列数值试验.数值试验结果表明,在非对称地震矩张量反演中,有必要引入S波进行P波与S波联合反演以提高反演的准确性和判定断层面的能力.      

Retrieval of source parameters of an event of the 2000 West Bohemia earthquake swarm assuming an anisotropic crust

We propose an inversion scheme for retrieval of characteristics of seismic point sources, which in contrast to common practice, takes into account anisotropy. If anisotropy is neglected during inversion, the moment tensors retrieved from seismic waves generated by sources situated in anisotropic media may be biased. Instead of the moment tensor, the geometry of the source is retrieved directly in our inversion; if necessary, the moment tensor can be then determined from the source geometry aposteriori. The source geometry is defined by the orientation of the slip vector and the fault normal as well as the strength of the event given by the size of the slip and the area of the fault. This approach allows direct interpretation of the source geometry in terms of shear and tensile faulting. It also makes possible to identify volumetric source changes that occur during rupturing. We apply the described algorithm to one event of the 2000 West Bohemia earthquake swarm episode. For inversion we use information of the direct P waves. The structure is approximated by three different models determined from travel-time observations. The models are inhomogeneous isotropic, inhomogeneous anisotropic, and homogeneous anisotropic. For these models we obtain seismic moments M_T = 3.2 − 3.8 × 10¹⁴ Nm and left-lateral near-vertical oblique normal faulting on a N-S trending rupture surface. The orientation of the rupture surface is consistent with fault-plane solutions of earlier studies and with the spatial distribution of other events during this swarm. The studied event seems to be accompanied by a small amount of crack opening. The amount of crack opening is slightly reduced when the inhomogeneous anisotropic model is assumed, but it persists. These results and additional independent observations seem to indicate that tensile faulting occurs as a result of high fluid pressure.     

Source process and tectonic implications of the Spanish deep-focus earthquake of March 29, 1954

The source process of the deep-focus Spanish earthquake of March 29, 1954 (m_b = 7.1, h = 630 km) has been studied by using seismograms recorded at teleseismic distances. Because of its unusual location, this earthquake is considered to be one of the most important earthquakes that merit detailed studies. Long-period body-wave records reveal that the earthquake is a complicated multiple event whose wave form is quite different from that of usual deep earthquakes. The total duration of P phases at teleseismic distances is as long as 40 s. This long duration may explain the considerable property damage in Granada and Malaga, Spain, which is rather rare for deep earthquakes. Using the azimuthal distribution of the differences between the arrival times of the first, the second and later P phases, the hypocenters of the later events are determined with respect to the first event. The focus of the second event is located on the vertical nodal plane of the first shock suggesting that this vertical plane is the fault plane. This fault plane which strikes in N2°E and dips 89.1°E defines a nearly vertical dip-slip fault, the block to the west moving downwards. The time interval and spatial separation between the first and the second events are 4.3 s and 19 km respectively, giving an apparent rupture velocity of 4.3 km/s which is about 74% of the S-wave velocity at the source. A third event occurred about 8.8 s after the first event and about 35.6 km from it. At least six to ten events can be identified during the whole sequence. The mechanism of some of the later events, however, seems to differ from the first two events. Synthetic seismograms are generated by superposition of a number of point sources and are matched with the observed signals to determine the seismic moment. The seismic moments of the later events are comparable to, or even larger than, that of the first. The total seismic moment is determined to be 7 · 10²⁷ dyn cm while the moments of the first and the second shocks are 2.1 · 10²⁶ dyn cm and 5.1 · 10²⁶ dyn cm, respectively. The earthquake may represent a series of fractures in a detached piece of the lithosphere which sank rapidly into the deep mantle preserving the heterogeneity of material property at shallow depths.     

Specific features of surface wave radiation by a shallow source

The uncertainty of the seismic moment tensor inversion from records of long period surface waves radiated by a shallow source is considered. The special cases of sources of pure thrust, normal, and strike-slip faults are examined in detail. Results of numerical modeling of Love and Rayleigh wave radiation patterns and frequency dependences of the radiation intensity are presented for these types of sources. It is shown that, in the general case, the focal mechanism and seismic moment of a shallow double couple can be uniquely determined from long period surface waves only if one of the nodal planes is subhorizontal. The threshold (maximum) value of the dip angle of this plane is determined by the source depth, the spectral range of observed surface waves, and the model structure in the vicinity of the source.     

Seismic moment tensor inversion using 3D velocity model and its application to the 2013 Lushan earthquake sequence

Source inversion of small-magnitude events such as aftershocks or mine collapses requires use of relatively high frequency seismic waveforms which are strongly affected by small-scale heterogeneities in the crust. In this study, we developed a new inversion method called gCAP3D for determining general moment tensor of a seismic source using Green's functions of 3D models. It inherits the advantageous features of the “Cut-and-Paste” (CAP) method to break a full seismogram into the Pnl and surface-wave segments and to allow time shift between observed and predicted waveforms. It uses grid search for 5 source parameters (relative strengths of the isotropic and compensated-linear-vector-dipole components and the strike, dip, and rake of the double-couple component) that minimize the waveform misfit. The scalar moment is estimated using the ratio of L₂ norms of the data and synthetics. Focal depth can also be determined by repeating the inversion at different depths. We applied gCAP3D to the 2013 M_s 7.0 Lushan earthquake and its aftershocks using a 3D crustal-upper mantle velocity model derived from ambient noise tomography in the region. We first relocated the events using the double-difference method. We then used the finite-differences method and reciprocity principle to calculate Green's functions of the 3D model for 20 permanent broadband seismic stations within 200 km from the source region. We obtained moment tensors of the mainshock and 74 aftershocks ranging from M_w 5.2 to 3.4. The results show that the Lushan earthquake is a reverse faulting at a depth of 13–15 km on a plane dipping 40–47° to N46° W. Most of the aftershocks occurred off the main rupture plane and have similar focal mechanisms to the mainshock's, except in the proximity of the mainshock where the aftershocks' focal mechanisms display some variations.     

Relative seismic moment tensor determination for Vrancea intermediate depth earthquakes

The method of relative seismic moment tensor determination proposed byStrelitz (1980) is extended a) from an interactive time domain analysis to an automated frequency domain procedure, and b) from an analysis of subevents of complex deep-focus earthquakes to the study of individual source mechanism of small events recorded at few stations.The method was applied to the recovery of seismic moment tensor components of 95 intermediate depth earthquakes withM _L=2.6–4.9 from the Vrancea region, Romania. The main feature of the obtained fault plane solutions is the horizontality ofP axes and the nonhorizontal orienaation ofT axes (inverse faulting). Those events with high fracture energy per unit area of the fault can be grouped unambiguously into three depth intervals: 102–106 km, 124–135 km and 141–152 km. Moreover, their fault plane solutions are similar to ones of all strong and most moderate events from this region and the last two damaging earthquakes (November 10, 1940 withM _W=7.8 and March 4 1977 withM _W=7.5) occurred within the third and first depth interval, respectively. This suggests a possible correlation at these depths between fresh fracture of rocks and the occurrence of strong earthquakes.     

The April 24, 2013 Changning M s4.8 earthquake: a felt earthquake that occurred in Paleozoic sediment

The dense broadband seismic network provides more high-quality waveform that is helpful to improve constraint focal depth of shallow earthquake. Many shallow earthquakes occurring in sediment were regarded as induced events. In Sichuan basin, gas industry and salt mining are dependent on fluid injection technique that triggers microseismicity. We adopted waveform inversion method with regional records to obtain focal mechanism of an M _s4.8 earthquake at Changning. The result suggested that the Changning earthquake occurred at a ESE thrust fault, and its focal depth was about 3 km. The depth phases including teleseismic pP phase and regional sPL phase shows that the focal depth is about 2 km. The strong, short-period surface wave suggests that this event is a very shallow earthquake. The amplitude ratio between Rayleigh wave and direct S wave was also used to estimate the source depth of the mainshock. The focal depth (2–4 km) is far less than the depth of the sedimentary layer thickness (6–8 km) in epicentral region. It is close to the depth of fluid injection of salt mining, which may imply that this event was triggered by the industrial activity.     

利用区域宽频地震数据反演2021年5月云南漾濞MS6.4地震震源破裂过程

本文基于有限断层模型反演方法,利用区域宽频带数据反演了2021年5月云南漾濞M_S6.4地震的震源破裂过程,结果显示:此次地震的发震断层走向为SE向,主要以右旋走滑为主.破裂主要发生在震源东南侧,最大错动量约为0.55 m,位于深度约9 km处,发生明显破裂的深度约达13 km.此次地震释放的标量地震矩为1.48×10¹⁸N·m,相当于矩震级M_W6.05.地震能量主要在前11 s释放.在深度为6~8 km处破裂速度有明显的变快,可能加剧了地表的震动.     

利用区域宽频地震数据反演2021年5月云南漾濞MS6.4地震震源破裂过程

本文基于有限断层模型反演方法,利用区域宽频带数据反演了2021年5月云南漾濞M_S6.4地震的震源破裂过程,结果显示:此次地震的发震断层走向为SE向,主要以右旋走滑为主.破裂主要发生在震源东南侧,最大错动量约为0.55 m,位于深度约9 km处,发生明显破裂的深度约达13 km.此次地震释放的标量地震矩为1.48×10¹⁸N·m,相当于矩震级M_W6.05.地震能量主要在前11 s释放.在深度为6~8 km处破裂速度有明显的变快,可能加剧了地表的震动.     

利用大地测量资料反演地震震源参数的若干问题

本文利用大地测量资料反演地震震源参数的一般理论,归纳为震源模型的建立、最优化计算的实施和观测资料的使用等三部分.文中着重研究了后两部分:提出了利用原始观测资料(地面长度、角度、高差、倾斜、应变的变化值)直接进行反演的严密方法;改进了单纯形最优化计算方法,加速了迭代的收敛并给出了收敛准则;按逐渐趋近法进行观测资料的标准化,保证了标准化后的各类资料都属于同一正态分布.文中根据地震前、后的大地测量资料,对1976年唐山7.8级地震的震源参数进行了反演,并对成果的稳定性和可靠性进行了检验,得到地震断层长度为112公里,走向为北东56,倾角为90,断层破裂至地面,断层面延深(宽)15公里,水平错距(右旋)2.51米,垂直错距(东南盘下降)0.70米,断层迹线中点位置为北纬3936'6、东经11811'4.由此计算得到的地震矩为1.451027达因厘米,应力降为37.2巴,应变降为5.610-5.这个结果与根据地震波及震时其它地球物理资料求得的结果相近,因此认为本次地震前(或后),发生相当于主震规模的断层蠕动的可能性不大.      

用应变地震观测求解震源矩张量的基本原理

震源机制解,即对地震矩张量的推断,对于地震研究具有至关重要的意义.应变地震观测是张量观测,与摆式地震仪的位移矢量观测不同,可以为地震研究提供新的数据源.本文讨论用应变地震观测求解震源矩张量的基本原理.在距离震源足够远的地方,地震波可以看成平面波,其性质决定于震源矩张量.假设平面地震波的应变张量可以由震源矩张量通过坐标变换计算得到,就可以通过观测应变地震波求解震源机制.这个假设至少对于双力偶震源机制是成立的.由此可以证明,在理想的无限介质中,只要有两个以上不同地点的应变地震波观测,就可以解出震源矩张量.这为解决震源机制问题提供了新的方法.目前的地震矩张量求解方法需要两方面的条件:或者需要很多观测点(例如体波反演),或者需要长周期地震波资料(例如面波反演).这些方法只适用于分析比较大的地震.对于小震,因为通常其震中周围不会有足够多的摆式地震仪观测点观测到其地震波,而地震波周期又短,难以利用传统方法给出可靠的震源机制解,所以只需少数观测点就能求解震源矩张量的新方法特别有意义.用应变地震观测求解震源机制,可以给出更为精确的结果.     

FOCAL MECHANISM OF THE AUGUST 8, M7.0 SICHUAN JIUZHAIGOU AND AUGUST 9, M6.6 XINJIANG JINGHE EARTHQUAKES OF 2017

The W-phase is a long period phase arriving between the P and S wave phases of a seismic source, theoretically representing the total near-and far-field long-period wave-field. Recent study suggests that the reliable source properties of earthquake with magnitude greater than ~M_W4.5 can be rapidly inverted by using the W-phase waveform data. With the advantage of W-phase, most of major earthquake research institutes in the world have adopted the W-phase based inversion method to routinely assess focal mechanism of earthquake, such as the USGS and GFZ. In this study, the focal mechanism of the August 8, 2017 M7.0 Sichuan Jiuzhaigou and August 9, 2017 M6.6 Xinjiang Jinghe earthquakes were investigated by W-phase moment tensor inversion technique using global seismic event waveform recordings provided by Incorporated Research Institutions for Seismology, Data Management Center. To get reliable focal mechanism, we strictly select raw waveform data and carry out inversion in stages. At first, we discard waveform without correct instrument information. Then we carry out an initial inversion using selected waveform data to get primary results. Using the preliminary results as input, we carry out grid-search based inversion to find the final optimal source parameters. The inverted results indicate that the August 8, M7.0 Sichuan Jiuzhaigou shock resulted from rupturing on a NW-trending normal fault with majority of strike-slip movement. The parameters of two nodal planes are strike 152.7°, dip 61.4°, rake -4.8° and strike 245.0°, dip 85.8°, rake -151.3° respectively, and focal depth is 14.0km. The August 9, Xinjiang Jinghe M6.6 shock resulted from rupturing on a south-dipping thrust fault with left-lateral strike-slip. The parameters of two nodal planes are strike 100.6°, dip 27.5°, rake 114.1° and strike 259.3°, dip 65.1°, rake 78.0°, and the focal depth is 16.0km. The direction of two nodal planes is consistent with regional seismotectonic background.     

Source parameters of some large earthquakes in Northern Aegean determined by body waveform inversion

Average source parameters for three large North Aegean events are obtained from body wave inversion for the moment tensor. The parameters for the events are as follows: The events exhibit dextral strike-slip faulting with theT axis striking NS and nearly horizontal, implying extension in this direction. The focal mechanisms obtained are in agreement with the seismotectonic regime of the North Aegean. It is known that the region is tectonically controlled by the existence of the strike-slip Anatolian fault and its westward continuation in the Aegean, as well as the NS extension the whole Aegean area undergoes.The components of the moment tensor show that the region is dominated by compression in the EW direction which is encompassed by extension in the NS direction. All the events were found to be shallow (10 km) with a source time function of approximately 8 s duration and small stress drop values.The teleseismic long period verticalP-waves exhibited distortions, that could be attributed to lateral inhomogeneities in the source structure or more probably to a nonflat water-crust interface.     

Tectonic analysis of mine tremor mechanisms from the Upper Silesian Coal Basin

Fault network of the Upper Silesian Coal Basin (USCB) is built of sets of strike-slip, oblique-slip and dip-slip faults. It is a typical product of force couple which acts evenly with the parallel of latitude, causing horizontal and anti-clockwise movement of rock-mass. Earlier research of focal mechanisms of mine tremors, using a standard fault plane solution, has shown that some events are related to tectonic directions in main structural units of the USCB. An attempt was undertaken to analyze the records of mine tremors from the period 1992–1994 in the selected coal fields. The digital records of about 200 mine tremors with energy larger than 1×10⁴ J (M _L >1.23) were analyzed with SMT software for seismic moment tensor inversion. The decomposition of seismic moment tensor of mine tremors was segmented into isotropic (I) part, compensated linear vector dipole (CLVD) part and double-couple (DC) part. The DC part is prevalent (up to 70%) in the majority of quakes from the central region of the USCB. A group of mine tremors with large I element (up to 50%) can also be observed. The spatial orientation of the fault and auxiliary planes were obtained from the computations for the seismic moment DC part. Study of the DC part of the seismic moment tensor made it possible for us to separate the group of events which might be acknowledged to have their origin in unstable energy release on surfaces of faults forming a regional structural pattern. The possible influence of the Cainozoic tectonic history of the USCB on the recent shape of stress field is discussed.