秦满忠 张元生 刘旭宙 李顺成

1) 中国兰州730000中国地震局兰州地震研究所 2) 中国兰州730000中国地震局地震预测研究所兰州科技创新基地 3) 中国北京100029中国地震局地质研究所地震动力学国家重点实验室



利用甘肃地震台网记录的sP前驱波研究汤加—斐济俯冲带“410”间断面深度

1) 中国兰州730000中国地震局兰州地震研究所 2) 中国兰州730000中国地震局地震预测研究所兰州科技创新基地 3) 中国北京100029中国地震局地质研究所地震动力学国家重点实验室

本文首先利用甘肃数字测震台网和甘东南野外观测流动台阵记录的汤加—斐济俯冲带深震震相sP， 使用倾斜叠加方法有效提取sP在近源区“410”速度间断面底界面反射的弱前驱波s410P； 其次基于CRUST1.0修正的IASP91模型， 通过震相sP与其前驱波s410P的到时差对比， 给出了汤加—斐济俯冲带“410”速度间断面底部反射点的深度约为398.5 km， 抬升约11.5 km， 这与俯冲带与近源区“410”间断面的相互作用有关； 最后得到了该区域新的地震波速度模型IASP91_Tonga. 本文结果将对更好地认识该俯冲带的复杂结构、 正确理解地球深部的动力学过程具有重要的参考价值．

sP震相 前驱波s410P 间断面 俯冲带 甘肃数字测震台网 甘东南野外观测流动台阵

引言

地震波(P， S， pP， sP等)在传播过程中， 遇到间断面会发生反射或转换， 所产生的新震相(P410P， P660P， S410S， S660S， S410P， S660P， s410P等)为间断面(“410”， “660”等)的深入研究提供了地震学基础(Flanagan， Shearer， 1998a； Deuss， 2009； Schmerr， Thomas， 2011)． 这些前驱震相在单条记录中的能量(幅度)往往很弱， 常淹没于噪声之中； 但是通过叠加大量的观测波形数据， 则可以压制噪声， 更加有效地提取与间断面相关的有用信号(Rost， Thomas， 2002； 臧绍先， 周元泽， 2002)．

很多研究者利用下行的转换震相在间断面的转换点深度确定该间断面的深度． 例如： Richards和Wicks(1990)使用SdP转换波, 利用非线性叠加方法对汤加地区下方的“670”间断面的深度和性质进行了深入研究； 谢彩霞等(2012)利用4次根倾斜叠加方法有效提取了离源下行的SdP次生转换震相， 进一步确认了汤加—斐济下方300 km附近速度间断面的存在； Li等(2008)通过叠加S--P转换波， 计算了中国东北地区“660”间断面的深度， 并探讨了俯冲带对“660”间断面的影响．

发生在俯冲带中的深震会在近源区速度间断面产生反射震相(p410P， s410P等)和转换震相(S410P， S660P等)(Flanagan， Shearer， 1998b)． 长周期震相(PP， SS)前驱波方法对狭窄俯冲带间断面的横向分辨具有很大局限性(Flanagan， Shearer， 1998a； Schmerr， Thomas， 2011)． 接收函数方法是研究地球内部间断面的一种有效方法， 可以通过P波与其上行转换波Pds的到时差来确定间断面的深度(Langston， 1979)， 但其主要适用于三分量地震台站下方壳、 幔间断面的研究．

在深源地震(震源深度h>300 km)分布密集的俯冲带(如汤加—斐济俯冲带)， 使用深震震相sP在近源区速度间断面底界面反射的震相s410P确定狭窄俯冲带间断面的横向变化具有很大的优势． 本文将通过倾斜叠加大量的观测波形数据， 有效提取sP在近源区“410”速度间断面底界面反射的弱前驱波s410P， 利用sP震相与其前驱波s410P震相的到时差计算汤加—斐济俯冲带“410”速度间断面的深度．

1 台网资料

