中美洲俯冲带南段上覆板块沉积物化学组分及其地质意义*

赵仁杰 鄢全树, 3 张海桃 关义立 石学法

1.山东科技大学地球科学与工程学院，青岛 266590 2.自然资源部第一海洋研究所，海洋地质与成矿作用重点实验室，青岛 266061 3.青岛海洋科学与技术试点国家实验室，海洋地质过程与环境功能实验室，青岛 266061

全球汇聚板块边缘根据上覆板块剥蚀物质(前缘剥蚀及底部剥蚀)和增生物质(前缘增生及底侵增生)之间的平衡关系，可以划分为剥蚀边缘和增生边缘(von Huene and Scholl, 1991; Cliftetal., 2009a; Straubetal., 2020; 赵仁杰等, 2020)。在剥蚀边缘，上覆板块弧前区域可能经历剥蚀和结构坍塌形成前缘剥蚀，同时其底部可能经历由水压致裂和磨蚀造成的底部剥蚀，这些过程将地壳物质带入俯冲带(Pichonetal., 1993; von Hueneetal., 2004; Cliftetal., 2009b; Scholl and von Huene, 2009; Saffer and Tobin, 2011)。通常认为俯冲剥蚀过程只将大陆地壳物质带入地幔中，然而随着对上覆板块弧前结构及组分研究的不断深入，弧前剥蚀的物质可能还包含增生洋壳以及增生沉积物，这些剥蚀再循环物质对俯冲带岩浆作用以及地幔不均一性作出重要的贡献(Plank and Langmuir, 1998; Workmanetal., 2004; Goss and Kay, 2006; Willbold and Stracke, 2006; Straubetal., 2015, 2020)。因此，确定俯冲剥蚀物质组分是研究剥蚀型汇聚板块边缘物质再循环的先决条件。

自早中新世以来，中美洲俯冲带已经从增生边缘演化成典型的剥蚀边缘(Vannucchietal., 2003)。目前中美洲海沟走向在哥斯达黎加南部发生向陆的弯曲，这可能与科科斯脊的俯冲有关(Vannucchietal., 2013; Lietal., 2018)。前人对哥斯达黎加南部弧前陆架层序重建的研究表明，科科斯脊俯冲造成弧前大量的俯冲剥蚀，并在初始的0.3Ma期间移除了体积达1.2×106km3的物质，形成典型的“弧前沉积(弧前由于底部被俯冲剥蚀移除，陆源沉积物直接在弧前沉积)”，接受了大量的陆源沉积(Vannucchietal., 2013, 2016a)。尽管从构造地质学和地球物理角度的研究指出该区域上覆板块基底存在明显的剥蚀现象，但是由于钻探取样技术有限，人们仍未获得中美洲俯冲带南部弧前的基底岩石样品。前人在俯冲再循环的研究中使用中美洲俯冲带南部(中美洲俯冲带南部指哥斯达黎加中部到巴拿马之间的区域，分界线见图1a中黄色实线)弧前露头中的蛇绿岩套组分代表剥蚀端元组分(Goss and Kay, 2006)，但该观点还存在争议。此外，在中美洲俯冲带的“俯冲物质浅部再循环产物”(大陆弧火山岩)中可能保存了俯冲剥蚀物质参与俯冲物质再循环的直接证据(Straubetal., 2020)。

图1 地质背景及样品描述

本文对综合大洋钻探计划344航次两个站位(U1380和U1413)所获得的沉积物样品进行研究。由于两个站位均未获得基底岩石，本文对钻孔中沉积物的粗碎屑层位开展了系统的主、微量元素及Sr-Nd-Pb-Hf同位素地球化学分析，将其与中美洲俯冲带南部不同年龄的大陆弧火山岩、中美洲俯冲带南部弧前基底(Goss and Kay(2006)选用的俯冲剥蚀端元)、加勒比大火成岩省基底组分及古加拉帕戈斯热点轨迹组分进行对比，确定中美洲俯冲带南部上覆板块基底组分。此外，本文还探讨了俯冲剥蚀物质对中美洲俯冲带南部大陆弧火山岩成因的可能影响。

1 地质背景和样品描述

