Influence of fluids and magma on earthquakes: seismological evidence

In this paper, we present seismological evidence for the influence of fluids and magma on the generation of large earthquakes in the crust and the subducting oceanic slabs under the Japan Islands. The relationship between seismic tomography and large crustal earthquakes (M=5.7-8.0) in Japan during a period of 116 years from 1885 to 2000 is investigated and it is found that most of the large crustal earthquakes occurred in or around the areas of low seismic velocity. The low-velocity zones represent weak sections of the seismogenic crust. The crustal weakening is closely related to the subduction process in this region. Along the volcanic front and in back-arc areas, the crustal weakening is caused by active volcanoes and arc magma resulting from the convective circulation process in the mantle wedge and dehydration reactions in the subducting slab. In the forearc region of southwest Japan, fluids are suggested in the 1995 Kobe earthquake source zone, which have contributed to the rupture nucleation. The fluids originate from the dehydration of the subducting Philippine Sea slab. The recent 2001 Geiyo earthquake (M=6.8) occurred at 50 km depth within the subducting Philippine Sea slab, and it was also related to the slab dehydration process. A detailed 3D velocity structure is determined for the northeast Japan forearc region using data from 598 earthquakes that occurred under the Pacific Ocean with hypocenters well located with SP depth phases. The results show that strong lateral heterogeneities exist along the slab boundary, which represent asperities and results of slab dehydration and affect the degree and extent of the interplate seismic coupling. These results indicate that large earthquakes do not strike anywhere, but only anomalous areas which can be detected with geophysical methods. The generation of a large earthquake is not a pure mechanical process, but is closely related to physical and chemical properties of materials in the crust and upper mantle, such as magma, fluids, etc.     

Seismic tomography and anisotropy in the source area of the 2008 Iwate–Miyagi earthquake (M 7.2)

In order to understand the generation mechanism of the Iwate–Miyagi earthquake (M 7.2) that occurred on 14 June 2008 in Northeast Japan, we determined high-resolution 3D seismic velocity (Vp, Vs) and Poisson's ratio (σ) structures as well as P-wave azimuthal anisotropy in and around the source area using 246,623 P-wave and 189,153 S-wave high-quality arrival times from 1941 aftershocks of the 2008 Iwate–Miyagi earthquake as well as 4878 shallow and intermediate-depth earthquakes in this area recorded by the dense High-Sensitivity Seismic Network (Hi-net) and the seismic stations of Tohoku University. Our results show that hypocenters of the mainshock and three large aftershocks are located in a boundary zone where both seismic velocity and Poisson's ratio change drastically in a short distance. A zone with pronounced low-velocity and high Poisson's ratio is revealed in the lower crust and uppermost mantle under the source area, which may reflect the arc magma and fluids ascending from the upper-mantle wedge. The P-wave azimuthal anisotropy is complex under the study area, which may also indicate the complex crustal structure there. These results indicate that the generation of the 2008 Iwate–Miyagi earthquake was influenced by the ascending arc magma and fluids associated with the dehydration reactions of the subducting Pacific slab under Northeast Japan.     

Three-dimensional Mapping of Magma Source and Transport Regions from Seismic Data: The Mantle Wedge beneath Northeastern Japan

—We investigate the distribution of partial melt in island arc using the seismic velocity structure of the mantle wedge beneath northeastern Japan. The comparison of the seismic tomography with laboratory velocity data on a partially-molten mantle rock yields estimates of melting zones in three dimensions. We employ experimental data on the degree of partial melt in hydrous peridotite to give constraints on the melt fraction and temperature. Melting and magma-rich zones derived from the velocity structure coincide with observed low Q zones. The results of the three-dimensional mapping indicate that the source of magma in island arc is diapir-like melting patches localized within the low velocity zones of the mantle wedge. Extensive volcanic activity along the volcanic front is due to the presence of vast magma-rich zones just beneath the Moho. Those melting zones in the uppermost mantle may, in turn, cause melting of lower crustal materials and produce felsic magma. Melt appears to stay at and beneath the Moho, where crystallization fractionation may proceed. Melt exists at greater depths in the back-arc region, which may correlate with across-arc variations of chemical compositions of the volcanic rocks observed in northeastern Japan. We suggest that magma migration in the ductile lower crust may cause low-frequency microearthquakes, and magma penetration into the brittle upper crust may produce mid-crustal S-wave reflectors.     

