扩散和富营养条件对弱混合水体中真菌寄生虫传播影响的建模

《Mathematical Biosciences》：Modelling the impact of diffusion and eutrophic conditions on fungal parasite transmission in poorly-mixed waters

　　

真菌寄生虫在水生生态系统中广泛分布，通过专性感染浮游动物宿主并诱导宿主死亡，在调节浮游生物种群中发挥重要作用。在本研究中，研究人员开发了一个传播动力学模型，用于描述未搅拌恒化器（unstirred chemostat）环境中浮游植物（phytoplankton）、浮游动物（zooplankton）和真菌寄生虫（fungal parasite）之间的相互作用。研究人员研究了模型的动力学性质，包括耗散性（dissipativity）、稳态（steady state）的存在性和稳定性，以及在主特征值（principal eigenvalue）和基本再生数（basic reproduction number, R₀）的适当条件下的一致持久性（uniform persistence）。数值模拟显示，高扩散率（diffusion rate）抑制了物种的长期共存。此外，研究人员的发现表明，富营养条件（eutrophic conditions）可以增强真菌寄生虫的传播并诱导种群波动。

真菌寄生虫在水生生态系统中广泛存在，通过专性感染浮游动物调节浮游生物种群，对水生群落结构具有重要影响。以往研究多集中在充分混合水体，忽略了自然水体中的空间异质性。为探究弱混合环境中扩散和营养梯度对真菌寄生虫传播的综合效应，研究人员在《Mathematical Biosciences》上发表了基于未搅拌恒化器（unstirred chemostat）的反应-扩散模型，系统分析了空间结构、营养水平与疾病传播的相互作用。研究得出以下关键结论：高扩散率或高洗出率抑制物种共存并降低寄生虫传播风险；富营养条件增强寄生虫传播并诱发种群振荡；低扩散或低洗出率促进寄生虫持续传播与复杂动态。该工作为理解弱混合水体中水生疾病的传播机制提供了理论框架，对生态管理具有重要指导意义。

研究人员采用的主要技术方法包括：1）构建一维空间反应-扩散偏微分方程模型（2.1）-（2.3），包含营养物、浮游植物、易感浮游动物、感染浮游动物和真菌孢子五个状态变量，利用Monod函数描述营养物限制的浮游植物生长，Holling II型功能反应刻画浮游动物摄食，并引入扩散系数D和洗出率ω模拟弱混合环境；2）通过线性化分析、特征值理论和比较原理，推导稳态存在性与稳定性条件，定义基本再生数R₀；3）采用有限差分法进行数值模拟，结合拉丁超立方抽样（Latin Hypercube Sampling, LHS）和偏秩相关系数（Partial Rank Correlation Coefficient, PRCC）进行全局敏感性分析，并利用蒙特卡洛模拟评估参数不确定性对R₀的影响。所有参数均基于文献假设，无实际样本队列来源。

研究结果部分：

**3.1 耗散性（Dissipativity）**：通过构造复合变量并利用比较原理，证明了系统（2.1）-（2.3）存在唯一全局正解，且解是耗散的，即所有状态变量最终有界，系统具有全局紧吸引子。

**3.2 稳态（Steady States）**：定义了四种稳态：营养物-only稳态E₁、营养物-浮游植物稳态E₂、无病稳态E₃（浮游植物与易感浮游动物共存）和地方病稳态E₄（所有物种共存）。通过特征值条件推导了E₁全局渐近稳定的条件（λ₁(D, g(N₁))<0）。当λ₁(D, g(N₁))>0时E₂存在唯一且全局吸引。进一步，当λ₁(D, ef(P₂))>0时E₃存在，且当R₀>1时系统一致持久，E₄存在。数值模拟显示，E₃和E₄可表现为稳态或空间非均匀周期解。

**3.3 数值模拟与敏感性分析（Numerical Simulations and Sensitivity Analysis of R₀）**：利用有限差分法模拟不同参数场景，发现所有物种生物量在恒化器入口处最高，向出口递减。全局敏感性分析显示，R₀与传播系数β、最大孢子产量q和营养输入浓度N₀呈显著正相关，与扩散系数D、洗出率ω和浮游植物细胞配额c呈显著负相关。蒙特卡洛模拟进一步验证了这些参数对R₀的驱动作用。

**4. 生态因素的影响（The effects of ecological factors）**：低扩散率（D<0.005）促进稳态共存；中等扩散率（0.0050.18）导致物种洗出。低洗出率ω促进浮游植物水华和周期振荡，高ω抑制水华并稳定系统。高营养输入N₀或低细胞配额c显著增加浮游植物和真菌孢子生物量，增强传播。高传播系数β（>0.42）或高孢子产量q（>34）使系统从无病稳态转向地方病周期振荡，增加动态复杂性。

讨论部分总结了模型的主要发现：扩散和洗出抑制寄生虫传播，富营养促进传播，且高传播强度诱发振荡。研究人员指出，模型假设所有扩散系数相等，若不等则分析困难，且未考虑自然死亡率、光照等因素。未来研究可纳入这些因素进一步探讨。研究结论翻译为：高扩散率抑制物种长期共存；富营养条件增强真菌寄生虫传播并诱导种群波动；低扩散或低洗出率有利于真菌寄生虫在浮游动物中持续传播并引发不稳定动态；高洗出率抑制寄生虫传播并降低水华概率。