广西科学  2018, Vol. 25 Issue (5): 532-543   PDF    
漓江干流水体主要离子化学变化特征及影响因素
张陶, 蒲俊兵, 李建鸿, 吴飞红, 袁道先     
中国地质科学院岩溶地质研究所,自然资源部/广西壮族自治区岩溶动力学重点实验室,广西桂林 541004
摘要: 【目的】 为保护漓江流域水环境及揭示区域碳酸盐岩风化过程提供基础数据及科学支撑。【方法】 通过收集一个水文年漓江干流典型断面的水样,并进行水化学分析来研究水体中离子化学变化特征及其控制因素。【结果】 漓江两个断面水化学类型都为Ca-HCO3型,河水主要离子组分受到水岩作用的控制。两个断面对比来看,上游大面圩(DM)断面受到部分降雨输入的影响,而下游普益(PY)断面受到更加强烈的碳酸盐岩风化作用的影响。漓江两个断面HCO3-、Ca2+离子浓度表现出雨季小于旱季的现象,主要受到雨季降雨和流量增大的影响,而SO42-、NO3-、Cl-离子浓度出现雨季大于旱季的现象,则主要受到农业活动、城市排污的影响。此外,由于人类活动的影响,漓江碳酸盐岩风化受到外源酸的影响,进而影响岩溶风化碳汇效应的评估。【结论】 漓江典型断面水体优势离子(HCO3-、Ca2+)变化特征主要受到地质条件的控制,岩溶作用强度表现出显著的空间差异性;SO42-、NO3-、Cl-离子则受到不同人类活动的影响,表现出显著的季节变化和空间变化。
关键词: 漓江     水化学     离子     碳酸盐岩风化     变化    
Chemical Characteristics Changes of Major Ion and Its Influencing Factors in the Main Stream of Lijiang River
ZHANG Tao , PU Junbing , LI Jianhong , WU Feihong , YUAN Daoxian     
Key Laboratory of Karst Dynamics, MNR & Guangxi, Institute of Karst Geology, Chinese Academy of Geological Sciences, Guilin, Guangxi, 541004, China
Abstract: 【Objective】 To provide basic data and scientific support for protecting the water environment of the Lijiang River basin and revealing the weathering process of the regional carbonate rock. 【Methods】 The characteristics of ion chemical changes and their controlling factors were studied by collecting water samples from typical sections of the Lijiang River in the hydrological year and performing hydrochemistry analysis. 【Results】 The hydrochemistry types of two sections in the Lijiang River were all Ca-HCO3 type, and the main ionic components of the river were controlled by water-rock interaction. In comparison of the two sections, the upstream Damianxu(DM) section was affected by partial rainfall input, and the downstream Puyi(PY) section was affected by more intense weathering of carbonate rocks. The HCO3- and Ca2+ ion concentrations in the two sections of the Minjiang River were affected by the increase of rainfall and flow in the rainy season, which made the two ions in rainy season less than that in the dry season. But due to the influence of agricultural activities and urban sewage, the SO42-、NO3- and Cl- ions in the rainy season in the Lijiang River was larger than that in the dry season. In addition, due to the influence of human activities, weathering of carbonate rocks in the Lijiang River was affected by exogenous acids, which might in turn affect the karst weathering carbon sinks evaluation. 【Conclusion】 The variation characteristics of the dominant ions (HCO3- and Ca2+) in the typical section of the Lijiang River are mainly controlled by the geological conditions, and the intensity of the karstification shows significant spatial differences; SO42-、NO3- and Cl- ions are influenced by different human activities, showing significant seasonal and spatial changes.
Key words: Lijiang River     hydrochemistry     ion     carbonate rocks weathering     variation    
0 引言

