生态学报  2017, Vol. 37 Issue (22): 7636-7646

文章信息

李凌宇, 于瑞宏, 田明扬, 胡海珠, 张笑欣, 冉立山, 杨现坤, 吕喜玺.
LI Lingyu, YU Ruihong, TIAN Mingyang, HU Haizhu, ZHANG Xiaoxin, RAN Lishan, YANG Xiankun, LÜ Xixi.
黄河二氧化碳逸出时空变化及其影响因素——以头道拐水文站为例
Spatial-temporal variations and influencing factors of carbon dioxide evasion from the Yellow River:An example of the Toudaoguai Gauging Station
生态学报. 2017, 37(22): 7636-7646
Acta Ecologica Sinica. 2017, 37(22): 7636-7646
http://dx.doi.org/10.5846/stxb201609091833

文章历史

收稿日期: 2016-09-09
网络出版日期: 2017-07-12
黄河二氧化碳逸出时空变化及其影响因素——以头道拐水文站为例
李凌宇1 , 于瑞宏1 , 田明扬1 , 胡海珠1 , 张笑欣1 , 冉立山2 , 杨现坤2 , 吕喜玺1     
1. 内蒙古大学生态与环境学院, 呼和浩特 010021;
2. 新加坡国立大学, 新加坡 17570
摘要: 河流连接着海洋碳库和陆地碳库,河流碳逸出是全球碳收支的重要组成部分。本文以黄河上游和中游分界点-内蒙古段头道拐断面为研究对象,采用Li-7000静态箱法监测了断面4个采样点在2013-2015年期间四季的CO2逸出通量(FCO2),并分析了FCO2时空变化规律。基于研究断面主要水文和水化学指标的野外监测和室内分析结果,探讨了FCO2的主要影响因素。结果表明断面FCO2介于14-186 mol m-2 a-1,平均值为84 mol m-2 a-1;水体CO2分压(pCO2)介于467-2101 μatm,平均值为995 μatm;DOC浓度介于2.7-13 mg/L。FCO2季节性差异明显:夏季FCO2为全年最大456 mmol m-2 d-1,冬季最小33 mmol m-2 d-1FCO2在4个采样点的空间差异显著:河道右岸S4点处最大为392 mmol m-2 d-1;河道中部S2和S3点基本相同;河道左岸S1点最小为86 mmol m-2 d-1FCO2与河道流速呈现较好的正相关关系,与pCO2中等相关,与pH负相关,与风速的相关性不明显,说明对于该研究断面河道流速较pCO2FCO2的贡献更大。本研究较为精细地探讨了头道拐断面的水体CO2逸出规律,表明即使在同一河道断面,FCO2也可能存在较大空间差异,流速较大处的FCO2较大,因此在野外监测FCO2时需要在河道断面选取具有代表性的采样点,特别是较大的河流。研究结果为黄河中上游CO2逸出量评价和河道断面FCO2监测点的布设提供了科学依据。
关键词: 二氧化碳逸出通量     时空变化     影响因素     黄河头道拐    
Spatial-temporal variations and influencing factors of carbon dioxide evasion from the Yellow River:An example of the Toudaoguai Gauging Station
LI Lingyu 1, YU Ruihong 1, TIAN Mingyang 1, HU Haizhu 1, ZHANG Xiaoxin 1, RAN Lishan 2, YANG Xiankun 2, LÜ Xixi 1     
1. School of Ecology & Environment, Inner Mongolia University, Hohhot 010021, China;
2. Department of Geography, National University of Singapore, 17570 Singapore
Abstract: The oceanic carbon pool and terrestrial carbon pool are connected by rivers. Carbon dioxide (CO2) evasion from rivers to the atmosphere represents a substantial flux in the global carbon cycle. The CO2 efflux (FCO2) and CO2 partial pressure (pCO2) in large rivers have been widely evaluated. Most studies concerning CO2 emission from the Yellow River, a typical river containing high sediment concentrations, focused on the lower reach and its estuary, but less is known about its upper and middle reaches. In this study, a river cross-section at the Toudaoguai Gauging Station in Inner Mongolia, the dividing point between the upper and middle reaches of the Yellow River was chosen as a study site. Evasion of CO2 was measured four times each year using Li-7000 static chamber method from 2013 to 2015 at four sampling points in a river cross-section. The spatial and temporal variations of FCO2 were analyzed. The relevant hydrological indexes, including water temperature, pH, and wind velocity as well as current velocity were measured at the four sampling points. The hydrochemical indicators, including ALK and DOC in water samples, were analyzed in the laboratory and pCO2 was estimated. The possible influential factors of FCO2 were further discussed using correlation analysis. The CO2 evasion from the river cross-section ranged from 14 to 186 mol m-2 a-1 and its average was 84 mol m-2 a-1. The pCO2 in the Yellow River at the Toudaoguai Gauging Station was within the range of 467-2101 μatm and the average value was 995 μatm. The concentration of DOC ranged from 2 to 13 mg/L. The FCO2 exhibited obvious seasonal variations, with the maximum FCO2 of 456 mmol m-2 d-1 occurring in summer and the minimum of 33 mmol m-2 d-1 occurring in winter. The FCO2 values were markedly different at sampling points, with the maximum value of 392 mmol m-2 d-1 at S4 near the right bank, similar values at S2 and S3 in the middle of the river section, and the minimum of 86 mmol m-2 d-1 at S1 near the left bank. The analysis of factors influencing FCO2 indicated that FCO2 was positively correlated with current velocity and pCO2, and negatively correlated with pH. There was no obvious correlation between FCO2 and wind speed. Results also showed that current velocity contributed more to FCO2 than to pCO2 in the river cross-section. In this study, evasion of CO2 from the Toudaoguai cross-section was determined on a relatively fine scale. The results suggested that a distinct spatial variation in FCO2 exists even at the level of river cross-section, with the maximum FCO2 found at the point with the highest current velocity. Thus, typical sampling points in a river cross-section should be chosen for FCO2 measurement. The study provided a scientific reference for both FCO2 evaluation in the upper and middle reaches of the Yellow River and FCO2 sampling in a river cross-section.
Key words: evasion of carbon dioxide     spatial-temporal variation     influencing factor     Toudaoguai Gauging Station of the Yellow River    

