生态学报  2018, Vol. 38 Issue (1): 291-297

文章信息

秦丽娟, 王万雄, 张锋, 赵有益.
QIN Lijuan, WANG Wanxiong, ZHANG Feng, ZHAO Youyi.
Allee效应对具有生境恢复的集合种群的影响
Impact of Allee effects on metapopulations with habitat restoration
生态学报. 2018, 38(1): 291-297
Acta Ecologica Sinica. 2018, 38(1): 291-297
http://dx.doi.org/10.5846/stxb201610082023

文章历史

收稿日期: 2016-10-08
修订日期: 2017-09-11
Allee效应对具有生境恢复的集合种群的影响
秦丽娟 , 王万雄 , 张锋 , 赵有益     
甘肃农业大学理学院数量生物学研究中心, 兰州 730070
摘要: Allee效应与种群的灭绝密切相关,其研究对生态保护和管理至关重要。Allee效应对物种续存是潜在的干扰因素,濒危物种更容易受其影响,可能会增加生存于生境破碎化斑块的濒危物种的死亡风险,因此研究Allee效应对种群的动态和续存的影响是必要的。从包含由生物有机体对环境的修复产生的Allee效应的集合种群模型出发,引入由其他机制形成的Allee效应,建立了常微分动力系统模型和基于网格模型的元胞自动机模型。通过理论分析和计算机模拟表明:(1)强Allee效应不利于具有生境恢复的集合种群的续存;(2)生境恢复有利于种群续存;(3)局部扩散影响了集合种群的空间结构、动态行为和稳定性,生境斑块之间的局部作用将会减缓或消除集合种群的Allee效应,有利于集合种群的续存。
关键词: 稳定性分析     局部扩散     元胞自动机模型     生境恢复    
Impact of Allee effects on metapopulations with habitat restoration
QIN Lijuan , WANG Wanxiong , ZHANG Feng , ZHAO Youyi     
Center for Quantitative Biology, College of Science, Gansu Agricultural University, Lanzhou 730070, China
Abstract: Understanding Allee effects is important for ecological conservation and management because it is strongly related to population extinction. The Allee effect is a destabilizing factor in species persistence, and endangered species are more susceptible to the effect as it will increase the risk of extinction associated with living in fragmented habitats. Many mechanisms may account for the Allee effect, such as organism-environment feedback, mate-finding limitations, cooperative defense, social dysfunction, inbreeding depression, and predator avoidance or defense. Because of various ecological mechanisms leading to the Allee effect, a large amount of experimental evidence in recent years has indicated that species have been subjected to multiple Allee effects. Therefore, it is necessary to study the influence of multiple Allee effects on population dynamics and persistence; however, few theoretical studies have been performed on multiple Allee effects. In this study, we began by assuming that metapopulations experience Allee effects caused by organism-environment feedback and also induced the Allee effect through other ecological mechanisms. We constructed a differential dynamical system and we used phase-plane analysis and numerical solutions to illustrate equilibriums and their system stability. Furthermore, we discussed and simulated the influences of the Allee effect on the spatial distribution of a metapopulation by using a cellular automata model based on a lattice model. The results obtained by theoretical analysis and computer simulations suggested that a strong Allee effect could markedly increase the extinction risk for a metapopulation with habitat restoration; (2) Habitat restoration is favorable for population persistence; and (3) local dispersal influences the spatial structure, dynamics, and stability of a population. Our simulations, which explicitly consider spatial structure, also demonstrated that local interaction among habitat patches can greatly mitigate the Allee effect and is favorable for the persistence of a metapopulation. This result suggests that limited local dispersal may be a stabilizing factor.
Key words: stability analysis     local dispersal     cellular automata model     habitat restoration    

