非岩溶水和硫酸参与溶蚀对湘南地区地下河流域岩溶碳汇通量的影响
黄奇波1,2,3, 覃小群2,3,*, 刘朋雨2,3, 张连凯2,3, 苏春田1,2,3
1.中国地质大学环境学院,湖北 武汉 430074
2.中国地质科学院岩溶地质研究所, 广西 桂林 541004
3.国土资源部广西岩溶动力学重点实验室,广西 桂林 541004
*通信作者:覃小群(1961-),女,广西梧州人,研究员,主要从事岩溶水文地质研究.E-mail:qxq@karst.ac.cn

作者简介:黄奇波(1982-),男,江西萍乡人,副研究员,主要从事岩溶水文地质、岩溶环境与全球变化研究.E-mail:qbohuang@karst.ac.cn

摘要

研究非岩溶水和硫酸参与溶蚀对地下河流域岩溶碳汇通量的影响,有助于提高岩石风化碳汇通量估算精度,对于推进地质作用与全球气候变化研究意义重大。选取湘南北江上游武水河流域内4条典型地下河为对象,通过水化学对比分析,揭示硅酸盐岩风化对流域地下水化学的重要影响。运用Galy方法计算流域非岩溶地层中的硅酸盐岩风化消耗大气/土壤CO2对岩石风化碳汇的重要贡献,并评价了H2SO4参与下碳汇通量的扣除比例。结果显示:①流域内有非岩溶地层的L01,L02地下河, Na+,K+和SiO2浓度明显高于纯碳酸盐L03和L04地下河,非岩溶地层中的硅酸盐的风化对地下河水中K+,Na+,SiO2浓度有一定贡献;② 4条地下河的[Ca2++Mg2+]/[HCO3-]当量比值为1.05~1.15,[Ca2++Mg2+]/[HCO3-+SO42-]的当量比值为0.99~1.08,Ca2++Mg2+相对于HCO3-过量,过量的Ca2++Mg2+与SO42-相平衡,证实硫酸参与流域碳酸盐岩的溶蚀;③L01和L02地下河岩石风化消耗的CO2通量中非岩溶地层中的硅酸盐风化消耗所占比例分别为3.36%和2.22%,而L03和L04地下河中硅酸盐风化消耗比例小于0.50%,表明有非岩溶地层存在的地下河流域,其岩石风化消耗的CO2通量中硅酸盐风化消耗占有一定比例;④在考虑硫酸参与碳酸盐岩溶蚀时,4条地下河的碳汇通量分别扣除4.84%,4.52%,6.20%和9.36%。

关键词: 地下河; 碳酸盐岩风化; 硅酸盐岩风化; 碳汇通量; 硫酸
中图分类号:P641.134 文献标志码:A 文章编号:1001-8166(2017)03-0307-12
The Influence of Allogenic Water and Sulfuric Acid to Karst Carbon Sink in Karst Subterranean River in Southern Hu’nan
Huang Qibo1,2,3, Qin Xiaoqun2,3, Liu Pengyu2,3, Zhang Liankai2,3, Su Chuntian1,2,3
1.Environmental School,China University of Geosciences, Wuhan 430074, China
2.Institute of Karst Geology, Chinese Academy of Geological Sciences,Guilin 541004, China
3.Karst Laboratory of Karst Dynamics, Ministry of Land and Resources & Guangxi Zhuang Autonomous Region, Guilin 541004, China

*Corresponding author:Qin Xiaoqun(1961-), female, Wuzhou City, Guangxi Zhuang Autonomous Region, Professor. Research areas include karst geo-hydrology.E-mail:qxq@karst.ac.cn

First author:Huang Qibo(1982-), male, Pingxiang City, Jiangxi Province, Associate Professor. Research areas include karst geo-hydrology and Karst environment with global chang.E-mail:qbohuang@karst.ac.cn

Abstract

Evaluating the impact of allogenic water and sulfuric acid on karst carbon sink not only helps to improve the accurate calculation of soil CO2 uptake by rock weathering, but also obtains a complete understanding of the global carbon cycle. Groundwater samples were collected from four karst subterranean rivers watershed within different lithology strata in Wushui Basin, upstream of Beijiang Basin, Hunan Province, for revealing the important impact of silicate weathering on hydrochemistry of groundwater. To estimate the contribution of soil CO2 uptake by silicate weathering to CO2 uptake by rock weathering, the Galy model was employed in this article. The important impact of sulfuric acid on CO2 uptake by carbonate weathering resulting from the substitution of carbonic acid by protons from sulfuric acid was investigated. Our results showed that the concentration of Na+, K+ and SiO2 in L01,L02 subterranean river with silicate strata in watershed were higher than that in L03,L04 subterranean river without silicate strata in watershed, which implied that the contribution of silicate weathering to Na+,K+ and SiO2 was very important in watershed within silicate strata . The changeable equivalent ratio between (Ca2++Mg2+) and HCO3- was 1.05 to 1.15, and the value of [Ca2++Mg2+]/[HCO3-+SO42-] was 0.99 to 1.08. The concentrations of Ca2+ and Mg2+ exceeded the equivalent concentrations of HC3-, and the excess of Ca2+ and Mg2+cations were compensated by SO42-, which suggested that sulfuric acid has an important influence on carbonate dissolution. The contribution of soil CO2uptake by silicate weathering to CO2 consumption in L01 and L02 subterranean river were 3.36% and 2.22%, respectively, whereas the contribution in L03, L04 subterranean river were less than 0.50%, indicating that the contribution of soil CO2 uptake by silicate weathering was important in the subterranean river basin within silicate strata. Due to the contributions made by sulfuric acid, the CO2 consumption in four subterranean rivers decreased by 4.84%, 4.52%, 6.20%, 9.36%, respectively.

Keyword: Karst subterranean rivers; Carbonate weathering; Silicate weathering; CO2 consumption flux; Sulfuric acid.

