环境科学  2018, Vol. 39 Issue (11): 5217-5226   PDF    
三峡库区典型河流水-气界面CO2通量日变化观测及其影响因素分析
罗佳宸1,2, 李思悦1     
1. 中国科学院重庆绿色智能技术研究院, 重庆 400074;
2. 中国科学院大学, 北京 100049
摘要: 为研究河流水-气界面CO2通量的季节和日变化特征;于2016年7月15~17日以及2017年11月4~6日对三峡库区嘉陵江支流竹溪河进行定点定时采集表层水样,并同步监测关键环境因子,采用亨利定律结合薄边界层模型计算其水-气界面CO2通量F(CO2).结果表明,竹溪河表层水CO2分压p(CO2)及界面CO2脱气通量呈现出显著的日间和季节变化,以及明显的日内变化特征:在上午09:00前后达到释放高峰,随后波动下降;水-气界面CO2通量日间均值分别为(100.9±31.6)、(78.6±12.1)、(83.9±29.7)、(137.5±42.1)、(147.6±34.0)、(132.4±21.7)mmol·(m2·d)-1;并表现出夏季表层水体CO2释放通量明显低于秋季,其均值分别为(87.8±27.5)mmol·(m2·d)-1和(139.2±34.0)mmol·(m2·d)-1;总体表现出大气CO2源的特征.竹溪河p(CO2)和F(CO2)受到诸多环境因子的影响,相关分析表明,pH、碱度、水温和气温是主要环境影响因子,CO2释放通量可以用pH和碱度预测.
关键词: CO2释放      水-气界面      日变化      季节变化      环境因子     
Daily Variation of CO2 Flux at Water-Air Interface and Analysis of Its Affecting Factors in a Typical River of the Three Gorges Reservoir
LUO Jia-chen1,2 , LI Si-yue1     
1. Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences, Chongqing 400074, China;
2. University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: Diurnal and seasonal characterization of CO2 partial pressure p(CO2) and CO2 areal flux F(CO2) at the water-air interface in an anthropogenic river in the Three Gorges Reservoir area was studied. A tributary of the Jialing River in Chongqing Municipality was chosen, and daily and seasonal samples were taken in summer and autumn, focusing on riverine p(CO2), F(CO2), and their associated controls. Henry's law combined with the thin boundary layer model was adopted to estimate the CO2 flux via the water-air interface. The results indicated that F(CO2) was not high on average, namely (87.8±27.5) mmol·(m2·d)-1 and (139.2±34.0) mmol·(m2·d)-1 in summer and autumn, respectively. The water-air interface F(CO2) showed significant hourly, daily, and seasonal variations. CO2 release peaked around 09:00 and then slightly decreased. We also found that pH, alkalinity, water, and temperature were significantly related to p(CO2) and F(CO2), whereas pH and alkalinity were the best predictors of F(CO2). This study aids understanding of the impacts of urbanization on CO2 emissions in the rivers and helps to re-evaluate local riverine CO2 budgets.
Key words: CO2 emission      water-air interface      daily change      seasonal change      environmental factors     

河流、湖泊、水库等陆地地表水域生态系统是大气CO2的重要来源[1~5]; Cole等[6]发表评论称, 全球河流对大气CO2的贡献量(以C计)约为0.35 Gt·a-1; 2013年, Raymond等[7]重新估算全球内陆水体CO2释放通量, 发现河流水-气界面CO2通量(以C计)为1.80 Gt·a-1.如此大的评估差异主要来自野外特别是小河流观察数据的缺乏; 因此, 开展典型河流的水-气界面碳释放研究, 量化河流水-气界面CO2脱气通量, 对评估河流碳循环过程及在区域碳循环中的作用具有重要的科学意义.

目前, 国外关于内陆水体水-气界面CO2通量的研究主要在南美的亚马逊河和北美的密西西比河等大江大河开展的定量工作[8~11]; 研究表明大多数河流溶解性CO2呈饱和状态; 并发现世界范围内河流在水体有机质浓度、pH、营养状况、水温和风速等方面存在很大的差别, 因此不同河流水-气界面CO2脱气通量F(CO2)差异较大.我国相比国外, 对水-气界面CO2通量的研究主要集中在长江、黄河等[12~15].例如, Li等[12]对长江上游龙川江观测, 结果表明水体p(CO2)为23.2~838.3 Pa, 平均值为124.2 Pa, 是大气CO2浓度的3~4倍. Ran等[13]报道了黄河上游水体的CO2分压, 介于20.2~3 030.0 Pa, 均值为283.8 Pa, 为大气CO2浓度的7倍.

