湖泊科学   2023, Vol. 35 Issue (6): 1853-1865.  DOI: 10.18307/2023.0601
0

专稿

引用本文 [复制中英文]

高扬, 罗飚, 沈迪, 贾珺杰, 陆瑶, 王朔月, 青藏高原水体碳源汇过程的重新认知与挑战. 湖泊科学, 2023, 35(6): 1853-1865. DOI: 10.18307/2023.0601
[复制中文]
Gao Yang, Luo Biao, Shen Di, Jia Junjie, Lu Yao, Wang Shuoyue. Recognition and challenges of the inland water carbon source and sink processes on the Qinghai-Tibet Plateau. Journal of Lake Sciences, 2023, 35(6): 1853-1865. DOI: 10.18307/2023.0601
[复制英文]

基金项目

国家杰出青年科学基金项目(42225103)和中国科学院稳定支持基础研究领域青年团队计划(YSBR-037)联合资助

通信作者

高扬, E-mail: gaoyang@igsnrr.ac.cn

文章历史

2023-06-19 收稿
2023-09-01 收修改稿

码上扫一扫

青藏高原水体碳源汇过程的重新认知与挑战
高扬1,2 , 罗飚1,3 , 沈迪4 , 贾珺杰1 , 陆瑶1,2 , 王朔月1,2     
(1: 中国科学院地理科学与资源研究所, 生态系统网络观测与模拟重点实验室, 北京 100101)
(2: 中国科学院大学资源与环境学院, 北京 100049)
(3: 中国石油天然气集团有限公司大港油田分公司第三采油厂, 沧州 061023)
(4: 兰州大学资源与环境学院, 兰州 730000)
摘要:陆地-大气和陆地-水体的碳输送与碳交换共同决定了陆地生态系统的碳平衡, 但长期以来青藏高原水体由于缺乏有效的观测数据, 导致青藏高原水体一直是全球碳循环研究中被忽视的碳汇功能区。本研究聚焦青藏高原水体, 重新认知水体碳侵蚀与碳沉积对青藏高原固碳的重要作用, 系统分析河湖库-大气间的碳运移与交换时空格局与驱动机制。通过归纳总结发现: (1)青藏高原河流岩石风化速率以及CO2消耗速率较高, 强烈的侵蚀作用对地球系统吸收和平衡CO2浓度起着重要的调控作用; (2)青藏高原湖泊是内陆生态系统的重要碳库, 12 ka以来湖泊沉积物碳储量约为73 kg C/m2; (3)青藏高原湖泊浮游植物和水生植物影响内陆水体C沉积, 年均初级生产力达(553±36) mg C/(m2·d)。(4)青藏高原湖泊目前为“碳源”, 碳排放量约为2.27 Tg C/a。此外, 本文还阐明青藏高原水体碳源汇过程评估的不确定性, 以及从方法学角度分析如何加强其水体的碳循环研究, 厘清全球变化下青藏高原水体的碳源汇机制, 并展望了青藏高原水体碳源汇功能未来研究所面临的挑战, 从而为深入理解全球变化下青藏高原生态系统在碳中和过程中的作用提供科学依据。
关键词碳迁移    碳转化    碳平衡    碳中和    碳交换    碳收支    青藏高原    
Recognition and challenges of the inland water carbon source and sink processes on the Qinghai-Tibet Plateau
Gao Yang1,2 , Luo Biao1,3 , Shen Di4 , Jia Junjie1 , Lu Yao1,2 , Wang Shuoyue1,2     
(1: Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, P.R. China)
(2: College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, P.R. China)
(3: China National Petroleum Corporation, Dagang Oilfield Company No. 3 Oil Extraction Plant, Cangzhou 061023, P.R. China)
(4: College of Earth and Environment Science, Lanzhou University, Lanzhou 730000, P.R. China)
Abstract: Carbon(C) transport and C exchange between land atmosphere and land water body jointly determine the C balance of terrestrial ecosystem. However, due to the lack of effective observation data for a long time, inland water body on the Qinghai-Tibet Plateau (QTP) are unheeded C source or sink in the study of global C cycle. This study focus on the inland water body on the Qinghai-Tibet Plateau, to reexplore that the importance of C erosion and C deposition to the C sequestration potential of terrestrial ecosystem, systematically explained the C transportation and exchange processes between rivers, lakes, reservoirs, and the atmosphere and associated driving mechanism. The results showed that: (1) the part of rivers on the QTP had high weathering rates and CO2 consumption rates, and strong erosion played an important regulatory role in the absorption and balance of CO2 concentration by the Earth system; (2) the lakes on the QTP were an important carbon reservoir for inland ecosystems, and the carbon storage of lake sediments since 12 ka was about 73 kg C/m2; (3) the phytoplankton and aquatic plants in lakes on the QTP affected the deposition of C in inland water bodies, with an average annual primary productivity of (553 ± 36) mg C/(m2·d); (4) the lakes on the QTP were currently the "carbon source" with a carbon emission of approximately 2.27 Tg C/a. In addition, this study clarified the uncertainty estimation on the C source and sink process of the inland water, analyzed how to strengthen the C cycle research from a methodological perspective, and clarified the C source and sink mechanism of the QTP water body under global change, and looked forward to the research challenge on C source and sink function the inland water body on the QTP in the future. All these provide a scientific basis for in-depth understanding of the role of terrestrial ecosystems in the C neutralization process under global change.
Keywords: C transport    C transformation    C balance    C neutrality    C exchange    C budget    Qinghai-Tibet Plateau    

二氧化碳(CO2)在陆-水界面、水-气界面和水-底泥界面通过无机碳(IC)库转化过程、有机碳(OC)库转化过程、无机与有机碳库之间的转化过程进行碳的固定和释放[1]。CO2从内陆水体转化到大气中,即所谓的“碳逃逸”,是全球碳循环的一部分[2],全球每年CO2逃逸量约为2.3 Pg C[3],其中15% 来自于内陆湖泊和水库,85% 来自河流。因此,内陆水体在全球碳收支和碳循环中都发挥着关键作用[4-5]。内陆水体的浮游植物和水生植物、异养生物以及微生物等固碳潜力巨大[6],它们主要通过IC与OC库转化、有机质的降解和碳的吸收或埋藏进行固碳[6-7]。此外,岩石矿物的化学风化长期以来一直被认为是去除大气中CO2的关键机制,在岩石化学风化过程中,流域通过水岩反应吸收大气CO2,产生流域碳汇作用,从而降低大气中CO2[8-9],但目前对定量理解化学风化如何响应侵蚀过程及其对全球碳循环的影响仍然有限。

近年来,随着全球变化的加剧,各环境要素通过影响水体的生物-物理-化学过程不断影响着内陆水体的碳交换过程[6]。例如,光照、温度、养分和金属改变水体生物的光合和呼吸过程[10-12],进而影响水-气之间的CO2交换以及水体碳源汇特征。降水量会通过对土壤的淋溶作用影响湖泊有机碳积累,进而影响湖泊沉积物对碳的固定。pH[13]和氧化还原电位[14]能够改变水体的物理化学环境进而调节水体碳酸盐体系,从而决定CO2的产生释放[15]。风速[16]和温度[17]也会通过控制气体传递速率和水中气体溶解度的变化驱动着水-气界面的CO2溶解和扩散过程。因此,准确定量水体碳收支是深入理解陆地生态系统碳循环和实现碳中和的重要前提。

