多尺度威胁塑造受胁水栖蝾螈(Necturus lewisi)的占域动态并揭示可能的灭绝债务(Extinction Debt)

《Ecology and Evolution》：Multiscale Threats Shape the Occurrence Dynamics of a Threatened Aquatic Salamander and Reveal a Possible Extinction Debt

　　

摘要：淡水生态系统受人为胁迫因子影响，导致约三分之一淡水动物面临灭绝威胁。Neuse河水生蝾螈（Necturus lewisi）是美国北卡罗来纳州东部Neuse河及Tar-Pamlico河流域特有的大型水栖蝾螈，于2021年被列入美国《濒危物种法》联邦受胁物种。生境退化被认为是驱动N. lewisi占域（occurrence）的主导威胁且效应可能存在滞后。美国鱼类及野生动物管理局（USFWS）草案恢复计划将调查该物种种群周转（colonization/extinction）的占域动态列为优先行动。研究人员假设：局地高质量河道内栖息地的存在会降低灭绝概率；集水区（contributing watershed）内扰动土地覆盖比例升高会增大灭绝概率；强干旱年份的灭绝概率亦升高。研究人员基于176个位点连续5年的检测/未检测（detection/non-detection）数据，在多尺度动态占域模型（dynamic occupancy modeling framework）中检验上述假设，并推导年占域率、周转率（turnover）及平衡占域率（equilibrium occupancy，稳定性）以判断空间上对威胁的响应是否滞后——即灭绝债务。最优模型支持位点特异性灭绝概率驱动假说，包括栖息地质量负向效应、流域内建设用地覆盖正向效应及年干旱强度正向效应。推导估计表明：乡村亚种群（即管理单元）年占域率最高，城市亚种群周转率最大，Neuse河流域多数亚种群平衡占域率低于维持稳定所需水平。占域动态及其衍生量暗示城市亚种群存在灭绝债务，随机干旱事件可能加速其兑现。本研究在灭绝债务背景下创新性应用动态占域模型框架，并为合作保护机构制定N. lewisi恢复方案提供指导信息。

论文解读——《Ecology and Evolution》发表：多尺度威胁塑造受胁水栖蝾螈(Necturus lewisi)的占域动态并揭示可能的灭绝债务(Extinction Debt)

研究背景

淡水生态系统支撑着与其面积不成比例的生物多样性，却承受生境退化、污染、水文改变及入侵物种等多重人为压力，致使约三分之一淡水动物濒临灭绝，依赖流水（lotic）系统的类群尤甚。美国东南部是水生生物多样性热点区域，特有种丰富但衰退显著。Neuse河水生蝾螈（Necturus lewisi）为北卡罗来纳州Neuse河与Tar-Pamlico河流域特有大型水栖蝾螈，历史广布于Piedmont及Coastal Plain区中大型溪流，2021年依《濒危物种法》列为受胁种。其分布退缩主因包括生境丧失、土地利用转变、污染、水坝及气候变化，自1980年至2015年占域位点减少约35%，上Neuse河流域受罗利—达勒姆都会区扩张影响尤为严重。已有研究表明局地生境质量（尤底质组成）是决定其占域的重要因子，景观尺度（沉积、土地利用变化）亦协同作用，但尚未有多尺度（局地+景观）叠加气候时间胁迫的综合动态占域研究。长寿命、低扩散率、狭域特有种在生境破坏后种群响应常具滞后性，即"灭绝债务（Extinction Debt，指生境不再适宜但物种暂存于退化生境，最终因无法维持而局地灭绝的延迟现象）"。传统灭绝债务评估需个体生命表矩阵，对难观测种不现实，而动态占域模型（Dynamic Occupancy Model，可估算初始占域ψ、年定殖概率γ、年灭绝概率ε并校正检测不完全性）可用较易获取的detect/non-detect数据推估斑块周转及平衡占域，为评估斑块水平灭绝债务提供替代途径。USFWS草案恢复计划明确需量化威胁大小及物种种群对威胁的响应，故有必要开展此研究。

研究人员利用2018—2023年五个冬季横跨176个位点的detect/non-detect调查数据，拟合多尺度协变量动态占域模型，量化局地生境质量、流域建设用地比例及年最大干旱等级对ε与γ的影响，并依最优模型参数推求各USFWS指定亚种群（管理单元）的年占域、平均周转率τ及平衡占域ψ_eq=γ/(γ+ε)，通过对比当前有限样本占域估计ψ^(fs)与ψ_eq判断非均衡（non-equilibrium，ψ^(fs)?ψ_eq）方向及潜在灭绝债务风险。

主要关键技术方法