自然界中不同形态的碳在岩石圈、陆地生态系统、大气和海洋4个碳库间的相互转换和运移过程称为全球碳循环[1]。河流作为全球碳循环的重要连接通道, 持续地从陆地向海洋输送碳, 在此过程中河流由水体向大气释放CO2, 是碳损失的重要途径之一[2-3]。已有研究表明黄河CO2逸出量与碳沉降量相当, 约占河流碳循环的1/3, 均大于向海洋的输送量[4]。全球河流系统的CO2年逸出量高达0.9—1.8亿t(GtC)[5], 与矿石燃料燃烧的碳排放通量和海洋-大气间的净碳通量属于同一量级[6-7]。因此, 开展河流水-气界面CO2通量研究, 对认识河流碳循环、区域碳循环乃至全球碳循环过程具有重要意义。

河流水体中的CO2在处于过饱和状态, 即水体中的CO2分压(pCO2)大于大气压时, 才能从水体逸出。因此, FCO2的研究通常以河流水体表面pCO2为基础。目前, 研究者在世界范围内的大江大河, 如南美的亚马逊河、北美的密西西比河、亚洲的长江、湄公河和黄河开展了pCO2FCO2的定量研究。例如, Li等[8]分析了湄公河下游42个水文站的数据得出水体表层pCO2均值为1090 μatm, FCO2为71 mol m-2 a-1;Li等[9]对长江龙川支流的研究发现该支流pCO2均值为1230 μatm, FCO2为27 mol m-2 a-1;密西西比河流域河流和湿地的CO2释放量可达1.2±0.3 MgC hm-2 a-1[10]

河流pCO2FCO2通常呈现出明显的时空差异, 导致差异的主要原因是河流水体中CO2的来源不相同, 并且相关的环境影响因素也有所不同。河流CO2的来源有外源和内源两种。外源主要有土壤中含碳有机物的矿化、土壤有机物的降解、植物根系的呼吸;内源主要有水体有机碳的原位呼吸作用、光化学降解作用以及碳酸盐的沉淀[11-12]。环境影响因素主要有温度、降水、流速、微生物以及人类活动等, 这些因素相互作用, 共同影响着河流CO2逸出[13]。研究表明流域降水给土壤提供了水分, 而微生物在适宜的温度条件下活跃度较高, 利于产生CO2进入河流, 从而增大了水体的pCO2。近些年人类活动引发的营养和热力学变化也可能增加河流pCO2[14]

关于FCO2的估算, 有静态箱法、模型估算法、涡度相关法以及气体示踪试验法等[15-17], 目前最常用的主要是前两种方法。静态箱法即红外水-气界面CO2分析仪法, 该法的优点是能直接在河流表面测得FCO2, 缺点是实验设备对水面风速和降雨可能存在扰动[18]。模型估算法又称梯度法, 表达式[17, 19]为:

(1)

式中, k为气体交换系数, pCO2水为河流表层的CO2分压, pCO2气为大气的CO2分压, 河流表层pCO2一般采用温度、pH结合碱度间接估算得到[20], 碱度则由滴定法得到。该模型中参数k的取值的不同, 结果有所差异, 因此存在争议。Peng采用模型估算法得到了乌江不同深度处的pCO2[21];Noriega和Araujo采用该法估算了巴西北部和东北部河口的pCO2FCO2 [22];Li等同样采用该法估算了湄公河下游的pCO2FCO2[8]

黄河是高含沙量河流的代表, 近些年有关其碳通量研究的报道越来越多。孙超[23]通过对黄河花园口和利津测站的观测, 研究水体pCO2与水沙的关系, 表明高泥沙含量对应较高的pCO2;王亮[24]观测了黄河花园口和利津测站水体中的各种形态的碳浓度, 并估算了河流碳的输运量和水-气界面通量。Ran等[25]对黄河干流碳逸出及其影响因素研究发现, 雨季FCO2大于旱季, 并且受人类影响, 输入海的碳通量在减少而FCO2在增加。以上研究表明目前对于黄河碳逸出的研究主要集中在黄河下游与河口地区, 而黄河中游和上游的相关研究较少。本文以黄河中游内蒙古段的头道拐断面为研究对象, 基于FCO2及其相关水文和水化学指标的监测分析, 探索了FCO2在年际、年内以及断面上的分布特征以及影响因素, 补充了黄河中上游CO2逸出规律的研究, 为黄河碳循环研究提供了科学参考。

1 研究区概况

