优化作物空间聚类以最小化病原体入侵：理论框架与木薯花叶病毒应用

本研究针对如何通过优化作物空间布局降低病原体初始入侵速率这一关键问题，开展了作物配置优化以最小化农业病原体入侵的研究。通过构建分析模型，确定了针对高斯、指数及幂律等不同传播核的最小化入侵速率的最佳聚类规模与间隔距离阈值，并应用于木薯褐条病毒(CBSV)的实际景观规划，为农业景观设计提供了理论依据与实用工具。

在全球化与气候变化的双重驱动下，农作物病害的传播速度和范围正以前所未有的规模扩大。尤其在小农经济占主导的非洲地区，木薯褐条病毒(CBSV)等毁灭性病害的暴发常导致区域性粮食危机。传统病害防控多依赖化学药剂或抗病品种，但农药残留、抗性丧失等问题日益突出。与此同时，作物空间配置这一“天然免疫屏障”的潜力长期被低估——尽管早有研究表明宿主聚集程度会影响病原传播效率，但“究竟何种空间布局能最有效阻滞疫情初期扩散”这一关键问题，始终缺乏定量化的理论支撑。

剑桥大学植物科学系的Yevhen F. Suprunenko与Christopher A. Gilligan团队在《Scientific Reports》发表的研究，首次建立了基于病原体传播核特征与宿主分布几何参数的解析模型，精准锁定了最小化初始侵染速率(r)的作物聚类方案。该研究突破了过去依赖复杂计算机模拟的局限，通过数学推导直接计算出关键阈值——当作物斑块边长小于L_H^*、间距大于Δ_H^*时，疫情暴发风险可降至理论最低值。更引人注目的是，团队将这一理论框架成功应用于非洲实地景观，为木薯种植区的CBSV防控提供了可直接落地的空间规划方案。

研究方法的核心在于构建了感染速率的双近似模型。针对高度聚集景观(聚类数少)，通过单聚类接触密度n̅₁推导出r≈β×n̅₁(L_H)；针对分散景观则采用全局密度近似r≈β×n。通过取二者最大值构建r_approximate=β×max[n̅₁(L_H), n]的合成模型，其预测结果与个体模型(IBM)仿真数据高度吻合(最大偏差≤22%)。基于此，通过求解n̅₁(L_H^*)=n的方程，解析出关键阈值L_H^*，并推导出对应间隔距离Δ_H^*=L_H^*×(√(A/A_H)-1)。应用阶段采用Szyniszewska(2020)提供的木薯分布栅格数据(1km分辨率)与Godding等(2023)估计的CBSV幂律传播核(α=3.75, p=0.12)，在喀麦隆与中非共和国边境区域(24×24km²)进行了案例计算。

人工景观与常用传播核的阈值特征

研究系统比较了高斯核、指数核与幂律核等常见传播核的阈值响应。对于短程传播主导的高斯核(尺度参数σ)，当宿主占比A_HA^-1≤0.1时，显式解析解L_H^*≈√(A_H/(b(0)A))与数值解高度吻合；而长尾的幂律核(b₀=3,4)因远程传播效应显著，显式解精度下降。值得注意的是，指数核(b₀=1)的L_H^*值始终大于高斯核，说明相同景观下高斯核传播的初始速率更高。

木薯景观与CBSV的应用实例

在A_HA^-1=0.011的案例区域，计算得CBSV的最佳聚类边长L_H^*=0.27km，间隔Δ_H^*=2.29km。为实现操作可行性，将景观聚合为4×4km²单元后生成空间化建议：在(4-8km,12-16km)单元内，推荐建设1×1km木薯集群并设置900m缓冲带；而(16-20km,12-16km)单元因宿主密度更低，集群尺寸需压缩至500×500m且缓冲带扩至1km。这种“聚类尺寸-缓冲带”的负相关关系为差异化规划提供了依据。

研究通过理论推导与实证验证相结合，揭示了景观空间结构对疫情发展的调控机制。双近似模型的存在使得阈值L_H^*实际为理论最小集群尺寸的上界，这种保守估计恰好为实际应用提供了安全边际。尽管分析聚焦于疫情初期动态(r与R₀)，但作者引用前期工作证实，该阈值对中期(6个月)疫情规模排序仍具预测力。值得注意的是，该方法突破了传统方形网格假设，在补充材料中推导的狭长地块适配公式，展现出对现实农田异质结构的扩展能力。

这项研究的深远意义在于将抽象的景观流行病学理论转化为可执行的规划工具。通过建立传播核尺度参数(a)与集群尺寸(L_H^*)的无量纲关系，使决策者能在病原体传播特征尚未完全明确时，仍能通过参数扫描提前优化布局。尤其对于资源匮乏地区，这种“以空间换时间”的防控策略，为降低农药依赖、实现可持续病害管理提供了新范式。未来研究可进一步整合抗病品种布局、药剂施用等控制措施的时空尺度，构建多目标优化的农业景观设计体系。