植物-蜜蜂互惠共生模型在沙尘暴效应下的动力学分析

《Scientific African》：Dynamical Analysis of a Plant–Honeybee Mutualism Model with Dust Storm Effects

　　

开花植物与蜜蜂（传粉者）之间的相互作用是生态系统中最常观察到的现象之一。研究通常关注简化的系统，聚焦于开花植物及其传粉者，以及针对其中一种或两种种群的捕食者。这些研究忽视了生态现象对相互关系的影响。基于上述观察，研究人员提出一个数学模型，用以考察沙尘暴对开花植物与蜜蜂种群之间互惠共生关系的影响。研究了解的正性、有界性和持久性；平衡点的可行性；局部和全局稳定性；以及某些分岔及其特征的基本结果。研究人员观察到，沙尘暴频率的增加导致共存平衡点处开花植物和蜜蜂的密度下降。数值结果表明，沙尘暴可能导致开花植物灭绝和蜜蜂种群衰退。开花植物对可持续生态系统至关重要，有助于生态系统内植物和蜜蜂种群的可持续性和稳定性。这些结果可为旨在保护蜜蜂和增强生物多样性的保护计划提供指导。

**论文解读文章**

**研究背景与问题**
生态系统中的开花植物与传粉者（如蜜蜂）之间的互惠共生关系是维持生物多样性和农业生产力的关键。全球约30万种植物依赖动物传粉者（包括蜜蜂）的传粉服务，这一过程通过异花授粉显著提升作物产量、遗传多样性和品质。然而，传粉者种群在全球范围内呈现下降趋势，主要归因于栖息地破坏、疾病、气候变化、农药使用和捕食等多重环境因素。沙尘暴作为其中一种不利因素，通过阻塞植物叶片气孔、阻碍光合作用、干扰蜜蜂嗅觉信号和飞行能力，对植物-蜜蜂互惠系统产生负面影响。尽管已有研究探讨了植物-传粉者模型中的捕食效应、蚂蚁和花蜜盗取等因素，但沙尘暴对植物-蜜蜂动态系统的综合影响尚未被充分建模和分析。因此，研究人员开展本研究，旨在阐明沙尘暴如何影响开花植物与蜜蜂的互惠关系，并评估其对生态系统平衡和可持续性的长期影响。

**研究内容与结论**
研究人员构建了一个由三个常微分方程组成的数学模型（PFH系统），分别描述沙尘暴浓度（P(t)）、开花植物生物量（F(t)）和蜜蜂种群数量（H(t)）的动态变化。模型考虑了沙尘暴的自然输入与耗散、植物对沙尘的清除作用（采用改进的Holling II型功能响应，引入风速因子(1+φ)）、植物的逻辑斯蒂增长、蜜蜂对植物的授粉益处、植物对蜜蜂的营养供给，以及沙尘暴对植物和蜜蜂的伤害效应。通过分析，研究人员证明了系统解的唯一性、正性、有界性和持久性，确定了所有可能的平衡点（沙尘暴点、无蜜蜂点、共存点），并利用Lyapunov方法分析了各平衡点的局部和全局稳定性条件。进一步，研究人员识别了横截分岔、鞍结分岔和Hopf分岔的发生条件，并通过数值模拟验证了理论结果。论文发表在《Scientific African》。

**主要关键技术方法**
研究人员采用常微分方程（ODE）建模，基于Holling II型功能响应描述沙尘清除与伤害过程，利用Lyapunov函数和Routh-Hurwitz判据进行稳定性分析，通过Sotomayor定理和Hopf分岔理论识别分岔类型，并采用经典四阶Runge-Kutta方法（步长h=0.001）进行数值模拟，初始条件设为(P₀,F₀,H₀)=(0.4,0.8,0.12)。模型参数来源于文献估计或合理假设（如风速φ=3，沙尘输入率Q=10 μ g m^-3 day^-1等）。

