好的，作为学术专家和评审专家，我将严格按照您的要求对这篇论文进行分析和解读。 中文标题 基于反应-扩散模型的入侵性林木害虫传播量化：英国橡树 processionary moth 扩张动态解析

《Mathematical Biosciences》：Reaction-diffusion models of invasive tree pest spread: quantifying the expansion of oak processionary moth in the UK

　　

专业术语翻译 英国林地、森林和城市树木景观正受到外来入侵物种的威胁，气候变化、贸易和运输加剧了这种威胁。入侵性林木害虫会削弱寄主并破坏森林生态系统，因此，对其进行定量建模和预测其传播至关重要。针对这一问题，本研究采用时空反应-扩散方程（reaction-diffusion equation）对一种入侵性害虫的传播进行建模，将空间分布表示为种群密度场（population density field）。研究人员对该难解方程进行了数值求解，并根据解确定了害虫到达场地中各位置的首达时间（first arrival times）。所采用的模型允许研究人员直接从模型参数中获得害虫传播的扩张速率（expansion rate），并通过马尔可夫链蒙特卡洛（Markov Chain Monte Carlo，MCMC）方案，在贝叶斯范式（Bayesian paradigm）下对这些参数进行推断。研究人员将该框架应用于英国橡树 processionary moth（Oak Processionary Moth，OPM）的持续蔓延，尽管采取了管理措施，但该虫害仍在不断扩张。研究表明，该方法有效地捕捉了害虫的传播过程，并且具有非恒定扩张速率的模型比恒定速率模型更符合数据。最后，研究人员讨论了表观非恒定扩张速率的可能原因，认为这可能是由于寄主物种和环境异质性的影响，而非内在扩张速率的真正变化。所提出的框架是定量模拟入侵性林木害虫传播的有力工具，可为未来的预测和管理方法提供支撑。

论文解读文章

研究背景与问题

英国乃至全球的林地、森林和城市树木景观正面临来自入侵物种的严峻生态、经济和流行病学威胁。据估计，这些入侵物种每年给英国经济造成超过40亿英镑的损失。英国境内已有超过100种入侵性林木害虫，其传入途径主要与国际贸易有关，而气候变化则被认为有利于这些害虫建立卫星种群并加速扩散。

其中，橡树 processionary moth（Oak Processionary Moth，OPM）是一个典型案例。该害虫于2006年随活橡木进口在英国南部成功定殖。OPM幼虫不仅会取食树叶导致树木衰弱，其体表的蜇毛（urticating setae）还会脱落并飘散在空气中，接触后可引起人和动物的皮炎，吸入后甚至会导致呼吸道刺激。自定殖以来，OPM的分布范围持续扩大。尽管政府自2011年起实施了管理策略，但疫情仍在增长。现有研究指出，到2050至2070年，英格兰北部和苏格兰的气候条件也将变得适宜OPM生存。

因此，建立一个稳健的框架来模拟入侵性林木害虫（尤其是OPM）的传播，以预测其扩张趋势并为政府的监测和管理计划提供依据，显得尤为迫切。然而，现有的研究方法存在局限性：例如，网络方法需要明确定义节点，难以扩展到新区域或模拟连续空间；传统的区室模型虽然能刻画种群动态，但无法直接给出从模型参数中推导出的扩张速率。此前对OPM扩张速率的最新估算采用了最大距离法，发现了两个不同速率的扩张阶段（2006-2014年为1.66公里/年，2015-2019年为6.17公里/年），但这种方法并未建立一个可用于未来预测的动态数学模型。

研究内容与结论

为解决上述问题，本研究开发了一种基于反应-扩散方程的数学框架来模拟OPM的传播，并从观测数据中推断模型参数。研究人员证明了该方法的有效性，发现一个允许扩张速率随时间变化的模型（双相模型）比恒定速率模型（单相模型）更能拟合观测数据。研究还表明，表观的非恒定扩张速率更可能是由寄主植物分布和环境异质性等外部因素造成的，而非害虫自身生物学特性的根本改变。这一框架为定量模拟入侵性林木害虫的传播提供了强大工具，有望为未来的预测和管理策略奠定基础。

关键技术方法

