喀斯特流域水土流失对水库碳汇效应的影响

王 翠, 周忠发, 李永柳, 张勇荣

(1.贵州师范大学 地理与环境科学学院/喀斯特研究院, 贵州 贵阳 550001;2.贵州省喀斯特山地生态环境国家重点实验室培育基地, 贵州 贵阳 550001)

水土流失影响生态安全、饮水安全、防洪安全、粮食安全，威胁并制约经济社会发展。近年来，随着全球碳循环研究的深入，水土流失导致的土壤碳迁移已被广泛关注。已开展降雨过程对土壤有机碳流失规律、有机碳在流失泥沙中的富集规律、泥沙含碳量与泥沙量的关系、泥沙有机碳富集比与雨强的关系等方面的研究[1-3]。还有水土保持措施减少土壤有机碳的输出[4]、水土流失影响土壤有机碳的分解转化、二氧化碳的释放和土壤碳汇、陆地生态系统碳循环的科学评估等[5]方面的研究。这些研究都是针对水土流失对土壤碳库的影响，对沉积区碳库的影响研究却相对较少。水库作为水土流失的沉积区，沉积物可揭示流域土壤侵蚀状态[6]，是土壤碳库与水库碳库产生交汇作用的重要区域，每年大概有1.50×1010～2.00×1010mg土壤有机碳及氮、磷、钾、硫等元素通过水土流失途径被输送到河流湖库中[7-8]。尽管我们知道水土流失带来的颗粒性有机碳使水库成为流失区碳的净汇并向大气释放温室气体[9-10]，但这个结论仍是对土壤碳库而言，水土流失对水库碳汇来说是碳源还是碳汇目前判断不一[11-14]。水库是陆地生态系统的重要组成部分，是全球“碳遗失汇”的重点关注区，前人因为对喀斯特作用的认识不足，无法科学计算和评估喀斯特地区水库碳汇效应，因此开展喀斯特区域水土流失对水库碳汇效应的影响研究，是对喀斯特地区水库碳汇效应研究的补充。平寨水库流域位于贵州省毕节市，喀斯特面积占区域总面积的73.34%。平寨水库是流域内唯一水库，是长江上游一级支流乌江上游的大型水库，也是贵州省首个大型跨地区、跨流域长距离水利调水工程的核心水源工程，其生态安全对贵州省实施西部大开发战略具有重要意义。目前在平寨水库流域开展的研究多为水质、水化学、生物碳泵等方面的研究[15-16]，尚未见流域水土流失对水库碳汇效应的研究报道。因此，本研究以平寨水库流域为研究区，开展流域水土流失对水库碳汇效应的影响研究，为区域水土流失治理，科学评估喀斯特地区河流生态系统碳循环提供科学依据。

1 材料与方法

1.1 研究区概况

研究区为平寨水库流域，是黔中水利枢纽工程的核心水源工程集雨区，流域面积833.77 km2。区域内碳酸盐岩与非碳酸盐岩交替出现，呈典型多层状碳酸盐岩组结构。平寨水库是研究区唯一水库，属典型喀斯特深水水库。水库主要干流为纳雍河，还有水公河、张维河、白水河、扈家河4条支流。纳雍河流域分布有煤炭资源和硫铁矿矿床[17-18]，其余区域未见大型矿产资源分布。研究区地形包括山地和丘陵两种，土地利用类型有耕地、林地、草地、建设用地、水域和未利用地6类，地表覆被类型以耕地和林地为主(图1)。水库2015年4月下闸蓄水，2018年12月通过蓄水验收并于月底开始向贵阳供水。

注：ZW，BS，HJ，PZ1，PZ2，SG1，SG2和NY分别为8个监测点的编号。

1.2 样品采集与分析

1.3 二氧化碳分压计算

二氧化碳分压根据亨利定律[20]进行计算，计算公式为：

(1)

(2)

(3)

(4)

式中：Ki为平衡常数,计算方程为：

pKCO2=-7×10-5T2+0.016T+1.11

(5)

pK1=1.1×10-4T2-0.012T+6.58

(6)

pK2=9×10-5T2-0.013 7T+10.62

(7)

式中：T为温度(℃)。

根据亨利定律，pCO2可根据以下公式计算：

(8)

α(H+)=10-[pH]

(9)

(10)

(11)

式中：α表示离子活度；I表示离子强度。

2 结果与分析

2.1 研究区水土流失分析