**研究结果**
- **模型建立（Model formulation）**：提出了包含沙尘暴、开花植物和蜜蜂的PFH系统，明确了各参数（如∝₁为植物清除沙尘率，∝₂为沙尘对植物伤害率，∝₃为沙尘对蜜蜂伤害率，β为蜜蜂对植物的益处率，γ为植物对蜜蜂的益处率等）的生态意义和单位。
- **模型分析（Model analysis）**：
- **唯一性（Uniqueness）**：系统满足Lipschitz条件，解唯一存在。
- **正性（Positivity）**：在正初始条件下，所有解保持正性。
- **有界性（Boundedness）**：当满足γL≤δ₃且s+βH₀>δ₂时，解一致有界。
- **持久性（Persistence）**：在满足sb>∝₂+bδ₂、γJ₁>∝₃(1+φ)P_Max/[1+c(1+φ)P_Max]+δ₃及∝₁(1+φ)F_Max>δ₁的条件下，系统一致持久。
- **平衡点分析（Equilibrium analysis）**：识别出沙尘暴点E₁=(P₁,0,0)、无蜜蜂点E₂=(P₂,F₂,0)和共存点E₃=(P₃,F₃,H₃)，并给出各平衡点存在的参数条件。
- **局部稳定性（Local stability）**：通过Jacobian矩阵特征值分析，得出E₁局部渐近稳定当s<∝₂(1+φ)P₁/[1+b(1+φ)P₁]+δ₂；E₂局部渐近稳定需满足γF₂<δ₃+∝₃(1+φ)P₂/[1+c(1+φ)P₂]等三个条件；E₃局部渐近稳定需满足Routh-Hurwitz判据S_i>0(i=1,3)及S₁S₂>S₃。
- **全局稳定性（Global stability）**：利用Lyapunov函数，E₁全局渐近稳定当(Lδ₁)(∝₁(1+φ)P/[1+b(1+φ)P])²≤s≤δ₂且(β+γ)F_Max≤δ₃；E₂全局渐近稳定需满足附加条件(24)；E₃全局渐近稳定需满足条件(25)和(26)。
- **分岔分析（Bifurcation）**：证明了在参数∝₂=∝₂*处发生横截分岔（TB）于E₁；在∝₁=∝₁*处发生鞍结分岔（SNB）于E₂；在∝₃=∝₃*处发生SNB于E₃；在δ₂=δ₂*处发生Hopf分岔（HB）于E₂；在∝₁=∝₁*处发生HB于E₃，并给出稳定极限环存在的条件（ρ₁₀<0）。
- **数值模拟（Numerical simulations and discussions）**：
- 沙尘输入率Q：低Q值（<1.49）时系统稳定于共存点；Q∈[1.49,9.4]时出现Hopf分岔（振荡）；Q>9.4时恢复稳定；Q=15时发生横截分岔，系统趋于无蜜蜂点。
- 植物清除沙尘率∝₁：低∝₁（<0.014）时系统趋于沙尘暴点；∝₁∈[0.014,0.18)时趋于无蜜蜂点；∝₁∈[0.18,0.27)时共存；∝₁=0.27时发生Hopf分岔，出现振荡。
- 风速φ：低风速（φ<15）时系统稳定共存；φ≥15时出现极限环（Hopf分岔）。
- 沙尘自然耗散率δ₁：δ₁∈(0,0.08)时共存；δ₁=0.08时出现Hopf分岔；δ₁≥1.43时恢复共存。
- 植物内在增长率s：s≤0.045时出现Hopf分岔（振荡）；s>0.045时稳定共存。
- 沙尘对植物伤害率∝₂：∝₂<0.51时共存；∝₂∈[0.51,0.57)时趋于无蜜蜂点；∝₂≥0.57时趋于沙尘暴点（两次横截分岔）。
- 沙尘对蜜蜂伤害率∝₃：∝₃<0.27时共存；∝₃≥0.27时趋于无蜜蜂点（横截分岔）。
- 蜜蜂对植物益处率β：β<0.07时趋于沙尘暴点；β≥0.07时共存（横截分岔）。
- 植物对蜜蜂益处率γ：γ<0.022时趋于无蜜蜂点；γ≥0.022时共存（横截分岔）。
- 植物自然耗散率δ₂：δ₂<0.06时共存；δ₂∈[0.06,0.14)时趋于无蜜蜂点；δ₂≥0.14时趋于沙尘暴点（两次横截分岔）。

**讨论与结论**
研究人员通过理论分析和数值模拟揭示了沙尘暴对植物-蜜蜂互惠系统的关键影响。沙尘暴强度的增加会降低共存平衡时植物和蜜蜂的密度，甚至导致系统失稳（出现振荡或灭绝）。风速和植物清除沙尘率等参数可显著改变系统动态，表明植物在缓解沙尘暴效应中起核心作用。研究结论指出：当系统达到唯一共存平衡时，沙尘暴效应增强会降低植物和蜜蜂密度，可能破坏系统稳定性；开花植物对可持续生态系统至关重要，支持植物和蜜蜂种群的稳定与恢复力。这些发现可为保护传粉者和促进生物多样性的保护计划提供理论依据。未来研究可进一步考虑湿度、随机沙尘事件、农药效应、疾病和温度等外部因素，以完善模型。