以历史分布区内按HUC10分层选取176个≥5阶（Strahler classification）永久溪流位点，2018—2023年连续5个冬季每位点布设10个Gee's Minnow Trap连续4夜（40 trap-nights），记录N. lewisi检出与否构建检测历史矩阵。于各点开展快速生境评估（Rapid Habitat Assessment，八项指标加权求和得0—100生境总分，区分Piedmont与Coastal Plain评分标准）；从2019 NLCD再分类土地覆被，用IDW-Plus工具HAiFLS（hydrologically active inverse flow length to stream）计算集水区各用地类型有效比例（遇坝排除上游贡献）；依据美国干旱监测（US Drought Monitor）赋予年最大干旱强度分级（0=None, 1=D0, 2=D1, 3=D2）；借邻近USGS站与流域面积比插补无测站位点日均流量并算TQ_mean表征洪枯变幅。采用R语言"unmarked"包拟合全局动态占域模型（初占ψ受生境分、底质、覆盖、LULC、TQ_mean、干旱指数影响；γ与ε受生境总分、LULC（开发/草场）、干旱分级影响；检测概率p受log(日均流量)及饵料老化影响），AIC_c选优，GOF检验无过离散。依最优模型参数算每位点年ψ^(fs)、γ、ε，进而对各亚种群平均得τ^(AVG)=γ?×(1?ψ?_t?1)+(1?ε?)×ψ?_t?1?ψ?_t?1（简言之反映年状态翻转倾向）及ψ_eq=γ?/(γ?+ε?)，参数不确定度由参数Bootstrap 500次获95%CI。

研究结果

3.1 Survey Results（调查结果）

五年间于Neuse及Tar河流域81个位点检出N. lewisi（n=691），原始占域率≈0.46。39个阳性位点跨年多季出现至少一年未检出。快速生境评估分值34—95（均值70.45±12.32），86%阳性位点总分≥60。

3.2 Occupancy Modeling Results（占域模型结果）

GOF无过离散，ΔAIC_c<2仅一不支持模型。最优模型：检测概率p受ln(日均流量)显著正向影响（β=0.311±0.069, p<0.00001）；初占ψ受生境总分正向但不显著（p=0.544）；定殖γ受流域开发用地正向边缘显著（β=4.91±2.983, p=0.099）；灭绝概率ε受生境总分负向（β=?1.290±0.458, p=0.005）、开发用地比例正向（β=9.18或表中10.23±3.904, p=0.016）、干旱分级正向（Drought 2: β=3.93±1.835, p=0.032; Drought 3: β=4.58±1.887, p=0.015）显著影响。无本地—景观交互作用入选。平均p=0.476，平均初占ψ=0.466，平均γ=0.185，平均ε=0.285。年有限样本占域五年波动微小（0.430—0.481）。亚种群比较：上Neuse最低ψ^(fs)=0.253，Trent最高ψ^(fs)=0.656；乡村亚种群（Contentnea, Lower Tar, Trent等）占域较高。

3.3 Turnover and Equilibrium Occupancy Results（周转率与平衡占域结果）

亚种群内年周转率变异小，亚种群间差异大。上Neuse τ^(AVG)=0.742（最高），Swift/Middle τ^(AVG)=0.567次之；其余多数≤0.29。平衡占域ψ_eq上Neuse与Flat/Eno低于当前ψ^(fs)（非均衡差?0.03与?0.06），暗示当前略低于维持稳定所需；Swift/Middle当前ψ^(fs)高于ψ_eq（差+0.04）但周转率极高，暗示尚未表现明显丧失的潜在灭绝债务。Fishing Creek、Little River等乡村亚种群ψ^(fs)>ψ_eq且τ低，近稳定。

讨论与结论翻译（总结）

最优模型证实N. lewisi位点灭绝概率受三重多尺度威胁驱动：局地河道生境质量（负效应，且在干旱年低质生境急剧抬高ε）、流域开发用地比例（正效应，>25%时ε显著上升，干旱年升幅更大）及年干旱强度（正效应）。初占ψ不受生境质量显著驱动，暗示成体可短期耐受劣质生境但长期稳定性敏感。生境质量对ε影响在无非干旱年微弱，干旱放大劣质生境灭绝风险，说明优质生境可缓冲干旱致局部灭绝。开发用地对ε正向效应可能与城市流域基底均质化、径流洪变（flashy flow）及底栖巢穴种丧失有关，TQ_mean未入选或因被开发用地变量吸收。草场/农牧用地未入选最优模型，未来转城镇用地方为主要关切。干旱历史上曾致Upper Tar亚种群丧失，气候变暖预期干旱频发将放大生境与城市化的灭绝促进效应。

周转率与平衡占域联合指示：Neuse河Piedmont区亚种群（Upper Neuse、Flat/Eno、Swift/Middle）周转率最高，Upper Neuse与Flat/Eno ψ^(fs)< />_eq表明灭绝债务正在兑现（已现占域下滑），Swift/Middle ψ^(fs)>ψ_eq但高τ暗示滞后兑现的潜在灭绝债务（相对更外围新区刚开发）；其他农村亚种群（Upper Tar、Lower Tar等）现ψ^(fs)接近或略低于ψ_eq，提示局部生境退化或历史干旱遗留的非均衡风险。整体五年观测占域微幅波动系长寿命成体延缓显现，不代表无威胁。

结论：Neuse及Tar-Pamlico河流域景观城市化与局地生境退化已令多处现/历史占域点不适于N. lewisi存续；占域动态推导暗示部分亚种群特别是Neuse河Piedmont城市亚种群存在灭绝债务，随机干旱可加速其偿还并引致分布收缩至最乡村残遗种群。动态占域模型框架结合集合种群理论可凭detect/non-detect数据初步评估长寿命隐蔽种的灭绝债务风险，为受胁物种恢复规划提供量化依据。本成果正纳入USFWS行将发布的正式恢复计划，指导生境修复与监测优先序设定。