Allee效应是低密度时种群的单位个体增长率与种群密度之间的正关系,由于直接关系到种群的灭绝,因此受到生态学和生物保护学的重视[1-3]。引起Allee效应的机制有很多,例如寻找配偶限制[4-5]、合作防御[6-7]、社会功能障碍[8]、近交衰退、捕食防御[9-10]以及生物与环境的反馈[9, 11]等。和局域种群类似,Allee效应也会出现在集合种群,其形成原因主要是扩散代价和侵占困难[8, 12-13]。由于各种各样的生态机制都会导致Allee效应产生,同一个种群可能会遭受两个或多个不同机制的Allee效应,因此两个或多个Allee效应在种群的动态和续存上的影响是很有必要的。Berec等[14]提出两个或多个Allee效应能影响一个种群的动态,并提供了很多的实验证据,目前对同一种群遭受两个或多个Allee效应的理论研究很少。

生物有机体与环境的正反馈(即生物体本身对环境的修复)是导致Allee效应产生的一个重要的生态学机制[9, 15]。生物的生长、存活和新陈代谢受物理环境(如温度、风和气流)和化学环境(氧气、毒素和荷尔蒙)的限制[16]。如果生物有机体在某种方式下能修复和改善自己的环境,较大集群的个体会生存的更好[17]。在生态系统中,生物有机体与环境之间的正反馈是重要的生态过程,种群的适合度通过环境作为中介间接地依赖于种群密度,从而导致产生Allee效应[18-19]。但是,当种群遭受其他机制引起的Allee效应时,生物有机体与环境的反馈对种群动态的影响仍不是很清楚。

对集合种水平上Allee效应的研究,都是基于经典的Levins模型[20-21]。Hui等研究表明,在集合种群水平上,当斑块侵占比例低于临界值时,集合种群将趋于灭绝[8]。McVinish和Pollett研究了集合种群的生境质量与Allee效应之间关系,当集合种群遭受Allee效应时,更需要保护生境免遭退化和破坏[22]。Zhang等在均匀场假设下研究了适宜斑块比例随时间变化的情况下,生境恢复和生境破坏对集合种群续存的影响[11],研究结果表明生物有机体对环境的修复是系统产生Allee效应(双稳态)的关键因素。本文提出了具有生境恢复的集合种群遭受其他机制引起的Allee效应时的均匀场模型,研究了这种Allee效应对集合种群的影响,并运用基于元胞自动机模型的计算机模拟,讨论了局部扩散对具有Allee效应的集合种群的影响。

1 模型建立 1.1 基本模型

经典的集合种群模型是建立在Levins斑块占据模型基础上的,Zhang等[11]人考虑到适宜斑块会因人为活动和自然因素的破坏而变为不适宜斑块,而不适宜斑块会因人类活动等因素的外部恢复和生物体自身内部修复变为适宜斑块,从而建立了下列模型:

(1)

式中,p表示占据斑块比例,h表示适宜斑块比例(1-h表示不适宜斑块比例),显然phc为侵占率,e为灭绝率,d为由于人类活动或自然灾害引起的斑块破坏率,λ为由生物体自身对不适宜斑块的内部修复率,μ表示外在因素(人为活动等)对不适宜斑块的恢复率,其他参数的含义和上面的相同,参数c, e, d, λμ均为非负常数。当λ=0时,方程(1)为经典的Levins模型[20-21]。对系统(1)通过理论分析,结果表明:生物体自身对环境的修复作用是系统产生Allee效应的关键因素,而且对集合种群的续存有很大影响,强调了生物有机体与环境的反馈在生物保护中的重要作用,。

1.2 具有Allee效应的均匀场模型

Hui等[8]明确指出集合种群与局域种群都存在Allee效应,当集合种群过小时,迁移个体数量也将很小,这些个体在迁移过程中由于死亡风险使得到达空斑块的数量降低,在斑块内就会发生传统的Allee效应。对系统(1),考虑到集合种群会遭受其他机制引起的Allee效应,得到下列模型

(2)

式中,为Allee效应,常数a≥0为Allee效应强度,a越大,Allee效应越明显。若a=0,则系统(2)与(1)相同。

易验证系统(2)有唯一边界平衡点E0(0, μ/(μ+d)),其雅可比矩阵的特征值λ1=-(e+d),λ2=-(μ+d)。因此,平衡点E0(0, μ/(μ+d))始终是局部渐近稳定的。进一步通过计算可知,系统(2)最多存在两个内部平衡点,用E1(p1*, h1*)和E2(p2*, h2*)(不妨设p1*p2*)表示,平衡点的精确解表达式很繁琐,可以通过Matlab求其数值解。下面讨论内部平衡点的局部稳定性。