岩石化学风化能够回收大气/土壤中的CO2产生碳汇效应, 积极参与全球碳循环[1~4]。由于方解石快速的溶解反应动力学机制, 碳酸盐岩风化碳汇是岩石风化碳汇的主要组成部分[5~9]。碳酸盐岩溶蚀形成的岩溶水中 HC O3-有一半来源于大气/土壤 CO2, 一半来源于碳酸盐岩[10, 11], 计算岩溶碳汇通量较为通用的方法是根据地下水流量和水中HC O3-浓度的水化学径流法[12~14]。完整的岩溶流域内经常存在非岩溶地层或夹层[15], 流经非岩溶地层或夹层的非岩溶水补给岩溶水的同时, 也带来了一部分HC O3-, 该部分HC O3-是碳酸溶蚀硅酸盐岩所产生的, 均源于大气/土壤 C O216, 17, 造成流域出口岩溶地下水中HC O3-来源于大气/土壤 CO2的量大于50%[15]。运用水化学径流法计算岩溶碳汇通量时由于忽略了非岩溶水带来的HC O3-, 这可能会造成碳汇效应被低估[18]。另一方面, 西南岩溶区大气酸沉降较为严重, 硫酸(H2SO4)参与碳酸盐岩溶蚀已被大量的研究所证实[19~22], 硫酸参与碳酸盐岩的溶蚀并不消耗大气或土壤中的CO2, 但会增加碳酸盐岩的溶蚀量并增加水体中HC O3-的含量[23, 24], 如何扣除地下水中来源于硫酸溶蚀成因的HC O3-也成为准确计算岩溶碳汇通量的关键[25~27]

Galy方法是Galy和France-Lanord用来估算Ganges-Brahmaputra流域的侵蚀速率和岩石风化消耗大气/土壤CO2通量的重要方法[28]。该方法将河流水化学离子组分归因于不同岩石端元的溶解, 根据河水中元素比值关系, 运用水化学计量方法来估算不同岩石端员风化对河水溶质的贡献。在我国云贵高原的乌江流域、南盘江、北盘江流域, 韩贵琳等[29]和刘丛强等[30]运用Galy方法计算了硅酸盐岩和碳酸盐岩风化速率及消耗大气/土壤CO2通量, 并评估了硫酸参与碳酸盐岩溶蚀的比例及对岩溶碳汇效应的影响程度。覃小群等[31]运用Galy方法计算出珠江流域平均岩石风化速率和由岩石风化作用引起的大气CO2消耗通量, 并证实由于碳酸盐岩风化溶蚀速率是硅酸盐岩的10.8 倍, 碳酸盐岩风化是流域碳汇过程及效应的主体。张连凯等[32]利用Galy方法研究了长江流域岩石化学风化与大气CO2消耗, 并估算了在考虑硫酸参与碳酸盐岩溶蚀时, 流域岩石风化碳汇通量应扣除约14%。然而, 将Galy方法用于计算典型地下河流域岩石风化碳汇通量的文献十分少见。

本研究选取湘南北江上游武水河流域内具有不同岩性条件的4条典型地下河为对象, 通过水化学组分对比分析, 揭示非岩溶地层中的硅酸盐岩风化对地下河水化学组成的重要影响。运用Galy方法计算岩溶流域内非岩溶地层中的硅酸盐岩风化消耗大气/土壤CO2通量对岩石风化碳汇通量的重要贡献, 并评价在硫酸的参与下, 流域碳汇通量的扣除比例。本研究为存在非岩溶水补给和外源酸参与溶蚀的地下河流域如何准确计算碳汇通量提供参考, 有助于提高岩石风化碳汇通量估算精度, 对于推进地质作用与全球气候变化研究意义重大。

1 研究区概况

研究区位于湖南省南部, 为长江流域与珠江流域分水岭地带, 归属珠江水系北江流域, 属亚热带季风性湿润气候区, 因南北气流受南岭山脉影响, 形成了气候温暖、四季分明、热量充足、雨水集中、春暖多变、夏秋多旱、严寒期短、暑期热长的气候特征。区内降水具有时空分布不均、年际变化大、雨量集中、前多而涝、后少而旱的特征, 上半年占全年降雨量的63.2%, 下半年为36.8%, 降水主要集中在4~6月, 占全年降水量的43.5%。湖南省是全国严重的酸雨污染区之一, 2004 年湖南省14 个城市的降水监测结果显示, pH≤ 5.6 的酸雨样占76.4%, pH≤ 4.5 的酸雨样占39.4%, 酸雨中硝酸的含量不及硫酸的10%, 主要为硫酸型酸雨[33]

研究区为基岩山区, 主要出露泥盆系的棋子桥组(Dq)、锡矿山组(Dx)、孟公坳组(Dm), 石炭系的马栏边组(Cm)、天鹅坪组(Ct)、石凳子组(Cs)、测水组(Cc)、梓门桥组(Cz), 二叠系的栖霞组(Pq)、孤峰组(Pg), 其中泥盆系的孟公坳组(Dm)、石炭系的测水组(Cc)及二叠系的孤峰组(Pg)岩性为页岩和砂页岩, 其他地层的岩性为灰岩夹少量白云岩, 碳酸盐岩面积占80%以上(图1)。山体海拔850~900 m, 地貌特征是北高南低、西高东低, 地表岩溶形态呈丘陵和洼地组合, 裸露基岩表层石芽、溶沟、溶槽广泛发育; 洼地底部浅平, 常有漏斗及落水洞发育, 以丘顶面和坡面发育较多, 为接受降水补给的通道; 受南北向断层和向斜构造控制, 地下河发育主要呈南北向和东西向; 南北向地下河规模大, 管道发育长10~12 km; 东西向地下河向东部金江河排泄, 规模较小, 长度1~3 km(图1)。

选取流域内的L01(邓家地下河)、L02(百里地下河)、L03(大泉地下河)、L04(夏家地下河)为研究对象(图1), 4条地下河的岩性及流域面积由野外实际调查圈定, 流量为建观测站后的实测流量, 采取每5天测1次流量, 遇雨天加密到1天1次的观测方法, 观测期为2015年5月至2016年6月。流域内有非岩溶地层出露的地下河为L01和L02, 其出露面积分别占流域面积比例的17.33%和13.00%, 而L03和L04地下河流域内均为碳酸盐岩地层。4条地下河的流量季节动态特征和岩性条件见表1

2 研究方法

