胶北地体荆山群大理岩碳氧同位素地球化学特征及其对Lomagundi-Jatuli事件的指示

引用本文

李凯月, 汤好书, 陈衍景, 薛莅治, 王玭, 孙之夫. 2020. 胶北地体荆山群大理岩碳氧同位素地球化学特征及其对Lomagundi-Jatuli事件的指示. 岩石学报, 36(4): 1059-1075, doi: 10.18654/1000-0569/2020.04.06

Li KY, Tang HS, Chen YJ, Xue LZ, Wang P and Sun ZF. 2020. Carbon and oxygen isotope geochemistry of marbles in the Jingshan Group, Jiaobei Terrane and its indication to the Lomagundi-Jatuli Event. Acta Petrologica Sinica, 36(4): 1059-1075(in Chinese with English abstract)

胶北地体荆山群大理岩碳氧同位素地球化学特征及其对Lomagundi-Jatuli事件的指示

李凯月¹, 汤好书², 陈衍景¹, 薛莅治¹, 王玭³, 孙之夫⁴

1. 北京大学造山带与地壳演化教育部重点实验室, 北京 100871;
2. 中国科学院地球化学研究所, 矿床地球化学国家重点实验室, 贵阳 550008;
3. 中国科学院南海海洋研究所边缘海与大洋地质重点实验室, 广州 510301;
4. 山东黄金集团有限公司, 莱州 261442

2019-04-16 收稿; 2019-11-25 改回.

基金项目: 本文受国家自然科学基金项目（U1906207、41890833、41872085）和国家重点基础研究发展计划（2012CB416602）联合资助

第一作者简介: 李凯月, 女, 1993年生, 博士生, 矿物学、岩石学、矿床学专业, E-mail:kaiyue_li@163.com

通讯作者: 陈衍景, 男, 1962年生, 博士, 教授, 矿床学研究方向, E-mail:yjchen@pku.edu.cn.

摘要: 古元古代大气圈快速大量充氧，地球表生系统发生一系列重要变化，其中之一是2.3~2.06Ga碳酸盐岩碳同位素正异常现象，称Lomagundi-Jatuli事件，颇受瞩目，成为当前国际研究前沿。胶北地体荆山群普遍经历高角闪岩相-麻粒岩相变质作用，其下部禄格庄组碳酸盐地层发育。本文研究了山后地区荆山群禄格庄组26件大理岩和7件白云质石英岩样品，获得δ¹³C_carb值在-0.8~+3.4‰（V-PDB），平均+1.4‰，其中大理岩δ¹³C平均值为+1.6‰，白云质石英岩δ¹³C平均值为+0.9‰。氧同位素值变化于9.3‰~20.2‰（V-SMOW），平均14.0‰。样品碳-氧同位素值之间具有很好的正相关性，且有13件样品Mn/Sr比值大于6，说明高角闪岩相的变质作用可能使δ¹³C_carb下降了3‰~5‰，原始沉积碳酸盐的δ¹³C_carb值高达+8.4‰，明显正异常，是全球性Lomagundi-Jatuli事件的响应。

关键词: Lomagundi-Jatuli事件大氧化事件碳同位素异常荆山群禄格庄组胶北地体

Carbon and oxygen isotope geochemistry of marbles in the Jingshan Group, Jiaobei Terrane and its indication to the Lomagundi-Jatuli Event

LI KaiYue¹, TANG HaoShu², CHEN YanJing¹, XUE LiZhi¹, WANG Pin³, SUN ZhiFu⁴

1. MOE Key Laboratory of Orogenic Belts and Crustal Evolution, Peking University, Beijing 100871, China;
2. State Key Laboratory of Ore Deposit Geochemistry, Institute of Geochemistry, Chinese Academy of Science, Guiyang 550008, China;
3. Key Laboratory of Ocean and Marginal Sea Geology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou 510301, China;
4. Shandong Gold Group Co., Ltd, Laizhou 261442, China

Abstract: The rapid rise of atmospheric oxygen content in Paleoproterozoic (i.e., Great Oxidation Event) resulted in many dramatic changes in earth's surficial system, including the 2.3~2.06Ma positive δ¹³C_carb excursion, termed Lomagundi-Jatuli Event, which has been a focus of international geological and geochemical studies. The Lugezhuang Formation, lower portion of the Jingshan Group, is a high amphibolite to granulite facies metamorphic carbonate association. Herewith we report the results obtained from isotope and element geochemical study of 26 marble and 7 dolomitic quartzite samples from the Lugzhuang Formation at Shanhou village. The δ¹³C_carb values of the samples range from -0.8‰ to +3.4‰ (V-PDB), with an average of +1.4‰. The δ¹³C_carb averages of marbles and dolomitic quartzites are +1.6‰ and +0.9‰, respectively. The δ¹⁸O_carb values range from 9.3‰ to 20.2‰ (V-SMOW), with an average of 14.0‰. The δ¹⁸O_carb and δ¹³C_carb values positively correlate, and 13 samples show Mn/Sr ratios of >6. This indicates that the δ¹³C_carb values of carbonates have been intensely modified by metamorphism, possibly causing a decrease of δ¹³C_carb values by 3‰ to 5‰. Hence, the δ¹³C_carb values originally sedimentary carbonates are possibly up to +8.4‰, showing positive δ¹³C_carb excursion characterizing the global Lomagundi-Jatuli Event.

Key words: Lomagundi-Jatuli Event Great Oxidation Event δ¹³C_carb anomaly Lugezhuang Formation of Jingshan Group Jiaobei Terrain

科学家共识地球形成初期大气圈是缺氧还原性的，现今大气圈富氧，氧气含量达21%(体积)。富氧大气圈是人类生存的必需条件，其形成机制、过程和时间自然备受关注。1980年之前，科学家普遍认为地球大气圈和水圈的氧化过程是缓慢的、渐变的，至少始于3.8Ga，主要发生在2.6~1.9Ga期间(Cloud, 1968; Schidlowski et al, 1975; Frakes, 1979)；1980年之后，受白垩纪末期恐龙灭绝事件研究的影响，陈衍景及合作者认识到古元古代水-气系统充氧事件及相关事件的突然性、短时性和系统性，2.3Ga前后突然大量出现红层、冰碛岩、含叠层石的厚层碳酸盐、菱镁矿、磷块岩、蒸发岩(石膏、硼酸盐等)、苏必利尔湖型条带状铁建造、石墨矿床以及沉积物负铕异常，提出了2.3Ga时地球表生环境由还原性突变为氧化性的观点(Chen, 1988；Chen et al, 1988; 陈衍景, 1987, 1990; 季海章和陈衍景，1990；陈衍景等, 1991, 1996)。

国际地层委员会前寒武纪分会高度重视2.3Ga环境突变事件及其分期意义，提出了增设2.3Ga为成铁纪(Siderian)与层侵纪(Rhyacian)界线的前寒武纪分期方案，获得了国际地球科学联合会的赞同(Cowie et al, 1989；陈衍景等，1994)。新界线增设引发了2.3Ga地质事件性质研究的热潮，学者们成功尝试了多种元素和C、S、N、Mo等同位素方法(Schidlowski, 1988; Baker and Fallick, 1989a, b；Karhu, 1993; Karhu and Holland, 1996；Bekker et al, 2004；Anbar et al, 2007; Guo et al, 2009)，更广泛、深刻、精确地揭示或证实了2.3Ga表生环境突变事件(Chen and Zhao, 1997; Melezhik et al, 1999；Pavlov and Kasting, 2002；Huston and Logan, 2004；唐国军等，2004；Tang et al, 2011, 2013a, b, 2016; Martin et al, 2013；Tang and Chen, 2013；Young, 2013, 2019；Chen and Tang, 2016; Condie, 2016; 陈衍景和汤好书，2018；汤好书和陈衍景, 2018a, b)。据Karhu and Holland (1996)估算，在2.4~2.2Ga期间(即2.3Ga前后)，大气圈自由氧含量从 < 10^-13PAL增至15%PAL (PAL=Present Atmosphere Level)，足见充氧量之大、速度之快。因此，Holland (2002)引入大氧化事件(Great Oxidation Event: GOE)的概念，专指2.3Ga左右大气圈成分由缺氧到富氧的变化(The change in the composition of the atmosphere from anoxygenic to oxygenated ca. 2.3Ga)。

大氧化事件概念的提出进一步激发了学者们对古元古代表生环境演化研究的兴趣，并提出了多种观点或模式。例如，Kump et al (2011)提出了脉动式氧化的观点，即：(1)2500Ma左右大气氧含量达到10^-8~10^-5PAL，能使黄铁矿氧化；(2)~2400Ma时或2320Ma前，非质量硫同位素分馏现象突然消失，表明大气圈氧含量急剧升高(Bekker et al, 2004；Kump, 2008；Guo et al, 2009)；(3)~2000Ma时大气氧含量达~1%PAL。再如，Murakami et al (2011)基于古土壤研究提出了古元古代大气氧含量逐渐上升的观点，认为2100~2000Ma期间快速上升至~1%PAL。最近，我国学者提出了大氧化事件的两阶段氧化模式，即2.5~2.3Ga水圈氧化和2.3~2.05Ga大气圈氧化(Tang and Chen, 2013; 陈威宇和陈衍景, 2018; 陈威宇等, 2018; 汤好书和陈衍景, 2018a, b；Chen et al, 2019)，得到了多数学者的赞同(Young, 2013, 2019；Condie, 2016; 翟明国等, 2018)。

在揭示2.3Ga环境突变或大氧化事件的研究中，碳酸盐碳同位素(δ¹³C_carb)发挥了关键作用。Schidlowski et al(1975, 1976)最早发现俄罗斯Karelia的Jatulian碳酸盐岩(2.65~1.95Ga)和津巴布韦Lomagundi省白云岩(~2.07Ga)的δ¹³C_carb分别高达8.8‰和13.6‰，提出大量藻类发育和堆埋，导致CO₂富集¹³C，沉积碳酸盐碳同位素正异常。Schidlowski et al(1975, 1976)的开创性研究直到1989年才引起重视，各大陆相继发现了古元古代δ¹³C_carb正异常现象(陈威宇和陈衍景，2018及其引文)，显示了事件的全球性(图 1)。由于该δ¹³C_carb正异常事件最初发现于Lomagundi群和Jatuli群，故被称为Lomagundi或Jatuli事件，或者称为Lomagundi-Jatuli事件(简称LJE)。LJE发生时间的厘定是个难题(Melezhik et al, 2013)，因为记录δ¹³C_carb正异常岩石都是沉积岩，缺乏合适的定年矿物，往往只能通过其它地质体定年结果限定沉积作用的最大和最小年龄，得到的年龄范围通常跨度较大。Fennoscandian地盾记录LJE事件的碳酸盐地层沉积于2.2~2.06Ga之间(Karhu, 1993; Karhu and Holland, 1996)；南美Marquette Range超群(Bekker et al, 2006；Vallini et al, 2006)、南非Transvaal超群(Hannah et al, 2004; Frauenstein et al, 2009)、北美Huronian超群和芬兰Kuusamo带Sericite Schist组(Melezhik et al, 2013)的年龄数据显示LJE始于2.31~2.22Ga(Melezhik et al, 2013)；Fennoscandian地盾Kuusamo带的结束时间是~2080Ma，Pechenga带和Imandra-Varzuga带则为~2058Ma(Melezhik et al, 2007)。Martin et al (2013)统计分析了最近20年获得的古元古代地层U-Pb和Re-Os年龄, 确认LJE发生在2306±9Ma至2057±1Ma之间，即2.3~2.06Ga。