研究意义】河流作为连接陆地与海洋两大生态系统的主要通道,是全球水循环的重要组成部分。河水溶解组分与流域岩性、地貌、气候、水文、土壤、植被覆盖和人类活动等要素密切相关。因此,通过分析河流的水化学组成,可以获取关于流域内元素的地球化学行为、岩石风化、人类活动等重要信息[1-3]。河流水化学的时空动态变化特征分析,则是了解外界环境包括自然因素和人为活动对河流影响程度、变化趋势以及变化幅度的一个重要手段[4]。【前人研究进展】我国的河流水化学研究始于20世纪60年代的乐嘉祥等[5]的研究。在这之后,关于河流化学特征的研究得到了迅速发展。陈静生等[6-7]对黄河和长江水化学变化的研究指出黄河离子化学具有显著的空间变化特征,且主要受到沉积岩风化作用和干旱气候影响下水中溶解盐的蒸发浓缩和结晶作用的控制;长江离子化学具有显著的时空变化特征,则主要受到岩石风化作用,尤其是碳酸盐矿物的溶解作用控制。王亚平等[8]表明长江水系的主要离子化学具有显著的空间变化特征,受到岩石风化作用的影响,且主要受到碳酸盐和蒸发岩矿物风化的贡献。韩贵琳等[9]在岩溶区河流的研究发现是流域盆地的地理岩性控制了河水的化学组成。于奭等[10]对珠江流域典型岩溶断面的研究表明岩溶地区河水化学特征主要受到碳酸盐岩溶解影响,同时受到人类活动不同途径的影响。孙平安等[11]对比两条地质条件差异显著的河流发现,两条河流主要离子都表现出冬高夏低的特征,且都受到流量变化的控制,此外两条河流离子化学主要受到河流地质条件的影响。前人研究表明, 河流水化学特征主要受到岩性、气候、水文特征、植被类型、地貌以及人类活动等的影响,尤其是岩性、气候的变化会使得河流水化学发生显著的时空变化特征。【本研究切入点】漓江流域属于亚热带季风气候,流域内表现出显著的旱雨季季节变化特征。漓江干流自北向南流,全长217 km,因此上下游气候差异较大。另外,漓江中游地区主要分布碳酸盐岩,其上游主要为花岗岩、碎屑岩,因而漓江流域基岩地质条件复杂,造成漓江流域特征出现巨大的空间异质性和时间动态性。【拟解决的关键问题】选择漓江干流上地质和气候条件差别显著的典型断面,利用两个断面冬夏(旱、雨)两个季节的水化学数据分析漓江流域水化学特征的时空差异性及其主要影响因素。

1 材料与方法 1.1 研究区概况

漓江是珠江流域西江水系一级支流桂江流域的中上游河段(图 1),其上游大榕江发源于广西壮族自治区桂林市兴安县西北部的越城岭主峰猫儿山,全长438 km,流域面积18 790 km2,河道平均坡降0.43‰,多年平均径流深1 033 mm,多年平均径流量1.44×1010 m3。漓江自北向南流,全长217 km,因此上下游气候差异较大。漓江流域属于亚热带季风气候,雨热同期,流域内降水比较丰沛,多年平均降水量达1 930.4 mm。流域降雨量年内分为明显的旱、雨季,其中雨季一般为每年的3-8月,约占全年降雨量的70%~80%,旱季为每年的11月至次年2月,其余时间为平水期。

图 1 研究区示意 Fig.1 Location map of the study site

漓江流域最发育的地层为泥盆系(D),出露面积为2 600 km2,约占流域范围的58%。其中上泥盆统融县组(D3r)灰岩是流域岩溶最发育的地层,分布较广,出露面积为643 km2,厚度为130~1 606 m,主要分布在桂林-葡萄-阳朔-兴坪-福利一带,同时区内还发育了少部分泥灰岩和硅质岩,主要分布在区内中南部的白沙和普益(图 2)。研究所选断面中上游大面圩断面(DM)位于桂林市区北部,其上游主要分布碎屑岩和碳酸盐岩夹层,而下游普益断面(PY)则主要流经较纯的碳酸盐岩区。

图 2 研究区岩性图(改自1:20万桂林幅水文地质普查报告) Fig.2 Lithologic map of the study area(Modified from the 1:200 000 Guilin geological survey report)
1.2 研究方法

于2016年12月到2017年11月在漓江上游大面圩断面和下游普益断面进行常规采样(图 1)。现场使用YSI ProDss进行水温(T)、pH值、电导率(Ec)、溶解氧(DO)、浊度(Tb)的监测,精度分别为0.1℃、0.01个pH单位、0.5%、0.1%和0.3 FNU。水样现场用直径50 mm,孔径为0.45 μm的醋酸纤维脂膜过滤。碱度在现场通过HCl滴定过滤水样进行测试,滴定剂为0.1 mmol/L HCl,滴定终点为pH值为4.5。用清洁的高密度聚乙烯瓶采集水样,采集之前先用水样润洗瓶子3遍,然后再装取水样。取600 mL过滤水样储存于高密度聚乙烯瓶中,用于阴离子和NH4+离子测定。取过滤水样装于事先用1:1体积浓度配置的HNO3溶液清洗过的50 mL聚乙烯瓶中,用于阳离子测定,取样后立即加入浓硝酸2~3滴,调节pH值< 2,用于阳离子测定。所有水样采集后立即放入便捷式冰袋保存,并于12 h内运至室内4℃冷藏。Cl-、NO3-、SO42-阴离子以及NH4+离子采用美国戴安公司ICS-900离子色谱仪分析,K+、Na+、Ca2+、Mg2+阳离子采用ICP-OES光谱仪分析(IRIS Intrepid Ⅱ XSP,Thermo Fisher Scientific,USA),测试精度为0.01 mg/L。阴阳离子分析误差 < 5%。以上测试工作均在中国地质科学院岩溶地质研究所岩溶地质资源环境监督检测中心完成。

