山东即墨马山粗面英安岩年代学与地球化学特征及其地质意义

朱晓青，侯方辉，刘洪滨，郭兴伟，孙天本，秦亚超，安郁辉，李凤春

1. 中国地质调查局青岛海洋地质研究所，青岛 266071

2. 青岛海洋科学与技术试点国家实验室海洋矿产资源评价与探测技术功能实验室，青岛 266237

3. 中国海洋大学海洋发展研究院，青岛 266100

4. 青岛市即墨区马山自然保护区管理处，青岛 266200

5. 中国冶金地质总局山东局测试中心，济南 250014

岩石圈增厚与下沉是大陆岩浆弧地区普遍存在的旋回性过程，典型例子如：中新生代时期的中国华北东部与美洲科迪勒拉等地区[1]。任纪舜等[2]根据中国东部大兴安岭-太行山-武陵山重力梯度带东、西两侧燕山造山期之后的差异性：西部保存良好的沉积盆地、广泛分布的中元古代—侏罗纪沉积盖层（中-上构造层）迥然不同于东部的变质基底大量出露（下构造层）特征，首次提出中国东部在印支-燕山旋回期间曾因挤压造山而隆升（成为高原）。邓晋福等[3]根据实验岩石学的证据，认为华北燕辽地区燕山期无负Eu异常的火山岩应形成于加厚的地壳底部（或山根带）。董树文等[4-5]提出燕山期华北周缘多板块汇聚的东亚汇聚模式，用于解释华北东部中生代东高西低古地势的动力学机制。张旗等[6]在上述研究基础上，首先提出“C型”埃达克岩的概念，将其用于描述中国东部燕山期富钾并具埃达克岩属性特征的岩浆岩，认为该类型岩浆岩源自增厚地壳的熔融，并根据中生代埃达克岩的分布大致圈定了东部高原的范围。尽管对于“C型”埃达克岩的成因及其所代表的大地构造意义仍存在不同的认识[6-12]，但越来越多的实例证明了“C型”埃达克岩的出现与地壳增厚具有对应关系[8]。

“中国大地构造演化和国际亚洲大地构造图编制（青岛所）”项目组在2017—2019年进行海区毗邻陆域的野外工作中，借助岩石主微量元素特征确定青岛即墨马山发育一套具埃达克岩地球化学特征的火山岩。该套中酸性火山岩因发育罕见的柱状节理而建设成为“马山石林”地质景观。1∶25万青岛市幅区域地质调查成果将该套火山岩归入青山群八亩地组火山岩地层中[13]。韩宗珠等[14]曾根据岩石主微量元素分析结果提出马山安山玢岩的分异结晶程度较高，为具火山弧属性的钙碱性火山岩，推测构造背景为大陆边缘弧环境，火山岩成生的动力学机制可能为扬子板块与华北板块在燕山期的推挤作用。最近，何登洋等[15]也对该套马山火山岩开展了锆石测年、全岩地球化学及锆石原位Lu-Hf同位素等工作，其认为该套火山岩来源岩浆形成于早白垩世（119.3±1.6 Ma）区域下地壳物质的熔融，与胶莱盆地早白垩世伸展变形时间对应，对应的动力学机制为古太平洋板片大规模回撤背景下的区域性伸展。上述研究者前期所做的部分样品显示出典型的“C型”埃达克岩地化特征，但未对其埃达克岩属性展开讨论[14-15]，更大的分歧在于对该套火山岩形成的构造背景存在挤压与伸展环境两种截然相反的认识[14-15]。

基于此，拟以该套马山柱状节理中-酸性火山岩为研究目标，利用其中的锆石U-Pb测年进一步限定其形成时代，借助岩石主微量元素分析与锆石原位Lu-Hf同位素特征探讨其来源岩浆性质，在前人研究的基础上，重点为下述3个科学问题的探讨提供约束与线索：①马山柱状节理火山岩的形成时间；②柱状节理火山岩的来源岩浆的地球化学特征；③该火山岩形成期对应的区域动力学机制等。

