湖泊科学   2024, Vol. 36 Issue (2): 499-511.  DOI: 10.18307/2024.0227
0

研究论文——生物地球化学与水环境保护

引用本文 [复制中英文]

刘纪明, 赵优, 肖启涛, 屈琦, 胡又心, 李郑杰, 易心钰, 卢伟志, 亚热带城市湖泊与河流CO2气体通量特征及其影响因素. 湖泊科学, 2024, 36(2): 499-511. DOI: 10.18307/2024.0227
[复制中文]
Liu Jiming, Zhao You, Xiao Qitao, Qu Qi, Hu Youxin, Li Zhengjie, Yi Xinyu, Lu Weizhi. Characteristics and influencing factors of CO2 flux in subtropical urban lakes and river. Journal of Lake Sciences, 2024, 36(2): 499-511. DOI: 10.18307/2024.0227
[复制英文]

基金项目

湖南省自然科学基金项目(2021JJ20012, 2021JJ41071)、湖南省教育厅项目(20B618)、湖南省水利科技项目(XSKJ2022068-35)、林业科技攻关与创新(XLKY202216)和湖南林业科技创新杰青培养项目(XLK202108-05)联合资助

通信作者

卢伟志, E-mail: weizhilu@csuft.edu.cn

文章历史

2023-07-03 收稿
2023-08-16 收修改稿

码上扫一扫

亚热带城市湖泊与河流CO2气体通量特征及其影响因素
刘纪明1,2 , 赵优1,2 , 肖启涛3 , 屈琦1,2 , 胡又心1,2 , 李郑杰1,2 , 易心钰4,5 , 卢伟志1,2     
(1: 中南林业科技大学生命科学与技术学院, 长沙 410004)
(2: 中南林业科技大学南方林业生态应用技术国家工程实验室, 长沙 410004)
(3: 中国科学院南京地理与湖泊研究所, 中国科学院流域地理学重点实验室, 南京 210008)
(4: 湖南省植物园, 长沙 410000)
(5: 湖南长株潭城市群森林生态系统国家定位观测研究站, 长沙 410116)
摘要:湖泊、河流等内陆水体是连接陆地生态系统和海洋的“长程碳环路”的重要节点,也是温室气体二氧化碳(CO2)排放源,在调节陆地、海洋间的碳迁移转换中发挥着重要作用。相对于自然水体,城市水体因面积小、水深浅且受监测方法限制,水-气界面碳通量经常被忽略。为探讨我国亚热带城市水体温室气体排放特征,本研究以湖南省长沙市典型城市水体,包括洋湖、西湖、松雅湖、月湖4个湖泊和湘江长沙段为研究对象,分别于2022年4和10月采用光化学反馈-腔增强吸收光谱法(OF-CEAS)和扩散模型法对水-气界面CO2通量进行对比测定。结果表明,长沙城市湖泊与河流春季为CO2排放源,秋季为吸收汇,河流水-气界面CO2通量呈显著季节差异。河湖之间CO2通量在春季表现为显著差异,秋季差异不显著。CO2通量与水体溶解氧、水体总氮浓度等呈显著正相关。2种方法的CO2通量对比测定在湖泊上显著相关,但对河流而言相关性不显著。研究揭示的城市湖泊与河流CO2气体的排放特征有利于深入探究城市水体碳的迁移转化,可对全面了解全球气候变化过程和河湖湿地温室气体减排和调控提供科学支撑。
关键词城市湖泊与河流    CO2通量    季节变化    影响因素    
Characteristics and influencing factors of CO2 flux in subtropical urban lakes and river
Liu Jiming1,2 , Zhao You1,2 , Xiao Qitao3 , Qu Qi1,2 , Hu Youxin1,2 , Li Zhengjie1,2 , Yi Xinyu4,5 , Lu Weizhi1,2     
(1: College of Life Science and Technology, Central South University of Forestry and Technology, Changsha 410004, P.R. China)
(2: National Engineering Laboratory for Applied Technology of Forestry & Ecology in South China, Central South University of Forestry and Technology, Changsha 410004, P.R. China)
(3: Key Laboratory of Watershed Geographic Sciences, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing 210008, P.R. China)
(4: Hunan Botanical Garden, Changsha 410000, P.R. China)
(5: Hunan Changsha-Zhuzhou-Xiangtan City Cluster National Research Station of Forest Ecosystem, Changsha 410116, P.R. China)
Abstract: Inland water bodies, such as lakes and rivers, are important nodes of the "long-range carbon loop" connecting terrestrial ecosystems and oceans, as well as sources of greenhouse gas carbon dioxide (CO2) emissions. Inland water bodies also play an important role in regulating carbon migration between land and sea. Compared with natural water bodies, carbon fluxes at the water-air interface are often overlooked in urban water due to their small sizes, shallow depths, and limitations in monitoring methods. In order to explore the greenhouse gas emission characteristics of such water bodies in China, this study took four typical urban water bodies in Changsha, Hunan Province, including Lake Yang, Lake West, Lake Songya and Lake Yue, and Changsha section of Xiangjiang River as research objects. The CO2 fluxes at the water-gas interface were measured by laser-based gas analyzer that used optical feedback-cavity-enhanced absorption spectroscopy (OF-CEAS) and diffusion model in April and October 2022, respectively. Results showed that lakes and rivers in Changsha city were CO2 sources in spring but CO2 sinks in fall, and there was a significant seasonal difference in CO2 fluxes at the river water-air interface. Moreover, the CO2 flux was positively correlated with the content of dissolved oxygen and total nitrogen. The comparative determination of CO2 fluxes between the two methods was significantly correlated on lakes, but not on rivers. The study reveals that the CO2 emission characteristics of urban lakes and rivers are beneficial to the in-depth investigation of carbon transport and transformation, which can provide scientific support for the comprehensive understanding of global climate change processes and the reduction and regulation of greenhouse gas emissions in inland wetlands.
Keywords: Urban lakes and river    CO2 flux    seasonal variability    influence factors    

