直接通量观测揭示冰岛冰川溪流的CO2汇：玄武岩风化与水文混合的主控作用

《Journal of Geophysical Research: Biogeosciences》：Direct Flux Measurements Reveal CO2 Sinks in Icelandic Glacial Streams Driven by Weathering and Hydrological Mixing

　　

冰川补给河流日益被视为高纬度碳循环的重要组成部分，但对前冰川溪流中CO₂交换的直接观测仍然有限，尤其是在年轻的玄武岩流域。此处，研究人员量化了源自Fallj?kull的冰川补给溪流Virkisá沿3 km河段的CO₂交换空间与季节模式，并以冰岛?r?faj?kull冰帽排出的四条冰川溪流测量作为补充。跨四个季节，研究人员在Virkisá沿线六个站点进行了154次直接CO₂通量测量，并结合水化学数据、溶解有机碳（DOC）浓度和基于稳定同位素的水源解析。该溪流始终表现为净大气CO₂汇，冰川出口附近吸收最强，向下游逐渐减弱。pH和距冰川距离是CO₂通量最强的预测因子。冰川出口附近升高的pH和CO₂吸收反映了新近侵蚀玄武岩沉积物的快速风化，生成CO₂不饱和（CO₂-undersaturated）融水，其中溶解CO₂主要转化为碳酸氢盐。随着水体逐渐与大气平衡，吸收向下游下降，但净吸收在整个河段持续存在。由冰川冰融水、雪融水、冰面径流和地下水混合塑造的融水化学特征进一步改变了水化学和CO₂交换。尽管流量组成存在季节变化，CO₂通量未显示显著季节或昼夜差异，尽管峰值吸收出现在春季融雪期间。因此，研究结果确定前冰川溪流为持久CO₂汇，并强调与玄武岩风化相关的水化学——而非生物活动——主要控制CO₂动态。

论文解读

研究背景与问题：冰川补给河流正被视为全球碳循环的动态组成部分，但其作为温室气体源或汇的作用仍受约束不足。内陆水体通常相对于大气CO₂过饱和，全球溪流和河流是大气CO₂的重要来源，近期估计为1.9—2.3 Pg C yr^?1，小型陡峭溪流具有较强气体交换效率。然而，冰川补给系统具有低温、高沉积物负荷、低植被覆盖和强水文变异性，可能偏离这一范式；加拿大、格陵兰、欧洲阿尔卑斯和青藏高原等地的研究显示部分冰川溪流可出现净CO₂吸收。空气—水CO₂交换受水—气CO₂浓度梯度和气体传输速度控制，后者受湍流和风速影响；吸收或逸散方向主要受CO₂浓度梯度控制，而该梯度受水文和生物过程改造。化学风化也是冰川系统CO₂通量的关键驱动：硅酸盐和碳酸盐风化消耗大气CO₂，硫化物氧化可在高径流期释放更多CO₂；冰川流域强烈物理侵蚀产生高反应性矿物表面，可维持快速化学反应并导致远离平衡状态。部分情况下，快速矿物溶解可驱动强CO₂不饱和并抵消微生物CO₂产生，尤其在可矿化有机碳浓度低时。因此，冰川流域究竟是大气CO₂净源还是净汇仍有争议。由于冰川系统采样在后勤和方法上的困难，直接原位通量测量仍然稀缺，许多研究依赖间接顶空方法且时间分辨率较粗，冰川CO₂汇的幅度、持久性和变异性仍约束不足；冰岛年轻玄武岩地形促进快速风化，但缺乏直接通量观测。研究人员因此沿Virkisá及?r?faj?kull冰帽排出的四条冰川溪流开展研究，以解析高纬度冰川补给河流系统CO₂交换的时空格局。

