Scattered Surface Wave Energy in the Seismic Coda

One of the many important contributions that Aki has made to seismology pertains to the origin of coda waves (Aki, 1969; Aki and Chouet, 1975). In this paper, I revisit Aki's original idea of the role of scattered surface waves in the seismic coda. Based on the radiative transfer theory, I developed a new set of scattered wave energy equations by including scattered surface waves and body wave to surface wave scattering conversions. The work is an extended study of Zeng et al. (1991), Zeng (1993) and Sato (1994a) on multiple isotropic-scattering, and may shed new insight into the seismic coda wave interpretation. The scattering equations are solved numerically by first discretizing the model at regular grids and then solving the linear integral equations iteratively. The results show that scattered wave energy can be well approximated by body-wave to body wave scattering at earlier arrival times and short distances. At long distances from the source, scattered surface waves dominate scattered body waves at surface stations. Since surface waves are 2-D propagating waves, their scattered energies should in theory follow a common decay curve. The observed common decay trends on seismic coda of local earthquake recordings particular at long lapse times suggest that perhaps later seismic codas are dominated by scattered surface waves. When efficient body wave to surface wave conversion mechanisms are present in the shallow crustal layers, such as soft sediment layers, the scattered surface waves dominate the seismic coda at even early arrival times for shallow sources and at later arrival times for deeper events.     

Source side seismic tomography(3STomo):A novelmethod to image the subsurface structure beneath seismically active region

We propose a novel seismic tomography method,Source Side Seismic Tomography(3STomo),which is designed particularly to image the subsurface structure beneath seismically active regions.Unlike the teleseismic tomography,in which the data are relative traveltime residuals between closely spaced stations for each teleseismic event,3STomo uses relative traveltime shifts between earthquakes within the study region for each distant station.Given the relatively evener distribution of global seismic stations,this method has unique advantages for imaging the structure beneath regions that have numerous earthquakes but lack of dense seismic stations,for example,some subduction zones and spreading ridges in the ocean.In addition,3STomo has potentially better vertical resolution at shallow depths than the traditional teleseismic tomography.The effect of the inaccurate source parameters on its resolution can be minimized by using depth phases and the technique of joint source and structure inversion.Numerical experiments and application to Luzon Island,Philippines show that 3STomo can be a valuable tool to investigate the subsurface structure beneath some areas where the traditional method cannot be applied to,or at least it can be used as a complementary component of conventional teleseismic tomography to obtain better back-azimuth coverage and achieve higher resolution at shallow depths in the inversion.     

南中国海区域海啸预警中心地震监测系统

简要介绍了南中国海区域海啸预警与减灾系统的建设和发展历程,同时重点阐述了地震监测系统构成及其基本功能。作为重要组成部分,地震监测系统通过地震数据的实时汇集、存储、自动处理和分析,并结合人机交互方式实现了地震定位、震源机制解和有限断层模型反演。实际应用表明,地震监测系统对全球6.0级以上地震定位时间不超过8 min,在震后10—15 min内完成W震相方法快速反演海底强震震源机制解,在震后短时间内完成有限断层模型反演,为海啸预警提供快速、准确、可靠的地震基本参数和震源特征参数。      

天祝-古浪地区双差地震层析成像与古浪地震构造探讨

利用观测走时、走时残差和波形互相关方法获得走时残差,采用双差地震层析方法对1996年6月1日天祝5.2级地震后布设的微震台网资料进行处理,获得了台网下方上地壳三维速度结构的基本特征和地震绝对位置的空间展布特征,并由此勾画出区域活动断裂在深部的延伸情况.地震层析结果表明,研究区上地壳速度结构在纵向和横向上具有明显的非均匀...     

一种低频地震事件的检测方法在汶川地震和芦山地震前的初步应用