以上研究表明, 虽然我国学者在河流水-气界面碳交换通量的观测上取得较大进展; 但碳循环伴生的水-气界面CO2通量研究主要建立在年或季这样大时间尺度上, 以日变化的小时间尺度上的定量观测十分有限.由于河流受人类活动例如城镇化影响扰动大, 相对于自然河流其温室气体排放更显著; 昼间尺度上的环境因子与水-气界面CO2通量可能会有不同的表现形式, 同时昼间通量观测可以辅助对河流的碳源/汇作用的评估.通过定量观测河流CO2释放通量, 以期评价人类活动对河流温室气体的影响[16].本研究于2016年7月15~17日以及2017年11月4~6日, 定点定时对位于重庆主城区段的嘉陵江支流竹溪河水-气界面CO2通量连续日间变化观测, 以及关键水化学和水质指标的监测; 分析pCO2F(CO2)的季节和日变化特征及其影响因素, 以期为今后深入认识小时间尺度河流水-气界面碳交换对大气CO2贡献的精确估算提供重要的数据支撑, 同时为河流界面CO2通量对该区域城镇化进程的响应提供基础数据.

1 研究地区

本研究根据前期考察, 观测点位设在嘉陵江支流竹溪河下游龙门大桥附近(106°32′19″E, 29°48′24″N, 图 1), 离汇入嘉陵江河口距离4 km.竹溪河源于重庆渝北区华泰乡境内, 全长65.05 km, 流经两江新区复兴片区全境, 最后汇入嘉陵江.目前, 竹溪河正面临城镇化进程, 加剧流域土地利用方式的转变, 引起诸如植被毁损和水土流失等生态环境问题.研究点位地处亚热带大陆性季风气候, 雨量充沛, 多年平均降水量为1071.5 mm; 年平均气温为18℃左右, 最热月多为8~9月, 最冷月多为12月~次年2月.

图 1 三峡库区竹溪河采样点示意 Fig. 1 Distribution of sampling sites of the Zhuxi River in the Three Gorges area

2 材料与方法 2.1 样品采集与分析方法

样品采集于2016年7月15~17日(夏季)和2017年11月4~6日(秋季), 采样期间天气晴好, 避免了降雨造成水面波动的影响.具体的采样时间分为2个时段包括08:00~12:00和14:00~18:00;每间隔10 min使用CyberScan PCD 650多参数水质分析仪(Eutech, 美国)测定pH值、水温、溶解氧(DO)和电导率(EC); 并同步使用Testo 410-1便携式叶轮风速仪(德图, 德国)测定气温和风速.每间隔10 min用200 mL聚乙烯瓶采集10 cm表层水样运回实验室进行化学分析.

于采样当天用0.70 μm Whatman GF/F玻璃纤维膜过滤水样, 将滤液分装于去离子水润洗的聚乙烯瓶中, 不留空隙, 用于溶解性总氮(TDN)、溶解性总磷(TDP)测定和溶解性有机碳(DOC)测定.立即取25.00 mL过滤水样用0.020 0 mol·L-1盐酸滴定碱度; 表层水体中p(CO2)根据碱度、pH及水温求得.

TDN和TDP采用国标法测定; TDN采用碱性过硫酸钾消解紫外分光光度法(GB 11894-89)测定, TDP采用钼酸铵分光光度法(GB 11893-89)测定[17]. DOC采用multi N/C 2100S总有机碳分析仪(德国耶拿, 德国)测定.

2.2 二氧化碳分压计算

由测得的pH、碱度和水温, 计算出水溶液中的CO2分压p(CO2)[12]:

(1)
(2)
(3)
(4)

式中, K0K1K2为平衡常数.由于水体的pH均大于7.7, HCO3-视为与碱度相当[12].

水-气界面的CO2脱气通量主要受大气和水体中二氧化碳分压差和风速等因素影响.根据Fick[18]定律, CO2界面通量可以用以下公式[18, 19]:

(5)
(6)
(7)
(8)
(9)

式中, Flux为水-气界面CO2交换通量[mmol·(m2·d)-1], K为气体扩散系数(cm·h-1), Kh为亨利常数, 即气体溶解度(mol·L-1·atm-1), pwater为CO2在水中分压(μatm), pair为现场温度及压力下CO2气体在空气中分压(μatm), Tk为水体绝对温度(K), K600是在20℃下CO2调节至Schmidt常数等于600的气体传输速度(cm·h-1), Sc为t℃下CO2的Schmidt常数, w为流速(cm·s-1)[19].

2.3 数据处理方案

使用SPSS 18.0进行配对T检验、Spearman非参数相关性分析和多元回归分析; 并利用SigmaPlot 13.0作图.

3 结果与分析 3.1 水体p(CO2)和CO2通量