龙门山断裂带多参数深部结构成像与地震成因研究

通过反演大量的纵、横波地震数据,获得了沿龙门山断裂带及周边区域的深部三维精细结构,结合前人二维大地电磁探测研究成果,提出龙门山断裂带地壳形变与深部速度结构和导电率不均匀性有关,探讨了2008年汶川和2013年芦山地震的诱发和产生与流体侵入及地壳形变的密切关系.本研究发现,2008年汶川地震发生在高速度、高泊松比和低电导率的区域,2013年芦山地震则位于高速度、低泊松比和低电导率的发震层.在上地壳中,四川前陆盆地的低速、低泊松比和低阻异常与松潘一甘孜地块的高速、高泊松比和高阻异常形成了鲜明的对比.在龙门山断裂带下方的两个低速和低阻块体,将龙门山断裂带分成南、中和北三段.我们的研究认为,这两个异常体与来自松潘甘孜地块的下地壳和(或)上地幔的局部熔融或流体侵入到龙门山断裂带的脆弱区有关.基于对汶川和芦山地震的余震分布特征及震源区的地震波速度、泊松比及电阻率参数分析,揭示了龙门山断裂带深部剧烈的地壳形变与流体应力积累对2008年汶川和2013年芦山地震的触发及其地震破裂过程具有重要的控制作用.     

华北地区壳内低速层与地壳流变的关系及其对强震孕育的影响

京津唐张地区普遍存在壳内低速层 ,鄂尔多斯块体内部没有发现低速层 ,壳内低速层的这种分布受新生代裂陷伸展的控制。华北地区中地壳下部和下地壳低速层是岩石塑性流变的结果 ,中地壳上部低速层是地壳裂陷伸展时形成的水平拆离带和韧性剪切带 ,岩石各向异性和流体作用可能是引起低速的原因。壳内软弱层 (低速和塑性流变层 )增强了块体层间的解耦作用 ,对地震孕育起着重要作用     

南北地震带地壳结构多参数成像及强震触发机制研究

通过反演由大量的纵、横波地震数据组成的综合数据集,获得了南北地震带地壳的多参数三维精细结构,探讨和分析了南北地震带的高地震活动性和强震频发的原因.成像结果表明,尽管1976年松潘一平武地震(M7.2)与2008年汶川地震(M8.0)以及2013年芦山地震(M7.0)均发生在高速、低泊松比异常区域,并且在其震源的下方均有一低速、高泊松比异常区域.我们认为,上述三个地震的触发与流体侵入导致的地壳形变之间有密切的联系.1955年炉霍地震(M7.4)和1973年康定地震(M7.1)均发生在鲜水河断裂带上,其震源中心区域表现为低速、高泊松比异常,可以解释为下地壳中的流体沿断层面上涌.在震源区的周边区域兼有高速、高泊松比异常,低速、高泊松比异常以及高速、低泊松比异常,可能分别与含流体的岩石、沿断裂带发育的变质岩以及坚硬的克拉通块体对应.流体的侵入不仅能够改变断层面上的应力情况,还能降低岩石骨架的岩石力学强度,进而触发地震.1970年云南通海大地震(M7.1)发生在哀牢山一红河断裂带附近的曲江断裂上,其震源处于高速度、低泊松比异常与低速度、高泊松比异常之间的边界区域,被认为是流体挤压后的应变能积累,最终导致脆性破裂,以至于发生地震.根据本次研究获得的多参数结构图像,结合前人的研究成果,我们认为南北地震带地壳强烈形变与流体侵入是造成该区域地震活动性较高及强震频发的两个主要因素.     

滇西地区壳幔解耦与腾冲火山区岩浆活动的深部构造研究