1 地质背景与采样

山东即墨马山地区在大地构造位置上处于中朝板块与扬子板块的结合带——苏鲁造山带北侧、胶莱盆地南部，因该地中性火山岩中发育规模壮观的柱状节理及其他地质景观，已建设成为省级地质公园。该地区的柱状节理中-酸性火山岩岩体直接覆盖于早白垩世莱阳群曲格庄组砂岩、粉砂岩之上，属于青山群八亩地组[13]（图1a）。

马山柱状节理火山岩岩体露头呈浅灰绿色，肉眼观察为斑状结构，发育充填方解石的晶洞（图1b，c），晶洞直径由几厘米到十几厘米不等。在马山地质公园的“石林”景观处（36°24′3″N、120°22′17″E）分散采集新鲜的柱状节理火山岩样品4块，每块约5 kg，分别用于薄片磨制、锆石挑选及主微量元素分析等。

图1 研究区地质概况、马山柱状节理火山岩露头与镜下特征a. 研究区地质简图（据文献[13]），b. 马山石林剖面，c. 粗面英安岩中的晶洞构造，d. 马山粗面英安岩镜下显微特征，左侧为单偏光（-），右侧为正交光（+）。Fig.1 Geological map of Mashan region with pictures of outcrops and photomicrographs of the volcanic rocks with columnar joints a. sketched geological map of study area （According to reference [13]）,b. section of columnar joints developed in Mashan,c. geode structure of the trachydacite,d. petrographic characteristics of Mashan trachydacite under microscope,with single polarized light （-）on the left and crossed light （+） on the right.

2 实验方法

2.1 锆石U-Pb、Lu-Hf及微量元素分析

选取样品MS-01用于锆石挑选、测年及薄片磨制。MS-01、MS-02、MS-03与MS-04用于主微量元素分析。锆石挑选、制靶与拍照均在河北廊坊诚信地质服务有限公司进行。按照标准的重矿物挑选流程对样品进行粉碎、重磁分选，双目镜下挑选获得用于定年的锆石颗粒。然后用环氧树脂将锆石固定并研磨抛光至锆石颗粒尺寸的1/2左右，以揭示锆石内部结构。在进行锆石U-Pb及Lu-Hf测试前，首先拍摄锆石的透射、反射及阴极发光图像，以保证在分析过程中避开锆石颗粒中的包裹体或裂隙。

锆石U-Pb测试在中国冶金局山东局测试中心完成，所用激光剥蚀系统为美国Conherent 公司生产的GeoLasPro 193 nm ArF准分子系统，电感耦合等离子体质谱仪的型号为ThermoFisher 公司生产的iCAPQ。具体实验流程为：以氦气为载气，束斑直径为30 μm，频率为10 Hz，能量密度约为8 J/cm2。单点采集时间模式为：30 s气体空白，50 s样品剥蚀，30 s冲洗；每10个样品点插入一组标样（锆石标样+成分标样）。采用91500（年龄为1 064±2 Ma）作为外标进行基体校正；成分标样采用NIST SRM610。其中29Si作为内标元素。锆石微量元素分析与UPb定年同步进行，以NIST 612作为外标计算未知微量元素的浓度，使用Pearce等[16]推荐值。利用29Si作为内标对分析结果进行标准化。轻稀土元素分析结果的平均不确定度为10%，其余元素为5%。

锆石原位Lu-Hf测试在北京科荟测试技术有限公司进行，利用NWR213nm 固体激光器对锆石进行剥蚀，激光剥蚀的斑束直径为55 μm，能量密度为7～8 J/cm2，频率为10 Hz，激光剥蚀物质以高纯He为载气送入Neptune Plus（MC-ICPMS）。以国际标样GJ-1作为监控样。采用179Hf/177Hf=0.732 5对Hf同位素比值进行指数归一化质量歧视校正，采用173Yb/172Yb=1.352 74对Yb同位素比值进行指数归一化质量歧视校正。εHf（t）计算采用衰变常数λ=1.867×10−11/a[17]，（176Lu/177Hf）CHUR=0.033 2，（176Hf/177Hf）CHUR=0.282 772[18]，（176Lu/177Hf）DM=0.038 4，（176Hf/177Hf）DM=0.283 25[19]；TCDM为二阶段Hf模式年龄；计算采用平均地壳176Lu/177Hf比值为0.015。

