青藏高原多年冻土碳释放临界点研究学术报告

本研究由中国科学院青藏高原研究所“青藏高原地球系统与资源环境”国家重点实验室（State Key Laboratory of Tibetan Plateau Earth System, Environment and Resources, Institute of Tibetan Plateau Research, Chinese Academy of Sciences）的Yuxi Wei、Juan Li、Jinzhi Ding等作者团队主导完成，合作单位包括兰州大学、北京大学等。该研究成果已于2026年发表在顶级学术期刊《自然·通讯》（Nature Communications）上，论文标题为《Permaffrost tipping point triggered by warming-driven loss of old carbon》。

一、研究学术背景

本研究属于全球变化生态学和生物地球化学循环交叉领域，核心关注点是多年冻土碳库的气候反馈机制。多年冻土储存了全球巨量的有机碳，其稳定性被认为是地球系统中关键的“气候引爆点”之一。气候变暖导致冻土融化，可能释放被封存数千年的“古老碳”，形成温室气体排放→加剧变暖→更多冻土融化的正反馈循环，从而严重放大气候变化。然而，关于多年冻土碳库在何种增温阈值下会跨越临界点、从碳汇转变为强劲的碳源，科学界仍存在重大不确定性。这主要源于三个关键科学问题：1）生态系统总初级生产力（Gross Primary Productivity, GPP）与生态系统呼吸（Ecosystem Respiration, Reco）对增温的响应轨迹是线性还是非线性？是否存在温度阈值？2）在增温下，生态系统呼吸（Reco）的增加多大程度上来源于深层古老碳的矿化释放？3）GPP与Reco的解耦（即光合固碳无法抵消呼吸释放）是否以及在多高的温升水平下发生，并驱动系统性的碳源转换？

针对这些问题，尽管已有一些高纬度冻土的研究，但对于全球最大的高山多年冻土区——青藏高原，仍缺乏长期、多梯度增温控制实验与深度分辨的碳通量观测相结合的实证研究。青藏高原储存了约47.36 Pg C的土壤碳，其变暖速率快于全球平均水平，活动层持续加深，使其碳库极为脆弱。因此，本研究旨在通过一个为期五年、包含四个增温梯度的野外控制实验平台，结合超过4万小时的近地表通量观测和土壤剖面CO2浓度及δ13C-CO2同位素剖面监测，定量揭示青藏高原多年冻土生态系统碳通量对渐进式增温的非线性响应，识别碳循环关键过程的温度阈值，并阐明古老碳释放在其中扮演的核心角色。

二、详细研究流程

本研究采用野外控制实验、高频连续监测、生物地球化学分析和数据建模相结合的综合方法，流程严谨且系统。

1. 实验平台建立与处理设置：研究地点位于西藏安多县（海拔4790米）的典型高山草甸多年冻土生态系统。研究团队构建了一个多梯度增温实验平台，采用精确反馈控制的红外加热系统，设置了四个处理：对照（环境温度）、+1°C、+2°C和+4°C目标增温。每个处理设4个重复样方，共16个样方，采用随机区组设计。增温系统通过实时监测空气温度并调整加热输出，确保增温目标的精确实现。实测土壤温度增温幅度随深度衰减，但在所有层次（10厘米至160厘米）均达到统计显著水平（例如，10厘米深度分别增温约+1.3°C, +1.6°C, +3.0°C；160厘米深度分别增温约+0.3°C, +0.4°C, +0.9°C）。

2. 生态系统碳通量高频监测：每个样方配备一个自动开闭透明箱室系统，与Picarro G2301气体分析仪相连。在2020年至2024年的五年间，系统以每1.5小时一次的频率自动测量净生态系统CO2交换量（Net Ecosystem Exchange, NEE），累积超过4万次有效观测。夜间NEE即为Reco。基于夜间Reco与土壤温湿度的关系，利用随机森林回归模型估算白天的Reco，进而计算出GPP（GPP = Reco - NEE）。该随机森林模型利用10厘米和20厘米深度的土壤温度和体积含水量作为预测因子，表现出高预测精度（R² 0.79-0.92），并通过随机掩蔽20%数据进行了严格验证。

3. 土壤剖面气体与同位素深度分辨采样：团队建立了永久性的深度特异性土壤气体剖面系统，在10、20、80、120和160厘米深度安装定制防水扩散室。在生长季每十天、非生长季每月采集各深度气体样品。使用Picarro G2201-i分析仪测定CO2浓度和稳定碳同位素比值（δ13C-CO2）。同时，结合传统的物理根际分离法（如环割、壕沟），在样方内设置不同处理的PVC collar（代表Reco、土壤呼吸Rs、异养呼吸Rh），定期用便携式土壤通量系统测量其通量和δ13C-CO2，为模型结果提供独立验证。