在东太平洋海隆东侧，两个次一级的构造板块(科科斯板块和纳兹卡板块)正沿着中美洲海沟向加勒比板块俯冲，是研究俯冲剥蚀边缘的典型区域(图1a)(Carretal., 2004; 李永祥等, 2013; 鄢全树和石学法, 2014)。沿中美洲海沟走向，输入板块的汇聚速率、火山弧地壳厚度、俯冲角度、年龄、组分等均存在显著的差异(Carretal., 1990; Lyleetal., 1995; Barckhausenetal., 2001; DeMets, 2001)。从危地马拉到哥斯达黎加，科科斯板块和加勒比板块之间的汇聚速率从60mm/yr增加到90mm/yr(DeMets, 2001)。弧下地壳厚度从危地马拉(48km)向尼加拉瓜(32km)逐渐减小，然后向南又逐渐增厚，在哥斯达黎加其厚度达到32～40km(Carretal., 1990)。俯冲角度从尼加拉瓜的60°减小到哥斯达黎加中部的40°(Lyleetal., 1995)，而在哥斯达黎加南部接近平板俯冲。从危地马拉到哥斯达黎加北部(中美洲俯冲带北部，分界线见图1a中黄色实线)，俯冲洋壳年龄大约为24.0Ma，是东太平海隆扩张形成的科科斯板块，其组分类似于大洋玄武岩(Barckhausenetal., 2001; Walker and Gazel, 2014)。然而，自哥斯达黎加中部到巴拿马(中美洲俯冲带南部，分界线见图1a中黄色实线)，俯冲板块被科科斯脊及西北侧海山群覆盖，其组分为洋岛玄武岩，年龄大约为13.0～14.5Ma(Carretal., 1990; Dziermaetal., 2011; Walker and Gazel, 2014)。科科斯脊及西北侧海山群是加拉帕戈斯热点与科科斯-纳兹卡扩张中心相互作用形成，科科斯脊长约1000km，宽200～300km，比周围海底高2km(Werneretal., 2003; Harppetal., 2005)。俯冲洋壳地形的差异造成中美洲海沟俯冲剥蚀速率的不同，中美洲俯冲带北部自早中新世正常洋壳俯冲以来，剥蚀速率为11.3～13.4km3/Myr/km；而中美洲俯冲带南部，自科科斯脊开始俯冲，短期的俯冲剥蚀速率可达1125km3/Myr/km，平均为108～123km3/Myr/km(Vannucchietal., 2013, 2016b)。同时，在中美洲火山弧南部，大陆弧火山岩在晚中新世发生显著变化，钙碱性的火山岩在～6Ma停止，同时出现碱性玄武岩和埃达克质岩，这些火山岩均具有加拉帕戈斯热点的同位素和微量元素特征，可能是由于俯冲样式变化造成的(Abratis and Wörner, 2001; Gazeletal., 2009, 2011; Morelletal., 2012; Morell, 2015)。最近根据火山岩的同位素地球化学数据和地震各向异性的综合分析，发现从哥斯达黎加到尼加拉瓜的弧下地幔楔中存在平行于海沟的物质流(Hoernleetal., 2008)。

为了探索大地震的成因机理与破坏过程，综合大洋钻探计划334和344航次在中美洲俯冲带南部哥斯达黎加西部奥萨半岛附近的海域实施了钻探，站位分布在俯冲板块和上覆板块(Harrisetal., 2013a)。其中站位U1380(8°35.99′N、84°4.39′W，水深502.7m)和U1413(8°44.46′N、84°6.80′W，水深540.0m)分别位于上覆板块陆坡的中部和上部。两站位的科学目标之一均为获得上覆板块基底岩石的性质、组分及物理特征，然而在实际实施过程中均未获得基底岩石(图1b)，两个站位详细的岩性单元见图1c(Harrisetal., 2013b, c)。本文选取U1380站位中5个层位：1个位于岩性单元ⅡA(深度为556.01m，岩性为含贝壳碎屑的粗砂)、2个位于岩性单元ⅡB(深度分别为701.61m和764.38m，岩性分别为粘土质粉砂和粉砂岩)、2个位于岩性单元Ⅲ中(深度为789.82m和796.98m，岩性均为粉砂质粘土岩)，对以上层位进行主、微量元素和Sr-Nd-Pb-Hf同位素测试；选取U1413站位中5个层位：2个位于岩性单元Ⅱ(深度为162.21m和170.78m，岩性均为粉砂质粘土岩)、3个位于岩性单元Ⅲ(深度分别为446.78m、513.56m和578.35m，岩性分别为粉砂岩，砂岩以及砂质粉砂岩)，其中上部的四个层位分别进行主、微量元素和Sr-Nd-Pb-Hf同位素测试，最底部的层位只进行Sr-Nd-Pb-Hf同位素测试。