2 结果与分析 2.1 漓江典型断面物理化学参数变化特征

DM断面旱季水温为11.9~26.8℃,平均值为16.6℃,雨季水温为12.1~25.3℃,平均值为20.4℃。PY断面旱季水温为13.1~29.9℃,平均值为18.4℃,雨季水温为13.6~28.3℃,平均值为23.2℃(表 1~2)。DM和PY断面水温都表现出雨季高于旱季(图 3),主要是由于研究区属于亚热带季风气候,雨热同期,且由于两个断面都是地表河流,对气温变化响应明显,变化较大。

表 1 DM断面物理化学参数 Table 1 Major physicochenmical indictors at DM transect

表 2 PY断面物理化学参数 Table 2 Major physicochenmical indictors at PY transect

图 3 DM和PY断面物理化学参数和离子月动态变化 Fig.3 Monthly dynamics of physicochenmical paremeters and ions at DM and PY transects

DM断面旱季pH值为7.58~8.04,平均值为7.87,雨季pH值为7.15~7.97,平均值为7.56。PY断面旱季pH值为7.77~8.85,平均值为8.33,雨季pH值为7.58~8.00,平均值为7.76(图 3)。DM和PY断面pH值都表现出雨季低于旱季,主要是由于研究区雨热同期,雨季土壤呼吸作用强烈,土壤中CO2浓度高,水-气交换活跃,雨水经土壤层淋滤溶解了大量CO2后输入河流的缘故。此外,研究区位于西南酸雨区重庆-贵阳-柳州一线,多年雨水pH值平均值约为4.5~5.0[12],因此,雨季pH值较低。

电导率(Ec)反映了河流中离子的强度,可以指示总离子组成及溶解态的无机物质组成。DM断面旱季Ec为88.1~113.3 μs/cm,平均值为100.9 μs/cm;雨季Ec为78.2~110.2 μs/cm,平均值为91.9 μs/cm。PY断面旱季Ec为185.6~229.9 μs/cm,平均值为201.7 μs/cm;雨季Ec为163.7~222.5 μs/cm,平均值为194.7 μs/cm。DM和PY断面都表现出雨季Ec小于旱季,且雨季波动范围更大,可能主要是雨水稀释作用和河流流量增大的影响,从图 3可以看出,两个断面雨季流量明显大于旱季。此外,DM断面和PY断面的Ec、Ca2+和HCO3-离子都表现出一致的季节变化,可能主要是受到岩溶作用的影响,而PY断面Ca2+和HCO3-离子都显著高于DM断面,说明PY断面流经大规模分布的碳酸盐岩地区后受到更强烈的岩溶作用的影响,这与两个断面地质背景一致。

DM断面旱季总溶解固体(TDS)为57.3~73.7 mg/L,平均值为65.6 mg/L,雨季TDS为50.7~72.0 mg/L,平均值为59.7 mg/L。PY断面旱季TDS为120.6~149.4 mg/L,平均值为130.0 mg/L,雨季TDS为106.4~145.0 mg/L,平均值为126.6 mg/L。DM和PY断面TDS都表现出旱季大于雨季,可能是由于旱季降雨量较少,水岩相互作用时间更长,使得旱季较雨季经受更强烈的化学侵蚀作用。从两个断面监测期间TDS的平均值来看,PY断面(128.8 mg/L)远高于DM断面(62.7 mg/L),而DM断面则略低于世界河流TDS的平均值,表明PY断面比DM断面更强烈的化学侵蚀作用[13-14]

2.2 漓江典型断面主要离子组成变化特征