青藏高原湖泊总面积达到42521 km2,占到全国湖泊总面积的57%,湖泊不仅在时间和空间上存在着剧烈的变化,而且自然地理环境也差异巨大;青藏高原湖泊跨越1000~6000 m海拔高度,其中4000~5000和5000~6000 m海拔高度范围内的湖泊数量分别为840和158个,合计占比达到88.5% [18]。青藏高原由于其独特的高寒气候,形成了以淡水、咸水和盐湖为主的较为原始和多样化的湖泊,不同类型的水体导致内陆水体碳交换量估算存在很大的不确定性[19-20]。例如,全球盐湖二氧化碳年排放量为0.11~0.15 Pg C/a,约占全球湖泊的39%~47% [21]。然而,青藏高原的盐湖作为碳汇,每年吸收了(10.28±1.65) Tg C的CO2[22]。此外,由于水生生物的代谢作用而产生的CO2,以及水体中有机碳经过分解与矿化形成无机碳脱气致使淡水湖成为已知的CO2逃逸到大气中的主要管道[21]。因此,本文旨在厘清全球变化下内陆水体的碳源汇过程及其驱动机制,重新认知青藏高原碳源汇功能与未来研究的挑战,为如何应对全球变化下青藏高原水体在碳中和过程中的作用提供重要科学依据。

1 青藏高原水体碳侵蚀与碳沉积特征 1.1 风化、侵蚀与CO2消耗

侵蚀过程中的化学作用主要表现为岩石的化学风化,而化学风化作用主要包括溶解、水解、水化、碳酸盐化、氧化和生物化学风化作用。根据风化和侵蚀作用的强度和性质,物理风化和化学风化所起的作用不同,化学风化的基本反应也有差别,使得形成地表物质的组成也不相同。侵蚀过程通过硅酸盐风化从土壤中调动陆地植被的OC,将其转移到河流和沉积物中,从而降低大气CO2;同时通过氧化含有机碳或硫化物的岩矿物也可释放CO2(图 1)。侵蚀、风化和碳循环之间耦合的前提是硅酸盐矿物和碳酸的反应,了解风化作用对碳中和的影响需要区分硅酸盐风化作用和碳酸盐风化作用(图 1)。硅酸盐矿物风化过程中吸收的大气-水体中的CO2,通过河流进入海洋后部分以碳酸盐的形式沉积(生物或化学沉积),被称为永久性无机碳汇(substantial inorganic carbon sink);碳酸盐风化吸收的CO2,通过河流进入水体后全部释放,被称为非永久性碳汇(insubstantial inorganic carbon sink)[23]

图 1 风化侵蚀的碳源汇过程 (紫色实心箭头代表风化和侵蚀碳过程,紫色虚线箭头代表风化作用产生的CO2流向) Fig.1 Carbon source and sink processes of weathering erosion (The solid purple arrows represent weathering and erosion carbon processes, and the dashed purple arrows represent the flow of CO2 from weathering)

在岩石风化过程中,如温度、降水和生物生产力等都会影响CO2通量[24-25]。适宜的温度是岩石风化的重要因素,温度过低不利于岩石风化,过高的温度使CO2更容易从水中析出。温度变化还会改变生物光合作用及活性,水生光合生物利用水体溶解无机碳形成有机碳的稳定碳汇,并在地表水和地下水中实现碳循环交换,即形成了生物光合作用和碳酸盐风化耦合的岩溶动力系统,进而影响岩溶碳汇[26]。径流量也是影响流域碳汇量的直接因素,径流量越大,岩石风化速率越大,从而形成的岩溶碳汇量越大[27]。此外,近年来由于大气沉降(酸雨),硫化物氧化而形成的硫酸会对碳酸盐进行溶解,进而也会形成一部分的碳,因此外源酸在岩石风化溶解中也发挥了极大的作用[28]

在全球尺度上,相比于通过硅酸盐风化的固碳速率,物理侵蚀速率更能剧烈影响生物圈POC在水域沉积物中的埋藏速率[29]。成岩OC通过河流输出与沉积物输出紧密相连,通过快速沉积并转移到河流系统中得以有效保存[30]。目前,大多数与侵蚀相关的净CO2收支模型中,没有考虑OC埋藏和氧化风化作用,这限制了定量预测由气候和构造运动导致的POC通量变化[9]。物理剥蚀速率和硅酸盐风化通量取决于侵蚀物质供应和化学反应速度,而碳酸盐风化与硅酸盐风化形成的溶质比例通常随侵蚀率而增加[31]。在低侵蚀率下,风化产物积聚在土壤中,保护下伏岩石免受与地表流体的相互作用,导致“供应有限”的风化[32-33]。在供应有限的条件下,随着侵蚀率的增加,风化和侵蚀呈线性扩展,暴露出更多的地表物质参与反应[34]。相反,在较高的侵蚀率下,矿物的供应速度与它们的反应速度不同,通量被“风化限制”(也被称为“动力学限制”)并由其他因素控制,如CO2浓度、流体流速和温度[35]。尽管侵蚀率与总风化通量之间存在明显的线性关系,但这种关系可能并不反映对长期CO2下降的侵蚀控制[36]

青藏高原是世界上独一无二的地域生态单元,在全球地球化学循环中发挥着重要的作用。其地质历史年轻,新构造运动活跃。新生代以来印度板块与亚洲板块碰撞造成青藏高原不断隆升[37],地势的不断隆升使高原地貌遭受强烈的内、外营利作用,地表处于不断的侵蚀、搬运和堆积过程中,大量岩石裸露加速了岩石化学风化过程与对大气CO2的消耗[38]。在岩石化学风化过程中,流域通过水岩反应吸收大气CO2,产生流域碳汇作用。在青藏高原地区,主要以硅酸盐与碳酸盐为主参与了流域的岩石风化作用。有研究显示青藏高原拉萨河流域碳酸盐风化速率为21.6~33.6 t/(km2·a),硅酸盐风化速率为3.61~6.62 t/(km2·a),对应的CO2净消耗通量分别为88.9×103和118×103 mol/(km2·a)[39]。尼洋河硅酸盐与碳酸盐的岩石风华速率分别为4.46和41.50 t/(km2·a),CO2净消耗通量分别为92×103和350×103 mol/(km2·a)[40]。此外,还有研究结果表明青藏高原南部的雅鲁藏布江流域全年岩石风化速率为42.03 t/(km2·a),流域岩石风化CO2消耗速率平均为407.13×103 mol/(km2·a)[41],而在全世界范围内,非喀斯特流域的平均CO2吸收通量约为158.89×103 mol/(km2·a)[42],因此,青藏高原的CO2消耗速率是高于世界河流的平均水平,青藏高原强烈的侵蚀作用对地球系统吸收和平衡CO2浓度起着重要的调控作用,进而影响着全球气候以及生态系统的变化。

1.2 沉积储碳