2 分析方法

2.1 主、微量元素分析

实验将冻干后的样品研磨成200目，在105℃的烘箱中烘干3小时，取0.05g的样品加入特氟龙熔样内胆中，然后加入1.50mL高纯HNO3和1.50mL高纯HF，加盖及钢套密闭，放入195℃的烘箱48小时，将冷却后的熔样内胆取出放在电热板上蒸干。然后再加入1mL HNO3并再次蒸干。最后加入3mL高纯HNO3和0.5mL Rh内标溶液，置于150℃的烘箱中24小时，冷却后，将提取液转移至聚酯瓶中，用Mill-Q水稀释待测。主量元素除了SiO2外，在自然资源部海洋地质与成矿作用重点实验室利用电感耦合等离子光学发射光谱法(ICP-OES)测定；SiO2是在山东第四地质矿产勘查院利用X射线荧光光谱法测试。微量元素在自然资源部海洋地质与成矿作用重点实验室利用电感耦合等离子质谱仪(ICP-MS)完成。对于含量>1.0%的主量元素，精度为±0.2%～2%，微量元素的精度<10%。同时测定了烧失量。标样为BHVO-2，测定值与推荐值在误差范围内是一致的(见表1)。

2.2 Sr-Nd-Pb-Hf同位素分析

称取50mg研磨至200目的样品于聚四氟乙烯(PTFE)熔样瓶中，加入2mL HF、1.5mL HNO3和0.2mL HClO4，拧紧瓶盖，将熔样瓶置于120℃电热板上加热约一周至其完全溶解。熔样完全后，开盖蒸干，再升温至180℃，去除残余HClO4。用于Sr、Nd同位素测试的样品蒸干后，用2.5mL/L的HCl将样品溶解，然后转移到离心管中，离心后吸取上清液，采用AG50W-X12树脂离子交换柱进行Sr和REE的化学分离，接收下来的REE再用P5O7萃淋树脂离子交换柱分离Sm和Nd；用于Pb同位素测试的样品蒸干后，用0.6mL/L的HBr溶解样品，离心后吸取上清液，采用AG1-X8树脂离子交换柱进行Pb同位素的化学分离。沉积物Hf同位素采用碱熔法溶样，称取0.5g研磨至200目的样品粉末与1g的Li2B4O7混合均匀，在铂金坩埚中用Rigaku高频全自动熔样机在1250℃熔融并冷却制成玻璃片，然后用高压压片机将玻璃片压碎，称取0.3～0.4g上述样品于样品瓶中，加入6～8mL 2.5M HCl，将样品并置于60℃的电热板上加热直至完全溶解。将溶解好的样品离心，以备上柱分离Hf。分离Hf和基体元素以及干扰元素采用HCl-单柱Ln-Spec提取色谱方法。

Sr、Nd、Pb同位素测试采用高精度多接收电感耦合等离子体质谱仪(MC-ICP MS)完成，该测试工作在自然资源部海洋地质与成矿作用重点实验室完成。在测试过程中Sr同位素测试的标样为NBS987，测定值为87Sr/86Sr=0.710268±7(2σ)，Nd同位素测试的标样为Shin Etsu JNdi-1，测定值为143Nd/144Nd=0.512120±3(2σ)，Pb同位素测试的标样为NBS981，测定值为206Pb/204Pb=16.938，207Pb/204Pb=15.493，208Pb/204Pb=36.725(见表2)。Hf同位素采用多接收电感耦合等离子体质谱仪(MC-ICP-MS)完成测试，该测试工作在中国科学院广州地球化学研究所同位素地球化学国家重点实验室完成，标样为BHVO-2，其测定值为176Hf/177Hf=0.282877±10(2σ)(见表2)。

表2 站位U1380和U1413各层位Sr-Nd-Pb-Hf同位素比值

3 分析结果

3.1 主、微量元素

站位U1380和U1413沉积物主量元素中含量最高的为SiO2，其平均含量为53.47%；其次从高到低依次为Al2O3、Fe2O3T、MgO、CaO、Na2O、K2O、TiO2、P2O5以及MnO；同时，两站位的烧失量均较高，平均为9.31%，数据详见表1。