甘肃地震台网由甘肃数字测震台网和甘东南野外观测流动台阵构成， 二者均处于青藏高原东北缘． 甘肃“十五”数字测震台网自2008年6月正式运行， 由兰州、 高台、 安西、 嘉峪关和天水等5个有人值守的国家数字地震台和39个区域遥测数字地震台组成(冯建刚等， 2012)， 目前拥有包括邻省(宁夏、 陕西、 四川、 青海、 内蒙)在内的76个数字测震台站． 甘东南野外观测流动台阵是中国地震局地质研究所和中国地震局兰州地震研究所共同合作在甘东南地区勘址架设， 共设有7条测线， 150个流动台站， 平均台间距约为10 km(图1)． 所有台站均使用REF TEK-130型数据采集器和Guralp公司生产的CMG-3ESPC宽频带地震计． 甘东南野外观测流动台阵运行期间(2009年11月—2011年12月)积累了连续、 可靠、 高质量的地震观测波形数据， 为深入研究地球内部物理结构提供了宝贵资料(秦满忠等， 2015)．

图1 甘肃数字测震台网和甘东南野外观测流动台阵分布

本文选取了甘肃数字测震台网和甘东南野外观测流动台阵共同记录的2011年9月15日发生在汤加—斐济俯冲带(21.61°S， 179.53°W，h=644 km，mB=7.3)的深远地震数字观测波形数据， 定位结果参考了USGS(2011)地震目录．

2 数据处理

深震震相sP与其前驱波s410P在地幔中的传播路径基本相同． sP是上行的S波在震中附近地表反射转换的P波； s410P是上行的S波在震中附近“410”间断面底界面反射转换的P波(图2)， 该震相在初至震相P与sP之间到达(图3)．

图2 射线(sP， s410P和P)传播示意图Fig.2 Ray paths of sP， s410P and P

首先对所选取的汤加—斐济深远地震观测波形数据作重采样、 去倾斜、 去均值处理;然后作带通滤波处理(0.1—0.01 Hz)， 同时去掉记录畸形、 信噪比差的观测波形； 最终选取了甘肃数字测震台网和甘东南野外观测流动台阵共同记录到的170条垂直分量观测波形数据． 如图3所示， 汤加—斐济深远震在台网记录中的震中距为88°—102°．

图3 甘肃地震台网记录到的2011年9月15日发生在汤加—斐济俯冲带的深远地震的数字观测波形sMohoP， sLABP和s410P震相分别表示sP在近源区莫霍面， 岩石圈和“410”的前驱波Fig.3 Waveforms of mB=7.3 event occurred in Tonga-Fiji on September 15， 2011， recorded by the Gansu Seismic Network The seismic phases sMohoP， sLABP and s410P are the precursors from near-source underside reflection off Moho， lithosphere-asthenosphere boundary，and “410”, respectively

通过倾斜叠加大量的观测波形数据， 可以压制噪声， 有效地提取与间断面相关的弱信号(Ritsemaetal， 1995； 臧绍先， 周元泽， 2002)． 在进行倾斜叠加的过程中， 选取震中距为94°的台站(红崖山台)作为参考台， sP为参考震相， 并以sP震相到时为零时刻(图3)． 信号i(s410P)在任意台s的到时为Δtir+ΔpΔDrs， 以红崖山台为参考台的倾斜叠加结果为

(1)

式中， Δtir为参考震中距记录的信号i(s410P)与参考震相sP之间的走时差， Δp为信号i在任意台s的水平慢度与信号i(震中距为参考台震中距)水平慢度之差， ΔDrs为任意台s与参考台的震中距之差，m为使用地震记录的台站数．

以sP为参考震相， 使用倾斜叠加方法获得了在震中附近“410”间断面底界面反射的前驱波信号s410P， 其反射点位置为(21.25°S， 179.99°W)． 该信号在倾斜叠加波形(图3中红色线条)中成像清晰， 与参考震相sP的观测走时差为137.43 s， 使用TauP软件结合理论IASP91全球速度模型得到的s410P与震相sP理论走时差为137.04 s； 通过使用基于CRUST1.0模型(Laskeetal， 2013)校正后的汤加—斐济地区速度模型IASP91_Tonga(图4)， 反算出“410”间断面底部反射点的深度约为398.5 km， “410”间断面抬升约11.5 km， 这与俯冲带与近源区“410”间断面的相互作用有关．

3 讨论与结论

甘肃数字测震台网和甘东南野外观测流动台阵记录的汤加—斐济深震震相pP和sP均具有尖锐、 清晰和较大振幅的记录特征， sP记录振幅甚至超过了直达P波， 这为我们使用前驱波(sMohoP， sLABP和s410P等)对近源区间断面(莫霍面， 岩石圈， “410”等)的深入研究提供了基础资料．

图4 全球速度模型IASP91以及修改后的汤加—斐济地区速度模型IASP91_Tonga Fig.4 IASP91 velocity model and modified velocity model (IASP91_Tonga) of the Tonga-Fiji zone in this paper

