气候对阿尔卑斯山草地生产力的直接与性状介导效应：Sentinel-2时间序列下的冠层功能属性、样地内光谱功能异质性与NIRv

《Global Change Biology》：Direct and Trait-Mediated Effects of Climate on Grassland Productivity Across the Alps

　　

理解气候如何通过能量限制和植被表达的功能属性变化塑造生态系统生产力仍具挑战性。这对山地草地中基于性状的生物多样性-生态系统功能研究至关重要，因为变暖与环境梯度相互作用。在此，研究人员区分欧洲阿尔卑斯山范围内气候对生产力的直接控制与通过冠层功能组成、生理状态及样地内光谱功能异质性起作用的间接路径。利用Sentinel-2时间序列（2017–2024），研究人员从色素投资和水分状态的光谱代理量化冠层功能结构，使用功能丰富度（functional richness，FR_ic）指标量化样地内光谱功能异质性，并使用植被近红外反射率（near-infrared reflectance of vegetation，NIR_v）量化植被生产力。研究人员结合面板数据模型与分段结构方程模型，将气候关联分解为直接成分和通过冠层功能属性起作用的性状介导成分，并使用变系数模型评估这些路径如何沿环境梯度变化。研究人员识别出两条连接气候与生产力的路径。更温暖的生长季直接提高生产力，但与较低的叶绿素与类胡萝卜素色素平衡相关，符合相对更高的类胡萝卜素投资，从而产生相反的间接路径。更大的冬季积雪主要通过直接路径降低生产力，而其冠层介导路径不确定。气候水量亏缺降低冠层水分状态和样地内光谱功能异质性，但其与生产力的直接和总关联与零无法区分。环境背景调节气候-植被关系。基线温度和水分可获得性强烈塑造气候对冠层性状和生产力的影响，而冠层属性-生产力关系在空间上变化很小。变暖对生产力的正向响应在更湿润和季节性变化更强的气候背景中更强，但在基线冬季积雪积累更大的地方减弱。总之，研究结果表明，山地草地短期生产力响应反映直接能量限制与通过整合冠层功能属性的间接路径之间的平衡，且这种平衡沿大气候梯度变化。

该论文发表在《Global Change Biology》，围绕气候如何通过能量限制与植被功能属性变化调控生态系统生产力这一核心问题。研究背景在于，气候通过能量、水分和季节性设定生物活动边界，但生态系统响应难以仅由直接气候效应解释；生物多样性可通过物种组成、结构和功能性状调节、缓冲或放大气候效应。山地草地沿陡峭环境梯度可从能量限制转向水分限制，气候变化既可直接通过生长季长度或雪动态影响生产力，也可间接通过性状组成和功能多样性改变资源利用。然而，直接气候控制与通过植被功能属性间接路径的相对贡献仍未解决，且这种平衡如何跨环境梯度变化尤其不确定。已有生物多样性-生态系统功能（biodiversity–ecosystem functioning，BEF）研究缺乏在自然系统中空间显式估计短期生产力变异中直接气候与间接功能属性的比例。研究人员因此整合遥感时间序列、面板数据和中介模型，在欧洲阿尔卑斯山天然草地中区分直接气候路径与通过冠层功能属性的间接路径。

关键技术方法不超过250字：研究队列为欧洲阿尔卑斯海拔>1000 m天然草地（CORINE Land Cover 2018 Class 3.2.1），基于Copernicus GLO30 DEM和阿尔卑斯弧多边形，随机抽取10,000个60 m半径圆形样地。利用Sentinel-2表面反射影像（2017–2024，10–20 m），限定生长季4–10月并保留分数植被覆盖度（fractional vegetation cover，FVC）≥0.60的像元；提取红边叶绿素指数（red-edge chlorophyll index，CI_re）、叶绿素-类胡萝卜素指数（chlorophyll–carotenoid index，CCI）和归一化差异红外指数（normalized difference infrared index，NDII）的80th percentile，计算FR_ic及生长季累计NIR_v。气候变量包括生长度日（growing degree days，GDD；ERA5-Land约9 km）、气候水量亏缺（climatic water deficit，CWD）和雪水当量（snow water equivalent，SWE；TerraClimate约4 km）。统计采用逆双曲正弦转换、线性混合模型、双向固定效应、分段结构方程模型、变系数混合模型，并以150 km空间块自助法1000次推断。