河流阴、阳离子三角图反映了河水的主要组成变化。DM和PY断面均以Ca2+、HCO3-离子为主,在阳离子中,DM和PY断面Ca2+离子所占比重最高,分别为80.5%~86.0%和75.1%~95.3%,平均值分别为82.8%和82.7%。Na+、K+和Mg2+一般都在20%以下。在阴离子中,DM和PY断面HCO3-离子所占比重最高,分别为74.3%~95.8%和72.0%~96.2%,平均值分别为83.9%和85.2%。Cl-、NO3-和SO42-一般都在20%以下(图 4)。说明漓江典型断面水化学类型为Ca-HCO3型,而这两种离子的浓度与岩溶作用强度有着密切联系,显示了碳酸盐岩溶蚀-沉淀过程对漓江典型断面水化学特征的控制作用。

图 4 DM和PY断面阴阳离子三角图 Fig.4 Triangular diagrams of the anions and cations at DM and PY transects

K+、Na+、Ca2+、Mg2+、Cl-、NO3-、SO42-、HCO3-是地下水中分布最广泛、含量最高的8种离子,在很大程度上决定了地下水水化学特征。通常岩溶水中Ca2+、Mg2+、HCO3-离子主要来源于岩溶作用,形成岩溶水体的优势离子[15],由图 3可知DM和PY断面的Ca2+、Mg2+、HCO3-离子都表现出雨季浓度偏小但波动较大,尤其在雨季出现极低值,可能主要受到大气降水的稀释作用和河流流量增大的影响,旱季浓度偏大但相对平稳,可能是由于旱季河流流量较小且稳定,同时长时间水岩作用的地下水补给河流,使得旱季离子浓度偏高。DM和PY断面的Ca2+、HCO3-离子之间存在较好的正相关关系(表 3~4),显示出两者的同源性即来源于岩溶作用。

表 3 DM断面阴阳离子相关关系矩阵 Table 3 Correlation matrix for anions and cations at DM transect

表 4 PY断面阴阳离子相关关系矩阵 Table 4 Correlation matrix for anions and cations at PY transect

Cl是保守元素,在地表水循环中也相对稳定,一般认为Cl主要来自海洋粒子,但也会受到人类活动的影响。海水中Cl-/Na+值为1.17,如果受到大气降水影响则Cl-/Na+值应该接近此值,而如果受到其他来源如人类活动影响则应该大于此值[16]。DM断面Cl-与Na+离子并没有表现出相关性,且Cl-/Na+平均值为0.93(表 3),要远小于大气降水和人为活动来源的比值,说明DM断面受到其它来源的Na+离子的干扰,DM断面上游为大溶江流域,广泛分布硅酸盐岩,其风化过程中产生Na+离子,而不产生Cl-,因而,DM断面受到上游硅酸盐岩风化的影响,应该是此部分Na+离子的主要来源[11],从而影响了大气降水和人为活动来源的比值。此外,DM断面Cl-、NO3-、SO42-相互之间都表现出正相关关系,表明可能具有较为一致的来源,可能主要受到人为活动的影响。PY断面虽然Cl-与Na+离子表现出一定的相关性,但是Cl-/Na+平均值高达1.62,说明有其它来源贡献了过量的Cl-离子。PY断面K+与Na+离子、Na+与Cl-、NO3-离子以及Cl-、NO3-、SO42-相互之间都表现出正相关关系,表明可能具有较为一致的来源,主要可能受到人为活动的影响。

2.3 漓江典型断面离子的来源特征

影响水体中溶质组成变化的主要因素是物质来源,但其物质来源常受自然条件(如岩性、土壤、大气沉降等)和人类活动(如工业、农业、城市等)的双重影响。Gibbs模型图是一种半对数坐标图,其纵坐标为对数坐标,代表水体的TDS,横坐标为普通坐标,表示水体的部分阴离子的比值c(Cl-)/c(Cl-+HCO3-)和部分阳离子的比值c(Na+)/c(Na++Ca2+),可以简单明了地显示自然水体中各种离子的起源机制(大气降水、水岩作用及蒸发浓缩效应)及其变化趋势过程。在Gibbs图中,中间部位表示水体水化学特征主要受到水岩相互作用控制;右上角表示主要受到蒸发结晶作用控制;右下角表示主要受到降雨作用控制[17]。从图 5中可以看出PY断面旱季和雨季河流水文地球化学特征变化不大,且都受到水岩相互作用的控制,大气降水的溶质输入较弱。DM断面旱雨季河流水文地球化学特征几乎没有变化,也都受到水岩相互作用的控制,但是相对于PY断面,DM断面明显也受到部分降雨作用的影响。