研究开展、结论与意义：研究人员在Virkisá沿冰川出口至下游3 km设置六个站点，跨四季开展154次直接CO₂通量测量，并补充前冰川湖湖岸及四条相邻冰川溪流测量；结合水化学、溶解有机碳（DOC）、稳定水同位素和机器学习增强的端元混合分析，量化水源贡献并识别控制因子。结果表明，Virkisá在整个研究河段持续为净大气CO₂汇，冰川出口附近吸收最强，向下游逐渐减弱但全程保持负通量；pH和距冰川距离是CO₂通量最强预测因子。新近侵蚀玄武岩沉积物的快速风化产生CO₂不饱和融水，使溶解CO₂主要转化为碳酸氢盐；下游随水体与大气逐渐平衡，吸收减弱。融水化学受冰川冰融水、雪融水、冰面径流和地下水混合进一步改造。尽管流量组成存在季节变化，CO₂通量无显著季节或昼夜差异，但春季融雪期出现峰值吸收。研究意义在于提供冰岛冰川溪流首批直接水—气CO₂通量观测，直接量化CO₂吸收幅度，并指出玄武岩风化相关水化学而非生物活动主要控制CO₂动态，为高纬度冰川区碳循环及冰川退缩情景下的预测提供经验基础。

关键技术方法（≤250字）：样本来源为冰岛东南部Virkisá（源出Fallj?kull和Virkisj?kull）及?r?faj?kull冰帽四条相邻冰川溪流，2023年3月至2025年8月覆盖四季，Virkisá六站点共154次测量。采用浮箱（floating chamber）连接CO₂传感器直接测水—气CO₂通量；同步测pH、电导率、温度、溶解氧，并分析溶解有机碳（DOC）、稳定水同位素（δ¹⁸O、δ²H）、阳离子（ICP-MS）和营养盐（CFA）。用机器学习（ML）增强的端元混合分析（end-member mixing analysis）估算冰融水、冰面水、降水、地下水与雪融水贡献，并以统计模型和广义可加模型（GAMs）分析通量与环境变量关系。

研究结果（原文结果部分无编号小标题，按原文内容归纳）：通过四季六站点直接通量测量发现，所有站点和整个研究期CO₂通量均为负，表明冰川末端下游至3 km范围持续净吸收大气CO₂。冰川出口春季出现最强瞬时吸收（?29.99 ± 17.93 mg C m^?2 h^?1）；年均最高吸收出现在前冰川湖下游的Site 2（?19.85 ± 11.00 mg C m^?2 h^?1），略高于冰川出口（?18.96 ± 13.59 mg C m^?2 h^?1）。随距冰川距离增加，CO₂通量逐渐趋近中性，Site 5和Site 6分别约为?4.48 ± 3.53和?4.85 ± 3.48 mg C m^?2 h^?1；前冰川湖平均通量为?3.42 ± 1.01 mg C m^?2 h^?1。沿程pH由出口8.2降至Site 6的7.18，平均水温由0.26°C升至1.41°C；电导率春季最高并总体向下游降低，溶解氧沿程变化较小。DOC在冰川出口最低（0.12 ± 0.06 mg L^?1），向下游增至3 km处0.23 ± 0.07 mg L^?1。通过通量—距离关系分析发现，春、夏、秋三季CO₂通量与距冰川距离呈显著线性关系，通量随距离增加而变得较不负，春季空间梯度最陡。通过水源贡献分析发现，春季雪融水占38%、冰融水占22%；夏季冰面水占39%、冰融水占31%；冬季雪融水占30%、冰融水占29%；但CO₂通量与降水、地下水或冰融水贡献无显著关系。通过季节和昼夜比较发现，出口最大吸收出现在融季初期，夏季最低（?11.63 ± 11.05 mg C m^?2 h^?1），冬季回升（?21.90 ± 8.65 mg C m^?2 h^?1）；多数站点无显著季节差异，仅Site 3夏季与冬季差异显著；所有站点未检测到显著昼夜差异。通过DOC、营养盐和水化学相关分析发现，DOC与CO₂通量无显著关系，NO₂^? + NO₃^?和NH₄⁺多低于0.015 mg L^?1检出限，磷酸盐在出口平均0.98 ± 0.87 mg L^?1但与通量无显著相关；水温和pH与CO₂通量呈强相关（r = 0.63，p < 0.001；r = ?0.62，p < 0.001），pH与电导率正相关（r = 0.58，p < 0.001），并与冰面水、地下水、雪融水贡献相关。通过主要离子比值分析发现，Ca/Na和Mg/Na分布指示以硅酸盐来源溶质为主、碳酸盐贡献较小，表明整个水文系统以硅酸盐风化为主。通过沿程关系分析发现，CO₂通量与pH呈明显反向关系：向下游通量变得较不负，pH趋近中性；CO₂通量与距离正相关（r = 0.47，p < 0.001），pH与距离负相关（r = ?0.55，p < 0.001）。通过四条相邻冰川溪流验证发现，同一冰帽排出的溪流均出现负CO₂通量，范围为Fjallsá的?28.7至Svinafellsá的?15.1 mg C m^?2 h^?1。