黄河是中国第二大河, 全长约5464 km[26], 水面落差4480 m, 流域面积约752443 km2, 平均含沙量达37 kg/m3[27], 是世界上含沙量最高的河流。本研究选取黄河上游与中游分界点—头道拐水文站(110°04′E, 40°16′N)为试验点。头道拐水文站是黄河干流重要的水文站点之一, 位于内蒙古托克托县河口镇, 位置如图 1所示。该站点所在区域属温带干旱半干旱气候, 年降水量为155—366 mm, 75%的降水集中于7—9月份。1956—2010年的水面蒸发量为1368 mm, 夏季蒸发量为全年最大[28]。站点所在流域内自产径流很少, 但过境流量很大, 流经该测站的径流量占全河的50%以上[29]。1987—2010年的黄河径流资料显示, 头道拐多年平均水量为258亿m3, 主要集中于汛期的7—10月[30]。黄河在巴彦高勒与头道拐区间, 主要接纳“十大孔兑”等支流汇入, 这些支流流域植被较差、土质疏松, 并且河道流程短、比降大, 因此汛期极易造成水土流失[31]。另外, 头道拐水文站的上游是河套灌区, 年引黄水量约50亿m3, 占黄河过境水量的七分之一[32]

图 1 研究断面位置示意图 Fig. 1 The position of the river cross-section
2 研究方法

本研究在头道拐水文站附近选取研究断面, 如图 2所示。为探究FCO2沿河道断面的空间分布规律, 在研究断面布设4个采样点, S1点靠近黄河左岸, S2和S3点分别位于河道中部, S4点靠近右岸, 分布如图 2所示。自2013年至2015年, 共实测得到近3年的FCO2数据, 采样时间分别为2013年7月、11月, 2014年和2015年的4月、6月、9月和12月。

图 2 研究断面采样点布设示意图 Fig. 2 The sampling points of the river cross-section

本研究运用静态箱法—红外水气CO2分析仪Li-7000(美国Li-Cor公司)监测水体pCO2[33]。由于Li-7000无法直接测量水体中的CO2浓度, 因此需要外接静态箱(图 3), 通过两根塑料管将静态箱与Li-7000气体分析仪连接, 形成空气闭合回路。本试验所用静态箱为体积0.018 m3, 与水面接触面积为0.09 m2的塑料箱, 为使塑料箱浮在水面, 在塑料外部套有12 cm厚泡沫板。测量时, 首先将静态箱举起, 使箱内CO2浓度与大气CO2浓度平衡, 之后将静态箱放置于河流表面, 保持其静止。Li-7000气体分析仪每间隔1 min记录一次pCO2数值, 每次监测持续20 min左右。Li-7000气体分析仪采用12V电池供电。水-气界面FCO2可通过下式计算[34]

(2)
图 3 Li-7000静态箱法原理示意图 Fig. 3 The schematic of Li-7000 static chamber method

式中, FCO2为CO2逸出通量(μmol m-2 s-1);dpCO2/dtpCO2在静态箱中积累量的斜率(μatm/s), 一般采用开始7—15 min之内pCO2呈线性升高趋势的数值计算;V为静态箱体积(m3);R为气体常数(m3 atm K-1 mol-1);T为静态箱内温度(K);S为静态箱覆盖水面面积(m2)。

野外试验过程中, 除了对每一采样点监测水-气界面CO2浓度变化, 还监测pH、流速、风速和水温等指标。水温和pH采用便携式水质监测仪(Multi3420)监测, 水温测量精度±0.2℃, pH测量精度±0.004。河道流速采用直读式流速仪(Global Water FP211)进行监测(m/s), 测量精度±0.03 m/s, 测量时将流速仪放置于每个采样点水面以下10 cm处, 分别测量3次取平均值。风速(m/s)采用标智风速仪(GM8901)进行监测, 精准度为±3%, 测量高度基本与静态箱采样时的水平高度一致, 每个采样点测量两次取平均值。取水样时, 在每个采样点水面以下5 cm处用超纯水润洗过的500 mL采样瓶进行水样采集。

冬季采样时, 分别在4个采样点处用冰镐在冰面上砸开略大于静态箱底面积的冰洞。12月的河面结冰厚度为20—40 cm, 河中央S2和S3处的冰层较厚。之后将静态箱放置在冰洞中进行监测, 步骤和方法同上。

取回水样后, 溶解有机碳(Dissolved Organic Carbon, DOC)使用德国VORE公司的TOC总碳总氮分析仪测定。碱度在采样后的5 h之内测定, 甲基橙做指示剂, 用0.1 mol/L的盐酸对采集的水样进行滴定, 将得到的碱度结合pH估算各采样点水体pCO2。河流系统中总无机碳(Dissolved Inorganic Carbon, DIC)是碳酸氢盐(HCO3-)、碳酸(H2CO3)、碳酸盐(CO32-)、溶解二氧化碳(CO2aq)的总和。无机碳中各组分构成比例随温度和pH值变化而变化, 温度与pH同时还影响大气中CO2与水体CO2aq的平衡。无机碳可以通过亨利定律来计算[35-36]

(3)
(4)
(5)
(6)

式中H2CO3* 是H2CO3与CO2aq的总和, K为给定温度的亨利常数。因此, 河流水体无机碳用以下公式进行计算[37]

(7)
(8)
(9)

河流水体的pCO2计算可以简单的表示为下式:

(10)

目前研究认为HCO3-代表了99%的碱度, 所以可以直接用做碱度进行计算[38]

3 结果 3.1 相关环境因子