图 5 DM和PY断面Gibbs图 Fig.5 Gibbs plot of river water at DM and PY transects

Gaillardet等[1]根据单一岩性小流域的水化学特征,描述了3大岩类端元,并对全球60条大河进行了分析。本研究基于此作出DM和PY断面Ca2+/Na+与HCO3-/Na+、Ca2+/Na+与Mg2+/Na+关系图(图 6)。两个断面Ca2+/Na+、Mg2+/Na+、HCO3-/Na+的值都比较高,基本都位于图的右上角,说明两个断面河水离子主要来源于碳酸盐岩风化,而受硅酸盐岩影响较小。

图 6 DM和PY断面Ca2+/Na+与HCO3-/Na+、Ca2+/Na+与Mg2+/Na+关系 Fig.6 Relationship of Ca2+/Na+ and HCO3-/Na+、Ca2+/Na+ and Mg2+/Na+ for river water at DM and PY transects

根据岩溶作用方程:

$ {\rm{CaC}}{{\rm{O}}_3} + {\rm{C}}{{\rm{O}}_2} + {{\rm{H}}_2}{\rm{O = C}}{{\rm{a}}^{{\rm{2 + }}}} + {\rm{HCO}}_3^ - , $ (1)
$ \begin{array}{l} {{\rm{C}}{\rm{a}}}{\rm{Mg}}{\left( {{\rm{C}}{{\rm{O}}_3}} \right)_2} + 2{\rm{C}}{{\rm{O}}_2} + 2{{\rm{H}}_2}{\rm{O = }}{{\rm{C}}{\rm{a}}}^{2 + } + \\ {\rm{M}}{{\rm{g}}^{{\rm{2 + }}}} + 4{\rm{HCO}}_3^ - 。\end{array} $ (2)

当纯方解石溶解平衡时,水体中的Mg2+/Ca2+为0;当纯白云岩溶解平衡时,水体中的Mg2+/Ca2+为1;当水体中方解石和白云岩溶解平衡时,Mg2+/Ca2+为0.5,水体中Mg2+/Ca2+值的变化主要受到地层中方解石和白云岩的比例控制[16],因此根据Mg2+/Ca2+值的变化可以来判断方解石和白云岩对地下河水化学的贡献。从图 7可以看出,DM和PY断面两个季节溶质组分都主要来源于方解石的溶解,表明方解石溶解-沉淀过程对两个断面河水水化学特征的控制作用。

图 7 DM和PY断面Mg2+/Ca2+ 和HCO3-关系 Fig.7 Relationship between Mg2+/Ca2+ and HCO3- at DM and PY transects

此外,根据上述岩溶作用方程,(Ca2++Mg2+)与HCO3-的摩尔比为0.5,而从下列公式:

$ \begin{array}{l} {\rm{C}}{{\rm{a}}_{{\rm{1 - }}x}}{\rm{M}}{{\rm{g}}_x}{\rm{C}}{{\rm{O}}_3} + {\rm{HN}}{{\rm{O}}_{\rm{3}}} = \left( {1 - x} \right){\rm{C}}{{\rm{a}}^{2 + }} + \\ x{\rm{M}}{{\rm{g}}^{{\rm{2 + }}}} + {\rm{HCO}}_3^ - + {\rm{NO}}_3^ - , \end{array} $ (3)
$ \begin{array}{l} 2\left( {{\rm{C}}{{\rm{a}}_{1 - x}}{\rm{M}}{{\rm{g}}_x}{\rm{C}}{{\rm{O}}_3}} \right) + {{\rm{H}}_2}{\rm{S}}{{\rm{O}}_4}{\rm{ = }}2\left( {1 - x} \right){\rm{C}}{{\rm{a}}^{2 + }} + \\ 2x{\rm{M}}{{\rm{g}}^{{\rm{2 + }}}} + {\rm{HCO}}_3^ - + {\rm{SO}}_4^{2 - }。\end{array} $ (4)

