鸟类稳定同位素（Stable Isotopes）的空间格局揭示气候与土地利用变化的生态影响

《Integrative Zoology》：Spatial Patterns of Avian Stable Isotopes Reveal Ecological Impacts of Climate and Land-Use Change

　　

稳定同位素分析为追踪动物食性来源与营养级结构提供了综合性研究途径，同时对环境变化与土地利用变化具有高度敏感性，因此成为评估动物资源利用与环境变异相互作用的有力工具。然而，跨越大量物种的空间同位素格局及其对气候变化和人为活动的响应仍知之甚少。本研究应用贝叶斯混合效应模型（Bayesian mixed-effects modeling）框架和同位素生态位分析（isotopic niche analyses），利用来自中国西南地区（1958–2007年）766只个体、涵盖200种鸟类（Aves）的羽毛样品，评估生态、气候和人为因素如何共同塑造跨类群的羽毛碳（δ¹³C）和氮（δ¹⁵N）同位素空间格局。结果显示，δ¹³C和δ¹⁵N值在不同食性和迁徙型功能群之间差异显著，反映了物种间资源利用策略和营养级组织的明显分化。羽毛δ¹⁵N值随牧草地覆盖度增加而升高，表明放牧景观中基线δ¹⁵N可能存在抬升效应。δ¹⁵N还与降水量呈正相关趋势，指示较湿润环境可能与较高的氮同位素基线相关。此外，同位素生态位分析揭示农田栖息地的同位素生态位宽度大于森林栖息地，尤其对植食性鸟类而言，凸显了农业扩张对资源利用的影响。总体而言，本研究展示了生态特征、气候因素和人为干扰如何共同塑造跨类群和跨尺度的鸟类营养格局，为理解不同生态系统的营养结构及土地利用变化对生态系统功能的后果提供了新见解。

本文发表于《Integrative Zoology》，研究人员以中国四川省及其周边地区1958年至2007年间收集的博物馆鸟类标本羽毛为研究对象，系统探讨了气候变化与土地利用变化对鸟类稳定同位素空间格局的生态影响。

研究背景与现存问题方面，全球气候变化和土地利用模式的持续转变深刻改变着陆生生态系统的碳循环和氮循环。人口快速增长与经济发展加剧了 land conversion（土地转换）、农业扩张和工业污染，导致栖息地破碎化及生态系统结构与功能的重大变化。这些转变改变了初级生产者的营养来源和同位素组成，并通过食物网传递，最终在更高营养级得以记录。稳定同位素分析因其能够保存反映膳食输入和环境条件的生物化学特征，已成为生态学和环境科学中量化能量流动、资源利用、营养循环和营养相互作用的核心分析工具，其中碳（δ¹³C）和氮（δ¹⁵N）同位素尤为关键：δ¹³C主要指示基部资源的光合作用途径（C₃、C₄或景天酸代谢[CAM]植物）及栖息地间的碳源异质性，而δ¹⁵N沿食物链逐步富集，是营养级位置和氮循环强度的可靠指标。在鸟类生态学中，基于δ¹³C和δ¹⁵N空间的同位素生态位分析日益被用于量化鸟类群落的资源分配和功能分化。然而，现有关于气候和人为驱动因素对鸟类稳定同位素影响的研究大多局限于少数物种或特定栖息地类型，对于环境异质区域内多样鸟类群落的空间同位素格局及其潜在生态机制仍缺乏理解。这一空白在四川省尤为突出——该地区高度异质的景观和陡峭的环境梯度支撑了广泛的生态系统类型，气候与土地利用相关过程可能共同导致鸟类同位素格局的空间变异。

为解决上述问题，研究人员开展了本项研究，旨在通过多物种、多因素的整合视角，评估生态特征、自然环境梯度和人为土地利用模式如何共同塑造鸟类羽毛δ¹³C和δ¹⁵N值的变化规律。

关键技术方法方面，研究以博物馆标本为基础获取样本来源，涉及1958–2007年收集的中国西部师范大学生物博物馆鸟类羽毛标本，覆盖四川省及邻近地区，涵盖高原、盆地、山地等多类生境。研究团队采用了贝叶斯混合效应模型（使用MCMCglmm软件包实现）分析生态与环境因素对同位素值的影响，以科（Family）为随机效应；利用SIBER（Stable Isotope Bayesian Ellipses in R）框架进行同位素生态位宽度分析，以标准椭圆面积校正小样本量（SEAc）表征核心同位素生态位，并通过贝叶斯后验分布（SEA.B）评估生态位宽度的不确定性；同时使用Z-score标准化处理连续预测变量，并对收集年份进行Suess效应校正以消除大气CO₂中¹³C含量持续下降的影响，还通过SIDER包估计组织-膳食分馏因子（TDFs）对实测同位素值进行膳食校正。

研究结果部分，在物种特征对稳定同位素值的影响方面，研究发现羽毛δ¹⁵N值在不同生态特征的物种间存在显著差异，尤其与食性和体质量相关。相对于肉食性鸟类（carnivores）作为参照组，杂食性（omnivores）和植食性鸟类（herbivores）的δ¹⁵N值显著更低，其中肉食性鸟类具有最高的δ¹³C和δ¹⁵N值，其次是杂食性，植食性最低。δ¹⁵N值与体质量呈正相关，这一关系在肉食性和杂食性鸟类中均保持显著。对于δ¹³C，食性、迁徙策略和体质量均为显著预测因子：杂食性鸟类的δ¹³C值显著低于肉食性鸟类，而植食性与肉食性之间的差异不显著；留鸟（sedentary）的δ¹³C值高于迁徙物种，该差异在肉食性类群中尤为明显；δ¹³C与体质量在整体模型及杂食性鸟类中均呈正相关。

