Application of receiver function method for characterizing the mantle transition zone structures in central and southern Alaska
-
摘要:
太平洋板块向北美板块俯冲,在北美洲阿拉斯加大陆边缘形成了碰撞造山带和俯冲构造带,前人对该区域阿拉斯加俯冲带地幔转换带的研究存在分歧,一种观点认为太平洋板块俯冲至阿拉斯加中南部下方的地幔转换带内部,引起地幔转换带增厚,另一种观点则认为地幔转换带增厚是由于库拉板块残留造成的。因此,研究阿拉斯加俯冲带区域地幔转换带的精确成像对于厘清太平洋板块俯冲机制和俯冲过程具有重要意义。本文利用美国地震联合学会(IRIS)在阿拉斯加区域布设的USArray部分固定台站接收到远震事件的波形资料,采用天然地震P波接收函数与共转换点道集叠加的方法,反演得到了阿拉斯加山脉和育空高原410 km、660 km间断面和地幔转换带的厚度结构。结果表明:由于太平洋板块俯冲至地幔转换带内部,使得410 km间断面在阿拉斯加山脉东部和育空高原与阿拉斯加山脉交界处的弧形区域出现了抬升现象,抬升幅度为0~20 km,进而导致在该区域出现不同程度的地幔转换带增厚的现象,验证了第一种观点。另外,该地区地幔转换带结构显示,660 km间断面在阿拉斯加山脉东部以及阿拉斯加山脉与育空高原交界的东部出现下沉现象,故此推测阿拉斯加中南部部分区域太平洋板块已经俯冲至地幔转换带底部。
Abstract:The subduction of the Pacific Plate towards the North American Plate has formed a collisional mountain belt and subduction zone at the edge of the continent in Alaska. Previous studies on the mantle transition zone beneath the Alaska subduction were in debate. It is believed that the thickening of the transition zone is due to the subduction of the Pacific Plate, or due to the thickening is caused by the residual Kula Plate. Therefore, accurate imaging of mantle transition zones in the Alaskan subduction zone is of great significance for understanding the subduction mechanism and courses of the Pacific Plate. This study provided waveform data from teleseismic events received in several stations of the US Array deployed by Incorporated Research Institutions for Seismology (IRIS) in the Alaska region, and adopts the natural earthquake P-wave receiver function and common conversion point stacking method to obtain 410 km and 660 km discontinuity sections and mantle transition zone thickness structure of the Alaska Range and Yukon Plateau. The new results show that due to the subduction of the Pacific Plate into the mantle transition zone, a 410 km long discontinuity is uplifted in the arc-shaped area of the eastern Alaska Mountains and the junction of the Yukon Plateau and the Alaska Mountains, with an uplift of 0~20 km, which in turn led to different degrees of thickening of the mantle transition zone in this area, which verifies the first view. In addition, the structure of the mantle transition zone in this area shows that a 660 km long discontinuity has subsided in the eastern part of the Alaska Mountain Range and the eastern part of the border between the Alaska Mountains and the Yukon Plateau. Therefore, it is assumed that the Pacific Plate in parts of south-central Alaska has subducted to the bottom of the mantle transition zone.
-
Keywords:
- receiver function /
- mantle transition zone /
- discontinuity /
- Alaska
-
为探究阿拉斯加俯冲带地幔转换带的构造特征,前人在该区域开展了广泛的研究,如全球、全区域尺度的体波走时层析成像研究[1-5]、接收函数研究[6-10]以及SS前驱波方法约束410 km和660 km界面起伏[11]等。
地震层析成像结果被用于讨论太平洋板块进入和穿过地幔转换带(Mantle Transtion Zone,MTZ,410~660 km)范围的问题。Qi等[12]在2007年的研究中发现400~450 km深度处板块运动速度为快速速度,同时在地幔转换带较低的位置也得到了板块运动的快速速度,因此他们认为地幔转换带较低的位置存在快速运动板块的残余。Van Stiphout等[13]应用层析成像技术研究得到的阿拉斯加俯冲带区域三维P波速度、b值(在地震学中,b值通常用于描述地震的震级分布特性)空间分布结果显示,深度为100 km的俯冲带下方活动火山呈现出高b值,且在深度为100 km的贝尼奥夫带顶部存在明显的高Vp/Vs异常,推测与板块脱水和产生岩浆有关;You等[14]使用层析成像技术对阿拉斯加南部进行了研究,发现该区域三维速度结构与P波各向异性在地壳及上地幔区域存在明显的速度不均匀性,活动火山下方沿着地幔楔的位置分布着明显的低速异常;Ferris等[7]和Ai等[15]采用接收函数方法对阿拉斯加地区进行了深入研究,但这些研究均未提供关于纵波、横波和泊松比深部结构的综合图像;Burdick等[4]在2017年构建的北美P波速度模型中得到了太平洋板块的末端似乎已经渗透到下地幔中,但在更深处模型的分辨率降低,不能对太平样板块的俯冲边界进行明确的计算;Maguire等[16]在2018年的研究模型表明,太平洋板块继续俯冲,俯冲至400 km以下,但对于该模型中同样存在一定深度(400 km)之下分辨率不高的问题,不能进行更精确的研究。Jiang等[17]在2018年构建的上地幔Vs模型中发现,俯冲板块在390 km深度处高速异常最为突出,且在地幔转换带中出现了较弱的高速异常。Gou等[18]的研究得到的P波速度模型显示了450~500 km深度处的太平洋板块,通过比较S波和P波速度情况,并通过“地下世界图谱”[5]确定了地幔转换带内的太平样板块末端。
不同层析模型的俯冲板块深度范围差异可能是由于地幔转换带模型分辨率不同造成的,因此需要使用不同方法来对阿拉斯加地区的地幔转换带和其可能范围及存在的板块进行成像。在成像研究中,接收函数法分析被用于研究全球地幔不连续性的地形,以便更深入地了解阿拉斯加下方的地幔转换带结构和俯冲板块的深度范围。
1. 区域地质背景
阿留申-阿拉斯加俯冲带东部(即阿拉斯加南部和阿拉斯加半岛)是由一系列堆积的地体构成[19],如雅库塔特地体、楚加奇地体、威廉王子地体等(图1)。这些地体起源于阿拉斯加南部北纬30°附近,并随着库拉板块持续地向极地方向移动[20],最终拼合到古老的陆核-北美克拉通之上。其中,弗兰盖利亚复合地体是一个大型地体,在100~55 Ma与阿拉斯加南部发生碰撞[21]。在55 Ma,楚加奇地体与弗兰盖利亚复合地体发生碰撞,导致楚加奇地体增生。在弗兰盖利亚复合地体增生后,阿拉斯加西南部地区经历了高达50°的逆时针旋转,将阿拉斯加和阿留申群岛的位置塑造成了我们今天所认识的阿拉斯加俯冲带的形状和位置。在白垩纪中期到古近纪早期,旋转发生之后,由于库拉板块的汇聚,沿着阿拉斯加南部形成了安第斯型的弧系。在24 Ma时,雅库塔特地体与太平洋板块耦合在一起,并随之俯冲至北美板块下方。
在板块俯冲的作用下,阿留申-阿拉斯加俯冲带成为北美构造运动最活跃的地区。该地区地震活动频发,但大部分为浅中源地震,震源深度超过300 km的深源地震不多。该俯冲带周围还孕育了比较多的火山,包括勘察加半岛6个火山、阿留申岛弧30个火山、阿拉斯加与库克湾19个火山。另外,在阿拉斯加地区存在较多的断层,如贯穿整个阿拉斯加州的德奈利断层、造成阿拉斯加中南部大部分地区发生火山和地震的阿留申巨型逆冲断层、将楚加奇地体和威廉王子地体分开的科恩台科特断层,它们都发育在古俯冲区、地质体内、裂谷沿线及两个板块相对运动摩擦的地区。
2. 数据与方法
2.1 数据来源
本研究主要使用USArray(https://earthquake.alaska.edu/)阵列中的AK台网的固定台站地震事件数据。其中AK地震台网的分布情况如图2中红色三角的排列所示,主要集中在阿拉斯加州育空高原和阿拉斯加山脉区域,在育空高原的排布较为工整,在阿拉斯加山脉南麓区域的台站分布更为密集;最终,本次研究共计使用了AK台网中103个地震台站数据(图2)。
前期通过在“国际地震中心ISC(International Seismological Centre(isc.ac.uk))”和“美国地质调查局(USGS)”的官网上搜索地震事件,选取了震中距为25°~100°、发震时刻为2015年1月1日至2023年1月31日之间且震级≥6.5级的地震事件。在此基础上,对地震目录中发震时间比较接近的地震进行筛选,确保2分钟内无重复性地震。经过筛选,本研究选取了201个地震事件,并以阿拉斯加研究区为中心绘制出震中与研究区域之间的相对位置关系(图3)。这些地震数目多,分布分散,足以保证台站附近的地震射线密度,从而进一步提高接收函数反演结果的可靠性。
图 3 本文涉及到的大于6.5级地震事件的分布及与研究区域的位置关系图图中黄色三角为研究区域,红色五角星为地震发生位置,小的红色圆圈为震中距25°的范围,大的红色圆圈为震中距100°的范围。Figure 3. The distribution of earthquakes with magnitude greater than 6.5 and the position of the study areaYellow triangle: the study area; red stars: earthquake epicenter; small red circles: the range of 25° from the epicenter, and large red circles: the range of 100° from the epicenter.2.2 方法原理
由于地球内部构造是各向异性的介质,因此地震波在地下传播的过程中会受到影响,而这种影响是多种因素共同作用的结果,通常对地震波传播受到的下方介质的影响特别关注,这样可以揭示地震台站下方的地质结构特征。通过分析这些影响,可以获得地震台站下方的地下结构速度间断面以及由此产生的转换波和多次反射波的信息。因此,接收函数通常被视为一种反演工具,可用于研究台站下方的地质结构[22-23]。
地震波形数据是由震源激发地震波经过一定的传播路径到达台站下方并被陆基台站或海底地震仪记录下来,包含震源响应S(t)、震源下方响应MS(t)、接收区(即台站下方的结构)响应MR(t)、传播路径响应Mr(t)、仪器响应I(t)的多重信息,因此地震波信号D(t)可以表示为上述5个响应的褶积形式[24]:
$$ D\left(t\right)=I\left(t\right)\mathrm{*}S\left(t\right)\mathrm{*}{M}_{\mathrm{s}}\left(t\right)\mathrm{*}{M}_{\mathrm{r}}\left(t\right)\mathrm{*}{M}_{\mathrm{R}}\left(t\right) $$ (1) 在上述5个响应中,只有仪器响应是已知的,而其余4种都是未知的,而且均与地下结构密切相关。为了将接收区介质的响应从其余响应中剥离出来,Langston提出等效震源假设[22],消除对远震P波波形的影响。S(t)、MS(t)、Mr(t)可以看作垂直分量和径向分量远震波形的共同因子。公式可在频率域中改写为:
$$ {D}_{\mathrm{z}}\left(\omega \right)=I\left(\omega \right)\mathrm{*}S\left(\omega \right)\mathrm{*}{M}_{\mathrm{s}}\left(\omega \right)\mathrm{*}{M}_{\mathrm{r}}\left(\omega \right)\mathrm{*}{M}_{\mathrm{R}}^{\mathrm{Z}}\left(\omega \right) $$ (2) $$ {D}_{\mathrm{R}}\left(\omega \right)=I\left(\omega \right)\mathrm{*}S\left(\omega \right)\mathrm{*}{M}_{\mathrm{s}}\left(\omega \right)\mathrm{*}{M}_{\mathrm{r}}\left(\omega \right)\mathrm{*}{M}_{\mathrm{R}}^{\mathrm{R}}\left(\omega \right) $$ (3) 根据公式(2)和(3)可得到去除接收区外其他所有响应,接收区介质在径向分量上的响应可以表示为:
$$ {M}_{\mathrm{R}}^{\mathrm{R}}\left(\omega \right)=\frac{{M}_{\mathrm{R}}^{\mathrm{Z}}\left(\omega \right)\mathrm{*}{D}_{\mathrm{R}}\left(\omega \right)}{{D}_{\mathrm{z}}\left(\omega \right)} $$ (4) 由于远震事件的入射角接近90°,垂直分量在接收区介质的响应MRZ(t)可以表示为一个尖脉冲δ(t),变换到频率域后值为1,因此公式变为:
$$ {M}_{\mathrm{R}}^{\mathrm{R}}\left(\omega \right)=\frac{{D}_{\mathrm{R}}\left(\omega \right)}{{D}_{\mathrm{z}}\left(\omega \right)} $$ (5) 再将MRR(t)经过傅里叶逆变换回到时间域,就得到了接收区介质响应的时间函数MRR(ω),被称为径向接收函数。
图4对径向接收函数的传播路径和震相进行了简单描述。莫霍面是地壳和地幔分界的速度间断面,当远震P波到达莫霍面时,大部分波形能量将继续作为P波(震相Pp)向上传播,而部分波形将转换为S波(震相Ps)穿透莫霍面向上传播。
3. 数据处理
3.1 数据预处理
编写Python脚本批量下载完成所需的AK网络的台站地震数据后,用Python软件的ObsPy工具包对地震波形数据进行了初步查看(图5)。其中图5A为AK台网中的D20K台站的数据,在各个分量都能够看到明显有效信号的同时,不同分量上的线性趋势也十分明显,进行数据接收函数计算时需要进行进一步的处理;图5B展示的台站数据是质量较好的波形信号,根据振幅的变化,可以很容易地判断出有效信号的到达时间。在对所下载的地震波形信号进行逐一检查后,发现AK台网的地震台站的原始信号也存在不同程度的趋势性。可以看出,预处理过程对获得高质量的地震波形数据和计算的接收函数结果至关重要。
通过图5可知,部分台站的数据存在信噪比较低的情况。为了使接收函数的计算以及波形反演出来的结果更加准确,进行计算前还需要对下载的数据进行预处理,以提升数据的信噪比和可靠性。本研究所涉及的预处理流程如下:
(1)地震仪器校正
由于不同台站的仪器响应可能不同,为使波形数据能相互比较,通过地震台站的仪器响应文件或者零极点文件对地震数据进行仪器响应的去除滤波,即对原始地震记录去平均与去线性趋势等处理和按照需求对波型做带通滤波处理,保留0.05~2.5 Hz频带地震记录,去除部分异常波动(图6)。
(2)旋转方位角
在地震事件的头段数据中添加地震事件的经纬度、震源深度以及地震台站的经纬度和高程信息,以便于计算反方位角,根据反方位角,将原始的ZNE三分量地震记录旋转到ZRT坐标;图7揭示了ZNE坐标系和ZRT坐标系两者之间的位置关系。
(3)截取地震波形
本文采用的方法是利用全球平均速度模型IASP91计算得到地震事件的预期走时,然后在地震记录上定位该事件并从中截取所需的波形数据。对于接收函数,所截取的地震记录为理论P波到时前10 s,以及P波到时后120 s地震数据,共计130 s(图8)。
3.2 接收函数计算及挑选
本文选用频率域迭代反卷积接收函数的方法对各个台站进行计算,分别获得了研究区域下载台站的径向接收函数和切向接收函数的结果。
为了获得台站高质量的接收函数数据,对提取的接收函数用Python的Seispy[25]工具包再次进行挑选,通过比对,挑选出P波信号突出、Ps震相在P波震相后3~5 s内有明显起跳、整体质量高的接收函数数据。在进行最终的接收函数汇总时,尽量保证接收函数为来自各个方位的地震事件计算而得,且每个台站都挑选了20条以上的有效地震记录,保障了结果的可靠性。在地壳和上地幔的边界处由于两边的物质组成元素相差过大,形成了强波阻抗界面,因此在速度结构上会有一个突然的增强,在莫霍面边界上会产生Ps、PpPs和PpPs+PpSs等震相,台站下方的莫霍面埋深和速度变化会影响到上述震相的到时,因此在选择数据时尽量选择这3个震相清晰的数据,用来提高后续共转换点叠加的准确性,获得可靠的410 km和660 km间断面形态分布特征。通过挑选,最终得到每个台站可使用的接收函数。图9给出了AK台网中BPAW台站和HAD台站下方按照后方位角排列的接收函数,图中可以看出莫霍界面的转换震相(Ps震相)清晰。
图 9 部分地震台站选取的接收函数显示A为BPAW台站下方按照后方位角排列的接收函数,B为HAD台站下方按照后方位角排列的接收函数。Figure 9. Display of selected receiver functions in some seismic stationsA represents the reception pattern beneath the BPAW station arranged according to the rear azimuth, and B represents the reception pattern beneath the HAD station arranged according to the rear azimuth.本文共选取了2个台站数据进行接收函数计算,2个台站的具体位置见图9,通过结果比对,可以很明显地观察到在径向接收函数的显示图中,每种震相相较于切向接收函数而言,在图中显示的更为明显,因此在本文的后续反演研究中,选择径向接收函数进行地壳厚度和440 km及660 km深度的叠加反演。同时,在径向分量的接收函数显示中,不仅直达P波的震相明显,而且能较为清晰地识别出3~5 s的Pms震相,这一特征为接下来的地壳深度的叠加反演提供了保障。
4. 结果
通过收集大量地震台站的地震事件信息来筛选接收函数,使用全球平均速度模型IASP91进行时深转换和共转换点道集的叠加,最终获得了阿拉斯加俯冲带区域阿拉斯加山脉及育空高原地区下方的410 km间断面和660 km间断面的形态特征,同时根据深度差获得过渡带厚度分布(图10)。
为了得到更为准确的共转换叠加结果,在进行叠加时,本文选择的圆形面元半径为150 km,并以100 km的间隔进行移动,最终得到阿拉斯加州中南部研究区域的叠加结果。通过三维的共转换点叠加的结果可以观察到410 km和660 km界面的区域起伏特征,其中410 km界面主要的埋深范围为380~450 km,660 km界面的主要埋深范围为630~680 km,剔除坏点之后,在410 km总共获得165个410 km界面的深度值,在660 km总共获得188个660 km界面的深度值。
图10A显示研究区域内410 km界面存在明显的纵向变化,从低纬度到高纬度,410 km的埋深整体呈现“深-浅-深”的状态分布,在64°N、156°~141°W的区域,410 km界面呈弧形上升,埋深从430 km上升到390 km左右,呈东西向展布的特征。410 km的界面主要在382~438 km之间变化,平均埋深为412.7 km。图10B展现出的660 km界面的分布特征大致与410 km界面的特征相似,但整体上升的区域较410 km界面向南北方向分布的更广,只有在65°N和60°~62°N、147°~141°W的区域出现下沉,分别从660 km下沉至680 km和670 km,整体上研究区域的660 km深度层埋深为634~688 km,660 km界面的平均埋深为663.6 km。
在进行共转换点道集叠加成像时,反演得到的数据会受到上地幔速度结构的影响,获取的叠加图显示410 km和660 km的埋深与实际埋深情况可能存在差异,因此使用地幔转换带的厚度来反映其自身结构更为合理。将反演得到的660 km界面的数据与410 km界面的数据点相减之后得到地幔过渡带的厚度,根据图10C可以直观地观测到阿拉斯加中南部区域地幔转换带厚度的变化情况。在本文的研究区中,地幔转换带的厚度在202~286 km之间变化,上述现象也显示出该区域地幔转换带结构在纵向上的不均一性,其中在64°N附近出现了弧状的厚度减小的现象,最薄的区域可达到202 km,小于全球地幔转换带的平均厚度,而在弧状区域的北侧,即65°~66°N的范围内,地幔转化带厚度出现明显增厚的现象,最厚处可达到286 km,大于全球地幔转换带的平均厚度。再次结合410 km和660 km界面的起伏情况进行研究,在弧状的地幔转换带变薄的区域,410 km界面并没有出现明显的抬升或下沉,界面埋深处于410~420 km;在660 km界面的显示中,可以看到弧状地幔转换带变浅区域有明显的向上抬升,界面约处于630~650 km。
通过三维共转换点道集叠加获取了研究区域410 km和660 km界面埋深的起伏情况,在两层界面上均发现了抬升和下沉,为了进一步了解和研究育空高原与阿拉斯加山脉北部连接处410 km和660 km界面的上突和地幔转换带厚度的减薄,本研究在形成三维共转换点成像的基础上,并考虑剖面上能够显示更多的台站数据,因此在研究区选择了5条测线进行了共转换点叠加计算并获取其剖面图(图11)。
测线A1沿151°W穿过阿拉斯加山脉和育空高原部分地区(60°~67°N)。在叠加剖面上,410 km间断面上震相清晰、可连续追踪,埋深上可以看到明显的起伏,大约在剖面图测线长度约320 km(62.8°N)的位置上,410 km间断面出现明显抬升,由原来的420 km左右的埋深抬升至390~400 km;同时,在间断面处于420 km的位置时,也能观察到剖面图在测线长度约200~300 km(61.8°~62.7°N)内410 km间断面的起伏,其余范围内的起伏情况较为平稳。整体来看,410 km间断面的埋深主要分布于390~420 km,埋深稍小于全球平均水平。与410 km间断面不同,660 km间断面震相清晰,埋深的横向起伏较为缓和,起伏位置主要发生在测线长度为300~400 km(62.7°~63.6°N),660 km间断面埋深主要在650~660 km,多数在660 km左右,与全球平均值接近。
测线A2沿145°W横穿育空高原部分地区(60°~68°N)。在整个叠加剖面上,410 km间断面上震相清晰、可连续追踪,埋深上在测线长度为0~550 km(60°~64.9°N)内十分平稳,埋深处在400 km左右的深度,在测线长度为550~750 km上410 km间断面发生起伏,且出现震相消失的情况,直到在测线长度760~880 km(66.9°~68°N)再次出现明显且连续的震相,且分别在390~440 km处形成两条震相清晰、可连续追踪的间断面,结合三维共转换点偏移叠加结果,由于66.9°~68°N区域处在共转换点叠加边缘,用于共转换点叠加的数量较少,结果存在不稳定性。与410 km间断面不同,660 km间断面震相清晰,埋深的横向起伏缓和,在测线长度700~880 km(66.5°~68°N)上出现下沉的现象,660 km间断面的整体埋深主要在660~680 km之间,稍大于全球660 km间断面的平均埋深水平。
测线B1沿65°N横穿育空高原部分区域(156°~1141°W)。由图11-B1可知,在整个叠加剖面上,410 km间断面上震相较为清晰、可连续追踪,在横向上可以清晰地看到410 km间断面的轻微起伏情况,在测线长度为260 km(150.5°W)处,410 km的间断面开始缓慢抬升,由原来的400 km左右抬升至380~390 km,整体来看,B1测线的410 km间断面埋深较浅且起伏缓和,主要分布于380~400 km。在后半段的410 km间断面的埋深稍小于全球410 km间断面的埋深水平。与410 km间断面不同,660 km间断面震相清晰,埋深在横向上出现明显的起伏,主要是在测线长度为380 km(147.9°W)的位置出现了间断面下沉的情况,从埋深660 km位置处下降约20 km到达680 km的埋深处,660 km间断面的整体埋深主要为650~680 km。
测线B2沿64°N横穿育空高原和阿拉斯加山脉交接地段(156°~141°W)。在叠加剖面上,410 km间断面上震相清晰、可连续追踪,埋深上在测线长度为200~720 km(141°~151.9°W)十分平稳,埋深处于约410 km的深度,在测线长度为0~200 km(151.9°~156°W)上410 km间断面呈现出向上抬升的现象,最浅的埋深可达到390 km,整体上,410 km间断面起伏缓和,大部分处于410 km左右的埋深。与410 km间断面不同,660 km间断面震相清晰,埋深在横向上几乎没有明显的起伏,整体埋深主要为650~660 km,符合660 km间断面埋深的全球平均值。
测线B3沿62°N,横穿阿拉斯加山脉部分区域(156°~141°W)。在叠加剖面上,410 km间断面震相清晰、可连续追踪,无明显激烈的横向变化,埋深由西向东逐渐恢复到全球正常埋深,从390 km逐渐恢复至410 km,并且在测线长度250~600 km(144.5°~151.2°W)间断面埋深为410 km深度,在测线长度为600~780 km(141°~144.5°W)时,间断面缓慢下沉至420 km深度处。660 km间断面的分布形式与410 km间断面的情况相似,但更为平缓。660 km间断面走势平缓,横轴上前600 km(144.5°~156°W)的范围内基本稳定在660 km的深度,之后逐渐下沉至680 km的深度。660 km间断面埋深稍大于全球平均水平。
当比正常地幔温度更低的地体(例如大洋板块)俯冲到地幔转换带时,且410 km的间断面具有正斜率,考虑到克拉布龙斜率的特征影响,深度为410 km的速度间断面会变浅;而具有负斜率的660 km速度间断面则会变得更深[26]。随着一系列变化的发生,410 km和660 km的速度间断面将会产生不同的反应:前者上升,后者下沉,地幔转换带的厚度也将增加(图12)。因为太平洋板块属于大洋板块,当太平洋板块俯冲至地幔转换带内部时,穿过410 km间断面,温度低的太平洋板块使410 km的速度间断面抬升,埋深变浅,进而使地幔转换带增厚。本文的研究结果表明部分区域存在660 km速度间断面下降、埋深变深的情况,根据克拉布龙斜率的特征影响,当低温板块与660 km速度间断面接触时,会使间断面下降,因此认为存在太平洋板块已俯冲至地幔转换带底部的情况。
5. 讨论
Van Stiphout等[13]利用接收函数共转换点叠加结果所显示的阿拉斯加中南部410 km间断面出现局部抬升至390~395 km深度处的现象,指出阿拉斯加中南部已经存在部分太平洋板块俯冲至地幔转换带内部的观点,不过并没有俯冲至转换带的底部。一些学者的研究也在一定程度上证实了这一结论,Gou等[18]在2019年采用三维P波速度与各向异性成像研究,指出阿拉斯加半岛下方存在一个高速异常体,该异常体从上地幔的顶部一直持续到深度大约450~500 km的位置。文章所提到高速异常体的走向基本上是东西方向,与阿留申俯冲带东段平行。该高速异常体一直向东延伸,到达阿拉斯加中南部地区下方。这一发现也验证北太平洋板块已俯冲到阿拉斯加半岛中南部下方转换带的观点。可以理解为,该高速异常体是与阿留申俯冲带东段平行的一条地质构造线,其延伸方向与阿拉斯加中南部地区下方的地质构造有关。
肖勇等[11]在2020年利用SS前驱波方法研究得到的结果认为阿拉斯加中南部下方地幔转换带增厚是由库拉板块残留引起,其研究在叠加过程中采用了2°的叠加半径,对于小尺度地质构造(板块残留)特征不利于分辨,同时他们认为接收函数和体波走时层析成像中探测的均是小尺度地质特征,如Dahm等[9]在2017年的研究显示,通过接收函数和体波走时层析成像,发现阿拉斯加中部地区下方410 km间断面上升幅度很小。Burdick等[4]在2017年所观测到的高速异常体仅在东西方向有大尺度的展布,而在南北方向的展布尺度大约只有100 km。因此他们认为在目前的研究中,阿拉斯加俯冲带的北太平洋板块并未俯冲到转换带深度范围内,且在该地区尚未出现大规模的地质构造变化[11]。
本研究在南北方向和东西方向分别设置了2条和3条测线,共转换点偏移叠加成像的结果显示,无论在南北方向还是在东西方向,都出现了地幔转换带层厚变厚的现象,且规模较大。Van Stiphout等[13]的研究显示仅有410 km间断面发生了抬升,抬升至390~395 km的深度范围,而在三维共转换点偏移叠加的图10中发现,不仅部分区域的410 km间断面有抬升情况,一些区域的660 km间断面同样存在下沉的现象,一般下沉至约680 km的深度范围,在A2、B1、B3的测线上存在200~400 km的区域显示660 km间断面的下沉。
基于上述分析,认为阿拉斯加俯冲带下方的太平洋板块已经俯冲至地幔转换带内部,考虑到研究区内的部分区域存在660 km间断面整体平缓下沉至大约680~690 km深度的现象,推测阿拉斯加中南部部分区域太平洋板块可能俯冲至地幔转换带底部,但是需要进一步的研究和佐证。
6. 结论
对阿拉斯加中南部区域内的103个地震台站所接收到的高质量地震事件信息进行接收函数计算,并对计算的接收函数进行二次挑选得到震相清晰的接收函数。进行三维共转换点偏移叠加时,使用圆形面元叠加,并设计面元半径为150 km,面元以55 km的步长移动,最终获得了阿拉斯加州中南部下方三维410 km和660 km间断面的起伏和分布形态;在测线上进行共转换点偏移叠加时,根据测线上台站分布情况,选择适应于测线的单次叠加距离和移动步长,得到了测线下方410 km和660 km间断面的起伏情况。研究表明:
(1)阿拉斯加山脉东部(60°~63°N、142°~149°W)以及阿拉斯加与育空高原交界的弧形区域(64°~66°N、141°~156°W)地幔转换带增厚,厚度为260~270 km。在阿拉斯加中南部区域,410 km间断面在育空高原与阿拉斯加山脉交界处的弧形区域(64°~65.8°N、141°~156°W)出现了抬升现象,抬升幅度为0~20 km;同时660 km间断面在阿拉斯加山脉与育空高原交界的东部(64.8°~66°N、141°~148°W)和阿拉斯加山脉东部部分区域(60°~63°N、142°~146°W)出现下沉现象,下沉幅度0~20 km,其他区域660 km间断面的埋深处于全球平均水平。
(2)阿拉斯加俯冲带下方的太平洋板块已经俯冲至地幔转换带内,考虑到研究区内的部分区域存在660 km间断面整体平缓下沉至670~680 km的现象,认为在阿拉斯加山脉东部(60°~63.5°N、142°~146°W)以及阿拉斯加山脉与育空高原交界的东部(64.8°~66°N、141°~148°W)太平洋板块可能俯冲至地幔转换带底部,但需要进一步的研究和佐证。
-
图 3 本文涉及到的大于6.5级地震事件的分布及与研究区域的位置关系图
图中黄色三角为研究区域,红色五角星为地震发生位置,小的红色圆圈为震中距25°的范围,大的红色圆圈为震中距100°的范围。
Figure 3. The distribution of earthquakes with magnitude greater than 6.5 and the position of the study area
Yellow triangle: the study area; red stars: earthquake epicenter; small red circles: the range of 25° from the epicenter, and large red circles: the range of 100° from the epicenter.
图 9 部分地震台站选取的接收函数显示
A为BPAW台站下方按照后方位角排列的接收函数,B为HAD台站下方按照后方位角排列的接收函数。
Figure 9. Display of selected receiver functions in some seismic stations
A represents the reception pattern beneath the BPAW station arranged according to the rear azimuth, and B represents the reception pattern beneath the HAD station arranged according to the rear azimuth.
-
[1] 柳存喜, 王志. 南阿拉斯加地壳及上地幔结构成像研究[J]. 地球物理学报, 2014, 57(7 doi: 10.6038/cjg20140708 LIU Cunxi, WANG Zhi. Structure imaging of the crust and upper mantle in south of Alaska [J]. Chinese Journal of Geophysics, 2014, 57(7). doi: 10.6038/cjg20140708
[2] Zhao D, Yamamoto Y, Yanada T. Global mantle heterogeneity and its influence on teleseismic regional tomography. Gondwana Research, 2013, 23(2): 595-616
[3] Martin-Short R, Allen R M, Bastow I D. Subduction geometry beneath south central Alaska and its relationship to volcanism. Geophysical Research Letters, 2016, 43(18): 9509-9517
[4] Burdick S, Vernon F L, Martynov V, et al. Model Update May 2016: Upper-Mantle Heterogeneity beneath North America from Travel-Time Tomography with Global and USArray Data. Seismological Research Letters, 2017, 88(2A): 319- 325
[5] van der Meer D G, van Hinsbergen D J J, Spakman W. Atlas of the underworld: slab remnants in the mantle, their sinking history, and a new outlook on lower mantle viscosity [J]. Tectonophysics, 2018, 723: 309-448. doi: 10.1016/j.tecto.2017.10.004
[6] Ai Y S, Zheng T Y, Xu W W, et al. A complex 660 km discontinuity beneath northeast China [J]. Earth and Planetary Science Letters, 2003, 212(1-2): 63-71. doi: 10.1016/S0012-821X(03)00266-8
[7] Ferris A, Abers G A, Christensen D H, et al. High resolution image of the subducted Pacific (?) plate beneath central Alaska, 50-150 km depth [J]. Earth and Planetary Science Letters, 2003, 214(3-4): 575-588. doi: 10.1016/S0012-821X(03)00403-5
[8] Rondenay S, Abers G A, van Keken P E. Seismic imaging of subduction zone metamorphism [J]. Geology, 2008, 36(4): 275-278. doi: 10.1130/G24112A.1
[9] Dahm H H, Gao S S, Kong F S, et al. Topography of the mantle transition zone discontinuities beneath Alaska and its geodynamic implications: constraints from receiver function stacking. Journal of Geophysical Research: Solid Earth, 2017, 122(12): 10352-10363
[10] van Stiphout A M, Cottaar S, Deuss A. Receiver function mapping of mantle transition zone discontinuities beneath Alaska using scaled 3-D velocity corrections [J]. Geophysical Journal International, 2019, 219(2): 1432-1446. doi: 10.1093/gji/ggz360
[11] 肖勇. 阿留申-阿拉斯加俯冲带和汤加-斐济俯冲带地幔过渡带结构研究[D]. 北京: 中国地震局地球物理研究所硕士学位论文, 2020. XIAO Yong. Study of the mantle transition zone structure beneath Alaska-Aleutian subduction zone and Tonga-Fiji subduction zone[D]. Master Dissertation of Institute of Geophysics, China Earthquake Administration, 2020.
[12] Qi C, Zhao D P, Chen Y. Search for deep slab segments under Alaska [J]. Physics of the Earth and Planetary Interiors, 2007, 165(1-2): 68-82. doi: 10.1016/j.pepi.2007.08.004
[13] van Stiphout T, Kissling E, Wiemer S, et al. Magmatic processes in the Alaska subduction zone by combined 3-D b value imaging and targeted seismic tomography [J]. Journal of Geophysical Research:Solid Earth, 2009, 114(B11): B11302.
[14] You T, Zhao D P. Seismic anisotropy and heterogeneity in the Alaska subduction zone [J]. Geophysical Journal International, 2012, 190(1): 629-649. doi: 10.1111/j.1365-246X.2012.05512.x
[15] Ai Y S, Zhao D P, Gao X, et al. The crust and upper mantle discontinuity structure beneath Alaska inferred from receiver functions [J]. Physics of the Earth and Planetary Interiors, 2005, 150(4): 339-350. doi: 10.1016/j.pepi.2004.12.002
[16] Maguire R, Ritsema J, Bonnin M, et al. Evaluating the Resolution of Deep Mantle Plumes in Teleseismic Traveltime Tomography. Journal of Geophysical Research: Solid Earth, 2018, 123(1): 384-400
[17] Jiang C X, Schmandt B, Ward K M, et al. Upper mantle seismic structure of Alaska from Rayleigh and S wave tomography [J]. Geophysical Research Letters, 2018, 45(19): 10350-10359. doi: 10.1029/2018GL079406
[18] Gou T, Zhao D P, Huang Z C, et al. Aseismic deep slab and mantle flow beneath Alaska: insight from anisotropic tomography [J]. Journal of Geophysical Research:Solid Earth, 2019, 124(2): 1700-1724. doi: 10.1029/2018JB016639
[19] Rossi G, Abers G A, Rondenay S, et al. Unusual mantle Poisson's ratio, subduction, and crustal structure in central Alaska [J]. Journal of Geophysical Research:Solid Earth, 2006, 111(B9): B09311.
[20] Haeussler P, Coe R S, Onstott T C. Paleomagnetism of the late Triassic hound island volcanics: revisited [J]. Journal of Geophysical Research:Solid Earth, 1992, 97(B13): 19617-19639. doi: 10.1029/92JB01361
[21] McClelland W C, Gehrels G E, Saleeby J B. Upper Jurassic-lower cretaceous basinal strata along the cordilleran margin: implications for the accretionary history of the Alexander-Wrangellia-Peninsular terrane [J]. Tectonics, 1992, 11(4): 823-835. doi: 10.1029/92TC00241
[22] Langston C A. Corvallis, Oregon, crustal and upper mantle receiver structure from teleseismic P and S waves [J]. Bulletin of the Seismological Society of America, 1977, 67(3): 713-724. doi: 10.1785/BSSA0670030713
[23] Vinnik L P. Detection of waves converted from P to SV in the mantle [J]. Physics of the Earth and Planetary Interiors, 1977, 15(1): 39-45. doi: 10.1016/0031-9201(77)90008-5
[24] Phinney R A. Structure of the Earth's crust from spectral behavior of long-period body waves. Journal of Geophysical Research, 1964, 69(14): 2997-3017.
[25] Xu M J, He J. Seispy: Python module for batch calculation and postprocessing of receiver functions [J]. Seismological Research Letters, 2023, 94(2A): 935-943. doi: 10.1785/0220220288
[26] 朱智昊. 利用接收函数方法研究中国中东部410和660千米间断面起伏[D]. 中国科学技术大学硕士学位论文, 2018 ZHU Zhihao. Topography of the 410-and 660-km discontinuities beneath the central and eastern China using receiver function[D]. Master Dissertation of University of Science and Technology of China, 2018.