河湖等内陆水体的二氧化碳(CO2)排放与吸收是全球碳循环的重要过程[1],其向大气排放的CO2量可相当于全球海洋和陆地的汇[2],而现如今内陆水体CO2通量的研究对象多集中在湖泊、水库、河流及河口[3-6]等大型自然水域,城市水体研究较少。已有研究表明内陆水体CO2排放高达2.1 Pg C/a[7],是陆地生态系统温室效应的重要贡献者[8],超过90% 来自湖泊与河流[1],其中,城市湖泊CO2排放高于非城市湖泊,城市河流的CO2排放是非城市河流的2倍[9-10],城市水库CO2排放量是非城市的3倍[11],可见城市水体CO2排放在区域和全球碳循环过程中至关重要。

相对于自然水体而言,城市水体多为缓流或静止水体、水深相对较浅,加上人类活动向城市水体输入大量营养盐和有机物导致其自净能力差[12-13]。一方面水体营养物质的大量输入会进一步刺激微生物活性和新陈代谢,导致水体中产生大量的CO2[14-16],所以城市水体常被认为是CO2的排放源;另一方面水体营养物质增加也有利于浮游植物的生长,通过光合作用,消耗和吸收CO2,减少了CO2排放[10, 17]。因此城市水体的碳源汇角色及其驱动机制需要更一步探究。

目前水-气界面CO2通量的测定包括扩散模型法、通量箱法、倒置漏斗法以及微气象法[18],其中通量箱法和扩散模型法由于价格低廉、操作简便、测定区域广阔,目前应用最广泛[17-19],且被认为估算结果更接近实际排放强度[20]。但扩散模型法受风速制约,研究结果通常变异较大且需进一步优化模型,而通量箱法每个点需要花费大量的时间采集气体,不便于进行大面积的温室气体监测,且受到许多因素的影响,例如箱体扰动、箱体内外温差、箱内气压变化、箱体大小以及箱内气体混合程度等[21-22]。已有对城市水体研究发现扩散模型法计算出来的CO2通量明显高于通量箱法[19, 23],而也有研究得出不同的结果[24]。在此背景下,本研究对通量箱法进行技术改进,耦合光化学反馈-腔增强吸收光谱技术形成新的通量箱测定方法OF-CEAS,该方法相较于传统通量箱法具有便携、现场监测、测量速度快、测量频次高等优点,形成了与扩散模型法可比较的空间覆盖的特点。

长沙作为长江中游地区重要中心城市,是典型的亚热带季风区城市之一,据2022年长沙统计年鉴显示,2021年长沙市降雨量1390.8 mm,年平均气温18.9℃,具有降水充沛和雨热同期的特点。长沙市河湖水系完整,河网密布,市区水域面积占比高达18%,同时湘江长沙段和长沙市众多湖泊具有典型亚热带城市水体的特点,水体受到农业源、工业源、生活源等污染,对温室气体总量的估算影响不可小视。因此本研究以长沙城市水体为例,采用OF-CEAS法和扩散模型法对城市湖泊与河流CO2通量进行对比研究,探究不同季节和2种方法间CO2通量差异,揭示CO2通量的时空变化特征及其控制要素。研究结果可揭示亚热带城市水体湖泊与河流的CO2产生和城市水体CO2通量的时空变化特征及其与水环境间的相互关系。

1 材料与方法 1.1 采样点布设

于2022年4和10月选择晴天对长沙城市湖泊和湘江长沙段水体CO2通量进行测定,每次采样和测定持续时间为3~4 d,采样点设置见图 1。通过权衡测定安全距离与样点的代表性,本研究所有采样点均离岸距离3 m以上。其中湖泊包括西湖(28°13′N, 112°56′E)、洋湖(28°08′N, 112°55′E)、月湖(28°14′N, 113°02′E)和松雅湖(28°16′N, 113°06′E),湘江长沙段(28°01′~28°18′N, 112°55′~112°57′E)穿长沙城而过,西湖和洋湖位于湘江西侧,月湖和松雅湖位于湘江东侧。本研究的湖泊面积范围在0.6~3 km2,水深在0.5~4 m之间,湖泊中有浮水、沉水、挺水和藻类等植物分布。每个湖泊各选择5个样点,湘江长沙段6个样点,共计26个采样点,每个样点采样3次作为重复,CO2通量测定与采样在9:00—16:00之间进行。

图 1 长沙市湖泊和河流采样点分布 Fig.1 Distribution of sampling sites of lakes and river in Changsha
1.2 CO2通量测定

CO2通量测定采用光化学反馈-腔增强吸收光谱法(OF-CEAS)和扩散模型法2种方法。在每个采样点现场通过OF-CEAS法进行直接测定获取水-气界面CO2通量(FCO2_OC)。同时使用采水器采集表层10 cm深度的湖水3 L,用60 mL注射器采取水样通过扩散模型法测定扩散通量(FCO2_DM)。

1.2.1 光化学反馈-腔增强吸收光谱法OF-CEAS

该方法是利用智能腔室连接CH4/CO2/H2O微量气体分析仪(美国/LI-COR,LI-7810),同时搭配漂浮板测量水体CO2通量FCO2_OC(附图Ⅰa),智能腔室体积为4244.10 cm3,测量时内无漂浮植物。采用光化学反馈-腔增强吸收光谱技术来测量空气中CO2浓度的变化,智能测量室可以通过连接微量气体分析仪测量CO2气体通量,CO2通量计算公式如下:

$ F_{\mathrm{c}}=\frac{10 V P_0\left(1-\frac{W_0}{1000}\right)}{\operatorname{RS}\left(T_0+273.15\right)} \frac{a c^{\prime}}{a t} $ (1)

式中,Fc是气体排放通量FCO2_OC(μmol/(m2·s)),V是体积(cm3),R是气体常数(8.314 Pa/(m3·K·mol)),P0是初始压强(kPa),W0是初始水蒸气摩尔分数(mmol/mol),S是水体表面积(cm2),T0是初始空气温度(℃);$\frac{a c^{\prime}}{a t}$是水分校正后的CO2摩尔分数初始变化速率(μmol/(mol·s)),整个计算过程可通过SoilFluxProTM(美国/LI-COR,version 5.2.0)软件完成。

1.2.2 扩散模型法

扩散模型法所测得的扩散通量FCO2_DM是通过表层水体和大气气体浓度差计算水-气界面交换速率而得到[25],每个采样点采集3个水体作为平行样品,水样用60 mL带有三通阀的医用注射器采集,每次采样前润洗3次注射器和采水器,在水中匀速缓慢抽取30 mL水样至注射器。本研究采用水-气界面CO2理论扩散模型计算扩散排放通量[25-27]