水体中的碳通过溶解泵和生物泵等过程导致的有机碳由水体表层向深层乃至水底的转移,将一部分碳封存在水底沉积物中长期不参与大气CO2循环, 起到“沉积物储碳”的作用[43-44]。对于许多水坝河流系统来说,持续的富营养化正在驱使水库趋向于自养,因为营养物质浓度的增加使得浮游生物群落的光合作用超过了呼吸作用,这种生产力变化的结果是增加了埋藏在沉积物中的碳[45]。内陆水体的沉积碳由无机碳和有机碳两部分组成,既有内源如湖泊内的浮游生物和大型水生植物等,也有外源如径流携带的陆源碎屑和土壤有机碳,一般可分为藻类、草本类、木质类和木炭4种类型[46];无机碳则来源于湖泊自生碳酸盐和外源碳酸盐[47]。外源碳酸盐是由湖盆流域母岩风化产生,由地表径流搬运至湖泊水体的碳酸盐,而湖泊自生碳酸盐包括湖水中无机化学沉淀产生的碳酸盐和生物壳体碳酸盐,以及少量沉积物埋藏后早期成岩作用产生的碳酸盐[48-49]。水生植物通过光合作用从大气中吸收CO2,以及水生生物死亡后降解沉降到底泥[50-51];自然界中的碳还会通过径流、土壤以及生物等陆地环境要素的传送作用进入湖泊,水体中POC会沉降到底部[52]。当大气中的pCO2pCH4大于水体中的pCO2pCH4时,CO2与CH4就会进入并溶解于水体,同时与沉积物进行碳交换[53]。因此,内陆水体的碳埋藏量是被忽视的一部分,存在巨大的不确定性(图 2)。

图 2 内陆水体碳源汇过程 Fig.2 Process of carbon source and sink in inland water

内陆水体沉积碳受多种因素的影响。在一些水热条件变化大的区域,温度[54]以及降水[55]对湖泊碳埋藏速率会产生较为显著的调控作用。温度的高低会影响水体中细菌对沉积物的分解以及影响湖泊初级生产力,而降水的多少不仅会影响湖泊水量,同时还通过对土壤的淋溶作用影响湖泊有机碳积累,进而影响湖泊沉积物的碳埋藏速率。一些湖泊的理化性质[56],例如湖泊的营养状况,湖泊水体类型以及湖泊中的溶解氧气等同样会影响沉积物碳积累。底层水氧气越少,越有利于有机碳储存,有机碳的埋藏效率也越高。此外,一些生物因素[57](如湖泊生物量、湖泊中生物类型)以及人类活动[58](土地利用类型的变换)也会加快或者减缓有机碳的埋藏速率。有研究表明在农田施肥后,部分化肥会随着地表径流进入湖泊,由此对湖泊中的营养多寡造成影响,进而影响湖泊有机碳储量。

青藏高原由于其复杂的地貌特征与气候条件,形成了世界上海拔最高、范围最广、数量最多,以盐湖与咸水湖集中为特色的高原湖泊群,而且青藏高原受人类活动等外界的扰动较小,在静水条件下,水体中的碳更容易沉降到湖泊底部,因此,青藏高原湖泊可以更大限度的将碳埋藏在沉积物中。目前,一些研究对青藏高原湖泊沉积物中的碳埋藏速率以及碳储量进行了估算。研究结果显示哈拉湖沉积物的碳埋藏速率为4.71 g C/(m2·a),尕海沉积物的碳埋藏速率为14.02 g C/(m2·a),兹格塘错沉积物的碳埋藏速率为6.80 g C/(m2·a)[59],青海湖沉积物的碳埋藏速率约为7.23 g C/(m2·a)[60],而欧洲一些湖泊沉积物碳埋藏速率大概为5.6 g C/(m2·a)[61],格陵兰岛湖泊沉积物碳埋藏速率为6 g C/(m2·a)[62],而美国一些地区的湖泊沉积物碳埋藏速率大概为5.63 g C/(m2·a)[63],可见,青藏高原湖泊沉积物的碳埋藏速率与全世界上其他地区的许多湖泊记录是可相提并论的。此外,还有研究结果表明青藏高原12 ka以来的湖泊碳储量约为73(16~274)kg C/m2[64]。因此,青藏高原湖泊沉积物中储存着大量有机碳,是内陆生态系统的重要碳库,在区域碳循环过程中扮演着重要角色。

2 青藏高原水体碳固定与碳逃逸特征 2.1 水体初级生产固碳

浮游植物是水体中最主要的初级生产者,虽然只占地球生物圈初级生产者生物量的0.2% [65-66],但浮游植物固碳对生物圈净初级生产(NPP)的贡献率约50%,是OC和IC储库间循环的最重要的环节[67-68]。同时,水生植物通过光合作用吸收大气CO2合成有机质,积累大量的OC和IC的过程也不容忽视[69-70]。这些浮游植物和水生植物的初级生产对营养时间和碳动力学产生了强烈的连锁效应,它将直接影响内陆水体碳沉积,控制水柱中O2的可用性(通过从主要生产者直接释放或通过漂浮的大型植物限制通入大气和通向大气的通量),从而间接影响水柱中的好氧矿化速率和磷浓度,因此在小型和浅水水生系统的碳沉积在全球碳循环中起着至关重要的作用[71]

2002—2020年间青藏高原湖泊初级生产力(GPP)分布格局如图 3所示[72]。2000年以来,青藏高原湖泊年均初级生产力在477~622 mg C/(m2·d)范围内不断波动,其平均值达到(553 ± 36) mg C/(m2·d)。大多数湖泊系统(72%)的GPP随时间呈显著下降趋势,尤其是位于高原内部和南部地区。大多数GPP显著增加的湖泊系统主要位于青藏高原东南部和东北部。在空间分布上,青藏高原湖泊初级生产力呈现出西高东低、北高南低的地理格局(图 3a)。位于高原西部海拔在4291 m的咸水湖曼冬错具有最高的年均GPP(5598 mg C/(m2·d)),而位于高原南部海拔4685 m的姆错丙尼具有最低GPP(106 mg C/(m2·d))。大多数GPP较高的湖泊系统位于青藏高原南部, 只有少数GPP高于平均值的湖泊系统位于青藏高原西部。此外,海拔较低或盐度较低的湖泊系统具有较高初级生产力。

图 3 青藏高原湖泊(a)总初级生产力、(b)CO2交换通量以及(c)CO2交换量空间分布格局 Fig.3 Spatial distribution pattern of (a) total primary productivity, (b)CO2 flux rates and (c) CO2 exchange flux in lakes on the Qinghai-Tibet Plateau
2.2 水体C逃逸

内陆水体与大气之间进行了数量巨大的CO2交换,但目前缺乏对全球尺度上内陆水体与大气之间基于CO2交换量空间分布的系统性研究[2]。在全球范围内CO2逃逸总量评估主要取决于3个参数:水体中CO2[73]、全球内陆水体(河流、湖泊、水库等)的总面积[74]以及气体交换速率(k,决定气体交换速率的物理学常数)[75-76]。内陆水体通常CO2过饱和,这反映了大气水量平衡,通过水体碱度和pH值计算得到的内陆水体中CO2含量显著高于大气中CO2含量。Butman等研究表明,在美国境内北纬25°至50°的温带地区河流每年排放的CO2数量大约为5亿吨[77],而对于内陆大型河流以及北极水体也被证明为CO2过饱和状态[78]