在环境与人为因素的影响方面，整体模型显示δ¹⁵N值与降水量存在显著正相关，同时与牧草地覆盖度呈显著正相关，而耕地覆盖度则呈现微弱的负相关趋势。这些模式在限制为留鸟数据集的敏感性分析中基本保持一致，但降水量的关联支持度减弱，牧草地的正效应仍保留，耕地保持微弱且非显著的负相关，δ¹³C则始终与环境及土地利用预测因子无关。

在生境类型间生态位宽度的差异方面，研究人员将数据集限制于留鸟后，对农田、灌丛和森林三种栖息地进行了同位素生态位比较。结果显示，植食性鸟类的栖息地相关差异最为显著：灌丛的核心同位素生态位最宽（SEAc = 44.753‰²），农田居中（20.073‰²），森林最窄（5.948‰²），所有成对比较的后验概率均超过0.99。肉食性鸟类呈现较弱模式，农田和灌丛的核心生态位面积均大于森林，但农田与灌丛之间的分离有限。杂食性鸟类的趋势相似但较不明显，灌丛的生态位面积略大于农田和森林，但栖息地对比弱于植食性鸟类，且后验支持度有限。等样本量重采样分析 broadly reproduced 了上述栖息地排序，证实了结果的稳健性。

讨论与结论部分，研究人员综合分析了上述发现的生态学意义。δ¹⁵N值从植食性到肉食性物种的逐步升高反映了沿摄食层次的营养富集模式，这与此前的研究一致。羽毛δ¹³C值在迁徙策略间的变异则可能反映了羽毛形成期间整合的碳源环境差异：留鸟和部分迁徙鸟类因在相对开阔的农田或农牧交错环境中觅食而具有较富集的δ¹³C值，而长途迁徙鸟类因关联于湿润森林或河岸生境而表现出更亏损的δ¹³C值，这些差异主要源于C₃作物（水稻、小麦）与C₄作物（玉米）共存形成的异质性农业基质。

体质量与羽毛稳定同位素组成的正相关关系可能由多种机制共同驱动：大型鸟类代谢周转较慢、羽毛生长期较长，整合了更长时间尺度的膳食输入；同时大型物种如大嘴乌鸦（Corvus macrorhynchos）和喜鹊（Pica pica）倾向于在受C₄作物和畜牧活动影响的开放栖息地觅食，而小型森林物种主要依赖C₃植被和昆虫猎物。

牧草地覆盖度与δ¹⁵N的正相关与耕地覆盖度的微弱负趋势形成对比，反映了不同土地利用类型间氮来源和转化途径的差异。放牧景观中有机氮输入（粪肥、堆肥）通常富集¹⁵N，而农田常用合成氮肥在同位素上较为亏损（接近0‰）。两种互补机制解释了这一模式：一是有机氮输入效应，牲畜密度增加促进粪肥沉积，通过氨挥发和部分反硝化作用提高土壤δ¹⁵N；二是放牧诱导的循环效应，持续践踏和植被移除改变土壤通气和水分状况，刺激微生物活动，加速氮矿化和反硝化。降水量与δ¹⁵N的正相关可能反映了区域水分条件对氮循环强度的影响，高降雨增强土壤微生物活性，促进矿化和反硝化，增加气态氮损失，优先移除¹⁴N而使剩余氮池富集¹⁵N。

同位素生态位分析进一步揭示了生境异质性和土地利用配置如何塑造鸟类资源利用和营养结构。植食性鸟类的栖息地分异最为清晰，灌丛保持最宽的核心同位素生态位，农田居中，森林最窄。灌丛复杂的植被结构和植物多样性提供了多样化的植物食源和庇护条件；农田马赛克（作物兜底、灌溉沟渠、草-作物混合斑块）支持了多营养级多样性，为鸟类提供了从谷物、种子到无脊椎动物的广泛食物资源；森林则可能支持更为均质的植物资源基础。杂食性鸟类的栖息地对比相对较弱，可能因其较高的膳食灵活性缓冲了生境特异性差异。肉食性鸟类中农田和灌丛的核心生态位均大于森林，但农田与灌丛之间的对比有限，可能反映了开放或边缘生境中更高的猎物多样性。

研究结论指出，稳定同位素分析为追踪鸟类营养结构和资源利用模式提供了有力工具，有助于增进对环境梯度间生态位分化的理解。通过将生态、气候和人为变量整合到综合性建模框架中，本研究揭示了牧草地覆盖度、降水量与羽毛δ¹⁵N之间的一致正相关关系，凸显了人类土地利用对食物网内氮循环的深刻影响。此外，食性和迁徙策略显著调控鸟类稳定同位素的变异，表明具有不同生态特征的物种在资源利用和能量传递中表现出不同的响应模式。总体而言，这一跨类群、多尺度的整合框架深化了对四川山地生态系统鸟类群落营养动态的理解，同时为评估人类活动对生态系统功能的长期影响提供了新视角。研究未观察到显著的时间格局，这可能归因于采样年份的不均匀分布；未来研究纳入更多时间段样本或构建更完整的时间序列，将有助于验证人为干扰和气候变化对鸟类稳定同位素动态的长期效应。