$ F=k\left(C_{\mathrm{w}}-C_{\mathrm{eq}}\right) $ (2)

式中,F为扩散排放通量FCO2_DM(μmol/(m2·s)),k为气体交换系数(cm/h),定义为风速和气体施密特数的函数[28]Ceq为大气中CO2浓度(μmol/L),Cw为CO2在表层水体中浓度(μmol/L),根据水-气平衡后采样管顶部CO2浓度和Bunsen系数计算。顶部CO2浓度通过水-气顶空平衡-气相色谱法测定,即向装有水样的注射器注入高纯氮气使注射器上方形成气室,在25℃条件下将注射器水平剧烈振荡20 min,然后静置2 h达到气-液两相平衡[29-30],将上层平衡气体注入密封真空瓶中(12 mL,英国),采用气相色谱仪(美国,安捷伦,Agilent 7890B)于24 h内测定完成。风速采用现场瞬时风速,Bunsen系数使用水温和盐度进行计算矫正,该方法测试示意图参照高洁[31]附图Ⅰb

1.3 水环境参数测定

现场原位测定环境参数包括水体温度(TW)、酸碱度(pH)、电导率(EC)、总固体溶解度(TDS)、电阻率(RT)、盐度(PSU)、氧化还原电位(ORP)和溶解氧(DO)等通过便携式多参数水质测定仪(意大利,HANNA,HI98194-04)测定;蓝绿藻-藻蓝蛋白(PC)、叶绿素a(Chl.a)、铵态氮(NH4+-N)通过YSI(美国, 赛莱默,EXO2多参数水质分析)测定;风速使用便携式风速仪测定。同时采集表层水样带回实验室后立即进行总有机碳(TOC)、总无机碳(TIC)、溶解性有机碳(DOC)、溶解性无机碳(DIC)、总氮(TN)和总磷(TP)浓度的测定。其中TOC、TIC、DOC、DIC浓度采用总有机碳分析仪(日本,岛津,TOC-V WP)测定;TN浓度采用过硫酸钾氧化-紫外分光光度法测定;TP浓度采用过硫酸钾氧化-钼蓝比色法测定。

1.4 数据分析

实验数据计算、处理与统计通过Micorsoft Excel 2010和IBM SPSS Statistics 27进行,采用相关性分析方法(Pearson correlations)和单因素方差分析(One-way ANOVA)方法分析数据之间的相关性和差异显著性,且符合检验要求;采用Rstudio(R4.2.1)软件“ggplot2”等包进行CO2通量数据汇总、绘制相关图件;采用ArcGIS 10.4进行采样点位分布图绘制。

2 结果分析 2.1 城市湖泊与河流水环境参数特征分析

湖泊与河流水质特征如表 1所示,4个城市湖泊和河流水体pH在7.6~8.9之间,均呈弱碱性。河湖中的DO浓度存在显著季节差异,且河流与除西湖外3个湖泊之间存在显著差异(P < 0.05),洋湖的DO浓度在秋季显著大于其他3个湖泊与河流。春季洋湖、西湖和湘江TP浓度显著高于松雅湖和月湖,秋季各湖泊与河流之间TP浓度无显著差异。洋湖的TN浓度在春、秋两季都显著高于其他3个湖泊和河流(P < 0.05)。所有湖泊的秋季TOC浓度均大于春季,其中西湖春秋两季之间存在显著差异(P < 0.05),且均高于其他湖泊和河流,西湖TOC浓度在春季显著高于松雅湖、月湖和河流(P < 0.05),在秋季,其浓度显著高于其他3个湖泊与河流(P < 0.05);河流秋季TOC浓度小于春季。湖泊和河流的TIC浓度范围在4.0~25.1 mg/L,均为秋季低于春季,洋湖、西湖以及河流的TIC在春秋季间存在显著差异(P < 0.05)。

表 1 4个城市湖泊和河流部分水质参数* Tab. 1 Water quality parameters of four urban lakes and river

为了进一步探讨CO2通量的影响因素,在秋季增测了PC、Chl.a、NH4+-N等指标。秋季湖泊PC、Chl.a、NH4+-N浓度高于河流(图 2),西湖的PC与Chl.a浓度最高,显著大于其他湖泊与河流(P < 0.05)(图 2ab)。湖泊整体的NH4+-N浓度显著大于河流(P < 0.05),其中洋湖的NH4+-N浓度最高且显著高于其他河湖(P < 0.05),西湖与月湖、湘江存在显著差异(P < 0.05),松雅湖与湘江存在显著差异(P < 0.05)(图 2c)。

图 2 4个城市湖泊和河流水体秋季PC(a)、Chl.a(b)和NH4+-N(c)浓度特征 (不同小写字母表示方差分析具有显著性差异,P < 0.05) Fig.2 Concentration characteristics of PC (a), Chl.a (b) and NH4+-N (c) in four urban lakes and river in autumn
2.2 城市湖泊与河流CO2通量特征

OF-CEAS法测定湖泊水体CO2通量春季均值为(0.10±0.48) μmol/(m2·s) (n=20),秋季为(-0.16±0.61) μmol/(m2·s) (n=20);河流水体CO2通量春季均值为(1.32±0.65) μmol/(m2·s) (n=5),秋季为(-0.23±0.10) μmol/(m2·s) (n=5)。湖泊的CO2通量表现为春季排放,秋季吸收;河流的CO2通量季节变化显著(P < 0.01),春季排放,秋季吸收(图 3a)。在春季,河湖之间CO2通量差异显著(P < 0.05),洋湖为排放源,其他3个湖泊为汇,且洋湖CO2与其他3个湖泊之间存在显著差异(P < 0.05),河流CO2通量与西湖、松雅湖、月湖之间存在显著差异(P < 0.05)(图 3b)。在秋季,河湖整体表现为吸收,无显著差异,但其中洋湖为排放源,其他3个湖泊为碳汇。洋湖CO2通量与西湖、月湖之间存在显著差异(P < 0.05),松雅湖CO2通量与月湖之间存在显著差异(P < 0.05)(图 3c)。对于湖泊来说秋季CO2通量低于春季都呈现降低的趋势。从季节波动上来看,河流FCO2_OC波动大于湖泊。