图 1 Lomagundi-Jatuli事件记录的碳酸盐岩剖面分布示意图(据Tang et al, 2011) Fig. 1 Distribution of the carbonates recording Lomagundi-Jatuli Event in the world (after Tang et al, 2011)

华北克拉通是全球最古老陆块之一，广泛发育了古元古代变质沉积建造，理应经历并记录了GOE、LJE等事件，但长期缺乏这些事件记录的报道，甚至缺乏专门研究。长期缺乏研究的主要原因是我国2.06Ga之前的地层变质较深，只有辽河群、滹沱群等少数地层单元为绿片岩相，其余为角闪岩相或麻粒岩相。在如此强烈变质改造的情况下，GOE、LJE等事件的标志性地质地球化学记录还能否被保存？即使保存下来，又能否被揭示出来？即使被揭示出来，又能否被学者们接受或认可？困难重重，使不少学者望而却步。一些关注GOE、LJE等事件的我国学者，也往往绕开华北克拉通而开展非洲、澳州或美洲大陆低变质或未变质古元古代沉积地层的研究，获得了容易被认可的GOE、LJE等的证据(Guo et al, 2009; Luo et al, 2016)。但是，华北克拉通GOE、LJE研究仍然被搁置，也仍然无法回答华北克拉通是否存在GOE、LJE记录等问题。显然，直接开展华北克拉通中高级变质地层的地质地球化学研究或探索，即使获得的结果或认识不如未变质区，也是回答上述问题，提升华北克拉通研究程度的必由之路。

根据华北克拉通石墨矿床的碳同位素研究资料，陈衍景等(2000)首先讨论了华北克拉通记录LJE的可能性，特别指出了胶东荆山群存在碳酸盐碳同位素正异常的信息，并推测荆山群原始沉积碳酸盐的δ¹³C_carb可能高于4.2‰。为检验这一推断，并探讨荆山群对古元古代大氧化事件及成矿大爆发事件的记录，我们对荆山群开展了专门研究，初步认为荆山群下部禄格庄组和野头组发育大量厚层大理岩，上部陡崖组大量发育石墨矿床，它们可能对应着记录了GOE的LJE和Shunga事件。本文旨在报道山后地区荆山群禄格庄组大理岩碳氧同位素及元素地球化学研究结果，探讨其对LJE的记录。

1 区域地质概况

胶北地体位于华北克拉通东部，北临渤海，西侧以郯庐断裂为界与鲁西地块相邻，东南端以烟台-五莲断裂为界与苏鲁超高压变质带相邻(杨进辉，2000；唐俊等，2004；刘平华等, 2011, 2013；谢士稳等，2014)。胶北地体主要出露新太古代胶东群、古元古代荆山群和粉子山群(孔庆友，2006)，以及大量中生代花岗岩(图 2)。

图 2 胶北地体区域地质图(据谢士稳等，2014修改) Fig. 2 Geological map of Jiaobei Terrain (modified after Xie et al, 2014)

荆山群是胶东半岛石墨矿、大理岩矿的主要赋矿层位(陈衍景等, 2000)，主要分布在胶北地体的平度仙山、明村，莱西南墅，莱阳荆山、旌旗山，牟平光山、祥山，海阳晶山、乳山午极等地。主体为一套碎屑岩夹碳酸盐岩的岩石组合，原岩恢复为一套多旋回的富铝泥沙质碎屑岩-钙镁质碳酸盐建造，间夹基性火山岩，是陆源浅海沉积产物，经历了高角闪岩相-麻粒岩相变质(季海章和陈衍景，1990；季海章等，1990; 卢良兆等，1996；孔庆友，2006)，自下而上可以分为禄格庄组、野头组和陡崖组(图 3)。禄格庄组分为下部的安吉村段和上部的光山段，安吉村段以(石墨、石榴、夕线)黑云片麻岩为主，间夹(石墨、黑云、石榴、透辉)变粒岩；光山段则以(蛇纹石、白云石)大理岩为主，夹有硅质大理岩、蛇纹石金云橄榄大理岩、角闪变粒岩和透辉变粒岩，偶夹有斜长角闪岩。野头组包括下部祥山段，以灰绿色透辉岩、透辉变粒岩为主，夹有斜长角闪岩、黑云斜长片麻岩以及长石石英岩；上部定国寺段，主要以蛇纹石大理岩、透辉石大理岩、方柱透闪大理岩等各类大理岩为主。陡崖组下部徐村段是石墨矿的主要含矿层位，主要岩性为石墨黑云变粒岩、石墨透辉变粒岩、透辉岩、黑云变粒岩、(石墨)黑云斜长片麻岩、石榴(夕线)二云(黑云)片岩，夹透闪石英岩、二云石英片岩等；上部水桃林段以石榴二云片岩、石榴黑云片岩为主，夹白云石英片岩等(张增奇和刘明渭，1996)。

图 3 荆山群地层柱(据李洪奎等，2013) Fig. 3 Stratigraphic units of the Jingshan Group (modified after Li et al, 2013)

前人针对荆山群地层做了许多年代学测试工作。Wan et al (2006)用SHRIMP锆石U-Pb年龄测试方法研究了烟台南部荆山群二云母夕线石榴片麻岩，获得最年轻的碎屑锆石²⁰⁷Pb/²⁰⁶Pb年龄为2175±16Ma。刘平华等(2011)获得莱西北墅南部夕线石榴斜长片麻岩的最年轻碎屑锆石²⁰⁷Pb/²⁰⁶Pb年龄为2100±10Ma。Tam et al (2011)对莱西南墅地区的麻粒岩和片麻岩碎屑锆石进行SHRIMP U-Pb年龄测试，获得最年轻的碎屑锆石年龄为2320±16Ma。Wan et al (2006)得到的荆山群二云母夕线石榴片麻岩的8颗变质锆石年龄在1859±11Ma~1902±8Ma之间，加权平均年龄为1882±12Ma。刘平华等(2011)获得栖霞小庄铺荆山群孔兹岩26颗变质锆石年龄为1847±8Ma~1879±5Ma，加权平均年龄为1868±3Ma。北墅夕线石榴斜长片麻岩的变质锆石年龄分为1853±8Ma~1922±12Ma和1833±9Ma~1844±11Ma两组，而由于最大年龄数据1922±12Ma的协和度小于90%，故将第二大年龄1901±20Ma作为这一组变质锆石的最大年龄。Tam et al (2011)测得栖霞太平庄麻粒岩21颗变质锆石年龄为1815±14Ma~1875±19Ma，加权平均年龄为1837±8Ma；莱西南墅陡崖组含石墨石榴石夕线石片麻岩的20颗变质锆石年龄为1790±19Ma~1853±27Ma，加权平均年龄1821±18Ma；莱西常家屯黑云片麻岩22颗变质锆石年龄1801±19Ma~1874±20Ma，上交点年龄为1836±68Ma；莱西南墅基性麻粒岩的4颗变质锆石年龄为1855±90Ma~2012±53Ma，加权平均年龄为1956±41Ma，但4颗锆石中的3颗年龄数据误差大于50Ma，可信度较低；平度双庙麻粒岩(08JB01-1-6)19颗变质锆石年龄为1633±150Ma~2125±130Ma，加权平均年龄为1884±24Ma。然而，双庙麻粒岩锆石年龄误差较大，最大可达230Ma，仅有7颗锆石年龄误差小于50Ma，且无协和度报道，本文仅做参考。

董春艳等(2010)对莱阳吕格庄地区侵入荆山群的变质闪长岩进行SHRIMP U-Pb年代学测试，获得14颗岩浆锆石年龄在1829±14Ma~1895±21Ma之间，加权平均年龄为1852±9Ma；莱西南墅变质辉长岩岩浆锆石年龄在1821±41Ma~1926±23Ma之间，11颗锆石的加权平均年龄为1865±11Ma。

荆山群碎屑锆石年龄分布较广，最大碎屑锆石年龄为2936±17Ma (Wan et al, 2006)，最小碎屑锆石年龄是2100±10Ma (刘平华等，2011)，说明物质来源复杂。变质锆石年龄范围相对集中，主要为1775~1875Ma，可能是较强变质作用时间，较可信的最大变质锆石年龄为1902±8Ma (Wan et al, 2006)。侵入荆山群的岩浆岩侵位年龄为~1.86Ga(董春艳等，2010)。本文将最年轻碎屑锆石年龄作为地层的最大沉积年龄，最大变质锆石年龄作为地层的最小沉积年龄，得到荆山群沉积年龄为2.10~1.90Ga。

2 样品和分析方法 2.1 样品及采样位置

本文研究的33件大理岩和白云质石英岩样品采自胶北地体北部的古元古代荆山群，位于平度市东北的南墅镇山后村附近(图 2、图 4)。山后村与南墅镇之间分布多个石墨矿床，以南墅石墨矿床最为知名，它们赋存于荆山群陡崖组地层。本文样品则采自位于荆山群下部的禄格庄组光山段，所有样品采自深度为725.2m的钻孔(ZK 1604)。

图 4 山后大理岩矿矿区地质图(山东正元地质资源勘查有限责任公司，2011^①) Fig. 4 Simplified geological map of the Shanhou marble deposit

① 山东正元地质资源勘查有限责任公司.2011.山东省平度市南墅镇山后大理岩矿床普查报告

33件样品分为大理岩和白云质石英岩两类。其中，26件为纯度较高的白云质大理岩，手标本呈乳白色(图 5a, c, d)，镜下白云石自形-半自形(图 5f-h)，含有少量石英颗粒(图 5g, i-l)；7件白云质石英岩(SiO₂>30%)样品呈灰白色(图 5b, c, e)，镜下石英颗粒呈不等粒状，碳酸盐矿物主要为白云石，次为方解石，方解石多呈网脉充填在石英颗粒间隙(图 5m-o)。详细描述见表 1。

图 5 荆山群禄格庄组样品地质特征 (a)大理岩, 样品SH07；(b)含钙石英岩, 样品SH26，靠近含碳角砾岩；(c)变粒岩和含钙石英岩, 样品SH30；(d)白色块状大理岩, 样品SH07；(e)石英岩, SH30；(f)粗粒白云石大理岩，发育石英网脉，SH09；(g)白云石大理岩，SH21；(h)大理岩遭受方解石化，SH03; (i)白云石和石英颗粒，SH01；(j)早期细粒白云石(DoI)大理岩，可见晚期重结晶粗粒白云石(DoII)和石英，SH05；(k、l)硅化蚀变, SH25, SH19; (m)石英岩，石英颗粒被方解石胶结，样品SH40；(n)角砾状构造石英岩, SH26; (o)硅化蚀变，白云母蚀变, SH26. Do-白云石; Qtz-石英; Cc-方解石; Ms-白云母 Fig. 5 Geological characteristics of samples from the Lugezhuang Formation, Jingshan Group

表 1 荆山群禄格庄组样品地质特征 Table 1 Geological characteristics of samples from the Lugezhuang Formation, Jingshan Group