将系统(2)中方程的右边分别记为f(p, h)和g(p, h),两条零等倾线方程分别为,则这两条等倾线上任一平衡点(p, h)处切线斜率分别为h1(p)=-fp/fhh2(p)=-gp/gh.在任一内部平衡点(p*, h*)的雅可比矩阵为

(3)

这里δ=(e+d)/c,易证上式中点(p*, h*)处fh>0, gp>0, gh<0。则有

由两条等倾线的性质,易知在E1(p1*, h1*)处h1<0, h2>0,在E2(p2*, h2*)处h1>h2>0且(p2*)2>(由h1(p2*)>0推出)。由于在E1(p1*, h1*)处det(J)<0,所以平衡点E1是不稳定的;而在E2(p2*, h2*)处det(J)>0, trJ<0平衡点E2是局部渐近稳定的。

1.3 基于元胞自动机模型的空间模拟

虽然动力系统是研究集合种群动态的经典方法,但它不能揭示集合种群的空间动态和分布模式,元胞自动机模型是研究集合种群空间动态和模式的一种重要方法。本文的模拟均在n×n的二维网格上进行,并采用冯.诺依曼邻体(4邻体)和摩尔邻体(8邻体),运用周期边界与同步更替,意味着二维网格像圆环面。每个斑块有3个状态:被物种占有斑块(用2表示)、适宜的空斑块(用1表示)以及不适宜的斑块(用0表示)。用表示第i斑块的邻体中被物种占有的斑块的比率,系统(2)的元胞自动机更替法则可以描述为:不适宜斑块以恢复率变为适宜的空斑块(0→1);适宜的空斑块以破坏率d变为不适宜的斑块(1→0),以侵占率变为被物种占有斑块(1→0);被占有的斑块以破坏率d变为不适宜的斑块(2→0),以自然灭绝率d变为适宜的空斑块(2→1)。具体同步更替法则见表 1

表 1 集合种群的元胞自动机更替法则 Table1 The updating rule of metapopulation Cellular Automata
状态State 比率Rate
0→1
1→0 d
1→2
2→0 d
2→1 e
  0:不适宜斑块,Unsuitable patch; 1:适宜空斑块,Suitable but empty patch; 2:占有斑块,Occupied patch; : 邻体斑块被物种占有的比率,The fraction of occupied patches in its neighborhood; λ:内部修复率,Internal habitat restoration rate; μ:外部修复率,External habitat restoration rate; c:侵占率,Colonization rate; e:灭绝率,Extinction rate; d:破坏率,Destruction rate; a: Allee效应,Allee effect
2 数值模拟与结果分析

首先,我们研究了Allee效应对具有生境修复的集合种群动态的影响。当集合种群遭受Allee效应影响(a>0)时,边界平衡点E0(0, μ/(μ+d))始终局部稳定,系统(2)存在两个内部平衡点(小的不稳定,大的稳定)或没有内部平衡点(图 1 A; 图 2A, C),增强Allee效应会导致种群灭绝, 不利于种群续存。从图 3可以看出,在相同的破坏比例下,生物体自身对环境的内部恢复率(λ)、人类对环境的外部恢复率(μ)以及Allee效应(a)对集合种群的续存都有很大的影响。生境恢复有利于集合种群的续存,而Allee效应不利种群的续存,而且集合种群的局部密度对其续存也有着决定性的影响。当种群局部密度低于某一阈值时,集合种群将无法续存,这表明Allee效应对集合种群的续存至关重要。同时,随着破坏率和Allee效应的增大,即使恢复率较大,也会导致种群灭绝(图 3)。

图 1 扩散距离影响下种群动态的相图 Fig. 1 Effect of dispersal distance on population dynamics shown by phase plane (A)为全局扩散(均匀场假设),(B)和(C)分别为邻体大小为8和4的局部扩散;灰线表示等倾线,黑线表示相平面轨线; 参数取值:λ=0.8, μ=0.1, c=0.8, e=0.1, d=0.1, a=0.1