表1 站位U1380和U1413各层位主量(wt%)、微量(×10-6)元素分析结果

根据站位U1380和U1413沉积物中各层位原始地幔标准化后的微量元素蛛网图(图2)可以看出，两站位中各层位微量元素变化范围相对较小，均呈现出Rb、Ba、U、K元素的富集，以及Nb、Ta、Th的亏损。相比于站位U1380，站位U1413更富集不相容元素，更亏损相容元素。

图2 站位U1380和U1413各层位原始地幔标准化微量元素蛛网图(标准化值据Sun and McDonough, 1989)

根据站位U1380和U1413沉积物中各层位球粒陨石标准化后的稀土元素配分图(图3)可以看出，两站位均呈现出轻微的轻稀土富集，其配分模式较为平坦。相比于站位U1380，站位U1413相对富集轻稀土，亏损重稀土，但是亏损和富集程度相对较低。

图3 站位U1380和U1413各层位球粒陨石标准化稀土元素配分图(标准化值据Sun and McDonough, 1989)

3.2 Sr-Nd-Pb-Hf同位素

各层位的同位素数据见表2。沉积物样品中，87Sr/86Sr(0.704679～0.707899)，206Pb/204Pb(18.858～19.255)，207Pb/204Pb(15.597～15.930)，208Pb/204Pb(38.588～39.408)的变化范围较大。但143Nd/144Nd(0.512883～0.512973)与176Hf/177Hf(0.283120～0.283170)的变化范围较小。沉积物Sr-Nd-Pb-Hf同位素相关性图解(图4)显示：沉积物组分与CLIP基底的Nd(143Nd/144Nd：0.512823～0.513104)、Pb(206Pb/204Pb：18.714～19.903；207Pb/204Pb：15.551～15.623；208Pb/204Pb：38.526～39.673)、Hf(176Hf/177Hf：0.283101～0.283153)同位素组分(Geldmacheretal., 2003; Hauffetal., 2000)基本一致。

图4 站位U1380和U1413各层位同位素图解

4 讨论

4.1 中美洲俯冲带南部上覆板块基底组分

研究剥蚀型汇聚板块边缘物质再循环的首要条件之一是确定上覆板块基底组分。中美洲俯冲带上覆板块为加勒比板块，该板块位于南美洲和北美洲之间，大部分被加勒比大火成岩省覆盖，其起源目前仍存在一定争议(Hoernleetal., 2002; James, 2005; Dürkefäldenetal., 2019)。最初的“原位模型”研究认为加勒比大火成岩省(CLIP)在南北美洲之间形成(Meschede and Frisch, 1998; James, 2005)。然而随着在小安的列斯岛弧(加勒比板块西侧)中发现了侏罗纪时期太平洋种属的放射虫，并且其年龄比南北美洲板块分离的时间更早，因此越来越多的学者支持加勒比板块起源于太平洋(Hoernleetal., 2002; Dürkefäldenetal., 2019)。该模型认为，加勒比大火成岩省是加拉帕戈斯地幔柱的头部在法拉隆板块上形成，然后与加勒比巨型弧碰撞，由于其具有浮力阻塞俯冲带，使俯冲极性发生由东向西的反转，随后插入南北美洲之间(Hoernleetal., 2002; Dürkefäldenetal., 2019)。

目前认为加勒比板块沿太平洋侧(中美洲俯冲带)岩浆杂岩体的起源较为复杂，可能由增生的太平洋海山组分、残余的古俯冲带物质、上升的加勒比大火成岩省基底组分以及增生的古加拉帕戈斯热点物质组成(Hauffetal., 2000; Hoernleetal., 2002)。因此，以上组分都可能作为上覆板块被俯冲剥蚀的物质。例如，Goss and Kay(2006)研究中美洲俯冲带南部物质再循环时，选择哥斯达黎加弧前布里卡(Burica)、戈尔菲托(Golfito)、奥萨(Osa)和克波斯(Quepos)处的蛇绿岩代表俯冲剥蚀物质。前人研究指出戈尔菲托和布里卡为加勒比火成岩省基底上升暴露的部分，而奥萨和克波斯则是由加拉帕戈斯热点形成的无震脊和海山(Hauffetal., 2000)。Goss and Kay(2006)的研究清晰表明中美洲哥斯达黎加处晚中新世以来埃达克质岩的形成受到俯冲剥蚀物质的影响。然而，随后的研究表明戈尔菲托可能并非加勒比大火成岩省基底组分，而是75～66Ma的岛弧组分(Buchsetal., 2010)。因此，Goss and Kay(2006)选择的端元是否能够真实代表弧前基底组分，前人并没有给出充分的证据。