2.2 岩石主微量元素分析

样品的主量元素与微量元素分析均在中国冶金地质总局山东局测试中心完成。4件样品的制备分为粗碎、中碎和细碎，每个阶段均包括破碎、过筛、混匀和缩分4道工序，实验室检测温度均保持25 ℃。主量元素分析采用X-射线荧光粉末熔片法，具体流程为：准确称取0.200 0 g样品，7.000 0 g混合熔剂放入瓷坩埚中，搅拌均匀，转入铂-金坩埚（wPt95%+wAu5%）内，加入2滴溴化锂溶液（4.2），将坩埚置于熔样机上以1 150 ℃熔融12 min，浇铸制备玻璃片。最后上机（ARL 9900XP型X射线荧光光谱仪）测试主量元素。主量元素的检测精度与准确度为：SiO2、Al2O3、MgO为0.1%，总铁TFe2O3、CaO为0.05%，TiO2、MnO、Na2O、K2O为0.01%，P2O5为0.005%。

微量元素测试采用电感耦合等离子质谱法：准确称取0.100 0 g样品于Teflon杯中，少量水润湿样品，依次加入2 mL HNO3和3 mL HF，加盖，放入钢套，烘箱185 ℃加热48 h后，冷却后拧开钢套，开盖，电热板130 ℃加热蒸至尽干，加入2 mL HNO3继续加热蒸干，以1.5 mL HNO3，1.5 mL MQ水冲洗杯壁，加盖，放入钢套，180 ℃加热12 h提取盐类，冷却后开盖，转移至100 mL 塑料刻度管中，稀释至刻度，摇匀，静置后利用ICP-MS测试微量元素。微量元素的检测精度与准确度为：U为0.003×10−6，Y、Nd、Hf为0.01×10−6，Cd为0.02×10−6，Mo、Be、Zr、Ta为0.05×10−6，Sc、Ni、W、Tl、Pb为0.1×10−6，Co、Cu、Ga、Sr、Sn为0.2×10−6，Cr、Ba为0.5×10−6，Th为0.8×10−6，Li、Rb为1.0×10−6，V、Zn为2.0×10−6。稀土元素的检出限除La、Ce、Pr、Nd、Sm、Gd为0.01×10−6外，其余均为0.003×10−6。

3 结果

3.1 岩相学

马山柱状节理火山岩的岩石薄片镜下分析显示：岩石主要由斑晶、基质组成。斑晶显示柱状，呈角闪石假象组成，绿鳞石化、少量绿泥石化等呈假象，粒度一般为0.2～0.8 mm，分布零散，少数聚斑状产出，部分具暗化边结构（图1d）。基质由长石、暗色矿物假象、石英组成，长石主要为斜长石，呈半自形微板条状，粒度一般＜0.1 mm，半定向或呈交织状分布，具交织结构，其间填隙暗色矿物（已蚀变为绿泥石等）及石英。杏仁体不规则状，大小0.2～1 mm零散分布，被石英、绿泥石等填充（图1d）。从岩石薄片分析即可发现石英含量较典型安山岩要高，薄片鉴定可将其归为英安岩类。

3.2 锆石U-Pb与Lu-Hf

马山柱状节理粗面英安岩中挑选获得的23颗锆石进行U-Pb测年（表1），选取协和度≥90%、结果显示最年轻一组12颗锆石U-Pb加权平均年龄为113.2 ± 1.3 Ma （MSWD=0.64，图2），对应早白垩世晚期Aptian与Albian期界线附近，其该组锆石的Th/U比值均大于0.4，表明锆石多为岩浆锆石[20]，可用于代表马山柱状节理粗面英安岩的形成时代。

