《Biological Invasions》：Threatening a keystone species: how climate and host phenology interact to trigger eruptive outbreaks of the invasive psyllid Acizzia uncatoides in a biodiversity hotspot
编辑推荐：
摘要：入侵性木虱Acizzia uncatoides通过定殖于严格特有的关键树种Acacia heterophylla，对留尼汪岛的生物多样性构成严重威胁。为揭示此次入侵的驱动因素，研究人员利用季节性自回归积分滑动平均模型（SARIMA/SARIMAX）分析了
摘要：入侵性木虱Acizzia uncatoides通过定殖于严格特有的关键树种Acacia heterophylla，对留尼汪岛的生物多样性构成严重威胁。为揭示此次入侵的驱动因素，研究人员利用季节性自回归积分滑动平均模型（SARIMA/SARIMAX）分析了跨越海拔梯度的五年（2017–2022）木虱丰度时间序列。木虱种群表现出极端的季节性和高振幅暴发（真空采样峰值超过每分钟35,000个个体），这与在其引入范围内缺乏有效的自上而下调控一致，这些位点未记录到天敌（Baujeu et al. 2025）也证实了这一点。在空间上，侵染严重程度支持了资源集中假说：高海拔的单作状林分维持了大规模暴发，而结构复杂的森林则提供了显著的关联抗性。建模识别出宿主物候（具体表现为呈现绿色嫩梢（幼嫩生长组织）的树木比例）是主要的自下而上驱动因子，表现出严格的营养同步性。相反，非生物因子起次要作用；温度施加了延迟的负反馈（2–3个月滞后），可能是通过加速叶片硬化而非直接的热胁迫。这些发现证实了该岛屿植物宿主对Acizzia uncatoides的高度易感性，并为实施经典生物防治以缓解该云雾森林生态系统的衰退提供了关键基线。虽然SARIMA模型为暴发预测提供了稳健的预测工具，但SARIMAX模型阐明了直接的生态机制。这些发现强调了包括生物防治在内的综合管理的迫切需求，以保护这一特有的关键树种。
论文解读文章
**研究背景与问题**
入侵性森林昆虫对树木健康和森林稳定性构成不成比例的威胁，其中植食性昆虫的暴发可导致大范围树木死亡。然而，目前对复杂高海拔热带生态系统中入侵昆虫种群动态的认知仍存在关键空白，尤其是在气候季节性、宿主物候与天敌释放相互作用方面。留尼汪岛作为马斯卡林群岛生物多样性热点核心区，其特有树种Acacia heterophylla（本地名“Tamarin des Hauts”）是热带山地云雾林的关键种，调控水文循环、防止水土流失并支持丰富的附生植物和特有鸟类。然而，入侵性木虱Acizzia uncatoides（原产澳大利亚）自2010年在留尼汪岛被发现后，专性侵染该特有树种，导致严重落叶和树木衰退。此前研究（Baujeu et al. 2025）已证实高海拔生境中几乎完全缺乏瓢虫捕食者和寄生蜂等天敌，表明该入侵物种处于“天敌逃逸”状态。因此，亟需解析其种群爆发机制，为实施经典生物防治提供基线。
**研究内容与结论**
研究人员基于2017–2022年沿海拔梯度（1350–2350 m）的月度监测数据，利用季节性自回归积分滑动平均模型（SARIMA/SARIMAX）分析了A. uncatoides的种群动态，并测试了三个假设：（i）种群呈爆发性“繁荣-崩溃”动态；（ii）宿主物候（绿色嫩梢比例）是暴发的主要驱动力；（iii）气候条件（温度、降雨）也起促进作用。结果表明：木虱种群表现出极端季节性和高振幅暴发（峰值超35,000个体/分钟），符合天敌逃逸特征；空间上，高海拔单作状林分（火山站点）的侵染程度比结构复杂的混合林（贝布尔森林）高出近100倍，支持资源集中假说；建模识别出绿色嫩梢比例是主要自下而上驱动因子，呈现严格营养同步性，而温度以2–3个月滞后产生负反馈（通过加速叶片硬化），累积降雨的直接作用可忽略。SARIMA模型表现出稳健的预测能力，而SARIMAX模型阐明了近端生态机制。该研究证实了该岛屿宿主的高易感性，并为计划引入的经典生物防治剂（瓢虫Harmonia conformis和跳小蜂Psyllaephagus sp. nov.）提供了关键基线。论文发表于《Biological Invasions》。
**主要关键技术方法**
1. **标准化真空采样**：使用改良园艺吸尘器（G-vac，Stihl SH56D）对每棵树进行1分钟吸虫，样本经冷冻后手动计数或重量法估算（基于1000个冻干成虫的平均质量0.111 mg）。
2. **宿主物候监测**：每月记录每站点10–15棵标记A. heterophylla树木的绿色嫩梢有无，计算比例作为外生变量。
3. **气候数据获取**：从METEOR平台提取日均温度和日累积降雨量，通过空间插值（距气象站>1 km的站点）获取，并计算月均温度和月累积降雨量。
4. **时间序列建模**：遵循Box-Jenkins框架，使用SARIMA（p,d,q）（P,D,Q）
12描述内生动态，并整合外生变量（绿色嫩梢比例、温度滞后、降雨滞后）构建SARIMAX模型；通过KPSS和OCSB检验确定差分阶数，Hyndman-Khandakar算法优化AICc，Ljung-Box检验验证残差白噪声；模型性能用AIC和RMSE评估，并计算精度增益（PG）。
**研究结果**
**季节性空间动态**：A. uncatoides丰度在各站点均呈显著年季节周期，峰值集中于10–12月（旱季末），随后1–2月快速下降。空间上，侵染程度与海拔正相关，但火山站点丰度远高于贝布尔站点：火山站点峰值平均丰度从VC_1400的3358个体（2018年11月）到VC_2350的35,259个体（2020年12月）；贝布尔站点峰值仅378（BB_1350）和668（BB_1550）个体（2018年12月）。年际变异表现为2017–2019年种群增加，2019–2020年下降（除最高位点外）。
**内生动态与SARIMA建模**：KPSS检验表明无需常规差分（d=0），OCSB检验确认需季节性差分（D=1）。最优SARIMA模型均包含季节性成分。非季节性部分：贝布尔1350 m、火山1400 m和1950 m站点选择一阶自回归项（AR1），表明种群惯性；贝布尔1550 m站点由移动平均项（MA1）描述，对随机冲击敏感；火山2350 m站点非季节性成分为零（p=0, q=0），动态几乎完全由年周期驱动。残差诊断显示前四个站点残差为白噪声，但火山2350 m站点仍有显著结构（p<0.05）。
**预测性能**：SARIMA模型对2021–2022年验证期的预测显示，预测值大致落在95%置信区间内，表明训练阶段识别的季节模式稳定。
**SARIMAX驱动因子分析**：变量选择表明环境驱动因子存在1–3个月滞后效应，无即时效应（Lag 0）。自回归项（AR1）在纳入环境预测因子后仍显著，表明种群内在惯性。绿色嫩梢比例：在火山站点呈1个月滞后显著正效应（VC_1950系数=2.478，VC_2350系数=1.404），在贝布尔站点呈2个月滞后效应。平均温度：在所有站点呈2–3个月滞后负效应（系数-0.151至-0.387）。累积降雨：在最终多元模型中系数接近零（<0.001），解释力可忽略。模型验证：除火山2350 m站点外，其余站点残差为白噪声。比较模型性能：环境“专家”模型在VC_1950和BB_1350站点精度增益分别达33.6%和26.2%，而VC_2350站点无增益（-0.3%），AIC支持环境模型。
**讨论与结论**
讨论部分指出：A. uncatoides的爆发性动态符合天敌逃逸假说（ERH），并因岛屿宿主易感性假说（Island Susceptibility Hypothesis）而加剧——岛屿植物长期与强烈植食作用隔离，缺乏有效防御机制，导致A. heterophylla对木虱的“生态天真的”防御不足。该木虱在留尼汪岛几乎专性选择A. heterophylla，而非其他已知替代宿主，支持了入侵害虫主动放弃共进化宿主以利用化学防御薄弱的岛屿物种的适应性转变。空间上，火山站点开放的单作状林分比贝布尔的结构复杂森林导致侵染程度高近100倍，支持资源集中假说，且结构复杂性提供关联抗性。SARIMAX模型确认绿色嫩梢比例是主要自下而上驱动因子，但种群崩溃发生在非生物条件仍有利时，表明存在“营养墙”机制：宿主进入生殖期后，营养分配转向生殖库，叶片硬化形成物理屏障，导致若虫高死亡率，而成虫可转向维持策略但抑制产卵。温度作为动力学调节因子而非绝对限制因子，通过加速叶片硬化缩短营养窗口；降雨通过机械冲刷作用影响若虫，但主要间接调控宿主物候。SARIMA模型提供了稳健的预测工具，而SARIMAX模型阐明了近端机制。研究人员建议采用双重框架：外生模型揭示机制，内生模型用于管理监测。
**研究结论部分翻译**：这些发现证实了该岛屿植物宿主对Acizzia uncatoides的高度易感性，并为实施经典生物防治以缓解该云雾森林生态系统的衰退提供了关键基线。虽然SARIMA模型为暴发预测提供了稳健的预测工具，但SARIMAX模型阐明了直接的生态机制。这些发现强调了包括生物防治在内的综合管理的迫切需求，以保护这一特有的关键树种。