湿润热带地区的单位面积CO2排放量最高,咸水湖的CO2排放量也不低,仅里海就贡献了20% 左右的全球湖泊碳排放,所有咸水湖则贡献了不到50% [21]。内陆水体中CO2逃逸是陆地呼吸的一部分,虽然从陆地到内陆水体的净初级生产力NPP输送仅占陆地NPP的4%,但每年异养呼吸和火灾与NPP的差值大约是15亿吨碳[79]。通过在英国流域的长期研究发现:约1.2%~2.2% 的陆地NPP从湖泊逃逸[80]。因此,忽略内陆水体碳排放可能造成区域尺度碳平衡的较大误差,而这部分在全球碳收支中十分重要[4]

青藏高原地区湖泊CO2交换通量形成了强差异化格局,从-16~140 mg CO2/(m2·h)不等(图 3b)。大部分湖泊CO2交换通量处在0~60 mg CO2/(m2·h)范围内。个别湖泊也具有固碳潜力,即CO2交换通量小于0。进一步统计了淡水湖和咸水湖的CO2平均交换通量(图 4),发现咸水湖CO2排放通量((38±34) mg CO2/(m2·h))高于淡水湖((33±28) mg CO2/(m2·h)),咸水湖的CO2吸收通量((-9±10) mg CO2/(m2·h))也高于淡水湖的吸收通量((-5±4) mg CO2/(m2·h))(绝对值相比)。因此,青藏高原咸水湖的平均净CO2排放通量((34±35) mg CO2/(m2·h))高于淡水湖((17±28) mg CO2/(m2·h))。原因可能是咸水湖盐度高,高盐毒害水生植物和浮游植物,导致较低的生物量和较弱的光合作用,但是较强的呼吸作用致使水中CO2分压高,CO2交换通量高,而淡水湖水生植物和浮游植物生长旺盛,较强的光合作用导致水中CO2分压低,CO2交换通量低[81]

图 4 青藏高原淡水湖与咸水湖CO2交换通量 Fig.4 CO2 flux rates between freshwater and saltwater lakes on the Qinghai-Tibet Plateau

青藏高原湖泊CO2交换量形成南低北高的格局,每年在-30万吨CO2到105万吨CO2范围内(图 3c)。整体上看,青藏高原湖泊是个巨大的碳源,每年平均净二氧化碳排放量达到(7±18)万吨,尤其是咸水湖高达(8±22)万吨,淡水湖低于咸水湖,但是每年平均净二氧化碳排放量也达到(3±4)万吨。咸水湖平均面积大于淡水湖,因此相比淡水湖来说,咸水湖既可能是一个巨大的碳汇,也可能是一个巨大的CO2排放源[81]

不同青藏高原湖泊的水温、盐度、碱度、pH等水环境要素的异质性制约了湖泊CO2交换通量,导致在利用有限采样数据升尺度评估青藏高原湖泊总体碳排放量时的高度不确定性。对青藏高原地区湖泊而言,本研究基于2003—2020年采样观测数据的评估结果为2.27 Tg C/a。Li等基于1998—2012年采样数据的估算结果为4.6 Tg C/a,置信区间范围为2.4~6.7 Tg C/a[18],而基于2014—2015年数据估算的CO2排放量达到10.1 Tg C/a[82];Ran等利用不同时期观测数据的估算结果为2.1 Tg C/a (1980s)和3.8 Tg C/a (2010s)[20](图 5)。以上研究也表明,青藏高原湖泊总体表现为碳源,每年CO2排放量范围在2~10 Tg C/a之间,而且其碳排放量随时间呈现显著上升趋势。除了各研究自身估算的不确定性,不同研究对青藏高原湖泊碳排放的估算结果也存在显著差异。造成显著差异的原因包括:1)不同研究所用数据的采样点位置不同;2)不同研究所用数据的采样时间不一致;3)受限于青藏高原高寒环境,采样观测数据量总体偏少。总体而言,当前研究结果表明青藏高原湖泊整体依然处于碳排放状态,但碳排放量估算仍然存在高度不确定性。

图 5 青藏高原湖泊CO2排放量不同研究估算结果对比 Fig.5 Comparison of CO2 emission estimates from different studies in lakes of the Qinghai-Tibet Plateau
2.3 水体温室气体排放

目前关于青藏高原河流CO2总体排放的研究当前还十分缺乏,1980s青藏高原河流在干季和湿季的总面积分别为(8902±4691)和(11091±5616) km2,2010s分别为(9662±5481)和(11368±5934) km2;平均CO2交换通量在1980s为(6656.4±7160.4) mg/(m2·d),到2010s为(5721.6±9362.4) mg/(m2·d);CO2总体排放量在1980s为(15.8±7.4) Tg C/a,到2010s上升至(18.6±7.9) Tg C/a[20]。据此推断,青藏高原河流CO2排放量是湖泊CO2排放量的4.9~7.5倍,占到青藏高原水体CO2排放量的83%~88%,是青藏高原水体排放的主要来源。相对于湖泊和河流,青藏高原水库的面积与数量较少。2010s统计结果共有水库239个,水库总面积(695.7±56) km2。青藏高原水库CO2排放量在1980s和2010s分别为(0.03±0.01)和(0.07±0.02) Tg C/a,仅占湖泊CO2排放量不到2% [20]

除了CO2以外,CH4也是水体碳交换的一个重要来源。青藏高原湖泊CH4的平均交换通量为5.2(0.0008~45.9) mmol/(m2·d)[82],淡水湖CH4净排放通量(4 mg CH4/(m2·h))普遍高于咸水湖(0.50 mg CH4/(m2·h))。青藏高原湖泊CH4总体排放量估算为1.6(0.45~2.7) Tg C/a, 约为青藏高原湖泊CO2排放量的35% [18]。青藏高原河流CH4的平均交换通量为11.9 mmol/(m2·d),约为湖泊CH4平均交换通量的2倍。Zhang等的长期观测分析结果指出,青藏高原东部主要河流的CH4排放量在0.37~1.23 Tg CH4/a[51]。因此,CH4也是青藏高原湖泊碳排放的重要组成部分。然而,由于缺乏足够的采样观测数据,青藏高原CH4排放量评估和CO2一样,还存在较大的不确定性,需要进一步开展相关研究。另一个主要的温室气体排放源是N2O。青藏高原湖泊N2O排放通量达到(2.16±1.77) mg N/(m2·d),年N2O总排放通量为(32.37±26.52) Gg N/a,占中国所有湖泊的88%,贡献了最多的中国N2O排放[18]。由于青藏高原海拔较高,湖泊缺氧加剧了湖泊沉积物中的硝酸盐反硝化过程形成N2O汇[83-84]。与CO2 (73.9%) 和CH4(25.7%)相比,湖泊N2O排放水平(0.4%) 极低[82]。总体而言,青藏高原湖泊的温室气体占中国湖泊总排放量的44%。

