青藏高原三江源草地碳汇对极端气候指数的响应量化：牧业农业对生态系统呼吸及碳管理的影响

《Agriculture》：Quantification of Extreme Climate Index Contributions to Grassland Carbon Sink of Sanjiangyuan, Tibetan Plateau: Effects of Pastoral Agriculture on Ecosystem Respiration and Carbon Management Implications

　　

青藏高原草地生态系统在全球碳循环中发挥着关键作用，然而极端气候事件对其碳汇动态定量影响的认识仍不充分。本研究以典型高寒牧区——三江源（Three-River-Source）地区为对象，采用净生态系统生产力（Net Ecosystem Productivity, NEP）估算、Theil–Sen趋势分析、变异系数（Coefficient of Variation, CV）、Hurst指数及岭回归（Ridge Regression）模型，量化了2000–2024年草地碳源/汇动态的时空特征及13项极端气候指数对其的贡献。结果表明：区域平均NEP以1.49 g C m?2a?1速率增加，90.78%的区域表现为碳汇，反映生态系统呼吸（Ecosystem Respiration, Reco）相对较弱且植被碳吸收占主导。但Hurst指数分析显示85.68%的区域呈反持续性（anti-persistence），暗示未来可能由碳汇转变为碳源，即生态系统呼吸增强且碳汇功能可能被碳源功能取代。岭回归分析表明，极端温度指数特别是暖昼日数（TX90p，相对贡献38.1%）、冷夜日数（TN10p）和暖期持续指数（Warm Spell Duration Index, WSDI）是NEP变化的主导驱动因子。基于此提出适应性管理策略：在变异性高且具反持续性的区域，应通过调控放牧强度、优化围栏管理及修复退化草地，减缓过度呼吸相关碳排放并维持碳汇稳定性，为高寒区日益加剧的气候极端事件下优化牧业农业碳管理及调节生态系统呼吸提供科学依据。

论文解读：青藏高原三江源草地碳汇对极端气候指数贡献的量化及牧业农业对生态系统呼吸与碳管理的启示

研究背景与立题依据

草地生态系统作为陆地生态系统的重要组成部分，其碳汇功能与生态系统呼吸（Ecosystem Respiration, Reco）紧密关联，而气候变化尤其是极端气候事件正显著改变草地碳循环过程。已有研究表明青藏高原草地净生态系统生产力（Net Ecosystem Productivity, NEP）总体呈东高西低、以碳汇为主，气候变暖和湿润有利于NEP提升，不同草地类型对增温响应存在差异，三江源地区气温显著升高已影响碳汇动态；但既有研究普遍忽视极端气候事件（极端干旱、极端低温、极端降水等）对草地生态系统碳循环的定量影响，相关量化研究十分匮乏。青藏高原作为全球气候变化敏感区和"放大器"，三江源作为典型高寒牧业农业区，其草地具有高土壤有机碳含量但对气候极端事件敏感，碳汇动态直接关系到区域碳平衡与可持续牧业发展。因此，量化极端气候指数对三江源草地NEP变化的贡献，对理解碳源/汇转换机制及制定适应性碳管理策略具有重要科学意义。本研究发表于《Agriculture》期刊。

主要技术方法概述

研究人员以2000–2024年三江源地区为研究区，整合区内23个气象站观测数据、ERA5-Land再分析逐日温湿降水数据（空间分辨率0.1°×0.1°）及NASA MODIS MOD17A3HGF年净初级生产力（Net Primary Productivity, NPP）产品（空间分辨率500 m），基于Pei等建立的气候驱动模型由NPP估算土壤异养呼吸（Heterotrophic Respiration, R_h）后计算NEP＝NPP－R_h（NEP＞0为碳汇，NEP＜0为碳源）。采用Theil–Sen Median斜率法分析NEP年际变化趋势，计算NEP像元尺度变异系数（Coefficient of Variation, CV）划分稳定等级，通过重标极差分析（Rescaled Range Analysis, R/S）求取Hurst指数判断未来NEP持续性（H＜0.5为反持续性，H＞0.5为持续性），结合二者划分NEP未来变化类型。基于逐日温湿数据计算13项极端气候指数（含5项极端温度指数：暖昼日数TX90p、冷夜日数TN10p、暖期持续指数WSDI、冷期持续指数CSDI、结冰日数ID0；8项极端降水指数：连续干旱日数CDD、连续湿润日数CWD、最大1日降水量RX1day、最大5日降水量RX5day、极湿日降水量R95p、特大湿日降水量R99p等），所有指数经极小-极大标准化后，采用岭回归（Ridge Regression）模型消除多重共线性，分解各极端气候指数对NEP变化趋势的相对贡献与绝对贡献。

研究结果

3.1 极端气候指数的分布特征

研究人员计算并分析了2000–2024年各极端温度与极端降水指数的年时序与空间趋势。极端温度指数整体呈暖化倾向——TX90p与WSDI增加，TN10p、CSDI及ID0减少；空间上TX90p全区普遍增多，WSDI增幅集中于东南部澜沧江源区，冷相关指数减幅在中东部及黄河源、通天河流域较明显。极端降水指数年时际波动大、拟合优度偏低——连续干旱日数CDD略降，连续湿润日数CWD略升，RX1day、RX5day、R95p、R99p均呈增加趋势（RX5day增幅强于RX1day），空间上RX系列与R95p、R99p增量集中于东南部，CDD多数区域略增且西南部更明显，CWD增幅最大处位于通天河流域。