结果部分保留原文小标题并浓缩如下。3.1 气候异常对草地生产力的直接与性状介导效应：通过样地内去均值线性混合模型和150 km空间块自助法，H1表明GDD升高降低CCI，CWD升高降低NDII和FR_ic，其余气候-冠层关联不确定；H2表明NIR_v随CI_re和CCI升高而增加，随NDII和FR_ic升高而下降；H3在纳入冠层路径后表明NIR_v随GDD升高而增加、随SWE升高而下降，CWD直接效应区间过零，CI_re和CCI仍正关联NIR_v，NDII和FR_ic仍负关联NIR_v。3.1.1 将路径整合：结构方程模型：分段结构方程模型识别出两条并行路径，即通过冠层光谱性状的间接路径与气候异常到生产力的直接路径。GDD间接路径主要经CCI，直接正效应更大，总效应为正；SWE直接和总效应为负，冠层介导路径不支持；CWD经NDII和FR_ic的间接点估计为正，但仅FR_ic介导的间接效应95%自助区间排除零，其余包括总效应不确定。3.2 环境调节与效应的空间变异：变系数混合模型显示环境背景调节气候-冠层和气候-生产力路径，最强调节出现在气候对冠层属性的H1；基线气候变量产生最广泛变异，地形也强烈调节部分路径，土壤效应一般较小。H1中GDD→CCI和CWD→NDII调节最强，GDD→CCI负关系随年均温和海拔升高而减弱、随温度季节性增强而增强；CWD→NDII负关系在长期水亏缺和温度季节性更大处增强、在更暖处减弱；CWD→FR_ic调节较弱。H2的环境调节总体较弱。H3中基线气候主导GDD→NIR_v，其正响应随长期降水和温度季节性增强而增强，随基线SWE增大而减弱；CWD→NIR_v亦受基线气候调节，尽管平均效应不确定。3.2.1 关注GDD→CCI的环境调节：基线环境解释19.7%的样地间变异，年均温、海拔和温度季节性是最清晰调节因子，空间投影为描述性。3.2.2 关注GDD→NIR_v的环境调节：静态环境解释约30%的样地间变异，支持年均降水、温度季节性和基线SWE三个调节因子，其余不确定，空间投影为描述性。

讨论部分总结：4.1 气候对生产力的直接与间接效应指出，暖生长季与较低CCI相关，符合更高的类胡萝卜素投资和光保护需求；CWD降低NDII和FR_ic，符合环境过滤和性状驱动理论预期；色素相关冠层属性与生产力正相关，符合质量比效应；NDII与生产力弱负相关，可能与能量-水共限有关；FR_ic负相关不能解释为个体或物种水平α功能多样性降低生产力，而是不利年份群落整合冠层状态空间异质性较高且生产力较低。直接GDD正效应和SWE负效应支持能量可用性与生长季长度限制；结构方程模型表明直接能量路径和间接冠层功能路径并存，有时方向相反。GDD经CCI的负间接效应未抵消直接正效应；CWD直接与性状介导路径相互抵消，总效应不确定。短期增温生产力增益可能以冠层功能完整性降低为代价。4.2 估计路径的环境调节强调，直接和间接气候效应具有背景依赖性；气候调节因子最广泛，地形调节部分路径，土壤效应较小；冠层属性-生产力关系弱调节，表明短期间优势性状机制可能压过背景依赖互补效应。GDD→NIR_v在长期降水多和温度季节性强处增强、在基线雪多减弱，说明水分可获得性允许暖年受益。4.3 视角与局限指出，估计路径方向的空间变异警告不应普适化BEF关系；需在其他山地系统重复；应扩展到稳定性、韧性和多功能性；研究时段较短且遗漏延迟群落和性状响应；气候效应可能非线性，需考虑阈值和饱和；遥感不确定性未完全传播；需链接光谱代理与物种性状数据库、地下过程及土地利用管理。4.4 结论译文：通过整合遥感与面板数据和中介模型，本研究提供了一个空间显式框架，用于将直接气候路径与通过冠层功能属性的间接路径分离。结果识别了短期生产力对变暖的响应与更高能量可用性相关的地点，以及冠层或水分限制可能限制这些响应的地点。这些条件性估计应在更长时间序列和其他山地系统中进行检验，之后才能推广到长期生态系统轨迹。