在2015年3月、6月、8月、11月对4条地下河进行了现场和取样回实验室分析测试。主要采用仪器现场测试、现场滴定和样品室内测试相结合的研究方法。野外现场利用德国WTW公司生产pH/Cond340i测试仪测试各样点水的pH值和水温, 分辨率分别为0.1 pH和 0.1 ℃; 使用德国Aquamerck公司生产的碱度计和硬度计测定河水的HC O3-和Ca2+质量浓度, 分辨率分别为6 mg/L和1 mg/L。样品在现场进行过滤(< 0.22 μ m Millipore滤膜)后再装入100 mL聚乙烯瓶中。用于阳离子(Ca2+, Mg2+, K+和Na+)测定的样品立即加入超纯盐酸酸化至 pH< 2, 密封于暗箱中保存; 用于测定阴离子的样品直接密封于暗箱中保存。样品带回实验室后, 用原子吸收光谱法 (Atomic Absorption Spectroscopy, AAS)测定阳离子(Mg2+, K+和Na+)含量。阴离子含量 (S O42+, Cl-和N O3-)用高效液相色谱(HPLC)仪测定。阴阳离子电荷平衡之差< 5%。

图1 研究区水地质图Fig.1 Gological sketch map of study area

采用钼酸盐光度法测量(硅钼黄法)水体中的SiO2浓度, 其原理:在酸性环境中, 钼酸铵与水中SiO反应, 生成黄色可溶的硅钼杂多酸络合物, 在一定深度范围内, 其黄色与SiO2的浓度成正比, 通过测定其吸光度并与硅标准曲线对照可求得SiO2浓度。测试精度为0.50 mg/L。测试结果如表2所示。

表1 4条地下河的流量和岩性特征 Table 1 The characteristic of lithology and discharge in four kart subterranean rivers
表2 4条地下河不同月主要离子组成 Table 2 Major ionic concentrations in different seasons from four karst subterranean rivers

图2 主要离子Piper图Fig.2 Piper diagram of ionic concentration in the samplers

3 结 果

表2可以看出, 4条地下河不同月份的pH值变化范围为6.99~7.87, 呈弱碱性, 显示出典型的岩溶地下水化学特征[34]。4条地下河6月和8月的水温大于3月和11月, 年均水温分别为17.7 ℃, 18.2 ℃, 17.9 ℃和17.4 ℃, 与当地年均气温基本一致。总体上, 4条地下河水温的季节动态变幅比当地气温的变幅小, 反映出岩溶地下水水温既受外部气温控制又对外部温度具有调蓄作用。4条地下河不同月份总溶解性固体(Total Dissolved Solids, TDS)含量的变化范围为193.00~259.21 mg/L, 平均分别为248.97, 215.89, 220.81, 220.91 mg/L, 远高于世界河流平均值(69 mg/L)[35], 显示了岩溶地下河流域剧烈的岩石化学侵蚀作用[36]

由水化学数据作出Piper三线图(图2), 可知4条地下河不同月份均以Ca2+, Mg2+和HC O3-为主要离子, 其他阴阳离子相对较少。水化学类型为低矿化度HCO3-Ca型水, 具有典型的“ 高钙富碳” 岩溶水特征[37]。阳离子中, 4条地下河均是Ca2+所占比重最高, Ca2+占阳离子总量的62.8%~90.7%; Mg2+次之, 占总量的7.4%~34.6%; K+和Na+离子含量很低, 两者仅占阳离子总量的1.5%~5.1%。阴离子中, 4条地下河均是HC O3-所占比重最高, 占阴离子总量的90.8%~96.2%; 其次为S O42-, 占总量的1.6%~5.2%; Cl-, N O3-含量低, 大多低于0.1 mmol/L, 因此本研究不考虑HNO3对碳酸盐岩的风化作用。4条地下河不同月份中的SiO2含量范围为4.12~7.16 mg/L, 平均浓度分别为 5.52, 6.68, 5.18, 4.96 mg/L。

利用FREECH软件计算出方解石饱和指数(SIc)和水体二氧化碳分压(pCO2)(表1)。发现4条地下河河水的SIc具有季节动态, 3月和6月SIc较大, 8月和11月较小, 其中L01, L02和L03地下河的SIc在3月和6月大于0, 8月和11月均小于0。4条地下河不同月份水体pCO2的变化范围为2 630~20 277 Pa, 显著高于大气二氧化碳分压(32 Pa), 这缘于大气降水下渗过程中溶解了更多的土壤CO2, 使岩溶地下水具有较高的HC O3-含量和水体pC O227, 38, 39。从图3可以看出, 4条地下河不同月份SIc和水体pCO2具有较强的负相关关系(R2=0.89, p< 0.05), 意味着二者的变化均受同一个因素即水中CO2浓度的控制[40]

4 讨 论
4.1 岩石风化端元鉴别

4.1.1 硅酸盐岩风化

4条地下河均是以大气降水为补给来源, 大气降水对岩溶地下水化学组分有一定影响。大气降水对地下水化学影响的研究, 通常选择 Cl-作为参照元素, 它是雨水中的主要成分且来源相对简单, 在水的循环过程中 Cl-的行为也相对保持稳定, 一般认为其主要通过降雨从海洋输送而来, 受人类活动的贡献相对较少[41]。4条地下河不同月份中Cl-和K+, Na+离子所占阴阳离子比例均小于10%, Na+与Cl-比值(Na+/Cl-)范围为0.41~1.47, 平均值为0.83, 接近标准海水的Na+/Cl-值(0.86)[42]; 且从图4可以看出, K+, Na+与Cl-存在一定的相关性, 意味着K+和Na+的来源与Cl-一致, 大部分来源于大气降水, 人类活动对4条地下河的K+和Na+影响较小; 这可从4条地下河不同月份水中S O42-和N O3-离子的浓度低得到证实。另一方面, 在0.05检测水平下, Na+与Cl-相关系数R2小于0.5, K+与Cl-的相关系数R2小于0.7, 这意味着K+和Na+的来源与Cl-又有所差异, 除大气降雨外, 应该还有其他来源。大量的研究结果显示, K+和Na+可来源于非岩溶地层中的硅酸盐岩矿物的风化[43], 如Na+可来源于石盐的溶解和硅酸盐风化; K+可来源于钾长石和云母的风化[42]。研究区L01和L02地下河流域内均有一定比例的非岩溶地层出露(表1), 其Na+(平均值为0.077 和0.052 mmol/L)和K+浓度(平均值为0.037 和0.023 mmol/L)明显高于纯碳酸盐岩L03和L04地下河(平均值为0.034, 0.031和0.017, 0.016 mmol/L)(表2), 表明地下河流域内非岩溶地层中硅酸盐的风化对地下河的K+和Na+具有一定贡献, 由于L01地下河流域中非岩溶地层的出露面积比例大于L02, 其K+和Na+浓度也明显高于L02。