根据青藏东部边缘的深部地球物理资料,分析了滇西地区壳幔耦合和腾冲火山区岩浆活动的深部构造特征,确认了地幔各向异性与上地幔速度结构(包括P波速度和S波速度)的内在联系,指出产生这一结果的原因与以腾冲火山区为中心的地幔热物质上涌有关:上地幔顶部平均温度升高导致介质强度降低,在印支块体的侧向挤压或印缅块体的向东俯冲作用下发生韧性变形,造成滇西地区地幔各向异性的快波方向与青藏东部地壳块体的旋转方向不一致.此外,鉴于中下地壳低速层的横向非均匀性,估计韧性流动并非贯通青藏高原的东部边缘,而是被不同的构造块体和边界断裂限定在局部地区.总体而言,滇西地区下地壳的地震波速度和电阻率偏低,具备发生韧性变形的构造条件.作为地壳和上地幔之间的解耦层,它使得青藏东部地壳块体旋转产生的构造应力未能传输至上地幔.腾冲火山区的地壳结构与不同时期的岩浆活动有关,火山区东侧的高速结构代表了上新世时期火山通道内冷凝固结的岩浆侵入体或难以挥发的高密度残留物质,火山区西侧的低速结构反映了更新世以来持续至今的岩浆活动,壳内岩浆源主要分布在10～20km的深度范围内,横向尺度约为15～20km,有可能通过地壳深部的断裂与上地幔岩浆源区相连,估计腾冲火山区下方的岩浆活动将持续进行.     

腾冲火山区的地震层析成像及其构造意义

利用滇西南临时台网和固定台站的地震数据反演了腾冲及邻近地区的P波速度结构,着重分析了腾冲火山区和龙陵7级地震震源区的地壳结构特点.研究结果表明,腾冲火山区下方10～20 km深度范围存在明显的低速区,其横向尺度大约在20～30 km之间;推测这一低速区代表了仍处于活动状态的壳内岩浆源,热流通道有可能通过腾冲断裂延伸至地...     

渤海海峡及周边区域地壳结构的层析成像特征

渤海海峡是连接中国东部山东半岛和辽东半岛的重要途径,其跨海通道的地壳稳定性研究受到高度关注.本文利用地震层析成像方法重建三维P波速度模型,揭示了渤海海峡及周边区域地壳和上地幔的构造特征.结果表明,渤海海峡的速度结构存在明显的非均匀性,海峡北部地壳速度较高,结构较为完整,断层活动不明显,与现今较弱的地震活动相吻合,但是地壳底部存在低速薄层,它有可能成为地壳和上地幔之间的滑脱带,需要开展进一步的研究加以确认.相比之下,海峡南部地壳速度偏低,附近区域地震活动频繁,与张家口—蓬莱断裂带通过于此有着密切的联系,该断裂持续不断的地震活动对海峡南部的地壳结构产生了较大的影响.在渤海南部,郯庐断裂带东、西两侧的地壳结构明显不同,西侧速度偏高,东侧至渤海海峡速度偏低,这一特征可能与此地区广泛发育的断层和地震活动有关.另外,受华北克拉通破坏及地幔上涌的影响,渤海地区地壳深部和上地幔速度偏低,郯庐断裂带及渤海海峡附近显示出深部热流的活动迹象,反映了岩石圈减薄和软流圈的局部抬升.     

华北强烈地震深部构造环境的探测与研究

20世纪六七十年代以来, 华北地区发生了一系列强烈地震． 强烈地震的孕育、 发生和发展与深部构造密切相关． 近50年来, 我国地震科学领域在强烈地震的地震构造和深部环境方面开展了大量的研究． 深部地球物理探测和地震层析成像结果揭示了华北地区地壳结构的基本特征, 并在强烈地震发生的深部构造环境等问题上取得了重要进展． 本文在回顾华北地区地壳上地幔结构探测的基础上, 对1966年邢台M_S7.2, 1976年唐山M_S7.8, 1975年海城M_S7.3和1679年三河—平谷M8.0地震的地震构造和深部构造环境进行评述． 深部地球物理数据的综合分析表明, 震源下方的低速异常带, 高角度超壳深断裂, 地壳深浅构造的不一致, 偏低的上地幔顶部速度和局部隆起的莫霍界面, 是华北伸展构造区深部孕震环境的共同特征．      

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.     

Seismic images under the Beijing region inferred from P and PmP data