讨论部分总结：4.1 冰岛前冰川溪流作为CO₂汇。该研究提供Virkisá首批直接水—气CO₂通量野外观测，解析冰川末端下游3 km及四季变化；Virkisá在整个研究河段包括前冰川湖均持续吸收CO₂，吸收幅度向下游下降但无显著季节和昼夜变化。观测值与加拿大高北极和Rongbuk冰川径流等报道的吸收重叠，幅度高于部分冰岛溪流、水库及大型冰岛河流的风化CO₂消耗估计；四条相邻冰川溪流的负通量支持该汇行为可能具有区域代表性。下游下降与初始CO₂不饱和水体逐渐与大气平衡、反应性矿物表面消耗及有机质输入增加等有关；pH与风化强度共同指示CO₂动态。春季纵向差异最强，可能与冬季积累的溶质和沉积物在融季初期释放、电导率升高及雪融水贡献与pH正相关有关；随融季推进，冰面水和冰融水稀释溶质与碱度；冬季低反应速率、低流量和部分冰盖可能限制气体交换和化学过程。4.2 时间动态有限。最高吸收出现在冰川出口春季融季初期，向下游季节差异减弱且全河段无昼夜差异，说明时间变异主要局限于冰川出口；浮箱测量已涵盖 prevailing 条件下流速和湍流对气体交换的影响。低DOC和低营养盐、无显著昼夜差异提示生物活动不是调查期CO₂交换的主导控制，化学自养过程可能参与但低DOC限制微生物代谢，物理化学过程更占主导。4.3 潜在物理化学机制。pH和距冰川距离是CO₂通量主要预测因子，pH更可能是系统地球化学状态指标而非直接驱动因子；玄武岩流域硅酸盐风化因新鲜矿物高反应性而高效，高碱度和高pH水体常使溶解CO₂低于大气平衡，出口pH约8.2利于CO₂转化为碳酸氢盐、降低CO₂分压并维持大气吸收；邻近Svínafellsj?kull冰下沉积物中硅酸盐和富铁矿物支持硅酸盐风化重要性。水文混合通过冰融水、雪融水、冰面水和地下水影响pH和水化学；冰与雪同位素特征相似限制端元区分，需谨慎解释混合结果，但混合分析与水文化学结果相互支持。前冰川湖未消除下游风化信号且湖面通量亦为负，Site 2吸收最高可能与高流量时富沉积物融水快速通过湖泊、反应性颗粒保持悬浮并继续风化有关，也可能与地下水或冰碛泉输入及冰下含水层和地热基融贡献有关。4.4 气候变化意义。升温、冰川退缩和高山景观绿化将通过影响呼吸、初级生产和矿物风化改变碳动态，并改变水文、沉积物供给和风化强度，进而影响pH、碱度和CO₂交换；冰岛冰川径流预计至少到本世纪中叶增加，前冰川溪流未来可能由大气CO₂汇转为源，仍需进一步研究，因为不同冰川区CO₂吸收机制可能不同。

研究结论翻译：研究人员的测量显示，Virkisá在四个季节以及沿冰川补给溪流数公里范围内始终是大气CO₂汇。来自同一冰帽排出的另外四条冰川补给溪流中测得的负通量表明，这种行为并不局限于Virkisá，尽管这些溪流仅在夏季采样。Virkisá沿程CO₂吸收的显著下降，加上有限的季节和昼夜变化，表明空间位置比采样时间更能一致地解释通量变异。尽管潜在机制无法得到确定性解析，但与pH和距离的观测关联将碳酸盐系统动态、大气平衡和风化识别为合理解释。直接、空间分辨的测量扩展了评估高纬度冰川补给水体CO₂交换的有限经验基础。