本文发展了一种单台检测低频地震事件的方法, 并用其检测2008年汶川M_S8.0地震、 2013年芦山M_S7.0地震震前的低频事件。 使用汶川地震、 芦山地震前周边的宽频带台站数月的垂直向连续波形, 经过2~8 Hz带通滤波、 窗长为10 s的包络平滑、 窗长为20 min的中值滤波, 去除近场天然地震及远场大震与低频噪声的影响后, 发现这两个大地震前有少数台站在连续数天的时间内, 出现持续时间长(数十分钟到数小时)、 包络中值明显大于背景噪声的现象, 其波形记录及包络特征与非火山震颤(NVT)的波形记录及包络特征具有较好的相似性, 我们认为这是汶川地震和芦山地震前的疑似慢地震事件。 初步应用表明, 本文所用的方法对于从宽频带连续波形资料中提取NVT信号是可行的。     

Tomographic imaging of P wave velocity structure beneath the region around Beijing

The three-dimension crustal and upper mantle structures in the region around Beijing were studied by seismic tomography. We used the P wave arrival times from local and teleseismic events. These events were recorded by 250 stations of the North China Seismic Array and 108 stations of the Beijing Telemetry Seismic Network. 118 869 P wave arrivals from 10 285 local events and 12 189 P wave arrivals from 107 teleseismic events were used in the inversion. We obtained the 3-D P wave velocity structure of the cru...     

多震相在近震体波层析成像中的应用

近震体波层析成像方法利用地震波走时来反演地下结构,这需要观测到时准确。在实际情况中,观测到时或多或少包含一定的误差,包括人为误差、仪器误差（如钟差）等。本文提出了一种新的近震体波层析成像方法,来消除随时间积累的误差对反演结果的影响,并验证了该方法的正确性和可行性。     

New advances of seismic tomography and its applications to subduction zones and earthquake fault zones: A review

Abstract There have been significant advances in the theory and applications of seismic tomography in the last decade. These include the refinements in the model parameterization, 3-D ray tracing, inversion algorithm, resolution and error analyses, joint use of local, regional and teleseismic data, and the addition of converted and reflected waves in the tomographic inversion. Applications of the new generation tomographic methods to subduction zones have resulted in unprecedentedly clear images of the subducting oceanic lithosphere and magma chambers in the mantle wedge beneath active arc volcanoes, indicating that geodynamic systems associated with the arc magmatism and back-arc spreading are related to deep processes, such as the convective circulation in the mantle wedge and deep dehydration reactions in the subducting slab. High-resolution tomographic imagings of earthquake fault zones in Japan and California show that rupture nucleation and earthquake generating processes are closely related to the heterogeneities of crustal materials and inelastic processes in the fault zones, such as the migration of fluids. Evidence also shows that arc magmatism and slab dehydration may also contribute to the generation of large crustal earthquakes in subduction regions.     

Crustal structure beneath Tien Shan orogenic belt and its adjacent regions from multi-scale seismic data

As one of the world's most active intracontinental mountain belts, Tien Shan has posed questions for researchers regarding the formation of different tectonic units and active shallow seismicity. Here, we used a huge data set comprising of 7094 earthquakes from local, regional and teleseismic seismic stations. We used waveform modeling and multi-scale double-difference earthquake relocation technique to better constrain the source parameters of the earthquakes. The new set of events provided us with better initial earthquake locations for further tomographic investigation. We found that reverse-faulting earthquakes dominate the whole study area while the fault plane solutions for earthquakes beneath the northwestern Tarim Basin and the Main Pamir Thrust are diverse. There is a low-velocity anomaly beneath Bashkaingdy at depth of 80 km, and high-velocity anomalies beneath central Tien Shan at shallower depths. These observations are the keys to understand the mechanism of Tien Shan's formation because of Tarim Basin northward and Kazakh Shield's southward subduction in the south and north respectively. Velocities beneath western Tien Shan are relatively high. We thus infer that the Western Tien Shan is relatively less deformed than the eastern Tien Shan primarily due to a relatively brittle mantle.     

雷琼火山与华夏块体的远震P波走时成像