In this study a new tomographic method is applied to over 1500 high-quality PmP (Moho reflected wave) travel-time data as well as over 38,500 high-quality first P-wave arrivals to determine a detailed 3D crustal velocity structure under Beijing and adjacent areas. Results of detailed resolution analyses show that the PmP data can significantly improve the resolution of the model in the middle and lower crust. After the PmP data are included in the tomographic inversion, our new model not only displays the tectonic feature appeared in the previous studies, but also reveals some new features. The Zhangjiakou-Bohai Sea fault zone (Zhang-Bo zone) is imaged as prominent and continuous low-velocity (low-V) anomalies in the shallower crust, while in the middle and lower crust it shows intermittent low-V anomalies extending down to the uppermost mantle. Furthermore, the pattern of low-V anomalies is different along the Zhang-Bo zone from the southeast to the northwest, indicating that there exist large differences in the dynamic evolution of Taihangshan and Yanshan uplifts and North China depression basin. Prominent low-V anomalies are visible under the source area of the 4 July 2006 Wen-An earthquake (M 5.1), suggesting that the occurrence of the Wen-An earthquake is possibly related to the effect of the crustal fluids probably caused by the upwelling of the hot and wet asthenospheric materials due to the deep dehydration of the stagnant Pacific slab in the mantle transition zone. The fluids in the lower crust may cause the weakening of the seismogenic layer in the upper and middle crust and thus contribute to the initiation of the Wen-An earthquake. This is somewhat similar to the cause of the 1695 Sanhe-Pinggu earthquake and the 1976 Tangshan earthquake in the region, as well as the 1995 Kobe earthquake in Japan and the 2001 Bhuj earthquake in India.     

A study on 3-D velocity structure of crust and upper mantle in Sichuan -Yunnan region, China

Many evidences indicate that the collision of two plates deformed strongly the crust of the SYR, and the deformation has been continued up to the present. In addition, the SYR is in the south segment of the South-North Seismic Zone of China, which is one of the regions in the Chinese mainland, where the seismic activity is very high, and the strong earthquakes frequently occurred. Since the 1970s, a series of large earthquakes with magnitude M>7.0 occurred in SYR, such as the 1970 Tongha…     

The influence of crustal structure on compositions of subduction-related magmas

A survey of Sr isotopic ratios and other compositional features of subduction-related magma suites reveals significant correlations between these averaged parameters and characteristics of the underlying crust (i.e., thickness, composition, and age). These observations lead to the conclusion that crust and(or) mantle rocks in the hanging walls of subduction zones are involved in modification of primary mafic magmas (typically basalt or basaltic andesite). It is proposed that mafic magmas will stagnate within the crust or uppermost mantle where they may differentiate and react with wall rocks. The extent to which such processes manifest themselves will depend upon details of the local crustal structure. In particular, the composition and age of the crust will strongly influence such parameters as Sr, Nd and Pb isotopic compositions. Such data strongly indicate the involvement of crustal rocks in locales underlain by old sialic crust (e.g., central Andes). Depending upon the level of magma stagnation and evolution within the crust, different trends in isotopic composition may result. These isotopic trends may be enhanced by partial melting of the wall rocks to produce relatively silicic anatectic magmas, and locally they may reflect subduction of continental sediments. Interpretation of the isotopic data may be more ambiguous in locales underlain by younger and more mafic continental crust (Cascades, E Eleutians) and those underlain by oceanic crust owing to the similarity in isotopic composition of primary magmas and the latter crustal materials. Yet some degree of crustal involvement in magmatic evolution seems highly probable even in these more primitive terranes. Consequently, most island arc magmas, and especially those more evolved than basalt, are probably not primary in the sense that they do not represent direct melts of the upper mantle. Studies of arc volcanic rocks may yield misleading conclusions concerning processes of magma generation related to subduction unless evolutionary processes are defined and their effects considered. It appears that modern volcanic arcs provide a poor analog for models of early crustal development because the modern mantle-derived magmatic components are more mafic in composition than average continental crust.     

Fine S wave velocity structure beneath Iwate volcano, northeastern Japan, as derived from receiver functions and travel times