图 2 Allee效应a对集合种群动态的影响 Fig. 2 Influence of the Allee effect aon the dynamics of metapopulation 参数取值:c=0.8, e=0.1, d=0.1, μ=0.1, (A)、(D)、(G)中λ=0, (B)、(E)、(H)中λ=0.4, (C)、(F)、(I)中λ=0.8

图 3 内部恢复率、外部恢复率与Allee效应对集合种群的续存的响应 Fig. 3 The response of internal restoration, external restoration and Allee effect under the persistence of metapopulation 参数取值:c=0.8, e=0.1, (A)μ=0.1, d=0.05; (B) μ=0.1, d=0.075; (C) μ=0.1, d=0.1; (D)λ=0.8, d=0.05; (E)λ=0.8, d=0.075; (F)λ=0.8, d=0.1

其次,为了研究局部扩散下Allee效应对集合种群动态的影响,我们在200×200的网格上采用元胞自动机模拟,并和均匀场模型的结果作比较。图 1中的(B)和(C)分别是邻体为8和4的元胞自动机模拟结果,结果表明局部扩散减弱或消除了集合种群的Allee效应,而且随着邻体大小的减小,达到平衡状态时的种群大小有明显的减小。图 2描述的是依赖于初始值的种群动态,模拟是在200×200的网格上将Allee效应a和种群大小p分为200×200份, 再分别运行3000步得到的,深灰色区域表示导致种群灭绝的初值,浅灰色区域表示导致种群续存的初值,(A)—(C)是均匀场假设下的情形,结果表明种群受到中等强度的Allee效应时,系统会存在两个平衡点(一个稳定,一个不稳定),当种群受到强Allee效应或种群初始密度较小时,种群趋于灭绝。图 2中(D)—(F)和(G)—(I)分别是8邻体和4邻体下元胞自动机模拟结果,结果表明:生物有机体对环境的修复作用有利于种群的续存,当物种遭受较强Allee效应时,较强的修复率可以防止种群灭绝;局部作用能减弱甚至消除集合种群的Allee效应,由于侵占率和生物体对斑块的恢复率是有限的,因此局部扩散下集合种群的密度低于均匀场情形下的密度,邻体大小不同时密度也有所不同(图 2)。不考虑统计随机性,适度的Allee效应也可导致空间异质性,强Allee效应会导致集合种群的灭绝(图 4)。

图 4 在200×200网格上,受Allee效应的集合种群的空间分布 Fig. 4 Spatial distribution of metapopulation subjected to Allee effect on a 200×200 lattice 时间t=3000,其中黑色、灰色和白色分别表示被占有的斑块、适宜的空斑和不适宜的斑块.参数取值:λ=0.8, μ=0.1, c=0.8, e=0.1, d=0.1
3 讨论

人类活动造成的生境破坏可以使广泛分布的物种退缩到狭小的破碎化生境中, 物种将因为Allee效应、自然灾害、流行病等趋向灭绝,生境破碎和丧失对自然生态系统的微小改变都将迫使人类付出更大的代价[23]。环境(生境)恢(修)复将有助于集合种群中物种的生存,可减小物种的灭绝风险[24]。生境恢复主要包括人为活动的外部恢复和生物有机体对自身环境的内部修复(即生物与环境的反馈)[11],生物有机体与环境的反馈在现实生活中是普遍存在的,特别在诸如干旱地区、潮间带等比较恶劣的环境[25-26],而且这种反馈可以引起物种复杂的空间分布模式[27-28]。因此,研究集合种群的续存与动态问题,不仅要考虑生境的破坏与破碎化,而且要考虑人类对种群环境的恢复和生物有机体对环境的修复。

集合种群是由局部种群组成的种群,每一个局部种群都存在灭绝的风险,种群的续存只能在集合种群上[29]。集合种群Allee效应产生的主要原因是建群困难和扩散损失。当集合种群过小时,发生迁移的个体数量也很小,这些个体在迁移过程中有死亡风险,使得到达适宜空斑块的数量降低,在斑块内就会发生传统意义上的Allee效应[30-31]。本文建立了生境恢复下具有Allee效应的集合种群均匀场模型和元胞自动机模型,通过平衡态分析和数值模拟发现,生境恢复有利于集合种群的续存,而Allee效应不利种群的续存,而且集合种群的局部密度对其续存也有着决定性的影响。