与研究断面水体FCO2相关的环境因子如表 1所示。水温的季节性变化明显, 在0—25℃范围内。pH的变化范围为8.01—8.45, 水体偏碱性。监测期间风速大部分为1—2 m/s, 仅在2014年12月风速最高达到6.3 m/s。研究断面的水体DOC浓度在2.7—13 mg/L, 其来源可能有多种, 分为内源和外源。内源主要是指水生植物通过光合作用吸收无机碳产生有机碳, 但是研究断面的水生植物较少[39], 而且夏季流速较快、泥沙含量较高和透光性较弱的特征限制了水生植物生长[40]。因此, 研究断面的DOC主要来自于外源输入, 包括土壤有机质的降解产物、人类生产和生活的排放物, 特别是工农业生产, 均是水体DOC不可忽略的来源。断面水体pCO2介于467—2101 μatm, 季节性变化显著, 丰水期的pCO2普遍大于枯水期, 最大值2101 μatm出现在2013年7月, 同时FCO2也达到最大510 mmol m-2 d-1

表 1 研究断面FCO2及其相关影响因素数据表 Table1 Data list of FCO2 and relative influencing factors at the river cross-section
日期Date水温/℃
Water
temperature
pH
风速
Wind velocity/
(m/s)
碱度
Alkalinity/
(μmol/L)
溶解有机碳
Dissolved
organic carbon/
(mg/L)
流速
Current
velocity/
(m/s)
pCO2
Partial pressure
of carbon
dioxide/(μatm)
FCO2
Carbon dioxide
evasion/
(mmol m-2 d-1)
2013-07258.011.332661.32101510
2013-1258.421.135530.569599
2014-04118.332.135175.40.8898152
2014-06238.372.233087.50.5946401
2014-09178.211.233206.01.51101504
2014-1208.326.345907.70.5115059
2015-04168.210.919037.10.8644145
2015-062513.00.5
2015-09168.301.833005.51.2957150
2015-1208.450.126732.70.246739
—表示未测得数据
3.2 FCO2的年际和年内季节变化

将每次采样时S1—S4点的FCO2取算术平均, 得到研究断面平均FCO2随时间变化过程, 如图 4所示。头道拐断面FCO2的年际变化总体趋同(图 4), 在2013—2015年期间, 冬季FCO2均为各年最低;春季FCO2在2014和2015年基本一致;2013和2014年的夏季FCO2均处于较高水平。2015年夏季FCO2缺测是因为河道在S1和S2点干涸。

图 4 研究断面的二氧化碳逸出通量(FCO2)年际和年内季节变化过程 Fig. 4 The annual and seasonal variations of carbon dioxide efflux (FCO2) at the river cross-section

研究断面的FCO2呈现出明显的年内季节性变化特征(图 4)。夏季FCO2为全年最高456 mmol m-2 d-1, 如果遇到枯水年份, 河道水量较低, FCO2也可能很低, 比如研究断面在2015年6月的河道流量极低。秋季次之为251 mmol m-2 d-1, 春季为148 mmol m-2 d-1, 而冬季为全年最低33 mmol m-2 d-1。一般认为冬季时, CO2会被冰层封于水体内, 因此水体可能储存了大量CO2, 但是该结果表明冬季冰层以下河流FCO2较低。

3.3 FCO2在研究断面的变化

将各采样点的FCO2在研究时段(2013年7月—2015年12月)内取算数平均值, 得到FCO2在头道拐断面上的空间分布(图 5)。结果表明研究断面的FCO2呈现出明显的空间差异性:靠近右岸的S4点处的FCO2为断面最大392 mmol m-2 d-1;左岸S1点处的最小为86 mmol m-2 d-1;位于河道中部S2和S3点处的FCO2基本相同, 分别为236 mmol m-2 d-1和237 mmol m-2 d-1, 介于S1和S4之间。

图 5 研究断面4个采样点的二氧化碳逸出通量(FCO2) Fig. 5 Carbon dioxide efflux (FCO2) at 4 sampling points at the river cross-section

综合研究断面的FCO2, 通过时间换算得到头道拐年均FCO2为84 mol m-2 a-1。该结果与其他河流的结果一致, 表明河流是CO2的源, 不断向大气释放CO2。相比于世界其他较大河流的FCO2:如湄公河下游的107.5 mol m-2 a-1、密西西比河的98.5 mol m-2 a-1、亚马逊河的69 mol m-2 a-1、哈德逊河的13.5 mol m-2 a-1以及长江河口的15.5—34.2 mol m-2 a-1, 该断面FCO2处于较高水平。对黄河而言, 位于上游的头道拐断面FCO2小于中游的125.6 mol m-2 a-1[5]

4 讨论

研究结果表明黄河头道拐断面是一个碳源, 而CO2逸出强度受诸多因素影响。从公式(1)可知, 大气的pCO2较稳定, 所以水体pCO2和气体交换系数k是控制FCO2的两个关键因素。

4.1 pCO2

根据温度、碱度和pH计算得到的头道拐断面水体pCO2普遍高于大气中的400 μatm:夏季pCO2为全年最大平均1524 μatm、秋季次之平均为1029 μatm、冬季和春季基本相同约为770 μatm, 同样表明黄河头道拐断面是一个碳源, 该结论与黄河其他站点的pCO2研究结果相一致[41-42]

