On “410” discontinuity depth of the Tonga-Fiji subduction zone by the sP precursor recorded by Gansu seismic network
-
摘要: 本文首先利用甘肃数字测震台网和甘东南野外观测流动台阵记录的汤加—斐济俯冲带深震震相sP, 使用倾斜叠加方法有效提取sP在近源区“410”速度间断面底界面反射的弱前驱波s410P; 其次基于CRUST1.0修正的IASP91模型, 通过震相sP与其前驱波s410P的到时差对比, 给出了汤加—斐济俯冲带“410”速度间断面底部反射点的深度约为398.5 km, 抬升约11.5 km, 这与俯冲带与近源区“410”间断面的相互作用有关; 最后得到了该区域新的地震波速度模型IASP91_Tonga. 本文结果将对更好地认识该俯冲带的复杂结构、 正确理解地球深部的动力学过程具有重要的参考价值.Abstract: 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.
-
引 言
地震波(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).
本文选取了甘肃数字测震台网和甘东南野外观测流动台阵共同记录的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).
图 3 甘肃地震台网记录到的2011年9月 15日发生在汤加—斐济俯冲带的深远 地震的数字观测波形sMohoP, sLABP和s410P震相分别表示sP在 近源区莫霍面, 岩石圈和“410”的前驱波 The seismic phases sMohoP, sLABP and s410P are the precursors from near-source underside reflection off Moho, lithosphere-asthenosphere boundary, and “410”, respectivelyFigure 3. Waveforms of mB=7.3 event occurred in Tonga-Fiji on September 15, 2011, recorded by the Gansu Seismic Network首先对所选取的汤加—斐济深远地震观测波形数据作重采样、 去倾斜、 去均值处理; 然后作带通滤波处理(0.1—0.01 Hz),同时去掉记录畸形、 信噪比差的观测波形; 最终选取了甘肃数字测震台网和甘东南野外观测流动台阵共同记录到的170条垂直分量观测波形数据. 如图 3所示,汤加—斐济深远震在台网记录中的震中距为88°—102°.
通过倾斜叠加大量的观测波形数据,可以压制噪声,有效地提取与间断面相关的弱信号(Ritsema et al,1995; 臧绍先,周元泽,2002). 在进行倾斜叠加的过程中,选取震中距为94°的台站(红崖山台)作为参考台,sP为参考震相,并以sP震相到时为零时刻(图 3). 信号i(s410P)在任意台s的到时为Δtir+ΔpΔDrs,以红崖山台为参考台的倾斜叠加结果为
式中,Δ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模型(Laske et al,2013)校正后的汤加—斐济地区速度模型IASP91_Tonga(图 4),反算出“410”间断面底部反射点的深度约为398.5 km,“410”间断面抬升约11.5 km,这与俯冲带与近源区“410”间断面的相互作用有关.
3. 讨论与结论
甘肃数字测震台网和甘东南野外观测流动台阵记录的汤加—斐济深震震相pP和sP均具有尖锐、 清晰和较大振幅的记录特征,sP记录振幅甚至超过了直达P波,这为我们使用前驱波(sMohoP,sLABP和s410P等)对近源区间断面(莫霍面,岩石圈,“410”等)的深入研究提供了基础资料.
鉴于前驱波震相在单条记录中的能量(振幅)往往很弱,本文通过倾斜叠加大量的观测波形数据,获得了在近源区“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”的影响是目前研究的一个热点(Li et al,2008). 该研究涉及到间断面的性质,又能反映地幔对流的性质和形式. 俯冲带是冷的下沉物质,其通过间断面时会引起间断面的形态变化. 若“410”为相变界面,俯冲带会使它抬升(Vidale,Benz,1992; Collier,Helffrich,1997); 若“660”为相变界面,俯冲板块中的冷物质将使它下沉. 因此,汤加—斐济俯冲带中冷的物质与“410”作用时,俯冲带及相邻地幔中的橄榄石-尖晶石的相变界面会上升(蒋志勇等,2003). 本文中得到的汤加—斐济俯冲带“410”间断面抬升约11.5 km与该结论一致.
本文使用倾斜叠加技术成功提取了sP前驱波(s410P,sLABP和sMohoP),获得了俯冲带间断面的深度,给出了汤加—斐济俯冲带新的地震波速度模型(IASP91_Tonga). 本文结果将为深入研究该区域的地球结构提供基础资料,对正确认识俯冲带的复杂结构以及深入理解地球深部的动力学过程具有重要的参考价值.
-
图 3 甘肃地震台网记录到的2011年9月 15日发生在汤加—斐济俯冲带的深远 地震的数字观测波形
sMohoP, sLABP和s410P震相分别表示sP在 近源区莫霍面, 岩石圈和“410”的前驱波 The seismic phases sMohoP, sLABP and s410P are the precursors from near-source underside reflection off Moho, lithosphere-asthenosphere boundary, and “410”, respectively
Figure 3. Waveforms of mB=7.3 event occurred in Tonga-Fiji on September 15, 2011, recorded by the Gansu Seismic Network
-
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.
doi: 10.1029/2000RG000100 USGS. 2011. Search earthquake archives[EB/OL]. [2015-04-26]. http://earthquake.usgs.gov/earthquakes/search/.
-
期刊类型引用(3)
1. 靳炎锡,秦满忠,李少华. 利用前驱波研究地幔间断面的若干进展. 地震工程学报. 2023(01): 209-219+243 . 百度学术
2. 崔辉辉,高雅健,周元泽. 伊豆-小笠原地区岩石圈软流圈边界地震学证据. 科学通报. 2017(07): 711-720 . 百度学术
3. 秦满忠,刘旭宙,邹锐,张元生,郭晓,魏从信,王亚红,孙点峰. 甘肃祁连山大容量气枪主动源最大探测范围. 地震工程学报. 2017(06): 1070-1075 . 百度学术
其他类型引用(3)