空间模式被很多生态学家关注,在物种的竞争与共存[32-33]、物种的多样性与续存[34]、传染病的传播[35-36]、空间格局的形成[37-38]等中起着关键的作用。空间模式对Allee效应的出现有很重要的意义。当生物体与环境的作用是局部有限时,集合种群的Allee效应被减弱。由局部作用产生的空间模式的出现可以解释Allee效应的这种减弱现象(图 4),因此空间模式的形成可以防止种群灭绝。

参考文献
[1] Allee W C. Animal aggregations, a study in general sociology. Chicago: University of Chicago Press, 1931.
[2] Birkhead T R. The effect of habitat and density on breeding success in the common guillemot (Uria aalge). Journal of Animal Ecology, 1977, 46(3): 751–764. DOI:10.2307/3638
[3] Lander R, Engen S, Sather B E. Stochastic Population Dynamics in Ecology and Conservation. Oxford: Oxford University Press, 2003.
[4] Courchamp F, Berec L, Gascoigne J. Allee effects in ecology and conservation. New York: Oxford University Press, 2008.
[5] Gordillo L F. Modeling ephemeral mating encounters in insects: The emergence of mate-finding Allee effects and applications to theoretical models of sterile release. Theoretical Population Biology, 2015, 104: 10-16.
[6] Zhang F, Hui C, Han X Z, Li Z Z. Evolution of cooperation in patchy habitat under patch decay and isolation. Ecological Research, 2005, 20(4): 461–469. DOI:10.1007/s11284-005-0072-7
[7] Zhang F, Hui C. Eco-evolutionary feedback and the invasion of cooperation in prisoner's dilemma games. PLoS One, 2011, 6(11): e27523. DOI:10.1371/journal.pone.0027523
[8] Hui C, Li Z Z. Distribution patterns of metapopulation determined by Allee effects. Population Ecology, 2004, 46(1): 55–63.
[9] Stephens P A, Sutherland W J, Freckleton R P. What is the Allee effect?. Oikos, 1999, 87: 185–190. DOI:10.2307/3547011
[10] Teixeira Alves M, Hilker F M. Hunting cooperation and Allee effects in predators. Journal of theoretical biology, 2017, 419: 13–22. DOI:10.1016/j.jtbi.2017.02.002
[11] Zhang F, Tao Y, Hui C. Organism-induced habitat restoration leads to bi-stability in metapopulations. Mathematical Biosciences, 2012, 240(2): 260–266. DOI:10.1016/j.mbs.2012.08.006
[12] Zhou S R, Wang G. Allee-like effects in metapopulations dynamics. Mathematical Biosciences, 2004, 189(1): 103–113. DOI:10.1016/j.mbs.2003.06.001
[13] Hammond D S. Post-dispersal seed and seedling mortality of tropical dry forest trees after shifting agriculture, Chiapas, Mexico. Journal of Tropical Ecology, 1995, 11(2): 295–313. DOI:10.1017/S0266467400008762
[14] Berec L, Angulo E, Courchamp F. Multiple Allee effects and population management. Trends in Ecology & Evolution, 2007, 22(4): 185–191.
[15] Ferson S, Burgman M A. The dangers of being few:demographic risk analysis for rare species extinction. New York State Museum Bulletin, 1990, 471: 129–132.
[16] Jones C G, Lawton J H, Shachak M. Organisms as ecosystem engineers. Oikos, 1994, 69(3): 373–386. DOI:10.2307/3545850
[17] Post D M, Palkovacs E P. Eco-evolutionary feedbacks in community and ecosystem ecology: interactions between the ecological theatre and the evolutionary play. Philosophical Transactions of the Royal Society. Series B: Biological Sciences, 2009, 364(1523): 1629-1640.
[18] Seto M, Iwasa Y. Regime shift and robustness of organism-created environments:A model for microbial ecosystems. Journal of Theoretical Biology, 2011, 269(1): 297–306. DOI:10.1016/j.jtbi.2010.10.035
[19] Wertheim B, Marchais J, Vet L E M, Dicke M. Allee effect in larval resource exploitation in Drosophila:an interaction among density of adults, larvae, and micro-organisms. Ecological Entomology, 2002, 27(5): 608–617. DOI:10.1046/j.1365-2311.2002.00449.x
[20] Levins R. Some demographic and genetic consequence of environmental eterogeneity for biological control. Bulletion of the Entomological Society of America, 1969, 15(3): 237–240. DOI:10.1093/besa/15.3.237
[21] Levins R. Extinction. Lecture Notes in Mathematics, 1970, 2(1): 75–107.
[22] McVinish R, Pollett P K. Interaction between habitat quality and an Allee-like effect in metapopulation. Ecological Modelling, 2013, 249: 84–89. DOI:10.1016/j.ecolmodel.2012.07.001
[23] 陈玲玲, 林振山, 梁仁君. 生境丧失对具有似Allee效应集合种群的影响及对策——以江苏盐城为例. 生态学报, 2007, 27(11): 4506–4515.
[24] 谢正磊, 林振山, 齐相贞, 管卫华. 基于栖息地恢复对群落不同种群演化影响的模拟. 生态学报, 2005, 25(6): 1397–1403.
[25] Mager D M, Hui C. A first record of biological soil crusts in the Cape Floristic Region. South African Journal of Science, 2012, 108(7/8): 1–4.
[26] Zarnetske P L, Hacker S D, Seabloom E W, Ruggiero P, Killian J R, Maddux T B, Cox D. Biophysical feedback mediates effects of invasive grasses on coastal dune shape. Ecology, 2012, 93(6): 1439–1450. DOI:10.1890/11-1112.1
[27] Hui C, Li Z Z, Yue D X. Metapopulation dynamics and distribution, and environmental heterogeneity induced by niche construction. Ecological Modelling, 2004, 177(1/2): 107–118.
[28] Han X Z, Hui C, Zhang YY. Effects of time-lagged niche construction on metapopulation dynamics and environmental heterogeneity. Applied Mathematics and Computation, 2009, 215(2): 449–458. DOI:10.1016/j.amc.2009.05.015
[29] Hui C, Zhang F, Han X Z, Li Z Z. Cooperation evolution and self-regulation dynamics in metapopulation:stage-equilibrium hypothesis. Ecological Modeling, 2005, 184(2/4): 397–412.
[30] 宋卫信, 张锋, 刘荣堂. 生境破坏的模式对集合种群动态和续存的影响. 生态学报, 2009, 29(9): 4815–4819.
[31] 刘会玉, 林振山, 温腾. 集合种群动态对栖息地毁坏时空异质性的响应. 生态学报, 2007, 27(9): 3711–3717.
[32] Tilman D. Competiton and biodiversity in spatially structured habitats. Ecology, 1994, 75(1): 2–16. DOI:10.2307/1939377
[33] Snyder R E. Spatiotemporal population distributions and their implications for species coexistence in a variable environment. Theoretical Population Biology, 2007, 72(1): 7–20. DOI:10.1016/j.tpb.2007.03.009
[34] Brandt A J, Seabloom E W. Seed and establishment limitation contribute to long-term native forb declines in California grasslands. Ecology, 2012, 93(6): 1451–1462. DOI:10.1890/11-0579.1
[35] Rass L. Pandemic bounds for an epidemic on an infinite lattice. Mathematical Biosciences, 2005, 195(2): 194–209. DOI:10.1016/j.mbs.2005.02.005
[36] Tuckwell H C, Toubiana L. Dynamical modeling of viral spread in spatially distributed populations:stochastic origins of oscillations and density dependence. Biosystems, 2007, 90(2): 546–559. DOI:10.1016/j.biosystems.2006.12.006
[37] Nowak M A, May R M. Evolutionary games and spatial chaos. Nature, 1992, 359(6398): 826–829. DOI:10.1038/359826a0
[38] Zhang F, Li Z Z, Hui C. Spatiotemporal dynamics and distribution patterns of cyclic competition in metapopulation. Ecological Modelling, 2006, 193(3/4): 721–735.