CO2和CH4净排放量从2000s的2.06和0.50 Tg C/a增加到2020s的4.58和0.61 Tg C/a,然而多数研究将夏季观测数据作为全年碳收支进行评估,这并未考虑冬季结冰期的碳收支情况。Li等的研究表明青藏高原所有盐湖吸收的CO2总量为(-10.28±1.65) Tg C/a,相当于青藏高原陆地生态系统净碳汇的1/3[20]。若考虑冬季季节性水化学过程和盐湖中湖冰对CO2的吸收效应,目前青藏高原CO2净排放量应低于此估值(4.58 Tg C/a),然而碳交换通量估算依然存在高度不确定性。自20世纪80年代初以来,青藏高原全面的气候变化如地表空气变暖和增湿、太阳变暗和风力减弱[85],显著影响了水体碳源汇效应。温度升高会导致水体生物物理作用加快以及碳排放增加[86],同时也会引起无冰期延长、冰川雪水融化[87-88]、湖泊数量面积增加,进而造成湖泊盐度降低和冰期变短,促进浮游植物生长,最终导致湖泊CO2吸收增强[22]。太阳变暗会引起水体自养生物的光合作用减弱进而减少对CO2的吸收[11]。风速减慢使气体扩散速率变慢以及气体交换通量变低[15]。未来青藏高原水体碳源汇向何方演变,仍需进一步加强结冰期的源汇观测评估。

3 青藏高原碳源汇过程研究挑战

水-土界面输入、水-气界面释放、光合作用、沉积作用等过程之间密切关联,河流碳输出通量是这些过程综合作用的结果,目前大多数研究主要针对单一过程独立开展,碳释放、累积、运输等过程的综合研究有待进一步加强[89]。随河流运移碳通量的估算方法主要有3类:通过碳含量与流量的乘积计算,利用水体中的碳浓度或通量观测数据与影响因素构建经验方程,以及生态水文与生物地球化学过程耦合模型[90-91]。直接利用碳含量与流量乘积的方法,受到观测数据时间序列或空间变化较少的限制,较难全面反应不同时间、空间尺度的碳通量[73]。生态水文与生物地球化学过程耦合模型通过链接陆-水-气界面的碳过程,综合观测数据与理论模型,可对不同气候、土地利用、土壤类型等条件下的水体碳传输通量进行实时动态、长时间序列、多空间位置的定量估算,如NICE-BGC模型[88-90]

解析陆-水、水-气、水-底泥等界面之间淋溶及侵蚀、光合/呼吸、溶解/沉淀、底泥释放/沉积等过程的碳通量是准确预测环境变化对水体碳通量与碳循环的影响,以及理解其控制机制的基础[53, 89]。区域尺度水土过程的碳流失通量较难直接观测估算得到,而生物地球化学过程模型综合考虑了不同条件下的物质循环过程,为定量估算水土过程的碳流失提供了有效工具[90]。水-气界面CO2通量直接准确测定较为困难,其定量计算方法有示踪法、箱式法、扩散模型法等[92-93]。当前内陆水体初级生产评估主要在特定区域利用实地采样数据来进行分析[6],而基于遥感观测与模型模拟的水体初级生产研究受制于地面采样与遥感反演技术,主要在局部区域开展[94],缺乏大范围、长时期、高精度的评估结果[82]。水体中碳沉积量估算主要是根据其沉积速率,并且水体中碳同位素组成的变化可以反应水体化学特征、水气界面CO2的交换程度、流域C3和C4植被变化情况、浮游植物光合/呼吸作用强弱等重要信息,同时可以结合C/N解析水体碳的不同来源[95]。此外,卫星遥感快捷、大范围和周期性的特点为陆地-河流-近海等生态系统碳循环研究提供了更好的研究手段,特别是随着拥有适合内陆水体水色波段的卫星发展,均可为开展基于卫星遥感的陆地-河流-近海的碳循环研究,提供了重要的技术支撑。

目前对于大多数水体碳源汇的研究,不论是基于实地采样或是通过遥感的方式来进行监测,基本都集中在开放水域的季节。而冰封期的水体由于在传统上被认为相对于开放水域季节的水体,在碳循环中并不那么重要,因此常常被排除在碳收支预算之外。但有研究发现,冰本身就可以从大气中吸收CO2。此外,随着冬季生态系统研究的深入,发现冰下栖息地并不像曾经假设的处于休眠状态,冰下生命活动仍然存在,碳交换在冰下持续进行着。因此,结冰期水体的碳交换量可能会对青藏高原内陆水体碳循环过程产生重要影响。

4 展望

未来如何加强青藏高原水体碳中和研究,首先需确定青藏高原水体碳汇(定义为CO2收支的剩余部分),对于青藏高原碳收支的闭合具有很重要的作用。青藏高原淡水湖、盐湖和咸湖随海拔梯度变化的巨大空间异质性,导致当前缺乏测量方法和足够数据来完全对青藏高原陆-水-气系统的异质性进行分类,以致不能将区域尺度的碳动态向全球尺度外推,而水-气CO2通量的不确定性将会对全球陆地生物圈与大气的碳交换及碳汇功能造成较大误差。在过去的几十年里,碳同位素方法提供了关于碳埋藏和风化速率的重要间接证据,未来研究需要运用新的技术与方法解决陆-水-气关键带碳收支中的有机和无机碳的去向问题,包括最新先进的碳来源、年龄和现代不稳定有机碳(如从浮游植物和现代维管植物)监测,以及新的CO2、pH和DOC的传感装置运用。同时提高耦合水动力的生物地球化学模型的参数敏感性分析,降低模型预测的不确定性,发展新一代耦合水动力的生物地球化学模型,从力学上考虑控制碳动力学的因素,如元素化学计量学和生物碳周转率等。

5 参考文献