图2 马山粗面英安岩锆石特征与测年结果a. 马山粗面英安岩的锆石阴极发光图像及对应的Pb206/U238年龄与εHf（t）值，b. 加权平均值，c. U-Pb年龄谐和图。Fig.2 Zircon characteristics and dating results of the trachyandesite in Mashana. cathodoluminescence （CL） images,yellow numbers=U-Pb ages and red number= εHf（t） values;b. weighted mean age,c. concordia diagrams of zircons from the Mashan trachydacite.

表1 马山粗面英安岩锆石U-Pb测试结果Table 1 Zircon U-Pb dating results of trachydacite at Mashan

最年轻一组锆石显示相似的左倾稀土配分特征（标准化值据McDonough和Sun[21]），重稀土富集，具δCe正异常及轻微的δEu负异常（图3）。碎屑锆石的稀土元素特征与形成时代关系密切，表明不同时代锆石的结晶条件存在差异（图3）。锆石原位Lu-Hf测试结果显示：所测锆石的εHf（t）值均为负值（−0.42～−19.48），对应的二阶段模式年龄均为古元古代—太古代（表2，图4）。

图4 马山粗面英安岩中锆石的εHf（t）值与对应的二阶段模式年龄Fig.4 εHf（t） values and two stage Hf model ages of zircons from trachydacite in Mashan

表2 马山粗面英安岩锆石原位Lu-Hf测试结果Table 2 Zircon in-situ Lu-Hf results of trachydacite at Mashan

图3 马山粗面英安岩中锆石的稀土元素配分曲线[21]Fig.3 Rare earth element distribution pattern of zircons from Mashan trachydacite[21]

3.3 主微量元素特征

马山柱状节理火山岩4件样品的主量元素含量如表3所示，SiO2含量较高，为65.44%～67.26%，平均为66.74%；Al2O3含量为13.37%～14.58%，平均为13.69%；K2O含量较高，为3.58%～3.98%，平均为3.79%；Na2O含量为3.77%～4.71%，平均为4.34%；全碱w（Na2O+K2O）含量为7.59%～8.24%，平均为8.13%；TiO2、MnO、P2O5含量均较低，平均含量分别为0.44%、0.09%、0.19%；MgO平均含量为2.14%；CaO平均含量为2.85%；全铁含量较高（3.61%～3.89%），平均含量为3.69%。岩石全碱TAS图解（据Le Bas[22]，图5）中，样品点均落于粗面英安岩区域，与镜下鉴定结果吻合。

图5 马山柱状节理粗面英安岩岩石类型划分TAS图解[22]Fig.5 TAS diagram of the trachydacite at Mashan[22]

表3 马山粗面英安岩全岩主量元素分析结果Table 3 Whole-rock major element compositions of trachydacite at Mashan %

4件马山粗面英安岩的微量元素（表4）特征为：Sr含量高（341×10−6～489×10−6），平均为413.5×10−6；Yb与Y含量低，均值分别为1.40×10−6、13.5×10−6；稀土元素总量较低（表5），ΣREE为119.2×10−6～131.1×10−6，平均为126.1×10−6，明显低于地壳岩浆岩平均值。稀土配分模式呈轻稀土显著富集、重稀土强烈亏损的右倾模式，（La/Yb）N为13.51～14.56，（Ce/Yb）N为9.31～9.98。δEu范围为0.77～0.93，平均值0.88，显示轻微负异常（图6）。此外，马山粗面英安岩明显富集大离子亲石元素Ba、Pb，亏损Th、Nb、Ti等高场强元素（图7）。

图6 马山粗面英安岩稀土元素配分曲线[21]Fig.6 REE distribution pattern of the trachydacite at Mashan[21]