4. 土壤与微生物群落分析：

   • 土壤理化性质与碳组分：实验开始前（2019年）和中期（2025年8月）采集不同深度土样，分析pH、土壤有机碳、总氮、有效磷。通过加速器质谱进行放射性碳（14C）测年，确认深层（80-160厘米）碳年龄为1845-3411年，属于古老碳库。2025年样品进一步分析了溶解性有机碳、颗粒有机碳和矿物结合有机碳。
   • 微生物群落表征：通过磷脂脂肪酸分析和宏基因组测序揭示微生物群落结构和功能变化。PLFA分析用于量化细菌、真菌生物量，并计算真菌/细菌比、革兰氏阳性/阴性菌比等指标。宏基因组测序揭示了微生物分类群和功能基因对增温的响应。

5. 呼吸源解析与数据统计分析：

   • 同位素混合模型：利用不同呼吸源（地上植物呼吸Rag、浅层根呼吸、浅层异养呼吸）和不同土壤深度碳库的δ13C特征值，应用贝叶斯稳定同位素混合模型（MixSIAR）量化了深层土壤呼吸对总土壤呼吸的贡献比例。模型运行了长链Markov Chain Monte Carlo，并通过Gelman-Rubin诊断、Geweke诊断、轨迹图和自相关分析确保了收敛。
   • 统计分析：使用Kruskal-Wallis和Wilcoxon检验分析不同增温处理间通量、气体浓度等的差异。通过三种方法量化温度敏感性：1）单位温度变化的通量变化量；2）基于时间序列指数拟合的时间温度敏感性Q10；3）利用广义可加混合模型控制湿度影响后推导的动态Q10。

三、主要研究结果

1. 净生态系统CO2交换（NEE）对增温的响应：五年监测表明，研究地点在2020-2024年本身就是一个净CO2源（+171.82 g CO2 m-2 yr-1）。增温显著增强了净CO2释放。与对照相比，+1°C、+2°C和+4°C处理使年平均NEE（净释放）分别增加了44%、80%和176%。NEE的温度敏感性随增温梯度呈非线性放大，例如在控制湿度后，其Q10值从对照的0.26增至+1°C的0.39、+2°C的0.63，并在+4°C时飙升至1.95。这种非线性响应源于GPP和Reco的差异化响应。

2. GPP与Reco的响应与解耦：

   • 生态系统呼吸（Reco）：Reco在所有增温处理下均持续增加，年通量在+1°C、+2°C和+4°C下分别增加15%、37%和40%。其中，+2°C增温下的响应最强，表明存在非线性温度敏感性。
   • 总初级生产力（GPP）：GPP表现出先促进后抑制的“热最适”响应。+1°C和+2°C增温使生长季GPP分别增加3%和21%，但+4°C增温则使GPP比对照降低14%。其温度敏感性系数也呈现先升后降的模式。
   • 解耦与主导作用：在低至中度增温（°C）下，Reco的增量是GPP增量的1-16倍，Reco的变化主导了NEE变化的73-87%。在极端增温（+4°C）下，虽然Reco仍是NEE变化的主要贡献者（78%），但GPP受到的热抑制（下降14%）协同作用，极大地削弱了生态系统的碳抵消能力，加剧了净碳损失。特别值得注意的是，从+2°C到+4°C的转变过程中，GPP出现了急剧下降（|δGPP| = 147.19 g CO2 m-2 yr-1），而Reco仅有微弱增加（|δReco| = 17.14 g CO2 m-2 yr-1），GPP下降的幅度远超Reco的增加。这标志着一个气候引爆点（2-4°C增温区间内）的出现，生态系统从“呼吸主导的弱碳源”转变为“光合抑制与呼吸持续共同作用的强碳源”。