图3 地下河水饱和指数(SIC)与pCO2的关系Fig.3 The correlation between SIC and pCO2 in karst subterranean river

图4 4条地下河中Cl-与Na+的关系(a)和与K+的关系(b)Fig.4 Variationsof Cl- vs. Na+ (a)and vs. K+ (b) in four karst subterranean river

地下水中SiO2主要来源于非岩溶地层中的硅酸盐的风化[41]。4条地下河不同月份的SiO2的含量与珠江[31]、乌江[26]及长江主河道[32]水体相比, 明显偏小, 主要原因是4条地下河流域主要为碳酸盐岩, 非岩溶地层比例最高的L01也不超过18%, 而珠江流域、乌江流域和长江流域非岩溶地层出露面积均大于50%。另一方面, 研究区Cl-主要来源于大气降水, 而K+和Na+除大气降雨外, 还来源于非岩溶地层中硅酸盐岩矿物的风化, 因此, 可用(Na++K+)/(Na++K++Cl-)值的大小来指示非岩溶地层中硅酸盐岩矿物风化对地下水贡献大小[34]。从图5可以看出, 4条地下河SiO2含量与(Na++K+)/(Na++K++Cl-)具有较好的正相关关系(R2=0.42, p< 0.05), 这是由于非岩溶地层中的硅酸盐岩矿物的风化使地下水中SiO2含量与(Na++K+)/(Na++K++Cl-)同步增加, 进一步证实了流域内非岩溶地层中硅酸盐岩矿物的风化对地下河的K+, Na+, SiO2具有一定贡献的认识。由于L01和L02地下河流域内有一定面积的非岩溶地层出露, 其SiO2的平均浓度(分别为5.52和5.68 mg/L)高于纯碳酸盐岩的L03和L04地下河流域(5.18和4.96 mg/L)。

图5 SiO2含量与(Na++K+)/(Na++K++Cl-)变化关系Fig.5 Variations of SiO2 vs. (Na++K+)/(Na++K++Cl-) in karst subterranean river

4.1.2 碳酸盐岩风化

Ca2+, Mg2+和HC O3-在4条地下河的阴阳离子中占有绝对优势, 碳酸盐岩的溶蚀对地下水化学组分具有控制作用[44]。开放的岩溶动力系统中, CO2是岩溶作用发生的重要驱动力[45], CO2进入岩溶动力系统促进碳酸盐岩的溶解, 表现为地下水中Ca2+, HC O3-浓度和水体pCO2增加, SIc减少。反之, 若CO2从系统逸出, 将发生碳酸盐岩的沉积, 水中Ca2+, HC O3-浓度和水体pCO2减少, SIc增加。因此, 通常情况下地下水中SIc和水体pCO2具有负相关关系[40], 图3很好地证实了这一点。从季节动态看, 4条地下河SIc的平均值和水体pCO2平均值也具有很强的对应关系(图6)。3月气温低, 降雨稀少, 土壤覆被的生产力和生物量相对较小, 土壤有机质含量较低, 土壤微生物含量较少, 土壤有机质分解产生CO2浓度低, 造成降雨下渗吸收的CO2含量少, 地下水的pCO2低, 水的侵蚀性低, SIc高; 6月随着气温回升, 降雨量增大, 土壤微生物活沃, 土壤覆被的生产力和生物量增加, 土壤有机质含量增高, 土壤CO2浓度升高, 地下水下渗能够吸收较多的土壤CO2, 水体pCO2增加, 对碳酸盐岩的溶蚀能力增加, SIc减小。8月气温最高, 生物活动最强, 地下水中pCO2达到最大, SIc最低; 11月随着气温回落, 生物活动减弱, 地下水中pCO2减小, SIc增加。

图6 地下水SIc和pCO2季节动态Fig.6 Temporal change of SIc and pCO2 in karst subterranean river

4.1.3 人为活动输入

人为活动产物可以通过废物直接排放和通过大气输入这2种途径输入水体[44]。其中, 以工业活动和农业活动的影响最为显著。地下水中N O3-主要来源于农业活动, 如大气输入、农肥、人畜的排泄物和“ 三废” [37], 而S O42-多来自于大气酸沉降和工矿业活动[46], 因此, 一般情况下, 水中的 N O3-和 S O42-浓度的变化主要反映来自农业活动和工业活动的影响。

4条地下河N O3-浓度为0.055~0.132 mmol/L, 平均值分别为0.107, 0.065, 0.098, 0.109 mmol/L, S O42-浓度为0.073~0.208 mmol/L, 平均值分别为0.122, 0.109, 0.124, 0.184 mmol/L。N O3-和 S O42-浓度远小于受酸雨和工业活动影响较大的乌江流域地表水和地下水[26], 也小于人类活动频繁的长江流域主河道水体[32]。另外, 所有样点的N O3-/Cl-与S O42-/Cl-具有较好的相关性(相关系数R2=0.69, p< 0.01)(图7), 其取值范围分别为1.1~2.6和 1.3~4.3, 也明显小于乌江上游段地表水体中的N O3-/Cl-与S O42-/Cl-[44], 表明4条地下河受农业和工业活动影响较轻。这与野外调查基本吻合, 4条地下河补给区为海拔850~900 m一级岩溶峰丛顶面, 岩溶发育, 以盲谷、落水洞发育较好(图1), 溶蚀洼地呈深盆状。人口相对稀少, 只在少量的洼地中种植玉米、花生, 因此人类生产、生活对地下水的影响有限。但研究区位于湖南南部, 湖南是受严重的酸雨污染区之一, 酸雨发生频率高, 但不同季节又存在明显差异的酸雨频率[47], 这可能是引起研究区4条地下河中不同季节S O42-差异的主要原因。

