Hydrogen-Oxygen-Sulfur Isotope Composition of Shihoushan Pyrite and Tungsten Polymetallic Deposit, Southern Jiangxi
-
摘要: 狮吼山硫铁-钨多金属矿床位于银坑-青塘整装勘查区北部,是赣南地区唯一大型硫铁矿床。磁黄铁矿-黄铁矿(-黄铜矿-白钨矿)矿体赋存于石炭系梓山组上段地层中含铁、含钙层位,主要形成于石英-硫化物阶段。本文通过分析原生矿石矿物中H-O-S同位素组成特征,结合Pb同位素和成矿年代测试结果,探讨成矿流体来源及成矿演化过程。矿石硫化物中δ34S组成特征(-5.50‰~-0.20‰,集中于-3.0‰~0.0‰)显示,硫源以岩浆硫为主,较宽的变化范围预示成矿流体遭受了叠加和改造作用。δD-δ18O同位素组成主要集中于岩浆水与变质水重叠区域(δD=-74.4‰~-48.0‰,δ18OH2O=3.76‰~10.86‰),说明成矿流体以岩浆水和变质水为主,后期有少量的天水混入。综合分析认为,该矿床成矿流体主要来自深部岩浆水,岩浆热液与含钙地层的接触交代作用形成大规模变质流体,再加上少量的天水混入,流体间的不混溶作用使成矿物质在岩体与含钙层位接触部位富集沉淀,形成热液充填交代型矿床。
-
关键词:
- 硫铁-钨多金属矿床 /
- H-O-S同位素 /
- 成矿流体演化 /
- 岩浆期后热液充填交代型矿床 /
- 赣南
Abstract:BACKGROUNDShihoushan Pyrite and Tungsten polymetallic deposit, located in the northern of Yinkeng-Qingtang Au-Ag polymetallic integrated exploration area, is the only large pyrite deposit in southern Jiangxi. The pyrrhotite-pyrite (-chalcopyrite-scheelite) orebody hosts in calciferous sandstone in the Zishan Formation of Carboniferous System, mainly formed in the quartz-sulfide stage. OBJECTIVESIn order to better understand the ore-forming material source and the evolution processes of Shihoushan deposit, the primary ores were selected as laboratory raw materials for stable isotope testing. METHODSH-O-S isotope composition of primary ore minerals was analyzed combined with the Pb isotope and metallogenic age results. The ore-forming fluid source and ore-forming evolution process are discussed. RESULTSThe δ34S values range from -5.50‰ to -0.20‰, which are mainly concentrated at -3.0‰-0.0‰ (n=11), show the typical signature of mantle S. The wide range of variation indicates that the ore-forming fluid has been subjected to superposition and modification. H-O isotope analyses show that δD=-74.4‰--48.0‰ (n=9), δ18OH2O=3.76‰-10.86‰ (n=9), indicating that the ore-forming fluid is composed mainly of magmatic water and metamorphic water, with minor meteoric water. CONCLUSIONSAccording to the comprehensive analysis, the ore-forming fluid of this deposit mainly comes from deep magmatic water. The contact between magmatic hydrothermal fluid and calcium-bearing strata forms a large-scale metamorphic fluid, mixed with a small amount of meteoric water. The fluid immiscibility makes the ore-forming materials precipitate in the contact between the rock mass and the calcium-bearing strata, forming a hydrothermal filling and metasomatic deposit. -
碳循环不仅是自然生态系统中物质和能量迁移转化的重要过程之一,而且对调节全球气候环境变化具有重要作用[1]。陆地碳酸盐化学风化作用(即岩溶作用)能够吸收固定大气CO2,形成溶解无机碳(DIC, dissolved inorganic carbon),且DIC随后被输送到河流、湖泊和海洋中,并被这些水体中的水生光合生物吸收利用,最终转化成有机碳,埋藏于水体底部,因此碳酸盐风化作用被认为是一种重要的陆地碳汇[2]。国内外研究表明,全球碳酸盐风化作用吸收固定的CO2通量巨大,每年可达0.11~0.80Pg碳,约占陆地碳汇总量的7%~25%[2-5]。刘再华等[5-6]认为,陆地碳酸盐风化-水生光合生物碳泵耦合过程至少在短时间尺度上能够代表一种净碳汇。由于陆地上碳酸盐岩分布较广,且碳酸盐岩溶蚀作用较硅酸盐岩强烈[7],其固碳能力巨大,在全球碳收支平衡中起着重要作用。另外,岩溶作用具有动态变化的特点[8-9],能够敏感地响应气候变化和地表植被覆盖/土地利用变化[10-11],故通过在岩溶地区增加植被或改变土地利用方式,有潜力为未来人工增汇、实现碳中和作出巨大贡献。
然而,目前岩溶作用碳汇效应还受到地质学家们的质疑,他们认为碳酸盐风化过程中吸收固定的CO2会通过水中CO2脱气和碳酸钙沉淀作用而重新返回大气[12]。这涉及岩溶碳汇的稳定性这一关键问题。在岩溶地区,岩溶泉或地下河水出露地表后,在水中和空气高二氧化碳分压(pCO2)梯度驱动下,常会发生CO2脱气作用。特别地,在高地势变化区、跌水和瀑布等地,CO2脱气作用尤为明显,并伴随有碳酸钙沉积[13-17]。碳酸钙沉积作用会将水中重碳酸(HCO3-)转变为CO2,释放进入大气,这会在一定程度上抵消了岩溶作用吸收固定的CO2量,从而给岩溶碳汇估算带来不确定性。因此,有必要对岩溶地下水出露地表后DIC变化及CO2脱气作用开展深入研究,以便更好地了解岩溶作用碳汇效应。前人研究表明,岩溶泉补给溪流系统CO2脱气具有时空变化,且受到了诸多因素影响,如源泉水水化学、流量、地形和生物活动等[13-17]。Hoffer-French等[13]研究了美国弗吉尼亚州Falling泉地表溪流,发现CO2脱气作用控制了溪水水化学演化,使得水中pH值和方解石饱和指数(SIc)沿流程逐渐升高,在跌水和瀑布处还产生了碳酸钙沉积。周小萍等[15-16]研究了重庆南川柏树湾表层岩溶泉溪流,也发现溪水DIC浓度沿流程逐渐降低,其降低的快慢反映了CO2脱气的强弱,他们发现脱气速率具有“枯水期>平水期>丰水期”的特点。然而,在岩溶泉补给的水池系统中,CO2脱气作用较弱,且水池中的水生光合生物能够吸收利用DIC,产生类似于海洋的“生物碳泵”效应[18]。其他研究表明,在地形平缓、流量很大的岩溶地下河补给的河流[19-20]或流经岩溶区的大型河流中[21],水生生物代谢作用也能吸收固定DIC,从而抑制CO2脱气。岩溶地下水的输入甚至能够抑制河流中生物呼吸作用产生的CO2排放[22]。尽管前人研究已经对岩溶地下水补给地表溪/河流有了一些认识,但是大多数研究只侧重于溪水地球化学时空变化特征,缺少对不同类型CO2脱气的定量性研究,对不同地形下溪流CO2脱气作用认识也还有限。
溪流水化学指标和溶解无机碳碳同位素组成(δ13CDIC)变化可以表征水体碳元素迁移转化,为研究溪流CO2脱气作用提供了重要手段[14, 20]。本文以广西桂林长流水表层岩溶泉补给的一条长约2.7km的地表溪流为例,利用水化学和同位素质谱测试技术,从泉口向下游,沿流程对溪水进行监测研究。根据溪水水化学指标和δ13CDIC变化,分析探讨溪流CO2脱气过程、通量及其影响因素,并估算CO2脱气通量及其相较于岩溶作用过程中吸收固定的大气CO2量的比例。根据长流水这一高地势、低流量且有碳酸钙沉积的环境下的溪流CO2脱气作用及影响因素,探讨和评估岩溶碳汇的稳定性,以期为岩溶作用碳汇效应提供更多证据。
1. 研究区概况
长流水岩溶泉位于广西桂林市雁山区大埠乡付中村,海拔约240m(图 1)。该区地处漓江西岸,属于中国南方喀斯特世界自然遗产核心区,地貌为典型的峰丛洼地。泉水补给区为一石峰环绕的岩溶洼地,面积约1.5km2。区内出露地层主要为上泥盆统亮晶颗粒灰岩,局部为白云岩化泥晶颗粒灰岩或夹细晶白云岩。泉水自洼地西侧的石峰西坡中部流出,为典型的表层岩溶泉,其流量季节变化较大,雨季较高,旱季较低[14]。泉水出露后,沿一较陡的渠道向坡下流动,而后进入峰林平原区,最终汇入良丰河的一条支流中(图 1),该段溪流总长度约2.7km。
2. 实验部分
2.1 样品采集
于2020年5月28日,自泉口向下游沿流程设置26个点(C1~C26;图 1)采集水样,样品编号为S1~S26。其中,在C1~C14段,相邻采样点间隔为15~30m。在C20和C24点,对支流水及混合前(上游5m处)和混合后(下游10m处)干流水进行了采样。在溪流横断面中部进行现场测试和采样。用预先酸洗过的聚乙烯瓶采集30mL水样,并加入2~3滴饱和HgCl2溶液,用于δ13CDIC测试。水样带回实验室后,置于冰箱中4℃保存,48h内进行测试。
2.2 实验测试
在野外,使用便携式多参数水质分析仪(德国WTW)现场测定溪水水温、pH和电导率(EC),其测量精度分别为0.01℃、0.01pH单位和1μS/cm。利用碱度和钙度测试盒(德国Merck公司)现场测定水中碱度和Ca2+浓度,精度分别为0.1mmol/L和2mg/L。在每个点进行两次滴定,取平均值。根据所测水温、pH、碱度和Ca2+浓度,使用PHREEQC 3.0程序计算水中pCO2、总无机碳(TIC)和HCO3-浓度以及SIc值。
水样δ13CDIC值采用GasBenchⅡ-IRMS在线磷酸分解法[23]进行测定。实验原理为:在25℃实验条件下,水样中的溶解无机碳与磷酸反应全部生成CO2,在载气(He)携带下经恒温色谱柱分离纯化,导入同位素比值质谱仪进行碳同位素比值测定。测试仪器为MAT253型同位素质谱仪(美国ThermoFinnigan公司),前处理装置为GasBenchⅡ(美国ThermoFisher公司)。仪器主要工作条件为:磁场加速电压9533kV;样品平衡温度25℃;色谱柱温度45℃;氦气流速2mL/min;样品采集时间630s;样品采集参考气峰5次,样品峰8次;系统稳定性≥5×10-6min;相对灵敏度≥0.2A/mbar。δ13C值测试结果相对于VPDB标准,精度优于±0.2‰。
2.3 坡度、流量和流速测量
在野外,利用数字水准仪(Trimble DiNi)测量了C1~C14渠段的高度,观测精度达0.3mm,安平精度±0.2″,并根据距离测量结果,换算成坡度。采样当天,在泉口处的矩形堰测定水位,换算为流量。在其他采样点,利用LSH10-1QC型便携式流速仪测定流速,测量范围为0.01~5.00m/s,准确度为1.0%±1cm/s,并根据渠宽和水深计算溪水流量。
3. 结果与讨论
3.1 沿流程溪水水化学变化
采样当天,长流水岩溶泉流量为27.9L/s,由于在采样前两天暴雨,该流量高于雨季平均值。泉水呈弱碱性,pH值为7.3。先前的研究已表明该泉水离子组成主要为Ca2+和HCO3-,而其他离子浓度都很低[15]。本研究中所测泉水Ca2+和HCO3-浓度分别为2.5mmol/L和4.9mmol/L,这个浓度与其他灰岩地区的岩溶泉水的浓度相近[9, 24-25],反映了流域和含水层内的灰岩溶蚀控制。泉水pCO2值为10-1.8 atm,这个值与桂林地区土壤pCO2值(10-1.6~10-1.4atm)[26]接近。水中SIc值约为0.2,显示系统内碳酸盐岩溶蚀作用已经达到化学平衡。
溪水水化学数据显示,自泉口(C1点)向坡下至C14点,各水化学指标呈现出明显的变化趋势(表 1;图 2)。总体上,水温、pH和SIc沿流程逐渐升高,分别升高了2.2℃、0.9个pH单位和0.9;而pCO2、EC、Ca2+和DIC浓度不断降低,分别下降了85%、34μS/cm、0.2mmol/L和0.7mmol/L。水温升高主要是由于外界气温较高,溪水逐渐与外界进行热量交换。Wang等[14]发现在冬天气温较低时,泉水出露后沿流程水温则会逐渐降低。溪水其他指标的变化则主要与溪水CO2脱气和碳酸钙沉积作用有关。具体解释如下。
表 1. 各渠段地球化学指标沿流程变化幅度Table 1. Amplitudes of downstream variation in geochemical indicators in different segments of the stream渠段 坡度
(°)高差
(m)距离
(m)ΔT
(℃)ΔpH ΔEC
(μS/cm)Δ[Ca2+]
(mmol/L)Δ[HCO3-]
(mmol/L)ΔSIc ΔlgpCO2 Δδ13CDIC
(‰)C1~C14 9.9 47 270 +2.2 +0.9 -34 -0.2 -0.7 +0.9 -0.9 +1.8 C14~C18 0.8 4 295 +1.1 -0.2 -9 -0.1 -0.3 -0.1 +0.2 -0.5 C18~C19 3.8 13 190 +0.6 +0.2 -24 -0.1 -0.2 +0.1 -0.2 +0.6 C21~C23 0.6 5 493 +0.6 -0.4 +18 +0.1 -0.1 -0.3 +0.4 -0.5 C25~C26 0.2 5 1450 +1.8 +0.3 +4 -0.1 -0.1 +0.3 -0.2 +1.0 注:Δ值为“+”,代表沿流程升高;Δ值为“-”,代表沿流程降低。 由于泉水pCO2值约是大气pCO2值(10-3.5atm)的40倍,故泉水出露地表后,在较高的pCO2梯度驱动下,水中的溶解态CO2很容易通过分子扩散形式逸出进入空气。可以看到,沿流程溪水中pCO2呈指数降低,至C13点下降到10-2.8atm以下(图 2)。最先逸出的CO2主要为泉水中来自土壤的物理溶解的CO2[14]。从以下岩溶作用反应式(1)和(2)可以看出,水中CO2不断逸出会逐渐消耗H+离子,从而使得pH值沿流程不断增高,在C13点达到8.2。此外,由于pCO2降低会使得水中化学平衡遭到改变,从而导致水中矿物饱和指数不断升高[27]。
CO2+H2O↔H2CO3 (1) H2CO3↔H++HCO−3 (2) CaCO3+CO2+H2O↔Ca2++2HCO−3 (3) 至离泉口约80m的C6点,溪水SIc值增大至0.8,在C12点达到1.1。已有研究表明,当SIc值达到或超过碳酸钙晶体成核的临界阈值0.5~0.8时,水流中就会发生碳酸钙沉淀[13]。这与野外观察一致,在C5~C9的渠道里可以看见大量钙华沉积。碳酸钙沉淀作用会消耗掉水中Ca2+和HCO3-[见反应式(3)的逆反应],从而使得溪水Ca2+和HCO3-浓度以及EC值沿流程逐渐降低。与此同时,沉淀过程会伴随CO2气体的生成和脱气,每沉积1mol的CaCO3,就会有1mol的HCO3-通过脱水/脱羟基作用释放出1mol的CO2气体[见反应式(3)],生成气态CO2很容易从水中逸出。
C1~C14渠段水化学出现较大变化的原因主要与渠道地形及其控制的水动力条件有关。自泉口(C1点)至C14点,渠道较陡,坡度为3°~23°(平均值为10°)。该渠道中水层较薄,多为2~5cm,且渠道内发育有多处边石坝、小跌水以及水坑,水流呈紊流。在这种水动力条件下,溪水中CO2脱气速率高,SIc值升高速度快,有助于在渠道中发生碳酸钙沉积[13-14]。
自C14至C18点,渠道坡度大大变缓,仅为0.8°,溪水水化学指标变化幅度较小(表 1;图 2)。在C17点,受部分稻田水的混入影响,溪水水化学有所变化。在C18~C19段,渠道转而变陡,坡度约为4°,溪水pH和SIc值又有所升高,而EC值和Ca2+浓度下降,指示发生了碳酸钙沉积。在C20点,有另一岩溶泉补给的溪流A汇入渠道中,溪流A流量为23.0L/s,其流程较短,所以溪水具有较低的温度、pH、SIc值和较高的pCO2、EC、Ca2+浓度(表 2)。该支流的混入也使得干流溪水水化学组成遭到了稀释。值得注意的是,混合后的溪水pH和SIc值相比混合前的两股水都要低。这主要是因为碳酸盐体系中的化学平衡呈指数变化,而溶质混合为线性变化[28],所以两股流量相近的水流混合会使得水流矿物饱和度变得更低。
表 2. 支流汇入后溪水地球化学变化Table 2. Changes in geochemistry of stream water caused by the mixing of tributaries样品编号 状态 HCO3-浓度
(mg/L)Ca2+浓度
(mg/L)EC
(μS/cm)pH pCO2
(10-3atm)SIc δ13CDIC
(‰)S19 汇合前 228.6 84 378 8.3 1.3 1.1 -13.0 S20 支流A 257.2 90 417 8.0 2.2 0.9 -14.7 S21 汇合后 248.8 86 397 7.9 2.7 0.8 -13.8 S23 汇合前 248.2 88 401 7.7 4.7 0.6 -14.0 S24 支流B 274.7 92 425 7.6 7.2 0.5 -15.7 S25 汇合后 255.3 88 409 7.6 5.5 0.5 -14.5 自C21至C22点,渠道坡度小,流速慢,溪水水化学变化幅度不大(图 2)。而在C22~C23段,溪水化学有较大改变,其原因也是有其他水流从两边水里流入溪流中。同C20处,C24处支流B的混入也使得溪水水化学有所改变,但是由于支流B流量(21.3L/s)相较于干流较小,并没有使得混合后的溪水pH和SIc值相比混合的水流都低。C25至C26点为平原区,在长约1.5km的渠段内,坡度仅为0.3°,溪水水化学变化幅度不大。
除了受物理化学过程控制外,溪水水化学的变化可能还受到水生生物代谢作用(如光合作用)的影响[18-19]。在C1~C14的渠道中,可以看见大量苔藓和藻类的生长。研究发现,藻类可以利用水中CO2和HCO3-进行光合作用[29],从而会改变水体水化学。但是,由于C1~C14渠道较陡,水流速度快(1~4m/s),且水中DIC浓度较大(>4mmol/L),生物代谢作用对水化学的影响应该有限。只有在流速缓慢、较深的渠段,水生生物代谢作用才有可能会对溪水水化学产生影响。
3.2 沿流程溪水δ13CDIC值的变化特征
长流水泉水的δ13CDIC值为-15.3‰,这个值与开放系统下碳酸盐溶蚀产生的δ13CDIC理论值(-16‰)接近[30]。支流B的δ13CDIC值也较低,为-15.72‰,说明该支流尚未发生强烈的CO2脱气作用。溪流水的δ13CDIC值变化范围为-15.2‰~-12.9‰。从泉口向下游,溪水δ13CDIC值变化趋势与pH、SIc一致,而与pCO2、EC、Ca2+浓度相反。自C1至C14点,溪水δ13CDIC值沿流程不断增高,增幅达1.8‰(图 2)。同样地,在C18~C19、C25~C26段,溪水δ13CDIC值也呈现不同程度的升高趋势(表 1)。具有不同δ13CDIC值的支流或其他水流的混入会显著地改变溪水δ13CDIC值。例如,支流A和B的汇入都使得溪水δ13CDIC有所偏负(表 2)。
溪水δ13CDIC演化与水温、CO2脱气、碳酸钙沉积和生物作用等有关[20, 31]。一般地,气态CO2和溶解态CO2之间的同位素分馏较小,仅为1.1‰[30],故水中CO2分子扩散对δ13CDIC的影响较小。然而,由于HCO3-脱水/脱羟基过程中的同位素分馏要大得多,在5~25℃时,δ13CDIC值为7.3‰~10.1‰,甚至更高,可达20‰[32]。因此,碳酸钙沉积引发的CO2脱气能够使得溪水δ13CDIC遭到很大改变。在长流水,δ13CDIC值和HCO3-浓度呈线性负相关(图 3),表明碳酸钙沉淀引起的CO2脱气可能是导致溪水中剩余DIC发生同位素富集的主要原因。由图 3也可知,受碳酸钙沉积影响的上游渠道,其δ13CDIC值和HCO3-浓度相比下游具有更高的斜率。生物代谢作用,如水生生物光合作用吸收DIC,也可能会改变δ13CDIC值[20, 33]。特别地,在流速缓慢、水较深之地,如C25~C26的某些渠段,生物代谢作用可能对溪水δ13CDIC升高有一定影响。
3.3 溪水CO2脱气通量及影响因素
3.3.1 溪水CO2脱气通量计算
根据溪水水化学和δ13CDIC沿流程变化可知,从泉口向下游存在两种不同类型的CO2脱气:一是,泉水中来自土壤的物理溶解的CO2脱气;二是,渠道里碳酸钙沉淀引发的CO2脱气。本文对这两种类型的CO2脱气进行定量分析,并估算岩溶泉水出露地表后沿溪流发生的CO2释放量占碳酸盐岩溶蚀过程中吸收固定的CO2量的比例。
一般地,可以根据上下游监测点溪水溶质浓度差来计算沿流程所发生的溶质质量转移。在此,使用PHREEQC 3.0程序获得的总无机碳(TIC)浓度沿流程的变化量来计算总的碳转移,即Δ[TIC],单位为mmol/L。由于HCO3-随pH变化,不能直接用来计算碳酸钙沉积过程中的CO2脱气。但从反应式(3)逆反应可以看出,每生成1mol的CO2,就会消耗1mol的Ca2+,所以对于碳酸钙沉淀引发的CO2脱气(Δ[CO2]p),可根据Ca2+浓度差值(Δ[Ca2+])计算。由于溪水Mg2+离子含量较低,仅为0.2mmol/L,且沿流程变化很小[14],故其对CO2脱气的贡献可以忽略不计。最后,利用公式(4)计算出物理溶解的CO2的脱气量(Δ[CO2]d)。
Δ[CO2]d=Δ[TIC]−Δ[CO2]p=Δ[TIC]−2×Δ[Ca2+] (4) 表 3列出各渠段的CO2转移量计算结果。可以看出自泉口C1点至C14点,在约270m的距离内,物理溶解的CO2脱气量最大,Δ[CO2]d=0.67mmol/L,同时碳酸钙沉积引发的CO2脱气量也较大,Δ[CO2]p=0.20mmol/L。而在下游平原区,如C25~C26段,流程约1.5km,两种类型CO2脱气量都很小。
表 3. 各渠段溪水CO2脱气变化幅度Table 3. Downstream change of CO2 outgassing flux along the stream water in different channel sections渠段 Δ[TIC]
(mmol/L)Δ[CO2]p
(mmol/L)Δ[CO2]d
(mmol/L)[DIC]s
(mmol/L)Pp Pd C1~C14 1.07 0.20 0.67 2.99 7% 22% C14~C18 0.08 0.10 - - - - C18~C19 0.16 0.05 0.06 - - - C21~C23 0 - - - - - C25~C26 0.12 0.05 0.02 - - - 本文进一步计算了C1~C14段溪水CO2脱气占泉水中DIC来自土壤CO2物理溶解和参与碳酸盐溶蚀的部分,即[DIC]s。从泉水pH、Ca2+浓度、pCO2和δ13CDIC等指标可以看出含水层内的碳酸盐岩溶蚀作用是在开放条件下进行的。因此,可以根据反应式(3)的平衡关系,用PHREEQC 3.0程序计算出的DIC浓度减去钙的浓度来确定[DIC]s值,约为2.99mmol/L,约占泉水总DIC的55%,这与蒲俊兵[34]计算的岩溶水体DIC中来自土壤CO2的比例一致。最后,计算出两种类型CO2脱气分别占[DIC]s的比例(Pd和Pp),见表 3。可以发现自泉口C1点至C14点,溪水中总的CO2脱气占了碳酸盐岩溶蚀过程中吸收固定CO2的29%,主要是物理溶解的土壤CO2返回了大气,而碳酸钙沉淀产生的CO2释放(主要消耗HCO3-)相对较小,仅造成了7%的CO2释放(表 3)。本研究还模拟了假设泉水补给溪流水在与大气pCO2(~10-3.5atm)达到完全平衡时(无碳酸钙沉淀),物理溶解的CO2沿流程脱气所占的比例。结果显示,最多会有48.9%的CO2会释放回大气中。但是,事实上这种扩散脱气达到完全平衡则需要很长时间和更长的距离。
3.3.2 溪水CO2脱气影响因素及对碳汇稳定性的启示
岩溶泉补给溪流CO2脱气受到许多因素的影响,如泉水水化学、温度、流量、流速、地形及生物作用等[13-17]。在长流水溪流约2.7km长的距离内,溪水CO2脱气量变化与地形决定的水动力条件密切相关。在C1~C14段,坡度较陡,水流厚度较薄,且多呈紊流,这样的水动力条件有利于CO2快速脱气,并使得溪水很快获得较高的碳酸钙饱和度,继而发生碳酸钙碳酸钙沉积。其原理为:水中CO2脱气时间常数受水流厚度控制[27],故水层越薄,CO2脱气越快。紊流状态(即气液界面大)有利于CO2快速脱气也同此理。Wang等[14]研究发现,当泉水具有较高的流量时,渠道中水层厚度增大,CO2脱气慢,从而导致碳酸钙饱和度升高也较慢,进而限制了碳酸钙沉积的发生。另外,前人研究发现,在流动水体中,固体表面和水溶液间存在有一层扩散边界层(DBL),该DBL厚度变化会影响Ca2+和HCO3-的传输,从而影响碳酸钙沉积速率;而较高的流速可以降低DBL的厚度,从而促进碳酸钙沉淀[35-36]。因此,在C1~C14段较高水流速度也可能是引发较多碳酸钙沉积的原因。而在平缓渠段,如C25~C26段,水流速度较慢,多呈层流,导致了CO2脱气较慢,碳酸钙沉积也很少。
这些发现也启示我们,在地势低缓地区、受流量很大的岩溶泉/地下河补给的溪/河流中,CO2脱气作用应该较弱,尤其是碳酸钙沉积诱发的CO2脱气有限。此外,我们还观测到当两股具有不同pCO2水流混合后,水化学会被稀释,碳酸钙饱和度会变得更低,这将会使得CO2脱气及碳酸钙沉积作用大大减弱。在长流水溪流汇入良丰河支流(其流量是溪流的数十倍),再到良丰河,最后进入漓江,水中DIC不断被稀释而变得更稳定[21, 37]。中国西南地区分布着许多流量非常大的岩溶泉水和地下河(流量通常>100L/s),其中很多地下河直接排入较大的河流、水库和湖泊中[19, 38]。在这种水文条件下,CO2脱气量相对较小。Zhang等[22]研究甚至发现了岩溶地下水的输入会对河流水体中的內源CO2脱气具有缓冲效应,从而减少了水体CO2排放。
在源头区,坡度较大,水流速度快,故水中生物代谢作用对水化学的影响有限,但是在下游较缓的河段,流速较慢,大量水生生物生长可能会对河水DIC产生影响。许多实验室模拟实验和野外研究结果表明,岩溶水体中水生光合生物可以吸收利用DIC(包括CO2和HCO3-)进行新陈代谢,并转化为有机碳[39-41],一些不易分解的惰性有机碳[42]能够被搬运或沉降在河流、水库、湖泊和海洋底部,这种效应被称为“生物碳泵”[43]。生物碳泵效应的存在也使得陆地碳酸盐风化作用至少在短时间范围内能够代表一种净碳汇[5-6, 44]。因此,在评估和模拟全球碳收支平衡时有必要考虑岩溶碳汇的贡献。
4. 结论
本文对广西桂林长流水一表层岩溶泉补给的地表溪流(约2.7km)进行了水化学和同位素分析测试,根据水化学和溶解无机碳碳同位素指标沿流程的变化,揭示了岩溶地下水出露地表后存在两种不同类型的CO2脱气:一是,地下水中来自土壤的物理溶解的CO2脱气;二是,渠道内碳酸钙沉淀引发的CO2脱气。这两种类型CO2脱气都显著受到了地形决定的水动力条件影响,在泉口下方的陡坡渠段,水层较薄、流速较快,溪水CO2脱气强烈,并伴随碳酸钙沉淀;而在下游平缓地段,水动力条件弱,使得溪水CO2脱气也较弱。此外,支流混合稀释会改变溪水地球化学特征,从而在一定程度上抑制CO2脱气。岩溶地下水补给溪流CO2脱气能够抵消部分岩溶作用吸收固定的大气CO2量,但是在长流水这一高地势、低流量且有碳酸钙沉积的环境下,其抵消的量也仅占29%,其中碳酸盐沉积仅造成了7%的CO2返回大气。
长流水溪流CO2脱气通量及影响因素研究结果启示我们,对于在低缓地区、受流量很大的岩溶泉/地下河补给的河流,其CO2脱气作用对岩溶碳汇的影响应该不大,尤其是碳酸钙沉淀引发的CO2脱气有限,加之受增强的水生光合生物固碳效应影响,岩溶碳汇具有很高的稳定性。这些发现为岩溶作用碳汇效应提供了更多证据,岩溶碳汇有潜力在全球碳收支平衡和未来陆地增汇中作出重要贡献。
要点
(1) 利用水化学和同位素技术分析了岩溶泉补给溪流二氧化碳脱气作用。
(2) 溪流二氧化碳脱气作用及通量主要受地形决定的水动力条件影响。
(3) 溪流碳酸钙沉积会造成岩溶碳汇不稳定,但在地势平缓区其影响有限。
HIGHLIGHTS
(1) Outgassing of carbon dioxide in a karst spring-fed stream was explored based on hydrochemical and isotopic analyses.
(2) Process and flux of carbon dioxide outgassing in streamwere mainly affected by topographically controlled hydrological conditions.
(3) Carbonate precipitation in streams caused instability of karst-related carbon sinks, but such impact was limited in streams from low-relief areas.
-
图 3 狮吼山硫铁多金属矿床成矿流体氢氧同位素组成图解(底图据路远发[22])
Figure 3.
表 1 狮吼山硫铁多金属矿床矿石硫化物的硫同位素组成
Table 1. Sulfur isotopic compositions of ore sulfides from the Shihoushan pyrite polymetallic deposit
样品编号 采样对象 测试矿物 采样位置 δ34SVCDT(‰) 数据来源 KD-b1 致密块状磁黄铁矿矿石 磁黄铁矿 160中段V1 -1.99 KD-b2 致密块状磁黄铁矿矿石 磁黄铁矿 160中段V1 -0.87 本文 KD-b4 团块状磁黄铁矿矿石 磁黄铁矿 160中段V1 -1.18 SHS-05 团块状磁黄铁矿矿石 磁黄铁矿 160中段V11 -1.30 SHS-06 团块状磁黄铁矿矿石 磁黄铁矿 160中段V11 -2.50 文献[10] SHS-07 团块状磁黄铁矿矿石 磁黄铁矿 160中段V11 -2.10 SHS2#-5 团块状磁黄铁矿矿石 磁黄铁矿 160中段V11 -1.90 SHS2#-2 磁黄铁矿-黄铁矿矿石 黄铁矿 160中段V11 -0.20 SHS2#-6a 磁黄铁矿-黄铁矿矿石 黄铁矿 160中段V11 -2.10 文献[10] SHS3#-2 磁黄铁矿-黄铁矿矿石 黄铁矿 160中段V11 -5.20 SHS3#-3 磁黄铁矿-黄铁矿矿石 黄铁矿 160中段V11 -5.50 表 2 狮吼山硫铁多金属矿床氢氧同位素组成
Table 2. Hydrogen and oxygen isotopic compositions of the Shihoushan pyrite polymetallic deposit
样品编号 矿石类型 测试矿物 采样位置 δDv-SNOW(‰) δ18Ov-pdb(‰) δ18Ov-SNOW(‰) t(℃) δ18OH2O(‰) 数据来源 KD-b4 团块状磁黄铁矿 石英 160中段V1 -74.4 -21.5 8.8 360 3.76 Shs-Q01 团块状磁黄铁矿 石英 160中段V1 -69.0 -18.7 11.6 360 6.56 Shs-Q02 团块状磁黄铁矿 石英 160中段V1 -71.3 -17.1 13.3 360 8.26 本文 Shs-Q03 团块状磁黄铁矿 石英 160中段V1 -63.7 -17.5 12.8 360 7.76 Shs-Q05 团块状磁黄铁矿 石英 160中段V1 -66.3 -14.5 15.9 360 10.86 SHS2#-2 团块状磁黄铁矿 石英 160中段V11 -55.0 / 9.4 360 4.36 SHS2#-4 团块状磁黄铁矿 石英 160中段V11 -49.0 / 13 360 7.96 文献[10] SHS2#-5 团块状磁黄铁矿 石英 160中段V11 -57.0 / 11.5 360 6.46 SHS3#-1 团块状磁黄铁矿 石英 160中段V11 -48.0 / 13.7 360 8.66 表 3 不同成因类型的硫铁矿床成矿流体特征对比
Table 3. Camparison of characteristics of ore fluids from different type pyrite deposits
特征 原生岩浆硫铁 变质流体 硫铁矿床 骆驼山式硫铁矿 向山式硫铁矿 狮吼山式硫铁矿 包裹体类型 液相和含子矿物包裹体 气液包裹体和CO2包裹体 液相为主,少量为气相CO2包裹体 液相包裹体 富气相包裹体和含子晶包裹体 流体体系 H2O-NaCl H2O-NaCl±CO2±CH4 H2O-NaCl-CO2 H2O-NaCl H2O-NaCl±CO2 温度(℃) >400 >200 370~520 90~270 180~500 盐度(wt%NaCl) 变化范围大 一般中低盐度 早期低→晚期高 8.01~21.1 1.05~4.24 δ34S(‰) -3~5 -20~22 0.24~6.46 10~23 -5.50~1.68 δ18O水(‰) 5~7 -16~25 -0.03~1.93 -11.01~1.60 3.76~10.86 δ18D水(‰) -80~-50 -140~-20 -85‰~-80‰ -51.6~-72.0 -74.4~-48.0 流体来源 岩浆水 变质水 以岩浆流体为主,晚期混入大气水 大气降水和热卤水 岩浆水与变质水混合成矿,晚期混入大气水 成矿类型 / / 岩浆热液交代型 热水沉积型 岩浆期后热液充填交代型 数据来源 徐兆文等[21]和张德会等[14] 徐兆文等[21]和张德会等[14] 梁新辉[23]和邢矿[24] 熊先孝等[25] 本文和赵正[10] -
[1] 王登红, 陈毓川, 陈郑辉, 等.南岭地区矿产资源形势分析和找矿方向研究[J].地质学报, 2007, 81(7):882-890. doi: 10.3321/j.issn:0001-5717.2007.07.002
Wang D H, Chen Y C, Chen Z H, et al.Assessment on mineral resource in Nanling Region and suggestion for further prospecting[J].Acta Geologica Sinica, 2007, 81(7):882-890. doi: 10.3321/j.issn:0001-5717.2007.07.002
[2] 王臻, 赵芝, 邹新勇, 等.赣南浅变质岩岩石地球化学特征及稀土成矿潜力研究[J].岩矿测试, 2018, 37(1):96-107. http://www.ykcs.ac.cn/article/doi/10.15898/j.cnki.11-2131/td.201710160167
Wang Z, Zhao Z, Zou X Y, et al.Petrogeochemical characteristics and metallogenetic potential of epimetamorphic rock in South Jiangxi Province[J].Rock and Mineral Analysis, 2018, 37(1):96-107. http://www.ykcs.ac.cn/article/doi/10.15898/j.cnki.11-2131/td.201710160167
[3] 陈毓川, 陈郑辉, 曾载淋, 等.南岭科学钻探第一孔选址研究[J].中国地质, 2013, 40(3):659-670. doi: 10.3969/j.issn.1000-3657.2013.03.001
Chen Y C, Chen Z H, Zeng Z L, et al.Research on the site selection of Nanling scientific drilling-1[J].Geology in China, 2013, 40(3):659-670. doi: 10.3969/j.issn.1000-3657.2013.03.001
[4] 曾载淋.南岭东段银坑地区主要金属矿产成矿规律研究与深部找矿探索[D].北京: 中国地质科学院, 2012.
Zeng Z L.Mineralization Regularity and Deep Prospecting of the Yinkeng Area in the East Nanling Metallognic Belt[D].Beijing: Chinese Academy of Geological Sciences, 2012.
[5] Feng C Y, Li D X, Zeng Z L, et al.Mineralization styles and genesis of the Yinkeng Au-Ag-Pb-Zn-Cu-Mn polymetallic orefield, Southern Jiangxi Province, SE China:Evidence from geology, fluid inclusions, isotopes and chronology[J].Acta Geologica Sinica, 2014, 88(3):825-844. doi: 10.1111/acgs.2014.88.issue-3
[6] 李伟, 刘翠辉, 贺根文, 等.赣南银坑矿田土壤与汞气地球化学特征与成矿预测[J].矿物岩石地球化学通报, 2016, 35(6):1195-1202. doi: 10.3969/j.issn.1007-2802.2016.06.010
Li W, Liu C H, He G W, et al.Geochemical anomalies and metallogenic prediction in the Yinkeng ore field, Southern Jiangxi[J].Bulletin of Mineralogy, Petrology and Geochemistry, 2016, 35(6):1195-1202. doi: 10.3969/j.issn.1007-2802.2016.06.010
[7] 赵正, 陈毓川, 曾载淋, 等.江西银坑W-Ag-Au多金属矿田成矿规律与找矿方向:兼论华南两个成矿系列叠加问题[J].地学前缘, 2017, 24(5):54-61. http://d.old.wanfangdata.com.cn/Periodical/dxqy201705006
Zhao Z, Chen Y C, Zeng Z L, et al.Jiangxi Yinkeng W-Ag-Au ore field's metallogenic regularity and prospecting direction:As well as the superposition of two metallogenic series in Southern China[J].Earths Science Frontiers, 2017, 24(5):54-61. http://d.old.wanfangdata.com.cn/Periodical/dxqy201705006
[8] Zhao Z, Zhou X P, Guo N X, et al.Superimposed W and Ag-Pb-Zn(-Cu-Au) mineralization and deep prospec-ting:Insight from a geophysical investigation of the Yinkeng orefield, South China[J].Ore Geology Reviews, 2018, 93:404-412. doi: 10.1016/j.oregeorev.2018.01.017
[9] 梁景时, 漆富勇, 范会虎, 等.江西省宁都县青塘硫铁矿床地质特征及成因探讨[J].西部探矿工程, 2014(6):169-173. doi: 10.3969/j.issn.1004-5716.2014.06.055
Liang J S, Qi F Y, Fan H H, et al.Geological characteristics and genesis of Qingtang pyrite deposits in Ningdu, Jiangxi[J].West-China Exploration Engineering, 2014(6):169-173. doi: 10.3969/j.issn.1004-5716.2014.06.055
[10] 赵正.南岭东段银坑矿田构造-岩浆活动与成矿规律研究[D].北京: 中国地质科学院, 2012: 124-136.
Zhao Z.Tectono-Magmatic Activity and Mineralization Regularity of Yinkeng Ore-field in the East Nanling Metallognic Belt[D].Beijing: Chinese Academy of Geological Sciences, 2012: 124-136.
[11] Ohmoto H.Systematics of sulfur and carbon isotopes in hydrothermal ore deposits[J].Economic Geology, 1972, 67(5):551-578. doi: 10.2113/gsecongeo.67.5.551
[12] 陕亮, 郑有业, 许荣科, 等.硫同位素示踪与热液成矿作用研究[J].地质与资源, 2009, 18(3):197-203. doi: 10.3969/j.issn.1671-1947.2009.03.007
Shan L, Zheng Y Y, Xu R K, et al.Review on sulfur isotopic tracing and hydrothermal metalogenesis[J].Geology and Resources, 2009, 18(3):197-203. doi: 10.3969/j.issn.1671-1947.2009.03.007
[13] 王艳娟, 胡援越, 申俊峰, 等.太行山南段北洺河铁矿S、Pb同位素组成及其对成矿物质来源的示踪[J].现代地质, 2011, 25(5):846-852. doi: 10.3969/j.issn.1000-8527.2011.05.003
Wang Y J, Hu Y Y, Shen J F, et al.Sulfur and lead isotope composition and tracing for sources of ore-forming materials in Beiming River iron deposits, Sourthern Taihang Mountains[J].Geoscience, 2011, 25(5):846-852. doi: 10.3969/j.issn.1000-8527.2011.05.003
[14] 张德会, 金旭东, 毛世德, 等.成矿热液分类兼论岩浆热液的成矿效率[J].地学前缘, 2011, 18(5):90-102. http://d.old.wanfangdata.com.cn/Periodical/dxqy201105010
Zhang D H, Jin X D, Mao S D, et al.The classification of ore-forming fluid and the efficiency of ore formation of magmatic hydrothermal solution[J].Earth Science Fronties, 2011, 18(5):90-102. http://d.old.wanfangdata.com.cn/Periodical/dxqy201105010
[15] 郭春影, 张文钊, 葛良胜, 等.氢氧同位素体系成矿流体示踪若干问题[J].矿物岩石, 2011, 31(3):41-47. doi: 10.3969/j.issn.1001-6872.2011.03.007
Guo C Y, Zhang W Z, Ge L S, et al.Several questions on tracing ore forming fluid by using hydrogen and oxygen isotope system[J].Journal of Mineralogy and Peteology, 2011, 31(3):41-47. doi: 10.3969/j.issn.1001-6872.2011.03.007
[16] Zhou J X, Huang Z L, Bao G P, et al.Sources and thermo-chemical sulfate reduction for reduced sulfur in the hydrothermal fluids, southeastern SYG Pb-Zn metallogenic province, SW China[J].Journal of Earth Science, 2013, 24(5):759-771. doi: 10.1007/s12583-013-0372-8
[17] Chaussidon M, Lorand J P.Sulphur isotope composition of orogenic spinel lherzolite massifs from Ariege (North-Eastern Pyrenees, France):An ion microprobe study[J].Geochimica et Cosmochimica Acta, 1990, 54(10):2835-2846. doi: 10.1016/0016-7037(90)90018-G
[18] 陈杰, 段士刚, 张作衡, 等.新疆西天山式可布台铁矿地质、矿物化学和S同位素特征及其对矿床成因的约束[J].中国地质, 2014, 41(6):1833-1852. doi: 10.3969/j.issn.1000-3657.2014.06.006
Chen J, Duan S G, Zhang Z H, et al.Geology, mineral chemistry and sulfur isotope geochemistry of the Shikebutai iron deposit in West Tianshan Mountains, Xinjiang:Constraints on genesis of the deposit[J].Geology in China, 2014, 41(6):1833-1852. doi: 10.3969/j.issn.1000-3657.2014.06.006
[19] 燕长海, 刘国印.豫西南铅锌多金属矿控矿条件及找矿方向[J].地质通报, 2004, 23(11):1143-1148. doi: 10.3969/j.issn.1671-2552.2004.11.016
Yan C H, Liu G Y.Metallogenic characteristics and ore-prospecting suggestions of lead-zinc polymetallic deposits in Southwestern Henan Province of China[J].Geological Bulletin of China, 2004, 23(11):1143-1148. doi: 10.3969/j.issn.1671-2552.2004.11.016
[20] 李占轲.华北克拉通南缘中生代银-铅-锌矿床成矿作用研究[D].北京: 中国地质大学, 2013: 30-31.
http://cdmd.cnki.com.cn/Article/CDMD-10491-1013354258.htm Li Z K.Metallogenesis of the Silver-Lead-Zinc Deposits along the Southern Margin of the North China Craton[D].Beijing: China University of Geoscience, 2013: 30-31.
[21] 徐兆文, 方长泉, 蒋少涌, 等.安徽青阳峙门口层状硫铁矿矿床地质特征及成因研究[J].地质找矿论丛, 2006, 21(1):10-14. doi: 10.3969/j.issn.1001-1412.2006.01.002
Xu Z W, Fang C Q, Jiang S Y, et al.Study on geological characteristics and genesis of stratified pyritic deposit in Shimenkou, Anhui Province[J].Contributions to Geology and Mineral Resources, 2006, 21(1):10-14. doi: 10.3969/j.issn.1001-1412.2006.01.002
[22] 路远发.Geokit:一个用VBA构建的地球化学工具软件包[J].地球化学, 2004, 33(5):459-464. doi: 10.3321/j.issn:0379-1726.2004.05.004
Lu Y F.Geokit-A geochemical toolkit for Microsoft excel[J].Geochimica, 2004, 33(5):459-464. doi: 10.3321/j.issn:0379-1726.2004.05.004
[23] 梁新辉.东秦岭骆驼山多金属硫铁矿地质特征与控矿模式探讨[D].北京: 中国地质大学, 2015: 40-42.
http://cdmd.cnki.com.cn/Article/CDMD-11415-1015385564.htm Liang X H.Geological Characteristics of the Luotuoshan Polymetallic Pyrite Deposit and Its Ore-controlling Model in Eastern Qinling, China[D].Beijing: China University of Geoscience, 2015: 40-42.
[24] 邢矿.豫西南栾川群地层特征及其与铅锌矿成矿关系研究[D].北京: 中国地质大学, 2005: 45-52.
http://cdmd.cnki.com.cn/Article/CDMD-11415-2005155143.htm Xing K.Study on the Stratigraphic Characters of the Luanchuan Group and the Relation to Pb-Zn Deposit Mineralization in the Southwest of Henan Province[D].Beijing: China University of Geoscience, 2005: 45-52.
[25] 熊先孝, 姚超美.宁芜向山式硫铁矿地质及矿床成因新探[J].化工矿产地质, 2001, 23(2):93-100. doi: 10.3969/j.issn.1006-5296.2001.02.005
Xiong X X, Yao C M.New thoughts on formation of Xiangshan-type pyrite deposit in Nanjing-Wuhu region[J].Geology of Chemical Minerals, 2001, 23(2):93-100. doi: 10.3969/j.issn.1006-5296.2001.02.005
-