本研究表明竹溪河表层水体p(CO2)和水-气界面CO2通量在小时(日内)、日间和季节尺度上均具有显著性差异(P < 0.01), 见表 1.

表 1 日内、日间和季节各指标配对T检验1) Table 1 Intra-day, daytime, and seasonal pair T-tests for each indicator

监测期间, 竹溪河表层水CO2分压均大于大气CO2分压, 表现出在上午09:00左右达最大值, 随后波动下降至下午18:00左右达到最小的日内变化规律(图 2). 2016年07月15~17日以及2017年11月4~6日河流表层水日间均值分别为(88.0±11.4)、(80.7±4.5)、(83.8±11.9)、(93.5±13.1)、(96.9±10.5)和(91.8±6.3) Pa.夏季表层水p(CO2)介于54.7~115.1 Pa之间, 平均值为(86.2±7.7) Pa; 秋季表层水体p(CO2)介于60.4~130.0 Pa之间, 平均值为(93.4±8.9) Pa; 总体来看, 秋季表层水体p(CO2)明显高于夏季, 在统计学上具有显著性差异(P < 0.01).

图 2 p(CO2)和CO2脱气通量日变化 Fig. 2 Daily changes in p(CO2) and CO2 flux

监测期间, 竹溪河水-气界面CO2脱气通量均为正值, 表现出上午09:00左右达最大值, 随后波动下降至下午18:00左右达到最小释放量的日内变化规律(图 2). 2016年7月15~17日以及2017年11月4~6日水-气界面CO2通量日间均值分别为(100.9±31.6)、(78.6±12.1)、(83.9±29.7)、(137.5±42.1)、(147.6±34.0)、(132.4±21.7) mmol·(m2·d)-1.夏季CO2释放通量在10.6~175.5 mmol·(m2·d)-1范围间波动, 平均值为(87.8±27.5) mmol·(m2·d)-1; 秋季水-气界面CO2通量在31.3~247.5 mmol·(m2·d)-1范围间波动, 平均值为(139.2±34.0) mmol·(m2·d)-1.总体来看, 秋季表层水体CO2脱气通量明显高于夏季(P < 0.01);可能是由于秋季表层水体温度相对较低, 太阳辐射较弱, 水体浮游植物光合作用减弱导致的.

3.2 主要环境因子 3.2.1 气象因子

监测期间的水温、气温和风速均具有明显的日内变化特征(图 3), 在小时、日间和季节尺度上均具有显著性差异(P < 0.01).竹溪河表层水温在持续12 h的监测中日内变化规律很一致, 总体表现出下午水温高于上午, 具有显著性差异(P < 0.01). 2016年7月15~17日以及2017年11月4~6日其表层水温日间平均值分别为(25.0±0.4)、(26.2±0.8)、(27.4±0.7)、(18.7±0.5)、(19.0±0.4)、(18.6±0.5)℃.夏季表层水温介于24.4~28.4℃之间, 均值为(26.2±1.1)℃; 秋季表层水温明显低于夏季, 具有显著性差异(P < 0.01), 变化范围为18.0~19.5℃, 平均值为(18.7±0.5)℃.

图 3 环境因子日变化 Fig. 3 Daily changes in main environmental factors

受太阳辐射的影响, 气温波动较大, 是一天中变化最明显的环境因子; 日内变化规律一致, 表现出先上升至下午15:00左右达最大值后波动下降的趋势. 2016年07月15~17日以及2017年11月4~6日其日间平均气温分别为(28.7±2.2)、(30.2±3.1)、(33.2±2.0)、(21.3±2.3)、(21.4±3.7)、(20.3±2.6)℃.夏季气温在24.9~35.8℃范围内波动, 平均为(30.7±3.1)℃, 秋季气温显著低于夏季(P < 0.01), 变化范围在15.6~28.6℃之间, 平均为(21.0±3.0)℃.

风速波动性较大, 夏季和秋季风速分别在0~2.2 m·s-1和0~1.0 m·s-1范围内变化, 其中受盛行季风影响夏季风速明显高于秋季, 具有显著性差异(P < 0.01).

3.2.2 水化学参数