本研究将弧前钻孔沉积物中的粗碎屑层位组分与弧前各种构造单元组分对比，以期获得真实的弧前基底组分。根据钻孔位置，研究推测钻孔中物质来源可能为中美洲火山弧组分、中美洲俯冲带南部弧前基底组分(Goss and Kay(2006)研究中选择的俯冲剥蚀端元)，古加拉帕戈斯热点组分以及加勒比大火成岩省基底组分。本文主要探讨各俯冲端元对中美洲俯冲带南部大陆弧岩浆作用的影响，所有测试结果为全岩样品数据，然而在沉积过程中，陆源物质和海洋自生组分共同沉积，因此，我们重点关注不受沉积过程影响的碱金属元素(Rb、K)以及高场强元素(Nb、Ta、Ti、Zr和Hf)(Plank and Langmuir, 1998; Plank, 2014)。从元素蛛网图(图2)中可以看出，两个站位的元素含量较为均一，表明其物质源区是相同的；沉积物粗碎屑层位中碱金属元素和高场强元素落在中美洲俯冲带南部>6Ma大陆弧火山岩、中美洲俯冲带南部<6Ma大陆弧火山岩以及加勒比大火成岩省基底组分的区域(图2a, b, d)；相比于中美洲俯冲带南部弧前基底和古加拉帕戈斯热点轨迹，沉积中粗碎屑层位具有较高的碱金属含量(图2c, e)。研究表明位于大陆边缘快速沉积形成的海洋沉积物稀土元素主要受陆源物质影响(Mclennanetal., 1990; Plank and Langmuir, 1998)。中美洲俯冲带南部弧前由于俯冲剥蚀形成了典型的“弧前沉积”，接受了大量的陆源沉积(Vannucchietal., 2013)，因此本研究中沉积物粗碎屑层位组分的稀土配分模式主要受陆源物质影响。从图3中我们排除中美洲俯冲带南部<6Ma大陆弧火山岩来源的可能性。我们进一步从Sr-Nd-Pb-Hf同位素角度进行对比，从图4中可以看出，沉积物Sr同位素变化范围较大，这是由于沉积物中碳酸盐对Sr同位素的影响(孟宪伟等, 2000)，因此在本文中Sr同位素并不用于源区判断；从图4a中，我们发现，中美洲俯冲带南部>6Ma大陆弧火山岩的Pb同位素明显低于沉积物粗碎屑层位的Pb同位素特征，结合Nd-Hf同位素特征(在Nd-Hf同位素图解中，沉积物未落在任何区域，这可能是由于前人缺乏对低Nd含量样品开展Hf同位素测试造成的)，我们认为沉积物中的粗碎屑可能来自出露在中美洲俯冲带南部的加勒比大火成岩省基底组分。综上，尽管上覆板块的钻孔未达到基底，然而上覆板块中沉积物中粗碎屑层位可代表上覆板块基底组分。

4.2 中美洲俯冲带南部物质循环模型

4.2.1 中美洲俯冲带南部俯冲再循环组分的可能端元