黄河水体表层pCO2FCO2的季节变化基本一致。夏季FCO2为全年最高, 表层水体pCO2平均值为1524 μatm, ΔpCO2达到1124 μatm, 促使水体中的CO2向大气释放, 而冬季FCO2为全年最低, 测得水体pCO2在2015年12月最小仅为467 μatm(表 1), ΔpCO2较小导致FCO2较小。与FCO2结果不同的是, pCO2在冬季和春季非常相近。该断面的FCO2pCO2总体呈中等相关水平(R2=0.24, P < 0.01), 如图 6所示。分析原因, 可能是因为FCO2除了与pCO2有关, 还与气体交换系数k有关。例如, 位于S3点的2014年秋季FCO2高达448 mmol m-2 d-1, 但是相应的pCO2仅为1181 μatm, 而其河道流速却高达1.84 m/s。

图 6 研究断面水体二氧化碳逸出通量(FCO2)与二氧化碳分压(pCO2)线性拟合结果 Fig. 6 The result of carbon dioxide efflux (FCO2) vs partial pressure of carbon dioxide (pCO2) with a linear fitting

研究表明水体pCO2主要由水温、pH、碱度和CO2浓度决定[43]。CO2浓度和碱度主要受光合作用、呼吸作用和钙化作用等生物过程控制[17]

4.1.1 水温

采样期间测得断面水温为0—25℃, 将水温和pCO2进行线性拟合, 表明河水温度与pCO2呈较好的正相关关系, 如图 7所示。夏季水温较高, 提高了水体中微生物的活性, 并降低了CO2溶解度, 因此FCO2相应增加, 冬季则反之[25, 26, 29]。Alin等[44]在研究亚马逊河时同样发现水温高会促进水-气界面的气体交换速率, 并增加CO2逸出强度。桂祖胜[45]在研究长江水-气界面CO2逸出季节变化时发现夏季FCO2为全年最高、秋季次之, 与本研究结果不同的是其冬季FCO2高于春季, 该差异可能是因为长江冬季无结冰现象。

图 7 研究断面水温、pH、溶解无机碳、溶解有机碳与二氧化碳分压的拟合结果 Fig. 7 The fitting results of water temperature, pH, dissolved inorganic carbon, dissolved organic carbon vs partial pressure of carbon dioxide at the river cross-section
4.1.2 pH

由水体中二氧化碳-碳酸盐体系的平衡过程可知(公式3)[46], pH升高会促使平衡向右移动, 更多的CO2生成碳酸盐溶解在水中, 导致水体pCO2降低, 水-气界面ΔpCO2减小, 所以FCO2减少。将采样点的pH和pCO2进行线性拟合, 得到pCO2与pH呈负相关性, 如图 7所示。研究表明当水体pH>8.3时, 河流相当于碳汇, 而当pH < 8.3时, 河流相当于碳源[9]。本研究断面超过一半水样的pH小于8.3, 同样说明该断面水体是CO2的源。

4.1.3 DIC和DOC

黄河水体中DIC主要由HCO3-构成, HCO3-又代表了99%的碱度[47]。黄河流域的黄土和灰钙土分布较广泛, 化学风化导致水体碱度较高[25], 研究断面的水体碱度最高可达4590 mg/L(表 1)。碱度直接影响水体pCO2, 但是图 7中的碱度和pCO2相关性并不好, 这点与前面FCO2pCO2呈中等相关的结果相符合。DOC和pCO2呈现出较好的正相关关系(R2=0.50, P < 0.001), 与相关研究结果符合, 即DOC的分解是控制水体pCO2的重要因素[42-43]

4.2 气体交换系数k

研究表明气体交换系数kFCO2的影响甚至大于水体pCO2, k值的大小主要由流速、风速控制, 同时还受到水深、河宽和河床坡度等多种环境因子的影响[48]

4.2.1 流速

头道拐断面的流速变化范围为0.2—1.5 m/s, S4点处的平均流速是断面最大1.04 m/s, S3点为0.99 m/s, S2点为0.81 m/s, S1点处最小仅为0.48 m/s。S4点流速较大因为该点位于河道的中泓线, 水深也最大。断面4个采样点处的平均流速大小分布规律与FCO2相似。将采样点的流速和FCO2进行线性拟合, 二者呈现较好的正相关关系, R2可达0.7, 如图 8所示。本研究结果说明头道拐断面河道流速是FCO2的关键影响因素, 大于pCO2对河道FCO2的贡献(R2=0.24)。这主要是因为流速较大有助于增大河流表面的湍流度和破碎度, 使水体与空气的接触面积变大, 从而加速了两者间气体的交换[49-50]

图 8 研究断面流速和风速与二氧化碳逸出通量的线性拟合结果 Fig. 8 The results of current velocity and wind speed vs carbon dioxide efflux with linear fitting at the study river section

夏季FCO2普遍较高, 降雨集中和流量较大是主要影响因素之一[38], 实测流量最高可达1010 m3/s(2014年8月21日)。冬季FCO2为全年最小, 与冰层以下较低的河水流速(0.13 m/s—0.3 m/s)有关, 较小的水体扰动减缓了CO2释放[50-51]

4.2.2 风速

已有研究表明风速是影响气体交换速率常数的主导因素, 对气体交换起着决定性作用。风会加快气体的交换速率, 带走河流水-气界面处的CO2, 从而增大河流水体与大气之间的ΔpCO2, 促使水体中的CO2逸出。张龙军等[52]在研究黄河碳通量影响因素时发现, 风速是影响FCO2的一个重要因素, 且两者表现出明显的正相关性。