可知,硫酸和硝酸参与对碳酸盐溶解的影响中,(Ca2++Mg2+)与HCO3-的摩尔比都为1,DM和PY断面两个季节的(Ca2++Mg2+)与HCO3-的摩尔比在监测期间绝大部分落在0.5~1,说明有外源酸参与了碳酸盐岩的溶解,从而带来了额外的Ca2+、Mg2+离子,在两个断面(Ca2++Mg2+)与HCO3-的摩尔比值上加上SO42-和NO3-后,两个断面(Ca2++Mg2+)与(HCO3-+SO42-)的摩尔比更接近0.5(图 8),表明硫酸参与了碳酸盐岩的风化,带来了额外的Ca2+、Mg2+离子。DM断面旱季(Ca2++Mg2+)与HCO3-的摩尔比在0.55~0.61,平均值为0.58;雨季在0.52~0.68,平均值为0.59。PY断面旱季在0.44~0.66,平均值为0.56,雨季在0.35~0.62,平均值为0.52。以上结果表明,DM断面硫酸和硝酸溶蚀对河水HCO3-的贡献率在旱季为16%,雨季为19%,PY断面该贡献率旱季为12%,雨季为5%。外源酸参与碳酸盐岩风化减少岩溶风化碳汇已有较多报道[18-20],毫无疑问,DM和PY断面由于外源酸的输入将会导致岩溶风化碳汇通量减少。

图 8 DM和PY断面(Ca2++Mg2+)和HCO3-、(Ca2++Mg2+)和(HCO3-+SO42-+NO3-)关系 Fig.8 Relationship between (Ca2++Mg2+) and HCO3-, (Ca2++Mg2+) and (HCO3-+SO42-+NO3-) at DM and PY transects

在所有的水样中,Cl-、NO3-、SO42-等阴离子的检出率为100%。河流中的NO3-主要来源于大气降水、土壤层中有机氮的转化和人类活动(包括施肥、排污等)[21-22]。从空间变化上,下游PY断面NO3-离子浓度明显大于上游DM断面,两个断面都是人类活动密集区域,旅游排污、农业活动密集,可能受到人为活动的影响。而两个断面NO3-离子浓度都表现出雨季大于旱季,可能是由于雨季是研究区施肥旺季,同时雨季气温较高,土壤和地下水中的微生物活动活跃,硝化作用进行快速从而加速了土壤氮的转化和运移,最后随着降雨进入河流中,造成河流中雨季NO3-较高[23-24]。河水中的SO42-多来自于硫化矿物的溶解、工矿业活动和大气沉降。在岩溶地下河系统中,Ca2+的来源主要是碳酸盐岩的溶解,较为单一,即使存在石膏等矿物的溶解,其阴离子也会相应的出现同步变化,因此可以用SO42-/Ca2+和NO3-/Ca2+的相互变化关系来指示河流溶质的来源特征[25]。DM和PY断面两个季节趋向两个方向的变化(图 9):一个是NO3-/Ca2+的比值低而SO42-/Ca2+的比值高,可能是受到工矿业活动对河流的影响;另一个是NO3-/Ca2+的比值高而SO42-/Ca2+的比值低,可能是受到农业活动或城市污水排放的影响。两个方向的对比发现,DM和PY断面在两个季节大部分时间的SO42-/Ca2+比值都分布在比值低的区域,只有PY断面(主要是雨季)部分分布在了SO42-/Ca2+比值高的区域,说明漓江典型断面主要受到农业活动和城市污水排放的影响,下游PY断面则在雨季受到部分工矿业的影响,这也与野外调查情况吻合。漓江流域主要分布在岩溶洼地,谷地以及峰林峰丛平原上,是主要的农业分布区和城镇聚集区,并且旅游业发达,因此一切的土地利用方式的改变以及污水排放都会引起河流溶质组分的改变,并且极易造成污染。

图 9 DM和PY断面NO3-/Ca2+和SO42-/Ca2+关系 Fig.9 Relationship between NO3-/Ca2+ and SO42-/Ca2+ at DM and PY transects

综上所述,虽然DM和PY断面水化学特征在监测期间受到水岩相互作用、人类活动以及降雨-流量稀释作用的影响,但HCO3-、Ca2+仍分别在两个断面旱、雨季阴、阳离子中占绝对优势,说明在监测期间河流主要还是受到水岩相互作用即方解石溶解-沉淀的影响。影响方解石溶解-沉淀最关键的因素就是水中CO2浓度的变化,且通过影响河水中pH值和pCO2来控制河流水化学的变化。从图 10中可以看出两个断面pH值和pCO2显示出很好的负相关性,说明两者的变化受到一个共同的因素即水体中CO2浓度的控制,进而驱动岩溶作用的进行。

图 10 DM和PY断面pH与pCO2关系 Fig.10 Relationship between pH and pCO2 at DM and PY transects
3 结论