俯冲剥蚀机制中，俯冲再循环物质一般包括俯冲洋壳，下伏岩石圈地幔、俯冲沉积物以及俯冲剥蚀物质(Scholl and von Huene, 2009; Straubetal., 2020)，这些再循环物质可能影响着俯冲带岛弧岩浆成因。因此，俯冲带岛弧火山岩成因的研究中，首先应明确可能的俯冲输入端元组分。前人研究表明，中美洲俯冲带南部大陆弧火山岩岩性在晚中新世时发生显著变化，岩性由钙碱性火山岩转变碱性玄武岩和埃达克质岩，其组分具有类似加拉帕戈斯热点的同位素和微量元素的特征，目前存在利用不同俯冲输入端元来解释该区域大陆弧火山岩岩性的变化(Abratis and Wörner, 2001; Goss and Kay, 2006; Gazeletal., 2009, 2011)。Abratis and Wörner(2001)认为厚的有浮力的无震脊(科科斯脊)俯冲在中美洲俯冲带南部形成构造板片窗，同时俯冲科科斯脊的前缘发生部分熔融对地幔楔产生影响。随后，Gazeletal.(2011)进一步提出板片拆沉模型，该模型认为科科斯脊俯冲造成先前俯冲的板片拆沉，致使受加拉帕戈斯热点影响的软流圈物质进入地幔楔，因此科科斯脊及西北侧海山群对中美洲俯冲带南部晚中新世大陆弧火山岩的形成起重要作用。以上模型均将中美洲俯冲带南部碱性火山岩和埃达克质岩的形成与科科斯脊俯冲联系起来，但是板块重建和低温热年代学研究指出，科科斯俯冲时间不可能早于3Ma(MacMillanetal., 2004; Morelletal., 2012)。近期，对334和344航次上覆板块钻孔(站位U1380和U1379)的沉积层序和古地磁研究将该时间进一步确定在1.8Ma(Lietal., 2018)。因此，Morelletal.(2012)和Morell(2015)认为中美洲俯冲带南部晚中新世火山岩岩性变化与俯冲样式的变化有关，而科科斯脊与海山俯冲不对该时期火山岩的变化产生影响。Goss and Kay(2006)从地球化学角度证明中美洲俯冲带南部埃达克质岩的形成与俯冲剥蚀有关。同时，Vannucchietal.(2013, 2016b)从构造地质学和地球物理的角度证明哥斯达黎加处确实存在大量的俯冲剥蚀。然而，目前造成上述争论的主要原因是人们目前仍未直接获得上覆板块基底组分。本文对上覆板块沉积物钻孔中粗碎屑层位的研究表明，上覆板块基底组分可能为加勒比大火成岩省基底组分。因此，目前中美洲俯冲带南部俯冲输入端元包含俯冲科科斯脊，俯冲海山，俯冲剥蚀地壳物质(加勒比大火成岩省基底)以及俯冲沉积物。

4.2.2 俯冲带物质再循环及其对中美洲俯冲带南部大陆弧火山岩成因的启示

在俯冲过程中，俯冲输入物质随着俯冲深度的增加会释放流体、熔体或者超临界流体，这些物质都会加入岛弧岩浆的源区，对岛弧火山岩的化学组分造成影响(Spandler and Pirard, 2013; Plank, 2014; 赵仁杰等, 2020; 张泽明等, 2020)。放射性同位素通常在岩浆演化过程中不发生分馏，因此被广泛应用于示踪地壳和地幔之间的物质循环(Hauffetal., 2003; Carpentieretal., 2008; Chauveletal., 2009; Yanetal., 2012)。我们根据所有可能的俯冲输入端元重新探讨中美洲俯冲带南部大陆弧火山岩晚中新世的岩性变化的原因，根据俯冲输入与大陆弧输出的同位素图解(图5)可以看出，中美洲俯冲带南部晚中新世以来的大陆弧火山岩受到了加勒比大火成岩省基底端元和海山端元的影响。根据元素在俯冲过程中的行为，由于轻稀土元素(Nd)一般只在熔体中迁移(Chauveletal., 2009)，我们推测俯冲剥蚀物质以熔体的形式影响该区域大陆弧岩浆作用。不同于Gazeletal.(2011)的研究，我们选择加勒比火成岩省基底代替俯冲科科斯脊的端元：首先科科斯脊俯冲时间存在争议，其俯冲时间可能为更晚的更新世；其次结合Yan and Shi(2016)对344航次钻孔获得的科科斯脊基底岩石同位素分析，以及前人发表的同位素数据(Hoernleetal., 2000; Werneretal., 2003)，我们发现新获得的样品Nd同位素范围更大，且具有更高的Sr同位素及低的208Pb/204Pb和206Pb/204Pb比值，使得科科斯脊和海山的同位素特征并不能解释中美洲俯冲带南部晚中新世以来火山岩的同位素组成。

图5 中美洲俯冲带南部俯冲输入与大陆弧输出物质循环同位素图解