本研究将每个采样点的风速和FCO2进行线性拟合, 表明FCO2与风速并没有呈现明显的相关性, 如图 8所示。可能是因为采样期间没有遇到大风天气, 大部分风速低于3.5 m/s。吕东珂[53]在研究水库FCO2时发现当风速小于3.5 m/s时, 对FCO2的影响较小, 而当风速大于3.5 m/s时, FCO2与风速就表现出明显的相关性。Raymond等[54]在哈德逊河的研究也发现FCO2未明显受到风速的影响, 同样可能是因为监测时风速较小, 即使出现风速较大的情况, 持续时间也并不长。

5 结论

本研究在2013—2015期间, 利用Li-7000静态箱法对黄河头道拐断面进行了FCO2野外监测试验, 分析了FCO2在年际、年内和断面上的变化规律, 并探讨了水体pCO2和气体交换系数kFCO2的影响。结果表明河流水体pCO2高于大气, 断面平均pCO2为995 μatm, 且表现出明显的季节性变化, 夏季平均pCO2最大1524 μatm, 秋季次之1029 μatm, 冬季与春季基本相同, 约为770 μatm。另外, DOC的分解是控制水体pCO2的重要因素, 而研究断面的DOC主要来自于外源。研究断面平均FCO2为84 mmol m-2 a-1, 呈现出与pCO2相似的年际变化规律和明显的季节性差异:夏季最高为456 mmol m-2 d-1, 与夏季温度高、流量大有关, 秋季次之251 mmol m-2 d-1, 春季148 mmol m-2 d-1, 冬季最低为33 mmol m-2 d-1, 可能是因为冬季温度低、微生物活性较低, 导致冬季冰层以下河流FCO2较低。研究断面FCO2表现出较强的空间差异性:靠近右岸S4点的FCO2为断面最大392 mmol m-2 d-1, 河段中部S2点和S3点的FCO2基本相同约为237 mmol m-2 d-1, 靠左岸S1点的FCO2最小86 mmol m-2 d-1。该结果表明即使在同一断面FCO2也可能存在较大差异, 研究河段流速最大处的FCO2也最大, 随机采样可能低估或高估水体CO2的逸出量, 因此在河道取样监测时需要选取具有代表性的采样点, 特别是在较大的河流中。研究断面的FCO2与流速和温度呈较好的正相关关系, 与pCO2的相关关系中等, 与pH呈负相关关系, 与风速并没有呈现明显的相关性。结果说明河道流速较水体pCO2FCO2的贡献更大。本研究为黄河中上游CO2逸出量的评价提供了科学参考, 但是仅一个断面不足以全面评价黄河的CO2释放能力, 需要在后续的试验中继续完善对黄河上中下游逸出CO2的研究。