图7 马山粗面英安岩微量元素蛛网图[21]Fig.7 Spider distribution pattern of the trachydacite at Mashan[21]

表4 马山粗面英安岩微量元素分析结果Table 4 Trace element compositions of trachydacite at Mashan 10−6

表5 马山粗面英安岩稀土元素分析结果Table 5 REE compositions of trachydacite at Mashan 10−6

4 讨论

4.1 马山粗面英安岩的埃达克岩属性、形成时代及区域对比

埃达克岩概念首先由Defant and Drummond于1990年在研究由年轻岩石圈俯冲熔融形成的现代弧岩浆时提出，系于现代弧岩浆由Kay[23]首先记录于阿留申群岛的埃达克岛（Adak Island），因之命名为Adakite[24]。埃达克岩初始定义用于描述由年轻的大洋岩石圈（≤25 Ma）俯冲形成的火山岩或侵入岩，具有SiO2含量≥56%，Al2O3≥15%（很少低于此值），Na2O＞3.5%，MgO含量通常＜3%（很少高于6%），相对岛弧性质的安山岩、英安岩和钠质流纹岩组合具有更低的Y和重稀土含量（Y≤18×10−6；Yb≤1.9×10−6），Sr含量高（很少＜400×10−6），高场强元素含量低等特征[24-25]。国内埃达克岩的研究，由张旗等在21世纪初引入并引发关于埃达克岩的广泛讨论[6-10,12,26-31]。张旗等将中国东部与板块消减无关、富K2O（Na2O/K2O=0.9～1.3）并具埃达克岩地化特征的陆相火成岩称之为“C型”埃达克岩，而将Defant 和Drummond[24]所定义的与年轻的大洋板块俯冲有关的埃达克岩称之为“O型”埃达克岩，并赋予“C”型埃达克岩明确的构造意义——代表增厚的地壳[6-8,31]。

本文所测4件马山粗面英安岩的主微量元素值，除Al2O3含量（平均值13.69%）略低外，各项指标均符合典型“C型”埃达克岩的上述特征[6-8,31]。Sr/Y-Y比值判别图解[24]与（La/Yb）N-YbN判别图解[32]的投点亦均落入典型埃达克岩范围内（图8）。4件样品的HREE与Y元素均表现为亏损，这可能是因为来源熔体的残留相中存在石榴石导致[24,33]。Eu的轻微负异常也表明来源岩浆源区很少或没有斜长石，对应的残留相为榴辉岩，表明岩浆很可能来自增厚地壳的熔融[6,8,12,33]。实验岩石学证据表明，Sr与Yb元素特征与压力密切相关[8,12]，马山粗面英安岩样品具有显著的高Sr（均值414×10−6）、低Yb（均值1.40×10−6）特征，可能代表了高压（增厚地壳）环境。此外，统计数据表明，对于安山岩类，在SiO2含量一定的情况下，K2O的含量越高对应来源岩浆所处的地壳厚度越大[34]，本文4件样品的K2O含量较高，Na2O/K2O平均值为1.15，亦暗示马山粗面英安岩的岩浆来源较深，可能形成于增厚地壳的熔融。4件马山粗面英安岩的Ce/Pb比值为0.97～1.73，均值1.31，与洋壳来源的Ce/Pb比值特征（14.38）截然不同[35]，表明来源岩浆源区不可能为洋壳。韩宗珠等[14]所测马山的“ms04”与“ms08”等2件样品、何登洋等[15]所测马山的“17MS01”、“17MS02”、“17MS03”、“17MS04”等4件样品的主微量元素特征与本文所获结果相似。综上，马山粗面英安岩为典型的“C型”埃达克岩，岩浆最可能来源于增厚地壳的熔融。

图8 马山粗面英安岩Sr/Y vs.Y[24]和(La/Yb)N-YbN[32]判别图解Fig.8 Sr/Y -Y and (La/Yb)N-YbN adakitic trace elemental discrimination diagrams for Mashan trachydacite[24,32]