利用广东及邻省的地震台网共94个地震台站记录的2009—2015年的256个远震波形数据,采用基于波形互相关的自适应迭代技术获得了18968个可靠的P波相对走时残差数据,进而应用FMTT(Fast Marching Teleseismic Tomography)地震层析成像方法确定了雷琼火山与华夏块体深至450km的三维P波速度结构.结果显示,雷琼火山下方100~450km深度显示明显低速异常,该低速异常随深度增加向华夏块体下方扩展.研究认为该低速异常由海南地幔柱热物质上涌造成,雷琼火山为板内热点火山,这与已有体波成像及接收函数成像等结果一致.根据该低速异常的整体特征,推测海南地幔柱向北东方向倾斜,并可延伸至450km以下.华夏块体下方的低速异常,可能是由海南地幔柱热物质上涌过程中不断向北东方向侵蚀而造成.结果同时显示研究区大陆至海域的岩石圈厚度可能具有由厚变薄的过渡特征.     

Lithospheric structure of Southeast Anatolia from joint inversion of local and teleseismic data

A joint tomographic inversion of local and teleseismic arrival times recorded at 41 seismic stations in southeast Anatolia is conducted to study the 3-D lithospheric velocity structure and its relation to the prevailing tectonic processes. A total of 21300 arrivals from local and teleseismic events are used in the final inversion. The tomographic model reveals prominent lower crustal/uppermost mantle low-velocity anomalies. High-velocity zones are imaged in the western part of the study area. The background seismic activity occurs mainly at the low-velocity areas and to a lesser extent in some high-velocity zones. Large crustal earthquakes occur in average velocity zones, but not in high-velocity areas that can resist stress. Results of the checkerboard resolution test indicate the reliability of the obtained images; while the large hit counts at most depth slices denote reasonable ray-path coverage for most parts of the study area. The obtained velocity anomalies are generally consistent with many previous geophysical measurements and give much deeper understanding of the current seismotectonic processes occurring in the region.     

用接收函数方法研究汶川和芦山地震之间未破裂段的地壳结构

芦山与汶川地震之间存在约40 km的地震空区.震源区和地震空区的深部构造背景的研究对深入了解中强地震的深部孕育环境及地震空区的地震活动性具有重要科学意义.利用本小组布设的15个临时观测地震台以及21个芦山科考台站和21个四川省地震局固定台站记录的远震数据,用H-K叠加方法得到各个台站的地壳厚度和平均泊松比,并构建了接收函数共转换点(CCP)偏移叠加图像以及反演得到台站下方的S波速度模型.我们的结果揭示了震源区和地震空区地壳结构特征差异:(1)汶川震源区的地壳平均泊松比为~0.28;芦山震源区为~0.29;而地震空区处于泊松比变化剧烈的区域;(2)汶川地震与芦山地震的震源区以西下方的Moho面呈现深度上的突变(这与前人的研究成果基本一致),分别从~44 km突变到~59 km,~40 km突变到~50 km,而地震空区地壳平均厚度呈现渐变性变化;(3)地震空区Moho面下凹且具有低速的上地壳.综合一维S波速度结构和H-k以及CCP的初步结果,这可能显示汶川地震的发震断裂在深部方向上向西倾斜并形成切割整个地壳的大型断裂;芦山地震则可能是由于上、下地壳解耦引起的;而地震空区处于两种地震形成机制控制区域的过渡带中.     

The influence of one-dimensional velocity sections on determining the key parameters of the earthquake sources in Azerbaijan

The paper addresses the construction of one-dimensional (1D) velocity models in the seismogenic regions of Azerbaijan taken individually and the analysis of implications of these models for estimating the key parameters of earthquake sources in Azerbaijan. We considered and analyzed the seismological data from the local earthquakes, the arrival times of the P-, P-g, Pn-, S-, Sg-, and Sn-waves recorded by the network of telemetry stations during the period from 2005 to 2014 with ml ≥ 2.5. For constructing the models, we used the VELEST program which calculates 1D velocity models from travel times of seismic waves. As a result, the 1D models were built for ten regions of Azerbaijan; the key parameters of the hypocenters of the earthquakes were recalculated; and the corrections to the body-wave arrival times at the observation stations were obtained, which increased the accuracy of locating the hypocenter of earthquakes.     

藏南地区中深源地震精确定深研究及其地球动力学含义