鉴于前驱波震相在单条记录中的能量(振幅)往往很弱， 本文通过倾斜叠加大量的观测波形数据， 获得了在近源区“410”间断面底界面反射的前驱波震相s410P， 并使用CRUST1.0模型校正后的汤加—斐济地区速度模型IASP91_Tonga， 反算出“410”间断面底部反射点的深度约为398.5 km．

图3中震相sLABP和sMohoP相对于参考震相sP的走时差分别为23.16 s和10.62 s. 我们同样对震源区的地壳及岩石圈结构作了修改(图4)， 使用TauP软件计算得到近源区岩石圈的深度约为46 km， 莫霍面深度约为10 km．

俯冲带对“410”和“660”的影响是目前研究的一个热点(Lietal， 2008)． 该研究涉及到间断面的性质， 又能反映地幔对流的性质和形式． 俯冲带是冷的下沉物质， 其通过间断面时会引起间断面的形态变化． 若“410”为相变界面， 俯冲带会使它抬升(Vidale， Benz， 1992； Collier， Helffrich， 1997); 若“660”为相变界面， 俯冲板块中的冷物质将使它下沉． 因此， 汤加—斐济俯冲带中冷的物质与“410”作用时， 俯冲带及相邻地幔中的橄榄石-尖晶石的相变界面会上升(蒋志勇等， 2003)． 本文中得到的汤加—斐济俯冲带“410”间断面抬升约11.5 km与该结论一致．

本文使用倾斜叠加技术成功提取了sP前驱波(s410P， sLABP和 sMohoP)， 获得了俯冲带间断面的深度， 给出了汤加—斐济俯冲带新的地震波速度模型(IASP91_Tonga)． 本文结果将为深入研究该区域的地球结构提供基础资料， 对正确认识俯冲带的复杂结构以及深入理解地球深部的动力学过程具有重要的参考价值．

冯建刚， 蒋长胜， 韩立波， 陈继锋. 2012. 甘肃测震台网监测能力及地震目录完整性分析[J]. 地震学报， 34(5): 646--658.

Feng J G， Jiang C S， Han L B， Chen J F. 2012. Analysis on the monitoring capability of seismic networks and completeness of earthquake catalogues in Gansu region[J].ActaSeismologicaSinica， 34(5): 646--658 (in Chinese).

蒋志勇， 臧绍先， 周元泽. 2003. 鄂霍次克海下间断面的起伏及俯冲带的穿透[J]. 科学通报， 48(4): 320--327.

Jiang Z Y， Zang S X， Zhou Y Z. 2003. Topographies of seismic velocity discontinuities and penetrations of subducting slabs beneath the Sea of Okhotsk[J].ChineseScienceBulletin， 48(4)： 320--327 (in Chinese).

秦满忠， 李顺成， 刘旭宙， 张淑珍， 张元生. 2015. 核幔界面衍射波Pdiff震相最大记录震中距探讨[J]. 地震学报， 37(2): 249--256.

Qin M Z， Li S C， Liu X Z， Zhang S Z， Zhang Y S. 2015. A discussion on the largest recorded epicentral distance for seismic phase Pdiffdiffracted on core-mantle boundary[J].ActaSeismologicaSinica， 37(2): 249--256 (in Chinese).

谢彩霞， 周元泽， 王卓君， 欧东新. 2012. 汤加—斐济地区300 km间断面的SdP转换波证据[J]. 地球物理学报， 55(5): 1591--1600.

Xie C X， Zhou Y Z， Wang Z J， Ou D X. 2012. Evidence of SdP conversion phases for the 300 km discontinuity beneath Tonga-Fiji region[J].ChineseJournalofGeophysics， 55(5): 1591--1600 (in Chinese).

臧绍先， 周元泽. 2002.N次根倾斜叠加方法在间断面研究中的应用[J]. 地球物理学报， 45(3): 407--415.

Zang S X， Zhou Y Z. 2002.N-th root slant stack and its application in study of mantle discontinuities[J].ChineseJournalofGeophysics， 45(3): 407--415 (in Chinese).

Collier J D， Helffrich G R. 1997. Topography of the “410” and “660” km seismic discontinuities in the Izu-Bonin subduction zone[J].GeophysResLett， 24(12): 1535--1538.

Deuss A. 2009. Global observations of mantle discontinuities using SS and PP precursors[J].SurvGeophys， 30(4): 301--326.