样品号	岩性特征	原岩岩性	孔深(m)	样品地质特征
SH01	大理岩	白云岩	90	变晶结构，块状构造；碳酸盐矿物85%，粒径0.01~1mm；石英15%，粒径多为0.1~1mm；发育宽约0.5mm石英脉(图 5i)
SH02	大理岩	白云岩	98.2	不等粒结构；碳酸盐矿物90%，粒径多为0.1~0.5mm；石英10%，粒径多小于0.5mm；硅化较强
SH03	大理岩	白云岩	103.5	不等粒结构；碳酸盐矿物95%，粒径0.05~0.5mm；石英5%，粒径小于0.1mm，多沿裂隙发育，伴有杂质沿裂隙混入(图 5h)
SH04	蚀变大理岩	泥质白云岩	105	碳酸盐矿物约35%，石英50%，透辉石约10%，白云母3%，赤铁矿化硫化物约2%；绢云母化和硅化强烈，发育碳酸盐细脉
SH05	大理岩	白云岩	108	变晶结构；碳酸盐矿物90%，石英9%，粒径都集中在 < 0.1mm和0.2~0.6mm两个范围(图 5j)；硫化物1%，粒径 < 0.1mm，发育石英脉
SH06	大理岩	白云岩	113	不等粒结构；碳酸盐矿物95%，粒径集中于1~4mm；石英5%，粒径约0.05~0.1mm，属填隙式石英脉，伴随弱硅化蚀变
SH07	大理岩	白云岩	115	不等粒结构；碳酸盐矿物95%，粒径0.5~5mm；石英5%，粒径 < 0.1mm；沿裂隙发育石英脉(图 5d)
SH08	钙硅质片麻岩	钙质泥岩	118	石英55%，粒径0.01~0.5mm；碳酸盐矿物44%，粒径 < 0.1mm；绢云母5%，硫化物1%；碳酸盐矿物发生硅化和绢云母化蚀变
SH09	石英大理岩	硅质白云岩	120.2	不等粒结构；碳酸盐矿物90%，粒径> 0.5mm；石英9%，少数粒径大于0.5mm，多数小于0.1mm；绢云母约1%；强硅化，弱绢云母化，石英沿裂隙充填(图 5f)
SH10	钙硅质片麻岩	钙硅质泥岩	122	石英55%，颗粒细小，多为蚀变产物。碳酸盐矿物43%，粒径多大于0.5mm，遭受强硅化蚀变；硫化物1%，绢云母1%
SH11	钙硅质片麻岩	钙硅质泥岩	123	石英42%，粒径较小；碳酸盐矿物56%，粒径0.1~1mm；白云母1%，黄铁矿1%；强硅化蚀变，沿裂隙有碳酸盐和少量石英充填
SH12	大理岩	白云岩	252.3	碳酸盐矿物95%，粒径大多小于0.1mm，部分粒径大于0.5mm；石英4%，多沿裂隙发育；少量硫化物(1%)
SH13	大理岩	白云岩	256.4	碳酸盐矿物95%，粒径0.2~1mm；石英5%，粒径0.2~0.5mm，沿裂隙发育石英脉
SH14	大理岩	白云岩	262	碳酸盐矿物95%，粒径0.1~1mm；石英5%，粒径0.05~0.5mm
SH15	硅化大理岩	白云岩	266	不等粒结构；碳酸盐矿物95%，粒径0.1~1mm；石英5%，粒径0.05~0.5mm；沿裂隙发育石英脉
SH16	硅化大理岩	白云岩	268	不等粒结构；碳酸盐矿物90%，粒径集中在 < 0.01mm和0.1~0.5mm两个范围；石英10%，多为硅化产物，粒径0.01~0.1mm
SH18	硅化大理岩	硅质白云岩	463	不等粒结构；碳酸盐矿物90%，石英10%，石英多为后期重结晶形成
SH19	硅化大理岩	硅质白云岩	486	不等粒结构；碳酸盐矿物78%，粒径差异大，多在0.1~1mm之间；石英20%，多为硅化产物，粒径小于0.1mm(图 5l)；白云母2%
SH20	硅化大理岩	硅质白云岩	498	不等粒结构；碳酸盐矿物90%，粒径差异大；石英10%，多为硅化产物，粒径小于0.1mm
SH21	硅化大理岩	硅质白云岩	505	碳酸盐矿物90%，粒径在1~5mm和0.01~0.2mm两个范围；石英10%，粒径0.05~0.5mm；颗粒较小的石英和碳酸盐矿物分布在颗粒较大的碳酸盐岩矿物边部，沿裂隙有石英和碳酸盐矿物充填(图 5g)
SH22	硅化大理岩	硅质白云岩	509	粗粒结构，块状构造；碳酸盐矿物96%，粒径> 5mm；石英4%，颗粒细小，主要沿裂隙充填
SH23	大硅化理岩	白云岩	515	不等粒结构；碳酸盐矿物85%，明显硅化；石英15%，部分颗粒约0.01mm，部分颗粒在0.1~0.3mm之间，水岩作用产物；风化较强
SH24	硅化大理岩	硅质白云岩	518	不等粒结构；碳酸盐矿物92%，大颗粒边部有重结晶的小颗粒碳酸盐和石英；石英8%，多沿裂隙呈脉状分布
SH25	大理岩	白云岩	525	碳酸盐矿物92%，粒径集中在0.1~1mm之间；石英8%，颗粒细小，多小于0.05mm，多沿裂隙发育(图 5k)
SH26	石英片麻岩	钙硅质泥岩	540	角砾状构造；石英70%，粒径达1mm，基质石英颗粒细小；碳酸盐矿物27%，粒径多小于0.2mm；白云母2%，少量硫化物1%。采自角砾岩附近(图 5n，o)
SH27	硅化大理岩	白云岩	548.5	不等粒结构；碳酸盐矿物90%，粒径差异大；石英10%，粒径多小于0.05mm；硅化明显
SH28	石英大理岩	硅质白云岩	581	不等粒结构；碳酸盐矿物85%，粒径多在0.1~0.6mm；石英15%，粒径多在0.1~0.8mm；水岩作用现象明显
SH29	大理岩	白云岩	587	等粒变晶结构；碳酸盐矿物95%，粒径多在0.3~0.6mm之间；石英5%，粒径约0.2~0.4mm
SH30	钙质石英岩	钙质硅质岩	592	不等粒结构(图 5e)；石英65%，粒径多小于0.03mm；碳酸盐矿物35%，粒径差异较大；硅化强烈
SH31	硅化大理岩	硅质白云岩	596.6	碳酸盐矿物90%，粒径多在0.2~0.5mm；石英10%，粒径在0.2~0.5mm和小于0.1mm两个范围；细粒石英为硅化产物
SH32	硅化大理岩	硅质白云岩	604	碳酸盐矿物75%，粒径多在0.02~0.5mm之间；石英24%，粒径多在0.01~0.5mm之间；硫化物约1%；蚀变明显
SH33	大理岩	硅质白云岩	609	碳酸盐矿物92%，粒径多小于0.5mm；石英8%，粒径多小于0.5mm；弱硅化，弱风化
SH40	钙质石英岩	钙质硅质岩	662	不等粒结构，石英75%，粒径差异较大；碳酸盐矿物24%，与石英共生或沿裂隙充填(图 5m)；白云母1%

表 1 荆山群禄格庄组样品地质特征 Table 1 Geological characteristics of samples from the Lugezhuang Formation, Jingshan Group

2.2 分析方法

碳-氧同位素分析在中国科学院地球化学研究所环境地球化学国家重点实验室完成。粉末样品首先通过德国Thermo Fisher Scientific的GasbenchⅡ在氦气流中冲洗8min，然后在真空系统中使用100%的磷酸与粉末样品在72℃恒温水浴中恒温反应4h以上，收集碳酸盐矿物释放出的CO₂用于质谱分析。使用Gasbench Ⅱ连续流进样，经Thermo Fisher Scientific MAT 253型气体质谱分析仪分析得出样品的δ¹³C值和δ¹⁸O值，碳同位素结果均采用V-PDB标准，氧同位素采用V-SMOW标准，测试精度为0.2%。

全岩主量元素测试在澳实分析检测(广州)有限公司完成。首先把包含硝酸锂在内的助溶剂加入待测样品中，待充分混合后高温熔融，把熔融物倒入铂金模子形成扁平玻璃片后，采用PANalytical Axios X射线荧光光谱仪进行分析。

全岩微量元素测试在中科院地球化学研究所矿床地球化学国家重点实验室完成。采用的仪器是四级杆型电感耦合等离子质谱仪，仪器型号为ELAN DRC-e ICP-MS，样品分析数据的相对精度优于10%，详细测试方法见Qi et al (2000)。

3 测试结果

荆山群禄格庄组33件样品碳-氧同位素分析结果，主量及Sr元素地球化学测试结果均列于表 2。大理岩的碳-氧同位素值在地层上的分布见图 6，频率分直方图见图 7。

表 2 荆山群禄格庄组样品主要氧化物(wt%)、Sr(×10^-6)及δ¹³C_V-PDB (‰)和δ¹⁸O_SMOW (‰)组成 Table 2 Composition of main oxides (wt%), Sr (×10^-6), δ¹³C_V-PDB (‰) and δ¹⁸O_SMOW (‰) of samples from the Lugezhuang Formation, Jiangshan Group