图 3 城市湖泊与河流CO2通量季节变化 (*表示方差分析检验结果P < 0.05, **表示P < 0.01) Fig.3 Seasonal CO2 fluxes of urban lakes and river

扩散模型法监测到湖泊水体CO2通量春季为(0.02±0.04) μmol/(m2·s) (n=20),秋季为(0.01±0.03) μmol/(m2·s) (n=20);河流水体CO2通量春季为(0.08±0.04) μmol/(m2·s) (n=5),秋季为(-0.0003±0.004) μmol/(m2·s) (n=5)。湖泊的CO2通量季节变化不显著,春秋两季整体表现为弱排放;河流的CO2通量季节变化显著(P < 0.05),春季弱排放,秋季弱吸收(图 3d)。在春季,河湖之间CO2通量差异显著(P < 0.05),洋湖CO2通量排放大于其它3个湖泊,湘江与西湖、松雅湖、月湖之间存在显著差异(P < 0.05)(图 3e)。在秋季,河湖之间CO2通量差异不显著,湖泊为弱源,河流为弱汇(图 3f)。洋湖、松雅湖为碳源,西湖、月湖为碳汇,且洋湖与西湖和月湖之间存在显著差异(P < 0.05)。对于4个湖泊来说秋季CO2通量都呈现降低的趋势。从季节波动上来看,河流FCO2_DM波动大于湖泊。

2.3 城市湖泊与河流CO2通量控制要素

2种方法测得的CO2通量与水环境参数相关性分析表明,CO2通量和水质因子的相关性表现出一致性。湖泊FCO2_OC和FCO2_DM,与EC、TDS、PSU、TIC、TN存在极显著正相关,与DO显著正相关(图 4a),与RT、pH存在极显著负相关,与TW显著负相关(图 4a);河流FCO2_OC和FCO2_DM,与TOC、DO、RT、TIC、TN存在极显著正相关(图 4b),与EC、TDS、PSU、TW存在极显著负相关(图 4b)。湖泊的CO2通量与NH4+-N、DOC、DIC存在正相关关系,其中与NH4+-N呈显著正相关,与PC、Chl.a存在负相关关系(表 2)。对河流来说FCO2_OC与PC、DOC、DIC存负相关关系,其中与DIC显著负相关,与Chl.a、NH4+-N呈正相关,河流FCO2_DM与PC存在显著正相关,与DOC存在显著负相关(表 2)。

图 4 春秋两季湖泊(a)和河流(b)通量与环境参数相关关系 (*表示方差分析检验结果P < 0.05, **表示P < 0.01, ***表示P < 0.001) Fig.4 Correlation significance analysis of parameters and fluxes of lakes (a) and river (b) in spring and autumn
表 2 秋季城市湖泊与河流CO2通量与部分水环境参数的Pearson相关性 Tab. 2 Pearson correlation between CO2 flux and water environmental parameters in urban lakes and river in autumn
2.4 OF-CEAS法与扩散模型法测定结果比较

OF-CEAS法和扩散模型法测定结果对比表明,2种方法所得CO2通量差异性显著。春、秋两季湖泊的FCO2_OC和FCO2_DM相关性极显著(春季R2=0.48, P < 0.001;秋季R2=0.77, P < 0.001,图 5a),且FCO2_OC大于FCO2_DM。河流的FCO2_OC和FCO2_DM两者之间的相关性不显著(P>0.05)(图 5b),且FCO2_OC大于FCO2_DM。2种监测方法在CO2通量结果估算上存在较大差异,无论是春季还是秋季,湖泊与河流的FCO2_DM均小于FCO2_OC。研究结果表明尽管OF-CEAS法和扩散模型法测定结果存在量级上的差异,但2种方法在揭示湖泊与河流源汇问题上有较好的一致性。

图 5 湖泊(a)与河流(b)OF-CEAS法和扩散模型法测定结果回归分析 Fig.5 Regression analysis of CO2 emission for OF-CEAS and diffusion model observations lakes(a) and river(b)
3 讨论 3.1 城市湖泊与河流CO2通量时空动态与驱动要素

本研究结果表明,长沙城市湖泊CO2通量交换较弱且时空特征复杂,整体呈现源汇交替状态。同时长沙作为亚热带城市,降水、气温等因素都可能会对通量排放或吸收造成影响。其中,洋湖在春、秋两季呈现较高的FCO2_OC排放,这可能与洋湖水体TIC、TN、NH4+-N浓度高于其他湖泊(表 1图 2c),富营养化程度较高有关(图 6)。其他3个湖泊FCO2_OC在春、秋两季都表现为吸收,这与以往研究表现一致[32]。城市河流CO2通量呈现出较明显的季节动态,春季河流CO2通量表现为释放,秋季表现为吸收状态,与南京市河流[33]、上海河流[23]在秋季呈现CO2通量为排放结果相反,这可能与人类活动强度大小、水体理化性质以及河道类型宽度等条件相关。长三角地区城市化和人类活动强度比长沙更高,内、外秦淮河和金川河全部为闸控河流,河道两边均为人筑石堤,河道宽度较窄;淀浦河、苏州河呈现出缺氧的特征,河道宽度分别为120和45 m,这些河流的宽度远小于湘江长沙段。同时水体周围用地环境也会显著影响气体通量,如湘江长沙段周边分布的居住小区、工业区、农作物区等;湖泊周边分布的居住小区、学校、医院、建筑施工的用地和待开发建设用地等。

城市河湖水-气界面CO2通量与水体碳(TOC、DOC、TIC、DIC)、氮(TN、NH4+-N)、磷(TP)营养盐及水质指标DO、pH相关性较高。本研究结果表明,长沙城市湖泊与河流水体均为富营养化状态(表 1)。营养水平高的湖泊往往会接纳更多的氮、磷和外源碳等物质[34]。一方面,水体氮、磷营养盐通过影响浮游植物初级生产力和水生生物新陈代谢过程间接影响水-气界面CO2排放通量变化[35-37]。水体氮、磷营养盐浓度增加可导致浮游生物呼吸作用释放的CO2增多,引起水中CO2浓度上升,促使水中CO2进入大气[38-39]。本研究中城市湖泊与河流CO2通量与TN浓度呈显著相关(图 6),湖泊中NH4+-N与CO2通量呈显著正相关(表 2),这些特征在洋湖表现得最为显著。另一方面,湖泊富营养化对碳通量的影响受内源藻类光合作用的影响,富营养化导致藻类增加,固碳能力加强,降低CO2排放[40],但城市频繁的人类活动导致输入大量的外源有机质,为CO2产生提供碳源,最终增加CO2的产生和排放[41-43]。对于河流而言,30% ~90% 的CO2可能源于陆域土壤无机碳的直接输入[44-45],且春季河流中TOC、TN、TIC浓度高于秋季(表 1),河流CO2通量与TOC、TN、TIC浓度呈显著正相关(图 4图 6),春季河流CO2通量表现为排放,秋季表现为吸收状态,这也说明长沙城市河流CO2通量与城市工业产业发展情况、碳氮等有机质浓度有关。