河流泥沙主要来源于流域水土流失，因此以干流纳雍河阳长水文站2014—2018年水文数据为代表进行流域水土流失分析(图2)。2014—2018年，纳雍河年流量最大1.66×109m3(2014年)，最小11.51 m3(2016年)。平寨水库流域水资源为降雨补给型，因此河流年流量呈自然起伏状态，降雨量多的年份流量大。泥沙通量在2.54×105～1.23×106t/a之间，最大值在2014年，最小值在2018年。年流量和年泥沙通量变化不同步，主要是由于2016年泥沙通量的突然升高。将平寨水库2014—2018年年际间的流量和泥沙通量进行相关分析，R2为0.486 4，无显著相关性。泥沙通量在2016年突然升高是由于下闸蓄水后水面抬升，淹没耕地和林地，淹没区的沉积物或土壤被水流不断侵蚀，短时间内形成大量泥沙进入河流。随着水库的调度运转，流速降低，河流两岸可侵蚀面积稳定，水库泥沙通量逐年降低。

图2 纳雍河2014-2018年水流量与泥沙通量变化

按月时间尺度对纳雍河流量和含沙量作对数、线性、指数、乘幂和多项式回归分析，含沙量与降雨的变化趋势分析，含沙量与降雨等的相关分析，分析结果见图3和表1。流量与泥沙含量回归方程R2最大为0.760 1，为指数回归模型，该结果与长江宜昌站分析结果不同[21]，表明研究区不同，水土流失机制和河流水文情势有差异，影响含沙量与流量关系。平寨水库流域降雨和含沙量表现为春夏季逐渐增大，秋冬季逐渐降低的变化周期，峰值出现在夏季。含沙量随降雨的变化而变化，峰值一般与降雨峰值同时出现。相关分析表明，研究区月时间尺度上含沙量与降雨、蒸发量、流量、流速都呈极显著正相关，流量和流速对含沙量的影响更大。因此在平寨水库流域，降雨是造成流域内水土流失的主要原因之一，短时强降雨对水库泥沙通量影响显著，春夏季水土流失剧烈，秋冬季减弱。

表1 研究区含沙量与降雨等参数的相关分析

图3 纳雍河月流量与含沙量回归分析与变化趋势分析

2.2 水库水化学分析

表2 水化学基本参数及水化学类型

图4 平寨水库水化学分析结果

2.3 水库碳汇效应分析

根据碳通量计算公式：

Flux=K(Cwater-Cair)

式中：Flux为CO2交换通量〔mg/(m2·h)〕；Cwater为CO2在水体中的浓度(mmol/L)；Cair为空气中CO2浓度(mmol/L)；K为CO2交换系数(cm/h)。

当K>0时,水库二氧化碳分压大于大气二氧化碳分压,水库向大气排放二氧化碳,为碳源,反之为碳汇。

表3 pH值、叶绿素a及与二氧化碳分压相关分析

图5 平寨水库二氧化碳分压

图6 平寨水库碳汇效应季节变化

3 讨 论

3.1 平寨水库流域水土流失变化及对水环境的影响

研究认为流量与泥沙含量一般规律为流量增加，水动力增强，水流挟沙力增大，泥沙含量升高[27]。平寨水库筑坝后水位上升，水动力减弱，河流泥沙通量逐年降低，与长江主要干支流水库建设对河流输沙量的影响研究结果一致[28]。

平寨水库流域水土流失与降雨、蒸发量、流速、流量在月时间尺度上呈极显著正相关，在年时间尺度上无相关性。一是用于统计分析的含沙量是水土流失的一部分，从年时间尺度上看，样本含量较少。二是由于水库筑坝破坏河流连续和洪水脉动规律，水流对泥沙的裹挟能力减弱，在此作用下，水土流失与降雨等影响因素失去年时间尺度上的相关性。三是影响水土流失的因素众多，各作用因子对水土流失的影响差异大[29-30]，在各因子的交互作用下，平寨水库流域降雨、流量、流速与水土流失在年际间较难形成相关关系，但在短时间尺度上，降雨是造成平寨水库流域水土流失的重要原因。

3.2 平寨水库流域水土流失对水库碳汇效应的影响

4 结 论

(1) 2014—2018年平寨水库流域干流纳雍河流量1.15×109～1.66×109m3/a，泥沙通量2.54×105～1.23×106t/a，泥沙通量和流量在年际间无相关性，在月时间尺度上极显著正相关，含沙量与流量是指数回归关系，与降雨峰值一般同时出现，短时强降雨是流域发生水土流失的主要原因之一。流域内春夏季水土流失剧烈，秋冬季减弱。

(3) 平寨水库二氧化碳分压具有时空异质性，表现为：冬季>秋季>春季>夏季，秋冬季为碳源，春夏季为碳汇。平寨水库流域水土流失与水库碳汇效应且具有相同的季节变化周期，水土流失会减弱水库的碳汇效应。