样品号	SiO₂	Al₂O₃	Fe₂O₃^T	MgO	CaO	Na₂O	K₂O	MnO	P₂O₅	TiO₂	LOI	Total	Sr	δ¹³C	δ¹⁸O	Mg/Ca	Mn/Sr
大理岩(n=26)
SH01	4.57	0.30	1.27	18.84	28.97	0.00	0.06	0.05	0.02	0.02	44.46	98.55	64.12	2.2	20.0	0.91	9.27
SH02	13.54	1.07	2.91	15.45	26.78	0.00	0.14	0.17	0.04	0.04	39.64	99.77	137.3	1.3	14.7	0.81	15.15
SH03	4.18	0.23	0.77	18.79	29.33	0.00	0.02	0.11	0.01	0.02	45.13	98.59	62.62	2.5	17.5	0.90	21.73
SH05	7.64	0.63	3.89	16.55	29.63	0.00	0.11	0.09	0.06	0.02	42.16	100.78	227.4	1.6	14.0	0.78	5.12
SH06	6.35	0.47	0.89	18.76	28.60	0.00	0.01	0.03	0.06	0.01	44.24	99.43	57.23	1.9	18.8	0.92	6.96
SH07	3.10	0.15	0.85	19.77	30.00	0.00	0.02	0.03	0.04	0.01	45.62	99.59	55.67	2.4	20.2	0.92	6.89
SH09	15.81	0.11	2.80	15.64	26.17	0.00	0.02	0.08	0.01	0.01	38.79	99.44	134.5	1.6	13.8	0.84	6.96
SH12	4.28	0.74	4.01	16.26	29.11	0.20	0.06	0.15	0.02	0.03	45.56	100.42	95.60	2.9	18.6	0.78	19.32
SH13	1.82	0.38	0.96	19.57	31.72	0.00	0.05	0.16	0.02	0.02	44.12	98.81	180.6	-0.8	10.8	0.86	10.83
SH14	5.29	0.59	2.00	17.15	28.80	0.00	0.05	0.12	0.02	0.03	44.16	98.21	160.9	0.4	12.8	0.83	9.34
SH15	4.73	0.61	1.03	18.44	28.83	0.00	0.07	0.08	0.02	0.02	44.43	98.26	288.7	2.8	13.3	0.90	3.37
SH16	5.10	1.20	1.30	18.03	28.39	0.02	0.09	0.09	0.02	0.04	44.11	98.39	429.5	2.6	12.8	0.89	2.48
SH18	16.91	1.80	4.05	14.17	24.75	0.00	0.17	0.17	0.01	0.08	36.71	98.83	90.49	3.4	10.6	0.80	23.52
SH19	25.97	2.68	2.07	12.87	21.06	0.00	0.71	0.06	0.01	0.09	32.99	98.51	260.1	-0.1	9.8	0.86	2.84
SH20	13.67	0.26	0.85	16.57	26.30	0.00	0.07	0.15	0.01	0.02	40.64	98.53	117.6	2.2	17.1	0.88	15.67
SH21	21.30	0.30	1.02	15.48	23.69	0.00	0.07	0.03	0.02	0.01	36.74	98.65	125.2	1.6	15.3	0.92	3.11
SH22	8.02	0.24	0.90	18.30	27.93	0.00	0.06	0.03	0.01	0.01	43.10	98.60	81.96	2.1	18.2	0.92	4.85
SH23	4.87	0.20	0.72	19.12	28.97	0.00	0.06	0.03	0.01	0.01	44.65	98.65	89.84	2.2	18.9	0.92	4.50
SH24	8.23	0.06	0.85	18.44	28.17	0.00	0.02	0.05	0.02	0.01	43.31	99.15	122.3	1.8	15.9	0.92	5.51
SH25	5.16	0.19	0.71	19.23	28.70	0.00	0.07	0.03	0.01	0.01	44.53	98.64	91.83	1.6	16.9	0.94	4.11
SH27	19.62	0.29	1.40	16.56	25.36	0.19	0.07	0.04	0.01	0.01	37.77	101.32	247.8	0.4	12.2	0.91	1.82
SH28	28.63	0.26	2.42	13.03	21.22	0.00	0.04	0.07	0.03	0.01	32.85	98.56	471.7	0.5	9.8	0.86	1.72
SH29	6.21	0.35	1.30	17.94	28.76	0.00	0.10	0.04	0.01	0.01	44.02	98.75	257.7	1.6	16.7	0.87	2.13
SH31	13.60	0.15	1.72	16.49	26.35	0.00	0.05	0.05	0.01	0.02	40.13	98.57	215.1	0.7	13.0	0.88	2.72
SH32	21.97	0.09	1.94	15.03	23.67	0.00	0.03	0.05	0.01	0.01	36.44	99.23	298.4	0.1	9.9	0.89	1.89
SH33	10.47	0.09	2.28	17.28	27.14	0.00	0.03	0.05	0.01	0.01	41.89	99.24	185.9	1.2	14.1	0.89	3.64
平均值	10.81	0.52	1.73	17.07	27.25	0.02	0.09	0.08	0.02	0.02	41.47	99.06	175.0	1.6	14.8	0.88	7.52
均方差	6.38	0.40	0.85	1.59	2.09	0.03	0.06	0.04	0.01	0.01	3.23	0.60	87.53	0.8	2.8	0.04	4.98
含钙石英岩(n=7)
SH04	36.82	13.17	7.53	5.25	11.08	0.00	3.21	0.11	0.13	1.07	20.01	98.37	108.8	0.0	11.2	0.66	12.00
SH08	40.11	15.28	6.58	5.07	9.06	0.40	3.75	0.09	0.04	0.59	18.32	99.29	254.2	0.6	9.3	0.78	4.17
SH10	50.22	8.59	3.68	5.85	10.23	0.29	1.48	0.07	0.06	0.45	18.19	99.11	348.2	2.4	10.3	0.80	2.55
SH11	33.14	12.81	4.52	7.34	13.35	0.34	1.61	0.09	0.06	0.71	24.80	98.76	460.0	2.5	10.7	0.77	2.35
SH26	56.95	8.72	3.57	4.17	8.89	0.17	0.58	0.08	0.03	0.01	16.34	99.49	215.5	0.1	10.7	0.66	4.42
SH30	73.95	0.47	0.60	4.51	7.93	0.03	0.11	0.02	0.05	0.02	11.04	98.73	85.41	0.3	12.9	0.80	2.22
SH40	67.18	0.32	3.08	4.33	9.38	0.00	0.06	0.09	0.01	0.01	14.07	98.52	158.9	0.2	10.6	0.65	6.63
平均值	51.20	8.48	4.22	5.22	9.99	0.17	1.54	0.08	0.05	0.41	17.54	98.90	233.0	0.9	10.8	0.73	4.91
均方差	12.71	4.62	1.70	0.80	1.34	0.14	1.13	0.02	0.02	0.34	3.19	0.34	103.8	0.9	0.7	0.06	2.52

图 6 禄格庄组地层柱状图及δ¹³C_V-PDB (‰)和δ¹⁸O_SMOW (‰)组成 Fig. 6 Stratigraphic column, δ¹³C (‰) and δ¹⁸O (‰) values of samples from the Lugezhuang Formation

图 7 禄格庄组碳酸盐样品δ¹³C (a)和δ¹⁸O (b)分布直方图 Fig. 7 Histograms for δ¹³C (a) and δ¹⁸O (b) values of samples from the Lugezhuang Formation

26件白云石大理岩样品具有较稳定的CaO和MgO含量，分别为21.06%~31.72%和12.87%~19.77%，平均分别为27.25±2.09%和17.07±1.59%。烧失量(LOI)和MgO/CaO(mol)较稳定，分别为32.85%~45.62%和0.78~0.94，平均分别为41.47±3.23%和0.88±0.04。SiO₂含量在1.82%~28.63%之间变化，平均10.81±6.38%。相较而言，7件含钙石英岩CaO、MgO含量较低，分别为7.93%~13.35%和4.17%~7.34%，平均分别为9.99±1.34%和5.22±0.80%。相应烧失量较低，为11.04%~24.80%，平均为17.54±3.19%。SiO₂含量在33.14%~73.95%之间，变化范围较大，平均51.20±12.71%。随着SiO₂含量升高，MgO/CaO值逐渐下降，为0.65~0.80，平均0.73±0.06(表 2)。

样品Al₂O₃含量多 < 1%，26件大理岩样品的Al₂O₃含量在0.06%~2.68%之间，平均0.52±0.40%。而7件含钙石英岩样品Al₂O₃含量较高，在0.32%~15.28%之间，平均8.48±4.62%。Fe₂O₃^T含量变化范围在0.60%~7.53%之间，大理岩样品平均为1.73±0.85%，含钙石英岩样品平均为4.22±1.70%，含量有所增加。所有样品的Na₂O含量 < 0.5%，K₂O的含量范围变化稍大，在0.01%~3.75%之间。MnO含量范围在0.02%~0.17%(154×10^-6~1725×10^-6)，大理岩样品和石英岩样品MnO的平均值均为0.08%。此外，样品中还含有少量的P₂O₅(0.01%~0.13%)和TiO₂(0.01%~1.07%)。微量元素实验数据显示，样品的Sr含量在55.67×10^-6~471.7×10^-6，大理岩样品平均值为175.0±87.53×10^-6，含钙石英岩样品平均值为233.0±103.8×10^-6。大理岩Mn/Sr范围为1.72~23.52，平均7.52±4.98；含钙石英岩Mn/Sr在2.22~12.00之间，平均值为4.91±2.52。

大理岩样品的δ¹³C值在-0.78‰~3.42‰之间，平均1.57±0.80‰(V-PDB)；δ¹⁸O值相对较高、变化范围大，为9.76‰~20.23‰，平均14.84±2.77‰(V-SMOW)。含钙石英岩样品δ¹³C值为0.02‰~2.50‰之间，平均0.88±0.90‰(V-PDB)；δ¹⁸O值较低且变化范围小，为9.28‰~12.86‰，平均10.82±0.70‰(V-SMOW)(表 2)。

4 讨论 4.1 碳酸盐岩碳氧同位素地球化学行为

在各种碳储库中，碳酸盐岩碳同位素组成最高，δ¹³C值平均为0.5‰ (表 3)。同时，碳酸盐岩的氧同位素组成也较高，δ¹⁸O_carb值通常大于20‰。据此可以推断，碳酸盐岩与其它碳储库进行碳同位素交换时，容易释放¹³C，汲取¹²C，其δ¹³C降低，同时，与之平衡反应的其它碳储库δ¹³C升高。事实上，学者们已经共识碳酸盐岩沉积之后的成岩、变质和流体作用均导致δ¹³C_carb和δ¹⁸O_carb降低，且δ¹³C_carb和δ¹⁸O_carb值随变质程度升高而降低(Veizer, 1983; Banner and Hanson, 1990; Jacobsen and Kaufman, 1999; Veizer et al, 1999; 陈衍景等, 2000; Melezhik et al, 2001a, b, 2005；汤好书等, 2008, 2009)。

表 3 地球主要碳储库的δ¹³C_V-PDB (‰) Table 3 The δ¹³C values (PDB‰) of major carbon reservoirs

