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

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

朱晓青, 侯方辉, 刘洪滨, 郭兴伟, 孙天本, 秦亚超, 安郁辉, 李凤春. 山东即墨马山粗面英安岩年代学与地球化学特征及其地质意义[J]. 海洋地质与第四纪地质, 2021, 41(6): 138-150. DOI: 10.16562/j.cnki.0256-1492.2021011801
引用本文: 朱晓青, 侯方辉, 刘洪滨, 郭兴伟, 孙天本, 秦亚超, 安郁辉, 李凤春. 山东即墨马山粗面英安岩年代学与地球化学特征及其地质意义[J]. 海洋地质与第四纪地质, 2021, 41(6): 138-150. DOI: 10.16562/j.cnki.0256-1492.2021011801
shu
ZHU Xiaoqing, HOU Fanghui, LIU Hongbin, GUO Xingwei, SUN Tianben, QIN Yachao, AN Yuhui, LI Fengchun. Geochronology and geochemistry of Mashan trachydacite, JiMo District, Shandong Province and their geological implications[J]. Marine Geology & Quaternary Geology, 2021, 41(6): 138-150. DOI: 10.16562/j.cnki.0256-1492.2021011801
Citation: ZHU Xiaoqing, HOU Fanghui, LIU Hongbin, GUO Xingwei, SUN Tianben, QIN Yachao, AN Yuhui, LI Fengchun. Geochronology and geochemistry of Mashan trachydacite, JiMo District, Shandong Province and their geological implications[J]. Marine Geology & Quaternary Geology, 2021, 41(6): 138-150. DOI: 10.16562/j.cnki.0256-1492.2021011801
shu

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

  • 1.

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

  • 2.

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

  • 3.

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

  • 4.

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

  • 5.

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

基金项目: 中国地质调查局海洋地质调查项目（DD20190365，DD20190377，DD20190205）
详细信息
    作者简介:

    朱晓青（1982—），男，博士，副研究员，主要从事构造地质与海洋地质研究，E-mail：xq_zhu0819@126.com

    通讯作者:

    侯方辉（1979—），男，博士，高级工程师，主要从事构造地质及海洋地质研究，E-mail：houfanghui3993@163.com

  • 中图分类号: P542

Geochronology and geochemistry of Mashan trachydacite, JiMo District, Shandong Province and their geological implications

  • 1.

    Qingdao Institute of Marine Geology, Qingdao 266071, China

  • 2.

    Laboratory for Marine Mineral Resources, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China

  • 3.

    Institute of Marine Development of Ocean University of China, Qingdao 266100, China

  • 4.

    Management Office of Mashan Nature Reserve, Jimo District, Qingdao 266200, China

  • 5.

    Testing Center of Shandong Bureau, China Metallurgy and Geology Bureau, Jinan 250014, China.