3.2 草地生态系统碳源/汇的时空分布

三江源草地NEP多年均值区域平均为130.16 g C m^?2a^?1，空间上东南高西北低；碳汇面积占总草地面积90.78%，碳源占9.22%。年际上区域年均NEP由2000年111.35 g C m^?2a^?1增至2024年165.22 g C m^?2a^?1，以1.49 g C m^?2a^?1速率上升，说明2000–2024年研究区整体以碳汇为主导且碳吸收能力波动增强。

3.3 草地生态系统碳汇功能的变异性

NEP的变异系数（CV）均值为0.22，空间上自东南向西北稳定性递减。稳定区（CV最小）分布于东部，较稳定区（占比约50%）分布于中部高寒草甸草原，不稳定过渡区分布于长江源东部及最北部，高度不稳定区分布于西部羌塘高寒荒漠草原。表明碳汇功能总体稳定但局部受水热波动及放牧/围封等人畜干扰而易失稳。

3.4 草地生态系统碳汇功能的持续性特征

Hurst指数介于0.10–0.86（均值0.39），表明NEP未来变化以强反持续性为主。"增长转下降"（Sen斜率＞0且H＜0.5）型占全区85.68%，预示多数区域将由当前碳汇主导转为潜在碳源主导，主因可能为2020年后增温加速冻土退化致土壤异养呼吸超过NPP；"持续增加"（Sen＞0且H＞0.5）型占10.77%，集中分布于黄河源东北部，受益于适度放牧抑制了呼吸波动；其余为零星散布的持续减少与"下降转增加"型（几无分布）。

3.5 极端气候指数对草地生态系统碳汇的贡献分析

岭回归显示极端温度指数累积相对贡献55.5%，极端降水指数44.5%。单项中暖昼日数（TX90p）相对贡献最高（38.1%），其次为冷夜日数（TN10p）与暖期持续指数（WSDI），最大1日降水量（RX1day）贡献最低（3.9%）。空间主导因子分布：TX90p为西北部主要影响因子（高温促进土壤微生物活性、增强Reco从而压低NEP），TN10p为黄河源及周边东部局部主导（冷夜减少延长呼吸时段、抬高根呼吸），WSDI主导南部，ID0显著影响西部，R95p为南部部分区主导降水因子。绝对贡献上TX90p正值最大（1.63 g C m^?2a^?1），CWD为降水因子中最大正值（0.22 g C m^?2a^?1），ID0呈负贡献；极端温度绝对贡献集中在东南部，极端降水绝对贡献集中在东北部。

3.6 草地农业管理模式与碳汇动态的关联分析

适度放牧可通过促进植被更新与养分循环维持植被生产力与生态系统呼吸间的平衡，使碳汇功能较稳定；禁牧初期植被固碳增强，但长期围封可致碳汇饱和甚至下降；高强度放牧压制光合、降低盖度并通过牲畜践踏破坏土壤结构、增强土壤呼吸，削弱碳汇。合理人为干扰应控制在生态系统自我调节弹性范围内。

讨论与结论总结

讨论部分指出：NEP变化是植被碳吸收与生态系统呼吸释放平衡的结果；极端温度指数（TX90p、TN10p、WSDI）通过影响高寒草地物候、生产力和土壤微生物呼吸主导NEP变异——适度增温可延长生长季促进固碳，但极端高温强化土壤呼吸、加速有机质分解，若植被生产力增幅不足以抵消呼吸增强则NEP下降；极端降水通过调节土壤水分与微生物活性影响NEP，但因年际波动大其线性拟合弱于温度因子。碳管理关键不仅是提升固碳量，更需在脆弱区维持现有碳汇稳定，过高放牧强度是重要人为碳损失因子，需控制草畜平衡与土壤扰动；围封管理应权衡恢复效益与生产需求，避免单一模式。建议按NEP变异性与Hurst持续性实施分区管理：高稳区采用围封+适度放牧混合模式；高变或反持续脆弱区严格控牧、修复退化草地并建立极端气候—牧业适应联动预警机制。

结论翻译如下：

(1) 2000–2024年三江源草地生态系统总体表现为碳汇。区域平均NEP以1.49 g C m^?2a^?1速率增加，碳汇面积占比90.78%；空间上东南高西北低。NEP持续增长表明总初级生产力（Gross Primary Productivity, GPP）增速超过生态系统呼吸（Reco），植被碳吸收增强主导碳汇变化。

(2) NEP稳定性存在显著空间分异，变异系数均值为0.22，区域碳汇功能总体稳定但局部不稳。Hurst指数均值0.39，"增长转下降"型区域占85.68%，提示当前NEP上升趋势在许多区域难以持续。

(3) 岭回归结果显示极端温度指数对NEP变化贡献大于降水相关指数（温度类累积相对贡献55.5% vs 降水类44.5%）；暖昼日数TX90p平均相对贡献最高（38.1%），其次为冷夜日数TN10p与暖期持续指数WSDI。极端温度事件加剧生态系统呼吸（Reco），危及净初级生产力（NEP）的长期稳定性。

(4) 三江源草地碳管理应考量空间异质性与对极端气候的敏感性，NEP变异性高或具反持续性趋势的区域在未来适应性管理中需重点关注，尤其在温度极端事件频发背景下。