年代学方面，何登洋等[15]近来借助LA-ICPMS锆石U-Pb定年，将马山粗面英安岩的形成时代定为119.3±1.6 Ma（n=9,MSWD=1.9）。但其用于确定喷出岩年龄的9颗锆石中年龄最小的2颗锆石年龄分别为112±1 Ma及115±1 Ma[13]，与本文所测的锆石U-Pb年龄在误差范围内一致，而迥然不同于其余7颗锆石年龄（约119～123 Ma）[15]。通过比较，可以推测马山粗面英安岩的实际形成时间可能更接近本文所测的113.2±1.3 Ma。

此外，马山地区该套早白垩世晚期的“C”型埃达克岩在区域上并非个例，鲁东与辽东地区亦存在大量同时期的埃达克岩，如：文登长山南花岗岩（114 Ma）和胶东的六度寺（115 Ma）、泰薄顶（114 Ma）、三佛山与辽东的古道岭（113 Ma）等[36]。

4.2 马山地区粗面英安岩来源岩浆形成的动力学机制初探

中国东部燕山期经历了由特提斯构造域/古亚洲洋构造域向太平洋构造域的转变[2,37-39]，晚中生代在东亚外侧形成壮观宏伟的燕山造山带的同时[2]，从俄罗斯楚科奇半岛到中国浙闽沿海发育了范围广大（长度超过5 000 km）的岩浆岩带[39]。对东亚大陆边缘晚中生代大规模岩浆活动的动力学机制一直存在较大分歧。胡受奚等[40]强调古太平洋板块向欧亚板块俯冲的弧后效应可能是主要原因[40]；王德滋等[41]根据晚侏罗世—早白垩世火山岩的地化特征，将长江断裂带以北、郯庐断裂带以东地区分布的火山岩归为橄榄粗安岩省，认为其生成于华北-扬子板块对接后的陆内造山与晚期伸展、古太平洋板块向欧亚板块的俯冲消减及与白垩纪后特提斯构造体制有关的碰撞-挤出作用，三大构造动力体制所导致的郯庐断裂带大规模左行平移营造出的张应力环境。董树文等[4-5]则强调燕山期东亚多块体的汇聚模式。万天丰等[42]提出圈层滑脱与高地温梯度的洋陆过渡型岩石圈沿构造断裂的活动可能是中国东部中生代岩浆起源的动力学机制。张旗[43]根据中国东部中生代缺少岛弧玄武岩和岛弧花岗岩及太平洋板块俯冲方向与大规模岩浆活动出现时代的不匹配等，认为中国东部中生代大规模岩浆活动与太平洋板块向西俯冲无关。

根据本文马山柱状节理粗面英安岩及胶辽地区广泛分布的早白垩世晚期“C型”埃达克岩，可以推测该时期马山地区粗面英安岩的岩浆来源于增厚地壳的熔融。但该时期熔融的增厚地壳是燕山运动早期的增厚地壳的残留（或范围缩减的高原）[31]还是由燕山运动多期次、间歇性幕式区域性挤压形成的增厚地壳？是来自华北板块下部地壳还是俯冲的扬子板块？从马山MS-01样品的锆石U-Pb与Lu-Hf测试结果分析，马山粗面英安岩来源岩浆中结晶锆石的二阶段模式年龄集中于古元古代—太古代（表2，图4）。古元古代—太古代的锆石组分在华北板块扬子北缘均广泛分布[44-45]。但根据本文及何登洋等[15]所测马山粗面英安岩样品的锆石UPb年龄分析，锆石中分别出现了代表扬子板块与苏鲁造山带的新元古代与晚石炭世组分，表明增厚地壳的熔融物质中含有扬子陆壳的组分[15]。此外，从马山粗面英安岩中具早白垩世晚期年龄锆石的Hf同位素特征看，并未呈现出在拉张环境下大陆弧地区产出锆石的Hf同位素受亏损地幔影响而导致εHf（t）值更向正值的趋势[1]。该结果亦暗示马山柱状节理粗面英安岩的岩浆来源更可能来源于区域挤压导致增厚地壳的熔融而非受控于太平洋板块俯冲的弧后效应。这与马山粗面英安岩的主微量元素特征所反映出来的“C型”埃达克岩所代表的构造含义一致。