1. 1.
   反应-扩散模型：采用Fisher-Kolmogorov-Petrovsky-Piskunov（FKPP）方程对OPM的时空传播进行建模。该方程结合了描述随机扩散过程的扩散项（由扩散系数 D 控制）和描述种群指数增长直至环境容纳量 K 的逻辑斯蒂增长项（由内禀增长率 r 控制）。
2. 2.
   数据预处理：使用了2006年至2023年间英国林业委员会和南安普顿大学提供的合并调查数据。通过凸包层分解（convex layer decomposition）方法处理原始观测点，剔除了波前异常值（如飞得过远的雄性个体）和波后内部点，最终保留了414个代表OPM首次到达年份的观测站点。
3. 3.
   贝叶斯推断与MCMC：采用贝叶斯范式来推断模型参数（如扩散系数 D 和增长率 r）。具体地，使用自适应马尔可夫链蒙特卡洛（MCMC）方案，从参数的先验分布出发，结合观测数据，计算参数的后验分布，从而得到最可能的参数值及其不确定性。
4. 4.
   模型比较：提出了两种模型——单相模型（假设恒定扩张速率）和双相模型（假设在某一时间点前后扩张速率发生变化）。通过计算似然函数并比较两种模型对数据的拟合优度，来确定哪个模型更优。

研究结果

3.1 单相模型（Monophasic Model）

通过MCMC对单相模型参数进行推断，得到了扩散系数 D 和增长率 r 的后验分布。由此计算出OPM的恒定扩张速率为 v = 2√(rD) = 4.67 km/yr（95%置信区间为 [4.06, 5.43] km/yr）。模型模拟的OPM前沿位置与观测数据总体吻合较好，但在早期年份（约2015年前）倾向于高估扩张速度，而在后期年份（约2015年后）则低估了扩张速度。

3.2 双相模型（Biphasic Model）

鉴于单相模型在不同时期的表现差异，研究人员引入了双相模型，允许在某个切换时间 τ 前后，扩张速率发生变化。MCMC推断结果表明，切换时间 τ 的后验均值约为2015年。在 τ 之前，扩张速率 v₁ 约为2.96 km/yr（95% CI: [1.42, 4.44] km/yr）；在 τ 之后，扩张速率 v₂ 显著增加到约6.82 km/yr（95% CI: [5.41, 8.47] km/yr）。通过比较边际似然（marginal likelihood），发现双相模型对数据的拟合效果明显优于单相模型，表明OPM的扩张确实经历了两个不同的阶段。

3.3 模型验证（Model Validation）

为了检验模型的预测能力，研究人员进行了留一法交叉验证（leave-one-out cross-validation）。对于每年的数据，用其余年份的数据训练模型，然后预测被剔除年份的OPM前沿位置。结果显示，模型预测的前沿位置与实际观测的前沿位置高度一致，验证了该反应-扩散模型框架的有效性和鲁棒性。

3.4 敏感性分析（Sensitivity Analysis）

研究人员还对初始条件的设定和波前异常值的处理进行了敏感性分析。结果表明，即使不剔除那些位于主波前前方的“异常”观测点，推断出的单相模型扩张速率也几乎没有变化。这说明模型的核心结论对数据预处理中的某些主观选择具有较强的稳健性。

总结与讨论

本研究成功地将反应-扩散模型（FKPP方程）应用于量化入侵性林木害虫OPM在英国的空间扩张。研究证实，一个允许扩张速率发生变化的双相模型比恒定速率模型能更好地描述OPM在2006至2023年间的传播过程。推断出的两个扩张阶段（~3 km/yr 前期和 ~7 km/yr 后期）与此前基于最大距离法的研究结果定性一致。

在讨论中，研究人员强调，观察到的扩张速率变化并不一定意味着OPM的内在生物学特性发生了改变（如扩散能力或繁殖率突变）。一个更合理的解释是，这种表观变化是由环境异质性驱动的。例如，OPM初期主要在伦敦西南部的城市公园（如Richmond Park和Bushy Park）中扩张，那里橡树密度高且集中；随着种群向外蔓延，进入了橡树密度较低、分布更分散的乡村地区，这可能减缓了扩张。而当种群进一步适应并跨越了某些地理或生态屏障后，又可能在新的适宜栖息地加速扩张。此外，管理措施的强度变化也可能影响观测到的扩张模式。总之，该研究提出的反应-扩散框架能够有效捕捉复杂的时空传播动态，并为未来在考虑环境异质性的情况下进行更精确的预测提供了坚实的基础。