监测期间的pH、EC在小时、日间和季节尺度上均具有显著性差异(P < 0.01, 图 3).竹溪河表层水体的pH均呈碱性, 在持续12 h的监测中上午表层水体pH明显低于下午, 具有显著性差异(P < 0.01). 2016年7月15~17日以及2017年11月4~6日竹溪河表层水体pH日间值分别在8.09~8.32、8.30~8.38、8.28~8.58、8.25~8.56、8.22~8.39和8.28~8.39范围内波动.夏季和秋季其pH分别在8.09~8.58和8.22~8.56范围内波动. EC表明了水体中的离子强度, 在持续12 h的监测中上午表层水体EC明显低于下午, 具有显著性差异(P < 0.01). 2016年7月15~17日以及2017年11月4~6日其表层水体EC日间值分别在354.0~376.0、411.0~442.0、494.0~519.0、660.90~669.9、667.4~688.8和667.5~677.4 μS·cm-1范围内波动.夏季和秋季EC分别在354.0~519.0和660.9~688.8 μS·cm-1范围内波动, 秋季表层水体EC远高于夏季, 具有显著性差异(P < 0.01);可能是由于夏季雨水充沛导致对河流的稀释作用造成的.

图 4可知, 竹溪河秋季表层水体TDN、TDP和DOC浓度分别在0.61~1.89、0.10~0.20和3.30~4.55 mg·L-1范围内波动; 其中营养物质波动性较大, 可能是由于微生物的代谢及外界营养盐物质的输入不同造成的; 从TDN、TDP和DOC营养物质来判断, 监测期间营养水平较低.

图 4 营养物质日变化 Fig. 4 Daily changes in riverine nutrients

4 讨论 4.1 竹溪河表层水体p(CO2)格局的影响因子分析

单因素方差分析表明, 在小时、日间和季节尺度下竹溪河表层水体p(CO2)均具有显著性差异(P < 0.01).本研究定点定时对位于重庆主城区段的嘉陵江支流竹溪河表层水体连续日间变化观测, 所处环境的风速、气压等都大致相同; 竹溪河正面临城镇化进程, 人为活动引起水体有机质浓度、水体营养状况、pH和水温等其他因素的差异可能是造成表层水在小时、日间和季节尺度下表层水体p(CO2)具有显著性差异的最重要原因.

有研究认为河流p(CO2)是流域侵蚀、有机物的降解及从水体到大气的CO2脱气等过程平衡的结果[8~11].水温、pH和溶解性CO2等因素共同影响水体p(CO2)分布, 一般河流水体中p(CO2)主要受到水中碳酸盐平衡及浮游植物的生化作用控制.在已有的河流水体CO2影响因素的研究中, 河流侵蚀带来的有机碳的呼吸作用是维持亚马逊河p(CO2)的主要机制[8], 英国泰晤士河表层水体p(CO2)受生物呼吸作用和光合作用控制[20].

本研究通过p(CO2)与环境因子的Spearman相关性分析(表 2), 探讨其主要环境因子, 研究表明其表层水p(CO2)与pH、水温和气温呈显著负相关关系.造成这种现象可能是监测期间pH的变化范围比较小, 且绝大分时间在弱碱环境下, 由于水体溶解性有机碳原位呼吸过程消耗氧气, 产生CO2, 进而降低pH; 同时, 水体中进行光合作用的生物及植物吸收CO2, pH上升.由于水体叶绿素含量极低, 光合作用消耗的CO2低于浮游生物呼吸作用产生的CO2量; 因此河流有机碳的原位呼吸作用是维持河流过饱和p(CO2)的主要因素[8~11].竹溪河表层水p(CO2)与水温呈明显的负相关, 夏秋两季相关系数分别为-0.312和-0.301, 这与许多学者的研究结果是一致的[21~23].水温升高, CO2在水体中的溶解度降低, 另外水温升高是光照强度增强的结果, 最终导致p(CO2)降低.

表 2 p(CO2)与环境变量的Spearman相关分析1) Table 2 Spearman correlation analysis between p(CO2) and key environmental variables

4.2 与其他水体CO2脱气通量的比较

众多研究表明世界范围内的水体主要表现出溶解性CO2过饱和现象, 总体表现出CO2“源”的特征(表 3)[1~5].如亚马逊河、Tay River、长江上游龙川江、黄河、三峡库区干流和重庆城市河流的研究表明其CO2脱气通量分别为559.3、31.7、131.5、33.0、126.0和446.5 mmol·(m2·d)-1, 总体而言处于温带地区的河流释放量最低, 亚热带河流次之, 热带河流在三者中最高[28~35].

表 3 世界河流CO2通量 Table 3 CO2 area flux in world rivers