从大区域上看，东亚在早白垩世发生过数次幕式的挤压事件[4-5,46-48]。与马山柱状节理粗面英安岩形成时代接近的挤压事件对应于Aptian晚期，Guo等[48]将该期区域性挤压事件的动力学机制解释为：华南板块与其东南部——由现在的民都洛岛、朗布隆岛、班乃岛、巴拉望岛、婆罗洲东部及苏拉威西北部等构成的微陆块发生碰撞所致。即墨马山相邻的灵山岛地区在早白垩世Aptian期也存在NWSE向挤压事件的明确记录，证据如：千层崖剖面西侧大型褶皱、平面X剪破裂所反映的古应力场特征等[49]。胶莱盆地西侧的临沂方城盆地八亩地组与上覆火山岩间存在沉积间断，甚至发育有古风化壳[50]。在燕山-辽西地区，该期碰撞事件导致早白垩世晚期的九佛堂组与阜新组之间形成区域性不整合[48,51-52]，向内陆西北方向甚至可影响至二连盆地，二连盆地内早白垩世晚期的腾格尔组与赛罕组之间存在不整合[48]等。

马山地区所处的胶莱盆地，目前绝大多数研究者认为其为拉张性盆地。但姜同海[53]根据①盆地发育于基底褶皱之上，盆地南部为一逆冲断裂；②白垩系下统的沉积速率与挤压型盆地一致（由慢变快）；③盆地边缘的莱阳群林寺山组底部发育磨拉石沉积，向内逐渐变为复理石沉积等特征，认为胶莱盆地属于挤压挠曲型盆地。特别是，早白垩世中晚期该地区构造应力场为NW-SE向挤压，与周缘地区记录的挤压事件吻合。如：Wu 等[54]利用高（锆石U-Pb）-中（白云母Ar-Ar）-低（锆石U-Th/He）3种温标对胶东地区早白垩世中晚期两处侵入岩体的热史进行约束，发现在117～110 Ma存在一期快速冷却事件。该期热事件可与区域上该期挤压隆升事件相对应。此外，Zhang 等[55]曾借助碳酸盐团簇同位素估算胶莱盆地晚白垩世期间的古海拔超过2 000 m，在约90 Ma之前甚至更高。但该时期高原或海岸山脉的形成是受控于鄂霍茨克板块与东亚大陆边缘的碰撞[56]还是从早白垩世晚期持续至晚白垩世早期的山脉或高原[57]，还需更多证据与线索。

综合上述马山早白垩世柱状节理粗面英安岩的“C型”埃达克岩属性、区域上同时期的构造、沉积响应及锆石U-Pb与Hf同位素特征，笔者更倾向于认为：早白垩世Aptian晚期，马山地区经历区域性挤压事件造成地壳增厚，进而导致原俯冲于华北板块之下的扬子板块部分熔融并沿同时期切割深度较大的断裂通道（如：牟平-即墨断裂）喷发形成范围广阔的火山岩盆地。

5 结论

（1）锆石U-Pb定年结果显示马山地区的柱状节理粗面英安岩的形成时代为早白垩世113.2 Ma，岩石全岩主微量元素分析结果显示其具有“C型”埃达克岩特征。

（2）马山粗面英安岩中所获锆石的二阶段模式年龄均为古元古代—太古代，暗示来源岩浆形成于古老基底的熔融。

（3）马山柱状节理粗面英安岩类的岩浆可能源于燕山期区域性NW-SE向挤压导致的大陆地壳增厚后的熔融，熔融物质包含来自俯冲于华北板块之下的扬子板块组分。