同时，由于特征元素在俯冲系统中具有不同的行为，也常被用于示踪俯冲物质循环(Pearceetal., 2005; Yanetal., 2019)。在俯冲过程中，Ba等大离子亲石元素在俯冲流体中迁移，而Th、Nb以及轻稀土则只在熔体中移动(Johnson and Plank, 1999)。因此，特征元素比值可以用来辨别俯冲组分(Pearceetal., 2005; Pearce and Stern, 2006; Zhangetal., 2019; Yanetal., 2019)。为了消除结晶分异的影响，利用Ba/Th示踪富水流体，利用Th/Nd来辨别含水熔体，利用Ba/Nb来计算俯冲组分释放的总流体(Pearceetal., 2005; Yanetal., 2019)。根据特征元素和同位素比值图解中可以看出，俯冲沉积物对中美洲俯冲带南部大陆弧火山岩产生影响(图6a, b)；同时，在晚中新世中美洲俯冲带南部大陆弧火山岩岩性发生变化时，其Th/Nd比值显著升高，该变化可能是由于俯冲熔体组分的变化引起的(Zhangetal., 2019; Yanetal., 2019)。在所有俯冲输入端元中，俯冲沉积物具有最高的Th/Nd比值，因此我们认为俯冲沉积物以熔体的形式对该地区火山岩产生影响，该结论与前人根据热力学模型计算的哥斯达黎加弧下俯冲沉积物发生部分熔融的认识相吻合(Peacocketal., 2005)。尽管前人通过Be和Tl同位素的研究指出，中美洲俯冲带南部俯冲沉积物对大陆弧火山岩的影响有限(Teraetal., 1986; Nielsenetal., 2017)，但是从特征元素比值可以看出俯冲沉积物的确参与大陆弧物质再循环过程。相对其他端元，俯冲沉积物中富集Ba、Th等特征元素，因此其加入很小的量足以影响该地区火山岩微量元素比值，但是其加入量不足以改变其同位素特征(Gazeletal., 2009; Nielsenetal., 2017; Zhaoetal.in print)。同时，我们并不能排除俯冲科科斯脊对近期中美洲俯冲带南部火山岩的影响，但是综合科科斯脊俯冲时间，我们推测晚中新世以来火山岩中类似加拉帕戈斯热点的物质主要受俯冲剥蚀物质的影响(Goss and Kay, 2006)。在晚中新世(～9Ma)，巴拿马三联点(Panama Triple Junction)位于尼克亚半岛的离岸处，由于俯冲倾角和汇聚速率变化，巴拿马三联点向东南迁移，由于俯冲角度的改变，俯冲板块释放的俯冲组分减少，造成中美洲俯冲带南部钙碱性火山岩的停止喷发(MacMillanetal., 2004; Morelletal., 2012; Morell, 2015)。同时，在其向东南迁移过程中造成明显的俯冲剥蚀(MacMillanetal., 2004; Saketal., 2009)，该过程将大量的弧前物质带入地幔楔，其与俯冲海山及俯冲沉积物共同影响中美洲俯冲带南部大陆弧火山岩的岩浆过程。尽管本文从地球化学角度初步证明了俯冲剥蚀对中美洲俯冲带南部大陆弧岩浆作用的影响，但是并未定量解释不同端元对中美洲俯冲带南部不同岩性火山岩的影响，未来的工作中，我们将结合非传统同位素(Li，Ba等)共同探讨中美洲俯冲带南部俯冲循环过程。

图6 中美洲俯冲带南部俯冲输入与大陆弧输出物质循环同位素及特征元素比值图解

5 结论

本研究对位于汇聚板块边缘上覆板块弧前钻孔中沉积物的粗碎屑层位进行了详细的全岩主、微量元素和Sr-Nd-Pb-Hf同位素分析研究，获得了如下三点新认识：

(1)研究中的粗碎屑层位组分可代表上覆板块基底组分，是剥蚀型汇聚板块边缘中物质俯冲再循环的一个重要端元，因此，弧前钻孔中沉积物的粗碎屑层位组分可为研究浅部和深部地幔地球化学变化提供重要线索。

(2)结合中美洲俯冲带南部大陆弧火山岩化学组分，我们识别出了俯冲剥蚀物质对该大陆弧火山岩成因的贡献，为俯冲剥蚀机制的存在提供了地球化学制约。

(3)自晚中新世以来，巴拿马三联点向东南的迁移引起中美洲俯冲带南部弧前大量的俯冲剥蚀，可能导致俯冲剥蚀物质、海山及俯冲沉积物共同进入了大陆弧地幔源区。

致谢感谢“决心号”科学钻探船上的全体船员、IODP-USIO技术人员以及中国IODP办公室。在研究期间，与哥伦比亚大学Susanne Staub教授进行了有益的讨论。同时，感谢两位匿名审稿人及俞良军老师对本文提出的宝贵修改意见。