本研究结果表明竹溪河CO2释放量均值为(113.5±40.2) mmol·(m2·d)-1, 表现出向大气释放CO2; 这与以往研究认为水生生态系统是大气的排放源的结论一致[1~5].对比三峡库区香溪河库湾[16]和闽江口养殖塘[35]日内水-气界面CO2通量观测, 表明CO2日内释放通量分别在10.9~52.7 mmol·(m2·d)-1和-11.6~40.4 mmol·(m2·d)-1范围内波动; 可知近年来, 包括城市化在内的人类活动影响下的河流CO2释放较明显, 可能是人为输入的有机质为原位呼吸创造了条件, 整体表现出向大气释放大量的CO2.本研究中, 其水-气界面CO2通量呈现出上午释放量大于下午的特征, 这与中午后光照强度增加, 有利于水生生物的光合作用密切相关.本研究与王亮等[27]人在三峡水库香溪河库湾夏季水-气界面温室气体通量24 h的连续观测的结果一致.

4.3 嘉陵江支流竹溪河水-气界面CO2通量变化的影响因子分析

受到人类活动的高度干扰, 竹溪河表层水水化学参数、营养盐和溶解碳的特征在不同时间尺度下有着显著的差异, 一般水体向大气中吸收或释放CO2, 主要是其内部新陈代谢的结果, 水体中生化作用对水-气界面CO2交换通量起着重要影响.本研究通过F(CO2)与关键环境因子的Spearman相关性(表 4)分析, 探讨其主要环境驱动因子, 本研究表明F(CO2)与pH、水温和气温呈极显著负相关关系.

表 4 水-气界面CO2通量与环境变量的Spearman相关分析1) Table 4 Spearman correlation analysis between CO2 flux and key environmental variables

4.3.1 气象因子

水-气界面CO2通量在季节变化、年变化等长时间尺度上的变异性与温度存在很大关系[21]; 但在以日变化的短时间尺度上与温度的依存关系较为复杂, 主要表现为:正相关性[23]、负相关性[27]或无相关性[22].通常情况下, 水温通过影响气体分子在水中的溶解度、微生物活性、水生植物光合、呼吸作用和气泡的形成影响CO2的产生过程[21~23].

本研究表明, 水温和气温与水-气界面CO2脱气通量均为极显著负相关, 其相关性分别可用回归方程F(CO2)=-7.639X水温+288.147 (R2=0.152, P=0.06)和F(CO2)=-3.870X气温+206.753 (R2=0.152, P < 0.01)来表示.受太阳辐射等因素的影响, 气温是日变化中最明显的环境因子; 其日内波动较大.通过相关性分析表明, 夏、秋季水-气界面CO2通量与气温相关系数分别为-0.466和-0.420, 有明显的相关性, CO2通量对气温具有响应.气温变化可直接引起水温升高或降低, 相关性分析表明, 气温和水温呈明显的正相关, 相关系数为0.894, 这与许多学者的研究结果是一致的[21, 23].

监测期间, 随着水温升高, 二氧化碳在水体中的溶解度降低; 另外水温升高主要是由于光照强度增强导致的, 以及该区域水流比较缓, 适合营养盐聚集, 利于水生生物的光合作用, 使得大气中的CO2进入水体.

风速是影响水-气界面CO2通量的重要因素, 通过影响气体交换速率以及水-气接触面积等多种途径, 进而影响水中碳物质的分解和二氧化碳的产生[33].本研究中, 水-气界面CO2通量与风速呈正相关关系, 这与文献[33]的研究结果是一致的.监测期间, 竹溪河表层水体上方风速在0~2.2 m·s-1范围内波动, 小于3.0 m·s-1, 低风速环境下对气体扩散强度、水体气泡释放及浮游植物进行光合作用等扰动较小, 使得水-气界面CO2与风速没有明显的相关关系.

4.3.2 水化学参数

pH能直接影响水体中碳酸盐体系(CO2、CO32-和HCO3-), 和水体溶存CO2的浓度密切相关, 进而影响CO2的吸收或排放[21~22, 28].当pH较高时, 水中游离的CO2易形成碳酸盐, 水体p(CO2)降低, 从而减少CO2的产生, 有利于大气环境中CO2进入水体; 反之, 较低的pH值会促进水体CO2释放到大气中[12].本研究表明, 夏、秋季竹溪河水-气界面CO2通量与pH呈极显著负相关, 相关系数分别为-0.564和-0.900;其相关性分别可用F(CO2)=-232.762XpH+20 024.118 (R2=0.627, P <0.01)和F(CO2)=-674.848XpH+5 775.944 (R2=0.853, P < 0.01)来表示.进一步通过逐步回归分析, 发现pH是影响竹溪河水-气界面CO2通量的最主要的环境因子.造成这种现象可能是由于水体进行光合作用的浮游植物较低, 河流有机碳的原位呼吸占主导地位, 使得:

平衡向左移动, 促进水体CO2释放, 进而降低pH.这与Li等[12]在长江上游龙川江的研究结果类似.