参考文献
[1] 姚冠荣, 高全洲. 河流碳循环对全球变化的响应与反馈. 地理科学进展, 2005, 24(5): 50–60. DOI:10.11820/dlkxjz.2005.05.006
[2] 孙会国, 汉景泰, 张淑荣, 吕喜玺. "05.06"西江特大洪水对河流碳输出通量的影响. 科学通报, 2006, 51(23): 2773–2779. DOI:10.3321/j.issn:0023-074X.2006.23.012
[3] 高全洲, 沈承德. 河流碳通量与陆地侵蚀研究. 地球科学进展, 1998, 13(4): 369–375.
[4] Ran L S, Wang S J, Fan X L. Channel change at Toudaoguai Station and its responses to the operation of upstream reservoirs in the upper Yellow River. Journal of Geographical Sciences, 2010, 20(2): 231–247. DOI:10.1007/s11442-010-0231-9
[5] 杨欢. 黄河中游pCO2的时空变化特征研究[D]. 呼和浩特: 内蒙古大学, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10126-1015364344.htm
[6] 田丽欣, 吴莹, 林晶, 朱卓毅, 张经. 冬季海南省东北部河流及近岸海区有机碳的分布特征. 热带海洋学报, 2010, 29(5): 136–141. DOI:10.11978/j.issn.1009-5470.2010.05.136
[7] 王敏. 长江主流碳的时空输运特征及三峡工程对其影响[D]. 青岛: 中国海洋大学, 2010. http://cdmd.cnki.com.cn/Article/CDMD-10423-1011027104.htm
[8] Li S Y, Lu X X, Bush R T. CO2 partial pressure and CO2 emission in the lower Mekong river. Journal of Hydrology, 2013, 504: 40–56. DOI:10.1016/j.jhydrol.2013.09.024
[9] Li S Y, Lu X X, He M, Zhou Y, Li L, Ziegler A D. Daily CO2 partial pressure and CO2 outgassing in the upper Yangtze River basin:a case study of the Longchuan River, China. Journal of Hydrology, 2012, 466-467: 141–150. DOI:10.1016/j.jhydrol.2012.08.011
[10] Richey J E, Melack J M, Aufdenkampe A K, Ballester V M, Hess L L. Outgassing from Amazonian rivers and wetlands as a large tropical source of atmospheric CO2. Nature, 2002, 416(6881): 617–620. DOI:10.1038/416617a
[11] Abril G, Martinez J M, Artigas L F, Moreira-Turcq P, Benedetti M F, Vidal L, Meziane T, Kim J H, Bernardes M C, Savoye N, Deborde J, Souza E L, Albéric P, de Souza M F L, Roland F. Amazon river carbon dioxide outgassing fuelled by wetlands. Nature, 2014, 505(7483): 395–398.
[12] Telmer K, Veizer J. Carbon fluxes, pCO2 and substrate weathering in a large northern river basin, Canada:carbon isotope perspectives. Chemical Geology, 1999, 159(1/4): 61–86.
[13] Ran L S, Lu X X, Sun H G, Han J T, Yu R H. Chemical denudation in the Yellow River and its geomorphological implications. Geomorphology, 2014, 231: 83–93.
[14] 于奭, 杜文越, 孙平安, 黄婕, 罗惠先, 李幼玲, 何师意. 亚热带典型河流水化学特征、碳通量及影响因素. 水文, 2015, 35(4): 33–41.
[15] 姚骁, 李哲, 郭劲松, 林初学, 陈永柏, 李翀. 水-气界面CO2通量监测的静态箱法与薄边界层模型估算法比较. 湖泊科学, 2015, 27(2): 289–296. DOI:10.18307/2015.0213
[16] 赵炎, 曾源, 吴炳方, 陈永柏, 王强, 袁超. 水库水气界面温室气体通量监测方法综述. 水科学进展, 2011, 22(1): 135–146.
[17] 丁虎, 刘丛强, 郎赟超, 李思亮, 李晓东, 汪福顺. 河流水-气界面碳交换研究进展及趋势. 上海大学学报:自然科学版, 2015, 21(3): 275–285.
[18] Ran L S, Lu X X, Sun H G, Han J T, Li R H, Zhang J M. Spatial and seasonal variability of organic carbon transport in the Yellow River, China. Journal of Hydrology, 2013, 498: 76–88. DOI:10.1016/j.jhydrol.2013.06.018
[19] Müller D, Warneke T, Rixen T, Müller M, Jamahari S, Denis N, Mujahid A, Notholt J. Lateral carbon fluxes and CO2 outgassing from a tropical peat-draining river. Biogeosciences Discussions, 2015, 12(13): 10389–10424. DOI:10.5194/bgd-12-10389-2015
[20] Hotchkiss E R, Hall Jr R O, Sponseller R A, Butman D, Klaminder J, Laudon H, Rosvall M, Karlsson J. Sources of and processes controlling CO2 emissions change with the size of streams and rivers. Nature Geoscience, 2015, 8(9): 696–699. DOI:10.1038/ngeo2507
[21] Peng X, Wang B L, Liu C Q, Liu X L, Wang F S. Diurnal variations of pCO2 in relation to environmental factors in the cascade reservoirs along the Wujiang River, China. Chinese Journal of Geochemistry, 2012, 31(1): 41–47. DOI:10.1007/s11631-012-0547-5
[22] Noriega C, Araujo M. Carbon dioxide emissions from estuaries of northern and northeastern Brazil. Scientific Reports, 2014, 4(4): 6164–6164.
[23] 孙超. 黄河调水调沙时期的碳输运特征及花园口站碳系统各分量的年变化[D]. 青岛: 中国海洋大学, 2007. http://cdmd.cnki.com.cn/article/cdmd-10423-2008012444.htm
[24] 王亮. 黄河干流碳输运及人类活动对其影响[D]. 青岛: 中国海洋大学, 2014. http://cdmd.cnki.com.cn/Article/CDMD-10423-1014369034.htm
[25] Ran L S, Lu X X, Yang H, Li L L, Yu R H, Sun H G, Han J T. CO2 outgassing from the Yellow River network and its implications for riverine carbon cycle. Journal of Geophysical Research:Biogeosciences, 2015, 120(7): 1334–1347. DOI:10.1002/2015JG002982
[26] Wang Z Yu, Xin Y Z, Gao D M, Li F M, Morgan J, Xing B S. Microbial community characteristics in a degraded wetland of the yellow river delta. Pedosphere, 2010, 20(4): 466–478. DOI:10.1016/S1002-0160(10)60036-7
[27] Ren M E. Sediment discharge of the Yellow River, China:past, present and future-a synthesis. Acta Oceanologica Sinica, 2015, 34(2): 1–8. DOI:10.1007/s13131-015-0619-6
[28] 丁磊. 黄河内蒙古段河道变迁和干流水量损失研究[D]. 呼和浩特: 内蒙古农业大学, 2012. http://cdmd.cnki.com.cn/Article/CDMD-10129-1013153776.htm
[29] 张建云, 王国庆, 贺瑞敏, 刘翠善. 黄河中游水文变化趋势及其对气候变化的响应. 水科学进展, 2009, 20(2): 153–158.
[30] 张永领, 董玉龙, 张东. 在三门峡水库影响下黄河有机碳的输送特征. 地理科学, 2015, 35(7): 912–918.
[31] 胡兴林, 王静, 蓝永超, 李文治, 赵昌瑞. 黄河内蒙段河道冲淤量与水沙输移过程的关系分析. 水文, 2012, 32(2): 44–48.
[32] Ran L S. Recent Riverine Carbon of the Yellow River:Fluxes, Outgassing and Burial[D]. Singapore:National University of Singapore, 2013. http://scholarbank.nus.edu.sg/handle/10635/48347?show=full
[33] Bergeron O, Margolis H A, Black T A, Coursolle C, Dunn A L, Barr A G, Wofsy S C. Comparison of carbon dioxide fluxes over three boreal black spruce forests in Canada. Global Change Biology, 2007, 13(1): 89–107. DOI:10.1111/gcb.2007.13.issue-1
[34] Frankignoulle M. Field measurements of air-sea CO2 exchange. Limnology and Oceanography, 1988, 33(3): 313–322. DOI:10.4319/lo.1988.33.3.0313
[35] Stumm W, Morgan J J. Aquatic chemistry:an introduction emphasizing chemical equilibria in natural waters.//Aquatic chemistry:An introduction emphasizing chemical equilibria in natural waters. 1981:795-795.
[36] Plummer L N, Busenberg E. The solubilities of calcite, aragonite and vaterite in CO2-H2O solutions between 0 and 90℃, and an evaluation of the aqueous model for the system CaCO3-CO2-H2O. Geochimica et Cosmochimica Acta, 1982, 46(6): 1011–1040. DOI:10.1016/0016-7037(82)90056-4
[37] Yao G R, Gao Q Z, Wang Z G, Huang X K, He T, Zhang Y L, Jiao S L, Ding J. Dynamics of CO2, partial pressure and CO2 outgassing in the lower reaches of the Xijiang River, a subtropical monsoon river in China. Science of the Total Environment, 2007, 376(1/3): 255–266.
[38] Hunt C W, Salisbury J E, Vandemark D. Contribution of non-carbonate anions to total alkalinity and overestimation of pCO2 in new England and new Brunswick rivers. Biogeosciences, 2011, 8(10): 3069–3076. DOI:10.5194/bg-8-3069-2011
[39] 袁永锋, 李引娣, 张林林, 余斌, 齐喜荣, 韩明轩, 张建军. 黄河干流中上游水生生物资源调查研究. 水生态学杂志, 2009, 2(6): 15–19.
[40] 徐雪梅. 秋季黄河主流碳的输运特征[D]. 青岛: 中国海洋大学, 2008. http://cdmd.cnki.com.cn/Article/CDMD-10423-2008174986.htm
[41] Li X L, Jia Q Y, Liu J M. Seasonal variations in heat and carbon dioxide fluxes observed over a reed wetland in northeast China. Atmospheric Environment, 2015, 127: 6–13.
[42] 张龙军, 徐雪梅, 温志超. 秋季黄河pCO2控制因素及水-气界面通量. 水科学进展, 2009, 20(2): 227–235.
[43] 黄文敏, 朱孔贤, 赵玮, 余博识, 袁希功, 冯瑞杰, 毕永红, 胡征宇. 香溪河秋季水-气界面温室气体通量日变化观测及影响因素分析. 环境科学, 2013, 34(4): 1270–1276.
[44] Alin S R, de Fátima F L, Rasera M, Salimon C I, Richey J E, Holtgrieve G W, Krusche A V, Snidvongs A. Physical controls on carbon dioxide transfer velocity and flux in low-gradient river systems and implications for regional carbon budgets. Journal of Geophysical Research:Biogeosciences, 2011, 116(G1): G01009.
[45] 桂祖胜. 长江干流及其河口水体pCO2的分布与影响因素的初步研究[D]. 青岛: 中国海洋大学, 2007. http://cdmd.cnki.com.cn/Article/CDMD-10423-2008011373.htm
[46] Larsen S, Andersen T, Hessen D O. The pCO2 in boreal lakes:organic carbon as a universal predictor?. Global Biogeochemical Cycles, 2011, 25(2): GB2012.
[47] Ran L S, Lu X X, Richey J E, Sun H G, Han J T, Yu R H, Liao S Y, Yi Q. Long-term spatial and temporal variation of CO2 partial pressure in the Yellow River, China. Biogeosciences, 2015, 12(4): 921–932. DOI:10.5194/bg-12-921-2015
[48] Wallin M B, Öquist M G, Buffam I, Billett M F, Nisell J, Bishop K H. Spatiotemporal variability of the gas transfer coefficient (KCO2) in boreal streams:implications for large scale estimates of CO2 evasion. Global Biogeochemical Cycles, 2011, 25(3): GB3025.
[49] Alin S R, Salimon C I, Richey J E, Gordon I S, Jeffrey E R, Gordon W H, Alex V K, Anond K. Physical controls on carbon dioxide transfer velocity and flux in low-gradient river systems and implications for regional carbon budgets. Journal of Geophysical Research Biogeosciences, 2011, 116(G1): G01009.
[50] Zappa C J, Raymond P A, Terray E A, McGillis W R. Variation in surface turbulence and the gas transfer velocity over a tidal cycle in a macro-tidal estuary. Estuaries, 2003, 26(6): 1401–1415. DOI:10.1007/BF02803649
[51] Wang F S, Wang Y C, Zhang J, Xu H, Wei X G. Human impact on the historical change of CO2 degassing flux in River Changjiang. Geochemical Transactions, 2007, 8: 7. DOI:10.1186/1467-4866-8-7
[52] 张龙军. 东海海-气界面CO2通量研究[D]. 青岛: 中国海洋大学, 2003. http://d.wanfangdata.com.cn/Thesis/Y521276
[53] 吕东珂. 哈尔滨周边泥炭型水库水-气界面CO2通量研究[D]. 哈尔滨: 东北林业大学, 2013. http://cdmd.cnki.com.cn/Article/CDMD-10225-1013359259.htm
[54] Raymond P A, Caraco N F, Cole J J. Carbon dioxide concentration and atmospheric flux in the Hudson River. Estuaries, 1997, 20(2): 381–390. DOI:10.2307/1352351