利用2004～2005年Hi-CLIMB计划架设在藏南日喀则附近由28台宽频带数字地震仪组成的二维台阵的地震记录,应用双差层析定位方法(TomoDD)对454个区域地震进行了精确重定位,共确定了340个地震的精确位置.发现区域震源深度呈明显的双层分布,其中有21 个地震的深度在50～80 km之间,与该地区的莫霍面深度相近.通过不同深度震源的理论地震图与实际地震图对比的方法,发现震源位置位于莫霍面上下地震图的震相存在明显差异,从而给出了存在地幔地震的直接观测证据.这一发现,对长期争论的地幔地震是否存在及大陆岩石圈流变结构的“三明治”模型有重要参考意义.     

Three-dimensional velocity inversion of teleseismic P-wave residuals from the telemetric network of the G.D.R.

Residuals of teleseismic P-wave arrival times at the CIPE network were used to estimate the three-dimensional velocity structure of the upper mantle by means of the inversion procedure described by Aki et al. (1977). Two models with different depth ranges (233 and 500 km, respectively) have been determined. They refer to the upper mantle beneath the Central European part of the Variscian fold-belt and its transition to the East European platform north of it. Both models reveal a lateral velocity increase within the upper mantle. It is directed to the northeast within the lithosphere and at least a part of the asthenosphere, and in a more northerly direction at greater depth. Furthermore, there seems to exist a local inhomogeneity west of Prague at a depth of some 200 km.     

应用地震层析成象技术研究深部速度结构

罗马尼亚Ｖａｎｃｅａ地震区是大陆上发生与板块磁撞和削减有关的中深部地震活动的地区之一。本文介绍了应用地震层析成象方法研究该地区深部速度结构成果。在研究中使用了地方和区域地震所记录的４３３个浅源和中部地震的到时资料反演求解深至２００ｋｍ的三维速度结构,在走时和射线路径的计算中利用了有效的三维射线跟踪技术,在反演中采用ＬＳＱＲ算法,高分辨率的地震层析图象揭示了速度结构的广泛不均匀性,结果表明,地震层析     

上海及其邻近地区三维地壳结构层析成像

利用上海及其邻近地区地震台网的地震记录及人工爆破资料,采用地震层析成像方法反演研究区地壳三维P波、S波速度结构模型。成像结果表明:不同深度的P波、S波速度扰动呈北西向展布特征;地震大多数分布在低速区的块体内,说明地震不仅与断裂活动有关,有可能还与物质结构和性质有关;上海及其邻近地区上地壳速度结构,断裂活动和物质性质三者之间存在密切关系。     

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.     

考虑到时误差的地震定位算法及其在四川地区2001-2008年地震定位的应用

精确的地震位置对于地震活动性、地震层析成像和地壳应力场反演具有相当重要的意义,对于地震速报也具有重要的应用价值。将观测到时的不确定性、台站高程、地震震源深度进行约束的同时,根据反演理论给出了地震震源位置精确估计和误差估计的方法。该算法联合考虑Pg波、Sg波、Pn波和Sn波的到时进行反演,数据量的增加可以增强地震位置的准确性,并可同时应用于地方震和区域地震。采用模拟数据对该地震定位算法进行检验发现,该算法在观测数据的不确定性不等时明显优于其他方法。将该算法应用于四川地区2001-2008年间的地震定位,得到的地震位置更加符合地震的丛集性并集中于断裂带附近。这些结果为四川地区的地震活动性、断层构造以及地震层析成像研究打下了基础,并且为汶川地震之前的地震活动前兆研究也提供了有益帮助。     

Pd-τc相容性检验方法在触发事件判别分析中的应用

结合福建地区地震预警系统的运行情况和处理能力, 利用日本KiK-net台网和福建地震台网记录的天然地震事件和人工爆破事件实测数据, 探讨了P_d-τ_c相容性检验方法在触发事件类型实时判别中应用的可行性． 结果表明, P_d-τ_c相容性检验方法能够较好地区分大震级地方震与误触发事件, 也能对远震事件进行粗略判别, 但对于小震级地方震和人工爆破事件则无法有效判别． 随后, 以P_d-τ_c相容性检验结果为基础, 本文设定强度参数P_d=10-3 cm为界线, 即可从众多触发事件中甄选出可能产生破坏的大震级地方震事件． 本文方法有助于提高触发事件类型判别结果的可靠性、 预警系统的处理效率和系统产出信息的准确度．