LA-ICP-MS zircon U-Pb dating of the undeformed granite in eastern Taiwan and its geological significance
-
摘要:
台湾源头山位于台湾岛东部太鲁阁带的北部,主体为花岗闪长岩-辉长岩"浆混岩"组合,发生了强烈变形变质,其中央分布有未变形的侵入体,具有典型岩浆岩特征,代表了变形变质期后的岩浆事件。通过开展野外地质调查和岩石学、同位素年代学研究,揭示了变形变质作用前后侵入岩在岩石特征及主量、微量元素方面存在明显差别,表明二者是不同构造岩浆事件的产物。未变形的侵入岩为石英闪长岩,具有高Pb、Rb、U、Th、Zr轻度富集特征,属典型的I型花岗岩,其LA-ICPMS锆石U-Pb年龄为82±1Ma,属晚白垩世坎潘期。这是台湾东部地区变形变质事件结束的标志,其后台湾东部进入了一个非挤压的新构造作用阶段。
-
关键词:
- 台湾东部 /
- 块状石英闪长岩 /
- LA-ICP-MS锆石U-Pb定年 /
- 晚白垩世坎潘期
Abstract:Yuantoushan is located in the north of Tailuko belt, eastern Taiwan. It is mainly composed of strongly deformed and metamorphosed granodiorite-gabbro "slurry mixed rock" combination. Deformed intrusions are in the center and show typical magmatic features, representing magmatic events after deformation and metamorphism. The authors carried out field geological survey and petrologic as well as isotopic geochronologic studies of metamorphic intrusive rocks and intrusive rocks. It is shown that there are significant differences in rock characteristics and main trace elements between the two kinds of rocks, indicating that the two kinds of rocks are products of different tectonic magmatic events. Undeformed intrusive rock is quartz diorite, with high Pb, Rb, U, Th, Zr, exhibiting mild enrichment characteristics of I type granite, and its LA-ICP-MS zircon U-Pb isotopic age is 82 ±1Ma, suggesting the Late Cretaceous Campanian period. This indicates that the deformation and metamorphism event in eastern Taiwan were ended, and the evolution then entered a non-squeezing neotectonic stage.
-
台湾源头山花岗岩体位于台湾东部的太鲁阁带北段,是太鲁阁带花岗岩类的重要岩体[1],长期以来一直作为变形变质地质体①,受到了广泛关注。地质研究中将片麻状侵入岩作为“开南冈组”进行填图,由于普遍发育矿物定向排列而成的面理和线理,曾被作为变质作用的典型代表,认为其为晚中生代片麻岩;近年来由于地质调查的不断深入,以及获取的锆石精确年龄等新资料,基本确定片麻岩是侵入岩体,时代为85~90Ma[2],属晚白垩世。Yui等[3-4]建立了台湾东部地区的构造演化模型,提出了3期构造,第一期为碰撞过程,卷入的地质体时代为晚二叠世-晚白垩世早期(大南澳片岩等),形成太鲁阁带,其形成时代为晚白垩世早期;第二期为俯冲过程,形成玉里带内的高压低温变质岩,可能发生于晚白垩世中期;第三期为碰撞过程,发生在上新世-更新世。第一期的变形变质作用十分强烈,其结束性标志的花岗岩侵入时代为85~90Ma。但已有测年的侵入岩均强烈变形变质,并卷入了一次十分强烈的构造过程,具强烈片麻状构造,此次构造作用常被归为第一期构造运动[5],并称为南澳运动。本文以位于源头山岩体的核心部之块状侵入岩为研究对象,采用野外地质调查,并对采集样品进行主量、微量元素分析及锆石LA-ICP-MS定年,将岩浆岩特征与构造作用相联系,以确定第一期构造作用结束的特征性标志。
1. 地质概况
台湾东部地区属于欧亚板块的边缘,北与太平洋板块,南与菲律宾洋板块接壤,以台东纵谷地壳对接带为界,西部属欧亚大陆板块,东北部为太平洋板块,东南部属菲律宾海板块(图 1)[1]。
图 1 台湾东部地区地质略图(据参考文献[1]修改)1-前古近纪(南澳片岩),变粒岩、云母片岩、绿泥石片岩、绿泥钠长片岩等;2-晚白垩世角闪片岩(枫树山岩体);3-源头山岩体(未变形石英闪长岩);4-源头山岩体(片麻状花岗闪长岩);5-绿泥片岩;6-黑云片岩和二云片岩;7-角闪片岩;8-片麻状花岗闪长岩;9-石英闪长岩;10-糜棱岩化;11-取样点Figure 1. Geological map of eastern Taiwan源头山岩体由颜沧波[6]创立,始称为源头山片麻岩,标准地为宜兰县南澳乡南澳村北北西方位的源头山,本次考察其建立地质单位的标准剖面。岩体以片麻状花岗闪长岩-辉长岩组成的浆混岩为主,已经强烈变形,主要表现为糜棱岩化、片麻理化等。在岩体中央部位分布未变形的石英闪长岩,与围岩呈突变关系,其中见有片麻状花岗闪长岩捕虏体,分布在接触带附近,以1~5cm的碎块状散布于边部。早期研究者多认为,未变形的块状石英闪长岩是变形的核心部分,保留了原岩面貌,并作为变形不均一的标志性特征①。本次观察到了未变形侵入岩体与变形的岩体存在明显差别,可能是不同期次侵入作用的产物。源头山岩体的北东部侵入至枫树山角闪岩体,枫树山角闪岩原岩可能为辉长岩,在南澳一带见有脉状体侵入至东澳片岩,后者由黑云片岩、绿泥石片岩、二云石英片岩等组成,原岩可能为含基性火山岩的陆源碎屑岩。
细粒石英闪长岩分布于源头山片麻岩体的中央部位,代号为δοK2,呈较小的岩体,岩石呈浅灰色,细粒花岗结构,块状构造(图版Ⅰ-a、d)。岩石由石英(10%)、斜长石(55%)、角闪石(20%)、黑云母(15%)组成。各矿物分布较均匀,杂乱排列,矿物粒径在0.2~2mm之间。斜长石为灰白色,半自形板状,聚片双晶、卡钠复合双晶、肖钠双晶发育。斜长石被少量绢云母矿物交代,晶体中包裹黑云母、角闪石。角闪石为普通角闪石,呈他形不规则状,切面2组解理完全,其夹角为124°与56°。角闪石Ng=褐黄色、Np=黄色,具双晶,斜消光,最高干涉色为二级蓝。黑云母呈黑色,鳞片状,镜下呈绿黄色,已被绿泥石(干涉色低)取代。石英呈等轴粒状或不规则状充填在其他矿物晶粒间,具波状消光。副矿物有钛铁矿、磷灰石。
源头山岩体的主体为片麻状-条带状花岗闪长岩,含辉长岩脉体或包体。因强烈变形,辉长岩呈长条状脉体或长椭圆形包体,含量约5%,原岩具有浆混杂岩特征(图版Ⅰ-b、e)。该特征与福建沿海的浆混杂岩十分相似,参与调查的两岸专家形成共识,认为源头山岩体由来自地幔的基性岩浆(辉长岩等)和来自地壳的花岗闪长质主体基本同时在相同空间形成。岩石受到强烈的糜棱岩化,呈条带状或片麻状构造,辉长岩呈脉体形成紧密褶皱或揉皱,并具糜棱岩特征,说明岩体形成之后经历了糜棱岩化过程。花岗闪长岩呈片麻状或条带状,由石英(20%)、钾长石(15%)、斜长石(55%)、角闪石(5%)、黑云母(10%)组成。石英遭受韧性剪切力作用,形态拉长,具定向排列,构成流状构造,形成粒径在0.03~0.3mm之间的复晶石英,波状消光。钾长石受强韧性剪切力作用,发生脆性破碎,形成碎斑(具格子状双晶)、糜棱基质。斜长石受强韧性剪切力作用,显脆性破碎,形成碎斑(具聚片双晶、卡钠复合双晶)、糜棱基质。黑云母为鳞片状,具定向排列,Ng=深黄色、Np=浅黄色。部分黑云母被白云母交代。角闪石为普通角闪石,呈柱状或他形。次生蚀变矿物有绢云母(30%)、绿帘石(3%)。辉长岩呈灰色,具脉状或条带状,脉宽3~10cm,主要由普通辉石(30%)、斜长石(30%)、角闪石(30%)和黑云母(10%)组成。
枫树山角闪(片)岩,岩石呈灰绿色,柱状变晶结构,片状构造(图版Ⅰ-c、f)。岩石由角闪石(70%)、斜长石(28%)和石英(2%)组成,少量绿帘石、十字石,副矿物为榍石等。柱状矿物角闪石平行定向排列,构成片理构造。角闪石为柱状变晶结构,呈柱状、长柱状,矿物粒径0.05~8mm,以0.3~ 2mm为主。角闪石横切面2组解理完全,其夹角为124°与56°。角闪石为绿黄色,多色性不明显,斜消光,中正突起,干涉色鲜艳。斜长石为不规则状外形,矿物粒径为0.1~2.6mm。斜长石中正突起,具聚片双晶、卡钠复合双晶,部分颗粒具环带结构,晶体中包裹有角闪石矿物。石英呈不规则状外形,矿物粒径在0.1~0.5mm之间。
2. 样品特征及测试方法
本次取样均在台湾南澳地区,包括源头山岩体(片麻状二长花岗岩和块状石英闪长岩)和枫树山岩体(角闪片岩)(表 1),位置分别为南澳乡的北溪和南澳附近的公路壁。
表 1 样品位置及岩性Table 1. Locations and lithology of chemically and chronologically analyzed samples样号 东经 北纬 岩石名称 采样地点 D0015 121°45′30″ 24°30′29″ 块状石英闪长岩 南澳北溪 D0017 121°46′09″ 24°30′12″ 片麻状糜棱岩化花岗闪长岩 南澳北溪 D0011-1 121°53′08″ 24°30′10″ 角闪片岩(枫树山岩体) 南澳 D0013-1 121°53′08″ 24°30′36″ 微晶角闪片岩(枫树山岩体) 南澳 取样是在台湾地质专家参与下完成的,属典型的岩石露头及剖面路线。D0015及D0017分别为块状石英闪长岩和片麻状花岗闪长岩,均采自河床的新鲜岩石;D0011和D0013-1取自南澳乡的海岸公路。经过详细的野外地质观察和镜下观察,并配套采集了全岩主量、微量和稀土元素分析样品,由福建省地质测试研究中心完成。全岩主量元素分析采用传统方法(符合GB/T14506.28-2010及DZG20.01-1991),其余按定量光谱分析(符合DD2005-01标准)。
由于片麻状花岗闪长岩已经有锆石定年成果[4],故本次只对块状石英闪长岩采锆石定年样。在野外考察和室内研究的基础上,挑选锆石、进行制靶、阴极发光(CL)、透射光和反射光照相,并圈出锆石测点,确定锆石测点的岩浆作用成因。锆石LA-ICP-MS定年测试工作在中国冶金地质总局山东局测试中心完成。
锆石微量元素含量和U-Pb同位素定年利用LA-ICP-MS分析完成。激光剥蚀系统为Newwave UP213,ICPMS为Bruker M90。激光剥蚀过程中采用氦气作载气、氩气为补偿气以调节灵敏度,二者在进入ICP之前通过一个Y型接头混合。每个时间分辨分析数据包括15~20s的空白信号和45s的样品信号。对分析数据的离线处理(包括对样品和空白信号的选择、仪器灵敏度漂移校正、元素含量及U-Th-Pb同位素比值和年龄计算)采用软件ICPMSDataCal[7-8]完成。
3. 分析结果
3.1 岩石化学及地球化学分析结果
本次分析结果包括源头山的块状石英闪长岩、强烈糜棱岩化花岗闪长岩、枫树山角闪片岩等岩石的主量、微量元素数据。其中2806、NA6A、NA8部分来自陈肇夏②的台湾变质岩的资料。枫树山岩体、源头山未变形岩体和强烈变形岩体主量元素、CIPW标准矿物及岩石化学指数均存在差别。①主量元素在上述三者之间存在显著差别:如SiO2分别为48.68%~54.7%、56.7%~58.25%和64.88%~65.52%;Ai2O3为15.43%~15.43%、16.76%~16.99%和16.64%~ 17.18%,属铝不饱和至过饱和;FeO为6.24%~6.38%、5.44%~6.16%和2.74~3.66%;MgO为4.13%~9.32%、2.16%~4.04%和1.89%~1.92%;K2O为0.17%~1.02%、1.62%~2.1%和2.7%~4.31%。②CIPW标准矿物石英Q、刚玉C、钾长石Or、钠长石Ab、钙长石An等均存在显著差别,如石英Q分别为0~15%、12%~14%和21%~30%、刚玉C为0、0~0.6%、3.6%~6.7%。③岩石化学指数的差别也较显著,如分异指数DI分别为53%~71%、77%~80%和87%~89%;A/CNK值分别为0.62~0.85、0.96~0.98和1.23~1.59,准铝质-过铝质;组合指数σ分别为0.88~1.05、1.63~1.88和1.33~ 1.88,但均小于4,属太平洋型;碱质及钾钠比趋于增高,即Na2O+K2O和K2O/Na2O指数均趋于增高。④在TAS图解[9-10](图 2-a)中,样品均投影在虚线下方,属亚碱性系列,枫树山岩体的样品投影在B(玄武岩)和O1(玄武安山岩)区,源头山未变形岩体样品投在O1(玄武安山岩)区和O2(安山岩),源头山强烈变形的岩石样品投影在O3(英安岩)区。硅-钾图解(图 2-b)中按邱家骧等[11]的分类图解,枫树山岩体的样品投影在玄武岩的低钾区和玄武安山岩的中钾区,源头山未变形岩体样品投在玄武安山岩区与安山岩区与中钾/高钾的界线附近,源头山强烈变形的岩石样品投影在高钾的英安岩区与流纹岩区内。
微量元素地球化学按李昌年[12]的划分方案,枫树山角闪(片)岩、源头山未变形石英闪长岩和强烈变形的花岗闪长岩,在微量元素组合上具有如下特征:①过渡元素V、Cr、Co、Ni趋于减少,如Co含量分别为32.14×10-6~48.28×10-6、17×10-6~25.65×10-6和5×10-6~11.32×10-6。②亲石元素中大离子亲石元素(LILE) Rb、Sr、Ba总体趋于增高,如Rb含量分别为5.6×10-6~60.18×10-6、65.87×10-6~89×10-6和114× 10-6~155×10-6;放射性生热元素U、Th也趋于显著增高,如U含量分别为0.09×10-6~0.45×10-6、1.03×10-6~2.1×10-6和2.8×10-6~3.38×10-6;不相容元素Zr和Hf呈成倍增长趋势,如Zr含量分别为62.7×10-6~85.6× 10-6、112×10-6~130.6×10-6和217.2×10-6~346×10-6,HREE(重稀土元素)也趋于增高,但变化较小;亲岩浆元素Th、La、Ce明显趋于增高,如La含量分别为3.26×10-6~7.66×10-6、20.6×10-6~22.62×10-6和37.4× 10-6~44.33×10-6;过亲岩浆元素Cr、Ni、Co明显趋于下降,如Co含量分别为32.14×10-6~48.28×10-6、17× 10-6~25.65×10-6和5×10-6~11.32×10-6。③岩石微量元素蛛网图(图 3-a)表明,枫树山角闪(片)岩、源头山未变形石英闪长岩和强烈变形的花岗闪长岩之间表现为由分散至集中的变化特点,Rb~Hf投影曲线呈分散状,形态相似,具Pb显著富集和Rb、Th、U弱富集,Nb、Ba弱亏损的特点;Sm~Lu近于平行,总体变化不大。
枫树山角闪(片)岩、源头山未变形石英闪长岩和强烈变形的花岗闪长岩的稀土元素含量具如下特征。①轻稀土元素(La~Eu)显著趋增,如Ce含量分别为8.58×10-6~17.48×10-6、28.1×10-6~46.61×10-6和79.7×10-6~90.9×10-6;②重稀土元素(Gd~Yb,Y)差别不大,如Y含量为15×10-6~30.9×10-6;③特征值差别显著,δEu趋于减少,枫树山角闪(片)岩为1,没有亏损,源头山未变形石英闪长岩为0.84~0.89,呈微弱亏损,而源头山强变形花岗闪长岩为0.58~ 0.66,为弱亏损。稀土元素总量及轻、重稀土元素比值显著增高;④枫树山角闪(片)岩的稀土元素配分模式(图 3-b)基本呈水平状,源头山未变形侵入岩和强烈变形侵入岩的稀土元素配分模式均呈左高右低的曲线。模式曲线近于直线,基本符合“I型”花岗岩特征,属壳-幔混合型。
3.2 锆石定年结果
锆石特征、测点位置及Th/U值,均说明锆石定年基础良好、结果可靠(图 5)。锆石呈浅橙黄色、淡玫瑰色,棱角-次棱角状、柱状、长柱状,金刚光泽,透明,晶体表面光洁明亮,部分表面较粗糙含铁染,部分含气、固体包裹体。粒径主要为0.08~0.15mm,次要为0.16~0.25mm,个别为0.26~0.35mm,延长系数主要为1: 2.5~1: 4,次要为1: 4~1: 5,个别1: 7~1: 9。主要晶体类型由柱面{110}不发育的{100}和锥面{311}{131}组成聚形,呈(可能是加工所致)棱角-次棱角状、柱状、长柱状;个别晶体类型由柱面{110}不发育的{100},锥面{111}{311}{131}组成聚形。CL图像(图 4)表明,锆石晶形不完整,多为碎屑状,显示锆石为长柱状等,振荡生长环带不发育,呈灰色-灰白交替出现的纹理,一般出现3条纹理,边界尚平直,连续性良好,总体呈单一特征,说明锆石的生长环境变化不大,为岩浆结晶形成。Th/U值为0.44~0.81,一般为0.6~0.7,均大于0.4,属岩浆成因锆石。不同锆石的特征相当均一,无明显差异,是由一次岩浆作用形成的矿物。
表 2 花岗岩主量、微量和稀土元素分析数据Table 2. Major, trace and rare earth elements of granites岩体 枫树山 源头山岩体 岩体 枫树山 源头山岩体 岩性 角闪片岩 微晶闪长岩(未变形) 细粒石英闪长岩(未变形) 片麻状花岗闪长岩 岩性 角闪片岩 微晶闪长岩(未变形) 细粒石英闪长岩(未变形) 片麻状花岗闪长岩 样号 D0011 D0013-1 D0015 2806 NA6A NA8 D0017 样号 D0011 D0013-1 D0015 2806 NA6A NA8 D0017 SiO2 48.68 54.7 56.7 58.25 65.52 64.88 63.25 Nd 7.15 10.96 23.41 24.5 34.7 34.8 37.28 TiO2 0.98 1.09 1.71 1.69 0.7 0.7 0.9 Sm 2.55 3.12 5 5.29 5.98 6.34 7.59 Al2O3 15.44 15.43 16.99 16.76 16.64 17.18 17.05 Eu 0.94 1.11 1.35 1.48 1.18 1.26 1.52 FeO 6.24 6.38 5.44 6.16 3.66 2.74 2.92 Gd 3.12 3.66 4.05 5.44 4.6 6.88 6.38 Fe2O3 2.1 3.45 1.66 1.06 1.1 1.03 1.86 Tb 0.6 0.64 0.64 0.84 0.83 0.63 1.05 MnO 0.079 0.11 0.124 0.15 0.11 0.03 0.137 Dy 3.24 3.44 3.04 4.95 MgO 9.32 4.13 4.04 2.16 1.89 1.89 1.92 Ho 0.76 0.84 0.7 1.12 CaO 11.46 7.43 5.8 5.37 2.28 1.1 2.37 Er 2.26 2.37 2.02 3.13 Na2O 2.27 2.18 3.11 3.26 2.78 2.52 3.49 Tm 0.32 0.33 0.3 0.44 K2O 0.17 1.02 1.62 2.1 2.7 4.31 3.53 Yb 2.25 2.28 1.96 2.94 4.2 1.18 2.95 P2O5 0.08 0.11 0.12 0.5 0.09 0.05 0.14 Lu 0.3 0.31 0.27 0.45 0.72 0.13 0.39 总和 99.59 99.62 99.64 99.66 Y 23.1 24.6 20.5 18 29 15 30.9 Co 48.28 32.14 25.65 17 10 5 11.32 AR 1.20 1.33 1.52 1.64 1.82 2.19 2.13 Ni 172.4 15.39 15.51 7 22 9 18.12 SI 46.4 24.1 25.5 14.7 15.6 15.1 14.0 Zn 68.93 100.9 91.09 116 141 97 89.12 FL 17.6 30.1 44.9 50.0 70.6 86.1 74.8 Rb 5.6 60.18 65.87 89 114 155 153.6 A/CNK 0.62 0.85 0.98 0.96 1.43 1.59 1.23 Sr 126 222 306.8 337 252 192 237.6 σ 1.05 0.88 1.63 1.88 1.33 2.13 2.43 Nb 2.23 4.17 9.07 18 15 18 15.31 Na2O+K2O 2.44 3.20 4.73 5.36 5.48 6.83 7.02 Cd 0.243 0.234 0.202 0.161 K2O/Na2O 0.07 0.47 0.52 0.64 0.97 1.71 1.01 Ba 34.03 129.2 303.2 396 546 949 616.5 DI 53.35 71.12 77.21 80.29 86.97 88.88 87.59 Hf 2.92 3.03 5.6 4.02 7.01 9.8 7.76 LREE/HREE 1.66 2.11 4.13 3.89 5.04 8.46 4.54 W 0.48 0.73 1.64 1.2 ΣREE 23.88 42.70 104.63 79.97 158.96 177.63 181.49 Pb 10.03 12.46 14.13 16 35 17 21.22 δEu 1.02 1.00 0.89 0.84 0.66 0.58 0.65 Th 1.04 1.98 7.86 6.9 15 19 17.74 Rb/Sr 0.04 0.27 0.21 0.26 0.45 0.81 0.65 U 0.09 0.45 1.03 2.1 3.7 2.8 3.38 Th/U 11.56 4.40 7.63 3.29 4.05 6.79 5.25 Cr 366.5 38.9 36.6 17 56 52 47.7 Ba/Sr 0.27 0.58 0.99 1.18 2.17 4.94 2.59 V 197.4 261 138.4 56 89 86 73.5 Q_CIPW 0.00 15.12 11.95 14.70 30.12 28.30 21.23 Zr 62.7 85.6 130.6 112 249 346 217.2 C_CIPW 0.00 0.00 0.00 0.57 5.35 6.73 3.60 La 3.26 7.66 22.62 20.6 37.4 44.33 42.27 Or_CIPW 1.04 6.28 9.84 12.74 16.37 26.42 21.38 Ce 8.58 17.48 46.61 28.1 79.7 90.9 83.15 Ab_CIPW 19.84 19.21 27.04 28.30 24.13 22.11 30.27 Pr 1.4 2.37 5.64 9.68 An_CIPW 32.47 30.52 28.37 23.98 11.00 5.32 11.12 注:主量元素含量单位为%,微量和稀土元素含量为10-6 锆石各测点的稀土元素及其配分模式表明,所测锆石的配分曲线总体形态相近,为左低、右高的重稀土元素相对富集的特征。Eu元素的亏损较显著,各样品的差异不明显,δEu为0.2~0.4。
综上所述,样品的锆石均属岩浆成因锆石,Th/U值大于0.4,稀土元素配分模式高度一致,锆石作为定年对象具有可靠的基础。
表 3表明,30个测点相对应的30组数据具较高的一致性,206Pb/238U表面年龄为77.2~85Ma,多数为80~83Ma,谐和度较高,分析数据的质量高,符合定年要求,即30组均属有效数据。采用图解的方法处理相关数据及求得年龄值。因年龄小于1000Ma,采用锆石206Pb/238U表面年龄并求其加权平均值(图 5)。
表 3 锆石U-Th-Pb同位素分析数据和计算的年龄值Table 3. Analytical data of zircon U-Th-Pb isotopes and calculated age values测点号 总Pb/10-6 232Th/10-6 238U/10-6 Th/U 同位素比值 表面年龄/Ma 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ D0015-1 6 259 319 0.81 0.0489 0.0026 0.0878 0.0047 0.0130 0.0002 142.7 124.1 85.4 4.4 83.2 1.3 2 4 160 222 0.72 0.0527 0.0034 0.0918 0.0057 0.0129 0.0003 322.3 146.3 89.1 5.3 82.7 1.7 3 13 435 753 0.58 0.0482 0.0019 0.0831 0.0031 0.0124 0.0002 109.4 90.7 81.0 2.9 79.8 1.1 4 5 201 343 0.59 0.0449 0.0023 0.0776 0.0041 0.0126 0.0002 75.9 3.8 80.6 1.6 5 8 334 490 0.68 0.0487 0.0021 0.0824 0.0037 0.0124 0.0003 131.6 101.8 80.4 3.4 79.3 1.8 6 8 354 436 0.81 0.0490 0.0019 0.0860 0.0034 0.0126 0.0002 150.1 88.9 83.8 3.2 80.9 1.2 7 4 122 220 0.56 0.0419 0.0034 0.0758 0.0058 0.0132 0.0003 error error 74.2 5.5 84.6 1.8 8 6 241 336 0.72 0.0477 0.0031 0.0841 0.0052 0.0129 0.0002 87.1 148.1 82.0 4.9 82.9 1.5 9 8 377 459 0.82 0.0483 0.0021 0.0855 0.0036 0.0129 0.0002 122.3 -93.5 83.3 3.4 82.7 1.3 10 6 171 334 0.51 0.0449 0.0024 0.0802 0.0040 0.0132 0.0002 error 78.4 3.8 84.8 1.5 11 5 229 298 0.77 0.0476 0.0023 0.0842 0.0040 0.0129 0.0002 76.0 114.8 82.1 3.8 82.5 1.5 12 4 149 220 0.68 0.0498 0.0028 0.0891 0.0047 0.0131 0.0002 183.4 126.8 86.7 4.4 84.1 1.5 13 6 142 349 0.41 0.0483 0.0032 0.0822 0.0052 0.0125 0.0002 122.3 142.6 80.2 4.9 79.8 1.5 14 7 293 391 0.75 0.0527 0.0028 0.0899 0.0045 0.0125 0.0002 322.3 122.2 87.4 4.2 79.8 1.3 15 9 404 549 0.74 0.0485 0.0026 0.0805 0.0041 0.0122 0.0003 124.2 131.5 78.6 3.8 77.9 2.1 16 5 178 280 0.64 0.0476 0.0031 0.0802 0.0052 0.0122 0.0002 76.0 157.4 78.3 4.9 78.0 1.5 17 5 170 273 0.62 0.0524 0.0027 0.0938 0.0046 0.0131 0.0003 301.9 114.8 91.0 4.2 84.1 1.6 18 3 126 196 0.64 0.0500 0.0034 0.0846 0.0051 0.0126 0.0003 194.5 157.4 82.5 4.8 81.0 1.6 19 4 160 258 0.62 0.0520 0.0032 0.0884 0.0051 0.0127 0.0003 283.4 110.2 86.0 4.8 81.1 1.6 20 6 223 364 0.61 0.0512 0.0029 0.0850 0.0048 0.0121 0.0002 250.1 126.8 82.8 4.5 77.2 1.4 21 8 380 448 0.85 0.0500 0.0026 0.0911 0.0046 0.0133 0.0002 194.5 122.2 88.6 4.3 84.9 1.3 22 8 343 477 0.72 0.0507 0.0029 0.0860 0.0052 0.0123 0.0003 233.4 125.0 83.8 4.9 78.6 1.9 23 4 133 221 0.60 0.0548 0.0030 0.0953 0.0056 0.0127 0.0002 466.7 124.1 92.5 5.2 81.1 1.6 24 3 87 154 0.56 0.0599 0.0055 0.1046 0.0085 0.0136 0.0004 598.2 200.0 101.1 7.8 87.2 2.6 25 5 152 343 0.44 0.0485 0.0029 0.0891 0.0049 0.0137 0.0002 120.5 133.3 86.7 4.5 87.4 1.6 26 4 175 266 0.66 0.0513 0.0033 0.0877 0.0054 0.0126 0.0002 257.5 148.1 85.4 5.0 80.5 1.6 27 6 273 356 0.77 0.0497 0.0033 0.0829 0.0050 0.0122 0.0003 183.4 147.2 80.9 4.7 78.5 2.0 28 5 115 305 0.38 0.0894 0.0399 0.1572 0.0693 0.0129 0.0003 1413.0 963.8 148.3 60.9 82.8 2.0 29 5 170 257 0.66 0.0482 0.0036 0.0868 0.0061 0.0133 0.0003 109.4 166.6 84.5 5.7 85.0 2.0 30 5 193 280 0.69 0.0483 0.0030 0.0849 0.0048 0.0131 0.0003 122.3 131.5 82.7 4.5 83.9 1.9 D0015的206Pb/238U年龄加权平均值为82±1Ma(MSWD=2.7),与谐和图(图 5-a)的交点年龄一致,定年结果可靠。源头山岩体未变形侵入岩的形成年龄为晚白垩世(坎潘期)。
4. 讨论
4.1 源头山块状石英闪长岩和强片麻状花岗岩-辉长岩(浆混岩)形成时代
林启文等①对苏澳幅的调查发现,在源头山岩体核部存在未变形的块状侵入岩,当时作为变形变质的核部,因未被变形变质影响,保留了原岩面貌。在野外实地考察时,台湾地学专家也持类似观点。本次取得的定年分析数据表明,其年龄为82± 1Ma(MSWD=2.7),属晚白垩世。Yui等[2]基于锆石LA-ICP-MS分析成果资料,认为大南澳一带的片麻状花岗岩类的年龄为87±1~89±1Ma,而源头山岩体(片麻状花岗闪长岩)为88±1Ma。这表明,块状石英闪长岩虽然位于源头山岩体的中心附近,但其年龄较新,为82±1Ma,与其受到强烈变形变质的主体部分存在约6Ma的时差,该时差虽小,但超过分析数据的误差范围(约1Ma),属具显著差别的不同年龄。地质事实也支持这种认识:①块状石英闪长岩与强片麻状花岗闪长岩-辉长岩呈突变关系,在内接触带还见有后者的捕虏体,表明强片麻状花岗闪长岩-辉长岩形成在先,块状石英闪长岩形成在后。②块状石英闪长岩未受变形变质作用影响,为花岗结构和块状构造;强片麻状花岗闪长岩受到强烈变形变质影响,形成糜棱岩化结构和条带状或片麻状构造,表明二者成岩后经历的构造过程存在明显差异。③岩石地球化学也存在显著差异,块状石英闪长岩与强片麻状花岗闪长岩存在主量、微量和稀土元素方面的差异,前者比后者呈低硅、铝、较低的亲石元素(LILE)Rb、Sr、Ba和稀土元素等。
综上,强片麻状花岗闪长岩与块状石英闪长岩是不同期次的侵入岩,前者岩石形成之后,还遭受了一次十分强烈的糜棱岩化过程,之后块状石英闪长岩侵入,形成源头山块状侵入岩体。
4.2 岩浆岩形成阶段
台湾东部的侵入岩主要分布在大南澳地区,经过许多地质学家的共同努力,已经成功地将原划分为片麻岩的地质体,作为正变质岩归为侵入岩体,主要是源头山岩体、开南岗岩体、溪畔岩体等,Yui等[2]据锆石同位素年龄成果,基本确定了岩浆岩形成的年龄为87±1~89±1Ma,属晚白垩世。该年龄无疑为台湾东部构造岩浆作用的重要依据。以糜棱岩化为界,片麻状花岗岩类属于构造作用之前的侵入体,而源头山岩体内的块状石英闪长岩是构造作用之后的产物,代表了后构造岩浆事件。即在晚白垩世大南澳地区经历了2个阶段的岩浆事件:第一阶段为87± 1~89±1Ma,主要岩性为花岗闪长岩-辉长岩(浆混岩);第二阶段为82±1Ma,主要岩性为石英闪长岩。相应的构造环境为张扭性构造和张性构造环境。
4.3 大南澳地区构造作用期次
本次结果表明,台湾东澳片岩由3个层次的构造层组成:东澳片岩、强烈变形的花岗闪长岩-辉长岩和未变形的石英闪长岩。东澳片岩常作为基底,由绿泥绿帘片岩、绿泥钠长片岩、云母石英片岩、石英片岩夹含石墨云母片岩组成,原岩为玄武岩、副长岩或碱性玄武岩及陆源碎屑岩,岩石主要为碱性系列,具低Ba、Sr和高Pb特征,稀土元素总量普遍较低,具有轻稀土元素富集、Eu无亏损的稀土元素配分模式,属陆源裂陷的火山(岩浆)-沉积作用产物,形成年龄约95.4±4.5Ma(火山碎屑锆石年龄)。强烈变形的花岗闪长岩-辉长岩属正变质岩类,以强烈糜棱岩化和片麻状构造发育为特征,时代为87± 1~89±1Ma,如前所述,岩石具有I型花岗岩特征,属钙碱性系列,是张扭性构造环境的产物,表明可能属地壳增厚下壳幔作用的结果。之后又经历了一次强烈的压性构造过程,形成强烈糜棱岩化,使花岗闪长岩-辉长岩(浆混岩)发生强烈糜棱岩化。至82±1Ma,又进入了张性构造环境阶段,形成未变形的块状石英闪长岩。
将上述内容总结如下:
(1)晚白垩世(塞诺曼期,约95.4±4.5Ma),陆缘裂陷构造环境,形成东澳片岩中的碱性火山岩系列。
(2)晚白垩世(康尼亚克期,约89±1Ma),张扭性构造环境,形成以花岗闪长岩-辉长岩为代表的侵入岩,属于钙碱性系列。
(3)岩浆岩形成之后,还经历了一次强烈的糜棱岩化过程。中深层次的构造作用,使包括花岗闪长岩在内的岩石发生了强烈糜棱岩化。
(4)晚白垩世(坎潘期,82±1Ma),转为张性构造环境,发生小规模的岩浆侵入,形成以块状石英闪长岩为特征的岩浆活动。
综上,晚白垩世(95.4~82Ma),大南澳地区经历了4次构造-岩浆与变形作用,即陆缘裂陷、张扭性构造-岩浆作用、压性构造(糜棱岩化)和张性构造作用4个阶段,这是太平洋板块多次活动的结果。
4.4 闽台构造岩浆事件对比
福建沿海地区的长乐-南澳断裂带活动可能始于晚三叠世-早侏罗世[14-17],经历了多次构造-岩浆事件[18-20],其结束时间约为100Ma。造山后构造岩浆阶段可能始于晚白垩世。长乐-南澳断裂带结束的标志性事件是强烈壳幔作用,发生地壳深部的物质重熔,形成酸性岩浆与来自地幔的基性岩浆上升。2种岩性同时同空间侵位,形成浆混杂岩体。此时,台湾东部地区发生(碱性)基性岩浆活动,并形成南澳绿片岩之原岩。在87~82Ma期间,长乐-南澳断裂带经历了2次构造岩浆事件,即87~84Ma的晶洞花岗岩(莆田东北部晶洞花岗岩为84.5± 1.2Ma)和未变形的变质的中酸性火山岩组合(莆田渡边英安岩约82.9±1.6Ma)。晶洞花岗岩作为典型的A型花岗岩,形成于张性构造环境。而中酸性火山岩组合(流纹岩-英安岩组合)也是张性构造环境的产物。
本次结果表明,台湾东部地区与福建地区具有一定的响应关系。
晚白垩世(塞诺曼期),台湾东部为陆缘裂陷的构造环境,形成东澳片岩中的碱性火山岩系列。而福建沿海地区形成地壳-地幔有关的浆混岩。
晚白垩世(康尼亚克期),台湾东部形成花岗闪长岩-辉长岩组合,如源头山主体、开南岗等。东南沿海则形成晶洞花岗岩、A型花岗岩和张性构造环境下的岩浆岩。
晚白垩世(坎潘期),台湾东部为张性构造环境,形成块状石英闪长岩等。福建沿海地区则见有莆田渡边一带的中酸性火山岩组合(流纹岩-英安岩组合)。
福建沿海地区的构造岩浆作用与台湾东部地区的岩浆活动十分吻合。显然,这不是巧合,是构造岩浆作用的必然,也是欧亚大陆板块边缘与太平洋板块交界处的频繁活动与内陆地区构造岩浆响应的关系。
5. 结论
(1)源头山未变形的块状石英闪长岩与其主体属不同期次构造岩浆活动的产物。
(2)LA-ICP-MS锆石U-Pb定年结果表明,源头山未变形的石英闪长岩的成岩年龄为82±1Ma,属晚白垩世(坎潘期)。
(3)变形变质作用有3个阶段,分别为晚白垩世的塞诺曼期、康尼亚克期和坎潘期。在台湾东部形成大南澳片岩、强片理化花岗闪长岩-辉长岩和块状石英闪长岩;相应的福建沿海形成“浆混杂岩”、晶洞花岗岩和流纹岩-英安岩组合。
致谢: 本文在调研工作中得到福建省地质调查研究院陈润生教授级高级工程师的帮助;台湾大学陈正宏、李寄瑀教授,台湾东华大学刘萤三教授,台湾师范大学叶恩肇教授等参加了野外考察及调查工作,审稿专家对论文提出了宝贵修改意见,在此表示衷心感谢。 -
图 1 台湾东部地区地质略图(据参考文献[1]修改)
1-前古近纪(南澳片岩),变粒岩、云母片岩、绿泥石片岩、绿泥钠长片岩等;2-晚白垩世角闪片岩(枫树山岩体);3-源头山岩体(未变形石英闪长岩);4-源头山岩体(片麻状花岗闪长岩);5-绿泥片岩;6-黑云片岩和二云片岩;7-角闪片岩;8-片麻状花岗闪长岩;9-石英闪长岩;10-糜棱岩化;11-取样点
Figure 1. Geological map of eastern Taiwan
表 1 样品位置及岩性
Table 1 Locations and lithology of chemically and chronologically analyzed samples
样号 东经 北纬 岩石名称 采样地点 D0015 121°45′30″ 24°30′29″ 块状石英闪长岩 南澳北溪 D0017 121°46′09″ 24°30′12″ 片麻状糜棱岩化花岗闪长岩 南澳北溪 D0011-1 121°53′08″ 24°30′10″ 角闪片岩(枫树山岩体) 南澳 D0013-1 121°53′08″ 24°30′36″ 微晶角闪片岩(枫树山岩体) 南澳 表 2 花岗岩主量、微量和稀土元素分析数据
Table 2 Major, trace and rare earth elements of granites
岩体 枫树山 源头山岩体 岩体 枫树山 源头山岩体 岩性 角闪片岩 微晶闪长岩(未变形) 细粒石英闪长岩(未变形) 片麻状花岗闪长岩 岩性 角闪片岩 微晶闪长岩(未变形) 细粒石英闪长岩(未变形) 片麻状花岗闪长岩 样号 D0011 D0013-1 D0015 2806 NA6A NA8 D0017 样号 D0011 D0013-1 D0015 2806 NA6A NA8 D0017 SiO2 48.68 54.7 56.7 58.25 65.52 64.88 63.25 Nd 7.15 10.96 23.41 24.5 34.7 34.8 37.28 TiO2 0.98 1.09 1.71 1.69 0.7 0.7 0.9 Sm 2.55 3.12 5 5.29 5.98 6.34 7.59 Al2O3 15.44 15.43 16.99 16.76 16.64 17.18 17.05 Eu 0.94 1.11 1.35 1.48 1.18 1.26 1.52 FeO 6.24 6.38 5.44 6.16 3.66 2.74 2.92 Gd 3.12 3.66 4.05 5.44 4.6 6.88 6.38 Fe2O3 2.1 3.45 1.66 1.06 1.1 1.03 1.86 Tb 0.6 0.64 0.64 0.84 0.83 0.63 1.05 MnO 0.079 0.11 0.124 0.15 0.11 0.03 0.137 Dy 3.24 3.44 3.04 4.95 MgO 9.32 4.13 4.04 2.16 1.89 1.89 1.92 Ho 0.76 0.84 0.7 1.12 CaO 11.46 7.43 5.8 5.37 2.28 1.1 2.37 Er 2.26 2.37 2.02 3.13 Na2O 2.27 2.18 3.11 3.26 2.78 2.52 3.49 Tm 0.32 0.33 0.3 0.44 K2O 0.17 1.02 1.62 2.1 2.7 4.31 3.53 Yb 2.25 2.28 1.96 2.94 4.2 1.18 2.95 P2O5 0.08 0.11 0.12 0.5 0.09 0.05 0.14 Lu 0.3 0.31 0.27 0.45 0.72 0.13 0.39 总和 99.59 99.62 99.64 99.66 Y 23.1 24.6 20.5 18 29 15 30.9 Co 48.28 32.14 25.65 17 10 5 11.32 AR 1.20 1.33 1.52 1.64 1.82 2.19 2.13 Ni 172.4 15.39 15.51 7 22 9 18.12 SI 46.4 24.1 25.5 14.7 15.6 15.1 14.0 Zn 68.93 100.9 91.09 116 141 97 89.12 FL 17.6 30.1 44.9 50.0 70.6 86.1 74.8 Rb 5.6 60.18 65.87 89 114 155 153.6 A/CNK 0.62 0.85 0.98 0.96 1.43 1.59 1.23 Sr 126 222 306.8 337 252 192 237.6 σ 1.05 0.88 1.63 1.88 1.33 2.13 2.43 Nb 2.23 4.17 9.07 18 15 18 15.31 Na2O+K2O 2.44 3.20 4.73 5.36 5.48 6.83 7.02 Cd 0.243 0.234 0.202 0.161 K2O/Na2O 0.07 0.47 0.52 0.64 0.97 1.71 1.01 Ba 34.03 129.2 303.2 396 546 949 616.5 DI 53.35 71.12 77.21 80.29 86.97 88.88 87.59 Hf 2.92 3.03 5.6 4.02 7.01 9.8 7.76 LREE/HREE 1.66 2.11 4.13 3.89 5.04 8.46 4.54 W 0.48 0.73 1.64 1.2 ΣREE 23.88 42.70 104.63 79.97 158.96 177.63 181.49 Pb 10.03 12.46 14.13 16 35 17 21.22 δEu 1.02 1.00 0.89 0.84 0.66 0.58 0.65 Th 1.04 1.98 7.86 6.9 15 19 17.74 Rb/Sr 0.04 0.27 0.21 0.26 0.45 0.81 0.65 U 0.09 0.45 1.03 2.1 3.7 2.8 3.38 Th/U 11.56 4.40 7.63 3.29 4.05 6.79 5.25 Cr 366.5 38.9 36.6 17 56 52 47.7 Ba/Sr 0.27 0.58 0.99 1.18 2.17 4.94 2.59 V 197.4 261 138.4 56 89 86 73.5 Q_CIPW 0.00 15.12 11.95 14.70 30.12 28.30 21.23 Zr 62.7 85.6 130.6 112 249 346 217.2 C_CIPW 0.00 0.00 0.00 0.57 5.35 6.73 3.60 La 3.26 7.66 22.62 20.6 37.4 44.33 42.27 Or_CIPW 1.04 6.28 9.84 12.74 16.37 26.42 21.38 Ce 8.58 17.48 46.61 28.1 79.7 90.9 83.15 Ab_CIPW 19.84 19.21 27.04 28.30 24.13 22.11 30.27 Pr 1.4 2.37 5.64 9.68 An_CIPW 32.47 30.52 28.37 23.98 11.00 5.32 11.12 注:主量元素含量单位为%,微量和稀土元素含量为10-6 表 3 锆石U-Th-Pb同位素分析数据和计算的年龄值
Table 3 Analytical data of zircon U-Th-Pb isotopes and calculated age values
测点号 总Pb/10-6 232Th/10-6 238U/10-6 Th/U 同位素比值 表面年龄/Ma 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ D0015-1 6 259 319 0.81 0.0489 0.0026 0.0878 0.0047 0.0130 0.0002 142.7 124.1 85.4 4.4 83.2 1.3 2 4 160 222 0.72 0.0527 0.0034 0.0918 0.0057 0.0129 0.0003 322.3 146.3 89.1 5.3 82.7 1.7 3 13 435 753 0.58 0.0482 0.0019 0.0831 0.0031 0.0124 0.0002 109.4 90.7 81.0 2.9 79.8 1.1 4 5 201 343 0.59 0.0449 0.0023 0.0776 0.0041 0.0126 0.0002 75.9 3.8 80.6 1.6 5 8 334 490 0.68 0.0487 0.0021 0.0824 0.0037 0.0124 0.0003 131.6 101.8 80.4 3.4 79.3 1.8 6 8 354 436 0.81 0.0490 0.0019 0.0860 0.0034 0.0126 0.0002 150.1 88.9 83.8 3.2 80.9 1.2 7 4 122 220 0.56 0.0419 0.0034 0.0758 0.0058 0.0132 0.0003 error error 74.2 5.5 84.6 1.8 8 6 241 336 0.72 0.0477 0.0031 0.0841 0.0052 0.0129 0.0002 87.1 148.1 82.0 4.9 82.9 1.5 9 8 377 459 0.82 0.0483 0.0021 0.0855 0.0036 0.0129 0.0002 122.3 -93.5 83.3 3.4 82.7 1.3 10 6 171 334 0.51 0.0449 0.0024 0.0802 0.0040 0.0132 0.0002 error 78.4 3.8 84.8 1.5 11 5 229 298 0.77 0.0476 0.0023 0.0842 0.0040 0.0129 0.0002 76.0 114.8 82.1 3.8 82.5 1.5 12 4 149 220 0.68 0.0498 0.0028 0.0891 0.0047 0.0131 0.0002 183.4 126.8 86.7 4.4 84.1 1.5 13 6 142 349 0.41 0.0483 0.0032 0.0822 0.0052 0.0125 0.0002 122.3 142.6 80.2 4.9 79.8 1.5 14 7 293 391 0.75 0.0527 0.0028 0.0899 0.0045 0.0125 0.0002 322.3 122.2 87.4 4.2 79.8 1.3 15 9 404 549 0.74 0.0485 0.0026 0.0805 0.0041 0.0122 0.0003 124.2 131.5 78.6 3.8 77.9 2.1 16 5 178 280 0.64 0.0476 0.0031 0.0802 0.0052 0.0122 0.0002 76.0 157.4 78.3 4.9 78.0 1.5 17 5 170 273 0.62 0.0524 0.0027 0.0938 0.0046 0.0131 0.0003 301.9 114.8 91.0 4.2 84.1 1.6 18 3 126 196 0.64 0.0500 0.0034 0.0846 0.0051 0.0126 0.0003 194.5 157.4 82.5 4.8 81.0 1.6 19 4 160 258 0.62 0.0520 0.0032 0.0884 0.0051 0.0127 0.0003 283.4 110.2 86.0 4.8 81.1 1.6 20 6 223 364 0.61 0.0512 0.0029 0.0850 0.0048 0.0121 0.0002 250.1 126.8 82.8 4.5 77.2 1.4 21 8 380 448 0.85 0.0500 0.0026 0.0911 0.0046 0.0133 0.0002 194.5 122.2 88.6 4.3 84.9 1.3 22 8 343 477 0.72 0.0507 0.0029 0.0860 0.0052 0.0123 0.0003 233.4 125.0 83.8 4.9 78.6 1.9 23 4 133 221 0.60 0.0548 0.0030 0.0953 0.0056 0.0127 0.0002 466.7 124.1 92.5 5.2 81.1 1.6 24 3 87 154 0.56 0.0599 0.0055 0.1046 0.0085 0.0136 0.0004 598.2 200.0 101.1 7.8 87.2 2.6 25 5 152 343 0.44 0.0485 0.0029 0.0891 0.0049 0.0137 0.0002 120.5 133.3 86.7 4.5 87.4 1.6 26 4 175 266 0.66 0.0513 0.0033 0.0877 0.0054 0.0126 0.0002 257.5 148.1 85.4 5.0 80.5 1.6 27 6 273 356 0.77 0.0497 0.0033 0.0829 0.0050 0.0122 0.0003 183.4 147.2 80.9 4.7 78.5 2.0 28 5 115 305 0.38 0.0894 0.0399 0.1572 0.0693 0.0129 0.0003 1413.0 963.8 148.3 60.9 82.8 2.0 29 5 170 257 0.66 0.0482 0.0036 0.0868 0.0061 0.0133 0.0003 109.4 166.6 84.5 5.7 85.0 2.0 30 5 193 280 0.69 0.0483 0.0030 0.0849 0.0048 0.0131 0.0003 122.3 131.5 82.7 4.5 83.9 1.9 -
福建省地质调查研究院.台湾省区域地质志[M].北京:地质出版社, 2015. Yui T F, Maki K, Lan C Y, et al. Detrital zircons from the Tananao metamorphic complex of Taiwan:Implications for sediment provenance and Mesozoic tectonics[J]. Tectonophysics, 2012, 541/543:31-42. doi: 10.1016/j.tecto.2012.03.013
Yui T F, Wu T W, Jahn B M. Geochemistry and plate-tectonic significance of the metabasites from the Tananao schist complex of Taiwan[J]. Journal of southeast Asian Earth Science, 1990, 4:357-368. doi: 10.1016/0743-9547(90)90006-Y
Yui T F, Lu C Y, Lo A.Tectonic evolution of the Tananao schist Complex of Taiwan[C]//Aubouin A, Bourgois J. Tectonics of Circum-Pacific Continental Margins.VSP, the Nethirlands, 1990:193-209.
Jahn B M, Martineau F, Peucat J J, et al. Geochronology of the Tananao schist complex, Tawan, and its regional tectonic significance[J]. Tectonophysics, 1986, 125:103-124. doi: 10.1016/0040-1951(86)90009-0
Yen T P. Some geological problems on the Tananao Schist[J]. Bull. Geol. Surv. Taiwan, 1954, 7:47-50. http://www.sciencedirect.com/science/article/pii/074395479090006Y
侯可军, 李延河, 田有荣. LA-MC-ICP-MS锆石微区原位U-Pb定年技术[J].矿床地质, 2009, 28(4):481-492. http://d.wanfangdata.com.cn/Periodical/kcdz200904010 Liu Y S, Gao S, Hu Z C, et al. Continental and oceanic crust recycling-induced melt-peridotite interactions in the Trans-North China Orogen:U-Pb dating, Hf isotopes and trace elements in zircons from mantle xenoliths[J]. Journal of Petrology, 2010, 51:537-571. doi: 10.1093/petrology/egp082
Le Maitre R W. A Classification of igneous rocks and glossary of terms[M]. Oxford:Blackwell Sci. Pub., 1989:1-193.
邓晋福, 刘翠, 冯艳芳, 等.关于火成岩常用图解的正确使用:讨论与建议[J].地质论评, 2015, 61(4):717-734. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzlp201504002&dbname=CJFD&dbcode=CJFQ 邱家骧, 林景仟.岩石化学[M].北京:地质出版社, 1991. 李昌年.火成岩微量元素岩石学[M].武汉:中国地质大学出版社, 1992. Masuda A, Nakamura N, Tanaka T. Fine structure of mutually normalized rare earth patterns of chondrites[J]. Geochim. Cosmochim. Acta, 1973, 49:2051-2067. http://www.sciencedirect.com/science/article/pii/0016703773901312
黄长煌.福建东山晚中生代花岗岩的LA-ICP-MS锆石同位素定年及其地质意义[J].福建地质, 2015, 34(4):261-271. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=fjdz201504001&dbname=CJFD&dbcode=CJFQ 黄长煌.福建东山变质岩LA-ICP-MS同位素年龄及其地质意义[J].中国地质, 2016, 43(3):738-750. doi: 10.12029/gc20160304 黄长煌.福建龙海深澳变质岩锆石U-Pb同位素年代学及地质意义[J].华东地质, 2016, 37(1):28-35. http://d.wanfangdata.com.cn/Periodical/hsdzykc201601004 黄长煌.福建港里地区变质岩与南园组火山岩地球化学特征及锆石U-Pb年龄对比[J].华东地质, 2017, 38(1):12-20. http://d.wanfangdata.com.cn/Periodical/hsdzykc201701002 冯艳芳, 邓晋福, 肖庆辉, 等.长乐-南澳构造带花岗岩类年代学岩石组合与构造演化[M].北京:地质出版社, 2013. 毛建仁, 厉子龙, 叶敏.华南中生代构造-岩浆活动研究:现状与前景[J].中国科学:地球科学, 2014, 44(12):2593-2617. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jdxk201412001&dbname=CJFD&dbcode=CJFQ 徐先兵, 张岳桥, 贾东等.华南早中生代大地构造过程[J].中国地质, 2009, 36(3):573-593. http://d.wanfangdata.com.cn/Periodical/zgdizhi200903007 林启文, 高铭建. 1: 5万台湾地质图说明书. 苏澳幅(第二版). 经济部中央地质调查所.2009. 陈肇夏. 台湾的变质岩. 经济部中央地质调查所. 1998.