A genetic algorithm inversion of receiver functions derived from a dense seismic network around Iwate volcano, northeastern Japan, provides the fine S wave velocity structure of the crust and uppermost mantle. Since receiver functions are insensitive to an absolute velocity, travel times of P and S waves propagating vertically from earthquakes in the subducting slab beneath the volcano are involved in the inversion. The distribution of velocity perturbations in relation to the hypocenters of the low-frequency (LF) earthquakes helps our understanding of deep magmatism beneath Iwate volcano. A high-velocity region (dV_S/V_S=10%) exists around the volcano at depths of 2–15 km, with the bottom depth decreasing to 11 km beneath the volcano’s summit. Just beneath the thinning high-velocity region, a low-velocity region (dV_S/V_S=−10%) exists at depths of 11–20 km. Intermediate-depth LF (ILF) events are distributed vertically in the high-velocity region down to the top of the low-velocity region. This distribution suggests that a magma reservoir situated in the low-velocity region supplies magma to a narrow conduit that is detectable by the hypocenters of LF earthquakes. Another broad low-velocity region (dV_S/V_S=−5 to −10%) occurs at depths of 17–35 km. Additional clusters of deep LF (DLF) events exist at depths of 32–37 km in the broad low-velocity zone. The DLF and ILF events are the manifestations of magma movement near the Moho discontinuity and in the conduit just beneath the volcano, respectively.     

前郭MS5.8和松原MS5.7地震震源区地壳速度结构与孕震环境

为了解2013年11月23日吉林前郭M_S5.8和2018年5月28日吉林松原M_S5.7地震震源区的地壳速度结构及其孕震环境,本研究收集了东北地区CEA固定台站、NECESSArray流动台站和五大连池WAVESArrary流动台站的连续波形数据,使用背景噪声成像方法,获得了前郭地震和松原地震震源区地壳S波速度结构.结果显示,这两个地震主要发生在相对高波速异常区,而震源区下方存在明显低波速异常,且该低速层呈北北东走向,并有向北延展的趋势.结果表明,震源区相对高波速异常区易于积累能量孕育地震,而震源区下方的低波速异常可能代表了流体作用.这种流体作用会降低断层面有效正应力从而触发地震.这种流体作用可能和太平洋板块深俯冲至我国东北地区（包括前郭地震和松原地震震源区）下方的地幔转换内形成"大地幔楔"结构与动力学密切相关.在"大地幔楔"结构中,由于地幔转换带中滞留板块脱水作用和地幔角流作用,容易形成热湿物质上涌,进而导致松辽盆地中北部岩石圈物质拆沉、地幔流体入侵至中下地壳、然后作用于断裂带,从而导致了两个震源区中强度地震的发生.     

华北地区地震和深部构造关系及其破裂机制探讨

华北地区的地震活动格外活跃, 其地壳深部构造的独特性决定了其地震活动性。 多条地学断面和地震剖面等地球物理研究表明, 华北地区的地震带(即活动构造带)一般均存在低速体或上地幔上隆; 华北强震区的地球物理详细勘测发现, 强震区也都存在低速体或上地幔上隆。 应用应力摩尔圆和库伦破裂准则, 利用上涌模型和地壳减薄模型研究地壳岩石破裂机制, 结果表明, 存在上地幔上隆和低速体的地区的地壳岩石比其它地区更易于发生破裂。 华北地区的地壳深部构造决定了地震的空间分布格局。     

川滇地区地壳上地幔三维速度结构研究

根据云南和四川地震台网174个台站记录的4625个区域地震初至Ｐ波和Ｓ波走时资料，并结合其它深部地球物理资料，确定了川滇地区地壳上地幔三维速度结构．在上地壳速度异常分布中，四川盆地为正异常，川西高原为负异常，龙门山断裂带为正、负异常的边界．龙门山断裂、鲜水河断裂以及红河断裂等，在下地壳和上地幔的速度异常中仍显示出构造分界特征，说明它们可能穿透了莫霍界面．腾冲火山区和攀西构造带在50km深度上呈现负速度异常，与上地幔温度和物质组成的差异相联系．川滇地区地壳结构的总体特征是:地壳和上地幔的低平均速度，地壳厚度变化剧烈，地壳和（或）上地幔存在高导层、高热流值．这些同印度板块与欧亚板块碰撞的构造背景有关．川滇菱形块体在地壳内总体上为正常或正异常速度，而其边界的深大走滑断裂存在负速度异常，它有助于地壳块体沿断裂的侧向挤出．在主要的地震带上，中下地壳的负速度异常与地震活动性相关．多数强烈地震发生在具有正速度异常或正常速度分布的上中地壳深度上，而其下方则通常是负速度异常带．      