漓江河水主要组分受到水岩作用的控制,上游DM断面主要流经碎屑岩和碳酸盐岩夹层,下游PY断面主要流经裸露的碳酸盐岩,但是两个断面主要化学类型都为Ca-HCO3型,这两种离子的浓度与岩溶作用强度有着密切联系,显示了碳酸盐岩溶蚀-沉淀过程对漓江典型断面水化学特征的控制作用,上游DM断面还受到部分降雨输入的影响,PY断面则受到更强烈的岩溶作用的影响。

虽然碳酸盐岩风化(主要石灰岩)控制了漓江河水的主要离子组分,硅酸盐风化对河水离子组分影响较小,但是从(Ca2++Mg2+)与HCO3-的摩尔比的变化可知,漓江碳酸盐岩风化受到外源酸的影响,进而可能影响岩溶风化碳汇。漓江河水由于受到雨季降雨和流量增大的影响,使得两个断面主要离子HCO3-、Ca2+的离子浓度都出现雨季小于旱季的情况。另外从SO42-/Ca2+和NO3-/Ca2+比值图可知,漓江河水主要受到农业活动、城市排污的影响,而PY断面受到部分工矿业的影响,尤其是在雨季,且由于降雨的影响,使得漓江SO42-、NO3-、Cl-离子出现雨季大于旱季的情况。

虽然漓江河水水化学特征受到水岩相互作用、人类活动以及降雨-流量稀释作用的共同影响,但HCO3-、Ca2+还是河流的优势离子,且受到水中CO2浓度变化的影响,其通过影响河水中pH值和pCO2,进而控制了河流水化学的变化。