Flanagan M P， Shearer P M. 1998a. Global mapping of topography on transition zone velocity discontinuities by stacking SS precursors[J].JGeophysRes， 103(B2): 2673--2692.

Flanagan M P， Shearer P M. 1998b. Topography on the 410-km seismic velocity discontinuity near subduction zones from stacking of sS， sP， and pP precursors[J].JGeophysRes， 103(B9): 21165--21182.

Langston C A. 1979. Structure under Mount Rainier， Washington， inferred from teleseismic body waves[J].JGeophysRes， 84(B4)： 4749--4762.

Laske G， Ma Z T， Masters G， Pasyanos M. 2013. A new global crustal model at 1×1 degrees[EB/OL]. [2015-04-26]. http:∥igppweb.ucsd.edu/～gabi/crust1.html#download.

Li J， Chen Q F， Vanacore E， Niu F L. 2008. Topography of the 660-km discontinuity beneath northeast China: Implications for a retrograde motion of the subducting Pacific slab[J].GeophysResLett， 35(1): L01302.

Richards M A， Wicks C W. 1990. S--P conversion from the transition zone beneath Tonga and the nature of the 670 km discontinuity[J].GeophysJInt， 101(1): 1--35.

Ritsema J， Hagerty M， Lay T. 1995. Comparison of broadband and short-period seismic waveform stacks: Implications for upper-mantle discontinuity structure[J].GeophysResLett， 22(23): 3151--3154.

Rost S， Thomas C. 2002. Array seismology: Methods and applications[J].RevGeophys， 40(3): 2-1--2-27. doi:10.1029/2000RG000100.

Schmerr N， Thomas C. 2011. Subducted lithosphere beneath the Kuriles from migration of PP precursors[J].EarthPlanetSciLett， 311(1/2): 101--111.

USGS. 2011. Search earthquake archives[EB/OL]. [2015-04-26]. http:∥earthquake.usgs.gov/earthquakes/search/.

Vidale J E， Benz H M. 1992. Upper-mantle seismic discontinuities and the thermal structure of subduction zones[J].Nature， 356(6371): 678--683.

On “410” discontinuity depth of the Tonga-Fiji subduction zone by the sP precursor recorded by Gansu seismic network

1)LanzhouInstituteofSeismology，ChinaEarthquakeAdministration，Lanzhou730000，China2)LanzhouBaseofInstituteofEarthquakeScience，ChinaEarthquakeAdministration，Lanzhou730000，China3)StateKeyLaboratoryofEarthquakeDynamics，InstituteofGeology，ChinaEarthquakeAdministration，Beijing100029，China

The sP seismic phase of the Tonga-Fiji subduction zone have a clear image in the Gansu Digital Seismic Network and southeastern Gansu temporary observation array. In this paper， we extracted the weak precursors s410P by slant stack method， and the s410P precursors are produced by deep focus earthquakes that result from near-source underside reflection off “410”. We compute the depth of “410” discontinuity by using the differential travel time between sP and its precursor s410P near Tonga-Fiji subduction zone based on the modified velocity model (IASP91_Tonga) in this zone based on CRUST1.0 model. The depth of “410” discontinuity is about 398.5 km， and uplifted about 11.5 km， which maybe relate to the interaction of the subduction zone and the “410” discontinuity in the near-source zone. This result is of great significance to recognize the complex structure of subduction zones and objectively understand the dynamic process of the deep Earth.

seismic phase sP; precursor s410P; discontinuity; subduction zone; Gansu Digital Seismic Network; southeastern Gansu temporary observation array

中国地震局地震预测研究所基本科研业务专项(2012IESLZ09)和国家自然科学基金(41504050)联合资助.

2015-06-06收到初稿， 2015-08-25决定采用修改稿.

e-mail: manzhongqin@eyou.com

10.11939/jass.2016.01.005

P315.63

A

秦满忠, 张元生, 刘旭宙, 李顺成. 2016. 利用甘肃地震台网记录的sP前驱波研究汤加—斐济俯冲带“410”间断面深度. 地震学报, 38(1): 53--58. doi:10.11939/jass.2016.01.005.

Qin M Z， Zhang Y S， Liu X Z， Li S C. 2016. On “410” discontinuity depth of the Tonga-Fiji subduction zone by the sP precursor recorded by Gansu seismic network.ActaSeismologicaSinica, 38(1): 53--58. doi:10.11939/jass.2016.01.005.