水体溶解氧含量决定了水体中有机物降解的途径和产物, 对水体生源要素循环有重要影响[34].在有氧情况下, 有机物主要转换为CO2, 反之在缺氧环境下主要生成CH4[34].本研究结果表明, 秋季竹溪河表层水DO浓度与水-气界面CO2通量为不明显正相关(r=0.02);其相关性可用F(CO2)=15.104XDO+68.553(R2=0.037, P=0.382)来表示.这主要是水体中DO浓度升高促进水体中有机物快速有氧分解, 释放出CO2; 另外监测期间天气多为多云转晴, 太阳辐射波动较大, 水体中叶绿素含量较低, 表层水中浮游植物进行光合作用强度较弱, 进而弱化了溶解氧与河流CO2释放的响应.

水生植物初级生产力及水生生物新陈代谢受到水体营养盐浓度的影响, 进而间接影响水-气界面CO2通量[35].本研究中, 秋季竹溪河水-气界面CO2通量日变化在绝大多数情况下与TDN和TDP的日时间尺度上的变化存在负相关关系, 分别为-0.070和-0.020.原因在于水体中氮磷等营养元素可为水生植物及浮游藻类提供营养成分, 使其活性增强, 增加光合作用, 但同时营养盐增加也促进细菌的活性.但河流水体溶存性CO2浓度是光合作用、呼吸作用综合作用的结果; 由于在监测期间, 竹溪河表层水营养盐及有机质含量较低, 不利于CO2产生及溶存; 这种相关性与杨平等[35]对闽江口养殖塘研究结果一致.

水体中有机质主要来源与系统内生物残体以及外源输入有关, 其中DOC为水体中CO2产生的直接碳源, 对水体CO2的产生与排放具有重要影响[36].通常情况下, 单位面积水体的CO2通量与水体中DOC浓度表现出正相关关系[36].本研究中, 秋季竹溪河表层水体DOC浓度在3.30~4.55 mg·L-1范围内波动; 通过相关性F(CO2)与DOC没有明显的相关性, 可能是水体可生物降解的有机质含量较低, 从而弱化了其对水-气界面CO2气体释放的贡献.

利用多元回归中的逐步回归分析法建立夏、秋两季F(CO2)与各环境因子间的最优回归方程, 分析各环境因子对CO2通量影响的贡献大小:

综上所述, 表层水pH和碱度是影响竹溪河水-气界面CO2通量的主要环境因子, 它们之间的负相关系数均达到显著性水平(夏季CO2R2=0.996, P < 0.001;秋季CO2R2=0.993, P < 0.001).

5 结论

(1) 三峡库区嘉陵江支流竹溪河p(CO2)和水-气界面CO2通量呈现显著的日内、日间及季节变化; 夏季和秋季p(CO2)分别为(86.2±7.7)Pa和(93.4±8.9)Pa, CO2通量分别为(87.8±27.5) mmol·(m2·d)-1和(139.2±34.0) mmol·(m2·d)-1; 并表现出在上午09:00左右达最大值, 随后波动下降至下午18:00左右达到最小释放量的日内变化规律.

(2) 通过相关性分析, CO2通量主要受到pH、碱度和水温的影响; 以及通过逐步回归分析, 发现pH是影响CO2交换通量的最主要的环境因子.