Seismic imaging of the entire arc of Tohoku and Hokkaido in Japan using P-wave, S-wave and sP depth-phase data

In order to better understand seismic structure and seismotectonics of the entire arc of Tohoku and Hokkaido in Japan, we combined arrival time data from earthquakes beneath Tohoku and Hokkaido land areas, and beneath the Pacific Ocean to determine the three-dimensional (3D) velocity structures (Vp and Vs) under the entire Northeast (NE) Japan-Kuril arc. We adopted 176,431 P-wave and 110,953 S-wave arrival times, from 5123 local earthquakes, and 2843 sP depth-phase data from 385 events that occurred beneath the Pacific Ocean. The 385 suboceanic events were accurately relocated by using P-wave, S-wave and sP depth-phase arrival time data jointly. The obtained results confirmed the major features delineated by previous studies and revealed some new features of the structural heterogeneity beneath NE Japan and the Kuril forearcs. High-velocity anomalies of the cold subducting Pacific slab and low-velocity anomalies in the hot mantle wedge were imaged clearly. Strong lateral heterogeneities were revealed on the upper boundary of the Pacific slab under the forearc region, which showed a good correlation with the spatial distribution of large interplate earthquakes. These results indicated that strong coupling sections (or asperities) and weak-coupled or decoupled patches might exist along the upper boundary of the Pacific slab. Widespread low-velocity anomalies were visible in the forearc mantle above the subducting Pacific slab, which might reflect serpentinization of the forearc mantle associated with the dehydration process of the subducting slab. Our results also showed a general tendency for seismic coupling in the asperities to be located around low-velocity areas on the slab boundary under the suboceanic region.     

A study on 3-D velocity structure of crust and upper mantle in Sichuan-Yunan region, China

Based on the first arrival P and S data of 4 625 regional earthquakes recorded at 174 stations dispersed in the Yunnan and Sichuan Provinces, the 3-D velocity structure of crust and upper mantle in the region is determined, incorporating with previous deep geophysical data. In the upper crust, a positive anomaly velocity zone exists in the Sichuan basin, whereas a negative anomaly velocity zone exists in the western Sichuan plateau. The boundary between the positive and negative anomaly zones is the Longmenshan fault zone. The images of lower crust and upper mantle in the Longmenshan fault, Xianshuihe fault, Honghe fault and others show the characteristic of tectonic boundary, indicating that the faults likely penetrate the Moho discontinuity. The negative velocity anomalies at the depth of 50 km in the Tengchong volcanic area and the Panxi tectonic zone appear to be associated with the temperature and composition variations in the upper mantle. The overall features of the crustal and the upper mantle structures in the Sichuan-Yunnan region are the lower average velocity in both crust and uppermost mantle, the large crustal thickness variations, and the existence of high conductivity layer in the crust or/and upper mantle, and higher geothermal value. All these features are closely related to the collision between the India and the Asia plates. The crustal velocity in the Sichuan-Yunnan rhombic block generally shows normal value or positive anomaly, while the negative anomaly exists in the area along the large strike-slip faults as the block boundary. It is conducive to the crustal block side-pressing out along the faults. In the major seismic zones, the seismicity is relative to the negative anomaly velocity. Most strong earthquakes occurred in the upper-mid crust with positive anomaly or normal velocity, where the negative anomaly zone generally exists below. Foundation item: National Scientific and Technological Development Program (95-973-02-02), the Climb Program (95-S-05-01) of National Scientific and Technological Ministry of China, and the State Natural Sciences Foundation of China (49874021). Contribution No. 02FE2004, Institute of Geophysics, China Seismological Bureau.