[1]
Barker PA, Hurrell ER, Leng MJ et al. Carbon cycling within an East African lake revealed by the carbon isotope composition of diatom silica: A 25-ka record from Lake Challa, Mt. Kilimanjaro. Quaternary Science Reviews, 2013, 66: 55-63. DOI:10.1016/j.quascirev.2012.07.016
[2]
Le Quéré C, Raupach MR, Canadell JG et al. Trends in the sources and sinks of carbon dioxide. Nature Geoscience, 2009, 2(12): 831-836. DOI:10.1038/ngeo689
[3]
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
[4]
Regnier P, Friedlingstein P, Ciais P et al. Anthropogenic perturbation of the carbon fluxes from land to ocean. Nature Geoscience, 2013, 6(8): 597-607. DOI:10.1038/ngeo1830
[5]
Wanninkhof R, Park GH, Takahashi T et al. Global ocean carbon uptake: Magnitude, variability and trends. Biogeosciences, 2013, 10(3): 1983-2000. DOI:10.5194/bg-10-1983-2013
[6]
Gao Y, Jia JJ, Lu Y et al. Determining dominating control mechanisms of inland water carbon cycling processes and associated gross primary productivity on regional and global scales. Earth-Science Reviews, 2021, 213: 103497. DOI:10.1016/j.earscirev.2020.103497
[7]
Jia JJ, Gao Y, Zhou F et al. Identifying the main drivers of change of phytoplankton community structure and gross primary productivity in a river-lake system. Journal of Hydrology, 2020, 583: 124633. DOI:10.1016/j.jhydrol.2020.124633
[8]
White AF, Brantley SL. The effect of time on the weathering of silicate minerals: Why do weathering rates differ in the laboratory and field?. Chemical Geology, 2003, 202(3/4): 479-506. DOI:10.1016/j.chemgeo.2003.03.001
[9]
Hilton RG, West AJ. Mountains, erosion and the carbon cycle. Nature Reviews Earth & Environment, 2020, 1(6): 284-299. DOI:10.1038/s43017-020-0058-6
[10]
Strzepek RF, Boyd PW, Sunda WG. Photosynthetic adaptation to low iron, light, and temperature in Southern Ocean phytoplankton. Proceedings of the National Academy of Sciences of the United States of America, 2019, 116(10): 4388-4393. DOI:10.1073/pnas.1810886116
[11]
Tanabe Y, Hori M, Mizuno AN et al. Light quality determines primary production in nutrient-poor small lakes. Scientific Reports, 2019, 9(1): 4639. DOI:10.1038/s41598-019-41003-9
[12]
Weisse T, Gröschl B, Bergkemper V. Phytoplankton response to short-term temperature and nutrient changes. Limnologica, 2016, 59: 78-89. DOI:10.1016/j.limno.2016.05.002
[13]
Natchimuthu S, Selvam BP, Bastviken D. Influence of weather variables on methane and carbon dioxide flux from a shallow pond. Biogeochemistry, 2014, 119(1): 403-413. DOI:10.1007/s10533-014-9976-z
[14]
Singh SN, Kulshreshtha K, Agnihotri S. Seasonal dynamics of methane emission from wetlands. Chemosphere-Global Change Science, 2000, 2(1): 39-46. DOI:10.1016/S1465-9972(99)00046-X
[15]
Li SY, Lu XX, 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. Journal of Hydrology, 2012, 466/467: 141-150. DOI:10.1016/j.jhydrol.2012.08.011
[16]
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(Pt A): 114433. DOI:10.1016/j.envpol.2020.114433
[17]
Palma-Silva C, Marinho CC, Albertoni EF et al. Methane emissions in two small shallow neotropical lakes: The role of temperature and trophic level. Atmospheric Environment, 2013, 81: 373-379. DOI:10.1016/j.atmosenv.2013.09.029
[18]
Messager ML, Lehner B, Grill G et al. Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 2016, 7: 13603. DOI:10.1038/ncomms13603
[19]
Li SY, Bush RT, Santos IR et al. Large greenhouse gases emissions from China's Lakes and reservoirs. Water Research, 2018, 147: 13-24. DOI:10.1016/j.watres.2018.09.053
[20]
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
[21]
Duarte CM, Prairie YT, Montes C et al. CO2 emissions from saline lakes: A global estimate of a surprisingly large flux. Journal of Geophysical Research, 2008, 113: G04041. DOI:10.1029/2007jg000637
[22]
Li XY, Shi FZ, Ma YJ et al. Significant winter CO2 uptake by saline lakes on the Qinghai-Tibet Plateau. Global Change Biology, 2022, 28(6): 2041-2052. DOI:10.1111/gcb.16054
[23]
Curl RL. Carbon shifted but not sequestered. Science, 2012, 335(6069): 655-655. DOI:10.1126/science.335.6069.655-a
[24]
Walker JCG, Hays PB, Kasting JF. A negative feedback mechanism for the long-term stabilization of Earth's surface temperature. Journal of Geophysical Research: Oceans, 1981, 86: 9776-9782. DOI:10.1029/jc086ic10p09776
[25]
Berner RA, Lasaga AC, Garrels RM. The carbonate-silicate geochemical cycle and its effect on atmospheric carbon dioxide over the past 100 million years. American Journal of Science, 1982, 283: 641-683. DOI:10.2475/ajs.283.7.641
[26]
Zeng SD, Xia J, Chen XD et al. Integrated land-surface hydrological and biogeochemical processes in simulating water, energy and carbon fluxes over two different ecosystems. Journal of Hydrology, 2020, 582: 124390. DOI:10.1016/j.jhydrol.2019.124390
[27]
Liu ZH. Contribution of carbonate rock weathering to the atmospheric CO2 sink. Carsologica Sinica, 2000, 19(4): 293-300. [刘再华. 碳酸盐岩岩溶作用对大气CO2沉降的贡献. 中国岩溶, 2000, 19(4): 293-300. DOI:10.3969/j.issn.1001-4810.2000.04.001]
[28]
Xie YC, Huang F, Yang H et al. Role of anthropogenic sulfuric and nitric acids in carbonate weathering and associated carbon sink budget in a Karst Catchment (Guohua), southwestern China. Journal of Hydrology, 2021, 599: 126287. DOI:10.1016/j.jhydrol.2021.126287
[29]
Galy V, France-Lanord C, Beyssac O et al. Efficient organic carbon burial in the Bengal fan sustained by the Himalayan erosional system. Nature, 2007, 450(7168): 407-410. DOI:10.1038/nature06273
[30]
Bouchez J, Galy V, Hilton RG et al. Source, transport and fluxes of Amazon River particulate organic carbon: Insights from river sediment depth-profiles. Geochimica et Cosmochimica Acta, 2014, 133: 280-298. DOI:10.1016/j.gca.2014.02.032
[31]
Jacobson AD, Blum JD. Relationship between mechanical erosion and atmospheric CO2 consumption in the New Zealand Southern Alps. Geology, 2003, 31(10): 865. DOI:10.1130/g19662.1
[32]
West AJ, Galy A, Bickle M. Tectonic and climatic controls on silicate weathering. Earth and Planetary Science Letters, 2005, 235(1/2): 211-228. DOI:10.1016/j.epsl.2005.03.020
[33]
Dixon JL, von Blanckenburg F. Soils as pacemakers and limiters of global silicate weathering. Comptes Rendus Geoscience, 2012, 344(11/12): 597-609. DOI:10.1016/j.crte.2012.10.012
[34]
Gabet EJ, Mudd SM. A theoretical model coupling chemical weathering rates with denudation rates. Geology, 2009, 37(2): 151-154. DOI:10.1130/g25270a.1
[35]
Maher K, Chamberlain CP. Hydrologic regulation of chemical weathering and the geologic carbon cycle. Science, 2014, 343(6178): 1502-1504. DOI:10.1126/science.1250770
[36]
Larsen IJ, Montgomery DR, Greenberg HM. The contribution of mountains to global denudation. Geology, 2014, 42(6): 527-530. DOI:10.1130/g35136.1
[37]
Raymo ME, Ruddiman WF, Froelich PN. Influence of late Cenozoic Mountain building on ocean geochemical cycles. Geology, 1988, 16(7): 649. DOI:10.1130/0091-7613(1988)0160649:iolcmb>2.3.co;2
[38]
Li SL, Chetelat B, Yue FJ et al. Chemical weathering processes in the Yalong River draining the eastern Tibetan Plateau, China. Journal of Asian Earth Sciences, 2014, 88: 74-84. DOI:10.1016/j.jseaes.2014.03.011
[39]
Xie YC, Yu S, Miao XY et al. Chemical weathering and CO2 consumption flux in Tibetan Plateau: A case of Lhasa River Basin. Earth Science Frontiers, 2023, 30(5): 510-525. [谢银财, 于奭, 缪雄谊等. 青藏高原流域岩石风化机制及其CO2消耗通量: 以拉萨河为例. 地学前缘, 2023, 30(5): 510-525.]
[40]
Liu X, Zhang D, Gao S et al. Chemical weathering and CO2 consumption flux in Tibetan Plateau: A case of Niyang River Catchment. Chinese Journal of Ecology, 2018, 37(3): 688-696. [刘旭, 张东, 高爽等. 青藏高原小流域化学风化过程及其CO2消耗通量: 以尼洋河为例. 生态学杂志, 2018, 37(3): 688-696. DOI:10.13292/j.1000-4890.201803.041]
[41]
Li BG, Bao YF, Hu MM et al. Hydrochemical characteristics and carbon sink calculation in Yarlung Zangbo River Basin. Yangtze River, 2023, 54(4): 49-56. [李伯根, 包宇飞, 胡明明等. 雅鲁藏布江流域水化学特征及碳汇计算. 人民长江, 2023, 54(4): 49-56.]
[42]
Li CJ, Wang SJ, Bai XY et al. Estimation of carbonate rock weathering-related carbon sink in global major river basins. Acta Geographica Sinice, 2019, 74(7): 1319-1332. [李朝君, 王世杰, 白晓永等. 全球主要河流流域碳酸盐岩风化碳汇评估. 地理学报, 2019, 74(7): 1319-1332. DOI:10.11821/dlxb201907004]
[43]
Jiao NZ, Herndl GJ, Hansell DA et al. Microbial production of recalcitrant dissolved organic matter: Long-term carbon storage in the global ocean. Nature Reviews Microbiology, 2010, 8(8): 593-599. DOI:10.1038/nrmicro2386
[44]
Jiao NZ. Carbon fixation and sequestration in the ocean, with special reference to the microbial carbon pump. Sci Sin Terrae, 2012, 42: 1473-1486. DOI:10.1360/zd-2012-42-10-1473
[45]
Dean WE, Gorham E. Magnitude and significance of carbon burial in lakes, reservoirs, and peatlands. Geology, 1998, 26(6): 535. DOI:10.1130/0091-7613(1998)0260535:masocb>2.3.co;2
[46]
Nowicka B, Ciura J, Szymańska R et al. Improving photosynthesis, plant productivity and abiotic stress tolerance-current trends and future perspectives. Journal of Plant Physiology, 2018, 231: 415-433. DOI:10.1016/j.jplph.2018.10.022
[47]
Lehmann J, Kleber M. The contentious nature of soil organic matter. Nature, 2015, 528(7580): 60-68. DOI:10.1038/nature16069
[48]
Webb JR, Santos IR, Maher DT et al. The importance of aquatic carbon fluxes in net ecosystem carbon budgets: A catchment-scale review. Ecosystems, 2019, 22(3): 508-527. DOI:10.1007/s10021-018-0284-7
[49]
Bogard MJ, Vachon D, St-Gelais NF et al. Using oxygen stable isotopes to quantify ecosystem metabolism in northern lakes. Biogeochemistry, 2017, 133(3): 347-364. DOI:10.1007/s10533-017-0338-5
[50]
Cotrufo MF, Soong JL, Horton AJ et al. Formation of soil organic matter via biochemical and physical pathways of litter mass loss. Nature Geoscience, 2015, 8(10): 776-779. DOI:10.1038/ngeo2520
[51]
Zhang H, Goll DS, Wang YP et al. Microbial dynamics and soil physicochemical properties explain large-scale variations in soil organic carbon. Global Change Biology, 2020, 26(4): 2668-2685. DOI:10.1111/gcb.14994
[52]
Tank SE, Fellman JB, Hood E et al. Beyond respiration: Controls on lateral carbon fluxes across the terrestrial-aquatic interface. Limnology and Oceanography Letters, 2018, 3(3): 76-88. DOI:10.1002/lol2.10065
[53]
Marx A, Dusek J, Jankovec J et al. A review of CO2 and associated carbon dynamics in headwater streams: A global perspective. Reviews of Geophysics, 2017, 55(2): 560-585. DOI:10.1002/2016rg000547
[54]
Dorrepaal E, Toet S, van Logtestijn RSP et al. Carbon respiration from subsurface peat accelerated by climate warming in the subarctic. Nature, 2009, 460(7255): 616-619. DOI:10.1038/nature08216
[55]
Hao ST, Zhou AF, Zhang XN et al. Progress of research on the burial efficiency of organic carbon and its influencing factors in lacustrine sediments. Journal of Earth Environment, 2017, 8(4): 292-306. [郝盛吞, 周爱锋, 张晓楠等. 湖泊沉积有机碳埋藏效率及其影响要素研究进展. 地球环境学报, 2017, 8(4): 292-306. DOI:10.7515/JEE201704002]
[56]
Downing JA, Prairie YT, Cole JJ et al. The global abundance and size distribution of lakes, ponds, and impoundments. Limnology and Oceanography, 2006, 51(5): 2388-2397. DOI:10.4319/lo.2006.51.5.2388
[57]
Sobek S, Durisch-Kaiser E, Zurbrügg R et al. Organic carbon burial efficiency in lake sediments controlled by oxygen exposure time and sediment source. Limnology and Oceanography, 2009, 54(6): 2243-2254. DOI:10.4319/lo.2009.54.6.2243
[58]
Van Oost K, Quine TA, Govers G et al. The impact of agricultural soil erosion on the global carbon cycle. Science, 2007, 318(5850): 626-629. DOI:10.1126/science.1145724
[59]
Yu LJ, Cheng Y, Wang B et al. Climate and vegetation codetermine the increased carbon burial rates in Tibetan Plateau Lakes during the Holocene. Quaternary Science Reviews, 2023, 310: 108118. DOI:10.1016/j.quascirev.2023.108118
[60]
Xu H, Lan JH, Liu B et al. Modern carbon burial in Lake Qinghai, China. Applied Geochemistry, 2013, 39: 150-155. DOI:10.1016/j.apgeochem.2013.04.004
[61]
Kastowski M, Hinderer M, Vecsei A et al. Long-term carbon burial in European lakes: Analysis and estimate. Global Biogeochem Cycles, 2011, 25: GB3019. DOI:10.1029/2010gb003874
[62]
Anderson NJ, D'Andrea W, Fritz SC. Holocene carbon burial by lakes in SW Greenland. Global Change Biology, 2009, 15(11): 2590-2598. DOI:10.1111/j.1365-2486.2009.01942.x
[63]
Munroe J, Brencher Q. Holocene carbon burial in lakes of the Uinta mountains, Utah, USA. Quaternary, 2019, 2(1): 13. DOI:10.3390/quat2010013
[64]
Wang M, Chen H, Yu ZC et al. Carbon accumulation and sequestration of lakes in China during the Holocene. Global Change Biology, 2015, 21(12): 4436-4448. DOI:10.1111/gcb.13055
[65]
Jiang YJ, He W, Liu WX et al. The seasonal and spatial variations of phytoplankton community and their correlation with environmental factors in a large eutrophic Chinese Lake (Lake Chaohu). Ecological Indicators, 2014, 40: 58-67. DOI:10.1016/j.ecolind.2014.01.006
[66]
Li T, Bai Y, Chen XY et al. Longtime variation of phytoplankton in the South China Sea from the perspective of carbon fixation. In: Remote Sensing of the Ocean, Sea Ice, Coastal Waters, and Large Water Regions 2017. September 11-14, 2017. Warsaw, Poland. SPIE, 2017. DOI: 10.1117/12.2278038.
[67]
Siegel DA, Behrenfeld MJ, Maritorena S et al. Regional to global assessments of phytoplankton dynamics from the SeaWiFS mission. Remote Sensing of Environment, 2013, 135: 77-91. DOI:10.1016/j.rse.2013.03.025
[68]
Sanders R, Henson SA, Koski M et al. The biological carbon pump in the North Atlantic. Progress in Oceanography, 2014, 129: 200-218. DOI:10.1016/j.pocean.2014.05.005
[69]
Brothers S, Hilt S, Meyer S et al. Plant community structure determines primary productivity in shallow, eutrophic lakes. Freshwater Biology, 2013, 58: 2264-2276. DOI:10.1111/fwb.12207
[70]
Brothers S, Kazanjian G, Köhler J et al. Convective mixing and high littoral production established systematic errors in the diel oxygen curves of a shallow, eutrophic lake. Limnology and Oceanography: Methods, 2017, 15(5): 429-435. DOI:10.1002/lom3.10169
[71]
Kazanjian G, Flury S, Attermeyer K et al. Primary production in nutrient-rich kettle holes and consequences for nutrient and carbon cycling. Hydrobiologia, 2018, 806(1): 77-93. DOI:10.1007/s10750-017-3337-6
[72]
Deng WQ, Sun K, Jia JJ et al. Evolving phytoplankton primary productivity patterns in typical Tibetan Plateau Lake systems and associated driving mechanisms since the 2000s. Remote Sensing Applications: Society and Environment, 2022, 28: 100825. DOI:10.1016/j.rsase.2022.100825
[73]
Buffam I, Turner MG, Desai AR et al. Integrating aquatic and terrestrial components to construct a complete carbon budget for a north temperate lake district. Global Change Biology, 2011, 17(2): 1193-1211. DOI:10.1111/j.1365-2486.2010.02313.x
[74]
Lehner B, Döll P. Development and validation of a global database of lakes, reservoirs and wetlands. Journal of Hydrology, 2004, 296(1/2/3/4): 1-22. DOI:10.1016/j.jhydrol.2004.03.028
[75]
McDonald CP, Rover JA, Stets EG et al. The regional abundance and size distribution of lakes and reservoirs in the United States and implications for estimates of global lake extent. Limnology and Oceanography, 2012, 57(2): 597-606. DOI:10.4319/lo.2012.57.2.0597
[76]
Alin SR, de Fátima F L, Rasera M, Salimon CI et al. Physical controls on carbon dioxide transfer velocity and flux in low-gradient river systems and implications for regional carbon budgets. Journal of Geophysical Research, 2011, 116(G1): G01009. DOI:10.1029/2010jg001398
[77]
Butman D, Raymond PA. Significant efflux of carbon dioxide from streams and rivers in the United States. Nature Geoscience, 2011, 4(12): 839-842. DOI:10.1038/ngeo1294
[78]
Frankignoulle M, Abril G, Borges A et al. Carbon dioxide emission from European Estuaries. Science, 1998, 282(5388): 434-436. DOI:10.1126/science.282.5388.434
[79]
Sabine CL. The global carbon cycle: Integrating humans. Scope, 2004, 62: 17-44.
[80]
Maberly SC, Barker PA, Stott AW et al. Catchment productivity controls CO2 emissions from lakes. Nature Climate Change, 2013, 3(4): 391-394. DOI:10.1038/nclimate1748
[81]
Jia JJ, Sun K, Lü SD et al. Determining whether Qinghai-Tibet Plateau waterbodies have acted like carbon sinks or sources over the past 20 years. Science Bulletin, 2022, 67(22): 2345-2357. DOI:10.1016/j.scib.2022.10.023
[82]
Yan FP, Sillanpää M, Kang SC et al. Lakes on the Tibetan Plateau as conduits of greenhouse gases to the atmosphere. Journal of Geophysical Research: Biogeosciences, 2018, 123(7): 2091-2103. DOI:10.1029/2017jg004379
[83]
Sasaki Y, Koba K, Yamamoto M et al. Biogeochemistry of nitrous oxide in Lake Kizaki, Japan, elucidated by nitrous oxide isotopomer analysis. Journal of Geophysical Research, 2011, 116(G4): G04030. DOI:10.1029/2010jg001589
[84]
Zhu D, Wu Y, Wu N et al. Nitrous oxide emission from infralittoral zone and pelagic zone in a shallow lake: Implications for whole lake flux estimation and lake restoration. Ecological Engineering, 2015, 82: 368-375. DOI:10.1016/j.ecoleng.2015.05.032
[85]
Yang L, Lu F, Zhou XP et al. Progress in the studies on the greenhouse gas emissions from reservoirs. Acta Ecologica Sinica, 2014, 34(4): 204-212. DOI:10.1016/j.chnaes.2013.05.011
[86]
Schrier-Uijl AP, Veraart AJ, Leffelaar PA et al. Release of CO2 and CH4 from lakes and drainage ditches in temperate wetlands. Biogeochemistry, 2011, 102(1): 265-279. DOI:10.1007/s10533-010-9440-7
[87]
Song KS, Wang M, Du JA et al. Spatiotemporal variations of lake surface temperature across the Tibetan Plateau using MODIS LST product. Remote sensing, 2016, 8(10): 854. DOI:10.3390/rs8100854
[88]
Zhang YY, Zhou YJ, Shao QX et al. Diffuse nutrient losses and the impact factors determining their regional differences in four catchments from North to South China. Journal of Hydrology, 2016, 543: 577-594. DOI:10.1016/j.jhydrol.2016.10.031
[89]
Vachon D, Sponseller RA, Karlsson J. Integrating carbon emission, accumulation and transport in inland waters to understand their role in the global carbon cycle. Global Change Biology, 2021, 27(4): 719-727. DOI:10.1111/gcb.15448
[90]
Nakayama T. Inter-annual simulation of global carbon cycle variations in a terrestrial-aquatic continuum. Hydrological Processes, 2020, 34(3): 662-678. DOI:10.1002/hyp.13616
[91]
Shi GH, Peng CH, Wang M et al. The spatial and temporal distribution of dissolved organic carbon exported from three Chinese Rivers to the China Sea. PLoS One, 2016, 11(10): e0165039. DOI:10.1371/journal.pone.0165039
[92]
Bastviken D, Sundgren I, Natchimuthu S et al. Technical Note: Cost-efficient approaches to measure carbon dioxide (CO2) fluxes and concentrations in terrestrial and aquatic environments using mini loggers. Biogeosciences, 2015, 12(12): 3849-3859. DOI:10.5194/bg-12-3849-2015
[93]
Lauerwald R, Laruelle GG, Hartmann J et al. Spatial patterns in CO2 evasion from the global river network. Global Biogeochemical Cycles, 2015, 29(5): 534-554. DOI:10.1002/2014gb004941
[94]
Deng YB, Zhang YL, Li DP et al. Temporal and spatial dynamics of phytoplankton primary production in Lake Taihu derived from MODIS data. Remote Sensing, 2017, 9(3): 195. DOI:10.3390/rs9030195
[95]
Ford WI, Fox JF. Isotope-based fluvial organic carbon (ISOFLOC) model: Model formulation, sensitivity, and evaluation. Water Resources Research, 2015, 51(6): 4046-4064. DOI:10.1002/2015wr016999