图7 4条地下河中N O3-/Cl-与S O42-/Cl-的关系Fig.7 Four variations of N O3-/Cl- vs. S O42-/Cl-in four karst subterranean river

4.2 H2SO4参与碳酸盐岩溶蚀

根据H2CO3风化碳酸盐岩的方程[13]:

Ca(1-x)MgxCO3+CO2+H2O =(1-x)Ca2++xMg2++2HC O3-

仅有H2CO3参与岩溶作用时, 溶蚀产物主要为Ca2+, Mg2+, HC O3-。此时, [Ca2++Mg2+]/[HC O3-]的当量比值为1, [S O42-]/[HC O3-]的当量比值为0。实测4条地下河的[Ca2++Mg2+]/[HC O3-]当量比值为1.05~1.15, 表明Ca2++Mg2+相对于HC O3-过量, 除HC O3-外还有其他阴离子来平衡Ca2++Mg2+图8表示了4条地下河中[Ca2++Mg2+]/[HC O3-+S O42-]的当量比值关系, 所有样品的[Ca2++Mg2+]/[HC O3-+S O42-]的当量比值为0.99~1.08, 平均值为1.04, 说明地下水中多余的Ca2+ 和Mg2+ 应当是S O42-离子来加以平衡的, 进而说明4条地下河中, 除了H2CO3参与碳酸盐岩的侵蚀以外, H2SO4也对碳酸盐岩的风化过程产生作用。

图8 地下河Ca2++Mg2+与HC O3-+ S O42-关系Fig.8 Ca2++Mg2+ vs.HC O3-+ S O42- in karst subterranean river

另一方面, 当H2SO4+H2CO3共同风化碳酸盐时, 溶蚀方程可表示为[31]:3Ca(1-x)MgxCO3+ H2SO4+H2CO3= 3(1-x)Ca2++3xMg2++S O42-+4HC O3-, 溶蚀产物主要为Ca2+, Mg2+, HC O3-, S O42-。此时, [Ca2++Mg2+]/HC O3-的当量比为1.5, [S O42-]/[HC O3-]的当量比为0.5。图9显示, 4条地下河均分布于H2CO3和H2SO4+H2CO3溶蚀碳酸盐岩2个端元组成之间, 说明除H2CO3外, H2SO4也参与了流域岩石化学风化过程。但所有样点更靠近H2CO3风化端元一侧, 意味着4条地下河主要受H2CO3风化控制, H2SO4风化贡献有限。

图9 [Ca2++Mg2+]/[HC O3-]与[S O42-]/[HC O3-]的当量比值关系Fig.9 Equivalent ratios of [ Ca2++ Mg2+]/[ HC O3-] vs.[ S O42-] /[ HC O3-] in karst subterranean rivers

4.3 不同端元风化贡献比例及CO2消耗通量

由以上分析可知, 4条地下河既受非岩溶地层中的硅酸盐岩风化影响, 又受H2SO4溶蚀的控制, 本研究利用地下河不同季度的流量(表1)和水化学数据(表2), 根据地下河元素比值之间的关系, 采用Galy 的模型[28], 并运用参考文献[31]中的公式(1)~(16)计算出4 条地下河中碳酸盐岩和非岩溶地层中硅酸盐岩风化消耗的CO2通量, 计算结果如表3所示。

表3可以看出, 4条地下河岩石风化消耗的CO2通量中硅酸盐岩风化消耗所占比例最大的为 L01地下河, 其硅酸盐岩风化消耗通量为27.93× 103 mmol/(km2· a), 占岩石风化CO2消耗通量(815.72× 103 mmol/(km2· a))的3.36%; 其次为L02地下河, 其硅酸盐岩风化消耗通量为16.13× 103 mmol/(km2· a), 占岩石风化CO2消耗通量(727.78× 103 mmol/(km2· a)) 的2.22%。最小为L03和L04地下河, 硅酸盐岩风化消耗通量分别为4.55× 103 mmol/(km2· a)和2.81× 103 mmol/(km2· a), 分别占岩石风化CO2消耗通量(1 109.09× 103 mmol/(km2· a)和859.27× 103 mmol/(km2· a))的0.41%和0.33%; 可见, 非岩溶地层中的硅酸盐岩风化消耗CO2通量和贡献比例大小受流域内非岩溶地层面积比例控制; 由于L01地下河流域内非岩溶地层面积比例最大, 为17.33%; 其次为L02, 非岩溶地层面积比例为13.00%; 2条地下河中硅酸盐风化消耗CO2通量所占比例大于2%。而L03和L04地下河流域内由于没有非岩溶地层, 其硅酸盐风化消耗CO2通量所占比例均小于0.5%。

仅考虑碳酸溶蚀作用时, L01, L02, L03和L04地下河岩石风化(碳酸盐岩和硅酸盐岩)消耗CO2速率分别为857.17× 103, 762.26× 103, 1 182.43× 103和947.98× 103 mol/(km2· a); 在硫酸与碳酸共同作用下, 4条地下河岩石风化吸收CO2速率分别为815.72× 103, 727.78× 103, 1 109.09× 103, 859.27× 103 mol/(km2· a)。可见在硫酸的参与下, 岩石风化消耗的CO2通量降低。在考虑硫酸参与碳酸盐岩溶蚀时, 4条地下河的碳汇通量分别扣除4.84%, 4.52%, 6.20%, 9.36%。但总体上, 4条地下河受硫酸影响较小, 其扣除的比例均小于10%。刘丛强等[30]的研究结果显示在硫酸参与下, 西南地区岩溶碳汇通量应扣除33%, 本研究的扣除比例明显小于该值, 这缘于本研究中的4条地下河流域受人类活动影响较小。

表3 4条地下河岩石风化CO2通量计算结果 Table 3 CO2consumption flux by rock weathering in four karst subterranean rivers
5 结 论