摘要

    摘要
  • 摘要: 山东即墨马山地区出露罕见的发育柱状节理的中酸性火山岩，对于其年代学及形成的大地构造背景鲜有讨论。岩石薄片与主量元素分析结果证明马山柱状节理火山岩属于粗面英安岩。利用LA-ICP-MS对马山柱状节理火山岩中的锆石进行U-Pb测年，确定其形成时间为113.2±1.3 Ma，属于早白垩世晚期。原位Lu-Hf分析获得锆石εHf(t)值范围为−19.5～−0.4，对应的二阶段模式年龄范围约为1.97～3.26 Ga，暗示其源区岩浆成生于古老地壳的熔融。马山粗面英安岩具有SiO2和Al2O3含量较高，MgO含量低及富含碱的主量元素的特征。其稀土元素总量较低，配分模式呈轻稀土富集、重稀土亏损的右倾模式，Eu轻微负异常，Sr含量高而Y和Yb含量低。马山粗面英安岩的主微量元素显示其具有埃达克岩特征，推测其来源岩浆可能形成于燕山期华北东部区域性地壳增厚后下部老地壳的熔融。
    Abstract: Intermediate-acid volcanic rocks with columnar joints expose in the Mashan area, Jimo, Shandong Province. The geochronology and tectonic setting of the volcanic rocks had been rarely studied. The petrographic characteristics and major element results show that the volcanic rocks at Mashan are trachydacite. Zircon LA-ICP-MS U-Pb dating of the trachydacite shows that the volcanic rocks erupted at 113.2 ± 1.3 Ma, that is, in the late Early Cretaceous. In-situ Lu-Hf analysis suggests that the zircon εHf (t) values range from −19.5 to −0.4 and two-stage model ages of zircons range from 1.97 Ga to 3.26 Ga, indicating that the magma was sourced from the melting of ancient crust. The Mashan trachyandesite is characterized by high content of SiO2 and Al2O3, low content of MgO, enriched Na2O and K2O. Moreover, all samples are featured by lower total contents of rare earth elements, enriched light rare earth elements, deficit heavy rare earth elements, slightly negative Eu anomaly, high Sr content, low Y and Yb contents. The major and trace elements of the trachydacite in Mashan are quite similar to adakite, and the adakitic magma could have been formed by the melting of the lower old crust after the regional crust thickening in Eastern Norch China Block during the Yanshanian.
    • HTML全文

  • 海洋沉积物是地球上最大的甲烷储库[1]。但因大部分甲烷在海洋环境中发生厌氧氧化作用（Anaerobic Oxidation of Methane，AOM），从而导致大气中只有不到2%的甲烷是由海洋排放而来[2]。通常，海洋沉积物中产生的甲烷在向上扩散过程中，与浅层沉积物孔隙水中向下扩散的硫酸盐相遇后，在甲烷氧化古菌和硫酸盐还原菌共同作用下发生AOM，导致硫酸盐和甲烷同时消耗，呈现出一个相反的垂直浓度变化剖面，并在沉积物层中形成硫酸盐—甲烷转换带（Sulfate Methane Transition Zone，SMTZ）[3-5]。

    以往研究大多认为，在SMTZ内仅发生AOM反应，而近期研究发现，海底沉积层中向下扩散进入SMTZ的硫酸盐通量大于向上扩散的甲烷通量[6-7]，表明SMTZ内部发生了有机质和甲烷共同竞争消耗硫酸盐的过程[8-9]。如在圣巴巴拉海盆就已发现SMTZ内部有35%～45%的SO42−是通过有机质的硫酸盐还原作用（Organoclastic Sulfate Reduction，OSR）所消耗[9]。另外，Yoshinaga等[10]研究发现，在SMTZ底部经常出现13CH4亏损，是因为SMTZ内还“隐藏”着某种产甲烷反应，如“AOM生成的DIC”二氧化碳还原产甲烷（Carbon Dioxide Reduction，CR）、有机质降解产甲烷（Methaogenesis from Organic matter degradation，ME）和AOM逆反应产甲烷（AOM back flux）等[11-13]。正是由于海洋沉积物中SMTZ内部存在AOM、OSR和“隐藏的”产甲烷反应等生物地球化学过程的交叉耦合，使其碳转化与碳循环过程十分复杂。Hong等[14]构建了一个基于化学计量学、扩散通量和同位素质量平衡的“箱式模型”，并成功利用孔隙水地球化学参数定量估算了AOM、OSR、CR和ME等反应过程对郁陵海盆沉积物中SMTZ内碳循环的相对贡献。

    杭州湾海底沉积物中蕴藏着大量的浅层生物气（主要成分是CH4），水深范围约0～25 m，作为一种海底天然气资源或一种海洋地质灾害因素已得到广泛关注[15]。但作为温室气体CH4的重要载体，海底浅层生物气又是一种潜在的气候致变因子，其释放、迁移和扩散行为可能引发一系列的负面环境效应，如温室效应的加剧、海洋生态环境恶化等。因此，研究杭州湾海底沉积物中的甲烷厌氧氧化以及相关的碳循环过程，对正确评估浅层生物气的生态环境效应具有重要的理论意义。

    本文对杭州湾YS6孔柱状沉积物样品进行了孔隙水部分阴阳离子（SO42−、Ca2+和Mg2+）、DIC浓度及δ13CDIC、CH4含量及δ13CCH4等地球化学分析，在此基础上利用“箱式模型”对硫酸盐—甲烷转换带内部CH4和DIC的碳转化与碳循环过程进行定量研究，以期为深刻认识杭州湾沉积物中甲烷的迁移与转化机制提供科学依据。

    1.   地质背景

    杭州湾位于浙江省东北部，西部湾顶连接钱塘江河口，湾口北部与长江口南翼毗邻，东、南部通过舟山群岛之间的潮汐水道与东海相通，是一个典型的喇叭型强潮河口湾，湾口宽约95 km，自口外向口内逐渐狭窄[16-17]。采样站位YS 6孔位于杭州湾内舟山群岛近岸海域（图1），水深约16 m，是长江、钱塘江等河流携带陆源碎屑物质入海的主要沉积场所，沉积有机质受陆源和海源物质输入综合影响[18]。研究区在第四纪的几次海进、海退过程中，交替沉积了数套富含有机质的淤泥层和砂层。淤泥层中的有机质经厌氧菌的生物化学作用产生生物气（主要成分是CH4），并呈垂向多层化发育，分布于海底百米深度以内的全新统和上更新统地层中[19-20]。第四纪以来，钱塘江沿海平原河口湾相沉积速率为4～5 mm/a，且长期处于稳定沉降过程中，这种快速沉积、持续沉降作用导致沉积物中浅层生物气充足[19]。同时，在地质历史时期的多次海进—海退旋回过程使得研究区沉积物孔隙水富含一定浓度的SO42−，为沉积物中生物气发生硫酸盐还原型AOM提供了物质基础。

    图  1  YS 6站位位置
    Figure  1.  Map showing the locations of Core YS 6 in Hangzhou Bay
    下载: 全尺寸图片 幻灯片

    2.   样品与方法

    YS6站位沉积物柱状样于2017年5月利用旋转钻井技术钻取获得，位于30°19'32.54"N、121°54'03.39"E，岩心柱总长60.6 m，直径108 mm，本文仅研究0～30 mbsf （meter below seafloor）层段。沉积物岩心钻获后，在甲板上现场进行分割、观察与描述，然后按2 m间隔进行沉积物孔隙水和顶空气等样品的采集，但18～30 mbsf之间未采集孔隙水样品。所有分析测试工作均在中国地质调查局青岛海洋地质研究所实验测试室完成。

    孔隙水样品利用Rhizon采样器进行采集。将其多孔聚合物管端插入沉积物中，另一端与抽成真空的、干净的20 mL医用注射器相连，大约1～2 h后，每个样品获得10～19 mL的孔隙水。随后，将获得的孔隙水样品分成两份，分别保存在两个10 mL的玻璃瓶中，并用氮气驱除玻璃瓶内顶部空间的空气。对其中一份水样注入10 μL饱和HgCl2溶液以消除微生物活性，用于溶解无机碳（Dissolved inorganic carbon，DIC）含量及其碳同位素分析。分析仪器为美国Thermo Fisher公司生产的MAT253型GasBench II-IRMS仪，分析方法参照杨涛等[21-22]，分析精度分别小于0.5%和0.10‰。对另一份注入0.1 mL HNO3（8 mol/L）以防止氧化还原和沉淀反应发生，用于阴、阳离子浓度分析。其中，SO42-浓度利用美国戴安公司生产的ICS-3000型离子色谱仪进行测定[23]，分析精度小于2%；Ca2+和Mg2+浓度利用美国Thermo Fisher 公司生产的IRIS INTREPID II XSP型全谱直读等离子体发射光谱仪（ICP-AES）进行测定[24]，分析精度小于3%。

    为最大限度地降低沉积物岩心在长时间放置、取样过程中的气体逸散损失，沉积物岩心钻取至甲板上以后，现场迅速对沉积物顶空气样品进行采集。顶空气样品取自孔隙水相同或相邻的沉积物岩心层位，利用头部被切割去除的一次性无菌医用注射器插取35 mL沉积物鲜样注入50 mL玻璃顶空瓶内，迅速加入10 mL饱和NaCl溶液后用丁基橡胶塞和铝盖密封，倒置、避光保存。顶空气样品在采样工作完成后，于一周内运抵实验室并完成气体成分及其碳、氢同位素组成的分析测试。分析仪器采用美国Thermo Fisher公司生产的Ultra Trace型气相色谱仪（GC-TCD）和MAT253型同位素比值质谱仪（GC-GC Isolink-IRMS），分析方法参照贺行良等[25-26]，测定精度分别小于3%和0.20‰。

    3.   碳循环模型与计算

    3.1   箱式模型

    沉积物中SMTZ内碳循环过程的研究采用“箱式模型”[14]，该模型包括5个组分，涉及2个固体组分（有机碳、碳酸盐岩）和5个溶解组分（SO42−、CH4、DIC、Ca2+、Mg2+）。基于在输入和输出速率、扩散通量之间的质量守恒，7个组分在SMTZ中参与碳循环的路径如图2所示。该模型的应用基于两个基本假设：第一，除AOM、OSR、CR、ME、以及碳酸盐矿物沉淀（Carbonate Precipitation，CP）作用外，其他反应对孔隙水中DIC的影响均忽略不计；第二，孔隙水Ca2+、Mg2+浓度剖面的变化完全是碳酸盐矿物沉淀与溶解的结果。“箱式模型”中所考虑的5个反应过程如下[14]：

    图  2  SMTZ内碳循环“箱式模型”图[14]
    溶解组分用黄色框表示，固体组分用灰色框表示；灰线箭头表示模型域之外的DIC、SO42−和CH4的输入/输出；5个主要反应路径AOM、OSR、CR、ME和CP分别用不同颜色的箭头标识。
    Figure  2.  Illustration of box model for SMTZ-internal carbon cycling[14]
    The dissolved components are shown in yellow and the solid components in gray; The gray arrow represents the input/output of DIC, SO42− and CH4 outside the box model; The reaction paths of AOM, OSR, CR, ME and CP are marked by different colored arrows.
    下载: 全尺寸图片 幻灯片

    （1）有机质的硫酸盐还原，OSR

    $$ {\rm{C}}{{\rm{H}}_2}{\rm{O}} + 0.5{\rm{SO}}_4^{2 - } \to {\rm{HCO}}_3^ - + 0.5{{\rm{H}}_2}{\rm{S}} $$ (1)

    （2）碳酸盐沉淀，CP

    $$ {\rm{C}}{{\rm{a}}^{2 + }} + {\rm{CO}}_3^{2 - } \to {\rm{CaC}}{{\rm{O}}_3};\;\;{\rm{M}}{{\rm{g}}^{2 + }} + {\rm{CO}}_3^{2 - } \to {\rm{MgC}}{{\rm{O}}_3} $$ (2)

    （3）甲烷厌氧氧化，AOM

    $$ {\rm{C}}{{\rm{H}}_4} + {\rm{SO}}_4^{2 - } \to {\rm{HCO}}_3^ - + {\rm{H}}{{\rm{S}}^ - } + {{\rm{H}}_2}{\rm{O}} $$ (3)

    （4）“AOM生成的DIC”产甲烷，CR

    $$ {\rm{C}}{{\rm{O}}_2} + 4{{\rm{H}}_2} \to {\rm{C}}{{\rm{H}}_4} + 2{{\rm{H}}_2}{\rm{O}} $$ (4)

    （5）有机质降解产甲烷，ME

    $$ {\rm{C}}{{\rm{H}}_2}{\rm{O}} + {{\rm{H}}_2}{\rm{O}} \to 0.5{\rm{C}}{{\rm{H}}_4} + 0.5{\rm{C}}{{\rm{O}}_2} + {{\rm{H}}^ + } $$ (5)

    需要说明的是，在箱式模型中OSR反应（式1）消耗了部分来自底层海水的SO42−产生DIC；AOM反应（式3）在SMTZ内同时消耗CH4和SO42−产生DIC，继而引发CP反应发生，消耗了大量来自底层海水的Ca2+和Mg2+（式2）。ME反应（式5）产生的DIC直接来自于有机物的降解，而CR反应（式4）利用的DIC则由AOM产生，与有机物的降解没有直接联系。FCH4.in表示自深部含甲烷沉积层扩散而来的CH4通量，与模型框架中考虑的任何反应无关。

    3.2   质量平衡

    考虑在SMTZ以上沉积层中CH4含量、以及SMTZ以下沉积层中SO42−含量均低于检出限，故实际计算中将假设SO42−和CH4在穿过SMTZ后均被完全消耗，即FSO4.out 和FCH4.out可以忽略。因此，根据图2碳循环框架图，在SMTZ范围内，SO42−、DIC和CH4的质量平衡关系可用以下3个方程式（式6—式8）来表述：

    $$ {{\rm{F}}_{{\rm{SO}}4.{\rm{in}}}} = {{\rm{R}}_{{\rm{OSR}} - {\rm{S}}}} + {{\rm{R}}_{{\rm{AOM}}}} $$ (6)
    $$ {{\rm{F}}_{{\rm{DIC}}.{\rm{in}}}} + {{\rm{R}}_{{\rm{ME}} - {\rm{DIC}}}} + {{\rm{R}}_{{\rm{OSR}} - {\rm{C}}}} + {{\rm{R}}_{{\rm{AOM}}}} = {{\rm{F}}_{{\rm{DIC}}.{\rm{out}}}} + {{\rm{R}}_{{\rm{CR}}}} + {{\rm{R}}_{{\rm{CP}}}} $$ (7)
    $$ {{\rm{F}}_{{\rm{CH}}4.{\rm{in}}}} + {{\rm{R}}_{{\rm{CR}}}} + {{\rm{R}}_{{\rm{ME}} - {\rm{CH}}4}} = {{\rm{R}}_{{\rm{AOM}}}} $$ (8)

    式中：“F”表示扩散通量，“R”表示反应速率，“in”和“out”指示SO42−、DIC和CH4扩散进出SMTZ的相对传输方向；RAOM为甲烷厌氧氧化反应速率、RCR为“AOM生成的DIC”产甲烷反应速率、RCP为碳酸盐沉淀反应速率、RME-DIC为有机质降解产甲烷作用中DIC的生成速率、ROSR-S和ROSR-C分别为OSR作用中硫酸盐的消耗速率和DIC的生成速率，根据反应式（1）中的化学计量关系可知：

    $$ {{\rm{R}}_{{\rm{OSR}} - {\rm{C}}}} = 2 \times {{\rm{R}}_{{\rm{OSR}} - {\rm{S}}}} $$ (9)

    假定孔隙水溶解组分的浓度梯度变化属于稳定状态，则孔隙水中经扩散通过SMTZ的5个溶解组分（Ca2+、Mg2+、DIC、CH4和SO42−）的扩散通量可根据Fick第一定律来计算[27-28]，计算公式如下：

    $${F_X} = - \varphi \cdot {D_X}\frac{{\partial {C_X}}}{{\partial Z}}$$ (10)

    式中，FX为扩散速率（mmol·cm−2·d−1），φ为沉积物孔隙度（平均约0.528，另文发表），DX为某组分X在沉积物中的扩散系数（cm2·d−1），CX为某组分X的浓度（mmol·L−1），Z为沉积物深度（m），$\dfrac{{\partial {C_X}}}{{\partial Z}}$为某组分X在沉积物中随深度变化的浓度梯度（mmol·L−1·m−1），根据浓度随深度变化的线性回归计算。扩散系数（DX）根据以下公式计算[28]：

    $${D_X} = \frac{{{D_0}}}{{1 - \ln({\varphi ^2})}}$$ (11)

    式中，D0为某组分在海水中的扩散系数，根据研究区海水温度取约20 ℃时的值[29]。

    扩散进出SMTZ内Ca2+和Mg2+的净通量用RCP表示，可根据公式（12）进行计算[30-31]：

    $$ {{\rm{R}}_{{\rm{CP}}}} = {{\rm{F}}_{{\rm{Ca}}}} + {{\rm{F}}_{{\rm{Mg}}}} = {{\rm{F}}_{{\rm{Ca}}{\rm{.in}}}} + {{\rm{F}}_{{\rm{Mg}}{\rm{.in}}}} - {{\rm{F}}_{{\rm{Ca}}{\rm{.out}}}} - {{\rm{F}}_{{\rm{Mg}}{\rm{.out}}}} $$ (12)

    3.3   碳同位素质量平衡

    各组分在扩散通过SMTZ时，除保持质量平衡外，轻、重碳同位素（12C和13C）质量平衡关系依然成立。根据式（7）和（8）有以下同位素质量平衡方程成立：

    $$\begin{split} &^{12}{{\rm{F}}_{{\rm{DIC}}{\rm{.in}}}}{ + ^{12}}{{\rm{R}}_{{\rm{OSR}} - {\rm{C}}}}{ + ^{12}}{{\rm{R}}_{{\rm{ME}} - {\rm{DIC}}}}{ + ^{12}}{{\rm{R}}_{{\rm{AOM}}}}=\\&{ ^{12}}{{\rm{F}}_{{\rm{DIC}}{\rm{.out}}}}{ + ^{12}}{{\rm{R}}_{{\rm{CR}}}}{ + ^{12}}{{\rm{R}}_{{\rm{CP}}}} \end{split}$$ (13)
    $$\begin{split} &^{13}{{\rm{F}}_{{\rm{DIC}}{\rm{.in }}}}{ + ^{13}}{{\rm{R}}_{{\rm{OSR }} - {\rm{C}}}}{ + ^{12}}{{\rm{R}}_{{\rm{ME}} - {\rm{DIC}}}}{ + ^{13}}{{\rm{R}}_{{\rm{AOM}}}}=\\&{ ^{13}}{{\rm{F}}_{{\rm{DIC}}{\rm{.out }}}}{ + ^{13}}{{\rm{R}}_{{\rm{CR}}}}{ + ^{13}}{{\rm{R}}_{{\rm{CP}}}} \end{split}$$ (14)
    $$ ^{12}{{\rm{F}}_{{\rm{CH}}4.{\rm{in}}}}{ + ^{12}}{{\rm{R}}_{{\rm{CR}}}}{ + ^{12}}{{\rm{R}}_{{\rm{ME}} - {\rm{CH}}4}} = 12{{\rm{R}}_{{\rm{AOM}}}} $$ (15)
    $$ ^{13}{{\rm{F}}_{{\rm{CH}}4.{\rm{in}}}}{ + ^{13}}{{\rm{R}}_{{\rm{CR}}}}{ + ^{13}}{{\rm{R}}_{{\rm{ME}} - {\rm{CH}}4}}{ = ^{13}}{{\rm{R}}_{{\rm{AOM}}}} $$ (16)

    另外，基于碳同位素质量平衡的考虑，结合测量得到的CH4和DIC等组分的碳同位素比值，可分别估算CH4和DIC的重碳（13C）和轻碳（12C）的扩散传输量。例如，可按方程式（17）来计算12CDIC和13CDIC的浓度分布：

    $$\begin{split} &{}^{12}{C_{DIC}} = \frac{{{C_T}}}{{\left( {1 + {r_{DIC}}} \right)}}\\ &{}^{13}{C_{DIC}} = {C_T} - {}^{12}{C_{DIC}} \end{split}$$ (17)

    式中：CT表示DIC的总浓度，rDIC是通过同位素值计算的DIC重碳（13C）和轻碳（12C）的浓度比值（13C/12C），根据公式（18）计算：

    $${r_{DIC}} = \left( {\frac{{\delta {}^{13}{C_{DIC}}}}{{1000}} + 1} \right) \times {r_{{\rm{std}}}}$$ (18)

    式中：rstd是碳同位素基准标准物质VPDB的重碳（13C）和轻碳（12C）的比值（13C/12C），取rstd = 0.011 237 2。其他组分的重碳（13C）和轻碳（12C）的浓度比值（如rCH4、rOM等）根据式（18）类推。

    “箱式模型”中的5个反应式（式1—式5）中的碳同位素分馏情况是理解SMTZ内碳循环过程的基础。5个组分之间通过包含不同程度碳同位素分馏的5个化学反应相关联，其中反应物和产物之间的碳同位素分馏效应可按公式（19）计算：

    $$\alpha = \frac{{{r_r}}}{{{r_p}}}$$ (19)

    式中，rr和rp分别为化学反应式中反应物和产物的碳同位素比率（13C/12C）。

    若5种化学反应按稳态一级反应考虑，则同位素分馏效应（α）还可用轻同位素的速率常数（kL）与重同位素的速率常数（kH）的比值来表示[32]：

    $${\rm{\alpha }} = \frac{{{{\rm{k}}^{\rm{L}}}}}{{{{\rm{k}}^{\rm{H}}}}}$$ (20)

    因此，各组分轻、重同位素的反应速率，将根据式（19）和式（20）计算如下[14]：

    $${}^{{\rm{12}}}{{\rm{R}}_{{\rm{OSR - S}}}} = \frac{1}{2}f \cdot {}^{{\rm{12}}}{{\rm{R}}_{{\rm{OM}}}}$$ (21)
    $$\begin{split} {}^{{\rm{12}}}{{\rm{R}}_{{\rm{OSR - C}}}} =& f \cdot {}^{{\rm{12}}}{{\rm{R}}_{{\rm{OM}}}} \\ {}^{{\rm{13}}}{{\rm{R}}_{{\rm{OSR - C}}}} =& {}^{{\rm{12}}}{{\rm{R}}_{{\rm{OSR - C}}}} \cdot {r_{OM}} \end{split} $$ (22)
    $$ \begin{split} & {}^{{\rm{12}}}{{\rm{R}}_{{\rm{AOM}}}} = {{\rm{F}}_{{\rm{SO4}}{\rm{.in}}}}{\rm{ - }}{}^{{\rm{12}}}{{\rm{R}}_{{\rm{OSR - S}}}} \\ &{}^{{\rm{13}}}{{\rm{R}}_{{\rm{AOM}}}} = {}^{{\rm{12}}}{{\rm{R}}_{{\rm{AOM}}}} \times \frac{{{r_{{\rm{CH4 - SMTZ}}}}}}{{{\alpha _{{\rm{AOM}}}}}} \\ \end{split} $$ (23)
    $$\begin{split} &{}^{{\rm{12}}}{{\rm{R}}_{{\rm{ME - CH4}}}} = {}^{{\rm{12}}}{{\rm{R}}_{{\rm{ME - DIC}}}} = \frac{{\rm{1}}}{{\rm{2}}}{}^{{\rm{12}}}{{\rm{R}}_{{\rm{ME}}}} = \frac{{\rm{1}}}{{\rm{2}}}(1{\rm{ - }}f) \cdot {}^{{\rm{12}}}{{\rm{R}}_{{\rm{OM}}}} \\ &{}^{{\rm{13}}}{{\rm{R}}_{{\rm{ME - CH4}}}} = {}^{{\rm{12}}}{{\rm{R}}_{{\rm{ME - CH4}}}} \cdot {r_{CH4 - bot}} \\ &{}^{{\rm{13}}}{{\rm{R}}_{{\rm{ME - DIC}}}} = {}^{{\rm{12}}}{{\rm{R}}_{{\rm{ME - DIC}}}} \cdot {r_{DIC - bot}} \\ \end{split} $$ (24)
    $$ \begin{split} & {}^{{\rm{12}}}{{\rm{R}}_{{\rm{CP}}}} = {F_{{\rm{Ca/Mg}}}}{\rm{/}}({\rm{1}} + {r_{DIC - SMTZ}}) \\ &{}^{{\rm{13}}}{{\rm{R}}_{{\rm{CP}}}} = {F_{{\rm{Ca/Mg}}}}{\rm{ - }}{}^{{\rm{12}}}{{\rm{R}}_{{\rm{CP}}}} \\ \end{split} $$ (25)

    式中，ROM表示有机质通过OSR和ME作用的降解速率：

    $$ {{\rm{R}}_{{\rm{OM}}}} = {{\rm{R}}_{{\rm{OSR}} - {\rm{C}}}} + {{\rm{R}}_{{\rm{ME}}}} $$ (26)

    f表示有机质通过OSR作用降解的质量分数；αAOM为AOM过程中的同位素分馏系数；r为有机质、DIC、CH4等组分在SMTZ处的13C与12C的比值，根据式（18）计算。

    结合式（13）与12ROM = 12ROSR-C + 12RME 可推导出：

    $$\begin{split} &{}^{{\rm{12}}}{{\rm{R}}_{{\rm{OM}}}} = \frac{{{}^{{\rm{12}}}{{\rm{F}}_{{\rm{DIC}}{\rm{.out}}}}{\rm{ - }}{}^{{\rm{12}}}{{\rm{F}}_{{\rm{DIC}}{\rm{.in}}}} + {}^{{\rm{12}}}{{\rm{R}}_{{\rm{CR}}}} + {}^{{\rm{12}}}{{\rm{R}}_{{\rm{CP}}}}{\rm{ - }}{}^{{\rm{12}}}{{\rm{R}}_{{\rm{AOM}}}}}}{{f + 0.5(1 - f)}} \\ &{}^{{\rm{13}}}{{\rm{R}}_{{\rm{OM}}}} = {}^{{\rm{12}}}{{\rm{R}}_{{\rm{OM}}}} \times {r_{OM}} \\[-15pt] \end{split} $$ (27)

    同时，为描述CR过程利用AOM产生的DIC生成甲烷（式4）的反应速率，假设在SMTZ内CR的速率与AOM的速率成正比，比例系数为b，则有方程式（28）成立：

    $$ \begin{split} &{}^{{\rm{12}}}{{\rm{R}}_{{\rm{CR}}}} = {\rm{b}}{}^{{\rm{12}}}{{\rm{R}}_{{\rm{AOM}}}}{\rm{ }} \\ &{}^{{\rm{13}}}{{\rm{R}}_{{\rm{CR}}}} = {\rm{b}}{}^{{\rm{13}}}{{\rm{R}}_{{\rm{AOM}}}} \times \frac{{\rm{1}}}{{{\alpha _{{\rm{CR}}}}}} \end{split} $$ (28)

    式中，b为比例系数，为0～1；αCR表示CR作用的同位素分馏系数，可根据沉积物岩心底层（此层碳同位素维持在一个较为固定的值，YS6孔约为20～30 mbsf层段，见图3）的DIC和CH4同位素值估算，即：

    图  3  YS6孔孔隙水部分溶解组分的垂直分布剖面
    阴影为SMTZ示意层；0 mbsf处为底层海水数据。
    Figure  3.  Pore water profiles of Core YS6
    Yellow shaded layer for SMTZ; Green shaded layer for MEZ; the data of 0 mbsf from bottom water.
    下载: 全尺寸图片 幻灯片
    $$ {{\rm{\alpha }}_{{\rm{CR}}}} = 1 + \left( {{{\rm{\delta }}^{13}}{{\rm{C}}_{{\rm{DIC}} - {\rm{bo}}{{\rm{t}}^ - }}}{{\rm{\delta }}^{13}}{{\rm{C}}_{{\rm{CH}}4 - {\rm{bot}}}}} \right) \times 1000 $$ (29)

    通常，在SMTZ处的甲烷浓度（13CCH4-SMTZ和12CCH4-SMTZ）接近于零，可视为远低于底部甲烷浓度（13CCH4-bot和12CCH4-bot）。因此，可以将轻、重甲烷的通量比（13FCH4.in/12FCH4.in）近似等同于底层轻、重甲烷的浓度比13CCH4-bot/12CCH4-bot，故根据方程式（15）和（16），可得以下关系式：

    $$\begin{split} \frac{{{}^{{\rm{13}}}{{\rm{F}}_{{\rm{CH4}}{\rm{.in}}}}}}{{{}^{{\rm{12}}}{{\rm{F}}_{{\rm{CH4}}{\rm{.in}}}}}} =& \frac{{\dfrac{{{}^{{\rm{13}}}{{\rm{C}}_{{\rm{CH4 - bot}}}} - {}^{{\rm{13}}}{{\rm{C}}_{{\rm{CH4 - SMTZ}}}}}}{L}}}{{\dfrac{{{}^{{\rm{12}}}{{\rm{C}}_{{\rm{CH4 - bot}}}} - {}^{{\rm{12}}}{{\rm{C}}_{{\rm{CH4 - SMTZ}}}}}}{L}}} \approx \frac{{{}^{{\rm{13}}}{{\rm{C}}_{{\rm{CH4 - bot}}}}}}{{{}^{{\rm{12}}}{{\rm{C}}_{{\rm{CH4 - bot}}}}}} \\ =& {r_{{\rm{CH4 - bot}}}} = \frac{{{}^{{\rm{13}}}{{\rm{R}}_{{\rm{AOM}}}} - {}^{{\rm{13}}}{{\rm{R}}_{{\rm{CR}}}} - {}^{{\rm{13}}}{{\rm{R}}_{{\rm{ME - CH4}}}}}}{{{}^{{\rm{12}}}{{\rm{R}}_{{\rm{AOM}}}} - {}^{{\rm{12}}}{{\rm{R}}_{{\rm{CR}}}} - {}^{{\rm{12}}}{{\rm{R}}_{{\rm{ME - CH4}}}}}} \end{split} $$ (30)

    式中，L为CCH4-SMTZ与CCH4-bot之间的深度。

    联合方程式（28）和（30），可推导得到b的计算公式如下：

    $$b \!=\! \frac{{\left( {{}^{{\rm{13}}}{{\rm{R}}_{{\rm{AOM}}}}\! -\! {}^{{\rm{13}}}{{\rm{R}}_{{\rm{ME \!-\! CH4}}}}} \right) \!-\! \left( {{}^{{\rm{12}}}{{\rm{R}}_{{\rm{AOM}}}} \!-\! {}^{{\rm{12}}}{{\rm{R}}_{{\rm{ME \!-\! CH4}}}}} \right) \cdot {r_{{\rm{CH4 \!-\! bot}}}}}}{{{}^{{\rm{13}}}{{\rm{R}}_{{\rm{AOM}}}}/{\alpha _{CR}} \!-\! {}^{{\rm{12}}}{{\rm{R}}_{{\rm{AOM}}}} \times {r_{{\rm{CH4 \!-\! bot}}}}}}$$ (31)

    4.   结果与讨论

    4.1   孔隙水溶解组分的垂向分布及对AOM的指示

    孔隙水部分阴阳离子（SO42−、Ca2+和Mg2+）、DIC与δ13CDIC、CH4与δ13CCH4、δ13CCO2等地化参数的垂直分布剖面如图3所示。从图3可知，研究区沉积层中SO42-和CH4呈现出一个相反的垂直浓度变化剖面，于海底约6～8 mbsf深度形成硫酸盐—甲烷过渡带（SMTZ，图3中阴影层）。尽管沉积物孔隙水和顶空气样品的采集分辨率较低，但仍可直观看出：0～8 mbsf沉积层中SO42−浓度从底层海水浓度（15.95 mM）呈近似线性梯度地降低至检出限（0.20 mM）以下，并在SMTZ以下始终小于检出限；6～10 mbsf沉积层中CH4含量从0 mM近线性增加至1.83 mM，且在10～30 mbsf层段均含有甲烷，含量为0.48～1.84 mM，为含甲烷沉积层。

    此外，从图3还可看出，孔隙水DIC含量随深度增加而呈“先升后降”的变化规律，在SMTZ层位对应着DIC含量的高值水平（约18 mM）；而δ13CDIC的垂直分布与DIC大致呈“镜像”关系，在SMTZ层位出现最低值（约−17‰，VPDB）。在10 mbsf以下沉积层中，DIC含量骤降至接近底层海水浓度（2.46 mM）、且δ13CDIC逐渐升高至超过底层海水水平（−2.66‰，VPDB），表明此层段可能存在着较强的二氧化碳还原途径的微生物产甲烷作用。Ca2+和Mg2+的垂直分布剖面大致呈“先降后升”的规律，在SMTZ层位以上Ca2+和Mg2+含量水平从底层海水浓度6.33和38.04 mM分别下降至约2.0和30 mM，SMTZ层位以下又表现出一定程度的升高趋势，这可能与含甲烷沉积层中因DIC生成甲烷导致碳酸盐矿物溶解、或因DIC含量低时碳酸盐矿物结晶慢等有关。

    从反应方程式（3）可知，AOM作用将会引起SMTZ内孔隙水HCO3−浓度升高，导致DIC含量增加、以及Ca2+和Mg2+生成自生碳酸盐矿物出现浓度负异常[33-34]。因此，根据YS6孔SMTZ层位出现DIC高值、δ13CDIC、Ca2+和Mg2+低值的分布特征，表明研究区海底约6～8 mbsf沉积层中发生了AOM。

    从δ13CCH4的垂直分布特征（图3）还可发现：在SMTZ层位，AOM作用达到最大值时，δ13CCH4呈现出最小值（约−90.3‰，VPDB），较深部含甲烷沉积层（约20～30 mbsf）的甲烷碳同位素值（平均−71.3‰，VPDB）偏负约20‰，即SMTZ处CH4表现出明显的亏损13C现象。这种现象在卡斯卡迪亚边缘、墨西哥湾、智利边缘、奥尔胡斯湾等海底沉积物中SMTZ处亦被发现[7,10]。然而，甲烷发生AOM时，由于优先消耗12CH4，氧化后剩余留存在SMTZ内的CH4应当富集13C，δ13CCH4值应变重[35-36]。Beulig等[7]和Yoshinaga等[10]研究认为，出现13CH4亏损的原因可能是因为在SMTZ内部或下边缘存在着某种“隐藏的”产甲烷过程，微生物产甲烷作用对碳同位素的分馏效应常常导致最小δ13CCH4值在SMTZ处出现。

    4.2   孔隙水溶解组分经SMTZ的扩散通量及其指示

    根据孔隙水Ca2+、Mg2+、DIC、CH4和SO42−浓度的垂直分布数据，按照Fick第一定律（式10和11）计算得到各溶解组分进出SMTZ的扩散通量列于表1中。从表1可知，扩散进入SMTZ内的SO42−通量为11.87 mmol·m−2·a−1、CH4通量为4.65 mmol·m−2·a−1，SO42−通量为CH4通量的2.5倍。此外，依表1数据还可计算得到在SMTZ内部消耗的Ca2+和Mg2+净通量（公式12）为4.35 mmol·m−2·a−1、DIC净通量（△FDIC = FDIC.in− FDIC.out）为9.49 mmol·m−2·a−1，两者的消耗量并不具质量守恒规律。可见，在SMTZ内部可能存在碳酸盐沉淀之外的某种DIC碳源的消耗过程（如CR等）。

    表  1  孔隙水部分溶解组分的扩散通量
    Table  1.  Diffusion fluxes of dissolved components in pore water
    组分扩散系数D0/（m2·s−1）扩散通量F/（mmol·m−2·a−1）符号
    SO42−8.91E-1011.87FSO4.in
    CH41.39E-094.65FCH4.in
    Ca2+6.72E-102.29FCa.in
    0.76FCa.out
    Mg2+5.91E-104.15FMg.in
    1.33FMg.out
    DIC9.89E-1016.72FDIC.out
    7.23FDIC.in
    下载: 导出CSV 
    | 显示表格

    4.3   SMTZ内碳循环反应速率及其指示

    根据SMTZ内碳循环过程中轻、重碳同位素的质量守恒原则，取甲烷发生AOM时碳同位素分馏系数的平均值αAOM =1.017（αAOM变化区间为1.004～1.030[14,35]）按照3.3节的计算方法分别计算了OSR反应质量分数（f）为0.00、0.50和1.00时RAOM、RCR、RCP、ROM、RME、RME-CH4、RME-DIC和ROSR-C等轻、重碳同位素反应速率（12C和13C）和总反应速率，结果如表2。

    表  2  SMTZ内轻、重碳同位素反应速率和总反应速率
    Table  2.  The reaction rates of light and heavy carbon isotopes in SMTZ and the total reaction rates
    项目参数计算结果
    轻、重DIC通量（mmol·m−2·a−1）12FDIC.out−16.53
    12FDIC.in−7.15
    13FDIC.out−0.18
    13FDIC.in−0.08
    碳同位素值（‰）δ13COM−23.32
    δ13CCH4-SMTZ−90.34
    δ13CDIC-SMTZ−17.04
    δ13CCH4-bottom−71.31
    δ13CDIC-bottom−0.42
    13C/12C比值rstd0.011 237
    rOM0.010 975
    rCH4-SMTZ0.010 222
    rDIC-SMTZ0.011 046
    rCH4-bottom0.010 436
    rDIC-bottom0.011 232
    分馏系数αCR1.070 9
    αAOM1.017 0
    反应分数f0.000.501.00
    b0.370.370.37
    轻同位素反应速率
    （mmol·m−2·a−1）    		12RAOM11.879.267.35
    12RCR4.353.392.69
    12RCP4.314.314.31
    12ROM12.3410.449.04
    12RME12.345.220.00
    12RME-CH46.172.610.00
    12RME-DIC6.172.610.00
    12ROSR-C0.005.229.04
    重同位素反应速率
    （mmol·m−2·a−1）    		13RAOM0.120.090.07
    13RCR0.040.030.03
    13RCP0.050.050.05
    13ROM0.140.110.10
    13RME0.140.060.00
    13RME-CH40.060.030.00
    13RME-DIC0.070.030.00
    13ROSR-C0.000.060.10
    总反应速率
    （mmol·m−2·a−1）    		RAOM11.999.357.42
    RCR4.393.422.72
    RCP4.364.364.36
    ROM12.4810.559.14
    RME12.485.280.00
    RME-CH46.242.640.00
    RME-DIC6.242.640.00
    ROSR-C0.005.289.14
    甲烷通量绝对差值（mmol·m−2·a−1)△FCH43.291.360.06
    　　注：△FCH4= FCH4.in−（RAOM−RCR−RME-CH4）。
    下载: 导出CSV 
    | 显示表格

    从表2计算结果不难发现：仅当f=1.00时，利用“箱式模型”经轻、重碳同位素（13C/12C）质量平衡（3.3节）和Fick第一定律（3.2节）计算得到的甲烷扩散通量（FCH4）几乎无差异，绝对差值△FCH4= 0.06 mmol·m−2·a−1（表1）。f=1.00时，参与SMTZ内碳循环的AOM、CR、CP、OSR和ME 5种反应的总反应速率分别为7.42、2.72、4.36、9.14和0.00 mmol·m−2·a−1。由此可见，YS6孔柱状沉积物中SMTZ内部的确存在着一个“隐藏的”产甲烷过程—即CR，SMTZ内碳循环主要是OSR、AOM、CR和CP 4种生物地球化学作用的耦合过程，而ME可视作未发生。

    4.4   SMTZ内DIC碳库的“源—汇”贡献分析

    利用5种反应（AOM、OSR、CR、ME和CP）对DIC的生成或消耗的总反应速率（RAOM、ROSR-C、RCR、RME-DIC和RCP）与DIC净扩散通量（△FDIC= FDIC.in− FDIC.out）的比值，评估了参与碳循环的各途径对SMTZ内部沉积物孔隙水DIC碳库的“源—汇”相对贡献率（图4）。其中，AOM、OSR和ME作用对孔隙水DIC生成的贡献OSR＞AOM＞ME，相对贡献率分别为0.78、0.96和0.00，而CR和CP作用对孔隙水DIC消耗的贡献CP＞CR，相对贡献率分别为0.29和0.46。

    图  4  5种反应对SMTZ内DIC碳库的相对贡献率
    Figure  4.  The relative contribution of each reaction to pore water DIC in SMTZ
    下载: 全尺寸图片 幻灯片

    正是因为SMTZ内部存在着一个“隐藏的”产甲烷过程—CR，该过程引起的碳同位素分馏效应远大于AOM过程[7,13,37]，从而导致SMTZ内部或边缘沉积层中CH4出现明显的13C亏损现象（δ13CCH4-SMTZ约为−90.3‰），较深部含甲烷沉积层亏损量Δδ13CCH4≈20‰（图3）。此外，为了平衡CR作用对DIC碳库的消耗，必然促进有机质的OSR降解作用进一步增强，导致OSR生成DIC的贡献率较高，ROSR-C表现出高值（图4）。

    4.5   SMTZ内碳循环过程的研究意义

    通常，沉积物顶空气样品在采集过程中可能会出现一定程度的脱气逸散，导致直接利用甲烷浓度剖面数据计算得到的CH4扩散通量被低估[14,38]。从图3可知，SMTZ下方δ13CCH4和δ13CCO2随深度分布曲线呈现出光滑的、渐变式的“凹凸”剖面，以及两者之间具有相对稳定的差值（Δ δ13CCH4-CO2= δ13CCH4− δ13CCO2≈−67‰），表明采样过程中并未因脱气作用导致明显的碳同位素分馏[39]。因此，3.3节碳同位素质量平衡法计算得到的甲烷通量基本可靠。根据轻、重碳同位素（12C和13C）质量平衡和Fick第一定律计算得到的甲烷通量之间几乎无差异（△FCH4= FCH4.in−（RAOM−RCR−RME-CH4）≈0，表2），表明YS6孔在沉积物顶空气样品采集过程中，甲烷气体发生逸散的程度基本可以忽略。因此，由于穿过YS6孔沉积物中SMTZ的硫酸盐通量大于甲烷通量，表明深部含甲烷气层向上扩散而来的CH4并不是驱动SMTZ内部SO42−还原的唯一电子供体。

    根据SMTZ内碳循环研究结果可发现，导致硫酸盐通量大于甲烷通量的原因有两种：一是，SMTZ内有机质直接通过OSR作用消耗了部分SO42−；二是，间接通过“隐藏的”产甲烷作用（CR）为SMTZ内部提供了额外的CH4，从而“隐蔽地”发生了额外的AOM反应消耗了部分SO42−。正如扩散进入SMTZ中的甲烷通量（FCH4）为4.65 mmol·m−2·a−1，硫酸盐通量（FSO4）为11.87 mmol·m−2·a−1，初看AOM消耗硫酸盐的比例为39%，而因为CR作用为SMTZ内部补充了额外的CH4（RCR = 2.72 mmol·m−2·a−1），使得实际发生的AOM作用增加至62%。因此，在定量评价SMTZ内部AOM作用的相对强弱时，SMTZ内可能存在的“隐藏的”产甲烷过程（如CR、ME等）生成的甲烷所介导的硫酸盐还原作用不容忽视。

    5.   结论

    （1）研究区在海底约6～8 mbsf沉积层中形成SMTZ。基于“箱式模型”，SMTZ内发生的碳循环过程主要是有机质的硫酸盐还原作用（OSR）、碳酸盐沉淀作用（CP）、甲烷厌氧氧化作用（AOM）和“AOM生成的DIC”产甲烷作用（CR）等4种生物地球化学过程的耦合，而有机质的降解产甲烷作用（ME）未发生。

    （2）AOM、OSR和ME对孔隙水DIC生成的贡献OSR＞AOM＞ME，贡献率分别为0.78、0.96和0.00，而CR和CP对孔隙水DIC消耗的贡献CP＞CR，贡献率分别为0.29和0.46。研究区SMTZ下边缘沉积层中CH4出现明显的13C亏损现象，亏损量Δδ13CCH4约为20‰，很可能与SMTZ内存在CR作用有关。

    （3）研究区由深部含甲烷沉积层向上扩散而来的CH4并不是驱动SMTZ内部SO42−还原的唯一电子供体。研究SMTZ内碳循环过程，揭示SMTZ内是否存在“隐藏的”产甲烷作用（如CR、ME等），在定量评估海底沉积物中AOM作用的相对强弱时不容忽视。

    致谢：样品测试、数据分析和绘图等工作得到了青岛海洋地质研究所张媛媛、辛文彩和宋晓芸、徐翠玲、李晶等的热心帮助，在此表示最诚挚的感谢。

  • 图  1   研究区地质概况、马山柱状节理火山岩露头与镜下特征

    a. 研究区地质简图（据文献[13]），b. 马山石林剖面，c. 粗面英安岩中的晶洞构造，d. 马山粗面英安岩镜下显微特征，左侧为单偏光（-），右侧为正交光（+）。

    Figure  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   马山粗面英安岩锆石特征与测年结果

    a. 马山粗面英安岩的锆石阴极发光图像及对应的Pb206/U238年龄与εHf（t）值，b. 加权平均值，c. U-Pb年龄谐和图。

    Figure  2.   Zircon characteristics and dating results of the trachyandesite in Mashan

    a. 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.

    下载: 全尺寸图片 幻灯片

    图  3   马山粗面英安岩中锆石的稀土元素配分曲线[21]

    Figure  3.   Rare earth element distribution pattern of zircons from Mashan trachydacite[21]

    下载: 全尺寸图片 幻灯片

    图  4   马山粗面英安岩中锆石的εHf（t）值与对应的二阶段模式年龄

    Figure  4.   εHf（t） values and two stage Hf model ages of zircons from trachydacite in Mashan

    下载: 全尺寸图片 幻灯片

    图  5   马山柱状节理粗面英安岩岩石类型划分TAS图解[22]

    Figure  5.   TAS diagram of the trachydacite at Mashan[22]

    下载: 全尺寸图片 幻灯片

    图  6   马山粗面英安岩稀土元素配分曲线[21]

    Figure  6.   REE distribution pattern of the trachydacite at Mashan[21]

    下载: 全尺寸图片 幻灯片

    图  7   马山粗面英安岩微量元素蛛网图[21]

    Figure  7.   Spider distribution pattern of the trachydacite at Mashan[21]

    下载: 全尺寸图片 幻灯片

    图  8   马山粗面英安岩Sr/Y vs.Y[24]和(La/Yb)N-YbN[32]判别图解

    Figure  8.   Sr/Y -Y and (La/Yb)N-YbN adakitic trace elemental discrimination diagrams for Mashan trachydacite[24, 32]

    下载: 全尺寸图片 幻灯片

    表  1   马山粗面英安岩锆石U-Pb测试结果

    Table  1   Zircon U-Pb dating results of trachydacite at Mashan

    点 号Th /10−6U/10−6Th/U比值年龄/Ma谐和度*/%
    207Pb/206Pb± 1σ207Pb/235U± 1σ206Pb/238U± 1σ207Pb/206Pb± 1σ207Pb/235U± 1σ206Pb/238U± 1σ
    MS-013.2370.20.010.049950.002970.133780.007410.019840.000391951341277127299
    MS-02165.0235.10.700.103040.003872.149010.078400.150080.002761679691165259011574
    MS-030.00.450.010.034650.016824.371172.079280.034120.01459----17073932169156
    MS-041327.3935.31.420.047820.001970.118100.004840.017850.00027100871134114299
    MS-05376.9294.51.280.049510.003290.119910.007490.017870.000401721561157114399
    MS-06100.5115.30.870.056430.005250.136870.012190.017670.0005347820513011113385
    MS-07899.6806.21.120.051670.002460.128890.006000.017870.000373331051235114292
    MS-0896.9148.00.650.116180.003705.079290.160910.312620.0057018985718332717542895
    MS-091121.3809.31.390.050880.002580.126930.006280.017950.000362351141216115294
    MS-1080.6120.90.670.081390.004430.717900.040650.063210.001671231140549243951067
    MS-110.06.20.000.090500.022780.450480.120010.032430.003611436496378842062341
    MS-12510.81110.00.460.051500.002060.125020.005060.017270.00030265931205110291
    MS-13167.72466.90.070.169040.004069.842730.238790.412520.0063325484024202222262991
    MS-14184.5166.11.110.052290.004710.129910.010390.017950.000462982101249115392
    MS-1569.1246.80.280.055780.002570.267340.011740.034360.000624431022419218490
    MS-16109.4315.40.350.052550.002370.241680.010640.032710.00060309992209208494
    MS-17937.5625.81.500.052140.002390.130700.005110.017940.000343001061255115291
    MS-18396.1555.50.710.053030.002220.130040.004860.017490.00032332941244112290
    MS-1952.91325.40.040.050070.001460.226080.006220.032090.00044198692075204398
    MS-20614.1290.22.120.046500.005310.114270.014120.017150.000483324311013110399
    MS-213467.11542.52.250.046980.001560.117370.004140.017660.0002956721134113299
    MS-221321.0881.11.500.047500.001890.120300.004820.017970.0003476891154115299
    MS-23412.3496.50.830.046820.002510.112490.005550.017280.00040391261085110397
    　　注：*谐和度=100${\rm{ \times \{ 1 - [}}{{\rm{(}}^{{\rm{206}}}}{\rm{Pb}}{{\rm{/}}^{{\rm{238}}}}{\rm{U}}\;{{\rm{ - }}^{{\rm{207}}}}{\rm{Pb}}{{\rm{/}}^{{\rm{235}}}}{\rm{U}}\;{\rm{)/(}}{{\rm{(}}^{{\rm{206}}}}{\rm{Pb}}{{\rm{/}}^{{\rm{238}}}}{\rm{U}}\;{\rm{}}{{\rm{ + }}^{{\rm{207}}}}{\rm{Pb}}{{\rm{/}}^{{\rm{235}}}}{\rm{U}}\;{\rm{)/2)]\} }}$
    下载: 导出CSV

    表  2   马山粗面英安岩锆石原位Lu-Hf测试结果

    Table  2   Zircon in-situ Lu-Hf results of trachydacite at Mashan

    点 号176Yb/177Hf±2σ176Lu/177Hf±2σ176Hf/177Hf±2σt (Ma)①T C ② DM±2σεHf(t)③±2σ
    MS-010.0107700.0001130.0004680.0000040.2823400.000014127197132−12.50.51
    MS-040.0926220.0021840.0034620.0000740.2822420.000021114220947−16.50.76
    MS-050.0800540.0015200.0030270.0000390.2822390.000021114221346−16.60.74
    MS-070.1451240.0042300.0052050.0001410.2822210.000020114226344−17.40.71
    MS-080.0049210.0002090.0001250.0000050.2812650.0000151898326333−11.10.54
    MS-090.1055450.0008740.0038230.0000370.2822140.000017115227238−17.50.61
    MS-120.1332180.0015710.0049970.0000550.2822320.000033110223972−17.01.17
    MS-130.0221710.0009790.0007140.0000250.2811770.0000272548308758−0.40.96
    MS-140.0640950.0016470.0025420.0000620.2821910.000029115231763−18.21.01
    MS-150.0224320.0005870.0008750.0000190.2821760.000016218228736−16.40.58
    MS-160.0137030.0001350.0006270.0000050.2822940.000012208203028−12.50.44
    MS-170.1102100.0003950.0041990.0000160.2821590.000017115239538−19.50.61
    MS-190.0241160.0005540.0011120.0000240.2821020.000012204245827−19.40.44
    MS-200.0828070.0014610.0033160.0000560.2822250.000025110224855−17.20.88
    MS-210.1500460.0041720.0052990.0001460.2822630.000023113217251−15.90.82
    MS-220.1322880.0024290.0051100.0000930.2822540.000022115219148−16.20.77
    MS-230.0741510.0040700.0028050.0001470.2821870.000018110233039−18.50.64
    　　注：①为锆石表观年龄；②为基于DM的二阶段模式年龄计算使用平均大陆地壳176Lu/177Hf值0.015。表中数据的误差为2σ（标准误差），不确定度为最后一位数字；③为εHf(t)值的计算使用锆石的表观年龄，该年龄为锆石结晶年龄的最小估计值，λ176Lu=1.867×10−11/a[17]，CHUR现今的176Lu/177Hf及176Hf/177Hf值分别为0.0332和0.282772[18]，DM现今的176Lu/177Hf及176Hf/177Hf值分别为0.0384和0.28325[19]。
    下载: 导出CSV

    表  3   马山粗面英安岩全岩主量元素分析结果

    Table  3   Whole-rock major element compositions of trachydacite at Mashan

    %
    样品号M-01M-02M-03M-04平均值
    SiO265.4467.2667.0667.1866.74
    TiO20.460.450.440.400.44
    Al2O314.5813.3713.4313.3813.69
    Fe2O33.613.893.653.623.69
    MnO0.080.100.080.090.09
    MgO2.101.952.142.362.14
    CaO3.402.672.762.582.85
    Na2O3.774.484.714.414.34
    K2O3.833.763.583.983.79
    P2O50.180.190.210.170.19
    L.O.I2.121.461.861.651.77
    下载: 导出CSV

    表  4   马山粗面英安岩微量元素分析结果

    Table  4   Trace element compositions of trachydacite at Mashan

    10−6
    样品号M-01M-02M-03M-04平均值
    Li37.541.436.438.638.5
    Be1.751.661.701.681.70
    Sc18.307.807.918.1010.50
    Ti28232700264024002641
    V53.155.156.154.454.7
    Cr68.563.565.681.669.8
    Co10.8011.0011.0011.4011.10
    Ni32.129.529.639.232.6
    Cu4.7645.8028.6028.2026.80
    Zn48.572.268.067.264.0
    Ga17.916.717.215.916.9
    Rb91.685.087.199.290.7
    Sr402341422489414
    Y13.413.613.913.013.5
    Zr130141134126133
    Nb9.511.010.810.610.5
    Mo0.610.880.580.660.68
    Cd0.0360.0610.0660.0380.051
    Cs0.880.620.651.000.79
    Ba23811461193917271877
    Hf3.734.354.404.004.12
    Ta0.731.201.001.151.02
    Tl0.380.380.380.440.40
    Pb30.354.639.243.641.9
    Th6.465.565.385.265.67
    U1.691.501.541.541.57
    下载: 导出CSV

    表  5   马山粗面英安岩稀土元素分析结果

    Table  5   REE compositions of trachydacite at Mashan

    10−6
    样品号M-01M-02M-03M-04平均值
    La29.629.930.827.829.5
    Ce52.453.154.649.552.4
    Pr6.346.186.345.776.16
    Nd22.923.324.021.823.0
    Sm3.963.863.923.643.85
    Eu0.941.081.101.001.03
    Gd3.533.283.383.103.32
    Tb0.500.470.480.460.48
    Dy2.502.532.602.462.52
    Ho0.500.510.520.480.50
    Er1.391.421.441.381.41
    Tm0.240.220.220.220.23
    Yb1.451.461.441.401.44
    Lu0.230.230.230.220.23
    Eu/Eu*0.770.930.920.910.88
    LaN/YbN13.8913.9414.5613.5113.98
    ΣREE126.48127.54131.07119.23126.08
    下载: 导出CSV
    • 参考文献(57)
  • [1]

    Cope T. Phanerozoic magmatic tempos of North China [J]. Earth and Planetary Science Letters, 2017, 468: 1-10. doi: 10.1016/j.jpgl.2017.03.022
    [2] 任纪舜, 陈廷愚, 牛宝贵, 等. 中国东部及邻区大陆岩石圈的构造演化与成矿[M]. 北京: 科学出版社, 1990: 1-205.

    REN Jishun, CHEN Tingyu, NIU Baogui, et al. Tectonic Evolution of the Continental Lithosphere and Metallogeny in Eastern China and Adjacent Areas[M]. Beijing: Science Press, 1990: 1-205.
    [3] 邓晋福, 刘厚祥, 赵海玲, 等. 燕辽地区燕山期火成岩与造山模型[J]. 现代地质, 1996, 10(2):137-148

    DENG Jinfu, LIU Houxiang, ZHAO Hailing, et al. Yanshanian igneous rocks and orogeny model in Yanshan-Liaoning area [J]. Geoscience, 1996, 10(2): 137-148.
    [4] 董树文, 吴锡浩, 吴珍汉, 等. 论东亚大陆的构造翘变: 燕山运动的全球意义[J]. 地质论评, 2000, 46(1):8-13 doi: 10.3321/j.issn:0371-5736.2000.01.002

    DONG Shuwen, WU Xihao, WU Zhenhan, et al. On tectonic seesawing of the East Asia Continent: Global implication of the Yanshanian Movement [J]. Geological Review, 2000, 46(1): 8-13. doi: 10.3321/j.issn:0371-5736.2000.01.002
    [5] 董树文, 张岳桥, 李海龙, 等. “燕山运动”与东亚大陆晚中生代多板块汇聚构造: 纪念“燕山运动”90周年[J]. 中国科学: 地球科学, 2018, 61(6):1888-1909

    DONG Shuwen, ZHANG Yueqiao, LI Hailong, et al. The Yanshan orogeny and late Mesozoic multi-plate convergence in East Asia: Commemorating 90th years of the “Yanshan Orogeny” [J]. Science China Earth Sciences, 2018, 61(6): 1888-1909.
    [6] 张旗, 钱青, 王二七, 等. 燕山中晚期的中国东部高原: 埃达克岩的启示[J]. 地质科学, 2001, 36(2):248-255 doi: 10.3321/j.issn:0563-5020.2001.02.014

    ZHANG Qi, QIAN Qing, WANG Erqi, et al. An East China Plateau in Mid-Late Yanshanian Period: Implication from adakites [J]. Chinese Journal of Geology, 2001, 36(2): 248-255. doi: 10.3321/j.issn:0563-5020.2001.02.014
    [7] 张旗, 许继峰, 王焰, 等. 埃达克岩的多样性[J]. 地质通报, 2004, 23(9-10):959-965

    ZHANG Qi, XU Jifeng, WANG Yan, et al. Diversity of adakite [J]. Geological Bulletin of China, 2004, 23(9-10): 959-965.
    [8] 张旗, 焦守涛. 埃达克岩来自高压背景: 一个科学的、可靠的、有预见性的科学发现[J]. 岩石学报, 2020, 36(6):1675-1683 doi: 10.18654/1000-0569/2020.06.02

    ZHANG Qi, JIAO Shoutao. Adakite comes from a high-pressure background: A scientific, reliable, predictable scientific discovery [J]. Acta Petrologica Sinica, 2020, 36(6): 1675-1683. doi: 10.18654/1000-0569/2020.06.02
    [9] 汪洋, 程素华. “C型埃达克岩”: 一个基于误解的概念?[J]. 矿物岩石地球化学通报, 2010, 29(3):284-292, 298 doi: 10.3969/j.issn.1007-2802.2010.03.012

    WANG Yang, CHENG Suhua. The “C-type Adakite”: A concept based on misunderstanding? [J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2010, 29(3): 284-292, 298. doi: 10.3969/j.issn.1007-2802.2010.03.012
    [10] 汪洋. 涉及埃达克岩实验的若干基本概念与问题辨析[J]. 岩石矿物学杂志, 2016, 35(1):162-176 doi: 10.3969/j.issn.1000-6524.2016.01.012

    WANG Yang. A discussion on some basic conceptions and problems related to the experimental study of adakite [J]. Acta Petrologica Et Mineralogica, 2016, 35(1): 162-176. doi: 10.3969/j.issn.1000-6524.2016.01.012
    [11]

    Martin H, Smithies R H, Rapp R, et al. D. An overview of adakite, tonalite–trondhjemite–granodiorite (TTG), and sanukitoid: relationships and some implications for crustal evolution [J]. Lithos, 2005, 79(1-2): 1-24. doi: 10.1016/j.lithos.2004.04.048
    [12] 王强, 郝露露, 张修政, 等. 汇聚板块边缘的埃达克质岩: 成分和成因[J]. 中国科学: 地球科学, 2020, 63(12):1992-2016 doi: 10.1007/s11430-020-9678-y

    WANG Qiang, HAO Lulu, ZHANG Xiuzheng, et al. Adakitic rocks at convergent plate boundaries: Compositions and petrogenesis [J]. Science China Earth Sciences, 2020, 63(12): 1992-2016. doi: 10.1007/s11430-020-9678-y
    [13] 山东省地质调查院. 1: 25万青岛市幅区域地质调查报告与地质图[R]. 2004.

    Shandong Institute of Geological Survey. 1: 250000 Qingdao Regional Geological Survey Report and Geological Map[R]. 2004.
    [14] 韩宗珠, 衣伟虹, 余芝华, 等. 青岛马山地区安山玢岩地球化学特征及成因机理研究[J]. 中国海洋大学学报, 2011, 41(S1):305-311

    HAN Zongzhu, YI Weihong, YU Zhihua, et al. Geochemistry and petrogenesis of andesite-prophyrite in Mashan area, Qingdao [J]. Periodical of Ocean University of China, 2011, 41(S1): 305-311.
    [15] 何登洋, 邱昆峰, 于皓丞, 等. 华北克拉通胶莱盆地马山地区早白垩世粗面英安岩岩石成因[J]. 岩石学报, 2020, 36(12):3705-3720 doi: 10.18654/1000-0569/2020.12.09

    HE Dengyang, QIU Kunfeng, YU Haocheng, et al. Petrogenesis of the Early Cretaceous trachy-dacite from Mashan in the Jiaolai Basin, North China Craton [J]. Acta Petrologica Sinica, 2020, 36(12): 3705-3720. doi: 10.18654/1000-0569/2020.12.09
    [16]

    Pearce N J G, Perkins W T, Westgate J A, et al. A compilation of new and published major and trace element data for NIST SRM 610 and NIST SRM 612 glass reference materials [J]. Geostandards Newsletter, 1997, 21(1): 115-144. doi: 10.1111/j.1751-908X.1997.tb00538.x
    [17]

    Söderlund U, Patchett P J, Vervoort J D, et al. The 176Lu decay constant determined by Lu-Hf and U-Pb isotope systematics of Precambrian mafic intrusions [J]. Earth and Planetary Science Letters, 2004, 219(3-4): 311-324. doi: 10.1016/S0012-821X(04)00012-3
    [18]

    Bichert-Toft J, Albarède F. The Lu-Hf isotope geochemistry of chondrites and the evolution of the mantle-crust system [J]. Earth and Planetary Science Letters, 1997, 148(1-2): 243-258. doi: 10.1016/S0012-821X(97)00040-X
    [19]

    Griffin W L, Pearson N J, Belousova E, et al. The Hf isotope composition of cratonic mantle: LAM-MC-ICPMS analysis of zircon megacrysts in kimberlites [J]. Geochimica et Cosmochimica Acta, 2000, 64(1): 133-147. doi: 10.1016/S0016-7037(99)00343-9
    [20] 吴元保, 郑永飞. 锆石成因矿物学研究及其对U-Pb年龄解释的制约[J]. 科学通报, 2004, 49(15):1554-1569 doi: 10.1007/BF03184122

    WU Yuanbao, ZHENG Yongfei. Genesis of zircon and its constraints on interpretation of U-Pb age [J]. Chinese Science Bulletin, 2004, 49(15): 1554-1569. doi: 10.1007/BF03184122
    [21]

    McDonough W F, Sun S S. The composition of the Earth [J]. Chemical Geology, 1995, 120(3-4): 223-253. doi: 10.1016/0009-2541(94)00140-4
    [22]

    Le Bas M J, Le Maitre R W, Streckeisen A, et al. A chemical classification of volcanic rocks based on the total alkali-silica diagram [J]. Journal of Petrology, 1986, 27(3): 745-750. doi: 10.1093/petrology/27.3.745
    [23]

    Kay R W. Aleutian magnesian andesites: Melts from subducted Pacific ocean crust [J]. Journal of Volcanology and Geothermal Research, 1978, 4(1-2): 117-132. doi: 10.1016/0377-0273(78)90032-X
    [24]

    Defant M J, Drummond M S. Derivation of some modern arc magmas by melting of young subducted lithosphere [J]. Nature, 1990, 347(6294): 662-665. doi: 10.1038/347662a0
    [25]

    Defant M J, Xu J F, Kepezhinskas P, et al. Adakites: some variations on a theme [J]. Acta Petrologica Sinica, 2002, 18(2): 129-142.
    [26] 王焰, 张旗, 钱青. 埃达克岩(adakite)的地球化学特征及其构造意义[J]. 地质科学, 2000, 35(2):251-256 doi: 10.3321/j.issn:0563-5020.2000.02.016

    WANG Yan, ZHANG Qi, QIAN Qing. Adakite: geochemical characteristics and tectonic significances [J]. Scientia Geologica Sinica, 2000, 35(2): 251-256. doi: 10.3321/j.issn:0563-5020.2000.02.016
    [27] 王强, 许继锋, 赵振华. 一种新的火成岩: 埃达克岩的研究综述[J]. 地球科学进展, 2001, 16(2):201-208 doi: 10.3321/j.issn:1001-8166.2001.02.010

    WANG Qiang, XU Jifen, ZHAO Zhenghua. The summary and comment on research on a new kind of igneous rock-Adakite [J]. Advance in Earth Sciences, 2001, 16(2): 201-208. doi: 10.3321/j.issn:1001-8166.2001.02.010
    [28]

    Rapp R P, Xiao L, Shim Izu N. Experimental constraints on the origin of potassium-rich adakites in eastern China [J]. Acta Petrologica Sinica, 2002, 18(3): 293-302.
    [29] 罗照华, 柯珊, 谌宏伟. 埃达克岩的特征、成因及构造意义[J]. 地球科学进展, 2002, 21(7):436-440

    LUO Zhaohua, KE Shan, CHEN Hongwei. Characteristics, petrogenesis and tectonic implications of adakite [J]. Geological Bulletin of China, 2002, 21(7): 436-440.
    [30] 翟明国. 埃达克岩和大陆下地壳重熔的花岗岩类[J]. 岩石学报, 2004, 20(2):193-194 doi: 10.3969/j.issn.1000-0569.2004.02.001

    ZHAI Mingguo. Adakite and related granitoids from partial melting of continental lower crust [J]. Acta Petrologica Sinica, 2004, 20(2): 193-194. doi: 10.3969/j.issn.1000-0569.2004.02.001
    [31] 张旗, 王元龙, 金惟俊, 等. 晚中生代的中国东部高原: 证据、问题和启示[J]. 地质通报, 2008, 27(9):1404-1430 doi: 10.3969/j.issn.1671-2552.2008.09.004

    ZHANG Qi, WANG Yuanlong, JI Weijun, et al. Eastern China Plateau during the Late Mesozoic: evidence, problems and implication [J]. Geological Bulletin of China, 2008, 27(9): 1404-1430. doi: 10.3969/j.issn.1671-2552.2008.09.004
    [32]

    Martin H. Effect of steeper Archean geothermal gradient on geochemistry of subduction-zone magmas [J]. Geology, 1986, 14(9): 753-756. doi: 10.1130/0091-7613(1986)14<753:EOSAGG>2.0.CO;2
    [33] 邓晋福, 苏尚国, 刘翠, 等. 关于华北克拉通燕山期岩石圈减薄的机制与过程的讨论: 是拆沉, 还是热侵蚀和化学交代?[J]. 地学前缘, 2006, 13(2):105-119 doi: 10.3321/j.issn:1005-2321.2006.02.009

    DENG Jinfu, SU Shangguo, LIU Cui, et al. Discussion on the lithospheric thinning of the North China craton: delamination? Or thermal erosion and chemical metasomatism? [J]. Earth Science Frontiers, 2006, 13(2): 105-119. doi: 10.3321/j.issn:1005-2321.2006.02.009
    [34]

    Best M G. Igneous and Metamorphic Petrology[M]. 2nd ed. Malden, Mass: Blackwell Science, 2003: 340-341.
    [35]

    Taylor S R, McClennan S M. The Continental Crust: Its Composition and Evolution[M]. Oxford: Blackwell Scientific Publication, 1985: 312.
    [36] 张旗, 王焰, 熊小林, 等. 埃达克岩和花岗岩: 挑战与机遇[M]. 北京: 中国大地出版社, 2008: 1-344.

    ZHANG Qi, WANG Yan, XIONG Xiaolin, et al. Adakite and Granite: Challenges and Opportunities[M]. Beijing: China Land Press, 2008: 1-344.
    [37] 赵越, 徐刚, 张拴宏, 等. 燕山运动与东亚构造体制的转变[J]. 地学前缘, 2004, 11(3):319-328 doi: 10.3321/j.issn:1005-2321.2004.03.030

    ZHAO Yue, XU Gang, ZHANG Shuanhong, et al. Yanshanian movement and conversion of tectonic regimes in East Asia [J]. Earth Science Frontiers, 2004, 11(3): 319-328. doi: 10.3321/j.issn:1005-2321.2004.03.030
    [38] 孟庆任. 中国东部燕山期沉积盆地演化与资源效应[J]. 矿物岩石地球化学通报, 2017, 36(4):567-569 doi: 10.3969/j.issn.1007-2802.2017.04.007

    MENG Qingren. Development of sedimentary basins in Eastern China during the Yanshanian period [J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2017, 36(4): 567-569. doi: 10.3969/j.issn.1007-2802.2017.04.007
    [39] 邵济安, 刘福田, 陈辉, 等. 大兴安岭—燕山晚中生代岩浆活动与俯冲作用关系[J]. 地质学报, 2001, 75(1):56-63 doi: 10.3321/j.issn:0001-5717.2001.01.006

    SHAO Ji’an, LIU Futian, CHEN Hui, et al. Relationship between Mesozoic Magmatism and Subduction in Da Hinggan-Yanshan Area [J]. Acta Geologica Sinica, 2001, 75(1): 56-63. doi: 10.3321/j.issn:0001-5717.2001.01.006
    [40] 胡受奚, 赵乙英, 胡志宏, 等. 中国东部中—新生代活动大陆边缘构造: 岩浆作用演化和发展[J]. 岩石学报, 1994, 10(4):370-381 doi: 10.3321/j.issn:1000-0569.1994.04.003

    HU Shouxi, ZHAO Yiying, HU Zhihong, et al. Evolution and development of tectonics and magmatism at the active continental margin of the East China (E106°) during Mesozoic and Cenozoic [J]. Acta Petrologica Sinica, 1994, 10(4): 370-381. doi: 10.3321/j.issn:1000-0569.1994.04.003
    [41] 王德滋, 任启江, 邱检生, 等. 中国东部橄榄安粗岩省的火山岩特征及其成矿作用[J]. 地质学报, 1996, 70(1):23-34

    WANG Dezi, REN Qijiang, QIU Jiansheng, et al. Characteristics of volcanic rocks in the Shoshonite province, Eastern China, and their metallogenesis [J]. Acta Geologica Sinica, 1996, 70(1): 23-34.
    [42] 万天丰, 赵庆乐. 中国东部构造-岩浆作用的成因[J]. 中国科学: 地球科学, 2012, 55(3):347-354 doi: 10.1007/s11430-011-4361-4

    WAN Tianfeng, ZHAO Qingle. The genesis of tectono-magmatism in eastern China [J]. Science China Earth Sciences, 2012, 55(3): 347-354. doi: 10.1007/s11430-011-4361-4
    [43] 张旗. 中国东部中生代岩浆活动与太平洋板块向西俯冲有关吗?[J]. 岩石矿物学杂志, 2013, 32(1):113-128 doi: 10.3969/j.issn.1000-6524.2013.01.010

    ZHANG Qi. Is the Mesozoic magmatism in eastern China related to westward subduction of the Pacific plate? [J]. Acta Petrologica ET Mineralogica, 2013, 32(1): 113-128. doi: 10.3969/j.issn.1000-6524.2013.01.010
    [44]

    Zhu X Q, Zhu W B, Ge R F, et al. Late Paleozoic provenance shift in the south-central North China Craton: Implications for tectonic evolution and crustal growth [J]. Gondwana Research, 2014, 25(1): 383-400. doi: 10.1016/j.gr.2013.04.009
    [45]

    Li S Z, Zhao S J, Liu X, et al. Closure of the Proto-Tethys Ocean and Early Paleozoic amalgamation of microcontinental blocks in East Asia [J]. Earth-Science Reviews, 2018, 186: 37-75. doi: 10.1016/j.earscirev.2017.01.011
    [46] 张岳桥, 李金良, 张田, 等. 胶莱盆地及其邻区白垩纪—古新世沉积构造演化历史及其区域动力学意义[J]. 地质学报, 2008, 82(9):1229-1257 doi: 10.3321/j.issn:0001-5717.2008.09.007

    ZHANG Yueqiao, LI Jinliang, ZHANG Tian, et al. Cretaceous to Paleocene Tectono-sedimentary evolution of the Jiaolai Basin and the contiguous areas of the Shandong Peninsula (North China) and Its geodynamic implications [J]. Acta Geologica Sinica, 2008, 82(9): 1229-1257. doi: 10.3321/j.issn:0001-5717.2008.09.007
    [47] 张田, 张岳桥. 胶北隆起晚中生代构造-岩浆演化历史[J]. 地质学报, 2008, 82(9):1210-1228 doi: 10.3321/j.issn:0001-5717.2008.09.006

    ZHANG Tian, ZHANG Yueqiao. Late Mesozoic tectono-magmatic evolution history of the Jiaobei Uplift, Shandong Peninsula [J]. Acta Geologica Sinica, 2008, 82(9): 1210-1228. doi: 10.3321/j.issn:0001-5717.2008.09.006
    [48]

    Guo Z X, Yang Y T, Zhao X Z, et al. Early Cretaceous tectonostratigraphic evolution of the Erlian Basin, NE China: A record of Late Mesozoic intraplate deformation in East Asia [J]. Marine and Petroleum Geology, 2019, 110: 539-564. doi: 10.1016/j.marpetgeo.2019.07.043
    [49] 吕洪波, 王俊, 张海春. 山东灵山岛晚中生代滑塌沉积层的发现及区域构造意义初探[J]. 地质学报, 2011, 85(6):938-946

    LÜ Hongbo, WANG Jun, ZHANG Haichun. Discovery of the Late Mesozoic slump beds in Lingshan Island, Shandong, and a pilot research on the regional tectonics [J]. Acta Geologica Sinica, 2011, 85(6): 938-946.
    [50] 孙天柱, 武斌. 临沂方城盆地中生代青山群火山岩地球化学特征及地质意义[J]. 山东国土资源, 2020, 36(6):14-22

    SUN Tianzhu, WU Bin. Geochemical characteristics and geological significance of volcanic rocks in Mesozoic Qingshan Group in Fangcheng Basin of Linyi City [J]. Shandong Land and Resources, 2020, 36(6): 14-22.
    [51] 王宇林, 郭强, 李喜海, 等. 辽西喀左盆地九佛堂组的划分及相关问题[J]. 地层学杂志, 2009, 33(2):147-154 doi: 10.3969/j.issn.0253-4959.2009.02.004

    WANG Yulin, GUO Qiang, LI Xihai, et al. Subdivision of the Jiufotang Formation and related issues in the Kazuo Basin [J]. Journal of Stratigraphy, 2009, 33(2): 147-154. doi: 10.3969/j.issn.0253-4959.2009.02.004
    [52] 刘为付, 刘双龙, 孙立新, 等. 包日温都构造带九佛堂组高分辨率层序结构及储层分析[J]. 地层学杂志, 2015, 39(3):310-318

    LIU Weifu, LIU Shuanglong, SUN Lixin, et al. High resolution sequence stratigraphy configurations and reservoir analysis of the Jiufotang Formation in the Baoriwendu tectonic belt [J]. Journal of Stratigraphy, 2015, 39(3): 310-318.
    [53] 姜同海. 胶莱盆地构造特征及演化[D]. 山东科技大学硕士学位论文, 2005: 1-45.

    JIANG Tonghai. The structural characteristics and evolution of Jiaolai Basin[D]. Master Dissertation of Shandong University of Science and Technology, 2005: 1-45.
    [54]

    Wu L, Monié P, Wang F, et al. Multi-phase cooling of Early Cretaceous granites on the Jiaodong Peninsula, East China: Evidence from 40Ar/39Ar and (U-Th)/He thermochronology [J]. Journal of Asian Earth Sciences, 2018, 160: 334-347. doi: 10.1016/j.jseaes.2017.11.014
    [55]

    Zhang L M, Wang C S, Cao K, et al. High elevation of Jiaolai Basin during the Late Cretaceous: Implication for the coastal mountains along the East Asian margin [J]. Earth and Planetary Science Letters, 2016, 456: 112-123. doi: 10.1016/j.jpgl.2016.09.034
    [56]

    Yang Y T. An unrecognized major collision of the Okhotomorsk Block with East Asia during the Late Cretaceous, constraints on the plate reorganization of the Northwest Pacific [J]. Earth-Science Reviews, 2013, 126: 96-115. doi: 10.1016/j.earscirev.2013.07.010
    [57]

    Chen P J. Paleoenvironmental changes during the Cretaceous in eastern China[M]//Okada H, Mateer N J. Developments in Palaeontology and Stratigraphy. Amsterdam: Elsevier, 2000, 17: 81-90.
    • 相关文章
    • 施引文献(5)

  • 期刊类型引用(4)

    1. 段晓勇，印萍，谢永清，曹珂，仇建东，李梅娜，李雪. 中国近海浅层气调查研究进展. 海洋地质与第四纪地质. 2024(03): 183-196 . 本站查看
    2. 李凤，段晓勇，贺行良，张媛媛，印萍，谢永清，杨磊，董超. 杭州湾沉积物中甲烷的垂向分布及控制因素. 海洋地质与第四纪地质. 2023(06): 112-121 . 本站查看
    3. 李晶，刘昌岭，吴能友，贺行良，孟庆国，许晓晴，陈烨. 海洋环境中甲烷好氧氧化过程的研究进展. 海洋地质与第四纪地质. 2021(05): 67-76 . 本站查看
    4. 常鑫，张明宇，谷玉，王厚杰，刘喜停. 黄、东海陆架泥质区自生黄铁矿成因及其控制因素. 地球科学进展. 2020(12): 1306-1320 . 百度学术

    其他类型引用(1)

    • 资源附件(0)
图(8)  /  表(5)
计量
  • 文章访问数:  2031
  • HTML全文浏览量:  514
  • PDF下载量:  48
  • 被引次数: 5
出版历程
  • 收稿日期:  2021-01-17
  • 修回日期:  2021-04-29
  • 网络出版日期:  2021-12-19
  • 刊出日期:  2021-12-27

目录

摘要
HTML全文

  • 1.   地质背景

  • 2.   样品与方法

  • 3.   碳循环模型与计算

  • 3.1   箱式模型

  • 3.2   质量平衡

  • 3.3   碳同位素质量平衡

  • 4.   结果与讨论

  • 4.1   孔隙水溶解组分的垂向分布及对AOM的指示

  • 4.2   孔隙水溶解组分经SMTZ的扩散通量及其指示

  • 4.3   SMTZ内碳循环反应速率及其指示

  • 4.4   SMTZ内DIC碳库的“源—汇”贡献分析

  • 4.5   SMTZ内碳循环过程的研究意义

  • 5.   结论

参考文献(57)
相关文章
施引文献(5)
资源附件(0)

/

返回文章
返回
微信公众号 微信公众号
订购码 订购码

版权所有 © 《海洋地质与第四纪地质》编辑部    备案号：鲁ICP备15036216号-7

地址：山东省青岛市即墨区观山路596号    邮编：266237    联系电话：0532-85755823    E-mail：hydzdsjdz@163.com

本系统由北京仁和汇智信息技术有限公司开发 技术支持： info@rhhz.net 百度统计