参考文献
[1]
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/2/3/4): 3-30.
[2]
HAN G, 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.
[3]
WU W. Hydrochemistry of inland rivers in the north Tibetan Plateau:Constraints and weathering rate estimation[J]. Science of the Total Environment, 2016, 541: 468-482. DOI:10.1016/j.scitotenv.2015.09.056
[4]
蒲俊兵.重庆市地下河发育、分布的控制机制及水文地球化学区域特征研究[D].重庆: 西南大学, 2011.
PU J B.Research on the controlling factors of formation and distribution of subterranean karst streams and its hydrogeochemistry regionality, Chongqing, China[D].Chongqing: Southwest University, 2011.
[5]
乐嘉祥, 王德春. 中国河流水化学特征[J]. 地理学报, 1963(1): 1-13.
LE J X, WANG D C. The river water chemistry characteristics in China[J]. Acta Geographica Sinica, 1963(1): 1-13. DOI:10.3321/j.issn:0375-5444.1963.01.001
[6]
陈静生, 王飞越, 何大伟. 黄河水质地球化学[J]. 地学前缘, 2006, 13(1): 58-73.
CHEN J S, WANG F Y, HE D W. Geochemistry of water quality of the Yellow River basin[J]. Earth Science Frontiers, 2006, 13(1): 58-73. DOI:10.3321/j.issn:1005-2321.2006.01.009
[7]
陈静生, 王飞越, 夏星辉. 长江水质地球化学[J]. 地学前缘, 2006, 13(1): 74-85.
CHEN J S, WANG F Y, XIA X H. Geochemistry of water quality of the Yangtze River basin[J]. Earth Science Frontiers, 2006, 13(1): 74-85. DOI:10.3321/j.issn:1005-2321.2006.01.010
[8]
王亚平, 王岚, 许春雪, 等. 长江水系水文地球化学特征及主要离子的化学成因[J]. 地质通报, 2010, 29(2/3): 446-456.
WANG Y P, WANG L, XU C X, et al. Hydro-geochemistry and genesis of major ions in the Yangtze River, China[J]. Geological Bulletin of China, 2010, 29(2/3): 446-456.
[9]
韩贵琳, 刘丛强. 贵州喀斯特地区河流的研究——碳酸盐岩溶解控制的水文地球化学特征[J]. 地球科学进展, 2005(4): 394-406.
HAN G L, LIU C Q. Hydrogeochemistry of rivers in Guizhou Province, China:Constraintson crustal weathering in karst terrain[J]. Advances in Earth Science, 2005(4): 394-406. DOI:10.3321/j.issn:1001-8166.2005.04.004
[10]
于奭, 孙平安, 杜文越, 等. 人类活动影响下水化学特征的影响:以西江中上游流域为例[J]. 环境科学, 2015, 36(1): 72-79.
YU S, SUN P A, DU W Y, et al. Effect of hydrochemistry characteristics under impact of human activity:A case study in the upper reaches of the Xijiang River basin[J]. Environmental Science, 2015, 36(1): 72-79.
[11]
孙平安, 于奭, 莫付珍, 等. 不同地质背景下河流水化学特征及影响因素研究:以广西大溶江、灵渠流域为例[J]. 环境科学, 2016, 36(1): 123-131.
SUN P A, YU S, MO F Z, et al. , Hydrochemical characteristics and influencing factors in different geological background:A case study in Darongjiang and Lingqu basin, Guangxi, China[J]. Environmental Science, 2016, 36(1): 123-131. DOI:10.3969/j.issn.1673-1212.2016.01.035
[12]
汤洁, 徐晓斌, 巴金, 等. 1992—2006年中国降水酸度的变化趋势[J]. 科学通报, 2010, 55(8): 705-712.
TANG J, XU X B, BA J, et al. , The change trend of precipitation acidity in China during 1992—2006[J]. Chinese Science Bulletin, 2010, 55(8): 705-712.
[13]
MEYBECK M. Global chemical weathering of surficial rocks estimated from river dissolved loads[J]. American Journal of Science, 1987, 287(5): 401-428. DOI:10.2475/ajs.287.5.401
[14]
GAO Q, TAO Z, HUANG X, et al. Chemical weathering and CO2 consumption in the Xijiang River basin, South China[J]. Geomorphology, 2009, 106(3/4): 324-332.
[15]
WHITE W B. Geomorphology and hydrology of karst terrains[M]. England: Oxgord University Press, 1988: 464.
[16]
MA R, WANG Y, SUN Z, et al. Geochemical evolution of groundwater in carbonate aquifers in Taiyuan, northern China[J]. Applied Geochemistry, 2011, 26(5): 884-897. DOI:10.1016/j.apgeochem.2011.02.008
[17]
GIBBS R J. Mechanisms controlling world water chemistry[J]. Science, 1970, 170(3962): 1088-1090. DOI:10.1126/science.170.3962.1088
[18]
PERRIN A S, PROBST A, PROBST J L. Impact of nitrogenous fertilizers on carbonate dissolution in small agricultural catchments:Implications for weathering CO2 uptake at regional and global scales[J]. Geochimica et Cosmochimica Acta, 2008, 72(13): 3105-3123. DOI:10.1016/j.gca.2008.04.011
[19]
JIANG Y, HU Y, SCHIRMER M. Biogeochemical controls on daily cycling of hydrochemistry and δ13C of dissolved inorganic carbon in a karst spring-fed pool[J]. Journal of Hydrology, 2013, 478: 157-168. DOI:10.1016/j.jhydrol.2012.12.001
[20]
ZHANG T, PU J, LI J, et al. Stable isotope and aquatic geochemistry of a typical subtropical karst subterranean stream in southwest China[J]. Carbonates and Evaporites, 2017, 32(3): 415-430. DOI:10.1007/s13146-017-0356-3
[21]
EDMUNDS W, SHAND P. Natural groundwater quality[M]. John Wiley & Sons.UK: Blackwell Publishing, 2008: 2-21.
[22]
沈照理, 朱宛华, 钟佐案. 水文地球化学基础[M]. 北京: 地质出版社, 1999: 189.
SHEN Z L, ZHU W H, ZHONG Z A. Hydrogeochemical basis[M]. Beijing: Geological Publishing House, 1999: 189.
[23]
VITÒRIA L, SOLER A, CANALS À, et al. Environmental isotopes (N, S, C, O, D) to determine natural attenuation processes in nitrate contaminated waters:Example of Osona (NE Spain)[J]. Applied Geochemistry, 2008, 23(12): 3597-3611. DOI:10.1016/j.apgeochem.2008.07.018
[24]
申宏岗, 曹建华, 潘根兴. 岩溶生态系统中土壤微生物量氮含量及其变化规律[J]. 生态环境, 2007, 16(6): 1728-1732.
SHEN H G, CAO J H, PAN G X. Contents of soil microbial biomass nitrogen and its change in karst ecosystem[J]. Ecology and Environment, 2007, 16(6): 1728-1732. DOI:10.3969/j.issn.1674-5906.2007.06.029
[25]
蒲俊兵, 袁道先, 蒋勇军, 等. 重庆岩溶地下河水文地球化学特征及环境意义[J]. 水科学进展, 2010, 21(5): 628-636.
PU J B, YUAN D X, JIANG Y J, et al. Hydrogeochemistry and environmental meaning of Chongqing subterranean karst streams in China[J]. Advances in Water Science, 2010, 21(5): 628-636.