流域内有非岩溶地层分布的L01和L02地下河其Na+ , K+ , SiO2浓度高于流域内为纯碳酸盐岩地层的L03和L04地下河, 地下河流域内非岩溶地层中硅酸盐的风化对地下河的K+, Na+, SiO2有一定贡献。4条地下河的[Ca2++Mg2+]/[HC O3-]当量比值为1.05~1.15, [Ca2++Mg2+]/[HC O3-+S O42-]的当量比值为0.99~1.08, 表明Ca2++Mg2+相对于HC O3-过量, 过量的Ca2++Mg2+与S O42-相平衡, 证实硫酸参与流域碳酸盐岩的溶蚀。

L01和L02地下河岩石风化消耗的CO2通量中非岩溶地层中的硅酸盐风化消耗所占比例分别为3.36%和2.22%, 而L03和L04地下河中硅酸盐风化消耗比例小于0.50%。有非岩溶地层存在的地下河流域, 其硅酸盐风化消耗CO2通量对岩石风化碳汇通量有一定贡献, 且贡献比例大小受流域内出露的非岩溶地层面积比例大小控制。在考虑硫酸参与碳酸盐岩溶蚀时, 4条地下河的碳汇通量分别扣除4.84%, 4.52%, 6.20%, 9.36%。证实硫酸参与岩石风化过程, 并减少岩石风化碳汇通量。

The authors have declared that no competing interests exist.