碳酸盐岩沉积之后的地质作用导致碳氧同位素降低的根本原因是：碳酸盐矿物的化学成分可表示为Mg、Ca、Fe、Mn的金属氧化物(MO)与CO₂的结合物，即MO(CO₂。根据软硬酸碱理论(参见戴安邦，1987；Chen and Zhao, 1997)，与O配合或结合时，C⁴⁺是比Mg²⁺、Ca²⁺、Fe²⁺、Mn²⁺更硬的酸，而¹⁸O是比¹⁶O更硬的碱，因此，在MO(CO₂中，¹⁸O倾向于分配到CO₂中。在碳酸盐岩经历成岩、变质作用而脱CO₂时，由于¹⁸O随CO₂逸失而导致剩余MO·CO₂的相对富集¹⁶O，即δ¹⁸O_carb降低。同样道理，¹³C比¹²C硬度大，更倾向于与¹⁸O结合，形成¹³C¹⁸O₂；因此，碳酸盐岩以任何形式脱CO₂时，都会造成残余碳酸盐岩δ¹³C_carb降低。然而，目前对δ¹³C_carb和δ¹⁸O_carb在后期地质作用中降低的幅度，尚缺乏定量估算公式或办法，只有一些统计性研究。

Schidlowski et al (1975)统计认为碳酸盐岩δ¹⁸O_carb值在重结晶或变质过程中降低约2‰~3‰，Veizer et al (1999)则认为碳酸盐岩δ¹⁸O_carb在成岩作用中降低~2‰。Bottinga (1969), Wada and Suzuki (1983), Schidlowski (1988)等研究表明，温度高于650℃的变质作用会使δ¹³C_carb降低幅度>3‰。就Fennoscandian地盾区而言，几乎没有经历蚀变的碳酸盐(白云岩和石灰岩)δ¹³C_carb为+5‰~+9.6‰(平均7.4±0.7‰)，经历了绿帘-角闪岩相变质的碳酸盐δ¹³C_carb值为+1‰到+5‰(Melezhik et al, 2013)。如果角闪岩相碳酸盐尚可保留初始沉积时¹³C富集特征的话，麻粒岩相碳酸盐则几乎不能保留¹³C富集信息，因其δ¹³C_carb值比角闪岩相降低5‰以上，δ¹⁸O_carb值降低幅度可达10‰以上，可从~25‰降低到~10‰，致使麻粒岩相碳酸盐平均δ¹⁸O_carb值一般不超过18‰(Baker and Fallick, 1989a, b; Melezhik et al, 2001a, b)。由上可见，变质等级越高，δ¹³C_carb和δ¹⁸O_carb值越低。

流体-岩石相互作用也会导致碳酸盐矿物δ¹³C和δ¹⁸O值下降(Melezhik et al, 1999; Kump et al, 2011)。例如，密西西比Burlington-Keokuk组白云岩较好记录了流体-岩石相互作用的历史(Banner and Hanson, 1990)，重结晶白云岩δ¹⁸O_carb值变化较大，δ¹³C_carb值变化较小，说明δ¹⁸O_carb相对于δ¹³C_carb更易发生变化。因此，Hudson (1977)建议δ¹⁸O_carb可作为岩石遭受流体作用的指示剂，冯伟民等(2003)和Aharon (2005)主张以δ¹⁸O_carb＝18‰作为判别碳酸盐岩发生热液蚀变作用的阀值，而Melezhik et al (2005)则主张δ¹⁸O_carb阀值为20‰。

LJE期间的碳酸盐碳、氧同位素值偏高。Schidlowski et al (1976)研究的LJE期间的δ¹⁸O_carb值平均为22.2±1.6‰，最大值达26.9‰。Veizer et al(1992a, b)认为早前寒武纪白云岩δ¹⁸O_carb较高，平均为26±2‰。Melezhik et al (1999)认为早前寒武纪白云岩δ¹⁸O_carb平均值在28‰左右，并将25‰~26‰作为岩石遭受蚀变或变质与否的阈值。

一些学者(Veizer, 1983; Valley, 1986; Banner, 1995；Jacobsen and Kaufman, 1999; Veizer et al, 1999; Ray, 2009)专门研究了碳酸盐岩δ¹⁸O较δ¹³C更易变化的原因。Banner and Hanson (1990)认为，在开放体系下，即使流体/岩石比值很低(F/R < 10)，碳酸盐矿物也能与富水流体之间达到氧同位素分馏平衡，但碳同位素分馏平衡则要在流体/岩石比值高达10³时才能实现，其原因是碳酸盐岩碳含量高，可以缓冲流体作用引起的δ¹³C_carb变化。

在流体作用下，碳酸盐岩δ¹³C_carb和δ¹⁸O_carb值降低方式可分为三类(汤好书等, 2008；Tang et al, 2011, 2013a)：(1) δ¹³C_carb和δ¹⁸O_carb同步降低(Veizer and Hoefs, 1976; Guerrera et al, 1997; Jacobsen and Kaufman, 1999)，发生在F/R比值较高的条件下；(2)近乎平行于δ¹⁸O_carb轴的极端变化趋势，发生在F/R比值较低的条件下(Banner and Hanson, 1990)，且流体富水贫碳质(如CO₂)；(3)近乎平行于δ¹³C_carb轴的极端变化趋势，是碳质流体作用的特定现象。

4.2 荆山群碳酸盐碳氧同位素组成

造山纪发生了全球性陆块汇聚、碰撞、造山作用，最终形成哥伦比亚超大陆。荆山群在造山纪哥伦比亚超大陆会聚过程中发生了强烈的变形变质，周喜文等(2001)用石榴子石-黑云母矿物对计算出南墅地区荆山群孔兹岩系的峰期变质温度在700~750℃之间，李凯月等(2018)通过石墨拉曼光谱研究估算张舍石墨矿所经历的最高变质温度在730~790℃之间。可见，胶北地体含石墨矿的孔兹岩系所经历的最高变质温度可能在700~800℃之间，变质等级达到高角闪岩相-麻粒岩相。那么，在如此强烈、复杂地质改造作用中，荆山群碳酸盐岩地球化学特征和碳氧同位素组成发生哪些改变呢？

4.2.1 禄格庄组碳酸盐岩流体作用及碳同位素变化

蚀变作用常导致碳酸盐岩Sr含量降低，Mn、Fe、Rb含量增加(Brand and Veizer, 1980; Derry et al, 1992; Melezhik et al, 2001a, b, 2008; Bekker and Kaufman, 2007)，Mn和Sr含量可以指示成岩作用的级别。Mn/Sr比值对碳酸盐岩的后期流体作用较为敏感(Kaufman et al, 1993; Bekker et al, 2008)。一般认为，Mn/Sr>6时，碳酸盐岩受后期流体作用影响较大(Kaufman and Knoll, 1995; Melezhik et al, 1999)。

在荆山群禄格庄组样品中，13件样品Mn/Sr比值超过6，但δ¹³C-Fe₂O₃^T、δ¹³C-Mn/Sr、δ¹⁸O-Mn/Sr之间的相关性不明显，R值分别为0.14、0.34、0.25(图 8a-c)，没有显示成岩作用对禄格庄组碳酸盐岩碳氧同位素的强烈影响。

图 8 禄格庄组样品元素协变图 Fig. 8 Correlations of elements in samples of the Lugezhuang Formation

早前寒武纪碳酸盐地层主要为白云岩，后期流体交代作用趋向于生成方解石(Melezhik et al, 2001a, b; 图 5j-l)。在图 8中，禄格庄组样品主要落入白云岩区域，Mg/Ca-Mn/Sr的相关系数R值为0.15，没有明显的相关性。δ¹³C-Mg/Ca、δ¹⁸O-Mg/Ca之间的相关系数R值分别为0.32和0.54(图 8e, f)，显示中等程度的正相关，表明白云岩的方解石化越强烈，δ¹³C和δ¹⁸O值越低。

Mg/Ca与SiO₂之间的相关系数为-0.71(图 8g)，显著负相关，说明硅化导致MgCO₃减少或者CaCO₃的相对增多(方解石化)。在图 8h中，SiO₂与LOI之间明显负相关，R=-0.98，原因之一是硅化导致岩石中CO₂等挥发分含量降低，原因之二是记录了原始沉积SiO₂与碳酸盐矿物之间的消长关系。从岩石薄片观测结果来看，石英交代、充填碳酸盐矿物的现象普遍、清楚，硅化证据充分。虽然无法排除SiO₂等组分是原始沉积形成，但排他性的岩相学依据不足够充分。9件样品SiO₂>25%，其δ¹⁸O值稳定在9.28‰~11.21‰，说明碳酸盐矿物与含SiO₂流体之间基本达到了氧同位素平衡(图 8i)；同时，这9件样品的δ¹³C值介于-0.06‰~2.50‰之间，明显低于其它研究样品。以上说明，流体作用使δ¹³C、δ¹⁸O值降低。

4.2.2 禄格庄组碳酸盐岩变质作用及氧同位素变化

禄格庄组33件样品δ¹⁸O值在9.28‰~20.23‰之间，即9.3‰~20.2‰之间，平均14.0‰。33件样品中，只有1件样品的δ¹⁸O值大于20‰，5件样品δ¹⁸O值大于18‰，表明沉积后经历了强烈的地质作用改造。如前所述，SiO₂与Mg/Ca、LOI呈负相关关系，指示样品发生了方解石化、硅化或钙硅酸盐化。图 8i中，SiO₂ < 30%时，即白云质大理岩，样品δ¹⁸O_carb值随着SiO₂含量增高而降低，直至δ¹⁸O_carb降至11‰左右；SiO₂>30%时，即白云质石英岩，δ¹⁸O_carb值稳定在11‰附近，不再随SiO₂含量变化而改变。这一现象表明，原始沉积白云岩可能经历了两次改造作用：(1)成岩或变质作用及其伴随的脱挥发份作用，造成δ¹⁸O_carb值降低和SiO₂含量相对增高；(2)以硅化为代表的交代作用，表现为SiO₂含量不断增高，而碳酸盐δ¹⁸O_carb值稳定在11‰附近，说明残余碳酸盐与流体之间的氧同位素分馏达到了平衡。

前述现象和认识与显微镜观察结果一致，显示样品后期硅化等热液作用较明显，主要依据是：(1)大理岩样品多由中-粗粒白云石组成，发育石英网脉(图 5f)；(2)虽然少部分石英呈现重结晶特点(图 5i)，但绝大部分为硅化形成(图 5j-l)；(3)强硅化样品，如白云质石英岩，多具有角砾状构造或碎裂结构，发育在角砾岩带附近(图 5b, n, o、图 6)；(4)山后大理岩矿区及附近广泛发育多种碎裂岩、糜棱岩(图 4)，并发育NNE或EW向含金石英脉。

4.3 荆山群碳同位素正异常对LJE的指示

在全球范围内，不同时代海相碳酸盐岩的δ¹³C均值为0.5±2.5‰。通常认为，δ¹³C_carb值明显高于或低于平均值的碳酸盐岩形成于特殊沉积环境(Melezhik et al, 1999, 2013; Sanyal et al, 2009)；当δ¹³C_carb值大于3‰时，即可被视为大于普通碳酸盐岩的异常沉积(Karhu, 1993)。在LJE期间，碳酸盐岩δ¹³C_carb值普遍大于4‰ (Melezhik et al, 2013)，集中在5‰~18‰之间，甚至出现δ¹³C_carb值达到28‰极端案例(Bekker et al, 2003)。Melezhik et al (1999)认为LJE由3或4个δ¹³C_carb值正漂移旋回组成，每个旋回的δ¹³C_carb值都正漂移至5‰以上，之后再回到0‰；在多旋回的循环往复中，δ¹³C_carb值始终大于0‰，充分证明当时海洋碳储库富集¹³C。

近年，一些学者探索性地开展了华北克拉通古元古代变质地层碳氧同位素研究，陆续发现多处古元古代碳酸盐岩碳同位素正向漂移现象。例如，辽北辽河群关门山组(2.3~1.9Ga)42件白云岩δ¹³C_carb值为4.6‰~5.9‰，平均为5.3‰(汤好书等，2008)。辽东辽河群大石桥菱镁矿区(2.3~1.85Ga)大石桥组白云岩δ¹³C_carb值为0.6‰~1.4‰，平均为1.2‰，大石桥组原始沉积白云岩的δ¹³C_carb值应在4.2‰左右(汤好书等，2009)。河南登封嵩山群(2.51~1.95Ga)碳酸盐岩δ¹³C_carb值为-1.0‰~4.2‰，平均2.5‰，具有碳同位素正异常现象(Lai et al, 2012)。山西五台地区滹沱群大石岭组175件δ¹³C_carb数据普遍为正值，变化于0.6‰~3.5‰，具有δ¹³C_carb正漂移特征(陈威宇和陈衍景, 2018)。考虑到高角闪岩相-麻粒岩相变质作用可使δ¹³C_carb值降低3‰以上(Wada and Suzuki, 1983)，陈衍景等(2000)提出胶东南墅石墨矿床院后矿区荆山群大理岩的原岩碳酸盐δ¹³C_carb值应大于4.5‰，因为兰心俨(1981)报道的2件大理岩δ¹³C_carb值分别为0.8‰和1.5‰，已经显示了δ¹³C_carb正异常。

本研究获得荆山群禄格庄组33件大理岩和石英岩之碳酸盐胶结物样品的δ¹³C_carb值，变化于-0.8‰~3.4‰之间，平均为1.4±0.9‰。在33件样品δ¹³C_carb值中，只有2件为负值，分别为-0.78‰和-0.06‰，其余为正值，最大为3.4‰，显示了正异常特征。图 9中，禄格庄组碳酸盐碳、氧同位素之间的相关系数(R)为0.55，其δ¹⁸O变化范围(9.3‰~20.2%)大于δ¹³C变化范围(-0.8‰~3.4‰)，δ¹⁸O-δ¹³C散点图略显平行于δ¹⁸O轴的趋势。禄格庄组部分样品与世界LJE期间未变质或浅变质海相碳酸盐相分布范围一致，多数样品类似于角闪岩相-麻粒岩相变质的~2.15Ga的Lofoten-Vesteralen大理岩特征(Baker and Fallick, 1989a, b; Martin et al, 2013)。也就是说，虽然总体属于麻粒岩相的荆山群变质程度高于同时期的Lofoten-Vesteralen大理岩，但其碳氧同位素特征仍像后者一样保留了LJE事件特征，甚至个别样品的碳氧同位素特征竟然像未变质的LJE时期碳酸盐岩，δ¹³C_carb正异常清楚。事实上，高级变质过程的脱水、脱挥发份、重结晶作用及流体作用导致禄格庄组碳酸盐碳氧同位素大幅度降低，部分样品达到了石英、碳酸盐、流体之间的氧同位素分馏平衡，δ¹⁸O_carb稳定在11‰附近。考虑到δ¹⁸O_carb降低如此强烈，可以推断禄格庄组大理岩δ¹³C_carb值至少低于其原始沉积碳酸盐3‰~5‰。现今禄格庄大理岩δ¹³C_carb值平均为1.4‰，初始沉积碳酸盐δ¹³C_carb平均值应不低于4.4‰，可能为6.4‰，最大值甚至达8.4‰以上，较好地记录了LJE事件。

图 9 荆山群禄格庄组碳酸盐δ¹³C-δ¹⁸O同位素与全球Lomagundi海相碳酸盐对比(底图据Tang et al, 2013a) Fig. 9 δ¹³C vs. δ¹⁸O of the Lugezhuang Formation, Jingshan Group compared with those of the 2.31~2.06Ga Lomagundi marine carbonates (base map and data source see Tang et al, 2013a and references therein for details)

5 结论

(1) 胶北地体荆山群禄格庄组发育大量碳酸盐地层，其33件样品δ¹³C_carb值为-0.8‰~3.4‰之间。其中，大理岩样品δ¹³C_carb值平均为1.6±0.8‰，白云质石英岩δ¹³C_carb值平均为0.9±0.9‰。碳氧同位素值之间具有很好的正相关性，部分样品Mn/Sr比值大于6。禄格庄组碳酸盐经历了重结晶-变质作用以及硅化作用，造成样品去碳酸盐化强烈。

(2) 高角闪岩相-麻粒岩相变质作用和流体交代作用是影响碳氧同位素值的主要因素，导致使禄格庄组大理岩碳同位素值至少下降了3‰~5‰，原始沉积碳酸盐δ¹³C平均值应该不低于4.4‰，正异常明显，记录了全球性的LJE事件。

致谢样品采集过程中得到山东黄金集团有限公司史先洪工程师的热情帮助；实验测试工作得到了中国科学院地球化学研究所矿床地球化学国家重点实验室的大力支持，微量元素由该实验室胡静、黄静高级工程师指导完成，碳氧同位素测试由该实验室谷静高级工程师帮助完成；三位审稿人提出了宝贵的修改意见；在此一并致以诚挚感谢！

参考文献

Aharon P. 2005. Redox stratification and anoxia of the Early Precambrian oceans: Implications for carbon isotope excursions and oxidation events. Precambrian Research, 137(3-4): 202-222

Anbar AD, Duan Y, Lyons TW, Arnold GL, Kendall B, Creaser RA, Kaufman AJ, Gorgon GW, Scott C, Garvin J and Buick R. 2007. A whiff of oxygen before the Great Oxidation Event?. Science, 317(5846): 1903-1909

Baker AJ and Fallick AE. 1989a. Heavy carbon in two-billion-year-old marbles from Lofoten-Vesteralen, Norway: Implications for the Precambrian carbon cycle. Geochimica et Cosmochimica Acta, 53(5): 1111-1115

Baker AJ and Fallick AE. 1989b. Evidence from Lewisian limestones for isotopically heavy carbon in two-thousand-million-year-old sea water. Nature, 337(6205): 352-354

Banner JL and Hanson GN. 1990. Calculation of simultaneous isotopic and trace element variations during water-rock interaction with applications to carbonate diagenesis. Geochimica et Cosmochimica Acta, 54(11): 3123-3137

Banner JL. 1995. Application of the trace element and isotope geochemistry of strontium to studies of carbonate diagenesis. Sedimentology, 42(5): 805-824

Bekker A, Karhu JA, Eriksson KA and Kaufman AJ.. 2003. Chemostratigraphy of Paleoproterozoic carbonate successions of the Wyoming Craton: Tectonic forcing of biogeochemical change?. Precambrian Research, 120(3-4): 279-325

Bekker A, Holland HD, Wang PL, Rumble III D, Stein HJ, Hannah JL, Coetzee LL and Beukes NJ. 2004. Dating the rise of atmospheric oxygen. Nature, 427(6970): 117-120

Bekker A, Karhu JA and Kaufman AJ. 2006. Carbon isotope record for the onset of the Lomagundi carbon isotope excursion in the Great Lakes area, North America. Precambrian Research, 148(1-2): 145-180

Bekker A and Kaufman AJ. 2007. Oxidative forcing of global climate change: A biogeochemical record across the oldest Paleoproterozoic ice age in North America. Earth and Planetary Science Letters, 258(3-4): 486-499

Bekker A, Holmden C, Beukes NJ, Kenig F, Eglinton B and Patterson WP. 2008. Fractionation between inorganic and organic carbon during the Lomagundi (2.22~2.1Ga) carbon isotope excursion. Earth and Planetary Science Letters, 271(1-4): 278-291

Bottinga Y. 1969. Calculated fractionation factors for carbon and hydrogen isotope exchange in the system calcite-carbon dioxide-graphite-methane-hydrogen-water vapor. Geochimica et Cosmochimica Acta, 33(1): 49-64

Brand U and Veizer J. 1980. Chemical diagenesis of a multicomponent carbonate system-1: Trace elements. Journal of Sedimentary Petrology, 50(4): 1219-1236

Chen WY and Chen YJ. 2018. Records of the great oxidation event in the Hutuo Group, Shanxi, China: A reassessment of the δ¹³C_carb data. Acta Petrologica Sinica, 34(12): 3709-3720 (in Chinese with English abstract)

Chen WY, Chen YJ, Li QG, Li JR, Li KY, Shu SP, Chen X and Tong ZD. 2018. Detrital zircon U-Pb ages of the Sijizhuang glacial diamictites of the Hutuo Group in Wutai Shan, Shanxi Province and implication for the Great Oxidation Event. Earth Science Frontiers, 25(5): 1-18 (in Chinese with English abstract)

Chen YJ. 1987. The geologic environment catastrophe at about 2300Ma. In: Pan LQ (ed.). Young Geologists in Nanjing University, 119-125 (in Chinese)

Chen YJ. 1988. Catastrophe of the geologic environment at 2300Ma. In: Abstracts of International Symposium on Geochemistry and Mineralization of Proterozoic Mobile Belts. Tianjin: 11

Chen YJ, Hu SX and Fu SG. 1988. The REE geochemical character and its significance of the crystalline basement in southern margin of North China platform. In: Abstracts of International Symposium on Geochemistry and Mineralization of Proterozoic Mobile Belts. Tianjin: 13

Chen YJ. 1990. Evidences for the catastrophe in geologic environment at about 2300Ma and the discussions on several problems. Journal of Stratigraphy, 14(3): 178-184 (in Chinese with English abstract)

Chen YJ, Ji HZ, Fu SG and Zhou XP. 1991. The challenge to traditional geological theory by the discover of catastrophe events in 2.3 billion years: A New understanding of some major geological problems. Advanced in Earth Science, 6(2): 63-68 (in Chinese)

Chen YJ, Ouyang ZY, Yang QJ and Deng J. 1994. A new idea of the Archean-Proterozoic boundary. Geological Review, 40(5): 483-488 (in Chinese with English abstract)

Chen YJ. 1996. Progresses in application of sedimentary trace element geochemistry to probe crustal composition and environmental change. Geology Geochemistry, (3): 12-30 (in Chinese)

Chen YJ and Zhao YC. 1997. Geochemical characteristics and evolution of REE in the Early Precambrian sediments: Evidence from the southern margin of the North China Craton. Episodes, 20(2): 109-116

Chen YJ, Liu CQ, Chen HY, Zhang ZJ and Li C. 2000. Carbon isotope geochemistry of graphite deposits and ore-bearing khondalite series in North China: Implications for several geoscientific problems. Acta Petrologica Sinica, 16(2): 233-244 (in Chinese with English abstract)

Chen YJ and Tang HS. 2016. The great oxidation event and its records in North China Craton. In: Zhai MG, Zhao Y and Zhao TP (eds.). Main Tectonic Events and Metallogeny of the North China Craton. Singapore: Springer, 281-304

Chen YJ and Tang HS. 2018. Determination of the great oxidation event in North China. In: Zhai MG, Zhang LC and Chen B (eds.). Major Geological Events and Metallogenesis in the Precambrian of North China Craton. Beijing: Science Press, 155-179 (in Chinese)

Chen YJ, Chen WY, Li QG, Santosh M and Li JR. 2019. Discovery of the Huronian glaciation event in China: Evidence from glacigenic diamictites in the Hutuo Group in Wutai Shan. Precambrian Research, 320: 1-12

Cloud PE Jr. 1968. Atmospheric and hydrospheric evolution on the primitive earth. Science, 160(3829): 729-736

Condie K C. 2016. Earth as An Evolving Planetary System. London: Elsevier

Cowie JW, Ziegler W and Remane J. 1989. Stratigraphic commission accelerates progress, 1984 to 1989. Episodes, 12(2): 79-82

Dai AB. 1987. Coordination Chemistry. Beijing: Science Press (in Chinese)

Derry LA, Kaufman AJ and Jacobsen SB. 1992. Sedimentary cycling and environmental change in the Late Proterozoic: Evidence from stable and radiogenic isotopes. Geochimica et Cosmochimica Acta, 56(3): 1317-1329

Dong CY, Wang SJ, Liu DY, Wang JG, Xie HQ, Wang W, Song ZY and Wan YS. 2010. Late Palaeoproterozoic crustal evolution of the North China Craton and formation time of the Jingshan Group: Constraints from SHRIMP U-Pb zircon dating of meta-intermediate-basic intrusive rocks in eastern Shandong Province. Acta Petrologica Sinica, 27(6): 1699-1706 (in Chinese with English abstract)

Faure G. 1986. Principles of Isotope Geology. 2^nd Edition. New York: Wiley, 589

Feng WM, Zheng YF and Zhou JB. 2003. Carbon and oxygen isotope geochemistry of marbles from the Dabie-Sulu orogenic belt. Acta Petrologica Sinica, 19(3): 468-478 (in Chinese with English abstract)

Frakes LA. 1979. Evolution and genesis of global climates. In: Fitzgerald DL (ed.). Climates throughout Geologic Time. Amsterdam: Elsevier

Frauenstein F, Veizer J, Beukes N, Van Niekerk HS and Coetzee LL. 2009. Transvaal Supergroup carbonates: Implications for Paleoproterozoic δ¹⁸O and δ¹³C records. Precambrian Research, 175(1-4): 149-160

Guerrera Jr A, Peacock SM and Knauth LP. 1997. Large ¹⁸O and ¹³C depletions in greenschist facies carbonate rocks, Western Arizona. Geology, 25(10): 943-946

Guo QJ, Strauss H, Kaufman AJ, Schröder S, Gutzmer J, Wing B, Baker MA, Bekker A, Jin QS, Kim ST and Farquhar J. 2009. Reconstructing earth's surface oxidation across the Archean-Proterozoic transition. Geology, 37(5): 399-402

Hannah JL, Bekker A, Stein HJ, Markey RJ and Holland HD. 2004. Primitive Os and 2316Ma age for marine shale: Implications for Paleoproterozoic glacial events and the rise of atmospheric oxygen. Earth and Planetary Science Letters, 225(1-2): 43-52

Hoefs J. 1997. Stable Isotope Geochemistry. 4^th Edition. Berlin: Springer-Verlag, 201

Holland HD. 2002. Volcanic gases, black smokers, and the Great Oxidation Event. Geochimica et Cosmochimica Acta, 66(21): 3811-3826

Hudson JD. 1977. Stable isotopes and limestone lithification. Journal of the Geological Society, 133(6): 637-660

Huston DL and Logan GA. 2004. Barite, BIFs and bugs: Evidence for the evolution of the Earth's early hydrosphere. Earth and Planetary Science Letters, 220(1-2): 41-55

Jacobsen SB and Kaufman AJ. 1999. The Sr, C and O isotopic evolution of Neoproterozoic seawater. Chemical Geology, 161(1-3): 37-57

Ji HZ and Chen YJ. 1990. Khondalite series and its mineral resources. Geology and Exploration, (11): 11-13 (in Chinese)

Ji HZ, Chen YJ and Zhao YY. 1990. Graphite ore deposits and Khondalite series. China Nonmetallic Minerals Industry, (6): 9-11 (in Chinese)

Karhu JA. 1993. Palaeoproterozoic evolution of the carbon isotope ratios of sedimentary carbonates in the Fennoscandian Shield. Geological Survey Finland Bulletin, 371: 1-87

Karhu JA and Holland HD. 1996. Carbon isotopes and the rise of atmospheric oxygen. Geology, 24(10): 867-870

Kaufman AJ, Jacobsen SB and Knoll AH. 1993. The Vendian record of Sr and C isotopic variations in seawater: Implications for tectonics and paleoclimate. Earth and Planetary Science Letters, 120(3-4): 409-430

Kaufman AJ and Knoll AH. 1995. Neoproterozoic variations in the C-isotopic composition of seawater: Stratigraphic and biogeochemical implications. Precambrian Research, 73(1-4): 27-49

Kong QY. 2006. Deposits in Shandong Province. Jinan: Shandong Science and Technology Press, 537-551 (in Chinese)

Kump LR. 2008. The rise of atmospheric oxygen. Nature, 451(7176): 277-278

Kump LR, Junium C, Arthur MA, Brasier A, Fallick A, Melezhik V, Lepland A, Crne AE and Luo G. 2011. Isotopic evidence for massive oxidation of organic matter following the Great Oxidation Event. Science, 334(6063): 1694-1696

Lai Y, Chen C and Tang HS. 2012. Paleoproterozoic positive δ¹³C excursion in Henan, China. Geomicrobiology Journal, 29(3): 287-298

Lan XY. 1981. Study on genesis of the Nanshu graphite deposit and on feature of Precambrian graphite-bearing rock association, Shandong. Journal of Changchun University of Earth Sciences, (3): 30-42 (in Chinese)

Li HK, Li YF, Geng K, Zhuo CY, Zhang YB, Liang TT and Wang F. 2013. Palaeoproterozoic tectonic setting in the eastern Shandong Province. Geological Survey and Research, 36(2): 114-130 (in Chinese with English abstract)

Li KY, Chen YJ, She ZB, Tang HS and Chen WY. 2018. Carbon isotope compositions and geochemical characteristics of the Zhangshe graphite deposit of the Jingshan Group, Jiaobei. Earth Science Frontiers, 25(5): 19-33 (in Chinese with English abstract)

Liu PH, Liu FL, Wang F and Liu JH. 2011. U-Pb dating of zircons from Al-rich paragneisses of Jingshan Group in Shandong peninsula and its geological significance. Acta Petrologica et Mineralogica, 30(5): 829-843 (in Chinese with English abstract)

Liu PH, Liu FL, Wang F, Liu JH and Cai J. 2013. Petrological and geochronological preliminary study of the Xiliu ~2.1Ga meta-gabbro from the Jiaobei terrane, the southern segment of the Jiao-Liao-Ji Belt in the North China Craton. Acta Petrologica Sinica, 29(7): 2371-2390 (in Chinese with English abstract)

Lu LZ, Xu XC and Liu FL. 1996. The Early Precambrian Khondalite Series in North China. Changchun: Changchun Publishing House, 1-273 (in Chinese)

Luo GM, Ono S, Beukes NJ, Wang DT, Xie SC and Summons RE. 2016. Rapid oxygenation of Earth's atmosphere 2.33 billion years ago. Science Advances, 2(5): e1600134

Martin AP, Condon DJ, Prave AR and Lepland A. 2013. A review of temporal constraints for the Palaeoproterozoic large, positive carbonate carbon isotope excursion (the Lomagundi-Jatuli Event). Earth-Science Reviews, 127: 242-261

Melezhik VA, Fallick AE, Medvedev PV and Makarikhin VV. 1999. Extreme ¹³C_carb enrichment in ca.2.0Ga magnesite-stromatolite-dolomite-'red beds' association in a global context: A case for the world-wide signal enhanced by a local environment. Earth-Science Reviews, 48(1-2): 71-120

Melezhik VA, Gorokhov IM, Fallick AE and Gjelle S. 2001a. Strontium and carbon isotope geochemistry applied to dating of carbonate sedimentation: An example from high-grade rocks of the Norwegian Caledonides. Precambrian Research, 108(3-4): 267-292

Melezhik VA, Gorokhov IM, Kuznetsov AB and Fallick AE. 2001b. Chemostratigraphy of Neoproterozoic carbonates: Implications for 'blind dating'. Terra Nova, 13(1): 1-11

Melezhik VA, Fallick AE and Pokrovsky BG. 2005. Enigmatic nature of thick sedimentary carbonates depleted in ¹³C beyond the canonical mantle value: The challenges to our understanding of the terrestrial carbon cycle. Precambrian Research, 137(3-4): 131-165

Melezhik VA, Huhma H, Condon DJ, Fallick AE and Whitehouse MJ. 2007. Temporal constraints on the Paleoproterozoic Lomagundi-Jatuli carbon isotopic event. Geology, 35(7): 655-658

Melezhik VA, Bingen B, Fallick AE, Gorokhov IM, Kuznetsov AB, Sandstad JS, Solli A, Bjerkgård T, Henderson I, Boyd R, Jamal D and Monize A. 2008. Isotope chemostratigraphy of marbles in northeastern Mozambique: Apparent depositional ages and tectonostratigraphic implications. Precambrian Research, 162(3-4): 540-558

Melezhik VA, Prave AR, Hanski EJ, Fallick AE, Lepland A, Kump LR and Strauss H. 2013. Reading the Archive of Earth's Oxygenation: Volume 3: Global Events and the Fennoscandian Arctic Russia-Drilling Early Earth Project. Berlin, Heidelberg: Springer-Verlag

Murakami H, Wang B and Kitoh A. 2011. Future change of Western North Pacific Typhoons: Projections by a 20-km-mesh global atmospheric model. Journal of Climate, 24(4): 1154-1169

Pavlov AA and Kasting JF. 2002. Mass-independent fractionation of sulfur isotopes in Archean sediments: Strong evidence for an anoxic Archean atmosphere. Astrobiology, 2(1): 27-41

Qi L, Hu J and Gregoire DC. 2000. Determination of trace elements in granites by inductively coupled plasma mass spectrometry. Talanta, 51(3): 507-513

Ray JS. 2009. Carbon isotopic variations in fluid-deposited graphite: Evidence for multicomponent Rayleigh isotopic fractionation. International Geology Review, 51(1): 45-57

Sanyal P, Acharya BC, Bhattacharya SK, Sarkar A, Agrawal S and Bera MK. 2009. Origin of graphite, and temperature of metamorphism in Precambrian Eastern Ghats Mobile Belt, Orissa, India: A carbon isotope approach. Journal of Asian Earth Sciences, 36(2-3): 252-260

Schidlowski M, Eichmann R and Junge CE. 1975. Precambrian sedimentary carbonates: Carbon and oxygen isotope geochemistry and implications for the terrestrial oxygen budget. Precambrian Research, 2(1): 1-69

Schidlowski M, Eichmann R and Junge CE. 1976. Carbon isotope geochemistry of the Precambrian Lomagundi carbonate province, Rhodesia. Geochimica et Cosmochimica Acta, 40(4): 449-455

Schidlowski M. 1988. A 3, 800-million-year isotopic record of life from carbon in sedimentary rocks. Nature, 333(6171): 313-318

Tam PY, Zhao GC, Liu FL, Zhou XW, Sun M and Li SZ. 2011. Timing of metamorphism in the Paleoproterozoic Jiao-Liao-Ji Belt: New SHRIMP U-Pb zircon dating of granulites, gneisses and marbles of the Jiaobei massif in the North China Craton. Gondwana Research, 19(1): 150-162

Tang GJ, Chen YJ, Huang BL and Chen CX. 2004. Paleoproterozoic δ¹³C_carb positive excursion event: Research progress on 2.3Ga catastrophe. Journal of Mineralogy and Petrology, 24(3): 103-109 (in Chinese with English abstract)

Tang HS, Chen YJ, Wu G and Lai Y. 2008. The C-O isotope composition of the Liaohe Group, northern Liaoning Province and its geologic implications. Acta Petrologica Sinica, 24(1): 129-138 (in Chinese with English abstract)

Tang HS, Wu G and Lai Y. 2009. The C-O isotope geochemistry and genesis of the Dashiqiao magnesite deposit, Liaoning Province, NE China. Acta Petrologica Sinica, 25(2): 455-467 (in Chinese with English abstract)

Tang HS, Chen YJ, Wu G and Lai Y. 2011. Paleoproterozoic positive δ¹³C_carb excursion in the northeastern Sino-Korean craton: Evidence of the Lomagundi Event. Gondwana Research, 19(2): 471-481

Tang HS and Chen YJ. 2013. Global glaciations and atmospheric change at ca.2.3Ga. Geoscience Frontiers, 4(5): 583-596

Tang HS, Chen YJ, Santosh M, Zhong H, Wu G and Lai Y. 2013a. C-O isotope geochemistry of the Dashiqiao magnesite belt, North China Craton: Implications for the Great Oxidation Event and ore genesis. Geological Journal, 48(5): 467-483

Tang HS, Chen YJ, Santosh M, Zhong H and Yang T. 2013b. REE geochemistry of carbonates from the Guanmenshan Formation, Liaohe Group, NE Sino-Korean Craton: Implications for seawater compositional change during the Great Oxidation Event. Precambrian Research, 227: 316-336

Tang HS, Chen YJ, Li KY, Chen WY, Zhu XQ, Ling KY and Sun XH. 2016. Early Paleoproterozoic metallogenic explosion in North China Craton. In: Zhai MG, Zhao Y and Zhao TP (eds.). Main Tectonic Events and Metallogeny of the North China Craton. Singapore: Springer, 305-328

Tang HS and Chen YJ. 2018a. Reconstruction of the Great oxidation Event sequence. In: Zhai MG, Zhang LC and Chen B (eds.). Major Geological Events and Metallogenesis in the Precambrian of North China Craton. Beijing: Science Press, 145-155 (in Chinese)

Tang HS and Chen YJ. 2018b. Enrichment of large amount of elements in supergene environment. In: Zhai MG, Zhang LC and Chen B (eds.). Major Geological Events and Metallogenesis in the Precambrian of North China Craton. Beijing: Science Press, 180-188 (in Chinese)

Tang J, Zheng YF, Wu YB, Zha XP and Zhou JB. 2004. Zircon U-Pb ages and oxygen isotopes of metamorphic rocks in the western part of the Shandong Peninsula. Acta Petrologica Sinica, 20(5): 1063-1086 (in Chinese with English abstract)

Valley JM. 1986. Stable isotope geochemistry of metamorphic rocks. Reviews in Mineralogy, 16(1): 445-489

Vallini DA, Cannon WF and Schulz KJ. 2006. Age constraints for Paleoproterozoic glaciation in the Lake Superior region: Detrital zircon and hydrothermal xenotime ages for the Chocolay Group, Marquette Range Supergroup. Canadian Journal of Earth Sciences, 43(5): 571-591

Veizer J and Hoefs J. 1976. The nature of O¹⁸/O¹⁶ and C¹³/C¹² secular trends in sedimentary carbonate rocks. Geochimica et Cosmochimica Acta, 40(11): 1387-1395

Veizer J. 1983. Chemical diagenesis of carbonates: Theory and application of trace element technique. In: Arthur MA, Anderson TF, Kaplan IR, Veizer J and Land LS (eds.). Stable Isotopes in Sedimentary Geology. V.10. Society of Economic Paleontologists and Mineralogists, 3.1-3.100

Veizer J, Clayton RN and Hinton RW. 1992a. Geochemistry of Precambrian carbonates: IV.Early Paleoproterozoic (2.25±0.25Ga) seawater. Geochimica et Cosmochimica Acta, 56(3): 875-885

Veizer J, Plumb KA, Clayton RN, Hinton RW and Grotzinger JP. 1992b. Geochemistry of Precambrian carbonates: V.Late Paleoproterozoic seawater. Geochimica et Cosmochimica Acta, 56(6): 2487-2501

Veizer J, Ala D, Azmy K, Bruckschen P, Buhl D, Bruhn F, Carden GAF, Diener A, Ebneth S, Godderis Y, Jasper T, Korte C, Pawellek F, Podlaha OG and Strauss H. 1999. ⁸⁷Sr/⁸⁶Sr, δ¹³C and δ¹⁸O evolution of Phanerozoic seawater. Chemical Geology, 161(1-3): 59-88

Wada H and Suzuki K. 1983. Carbon isotopic thermometry calibrated by dolomite-calcite solvus temperatures. Geochimica et Cosmochimica Acta, 47(4): 697-706

Wan YS, Song B, Liu DY, Wilde SA, Wu JS, Shi YR, Yin XY and Zhou HY. 2006. SHRIMP U-Pb zircon geochronology of Palaeoproterozoic metasedimentary rocks in the North China Craton: Evidence for a major Late Palaeoproterozoic tectonothermal event. Precambrian Research, 149(3-4): 249-271

Xie SW, Wang SJ, Xie HQ, Liu SJ, Dong CY, Ma MZ, Liu DY and Wan YS. 2014. SHRIMP U-Pb dating of detrital zircons from the Fenzishan Group in eastern Shandong, North China Craton. Acta Petrologica Sinica, 30(10): 2989-2998 (in Chinese with English abstract)

Yang JH. 2000. Age and metallogenic dynamics of gold mineralization in Jiaodong Peninsula, eastern China: Constraints in the interaction of mantle/crust and metallogenesis/lithospheric evolution. Ph. D. Dissertation. Beijing: Institute of Geology and Geophysics, Chinese Academy of Science, 1-133 (in Chinese with English summary)

Young GM. 2013. Precambrian supercontinents, glaciations, atmospheric oxygenation, metazoan evolution and an impact that may have changed the second half of Earth history. Geoscience Frontiers, 4(3): 247-261

Young GM. 2019. Aspects of the Archean-Proterozoic transition: How the great Huronian Glacial Event was initiated by rift-related uplift and terminated at the rift-drift transition during break-up of Lauroscandia. Earth-Science Reviews, 190: 171-189

Zhai MG, Zhang LC and Chen B. 2018. Precambrian Great Geologic and Metallogenic Events in North China Craton. Beijing: Science Press, 1-29 (in Chinese)

Zhang ZQ and Liu MW. 1996. Rocks and Strata of Shandong Province. Wuhan: China University of Geosciences Press, 1-310 (in Chinese)

Zhou XW, Dong YS and Wei CD. 2001. The genesis and evolution of the metamorphic minerals of khondalite series in Nanshu district of Shandong Province. Journal of Changchun University of Science and Technology, 31(2): 116-121 (in Chinese with English abstract)

陈威宇, 陈衍景. 2018. 大氧化事件在山西滹沱群中的记录:碳酸盐岩碳同位素资料分析. 岩石学报, 34(12): 3709-3720.

陈威宇, 陈衍景, 李秋根, 李建荣, 李凯月, 疏孙平, 陈西, 佟子达. 2018. 山西五台山滹沱群四集庄冰碛岩碎屑锆石年龄及其对大氧化事件研究意义. 地学前缘, 25(5): 1-18.

陈衍景. 1987.论23亿年前地质环境的突变.见: 潘龙泉编.南大青年地质学家, 119-125

陈衍景. 1990. 23亿年地质环境突变的证据及若干问题的讨论. 地层学杂志, 14(3): 178-184.

陈衍景, 季海章, 富士谷, 周小平. 1991. 23亿年灾变事件的揭示对传统地质理论的挑战——关于某些重大地质问题的新认识. 地球科学进展, 6(2): 63-68.

陈衍景, 欧阳自远, 杨秋剑, 邓健. 1994. 关于太古宙-元古宙界线的新认识. 地质论评, 40(5): 483-488.

陈衍景, 杨秋剑, 邓建, 季海章, 富士谷, 周小平, 林清. 1996. 地球演化的重要转折——2300Ma时地质环境灾变的解释及其意义. 地质地球化学, (3): 106-125.

陈衍景, 刘丛强, 陈华勇, 张增杰, 李超. 2000. 中国北方石墨矿床及赋矿孔达岩系碳同位素特征及有关问题讨论. 岩石学报, 16(2): 233-244.

陈衍景, 汤好书. 2018.华北大氧化事件的确定.见: 翟明国, 张连昌, 陈斌编.华北克拉通前寒武纪重大地质事件与成矿.北京: 科学出版社, 155-179

戴安邦. 1987. 配位化学. 北京: 科学出版社.

董春艳, 王世进, 刘敦一, 王金光, 颉颃强, 王伟, 宋志勇, 万渝生. 2010. 华北克拉通古元古代晚期地壳演化和荆山群形成时代制约——胶东地区变质中-基性侵入岩锆石SHRIMP U-Pb定年. 岩石学报, 27(6): 1699-1706.

冯伟民, 郑永飞, 周建波. 2003. 大别-苏鲁造山带大理岩碳氧同位素地球化学研究. 岩石学报, 19(3): 468-478.

季海章, 陈衍景. 1990. 孔达岩系及其矿产. 地质与勘探, (11): 11-13.

季海章, 陈衍景, 赵懿英. 1990. 孔达岩系与石墨矿床. 建材地质, (6): 9-11.

孔庆友. 2006. 山东矿床. 济南: 山东科学技术出版社, 537-551.

兰心俨. 1981. 山东南墅前寒武纪含石墨建造的特征及石墨矿床的成因研究. 长春地质学院学报, (3): 30-42.

李洪奎, 李逸凡, 耿科, 禚传源, 张玉波, 梁太涛, 王峰. 2013. 鲁东地区古元古界形成的大地构造环境探讨. 地质调查与研究, 36(2): 114-130.

李凯月, 陈衍景, 佘振兵, 汤好书, 陈威宇. 2018. 胶北荆山群张舍石墨矿碳同位素特征及其地质意义. 地学前缘, 25(5): 19-33.

刘平华, 刘福来, 王舫, 刘建辉. 2011. 山东半岛荆山群富铝片麻岩锆石U-Pb定年及其地质意义. 岩石矿物学杂志, 30(5): 829-843.

刘平华, 刘福来, 王舫, 刘建辉, 蔡佳. 2013. 胶北西留古元古代~2.1Ga变辉长岩岩石学与年代学初步研究. 岩石学报, 29(7): 2371-2390.

卢良兆, 徐 .学纯, 刘福来. 1996. 中国北方早前寒武纪孔兹岩系. 长春: 长春出版社, 1-273.

唐国军, 陈衍景, 黄宝玲, 陈从喜. 2004. 古元古代δ¹³C_carb正向漂移事件: 2.3Ga环境突变研究的进展. 矿物岩石, 24(3): 103-109.

汤好书, 陈衍景, 武广, 赖勇. 2008. 辽北辽河群碳酸盐岩碳-氧同位素特征及其地质意义. 岩石学报, 24(1): 129-138.

汤好书, 武广, 赖勇. 2009. 辽宁大石桥菱镁矿床的碳氧同位素组成和成因. 岩石学报, 25(2): 455-467.

汤好书, 陈衍景. 2018a.大氧化事件序列重建.见: 翟明国, 张连昌, 陈斌编.华北克拉通前寒武纪重大地质事件与成矿.北京: 科学出版社, 145-155

汤好书, 陈衍景. 2018b.表生环境巨量元素富集.见: 翟明国, 张连昌, 陈斌编.华北克拉通前寒武纪重大地质事件与成矿.北京: 科学出版社, 180-188

唐俊, 郑永飞, 吴元保, 查向平, 周建波. 2004. 胶东地块西部变质岩锆石U-Pb定年和氧同位素研究. 岩石学报, 20(5): 1063-1086.

谢士稳, 王世进, 颉颃强, 刘守偈, 董春艳, 马铭株, 刘敦一, 万渝生. 2014. 华北克拉通胶东地区粉子山群碎屑锆石SHRIMP U-Pb定年. 岩石学报, 30(10): 2989-2998.

杨进辉. 2000.胶东地区金矿床成矿时代及其成矿地球动力学背景——兼论壳幔相互作用与成岩成矿.博士学位论文.北京: 中国科学院地质与地球物理研究所, 1-133

翟明国, 张连昌, 陈斌. 2018. 华北克拉通前寒武纪重大地质事件与成矿. 北京: 科学出版社, 1-29.

张增奇, 刘明渭. 1996. 山东省岩石地层. 武汉: 中国地质大学出版社, 1-310.

周喜文, 董永胜, 魏存弟. 2001. 山东南墅地区孔兹岩系变质矿物的成因及演化. 长春科技大学学报, 31(2): 116-121.


岩石学报 2020, Vol. 36 Issue (4): 1059-1075, doi: 10.18654/1000-0569/2020.04.06	PDF