图 6 DO、TIC、TN与FCO2_OC及FCO2_DM相关性分析 Fig.6 Correlation analysis of DO, TIC, TN with FCO2_OC and FCO2_DM diffusion flux

研究还表明长沙湖泊与河流水体均呈碱性(表 1),湖泊水体CO2通量与pH呈显著负相关,pH可通过改变水体中碳酸盐的动态平衡与形态直接影响CO2的排放,当pH值较高时,碳在水中主要以碳酸盐形式存在,水中CO2浓度减小,进而CO2排放减小。DO在一定程度上反映了水体光合作用和呼吸作用的强度,一方面水体中DO浓度越高,表明水体中水生和浮游植物的光合作用越强,溶存的CO2越少,CO2通量越小,当其通量为负值时表现为对CO2的吸收就越强;另一方面,随着DO浓度的升高,水生生物的呼吸作用加强,产生和排放的CO2就增加。本研究中,城市水体CO2通量与DO浓度呈显著正相关(图 4图 6),除西湖外其他3个湖泊的DO浓度在秋季增加,且对CO2通量吸收增强,推测在长沙城市湖泊中水生和浮游植物的光合作用可能占主导。而河流春季排放CO2,秋季吸收CO2,且春季水温偏低,DO浓度大于秋季,推测在春季河流可能主要以水生生物呼吸作用占主导,秋季以水生植物光合作用占主导。

3.2 水体CO2通量测定方法对比分析

本研究改进传统的静态箱法,构建基于光化学反馈-腔增强吸收光谱静态箱法与扩散模型法共同监测CO2通量。OF-CEAS法在传统的通量箱基础上进行改进,主要包括:(1)风箱控制的闭合机制,保证气室的放置对每个测量样点的完整性;(2)气室通风口能保证在有风的环境下也能使腔室内外的压强保持平衡;(3)安装的水体环使得水面大小一致确保水体的可比较性。OF-CEAS法通过改进设计允许在多个地点进行快速测量,也可用于对特定应用或实验类型进行长期测量,实现单个现场的多次重复。研究中OF-CEAS法因其腔室设计测定的是水-气界面总通量,而扩散模型法仅估算气体扩散通量,未考虑其他因素,例如植物分布和传输的影响,故OF-CEAS法估算通量高于扩散模型法,这与之前研究结果相似[46-47]。Duchemin等[20]在加拿大魁北克北部2个水库的研究中发现,漂浮箱法测定的CO2通量比扩散模型法高1.6~4倍;Vachon等[47]认为箱法造成的通量高估会随着风速的增加而降低,当风速大于4 m/s时,两者估算结果趋于一致。本研究区位于城市,城市周边的建筑物导致风速大小有很大的不确定性,加之浮游植物快速繁殖导致水体碳酸盐平衡体系向着碱性方向转变[48],因此基于碳酸盐平衡体系和风速模型的扩散模型估算的结果偏低。

本研究中湖泊上CO2通量2种方法的相关性显著强于河流水体(图 5)。对于湖泊来说,由于面积较小,水体属于静态水体受风浪影响较小,所以2种方法研究结果存在一致性,表明2种方法在城市湖泊CO2气体排放中具有较好的适用性;对于湘江这种具有通航功能的河流来说,水体属于缓流动态水体,加之受到风浪或航运的影响极大地改变了河道的水动力学,从而影响CO2通量结果,2种方法研究结果一致性不高,当然其他因素譬如水生植物分布等也是影响两者不一致的重要原因。

3.3 国内外湖泊与河流水-气界面CO2通量对比

通过对长沙城市湖泊与河流CO2通量的研究发现,城市水体不同湖泊与河流之间存在较大的空间变异性。由表 3可知,不同的气候地区和不同的监测方法都会导致CO2通量估算值存在差异。2种方法测得的长沙城市湖泊与河流CO2通量比同为亚热带地区的巢湖流域的河流排放低[9],略小于属于热带草原气候的赞比西盆地的河流,高于寒带南极洲的湖泊。OF-CEAS法测得的CO2通量与亚热带上海淀浦河和苏州河相近,与温带的内蒙古,东北地区的湖泊相近。扩散模型法测得的CO2通量值与亚热带阿根廷中东部湖泊、武汉东湖、加拿大Wascana Lake、西班牙南部Guadalcacín Lake,以及温带长春的南湖、北湖、地理所内湖、美国威斯康星州Mendota Lake值相近。以往对不同类型的湖泊与河流水面CO2通量研究表明,植被覆盖、气候和人类活动对城市地区的CO2通量有显著的影响[49]。同时湖泊与河流的海拔、面积大小、水深、水表湍流和富营养化等影响也不容忽视,对东北地区95个湖泊的研究表明城市湖泊CO2通量高于非城市地区且CO2的平均排放量随着湖泊面积的增加而下降。在巢湖流域,城市河流CO2排放大于非城市地区的河流。由表 3可知,整体来说富营养化程度越高会导致水体CO2排放增加,但同时也有富营养化水体对CO2表现出吸收的状态,这可能与水中水生生物类型和初级生产强度等水体环境差异有关。因此,对亚热带乃至全球内陆水体温室气体的评估不仅要考虑河湖等水体本身所处区域的不同以及面临的环境问题,还需要考虑监测方法的准确性、普适性、监测频次的统一以及更大更多区域尺度的研究。

表 3 不同湖泊与河流水-气界面CO2通量值 Tab. 3 CO2 flux at the water air interface of different lakes and rivers
4 结论