参考文献
[1] Amiotte-Suchet P, Probst J L, Ludwig W. Worldwide distribution of continental rock lithology: Implications for the atmospheric/soil CO2 uptake by continental weathering and alkalinity river transport to the oceans[J]. Global Biogeochemical Cycles, 2003, 17(2): 1 038. [本文引用:1]
[2] Jacobson A D, Blum J D, Chamberlain C P, et al. Climatic and tectonic controls on chemical weathering in the New Zealand Southern Alps[J]. Geochimica et Cosmochimica Acta, 2003, 67(1): 29-46. [本文引用:1]
[3] Amiotte-Suchet P, Probst J L. A global model for present-day atmospheric/soil CO2consumption by chemical erosion of continental rocks (GEM-CO2)[J]. Tellus, 1995, 47(1): 273-280. [本文引用:1]
[4] Luo Weijun, Wang Shijie, Liu Xiuming. Research progresses and prospect of chimney effect about carbon cycle in the Karst cave system[J]. Advances in Earth Science, 2014, 29(12): 1 333-1 340.
[罗维均, 王世杰, 刘秀明. 喀斯特洞穴系统碳循环的烟囱效应研究现状及展望[J]. 地球科学进展, 2014, 29(12): 1 333-1 340. ] [本文引用:1]
[5] Liu Z, Dreybrodt W. [J]. Geochimica et Cosmochimica Acta, 1997, 61(14): 2 879-2 889. [本文引用:1]
[6] Liu Z, Dreybrodt W, Wang H. A new direction in effective accounting for the atmospheric CO2 budget: Considering the combined action of carbonate dissolution, the global water cycle and photosynthetic uptake of DIC by aquatic organisms[J]. Earth-Science Reviews, 2010, 99(3/4): 162-172. [本文引用:1]
[7] Liu Zaihua. New progress and prospects in the study of rock-weathering-related carbon sinks[J]. China Science Bulletin, 2012, 57(2/3): 95-102.
[刘再华. 岩石风化碳汇研究的最新进展和展望[J]. 科学通报, 2012, 57(2/3): 95-102. ] [本文引用:1]
[8] Pu Junbing, Jiang Zhongcheng, Yuan Daoxian, et al. Some opinions on rock-weathering-related carbon sinks from the IPCC fifth assessment report[J]. Advances in Earth Science, 2015, 30(10): 1 081-1 090.
[蒲俊兵, 蒋忠诚, 袁道先, . 岩石风化碳汇研究进展: 基于 IPCC 第五次气候变化评估报告的分析[J]. 地球科学进展, 2015, 30(10): 1 081-1 090. ] [本文引用:1]
[9] Liu Zaihua, Dreybrodt W, Liu Huan. Atmospheric CO2 sink: Silicate weathering or carbonate weathering?[J]. Quaternary Sciences, 2011, 31(3): 426-430.
[刘再华, Dreybrodt W, 刘洹. 大气CO2汇: 硅酸盐风化还是碳酸盐风化的贡献?[J]. 第四纪研究, 2011, 31(3): 426-430. ] [本文引用:1]
[10] Yuan D. The carbon cycle in karst[J]. Zeitschrift fur Geomorphologie Neue Folge, 1997, 108(Suppl. ): 91-102. [本文引用:1]
[11] Jiang Z C, Yuan D X. CO2 source-sink in karst processes in karst areas of China[J]. Episodes, 1999, 22(1): 33-35. [本文引用:1]
[12] Xu Shengyou, Jiang Zhongcheng. Preliminary estimate the relationship between Karstification and CO2 sources and sink[J]. Chinese Science Bulletin, 1997, 42(9): 953-956.
[徐胜友, 蒋忠诚. 我国岩溶作用与大气温室气体CO2源汇关系的初步估算[J]. 科学通报, 1997, 42(9): 953-956. ] [本文引用:1]
[13] Yuan Daoxian, Zhang Cheng. Karst dynamics theory in China and its practice[J]. Acta Geoscientica Sinica, 2008, 29(3): 355-365.
[袁道先, 章程. 岩溶动力学的理论探索与实践[J]. 地球学报, 2008, 29(3): 355-365. ] [本文引用:1]
[14] Huang Qibo, Liu Pengyu, Qin Xiaoqun, et al. The characteristics of karst carbon sink in the Guijiang Catchment[J]. Carsologica Sinica, 2011, 30(4): 437-442.
[黄奇波, 刘朋雨, 覃小群, . 桂江流域岩溶碳汇特征[J]. 中国岩溶, 2011, 30(4): 437-442. ] [本文引用:1]
[15] Huang Fen, Tang Wei, Wang Jinliang, et al. The influence of allogenic water on karst carbon sink: A case study in the Maocun Subterranean River in Guilin, China[J]. Carsologica Sinica, 2011, 30(4): 417-421.
[黄芬, 唐伟, 汪进良, . 外源水对岩溶碳汇的影响——以桂林毛村地下河为例[J]. 中国岩溶, 2011, 30(4): 417-421. ] [本文引用:2]
[16] Amiotte-Suchet P, Probst A, Probst J L. Influence of acid rain on CO2 consumption by rock weathering: Local and global scales[J]. Water Air and Soil Pollution, 1995, 85(3): 1 563-1 568. [本文引用:1]
[17] Gaillardet J, Dupré B, Louvat P, et al. Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers[J]. Chemical Geology, 1999, 159(1/4): 3-30. [本文引用:1]
[18] Kang Zhiqiang, Yuan Daoxian, Chang Yong, et al. The main controlling factor of karst carbon sequestration: About water cycle[J]. Journal of Jilin University (Earth Science Edition), 2011, 41(5): 1 542-1 547.
[康志强, 袁道先, 常勇, . 岩溶碳汇的主控因子——水循环[J]. 吉林大学学报: 地球科学版, 2011, 41(5): 1 542-1 547. ] [本文引用:1]
[19] Han G L, Liu C Q. Water geochemistry controlled by carbonate dissolution: A study of the river waters draining karst-dominated terrain, Guizhou Province, China[J]. Chemical Geology, 2004, 204(1/2): 1-21. [本文引用:1]
[20] Li S L, Calmels D, Han G L, et al. Sulfuric acid as an agent of carbonate weathering constrained by δ13CDIC: Examples from Southwest China[J]. Earth Planetary Science Letters, 2008, 270(3/4): 180-199. [本文引用:1]
[21] Li X D, Liu C Q, Harue M, et al. The use of environmental isotopic (C, Sr, S) and hydrochemical tracers tocharacterize anthropogenic effects on karst groundwater quality: A case study of the Shuicheng Basin, SW China[J]. Applied Geochemistry, 2010, 25(12): 1 924-1 936. [本文引用:1]
[22] Li X D, Liu C Q, Liu X L, et al. Identification of dissolved sulfate sources and the role of sulfuric acid in carbonate weathering using dual-isotopic data from the Jialing River, Southwest China[J]. Journal of Asian Earth Sciences, 2011, 42(3): 370-380. [本文引用:1]
[23] Jiang Yongjun, Yuan Daoxian. Geochemical tracers to characterize effects of urbanization on karst groundwater quality from Nanshan underground river system, SW China[J]. Quaternary Science, 2014, 34(5): 1 044-1 053.
[蒋勇军, 袁道先. 城市发展对岩溶地下水质影响的地球化学示踪——以重庆南山老龙洞地下河系统为例[J]. 第四纪研究, 2014, 34(5): 1 044-1 053. ] [本文引用:1]
[24] Zhang Xingbo, Jiang Yongjun, Qiu Shulan, et al. Agricultural activities and carbon cycling in karst areas in Southwest China: Dissolving carbonate rocks and CO2 sink[J]. Advances in Earth Science, 2012, 27(4): 466-476.
[张兴波, 蒋勇军, 邱述兰, . 农业活动对岩溶作用碳汇的影响: 以重庆青木关地下河流域为例[J]. 地球科学进展, 2012, 27(4): 466-476. ] [本文引用:1]
[25] Jiang Y G. The contribution of human activities to dissolved inorganic carbon fluxes in a karst underground river system: Evidence from major elements and δ13CDIC in Nand ong, Southwest China[J]. Journal of Contaminant Hydrology, 2013, 152: 1-11, doi: 10.1016/j.jcohyd.2013.05.010. [本文引用:1]
[26] Huang Qibo, Qin Xiaoqun, Liu Pengyu, et al. Influence of sulfuric acid to karst hydrochemical and δ13CDIC in the upper and middle reaches of Wu River[J]. Environmental Science, 2015, 36(9): 102-111.
[黄奇波, 覃小群, 刘朋雨, . 硫酸对乌江中上游段岩溶水化学及δ13CDIC的影响[J]. 环境科学, 2015, 36(9): 102-111. ] [本文引用:2]
[27] Huang Qibo, Qin Xiaoqun, Liu Pengyu, et al. The impact of human activities to δ13CDIC of karst groundwater and carbon sink in the upper and middle reaches of of Wu River[J]. Quaternary Science, 2016, 36(6): 1 358-1 369.
[黄奇波, 覃小群, 刘朋雨, . 人为活动对乌江中上游段岩溶地下水δ13CDIC及碳汇效应的影响[J]. 第四纪研究, 2016, 36(6): 1 358-1 369. ] [本文引用:1]
[28] Albert Galy, Christian France-Lanord. Weathering processes in the Ganges-Brahmaputra Basin and theriverine alkalinity budget[J]. Chemical Geology, 1999, 159(1/4): 31-60. [本文引用:2]
[29] Han Guilin, Liu Congqiang. Hydrogechemistry of rivers in Guizhou Province China: Constraints on crustal weathering in karst terrain[J]. Advances in Earth Science, 2005, 20(4): 394-406.
[韩贵琳, 刘丛强. 贵州喀斯特地区河流的研究——碳酸盐岩溶解控制的水文地球化学特征[J]. 地球科学进展, 2005, 20(4): 394-406. ] [本文引用:1]
[30] Liu Congqiang, Jiang Yingkui, Tao Faxiang, et al. Chemical weathering of carbonate rocks by sulfuric acid and the carbon cycling in Southwest China[J]. Geochimica, 2008, 37(4): 404-414.
[刘丛强, 蒋颖魁, 陶发祥, . 西南喀斯特流域碳酸盐岩的硫酸侵蚀与碳循环[J]. 地球化学, 2008, 37(4): 404-414. ] [本文引用:2]
[31] Qin Xiaoqun, Jiang Zhongcheng, Zhang Liankai, et al. The difference of the weathering rate between carbonate rocks and silicate rocks and its effects on the atmospheric CO2 consumption in the Pearl River Basin[J]. Geological Bulletin of China, 2015, 34(9): 1 749-1 757.
[覃小群, 蒋忠诚, 张连凯, . 珠江流域碳酸盐岩与硅酸盐岩风化对大气CO2汇的效应[J]. 地质通报, 2015 , 34(9): 1 749-1 757. ] [本文引用:3]
[32] Zhang Liankai, Qin Xiaoqun, Liu Pengyu, et al. Chemical denudation rate and atmospheric CO2 consumption by H2CO3 and H2SO4 in the Yangtze River Catchment[J]. Geological Bulletin of China, 2016, 90(8): 1 933-1 943.
[张连凯, 覃小群, 刘朋雨, . 硫酸参与的长江流域岩石化学风化与大气CO2消耗[J]. 地质学报, 2016, 90(8): 1 933-1 943. ] [本文引用:3]
[33] Zhu Xuanxiang. The present situation, cause of formation and countermeasure research of the acid rain in Hunan Province[J]. Technological Development of Enterprise, 2006, 25(10): 28-30.
[朱选祥. 湖南省酸雨现状、成因及对策研究[J]. 企业技术开发, 2006, 25(10): 28-30. ] [本文引用:1]
[34] Chen Jingsheng, He Dawei. Chemical characteristics and genesis of major ions in the Pearl River Basin[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 1999, 35(6): 786-793.
[陈静生, 何大伟. 珠江水系河水主要离子化学特征及成因[J]. 北京大学学报: 自然科学版, 1999, 35(6): 786-793. ] [本文引用:2]
[35] Meybeck M. Pathways of major elements from land to ocean through rivers[M]∥Martin J M, Burton J D, Eisma D, eds. River Inputs to Ocean Systems. New York: United Nations Press, 1981: 18-30. [本文引用:1]
[36] Li Jun, Liu Congqiang, Li Longbo, et al. The impacts of chemical weathering of carbonate rock by sulfuric acid on the cycling of dissolved inorganic carbon in Changjiang River water[J]. Geochimica, 2010, 39(4): 305-313.
[李军, 刘丛强, 李龙波, . 硫酸侵蚀碳酸盐岩对长江河水DIC循环的影响[J]. 地球化学, 2010, 39(4): 305-313. ] [本文引用:1]
[37] Pu Junbing, Yuan Daoxian, Jiang Yongjun, et al. Hydrogeochemistry and environmental meaning of Chongqing subterranean karst streams in China[J]. Advances in Water Science, 2010, 21(5): 628-636.
[蒲俊兵, 袁道先, 蒋勇军, . 重庆岩溶地下河水文地球化学特征及环境意义[J]. 水科学进展, 2010, 21(5): 628-636. ] [本文引用:2]
[38] Huang Qibo, Qin Xiaoqun, Liu Pengyu, et al. The impact of acid rain to δ13CDIC of karst groundwater and carbon sink in dry season in Guilin[J]. Earth Science—Journal of China University of Geosciences, 2015, 40(7): 1 237-1 247.
[黄奇波, 覃小群, 刘朋雨, . 酸雨对桂林枯水期岩溶地下水δ13CDIC及碳汇效应的影响[J]. 地球科学——中国地质大学学报, 2015, 40(7): 1 237-1 247. ] [本文引用:1]
[39] Huang Q B, Qin X Q, Yang Q Y, et al. Identification of dissolved sulfate sources and the role of sulfuric acid in carbonate weathering using δ13CDIC and δ34S in karst area, Northern China[J]. Environmental Earth Science, 2016, 75(1): 530-539. [本文引用:1]
[40] Pu Junbing. Dissolved inorganic carbon and stable carbon isotope in karst subterranean streams in Chongqing, China[J]. Carsologica Sinica, 2013, 32(2): 123-132.
[蒲俊兵. 重庆地区岩溶地下河水溶解无机碳及其稳定同位素特征[J]. 中国岩溶, 2013, 32(2): 123-132. ] [本文引用:2]
[41] Jiemei, An Yanling, Wu Qixin, et al. Hydrochemical characteristics and sources of Qingshuijiang River Basin at wet season in Guizhou Province[J]. Environmental Science, 2015, 36(5): 1 565-1 572.
[吕婕梅, 安艳玲, 吴起鑫, . 贵州清水江流域丰水期水化学特征及离子来源分析[J]. 环境科学, 2015, 36(5): 1 565-1 572. ] [本文引用:2]
[42] Wang Peng, Shang Yingnan, Shen Licheng, et al. Characteristics and evolution of hydrochemical compositions of freshwater lake in Tibetan Plateau[J]. Environmental Science, 2013, 34(3): 874-881.
[王鹏, 尚英男, 沈立成, . 青藏高原淡水湖泊水化学组成特征及其演化[J]. 环境科学, 2013, 34(3): 874-881. ] [本文引用:2]
[43] Chen Jingsheng, Wang Feiyue, Xia Xinghui. Geochemistry of water quality of the Yangtze River Basin[J]. Earth Science Frontiers, 2006, 13(1): 74-85.
[陈静生, 王飞越, 夏星辉. 长江水质地球化学[J]. 地学前缘, 2006, 13(1): 74-85. ] [本文引用:1]
[44] Huang Qibo, Qin Xiaoqun, Liu Pengyu, et al. Major ionic features and their controlled factor in the Upper- Middle Reaches of Wu River[J]. Environmental Science, 2016, 37(5): 1 779-1 787.
[黄奇波, 覃小群, 刘朋雨, . 乌江中上游段河水主要离子化学特征及控制因素[J]. 环境科学, 2016, 37(5): 1 779-1 787. ] [本文引用:3]
[45] Liu Zaihua. Two important atmospheric CO2 sinks[J]. Chinese Science Bulletin, 2000, 45(21): 2 348-2 351.
[刘再华. 大气CO2两个重要的汇[J]. 科学通报, 2000, 45(21): 2 348-2 351. ] [本文引用:1]
[46] Zhang Dong, Huang Xingyu, Li Chengjie. Sources of riverine sulfate in Yellow River and its tributaries determined by sulfur and oxygen isotopes[J]. Advances in Water Science, 2013, 24(3): 418-426.
[张东, 黄兴宇, 李成杰. 硫和氧同位素示踪黄河及支流河水硫酸盐来源[J]. 水科学进展, 2013, 24(3): 418-426. ] [本文引用:1]
[47] Xu Xinhui, Gao Hongwen. Analysis of the distribution and causes of acid rain in southern China[J]. Sichuan Environment, 2011, 30(4): 135-139.
[许新辉, 郜洪文. 中国南方酸雨的分布特征及其成因分析[J]. 四川环境, 2011, 30(4): 135-139. ] [本文引用:1]