Are there links between the localized deformation in the northeastern section of Wenchuan earthquake source zone and the May 25，2008 Qingchuan strong aftershock?
-
摘要: 2008年汶川地震之后，通过InSAR观测到青川县木鱼镇附近存在一个长约为15 km、宽约为10 km、地表位移数十厘米的局部形变区。前人分析认为，该形变区是由MS6.4的青川强余震造成的，但拟合地表形变数据所采用的震源深度和震源机制解与地震学反演的结果具有较大差异。本文利用远震体波和瑞雷波振幅谱进一步测定了青川强余震的震源深度和震源机制解，计算了此次事件造成的地表位移场，认为青川强余震并非造成木鱼镇地区局部形变的直接原因，并讨论该局部形变区可能的成因。Abstract: A localized deformation zone near Muyu town of Qingchuan county was observed on the InSAR interferograms of the great 2008 Wenchuan earthquake. The Muyu deformation zone was about 15 km long and 10 km wide, with line of sight displacement up to dozens of centimeters. Some researchers proposed that the deformation zone is caused by a strong aftershock （MW6.1） in Qingchuan county on May 25, 2008, and they adopted seismic source parameters very different from seismic inversions. To verify reliability of the seismic source inversions, we use teleseismic body wave and local Rayleigh wave spectra to constrain focal depth and fault plane of the strong aftershock. Thereafter, we compute ground displacement and proposed that the aftershock was not the direct cause of the localized deformation zone. Then we discuss some candidate mechanisms which might explain the localized deformation.
-
Keywords:
- focal depth /
- source mechanism /
- ground deformation /
- InSAR
-
引言
夏威夷群岛的火山喷发记录始于1790年，其中的大部分火山事件，都可以定位在群岛的东南方向，也就是地质年代形成最晚的夏威夷大岛（Mooney，2018）。该岛共包含五座盾状火山，其中基拉韦厄火山的持续喷发导致其成为全球火山研究的热点。Nettles和Ekström （2004）重新分析了1975年卡拉帕纳（Kalapana）地震的长周期数据，观察到该地震的震源机制解显示为向陆地方向倾覆的冲断层；Wright和Klein （2006）通过研究夏威夷活火山的岩浆运移通道，认为岩浆通道在30 km以下是垂直的，沿海洋地幔边界是近似水平的，在基拉韦厄顶峰下接近垂直（Wright，Klein，2006）；Benz等（2002）和Lin等（2014）在基拉韦厄火山地区通过反演地震速度结构，发现在基拉韦厄火山埋深4—12 km处孕育了一个8 km宽的复杂高波速管状岩体，岩浆通道近乎垂直，火山东部裂谷带下方8—11 km处存在异常体，具有低速纵波、低速横波和高纵横波速比的特点，认定其为岩浆囊（Lin et al，2014）；Lei和Zhao （2006）将穿过外核的震相用于新的全球三维层析成像，观察到夏威夷的地幔柱成像结果表现为地表到核幔边界间的连续低速异常。
在夏威夷深部构造研究方面，科研人员大多采用地震学手段，然而，大地电磁测深方法在深部构造探测特别是火山研究中具有独特的优势。地下电阻率对介质的温度、含水饱和度和熔融状态非常敏感，而地下电性的变化通常是由岩浆的喷出和补充等循环变化造成的，因此电磁监测对研究低阻岩浆囊的动态变化较为灵敏。Constable和Heinson （2004）在海底布设了7个海洋大地电磁仪器，反演结果显示，在地幔柱顶部150 km处出现一个约为10 Ω·m的低阻体，且地下的导电性是各向异性的，垂直电导率高于水平电导率；汤吉等（1997，2001）观测到在天池火山区12 km深处存在岩浆囊，并且在不同剖面上表现出不同形态，其团队之后在阿尔山活火山区进行大地电磁观测，反演结果显示，新、老两条火山带的岩浆通道在深部可能是同源的（汤吉等，2005）；仇根根等（2014）在天池火山区同一地点进行大地电磁反演，得到了浅部发育的多个岩浆囊的位置及其发展趋势；李世文等（2020）通过大地电磁反演得出东北地区新生代火山区的三维电性结构，认为在高阻岩石圈下，存在多个与新生代火山相对应的局部低阻异常；Aizawa等（2011）在日本樱岛进行连续大地电磁测量，认为岩浆挥发物与地下水之间的混合作用可能会导致电阻率在海平面附近的某个深度发生变化；Díaz等（2012）通过分析大地电磁数据和岩石特征，推断出拉斯卡尔火山和普那火山下存在向南部延伸的大规模高电导率带，安第斯山脉南部存在两个不同深度和组成的岩浆储层；Abdallah等（2020）对库车火山群进行大地电磁测量，得到了该火山群下方的三维电阻率结构以及浅部热液系统与深部岩浆通道之间的联系。这些研究成果充分证明了大地电磁方法在探测低阻火山岩浆囊方面拥有的独特优势。
因此，大地电磁测深方法可用于探测火山深部电性结构，了解地下岩浆囊的构造赋存形态以及岩浆上升通道及深度展布范围，为岩浆活动和火山的危险性预测提供电性数据支持。然而，我们对夏威夷岛的深部构造特别是电性结构认识还很不充分，因此本文拟对夏威夷地区的深部电性结构进行研究，以期为认识夏威夷岛的深部结构提供电性依据。
1. 野外数据处理与分析
1.1 数据来源
文中的大地电磁数据来源于加拿大凤凰公司在夏威夷岛上布设的长期观测站（图1），位置为（19.8°N，155.6°W），测点海拔为1 665 m。台站采集到的自然发生的地磁场和感应地球电流的频率范围为10 400—0.005 5 Hz，通过趋肤深度公式进行粗略计算，可知反演深度可达100 km，100 km深度以下的地层电性信息由于超低频观测数据质量太差而参考价值不高。
![]()
图 1 夏威夷岛的构造背景及测点位置Figure 1. The tectonic background and distribution of the MT station of Hawaii Island1.2 近期地震与火山事件
夏威夷群岛地震及火山活动频发，参考中国地震台网中心国家地震科学数据中心（2021）发布的全球地震动态，可知自2018年5月至2019年4月大约一年的时间里，夏威夷大岛发生的M≥5.0的火山喷发地震高达54次。在大岛上标注出54次地震的震中位置（图2），可以看出地震非常频繁并且集中，震中大多汇聚在基拉韦厄火山附近。表1为部分较大震级的地震群的活动记录，这些震群的地震类型可能包括了火山构造地震、低频地震以及颤动（陈洪洲等，2004）。这些地震事件或发生在裂谷系统下的浅部（深约1—4 km），或沿着火山口下垂直的岩浆通道附近发生，且多以震群形式出现，有空间集中和震级较小的特征（傅征祥，1991）。
![]()
图 2 2018年5月至2019年4月地震的震中位置Figure 2. The epicenters of earthquakes from May 2018 to April 2019表 1 夏威夷岛部分较大震级的地震事件Table 1. Some earthquakes with large magnitude on Hawaii Island发震日期 北纬/° 西经/° 深度/km mb 年-月-日 时:分 2018-05-03 20:30 19.34 155.07 6 5.1 2018-05-04 22:32 19.31 155.00 2 6.9 2018-06-08 12:44 19.41 155.28 −1 5.2 2018-07-02 11:24 19.39 155.27 −1 5.3 2019-03-13 10:55 19.33 155.20 7 5.5 2019-04-14 03:09 19.47 155.79 13 5.3 1.3 数据质量分析
本文中数据的处理及反演使用的是MTPioneer软件。2018年7月15日的原始资料2018715M，2018715N，2018715O的视电阻率曲线如图3所示。其中：12个月份对应记为1—9，A，B，C；日期最后一位字母A-X分别对应了一天中的24个小时。2018年7月15日13时26分，夏威夷岛上发生了mb5.2地震，图3b的2018715N表示地震发生时记录到的视电阻率Rxy，低频段的数据质量很差，不能参与反演，图3a，c中2018715M和2018715O分别表示震前、震后记录到的视电阻率Rxy，其中中频部分的数据质量较好，低频的部分频点出现了飞点现象，高频部分出现整段的上移。由于本文主要关注中深层的电性信息，因此选择了频率小于629 Hz的频点进行反演。反演之前对部分受干扰严重的数据进行了编辑、圆滑、剔除等操作，使其达到反演要求。
![]()
图 3 原始数据中视电阻率Rxy的质量日期最后一位字母A−X分别对应一天中的24时，下同Figure 3. The quality of resistivity Rxy for the original dataThe last letters A-X of the date correspond to 24 hours in a day, the same below（a） 2018715M；（b） 2018715N；（c） 2018715O2. 大地电磁一维反演与讨论
对大地电磁数据进行自适应正则化一维反演，获得了观测台站下方深部电性结构及其长期变化特征。本文反演所用的阻抗数据均为Rxy，反演过程中，设置截断误差为1，最大迭代次数为50，反演初始模型为均匀半空间。
2.1 地震及火山活动前后电性特征分析
许多研究表明，夏威夷的地震活动直接或间接地与火山活动过程有密切关系。图4给出了三种典型火山地震发生前后的地下电性变化曲线，深度范围大约为0.5—130 km。给视电阻率加5%的误差，可以看出地震前后两条曲线基本无重叠部分，证明了视电阻率随时间变化的可行度。由图4可见，随着深度的增加，电阻率从地表向深处开始快速下降，到了2—3 km深度处，电阻率降低到最小值约20 Ω·m，与实测的熔融玄武岩熔岩的电阻率（大约为40 Ω·m）相近（Bartel et al，1983），推测其可能为低阻岩浆囊位置之一。在3 km深度以下，地震前后的电性差异出现明显变化，而且变化趋势相似，即电阻率先升高后降低，到了约100 km深度处逐渐趋于稳定的低阻异常，推测该深度下可能存在复杂的岩浆结构。
![]()
图 4 火山地震发生前后地下电性变化示意图图（b）−（d）线条上的短线为误差棒，误差棒的长短表示误差大小（a） 无地震发生；（b） 2018年5月3日；（c） 2018年5月4日；（d） 2018年6月8日Figure 4. Schematic diagram of underground electrical changes before and after volcanic earthquakesThe short line on the curves of Figs.（b）−（d） is error bar，and the length of the bar represents the error size（a） No earthquake；（b） May 3，2018；（c） May 4，2018；（d） June 8，2018火山活动与深部岩浆运输和物质熔融存在着密不可分的关系，火山喷发和强烈的浅层地震活动诱发了地下电阻率的变化，这是因为裂缝的发育会导致岩石中流体含量及其连通性发生变化，而岩浆和气体上升迁移又会导致岩石破裂。一般来说，熔融或部分熔融后的岩浆的电阻率值很低，地下电性变化通常对应着岩浆物质的喷出、补充等构造变化，但岩浆的上升通道和储存路径等细节却不明确。图4a表示2019年9月1日在没有地震发生的情况下，相近三个小时内的地下电性变化曲线，电阻率有变化但不明显，由此可证实电阻率变化的确与地震相关；图4b表示2018年5月3日地震发生前后深部电性变化，推测可能是由于火山型地震改变了岩浆输送通道，部分岩浆被挤压进入其它通道，导致地震后的电阻率大于地震前的电阻率；图4c表示2018年5月4日基拉韦厄南翼mb6.9地震前后的地下深部电性变化，可见其电性变化特征从浅至深在不同深度上都与2018年5月3日有所差异，随着深度的增加，电性差异逐渐增大，特别到了20 km以下，火山喷发后的电阻率明显高于喷发前的电阻率，可能是由于该层位大量的岩浆被抽送或运移至东南方向，导致了电阻率的升高；图4d表示2018年6月8日地震发生前后的电阻率变化，观测点下方100 km深度处，地震前后的电阻率变化差异不大，可能是此次地震对地下地层结构的影响较小，岩浆的喷出和补充变化不明显或者喷出与补充达到了一种动态平衡。Aizawa等（2011）通过对比在日本樱岛采集的视电阻率时间序列和多种地球物理数据，认为视电阻率的长期变化与潮汐、降水量之间无显著关联，故本文不考虑潮汐和降水对夏威夷岛地下视电阻率变化的影响。
2.2 夏威夷地区深部长期电性变化特征
为分析大岛地区较长时期内的深部电性变化特征，本文选择了连续18个月（2018年4月至2019年9月中每月第4日，或距4日最近一天）采集的大地电磁数据进行反演，选择典型的7条反演曲线进行展示（图5）。Díaz等（2020）认为复杂的岩浆系统具有多层次的储集能力，可能由不同深度和组分的储层组成。夏威夷地区的反演结果显示出两处明显的低阻异常，在观测点下方2—3 km处存在一个薄岩浆层，电阻率大约为20 Ω·m，该层位正对应了大岛裂谷系统中地震频发的深度，这些震群是由于岩浆囊在通道中的体积变化或压强变化引起应力集中的结果；在3—100 km深度范围内，视电阻率曲线随深度增加整体呈现先增大再逐渐降低的特征，大约在30—300 Ω·m之间变化，且在10—100 km深度范围内电性变化较为明显，推测为岩浆通道。该深度范围与Matoza等（2021）探测到的深度在20—50 km之间的横向连续地震活动扩散带相对应。而100 km深度左右电阻率变化范围在20—70 Ω·m，结合Constable和Heinson （2004）的反演认为地幔柱在40—400 km存在一个低电阻体，电阻率约为10 Ω·m。Witze （2013）认为在地下约110—155 km处存在一个暖池，该暖池并非位于夏威夷主岛的正下方，而是在主岛西侧100 km处，推测100 km以下可能存在大型的复杂岩浆储层，为上方岩浆囊提供深部热源（张炯等，2017）。
![]()
图 5 2018年5月至2019年9月测点下方的电阻率变化示意图Figure 5. Schematic diagram of resistivity change under the station from May 2018 to September 2019结合夏威夷岛的地震动态和每月的电性变化特征进行分析，可以看出在地震和火山活动频发的2018年5—7月份，电阻率呈上升趋势（图6）。地震的电阻率效应源于微裂隙内孔隙流体的变化，故推测地震余震期间电阻率的增加趋势是由震后应变释放区应力恢复和孔隙度恢复所致（王立凤等，2017）。
![]()
图 6 2018年4月至8月测点下方的电阻率变化示意图Figure 6. Schematic diagram of resistivity change under the station from April to August，20183. 结论
基于加拿大凤凰公司在夏威夷大岛上布设的长期大地电磁观测台站所采集的数据，对该地区开展了大地电磁反演研究，获取了其深部电性结构，单个台站得到的结果仅对于该测点下方是合理的，可能不能代表远离该测点下方的地下结构特征。结合该地区一年来发生的多次火山地震事件，发现观测点下方约2—3 km处存在明显的低阻岩浆囊，100 km以下可能是热源储层。研究表明，使用长期大地电磁观测台站可以对火山下方的电性结构特征进行长期的动态监测，大地电磁勘探对于认识和掌握火山区深部构造及变化特征具有重要的参考价值。本文仅得到了对夏威夷大岛的单点初步探测结果，要进一步了解大岛的深部电性结构及其动态变化特征，仍需进行更大范围和更长期的观测和研究。
本研究数据来自加拿大凤凰地球物理公司的公开数据，作者在此表示感谢。
-
![]()
图 1 汶川地震震源区东北端InSAR位移场及余震分布(修改自Hashimoto et al，2009 )
Figure 1. Deformation observed by InSAR and after-shocks (circles) in the northeastern section of Wenchuan earthquake source zone（modified from Hashimoto et al, 2009 ）
![]()
图 2 青川强余震及用于震源参数反演的远震台站分布（a）和波形拟合残差随深度分布（b）
Figure 2. Teleseismic stations and the Qingchuan strong aftershocks （a） as well as waveforms misfit versus depth （b）
![]()
图 3 青川强余震震源机制解及CAPtel反演后的远震体波观测波形（黑色）与理论波形（红色）对比图
台网及台站代号标记在波形左侧上方，台站下方分别为震中距/方位角，波形下方的两行数字分别为波形对齐所需的时移 （上行）和相关系数百分值 （下行）
Figure 3. Comparison between observed (black) and synthetic (red) seismograms after the CAPtel inversion of the Qingchuan strong aftershock.
Seismic network and stations are labeled to the left of the traces，while epicentral distance and azimuth are displayed below the station code，time shift （above） and percentile of cross correlation coefficient （below） are displayed under the waveforms
![]()
图 4 ABKT地震台的远震直达P波及其深度震相观测与理论波形对比
P波按照实际到时对齐，理论地震图所使用的深度标记于波形的下方
Figure 4. Observed and synthetic P，pP and sP for station ABKT
The waveforms are aligned on onset of P wave，and the focal depth used in synthetics are displayed below each trace
![]()
图 5 恩施地震台的三分量波形对比
Figure 5. Comparison among radial，tangential and vertical waveform recorded by station ENH
![]()
图 6 恩施台的振幅谱 （a），频散曲线测量图 （b） 以及时间域地震波形图 （c）
Figure 6. Spectral amplitude versus period （a），group velocity versus period （b） and seismograms （c） at station ENH
![]()
图 7 恩施台的瑞雷波频谱极小点周期与震源深度的关系
Figure 7. The dependence of period at spectral null and focal depth
![]()
图 8 不同质心深度时的InSAR位移场
断层位错上端点深度分别为 6.02 km （a）， 11.02 km （b）， 16.02 km （c）及21.02 km （d）
Figure 8. Prediction InSAR deformation with different depths
The depth of the upper edge of the ruptured fault is placed at 6.02 km （a），11.02 km （b）， 16.02 km （c） and 21.02 km （d），respectively
表 1 汶川强余震目录
Table 1 Catalog of strong aftershocks of the Wenchuan earthquake
序号 发震时刻 （UTC） 震中位置 震源
深度/kmMS MS7 mL mb mB MW 年-月-日 时：分：秒 东经/° 北纬/° 1 2008-05-12 06：43：14 103.82 31.27 14 6.3 5.9 5.7 6.2 6.3 2 2008-05-12 06：54：16 103.59 31.26 13 5.8 5.8 5.3 5.8 6.2 3 2008-05-12 11：11：01 103.67 31.26 14 6.3 6.2 5.8 5.8 6.2 6.1 4 2008-05-13 07：07：08 103.42 30.95 14 6.1 5.9 5.7 5.6 6.0 5 2008-05-17 17：08：24 105.08 32.20 13 6.1 5.9 5.9 5.5 6.1 6 2008-05-25 08：21：47 105.48 32.55 14 6.4 6.2 6.3 5.7 6.1 6.1 7 2008-07-24 07：09：27 105.61 32.76 10 6.0 5.6 5.7 5.7 5.7 8 2008-08-01 08：32：41 104.85 32.02 14 6.2 6.0 5.9 5.6 6.1 9 2008-08-05 09：49：15 105.61 32.72 13 6.5 6.3 6.0 5.8 6.3 6.0 注：表中序号为3，6，9的数据来自GCMT （2008），其余来自中国地震台网中心 （2008）。 
下载: 导出CSV
表 2 2008年5月25日青川MW6.1余震震源参数对比
Table 2 The source parameters of the May 25，2008 MW6.1 Qinchuan aftershock from various authors or agencies
数据来源 震中位置 深度/km 节面Ⅰ 节面Ⅱ 东经/° 北纬/° 走向/° 倾角/° 滑动角/° 走向/° 倾角/° 滑动角/° GCMT （2008） 105.45 32.57 27 59 84 178 149 88 6 郑勇等 （2009） 105.39 32.62 18 63 64 −171 329 82 −26 郭祥云等 （2010） 251 74 170 344 80 16 吕坚等 （2008） 105.37 32.62 19* 中国地震台网中心 （2008） 105.48 32.55 14* 注：深度为质心深度，标*为破裂起始深度。
下载: 导出CSV
表 3 理论地震图及地表形变使用的速度结构模型
Table 3 Crustal structure used in this study
底部深度/km vP/(km·s−1） vS/（km·s−1） 密度/（103 kg·m−3） 1 2.50 1.00 1.20 2 4.00 2.10 2.40 22 6.10 3.50 2.75 42 6.30 3.60 2.80 46 7.20 4.00 3.10 ∞ 8.00 4.47 3.35
下载: 导出CSV
-
郭祥云, 陈学忠, 李艳娥. 2010. 2008年5月12日四川汶川8.0级地震与部分余震的震源机制解[J]. 地震, 30(1): 50-60. Guo X Y, Chen X Z, Li Y E. 2010. Focal mechanism solutions for the 2008 MS8.0 Wenchuan earthquake and part of its aftershocks[J]. Earthquake, 30(1): 50-60 （in Chinese）.
黄媛, 吴建平, 张天中, 张东宁. 2008. 汶川8.0级大地震及其余震序列重定位研究[J]. 中国科学: D辑, 38(10): 1242-1249. Huang Y, Wu J P, Zhang T Z, Zhang D N. 2008. Relocation of the M8.0 Wenchuan earthquake and its aftershock sequence[J]. Science in China: Series D, 51(12): 1703-1711.
罗艳, 倪四道, 曾祥方, 郑勇, 陈棋福, 陈颙. 2010. 汶川地震余震区东北端一个余震序列的地震学研究[J]. 中国科学: 地球科学, 40(6): 677-687. Luo Y, Ni S D, Zeng X F, Zheng Y, Chen Q F, Chen Y. 2010. A shallow aftershock sequence in the north-eastern end of the Wenchuan earthquake aftershock zone[J]. Science China : earth Sciences, 53(11): 1655-1664.
吕坚, 苏金蓉, 靳玉科, 龙锋, 杨雅琼, 张致伟, 汤兰荣, 李超. 2008. 汶川8.0级地震序列重新定位及其发震构造初探[J]. 地震地质, 30(4): 917-925. Lü J, Su J R, Jin Y K, Long F, Yang Y Q, Zhang Z W, Tang L R, Li C. 2008. Discussion on relocation and seismo-tectonics of the MS8.0 Wenchuan earthquake sequences[J]. Seismology and Geology, 30(4): 917-925 （in Chinese）.
钱韵衣, 倪四道. 2016. 核幔边界反射震相ScS对远震体波反演震源参数精度影响[J]. 地球物理学报, 59(6): 2014-2027. Qian Y Y, Ni S D. 2016. The effects of the core-reflected wave ScS on source parameters in inversion with teleseismic body waves[J]. Chinese Journal of Geophysics, 59(6): 2014-2027 （in Chinese）.
秦刘冰, 陈伟文, 倪四道, 韩立波, 罗艳. 2014. 基于相对质心震中的地震破裂方向性测定方法研究: 以2008年云南盈江MS6.0地震为例[J]. 地球物理学报, 57(10): 3259-3269. Qin L B, Chen W W, Ni S D, Han L B, Luo Y. 2014. A method of resolving earthquake rupture directivity with relative centroid location and it’s application to the 2008 Yingjiang MS6.0 earthquake[J]. Chinese Journal of Geophysics, 57(10): 3259-3269 （in Chinese）.
孙建宝, 梁芳, 沈正康, 徐锡伟. 2008. 汶川MS8.0地震InSAR形变观测及初步分析[J]. 地震地质, 30(3): 789-795. Sun J B, Liang F, Shen Z K, Xu X W. 2008. InSAR deformation observation and preliminary analysis of the MS8 Wenchuan Earthquake[J]. Seismology and Geology, 30(3): 789-795 (in Chinese).
王家庆, 张国宏, 单新建, 张迎峰. 2016. 2008年西藏改则地震多视角InSAR成果的三维形变解算及初步分析[J]. 地震地质, 38（4）:978-986. Wang J Q, Zhang G H, Shan X J, Zhang Y F. 2016. Three-dimensional deformation of the 2008 Gaize earthquakes resolved from InSAR measurements by multiple view angles and its tectonic implications[J]. Seismology and Geology, 38(4):978-986 （in Chinese）.
王德才, 倪四道, 李俊. 2013. 地震烈度快速评估研究现状与分析[J]. 地球物理进展, 28(4): 1772-1784. Wang D C, Ni S D, Li J. 2013. Research status of rapid assessment on seismic intensity[J]. Progress in Geophysics, 28(4): 1772-1784 （in Chinese）.
王卫民, 赵连锋, 李娟, 姚振兴. 2008. 四川汶川8.0级地震震源过程[J]. 地球物理学报, 51(5): 1403-1410. Wang W M, Zhao L F, Li J, Yao Z X. 2008. Rupture process of the MS8.0 Wenchuan earthquake of Sichuan, China[J]. Chinese Journal of Geophysics, 51(5): 1403-1410 （in Chinese）.
韦生吉. 2009. 稀疏台网震源参数方法研究[D]. 合肥: 中国科学技术大学: 1–152. Wei S J. 2009. Constraining Source Parameters with Sparse Network[D]. Hefei: University of Science and Technology of China: 1–152 (in Chinese).
吴忠良, 臧绍先. 1991. 用体波合成地震图方法确定渤海、永善两大地震的震源参数[J]. 地震学报, 13(1): 1-8. Wu Z L, Zang S X. 1991. Source parameters of Bohai earthquake of July 18, 1969 and Yongshan earthquake of May 11, 1974 from synthetic seismogram of body waves[J]. Acta Seismologica Sinica, 13(1): 1-8 （in Chinese）.
张培震, 徐锡伟, 闻学泽, 冉勇康. 2008. 2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因[J]. 地球物理学报, 51(4): 1066-1073. Zhang P Z, Xu X W, Wen X Z, Ran Y K. 2008. Slip rates and recurrence intervals of the Longmen shan active fault zone, and tectonic implications for the mechanism of the May 12 Wenchuan earthquake, 2008, Sichuan, China[J]. Chinese Journal of Geophysics, 51(4): 1066-1073 （in Chinese）.
张勇, 许力生, 陈运泰, 冯万鹏, 杜海林. 2008. 2007年云南宁洱MS6.4地震震源过程[J]. 中国科学: D辑, 38(6): 683-692. Zhang Y, Xu L S, Chen Y T, Feng W P, Du H L. 2009. Source process of MS6.4 earthquake in Ning’er, Yunnan in 2007[J]. Science in China: Series D, 52(2): 180-188.
郑勇, 马宏生, 吕坚, 倪四道, 李迎春, 韦生吉. 2009. 汶川地震强余震(MS≥5.6)的震源机制解及其与发震构造的关系[J]. 中国科学: 地球科学, 39(4): 413-426. Zheng Y, Ma H S, Lü J, Ni S D, Li Y C, Wei S J. 2009. Source mechanism of strong aftershocks (MS≥5.6) of the 2008-05-12 Wenchuan earthquake and the implication for seismotectonics[J]. Science in China: Earth Sciences, 2009, 52(6): 739-753.
中国地震台网中心. 2008. 中国地震目录[EB/OL]. [2018–02–01]. http://data.earthquake.cn. China Earthquake Network Center. 2008. China earthquake catalogue[EB/OL]. [2018–02–01]. http://data.earthquake.cn (in Chinese).
Aki K, Richards P G. 2002. Quantitative Seismology[M]. 2nd ed. California: University Science Books: 34.
Chen W W, Ni S D, Kanamori H, Wei S J, Jia Z, Zhu L P. 2015. CAP joint, a computer software package for joint inversion of moderate earthquake source parameters with local and teleseismic waveforms[J]. Seismol Res Lett, 86(2A): 432-441.
Elliott J R, Walters R J, England P C, Jackson J A, Li Z, Parsons B. 2010. Extension on the Tibetan Plateau: Recent normal faulting measured by InSAR and body wave seismology[J]. Geophys J Int, 183(2): 503-535.
Frohlich C, Davis S D. 1999. How well constrained are well-constrained T, B, and P axes in moment tensor catalogs?[J]. J Geophys Res, 104(B3): 4901-4910.
GCMT. 2008. Monthly CMT solutions [EB/OL]. [2017–10–22]. http://www.ldeo.columbia.edu/~gcmt/projects/CMT/catalog/NEW_MONTHLY/2008/mar08.ndk.
Hao K X H, Si H J, Fujiwara T, Ozawa T. 2009. Coseismic surface-ruptures and crustal deformations of the 2008 Wenchuan earthquake MW7.9, China[J]. Geophys Res Lett, 36(11): L11303.
Hashimoto M, Enomoto M, Fukushima Y. 2010. Coseismic deformation from the 2008 Wenchuan, China, earthquake derived from ALOS/PALSAR images[J]. Tectonophysics, 491(1/2/3/4): 59-71.
He X, Ni S. 2017. Rapid rupture directivity determination of moderate dip-slip earthquakes with teleseismic body waves assuming reduced finite source approximation[J]. J Geophys Res, 122(7): 5344-5368.
He X H, Ni S D, Liu J. 2015. Rupture directivity of the August 3rd, 2014 Ludian earthquake (Yunnan, China)[J]. Science China: Earth Science, 58(5): 795-804.
Huang M H, Tung H, Fielding E J, Huang H H, Liang C R, Huang C, Hu J C. 2016. Multiple fault slip triggered above the 2016 MW 6.4 Meinong earthquake in Taiwan[J]. Geophys Res Lett, 43(14): 7459-7467.
IRIS. 2008. Earthquake event waveform data[EB/OL]. [2018–02–01]. http://ds.iris.edu/wilber3/find_event.
Jia Z, Ni S D, Chu R S, Zhan Z W. 2017. Joint inversion for earthquake depths using local waveforms and amplitude spectra of Rayleigh waves[J]. Pure Appl Geophys, 174(1): 261-277.
Kikuchi M, Kanamori H. 1991. Inversion of complex body waves: III[J]. Bull Seismol Soc Am, 81(6): 2335-2350.
Langston C A. 1980. A note on spectral nulls in Rayleigh waves[J]. Bull Seismol Soc Am, 70(4): 1409-1414.
Li Z H, Elliott J R, Feng W P, Jackson J A, Parsons B E, Walters R J. 2011. The 2010 MW 6.8 Yushu (Qinghai, China) earthquake: Constraints provided by InSAR and body wave seismology[J]. J Geophys Res, 116(B10): B10302.
Okada Y. 1985. Surface deformation due to shear and tensile faults in a half-space[J]. Bull Seismol Soc Am, 75(4): 1135-1154.
Qian Y Y, Ni S D, Wei S J, Almeida R, Zhang H. 2017. The effects of core-reflected waves on finite fault inversions with teleseismic body wave data[J]. Geophys J Int, 211(2): 936-951.
Shen Z K, Sun J B, Zhang P Z, Wan Y G, Wang M, Bürgmann R, Zeng Y H, Gan W J, Liao H, Wang Q L. 2009. Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake[J]. Nat Geosci, 2(10): 718-724.
Stein S, Wiens D A. 1986. Depth determination for shallow teleseismic earthquakes: Methods and results[J]. Rev Geophys, 24(4): 806-832.
Tsai Y B, Aki K. 1970. Precise focal depth determination from amplitude spectra of surface waves[J]. J Geophys Res, 75(29): 5729-5744.
Tong X P, Sandwell D T, Fialko Y. 2010. Coseismic slip model of the 2008 Wenchuan earthquake derived from joint inversion of interferometric synthetic aperture radar, GPS, and field data[J]. J Geophys Res, 115(B4): B04314. doi: 10.1029/2009JB006625.
Wang R J, Martıń F L, Roth F. 2003. Computation of deformation induced by earthquakes in a multi-layered elastic crust: FORTRAN programs EDGRN/EDCMP[J]. Comput Geosci, 29(2): 195-207.
Zhan Z W, Helmberger D, Simons M, Kanamori H, Wu W B, Cubas N, Duputel Z, Chu R S, Tsai V C, Avouac J P, Hudnut K W, Ni S D, Hetland E, Culaciati F H O. 2012. Anomalously steep dips of earthquakes in the 2011 Tohoku-Oki source region and possible explanations[J]. Earth Planet Sci Lett, 353/354: 121-133.
Zhao L S, Helmberger D V. 1994. Source estimation from broadband regional seismograms[J]. Bull Seismol Soc Am, 84(1): 91-104.
计量
- 文章访问数: 2312
- HTML全文浏览量: 1271
- PDF下载量: 130