3. 呼吸来源的转变与古老碳的释放：

   • 剖面特征：土壤CO2浓度随深度增加，即使在有机碳含量极低（%表层）的160厘米深处，非生长季浓度仍可达表层的3-5倍。深层δ13C-CO2值全年保持相对稳定且偏负（约-18.34‰至-18.01‰），与古老碳分解的信号一致。氧气剖面分析排除了全年大部分时间内扩散受限导致CO2积累的主因，支持了深层存在持续生物生产的结论。
   • 来源解析：MixSIAR模型和同位素证据共同表明，增温显著增强了生长季深层土壤碳对土壤呼吸的贡献。在对照条件下，深层呼吸贡献占生长季土壤呼吸的67%，而在增温条件下，这一比例上升至76%。这意味着，增温导致的额外CO2释放，主要来源于生长季对深层、古老有机碳库的活化分解。
   • 微生物与碳库机制：宏基因组和PLFA分析显示，在+4°C极端增温下，微生物群落向寡营养型、资源保守型策略转变（如Candidatus_Dormibacterota和Verrucomicrobiota富集），真菌/细菌比下降，表明碳有效性降低而非微生物周转加快。土壤碳组分分析发现，表层活性碳库（颗粒有机碳）和传统认为较稳定的碳库（矿物结合有机碳）均显著减少，深层溶解性有机碳也下降26%。这证实了强烈增温首先消耗活性碳，进而驱动微生物群落演替以利用更顽固的古老碳库，从而维持较高的Reco水平。这种深层碳库的“代谢阀门”机制在+4°C增温下被显著放大。

四、研究结论与意义

本研究通过五年多梯度增温实验，首次在青藏高原高山多年冻土生态系统中直接观测并证实了由变暖驱动的古老碳释放所触发的碳循环临界点。主要结论如下： 1. 存在一个2-4°C的增温临界区间：在此区间内，生态系统碳平衡发生质变。当增温超过光合作用的温度阈值时，GPP开始下降，而Reco则通过活化深层古老碳库得以维持甚至缓慢增加，导致GPP与Reco严重解耦。 2. 深层古老碳释放是关键机制：增温不仅加速了表层碳循环，更重要的是打开了深层、古老碳库的“代谢阀门”。在持续+4°C增温下，深层古老碳的分解贡献了约76%的土壤呼吸，相当于约59%的生态系统呼吸。这挑战了关于深层土壤碳稳定性的传统认知。 3. 定量预估碳释放潜力：基于实验结果外推，到本世纪末（2081-2100），青藏高原多年冻土区预估的增温（2.69°C）可能导致古老碳释放增加24-47 g CO2 m-2 yr-1。在高纬度碳更丰富的多年冻土区，等效变暖可能引发更高的碳排放。

本研究的科学价值在于：为多年冻土碳-气候正反馈理论提供了关键的、基于控制实验的实证证据；首次定量给出了高山多年冻土生态系统碳通量关键过程的非线性响应阈值；强调了在评估冻土碳反馈时必须考虑深度分辨的碳动态和古老碳的动员过程。其实践价值在于：研究结果对改进地球系统模型至关重要，模型需要纳入这些非线性阈值行为和深层碳动态，才能更准确预测未来气候。此外，研究警示，即使实现《巴黎协定》2°C温控目标的下限，部分地区（如快速变暖的青藏高原）的多年冻土生态系统也可能已接近或跨越碳源转换的临界点，对全球碳预算和气候治理目标构成严峻挑战。

五、研究亮点

1. 创新的实验设计：建立了全球首个在高山多年冻土区结合多梯度精确增温控制、高频通量监测和土壤剖面同位素深度解析的综合性实验平台，填补了方法学空白。
2. 揭示关键阈值与机制：不仅观测到通量的非线性响应，更重要的是揭示了“光合热抑制”与“古老碳活化释放”协同驱动碳源转换的双重机制，并首次在野外实验中明确指出了2-4°C的潜在引爆点区间。
3. 多学科证据链完整：融合了生态系统生理学、稳定同位素地球化学、微生物生态学和生物地球化学等多学科手段，从通量、同位素信号、微生物群落结构到土壤碳库变化，形成了相互印证、逻辑严密的完整证据链，强有力地支持了“古老碳释放触发临界点”的核心假设。
4. 重要的区域与全球意义：聚焦于对气候变化高度敏感且碳储量巨大的青藏高原，其研究结论对于理解全球高山及高纬度多年冻土碳库的脆弱性具有重要的参考价值和预警意义。

六、其他有价值内容

研究也坦陈了其局限性并指出了未来方向：1）五年实验周期虽能揭示清晰的转变趋势，但不足以评估长期（如几十年尺度）的适应性响应（如植被群落演替、微生物功能适应）。2）离散的增温处理（+1, +2, +4°C）只能界定临界点的大致区间（2-4°C），而非精确温度值，未来需要更精细的温度增量实验来进一步约束。3）未来研究需结合长期监测、植被动态和微生物基因组分析，以检验潜在的适应性过程是否会改变所观测到的碳损失的持久性。这些思考体现了研究的严谨性和前瞻性。