参考文献
[1] 李哲, 白镭, 郭劲松, 等. 三峡水库两条支流水-气界面CO2、CH4通量比较初探[J]. 环境科学, 2013, 34(3): 1008-1016.
Li Z, Bai L, Guo J S, et al. Comparative study on water-air CO2, CH4 flux in two tributaries in the Three Gorges Reservoir, China[J]. Environmental Science, 2013, 34(3): 1008-1016.
[2] Knoll L B, Vanni M J, Renwick W H, et al. Temperate reservoirs are large carbon sinks and small CO2sources:Results from high-resolution carbon budgets[J]. Global Biogeochemical Cycles, 2013, 27(1): 52-64. DOI:10.1002/gbc.20020
[3] Erkkilå K M, Ojala A, Bastviken D, et al. Methane and carbon dioxide fluxes over a lake:comparison between eddy covariance, floating chambers and boundary layer method[J]. Biogeosciences, 2018, 15(2): 429-445. DOI:10.5194/bg-15-429-2018
[4] Müller D, Warneke T, Rixen T, et al. Lateral carbon fluxes and CO2 outgassing from a tropical peat-draining river[J]. Biogeosciences, 2015, 12(20): 5967-5979. DOI:10.5194/bg-12-5967-2015
[5] Huang W J, Cai W J, Wang Y C, et al. The carbon dioxide system on the Mississippi river-dominated continental shelf in the northern Gulf of Mexico:1. Distribution and air-sea CO2 flux[J]. Journal of Geophysical Research Oceans, 2015, 120(3): 1429-1445. DOI:10.1002/2014JC010498
[6] Cole J J, Prairie Y T, Caraco N F, et al. Plumbing the global carbon cycle:integrating inland waters into the terrestrial carbon budget[J]. Ecosystems, 2007, 10(1): 172-185. DOI:10.1007/s10021-006-9013-8
[7] Raymond P A, Hartmann J, Lauerwald R, et al. Global carbon dioxide emissions from inland waters[J]. Nature, 2013, 503(7476): 355-359. DOI:10.1038/nature12760
[8] Sawakuchi H O, Neu V, Ward N D, et al. Carbon dioxide emissions along the lower amazon river[J]. Frontiers in Marine Science, 2017(4): 76.
[9] Li S Y, Lu X X, Bush R T. CO2 partial pressure and CO2 emission in the lower Mekong river[J]. Journal of Hydrology, 2013, 540: 40-56.
[10] Oliveira A P, Cabeçadas G, Mateus M D. Inorganic carbon distribution and CO2 fluxes in a large European estuary (Tagus, Portugal)[J]. Scientific Reports, 2017, 7(1): 7376. DOI:10.1038/s41598-017-06758-z
[11] Almeida R M, Pacheco F S, Barros N, et al. Extreme floods increase CO2 outgassing from a large Amazonian river[J]. Limnology and Oceanography, 2017, 62(3): 989-999. DOI:10.1002/lno.10480
[12] Li S Y, Lu X X, He M, et al. Daily CO2 partial pressure and CO2 outgassing in the upper Yangtze River basin:a case study of the Longchuan River, China[J]. Journal of Hydrology, 2012, 466-467: 141-150. DOI:10.1016/j.jhydrol.2012.08.011
[13] Ran L, Lu X X, Richey J E, et al. Long-term spatial and temporal variation of CO2 partial pressure in the Yellow River, China[J]. Biogeosciences, 2015, 12(4): 921-932. DOI:10.5194/bg-12-921-2015
[14] 翟惟东.南海北部与珠江河口水域CO2通量及其调控因子[D].厦门: 厦门大学, 2003. 5-6.
Zhai W D. Air-sea fluxes of carbon dioxide and upper ocean biogeochemical processes in the Northern South China Sea and the Pearl river estuary[D]. Xiamen: Xiamen University, 2003. 5-6.
[15] 张陶, 李建鸿, 蒲俊兵, 等. 桂江流域夏季水-气界面CO2脱气的空间变化及其影响因素[J]. 环境科学, 2017, 38(7): 2773-2783.
Zhang T, Li J H, Pu J B, et al. Spatial variations of CO2 degassing across water-air interface and its impact factors in summer in Guijiang river, China[J]. Environmental Science, 2017, 38(7): 2773-2783. DOI:10.3969/j.issn.1000-6923.2017.07.044
[16] 黄文敏, 朱孔贤, 赵玮, 等. 香溪河秋季水-气界面温室气体通量日变化观测及影响因素分析[J]. 环境科学, 2013, 34(4): 1270-1276.
Huang W M, Zhu K X, Zhao W, et al. Diurnal changes in greenhouse gases at water-air interface of Xiangxi river in Autumn and their influencing factors[J]. Environmental Science, 2013, 34(4): 1270-1276.
[17] 国家环境保护总局. 水和废水监测分析方法[M]. (第四版).北京: 中国环境科学出版社, 2002: 239-284.
[18] UNESCO/IHAGHG. Greenhouse gas emissions related to freshwater reservoirs[R]. The World Bank, 2010. 64-127.
[19] Alin S R, de Fátima F L, Rasera M, Salimon C I, et al. Physical controls on carbon dioxide transfer velocity and flux in low-gradient river systems and implications for regional carbon budgets[J]. Journal of Geophysical Research, 2011, 116(G1): G01009.
[20] Neal C, Harrow M, Williams R J. Dissolved carbon dioxide and oxygen in the river Thames:Spring-summer 1997[J]. Science of the Total Environment, 1998, 210-211: 205-217. DOI:10.1016/S0048-9697(98)00013-8
[21] Stadmark J, Leonardson L. Emissions of greenhouse gases from ponds constructed for nitrogen removal[J]. Ecological Engineering, 2005, 25(5): 542-551. DOI:10.1016/j.ecoleng.2005.07.004
[22] Schrier-Uijl A P, Veraart A J, Leffelaar P A, et al. Release of CO2 and CH4 from lakes and drainage ditches in temperate wetlands[J]. Biogeochemistry, 2011, 102(1-3): 265-279. DOI:10.1007/s10533-010-9440-7
[23] Natchimuthu S, Selvam B P, Bastviken D. Influence of weather variables on methane and carbon dioxide flux from a shallow pond[J]. Biogeochemistry, 2014, 119(1-3): 403-413. DOI:10.1007/s10533-014-9976-z
[24] Richey J E, Melack J M, Aufdenkampe A K, et al. Outgassing from Amazonian rivers and wetlands as a large tropical source of atmospheric CO2[J]. Nature, 2002, 416(6881): 617-620. DOI:10.1038/416617a
[25] Hélie J F, Hillaire-Marcel C, Rondeau B. Seasonal changes in the sources and fluxes of dissolved inorganic carbon through the St. Lawrence River-Isotopic and chemical constraint[J]. Chemical Geology, 2002, 186(1-2): 117-138. DOI:10.1016/S0009-2541(01)00417-X
[26] 李双, 王雨春, 操满, 等. 三峡库区库中干流及支流水体夏季二氧化碳分压及扩散通量[J]. 环境科学, 2014, 35(3): 885-891.
Li S, Wang Y C, Cao M, et al. Partial pressure and diffusion flux of dissolved carbon dioxide in the mainstream and tributary of the central Three Gorges Reservoir in Summer[J]. Environmental science, 2014, 35(3): 885-891.
[27] 王亮, 肖尚斌, 刘德富, 等. 香溪河库湾夏季温室气体通量及影响因素分析[J]. 环境科学, 2012, 33(5): 1471-1475.
Wang L, Xiao S B, Liu D F, et al. Fluxes of greenhouse gases from Xiangxi River in summer and their influencing factors[J]. Environmental science, 2012, 33(5): 1471-1475.
[28] 蒋滔, 郭劲松, 李哲, 等. 三峡水库不同运行状态下支流澎溪河水-气界面温室气体通量特征初探[J]. 环境科学, 2012, 33(5): 1463-1470.
Jiang T, Guo J S, Li Z, et al. Air-water surface greenhouse gas flux in Pengxi River at different operational stages of the Three Gorges Reservoir[J]. Enviromental science, 2012, 33(5): 1463-1470.
[29] Huang W M, Bi Y H, Hu Z Y, et al. Spatio-temporal variations of GHG emissions from surface water of Xiangxi River in Three Gorges Reservoir region, China[J]. Ecological Engineering, 2015, 83: 28-32. DOI:10.1016/j.ecoleng.2015.04.088
[30] Zhao Y, Wu B F, Zeng Y. Spatial and temporal patterns of greenhouse gas emissions from Three Gorges Reservoir of China[J]. Biogeosciences, 2013, 10(2): 1219-1230. DOI:10.5194/bg-10-1219-2013
[31] Wang X F, He Y X, Yuan X Z, et al. pCO2 and CO2 fluxes of the metropolitan river network in relation to the urbanization of Chongqing, China[J]. Journal of Geophysical Research Biogeosciences, 2017, 122(3): 470-486. DOI:10.1002/2016JG003494
[32] Li S L, Chetelat B, Yue F J, et al. Chemical weathering processes in the Yalong River draining the eastern Tibetan Plateau, China[J]. Journal of Asian Earth Sciences, 2014, 88: 74-84. DOI:10.1016/j.jseaes.2014.03.011
[33] Guérin F, Abril G, Serça D, et al. Gas transfer velocities of CO2 and CH4 in a tropical reservoir and its river downstream[J]. Journal of Marine Systems, 2007, 66(1-4): 161-172. DOI:10.1016/j.jmarsys.2006.03.019
[34] 程炳红, 郝庆菊, 江长胜. 水库温室气体排放及其影响因素研究进展[J]. 湿地科学, 2012, 10(1): 121-128.
Cheng B H, Hao Q J, Jiang C S. Research progress on the emission of greenhouse gases from reservoir and its influence factors[J]. Wetland Science, 2012, 10(1): 121-128. DOI:10.3969/j.issn.1672-5948.2012.01.017
[35] 杨平, 仝川, 何清华, 等. 闽江口鱼虾混养塘水-气界面温室气体通量及主要影响因子[J]. 环境科学学报, 2013, 33(5): 1493-1503.
Yang P, Tong C, He Q H, et al. Greenhouse gases fluxes at water-air interface of aquaculture ponds and influencing factors in the Min river estuary[J]. Acta Scientiae Circumstantiae, 2013, 33(5): 1493-1503.
[36] Jacinthe P A, Filippelli G M, Tedesco L P, et al. Carbon storage and greenhouse gases emission from a fluvial reservoir in an agricultural landscape[J]. Catena, 2012, 94: 53-63. DOI:10.1016/j.catena.2011.03.012