对长沙市城市湖泊与河流水-气界面CO2通量研究表明,长沙城市湖泊与河流之间CO2通量春季差异显著(P < 0.05),且城市河流的水-气界面CO2通量呈显著季节差异(P < 0.05)。整体来说,河流与湖泊春季为排放源(FCO2_OC=(0.34±0.72) μmol/(m2·s)) (n=25),秋季为吸收汇(FCO2_OC=(-0.17±0.55) μmol/(m2·s))(n=25)。CO2通量与水体溶解氧、水体总氮、总无机碳浓度等呈显著正相关(P < 0.05)。2种方法的对比研究在湖泊具有显著相关(P < 0.05),但在河流相关并不显著,表明未来需要开展不同观测方法的比较以及不确定性来源解析,研究结果对未来内陆水体的碳排放研究具有重要意义。

5 附录

附图Ⅰ见电子版(DOI: 10.18307/2024.0227)。

附图Ⅰ 水-气界面CO2通量测定方法示意:(a)光化学反馈-腔增强吸收光谱法(OF-CEAS);(b)扩散模型法测定表层水中CO2浓度操作示意图 AttachedFigI Determination of CO2 flux at the water-air interface: (a) OF-CEAS method device schematic diagram; (b) Operation diagram for determination of CO2 concentration in surface water
6 参考文献

[1]
Ran LS, Butman DE, Battin TJ et al. Substantial decrease in CO2 emissions from Chinese inland waters due to global change. Nature Communications, 2021, 12: 1730. DOI:10.1038/s41467-021-21926-6
[2]
Borges AV, Darchambeau F, Teodoru CR et al. Globally significant greenhouse-gas emissions from African inland waters. Nature Geoscience, 2015, 8(8): 637-642. DOI:10.1038/ngeo2486
[3]
Baliña S, Sánchez ML, Izaguirre I et al. Shallow Lakes under alternative states differ in the dominant greenhouse gas emission pathways. Limnology and Oceanography, 2023, 68(1): 1-13. DOI:10.1002/lno.12243
[4]
Deng HG, Zhang ZB, Liu T et al. Dissolved greenhouse gas concentrations and the influencing factors in different vegetation zones of an urban lake. J Lake Sci, 2019, 31(4): 1055-1063. [邓焕广, 张智博, 刘涛等. 城市湖泊不同水生植被区水体温室气体溶存浓度及其影响因素. 湖泊科学, 2019, 31(4): 1055-1063. DOI:10.18307/2019.0409]
[5]
Jiang XY, Zhang L, Yao XL et al. Greenhouse gas flux at reservoirs of Jiangxi Province and its influencing factors. J Lake Sci, 2017, 29(4): 1000-1008. [姜星宇, 张路, 姚晓龙等. 江西省水库温室气体释放及其影响因素分析. 湖泊科学, 2017, 29(4): 1000-1008. DOI:10.18307/2017.0424]
[6]
Wang XF, Yuan XZ, Chen H et al. Review of CO2 and CH4 emissions from rivers. Environmental Science, 2017, 38(12): 5352-5366. [王晓锋, 袁兴中, 陈槐等. 河流CO2与CH4排放研究进展. 环境科学, 2017, 38(12): 5352-5366. DOI:10.13227/j.hjkx.201701049]
[7]
Raymond PA, Hartmann J, Lauerwald R et al. Global carbon dioxide emissions from inland waters. Nature, 2013, 503(7476): 355-359. DOI:10.1038/nature12760
[8]
Yang P, He QH, Huang JF et al. Fluxes of greenhouse gases at two different aquaculture ponds in the coastal zone of southeastern China. Atmospheric Environment, 2015, 115: 269-277. DOI:10.1016/j.atmosenv.2015.05.067
[9]
Zhang WS, Li HP, Xiao QT et al. Urban Rivers are hotspots of riverine greenhouse gas (N2O, CH4, CO2) emissions in the mixed-landscape Chaohu Lake basin. Water Research, 2021, 189: 116624. DOI:10.1016/j.watres.2020.116624
[10]
Wen ZD, Song KS, Zhao Y et al. Carbon dioxide and methane supersaturation in lakes of semi-humid/semi-arid region, Northeastern China. Atmospheric Environment, 2016, 138: 65-73. DOI:10.1016/j.atmosenv.2016.05.009
[11]
Wang XF, He YX, Yuan XZ et al. Greenhouse gases concentrations and fluxes from subtropical small reservoirs in relation with watershed urbanization. Atmospheric Environment, 2017, 154: 225-235. DOI:10.1016/j.atmosenv.2017.01.047
[12]
Havens KE, Fukushima T, Xie P et al. Nutrient dynamics and the eutrophication of shallow Lakes Kasumigaura (Japan), Donghu (PR China), and Okeechobee (USA). Environmental Pollution, 2001, 111(2): 263-272. DOI:10.1016/S0269-7491(00)00074-9
[13]
Sun HY, Lu XX, Yu RH et al. Eutrophication decreased CO2 but increased CH4 emissions from lake: A case study of a shallow Lake Ulansuhai. Water Research, 2021, 201: 117363. DOI:10.1016/j.watres.2021.117363
[14]
胡茂俊. 富营养化自然水体中气体释放特征研究[学位论文]. 南京: 南京农业大学, 2016.
[15]
Wang H, Yu LF, Zhang ZH et al. Molecular mechanisms of water table lowering and nitrogen deposition in affecting greenhouse gas emissions from a Tibetan alpine wetland. Global Change Biology, 2017, 23(2): 815-829. DOI:10.1111/gcb.13467
[16]
Schimel JP, Gulledge J. Microbial community structure and global trace gases. Global Change Biology, 1998, 4(7): 745-758. DOI:10.1046/j.1365-2486.1998.00195.x
[17]
Yang P, Tong C. Emission paths and measurement methods for greenhouse gas fluxes from freshwater ecosystems: A review. Acta Ecologica Sinica, 2015, 35(20): 6868-6880. [杨平, 仝川. 淡水水生生态系统温室气体排放的主要途径及影响因素研究进展. 生态学报, 2015, 35(20): 6868-6880. DOI:10.5846/stxb201406231298]
[18]
杜鹏程. 水库水-气界面温室气体通量监测方法对比研究[学位论文]. 宜昌: 三峡大学, 2019.
[19]
Xie H, Long L, Mu XH et al. Comparison of floating static chamber and boundary layer equation for measuring CO2 and CH4 emission from urban river. Journal of China Three Gorges University: Natural Sciences, 2019, 41(5): 79-83. [谢恒, 龙丽, 穆晓辉等. 城市水体CO2和CH4通量监测的静态箱法与薄边界层模型估算法比较. 三峡大学学报: 自然科学版, 2019, 41(5): 79-83.]
[20]
Duchemin E, Lucotte M, Canuel R. Comparison of static chamber and thin boundary layer equation methods for measuring greenhouse gas emissions from large water bodies. Environmental Science & Technology, 1999, 33(2): 350-357. DOI:10.1021/es9800840
[21]
Liang NS, Nakadai T, Hirano T et al. In situ comparison of four approaches to estimating soil CO2 efflux in a northern larch (Larix kaempferi Sarg.) forest. Agricultural and Forest Meteorology, 2004, 123(1/2): 97-117. DOI:10.1016/j.agrformet.2003.10.002
[22]
Pumpanen J, Kolari P, Ilvesniemi H et al. Comparison of different chamber techniques for measuring soil CO2 efflux. Agricultural and Forest Meteorology, 2004, 123(3/4): 159-176. DOI:10.1016/j.agrformet.2003.12.001
[23]
Chang SQ, Wang DQ, Yu L et al. Greenhouse gas emission characteristics from urban rivers in Shanghai. Research of Environmental Sciences, 2015, 28(9): 1375-1381. DOI: 10.13198j.issn.1001-6929.2015.09.06. [常思琦, 王东启, 俞琳等. 上海城市河流温室气体排放特征及其影响因素. 环境科学研究, 2015, 28(9): 1375-1381. ]
[24]
Huang T, Wang XF, Liu TT et al. Spatiotemporal variations and influencing factors of CO2 and CH4 emissions from urban small landscape waters. Acta Ecologica Sinica, 2021, 41(15): 6024-6037. [黄婷, 王晓锋, 刘婷婷等. 城市小型景观水体CO2与CH4排放特征及影响因素. 生态学报, 2021, 41(15): 6024-6037. DOI:10.5846/stxb202005141226]
[25]
Wanninkhof R. Relationship between wind speed and gas exchange over the ocean. Journal of Geophysical Research: Oceans, 1992, 97(C5): 7373-7382. DOI:10.1029/92JC00188
[26]
Cole JJ, Caraco NF. Atmospheric exchange of carbon dioxide in a low-wind oligotrophic lake measured by the addition of SF6. Limnology and Oceanography, 1998, 43(4): 647-656. DOI:10.4319/lo.1998.43.4.0647
[27]
Crusius J, Wanninkhof R. Gas transfer velocities measured at low wind speed over a lake. Limnology and Oceanography, 2003, 48(3): 1010-1017. DOI:10.4319/lo.2003.48.3.1010
[28]
Raymond PA, Cole JJ. Gas exchange in rivers and estuaries: Choosing a gas transfer velocity. Estuaries, 2001, 24(2): 312-317. DOI:10.2307/1352954
[29]
杨平. 福建主要河口滨海围垦养虾塘养殖期CO2和CH4动态研究[学位论文]. 福州: 福建师范大学, 2017.
[30]
Liang JH, Tian LL, Zhou ZY et al. Characteristics and drivers of dissolved carbon dioxide and methane concentrations in the nantiaoxi river system in the upper reaches of the Taihu Lake basin during summer-autumn. Environmental Science, 2021, 42(6): 2826-2838. [梁佳辉, 田琳琳, 周钟昱等. 太湖流域上游南苕溪水系夏秋季水体溶存二氧化碳和甲烷浓度特征及影响因素. 环境科学, 2021, 42(6): 2826-2838. DOI:10.13227/j.hjkx.202009095]
[31]
Gao J, Zheng XH, Wang R et al. Preliminary comparison of the static floating chamber and the diffusion model methods for measuring water-atmosphere exchanges of methane and nitrous oxide from inland water bodies. Climatic and Environmental Research, 2014, 19(3): 290-302. [高洁, 郑循华, 王睿等. 漂浮通量箱法和扩散模型法测定内陆水体CH4和N2O排放通量的初步比较研究. 气候与环境研究, 2014, 19(3): 290-302.]
[32]
Wen ZD, Song KS, Zhao Y et al. Seasonal variability of greenhouse gas emissions in the urban lakes in Changchun, China. Environmental Science, 2016, 37(1): 102-111. [温志丹, 宋开山, 赵莹等. 长春城市水体夏秋季温室气体排放特征. 环境科学, 2016, 37(1): 102-111. DOI:10.13227/j.hjkx.2016.01.015]
[33]
韩洋. 南京市河流温室气体排放通量及其影响因素的研究[学位论文]. 南京: 南京信息工程大学, 2013.
[34]
Xiao QT, Duan HT, Qi TC et al. Environmental investments decreased partial pressure of CO2 in a small eutrophic urban lake: Evidence from long-term measurements. Environmental Pollution, 2020, 263: 114433. DOI:10.1016/j.envpol.2020.114433
[35]
Huttunen JT, Lappalainen KM, Saarijärvi E et al. A novel sediment gas sampler and a subsurface gas collector used for measurement of the ebullition of methane and carbon dioxide from a eutrophied lake. Science of the Total Environment, 2001, 266(1/2/3): 153-158. DOI:10.1016/S0048-9697(00)00749-X
[36]
Yang P, Tong C, He QH et al. Greenhouse gases fluxes at water-air interface of aquaculture ponds and influencing factors in the Min River Estuary. Acta Scientiae Circumstantiae, 2013, 33(5): 1493-1503. [杨平, 仝川, 何清华等. 闽江口鱼虾混养塘水-气界面温室气体通量及主要影响因子. 环境科学学报, 2013, 33(5): 1493-1503. DOI:10.13671/j.hjkxxb.2013.05.018]
[37]
Peng X, Wang BL, Liu CQ et al. 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
[38]
Davies JM, Hesslein RH, Kelly CA et al. PCO2 method for measuring photosynthesis and respiration in freshwater lakes. Journal of Plankton Research, 2003, 25(4): 385-395. DOI:10.1093/plankt/25.4.385
[39]
Yang K, Liu LX, Han YW et al. Effects of eutrophication on greenhouse gas emissions from urban water. Environmental Science & Technology, 2023, 46(S1): 231-236. [杨科, 刘丽香, 韩永伟等. 富营养化对城市水体温室气体排放的影响. 环境科学与技术, 2023, 46(S1): 231-236.]
[40]
Qi TC, Xiao QT, Miao YQ et al. Temporal and spatial variation of carbon dioxide concentration and its exchange fluxes in Lake Chaohu. J Lake Sci, 2019, 31(3): 766-778. [齐天赐, 肖启涛, 苗雨青等. 巢湖水体二氧化碳浓度时空分布特征及其水-气交换通量. 湖泊科学, 2019, 31(3): 766-778. DOI:10.18307/2019.0315]
[41]
Hope D, Kratz TK, Riera JL. Relationship between and dissolved organic carbon in northern Wisconsin Lakes. Journal of Environmental Quality, 1996, 25(6): 1442-1445. DOI:10.2134/jeq1996.00472425002500060039x
[42]
Han Y, Zheng YF, Wu RJ et al. Greenhouse gases emission characteristics of Nanjing typical waters in Spring. China Environmental Science, 2013, 33(8): 1360-1371. [韩洋, 郑有飞, 吴荣军等. 南京典型水体春季温室气体排放特征研究. 中国环境科学, 2013, 33(8): 1360-1371.]
[43]
Huttunen JT, Alm J, Liikanen A et al. Fluxes of methane, carbon dioxide and nitrous oxide in boreal lakes and potential anthropogenic effects on the aquatic greenhouse gas emissions. Chemosphere, 2003, 52(3): 609-621. DOI:10.1016/S0045-6535(03)00243-1
[44]
Richey JE, Melack JM, Aufdenkampe AK et al. Outgassing from Amazonian rivers and wetlands as a large tropical source of atmospheric CO2. Nature, 2002, 416(6881): 617-620. DOI:10.1038/416617a
[45]
Abril G, Martinez JM, Artigas LF et al. Amazon River carbon dioxide outgassing fuelled by wetlands. Nature, 2014, 505(7483): 395-398. DOI:10.1038/nature12797
[46]
Matthews CJD, St Louis VL, Hesslein RH. Comparison of three techniques used to measure diffusive gas exchange from sheltered aquatic surfaces. Environmental Science & Technology, 2003, 37(4): 772-780. DOI:10.1021/es0205838
[47]
Vachon D, Prairie YT, Cole JJ. The relationship between near-surface turbulence and gas transfer velocity in freshwater systems and its implications for floating chamber measurements of gas exchange. Limnology and Oceanography, 2010, 55(4): 1723-1732. DOI:10.4319/lo.2010.55.4.1723
[48]
Gong XJ, Wang XF, Yuan XZ et al. Effects of field towns development on the dissolved and diffusion fluxes of greenhouse gases in Heishuitan River basin, Chongqing. Acta Ecologica Sinica, 2019, 39(22): 8425-8441. [龚小杰, 王晓锋, 袁兴中等. 城镇快速发展对河流温室气体溶存及扩散通量的影响——以重庆市黑水滩河流域场镇为例. 生态学报, 2019, 39(22): 8425-8441. DOI:10.5846/stxb201810112205]
[49]
Ramamurthy P, Pardyjak ER. Toward understanding the behavior of carbon dioxide and surface energy fluxes in the urbanized semi-arid Salt Lake Valley, Utah, USA. Atmospheric Environment, 2011, 45(1): 73-84. DOI:10.1016/j.atmosenv.2010.09.049
[50]
Morales-Pineda M, Cózar A, Laiz I et al. Daily, biweekly, and seasonal temporal scales of pCO2 variability in two stratified Mediterranean reservoirs. Journal of Geophysical Research: Biogeosciences, 2014, 119(4): 509-520. DOI:10.1002/2013JG002317
[51]
Finlay K, Leavitt PR, Patoine A et al. Magnitudes and controls of organic and inorganic carbon flux through a chain of hard-water lakes on the northern Great Plains. Limnology and Oceanography, 2010, 55(4): 1551-1564. DOI:10.4319/lo.2010.55.4.1551
[52]
Loken LC, Crawford JT, Schramm PJ et al. Large spatial and temporal variability of carbon dioxide and methane in a eutrophic lake. Journal of Geophysical Research: Biogeosciences, 2019, 124(7): 2248-2266. DOI:10.1029/2019jg005186
[53]
Li DP, Jiang JBL, Xiang Y et al. Greenhouse gas emissions of typical waterbody in summer in Urumqi city. Environmental Protection Science, 2020, 46(6): 122-128. [李典鹏, 蒋靖佰伦, 向殷等. 乌鲁木齐市典型水体夏季温室气体排放初探. 环境保护科学, 2020, 46(6): 122-128.]
[54]
Xing YP, Xie P, Yang H et al. Methane and carbon dioxide fluxes from a shallow hypereutrophic subtropical Lake in China. Atmospheric Environment, 2005, 39(30): 5532-5540. DOI:10.1016/j.atmosenv.2005.06.010
[55]
Campeau A, Lapierre JF, Vachon D et al. Regional contribution of CO2 and CH4fluxes from the fluvial network in a lowland boreal landscape of Québec. Global Biogeochemical Cycles, 2014, 28(1): 57-69. DOI:10.1002/2013gb004685
[56]
Teodoru CR, Nyoni FC, Borges AV et al. Dynamics of greenhouse gases (CO2, CH4, N2O) along the Zambezi River and major tributaries, and their importance in the riverine carbon budget. Biogeosciences, 2015, 12(8): 2431-2453. DOI:10.5194/bg-12-2431-2015
[57]
Zhu RB, Liu YS, Xu H et al. Carbon dioxide and methane fluxes in the littoral zones of two lakes, east Antarctica. Atmospheric Environment, 2010, 44(3): 304-311. DOI:10.1016/j.atmosenv.2009.10.038
[58]
Crawford JT, Striegl RG, Wickland KP et al. Emissions of carbon dioxide and methane from a headwater stream network of interior Alaska. Journal of Geophysical Research: Biogeosciences, 2013, 118(2): 482-494. DOI:10.1002/jgrg.20034