莫昌健, 韓曉卓

生態(tài)位構(gòu)建作用對捕食-競爭系統(tǒng)的影響

莫昌健1,2, 韓曉卓1,*

1. 廣東工業(yè)大學(xué), 應(yīng)用數(shù)學(xué)學(xué)院, 廣州 510512
2. 廣東松山職業(yè)技術(shù)學(xué)院, 基礎(chǔ)部, 韶關(guān) 512000

生態(tài)位構(gòu)建是指有機體所具有的改變其生存環(huán)境以促進自身生存與發(fā)展的能力。文章通過數(shù)值模擬競爭環(huán)境下具有捕食關(guān)系的各物種間的相互作用, 探討了生態(tài)位構(gòu)建作用對于捕食-競爭系統(tǒng)中物種數(shù)量變化及其空間分布模式的影響。研究結(jié)果表明: 物種的生態(tài)位構(gòu)建行為不會破壞生態(tài)環(huán)境的穩(wěn)定性, 但能有效提高構(gòu)建物種的競爭力, 維持其在不利環(huán)境中與競爭物種共處的平衡狀態(tài); 同時, 生態(tài)位構(gòu)建行為加強了構(gòu)建物種的空間擴散能力, 促使構(gòu)建物種迅速遷移至有利于自身生存的空間生境中, 成為捕食-競爭系統(tǒng)中的穩(wěn)定生存物種; 而由于可供分配的資源有限, 有機體的構(gòu)建行為與擴張行為之間存在的trade-off關(guān)系使其空間分布區(qū)域維持在一定的平衡狀態(tài), 不會出現(xiàn)無限制的空間擴散。

生態(tài)位構(gòu)建; 反應(yīng)擴散方程; 妥協(xié); 平衡狀態(tài)

1 前言

生態(tài)位構(gòu)建是指有機體通過新陳代謝、活動和選擇, 部分的創(chuàng)建和部分的毀滅自身生態(tài)位、改變環(huán)境, 進而改變其環(huán)境中生物與非生物的自然選擇源的一種能力[1–5]。自生態(tài)位構(gòu)建理論提出以來, 一系列關(guān)于生態(tài)位構(gòu)建的研究成果的發(fā)表, 使得該理論逐漸引起生態(tài)學(xué)家與進化生物學(xué)家的興趣[6–8]。2004年惠蒼[9]首先將生態(tài)位構(gòu)建作用架構(gòu)到集合種群水平上; 2006年韓曉卓等人[10–11]討論了具有生態(tài)位構(gòu)建的有機體或種群對破壞生境的重新修復(fù)作用以及相應(yīng)的集合種群動態(tài)與空間分布的變化。Kylafis和Loreau于2008年[12]研究了單一物種的生態(tài)位構(gòu)建對于環(huán)境及進化的影響, 并于2011年[13]以Leibold[4]所提出的模型為基礎(chǔ), 構(gòu)建了具有生態(tài)位構(gòu)建行為的捕食-競爭系統(tǒng)模型。2014年韓曉卓和惠蒼[14]研究了環(huán)境梯度下生態(tài)位構(gòu)建行為對于物種基因分布及適應(yīng)性景觀的影響, 同年Buser[15]研究了生態(tài)位構(gòu)建行為對于互利共生系統(tǒng)下生物進化的影響。經(jīng)過學(xué)者的不斷補充完善, 生態(tài)位構(gòu)建理論已在種的適合度測定[16]、優(yōu)勢度分析[17]、種間競爭與共存[18–19]、群落結(jié)構(gòu)與多樣性維持[20–23]、物種生態(tài)位構(gòu)建與進化對策[24–25]等研究中廣泛應(yīng)用, 關(guān)于生態(tài)位構(gòu)建的定量研究[26]亦逐漸受到重視。

本文以Kylafis和Loreau提出的捕食-競爭系統(tǒng)為基礎(chǔ), 建立了捕食-競爭系統(tǒng)的反應(yīng)擴散模型。通過數(shù)值模擬研究, 探討了生態(tài)位構(gòu)建對于捕食-競爭中各物種變化趨勢及其空間擴散的影響。研究結(jié)果表明, 生態(tài)位構(gòu)建行為在不破壞系統(tǒng)穩(wěn)定性的同時,明顯提升了構(gòu)建物種的競爭力, 使其在不利環(huán)境中依然能與競爭物種共處; 同時, 生態(tài)位構(gòu)建行為加強了物種的空間擴散能力, 使得物種能夠迅速遷徙到更適合自己生存的空間, 成為捕食-競爭系統(tǒng)中的穩(wěn)定生存物種。這一研究結(jié)果從數(shù)值角度闡釋了生態(tài)位構(gòu)建對于物種數(shù)量以及空間擴展的影響, 拓展了Kylafis和Loreau的研究結(jié)果。

2 模型

2.1捕食-競爭系統(tǒng)

生態(tài)位構(gòu)建是指有機體對其環(huán)境上自然選擇源的改變, 而在生態(tài)時間尺度上, 其重點強調(diào)物種的構(gòu)建行為對其自身所在環(huán)境的改變和影響。在物種的實際生存中, 物種數(shù)量主要受到兩個因素的限制:分別是有限的自然資源和捕食者自身對于物種的捕食現(xiàn)象。據(jù)此, 可以建立一個基本的捕食-競爭系統(tǒng)(系統(tǒng)1): 該系統(tǒng)包含A和B兩個競爭物種和一個捕食者P。為了引入生態(tài)位構(gòu)建的作用, 在該模型中同時考慮自然資源R的動態(tài), 其中兩競爭物種通過攝取自然資源R而生存; 捕食者P則以捕食兩競爭物種作為自身生存的前提。模型具體表示如下(系統(tǒng)1):

其中i=A, B,im和ic分別表示捕食者對物種i 的捕食效率和轉(zhuǎn)化效率,if和ia表示競爭物種i對資源的攝取效率及轉(zhuǎn)化效率,pd,id則代表捕食者P和競爭物種i的自然死亡率, k代表自然資源自身的轉(zhuǎn)化效率, 而S則代表自然資源的供應(yīng)總量。

生態(tài)位構(gòu)建的作用可以通過兩方面來體現(xiàn), 一種是直接的方式, 即物種通過構(gòu)建有利于自身生存的環(huán)境以加強自身對于環(huán)境的適應(yīng)性, 另一種則是間接的方式, 即物種通過構(gòu)建抑制競爭物種生存的環(huán)境以降低競爭物種對于環(huán)境的適應(yīng)性。本文重點考慮前者, 在系統(tǒng)(1)的基礎(chǔ)上針對物種A, 通過增加其抵抗捕食者的被捕食能力, 建立具有生態(tài)位構(gòu)建作用的捕食-競爭系統(tǒng), 模型表示如下(系統(tǒng)2):

其中參數(shù)Pn為競爭物種A對捕食者的生態(tài)位構(gòu)建強度, 當(dāng)pn=0時, 系統(tǒng)(2)退化為系統(tǒng)(1), 其他參數(shù)含義同系統(tǒng)(1)。

為簡化模擬與討論, 考慮到物種對資源的轉(zhuǎn)化效率對于其在生存競爭中所產(chǎn)生的影響甚微, 假設(shè)各競爭物種對于自然資源R的轉(zhuǎn)化效率以及捕食者P對于各競爭物種的轉(zhuǎn)化效率均為理想狀態(tài), 即設(shè)定參數(shù)ci, ai, k, s均取值為1; 同時, 為了更好的反映各物種間的關(guān)系, 進一步假設(shè)捕食-競爭系統(tǒng)中競爭物種A和物種B具有如下層次特征: 即物種A相對于B而言, 具有較強的對自然資源的攝取能力,但更易受到捕食者P的捕食; 這兩者的競爭特征反映在系數(shù)上則體現(xiàn)為: mA＞m B, fA＞fB。在數(shù)值模擬中, 通過固定參數(shù)fA、fB及mB, 選擇變化mA和nP的方法, 對上述兩個系統(tǒng)的穩(wěn)定性及各物種數(shù)量的變化趨勢進行對比分析, 以清楚地了解生態(tài)位構(gòu)建作用對捕食-競爭系統(tǒng)的影響。

2.2捕食-競爭系統(tǒng)反應(yīng)擴散模型

從空間生態(tài)學(xué)角度上講, 環(huán)境的異質(zhì)性是導(dǎo)致物種形成及分布的重要因素, 也是空間模擬的主要原因。其中, 偏微分方程及反應(yīng)擴散方程一直是數(shù)學(xué)生態(tài)學(xué)研究的主要手段。經(jīng)典的常微分方程描述生物系統(tǒng)的時變動態(tài), 而偏微分方程則可以給出空間上每一點的時間動態(tài)及密度分布[27]。因此, 為了進一步探討生態(tài)位構(gòu)建作用對物種空間分布動態(tài)的影響, 我們在上述系統(tǒng)(1)和(2)的基礎(chǔ)上分別建立了具有空間擴散項的偏微分方程。其中, 反應(yīng)項即為上述系統(tǒng)(1)和(2)的主體部分, 描述物種自身的增長動態(tài); 而擴散項則由布朗隨機運動描述, 視有機體在空間生境中做布朗隨機運動且運動的速率是恒定的無方向差別量。依據(jù)假設(shè), 相應(yīng)的擴散模型分別為基本捕食-競爭系統(tǒng)的反應(yīng)擴散模型(系統(tǒng)3)和具有生態(tài)位構(gòu)建的捕食-競爭系統(tǒng)的反應(yīng)擴散模型(系統(tǒng)4), 表示如下:

其中PD、AD、BD分別表示物種P、A、B物種的擴散系數(shù)。首先, 為了單純觀察生態(tài)位構(gòu)建對于物種擴散的影響, 將擴散系數(shù)統(tǒng)一設(shè)置為1; 其次, 進一步考慮到有機體自身用于生長、生命維持以及繁殖的時間、能量和其他資源總是有限的, 其必須在這些可能的需求方面分配有限資源[28]。對于該系統(tǒng),需求方面的沖突會導(dǎo)致生態(tài)位構(gòu)建能力與擴散能力之間的妥協(xié)。因此, 假設(shè)AD+Pn=1。選取一個大小為nn×(n = 300)的空間顯含模式的正方形網(wǎng)格區(qū)域,其邊界條件為零流量輸出邊界條件。通過空間模擬方法, 對兩種不同擴散因子假設(shè)下的模型中分別進行對照分析, 從而了解生態(tài)位構(gòu)建對于捕食-競爭系統(tǒng)中各物種空間擴散的影響。

3 結(jié)果與分析

3.1捕食-競爭系統(tǒng)的穩(wěn)定性分析

為了更好的反映生態(tài)位構(gòu)建對于捕食—競爭系統(tǒng)中各物種數(shù)量變化趨勢的影響, 首先利用常微分方程的穩(wěn)定性理論分析系統(tǒng)的非平凡穩(wěn)定平衡點及平衡條件, 以確保后續(xù)的討論是在系統(tǒng)中捕食者P和競爭物種A、B共存的條件下進行的。

類似的, 通過求解系統(tǒng)(2)的平衡點, 不難得出系統(tǒng)(2)穩(wěn)定平衡的限制條件為: (2.1)。通過對比兩系統(tǒng)的限制條件, 發(fā)現(xiàn)相較于系統(tǒng)(1)的穩(wěn)定條件, 系統(tǒng)(2)的差異在于構(gòu)建系數(shù)Pn的引入, 因此,下文將著重分析該參數(shù)對系統(tǒng)穩(wěn)定性及物種間數(shù)量關(guān)系的影響。

3.1.1 捕食-競爭系統(tǒng)(系統(tǒng)1)的時間動態(tài)

數(shù)值模擬結(jié)果顯示: 捕食者P, 競爭物種A、B動態(tài)呈現(xiàn)阻尼振蕩趨勢(見圖1a), 三者的相空間曲線最終螺旋匯聚于內(nèi)部平衡點(見圖1b)。進一步觀察系統(tǒng)(1)的動態(tài)變化趨勢(見圖2)發(fā)現(xiàn): 當(dāng)mA=0.21時, 物種A處于優(yōu)勢并漸趨穩(wěn)定, 物種B處于劣勢并漸趨滅亡; 當(dāng)mA=0.2119時, 兩競爭物種維持平衡狀態(tài); 當(dāng)mA=0.4時, 競爭物種A迅速消亡, 物種B取代物種A成為競爭中的穩(wěn)定生存物種。以上動態(tài)變化趨勢顯示與的接近程度影響著兩物 種的數(shù)量變化, 當(dāng)與十分接近時, 物種A在競爭中處于優(yōu)勢, 當(dāng)遠大于時, 物種B 取代物種A, 成為競爭中的穩(wěn)定生存物種。

圖1 系統(tǒng)(1)的三物種動態(tài)圖(圖1a)及相空間圖(圖1b)Fig. 1 The diagrams of species temporal dynamics (a) and phase plane (b) of system (1)

圖2 系統(tǒng)(1)的物種數(shù)量動態(tài)圖, (a)、(b)、(c)所對應(yīng)的mA分別為0.21, 0.2119和0.4Fig. 2 The temporal dynamics of three species in system (1); the parameters of mAin (a), (b) and (c) are 0.21, 0.2119 and 0.4.

3.1.2 具有生態(tài)位構(gòu)建項模型(系統(tǒng)2)的時間動態(tài)

數(shù)值模擬結(jié)果顯示, 系統(tǒng)(2)的相空間圖(見圖3)及物種數(shù)量動態(tài)趨勢圖(見圖4)與系統(tǒng)(1)類似, 均為圍繞內(nèi)部平衡點的阻尼振蕩動態(tài)。仔細觀察, 不難發(fā)現(xiàn): 當(dāng)Pn=0.1885時, 物種A和B的數(shù)量水平相當(dāng); 而當(dāng)構(gòu)建系數(shù)繼續(xù)增加時, 系統(tǒng)(2)出現(xiàn)了異于系統(tǒng)(1)的情形: 物種A取代物種B成為競爭中的優(yōu)勢方, 物種B漸趨消亡。

對比上述兩組結(jié)果, 可以看出生態(tài)位構(gòu)建系數(shù)nP的引入并未改變捕食-競爭系統(tǒng)的穩(wěn)定性, 但是生態(tài)位構(gòu)建行為提升了構(gòu)建物種A的競爭力, 使該物種不僅能在不利環(huán)境中(Am的增加表明競爭物種A所面臨的被捕食壓力越來越大)維持與物種B共處的平衡狀態(tài), 而且在不斷增加構(gòu)建系數(shù)的情況下,取代物種B成為競爭中的穩(wěn)定生存物種。

3.2捕食-競爭系統(tǒng)的空間動態(tài)

對于反應(yīng)擴散模型的模擬, 根據(jù)其數(shù)值模擬對于時間空間步長的依賴性和敏感度, 我們合理的選擇了較小的時空步長以盡可能降低模擬誤差。結(jié)果顯示, 對于系統(tǒng)3, 當(dāng)其不具有構(gòu)建因素時, 捕食者以及競爭物種主要聚集在初始生存區(qū)域, 并未發(fā)生明顯的空間擴散行為(圖5)。當(dāng)競爭物種A具有生態(tài)位構(gòu)建行為時, 各物種的生存區(qū)域發(fā)生了明顯變化(圖6): 特別是具有構(gòu)建行為的物種A, 因為構(gòu)建行為提升了其躲避捕食者P的能力, 所以其表現(xiàn)出離開原有的生存區(qū)域、并遷移至生存風(fēng)險較低的生存區(qū)域的空間擴張行為; 對于不具有構(gòu)建行為的競爭物種B, 其生存環(huán)境與能力并未發(fā)生改變, 但其相對于構(gòu)建物種A, 其生存能力在減弱, 因此其生存區(qū)域逐漸收縮; 而捕食者P則因為競爭物種A的遷移也相應(yīng)的擴大了自己的捕食范圍, 進而擴大了自身的生存區(qū)域。

圖3 系統(tǒng)2的相空間圖, (a)、(b)、(c)所對應(yīng)的nP分別為0.1, 0.1885和0.2Fig. 3 The diagram phase plane of system (2); the parameters of npin (a), (b) and (c) are 0.1, 0.1885 and 0.2.

圖4 系統(tǒng)2的物種數(shù)量動態(tài)圖, (a)、(b)、(c)所對應(yīng)的nP分別為0.1, 0.1885和0.2Fig. 4 The temporal dynamics of three species in system (2); the parameters of npin (a), (b) and (c) are 0.1, 0.1885 and 0.2.

圖5 模型3中各物種的空間分布圖, 參數(shù)為: fA=0.2, mA=0.4, fB=0.1, mB=0.1, T=3000.Fig. 5 The spatial distribution of three species in model (3); the parameters are fA=0.2, mA=0.4, fB=0.1, mB=0.1, T=3000.

可以注意到, 在其他環(huán)境參數(shù)不變的情況下,一方面, 若不考慮生態(tài)位構(gòu)建與擴散能力間的妥協(xié)關(guān)系, 構(gòu)建系數(shù)的變化也會影響系統(tǒng)中各物種的空間分布變化。模擬結(jié)果(圖7,圖8和圖9)表明: 當(dāng)構(gòu)建系數(shù)逐漸加強時, 在物種B的生存區(qū)域逐漸收縮的同時, 構(gòu)建物種A和捕食者P的生存區(qū)域及空間密度發(fā)生顯著變化.其中構(gòu)建物種A的生存區(qū)域逐漸擴大, 空間密度也在逐漸增加, 當(dāng)達到一定閾值(Pn=0.1885)后, 物種A出現(xiàn)了回溯遷移的空間擴展行為; 另一方面, 對于捕食者P而言, 當(dāng)構(gòu)建系數(shù)較小時, 其生存區(qū)域以及空間密度都逐漸提升, 但隨著構(gòu)建系數(shù)的逐漸增加, 其生存區(qū)域開始出現(xiàn)收縮, 邊界空間密度開始降低。

圖6 模型4中各物種的空間分布圖, 參數(shù)值Pn=0.1, 其他參數(shù)同圖5.Fig. 6 The spatial distribution of three species in model (4);Pn=0.1, other parameters are same as those in Fig.5.

圖7 模型4中物種P的空間區(qū)域在不同Pn下的快照, (a)、(b)、(c)中的Pn分別為0.1, 0.1885, 0.3Fig. 7 Snapshots of species P distribution with differentPn; the parameters in (a), (b) and (c) are 0.1, 0.1885, 0.3.

圖8 模型4中物種A的空間區(qū)域在不同Pn下的快照,Pn參數(shù)同圖7.Fig. 8 Snapshots of species A distribution with differentPn; the parameters are same as those in Fig.7.

圖9 模型4中物種B的空間區(qū)域在不同Pn下的快照,Pn參數(shù)同圖7.Fig. 9 Snapshots of species B distribution with differentPn; the parameters are same as those in Fig.7.

同時, 物種的構(gòu)建行為對于其空間擴散的作用會受到trade-off效應(yīng)的制約, 從而影響著系統(tǒng)中各物種的空間分布, 模擬結(jié)果(見圖10, 圖11和圖12)顯示: 當(dāng)模型4受到trade-off效應(yīng)影響時, 對于不具有構(gòu)建行為的物種B而言, 其空間區(qū)域及空間密度較之原模型(見圖7)明顯減弱; 而對于構(gòu)建物種A及捕食者P而言, 雖然二者的空間密度并未發(fā)生變化,但它們的空間分布區(qū)域基本維持在Pn=0.5時的區(qū)域面積, 并未出現(xiàn)進一步的擴散。由此可以看出, 構(gòu)建行為雖然會促進物種空間擴散, 但由于物種自身能力所存在的trade-off效應(yīng), 這種促進作用會隨著構(gòu)建系數(shù)的增加而逐漸降低, 使得物種的空間區(qū)域維持在一定的平衡狀態(tài), 不會出現(xiàn)無限制的空間擴散。

圖10 具有trade-off效應(yīng)的模型4中物種P的空間區(qū)域在不同Pn下的快照,Pn參數(shù)同圖7.Fig. 10 Snapshots of species P distribution with trade-off effects; the parametersPn are same as those in Fig.7.

圖11 具有trade-off效應(yīng)的模型4中物種A的空間區(qū)域在不同Pn下的快照,Pn參數(shù)同圖7.Fig. 11 Snapshots of species A distribution with trade-off effects; the parameters are same as those in Fig.7.

圖12 具有trade-off效應(yīng)的模型4中物種B的空間分布在不同Pn下的快照,Pn參數(shù)同圖7.Fig. 12 Snapshots of species B distribution with trade-off effects; the parameters are same as those in Fig.7.

4 討論與總結(jié)

本文通過對具有生態(tài)位構(gòu)建行為的捕食-競爭系統(tǒng)和考慮到空間擴散因素而建立的反應(yīng)擴散模型的數(shù)值模擬和比較分析, 研究了物種的構(gòu)建行為對捕食-競爭系統(tǒng)演化動態(tài)和空間分布模式的影響。首先, 當(dāng)物種不存在生態(tài)位構(gòu)建行為時, 對自然資源的攝取能力較強但更易受到捕食的物種A的數(shù)量變化趨勢是前期優(yōu)勢、后期劣勢并漸趨消亡; 而攝取能力較弱且不易被捕食的物種B的數(shù)量變化趨勢則與之恰好相反, 并最終取代競爭物種A; 當(dāng)競爭物種A具有生態(tài)位構(gòu)建行為時, 系統(tǒng)的穩(wěn)定性并未發(fā)生改變, 但明顯的, 物種A通過構(gòu)建行為提高了自身競爭力, 使得其在不利環(huán)境中依舊能與競爭物種B共存, 甚至取而代之成為系統(tǒng)的穩(wěn)定生存物種。因此,生態(tài)位構(gòu)建行為也是導(dǎo)致物種穩(wěn)定共存的因素之一,這一結(jié)果得到之前多項研究的支持[30–31], 特別是惠蒼等人2004年[9]的研究表明生態(tài)位構(gòu)建的機制可能是解釋許多群落中競爭共存的潛在模式。本文的研究結(jié)果再一次說明有機體的生態(tài)位構(gòu)建行為理應(yīng)在生態(tài)學(xué)研究中備受重視[7]。

本文的重點在于考慮了具有擴散因子的生態(tài)位構(gòu)建物種在競爭與捕食雙重作用下的反應(yīng)擴散方程。嚴(yán)格意義上來講, 這是首次在生態(tài)位構(gòu)建的研究工作中使用偏微分方程這一工具。本研究的模擬結(jié)果充分展示了構(gòu)建行為對具有擴散因子的物種的空間分布動態(tài)的影響。其中, 具有生態(tài)位構(gòu)建作用的物種A會離開原有的生存區(qū)域, 并因擴散因素而遷移至生存風(fēng)險更低的空間區(qū)域; 捕食者P則由于物種A的遷徙而相應(yīng)的擴大了自身的捕食區(qū)域, 進而也擴展了自身生存區(qū)域。尤其是隨著構(gòu)建系數(shù)的增加, 物種A和捕食者P的生存區(qū)域及空間密度發(fā)生了明顯的變化。其中構(gòu)建物種A的生存區(qū)域逐漸擴大, 當(dāng)達到一定閾值后, 物種A出現(xiàn)了回溯遷移的空間擴展行為, 此時, 捕食者P生存區(qū)域也開始出現(xiàn)收縮, 且邊界空間密度降低。在這一結(jié)果中, 閾值的出現(xiàn)和空間擴張行為的回溯現(xiàn)象是源于具有構(gòu)建行為的有機體與環(huán)境之間的反饋機制。這也是生態(tài)位構(gòu)建理論的核心, 即有機體與環(huán)境之間相輔相成、相互影響與制約的正反饋作用能夠產(chǎn)生豐富的進化現(xiàn)象與生態(tài)學(xué)動態(tài)[5]。

Trade-off關(guān)系是指有機體由于可供分配的資源有限, 特定兩個特征間的相互作用, 即提高一種特征的優(yōu)勢的時候, 另一種特征的優(yōu)勢將降低?，F(xiàn)已開展了諸多關(guān)于trade-off關(guān)系的理論研究, 例如均質(zhì)斑塊環(huán)境上的集合種群競爭共存的關(guān)鍵是侵占能力與競爭能力之間的妥協(xié)[32]。Tilman的模型提供了一個由固著個體組成的許多群落中高物種豐富度的完美解釋。另外, 到目前為止, 促進穩(wěn)定共存的因子已經(jīng)被很好地研究和討論過。例如, 空間細分和局部擴散會導(dǎo)致奇異的空間動態(tài)并使捕食系統(tǒng)動態(tài)穩(wěn)定[31]。侵占-競爭妥協(xié)會形成競爭者與逃逸者的穩(wěn)定共存[32]?；萆n等人2004年[9]的研究則說明空間尺度的物種共存可能得益于競爭力、侵占力以及生態(tài)位構(gòu)建力之間巧妙的妥協(xié)。而本文基于生態(tài)位構(gòu)建能力和擴散能力間妥協(xié)關(guān)系的假設(shè), 研究結(jié)果表明:生態(tài)位構(gòu)建行為對于物種空間擴散的促進作用并不會導(dǎo)致物種無限制的空間擴張, 而是隨著構(gòu)建系數(shù)的增加而逐漸降低, 進而維持物種空間區(qū)域的平衡狀態(tài)。

Kylafis和Loreau于2011年[13]曾對捕食-競爭系統(tǒng)進行了初步的理論研究, 他們認為: 生態(tài)位構(gòu)建有助于增加構(gòu)建物種抵抗捕食者的能力, 維持構(gòu)建物種在不利條件下與競爭物種共處的平衡狀態(tài), 這與本文的數(shù)值模擬結(jié)果一致。進一步對系統(tǒng)中各物種添加擴散項后的模擬結(jié)果證實, 構(gòu)建行為還有助于提升構(gòu)建物種的空間擴展能力, 擴展其生存區(qū)域。這一研究將原有的研究從之前單一的物種數(shù)量變化擴展至生存區(qū)域的變化, 對于探討生態(tài)位構(gòu)建行為對物種生態(tài)分布及種群變化趨勢的影響有著重要意義。由于本文的模型是在理想化情況下構(gòu)建的,并未考慮實際物種共存中存在的Allee效應(yīng)等問題, Allee效應(yīng)會對具有構(gòu)建行為的捕食-競爭系統(tǒng)中各物種的數(shù)量變化及其空間擴展造成怎樣的影響有待進一步研究。

致謝本文得到國家自然科學(xué)青年基金資助項目-生態(tài)位構(gòu)建對環(huán)境梯度下物種分布格局的作用(31100308)的資助。

參考文獻

[1] LALAND K N, ODLING-SMEE F J, FELDMAN M W. The evolutionary consequences of niche construction: a theoretical investigation using two-locus theory[J]. Journal of Evolutional Biology, 1996, 9: 293–316.

[2] DAY R L, LALAND K N, ODLING-SMEE F J. Rethinking Adaptation. the niche-construction perspective[J]. Perspectives in Biology and Medicine, 2003, 46: 80–95.

[3] ODLING-SMEE F J, LALAND K N, FELDMAN M W. Niche construction: the neglected process of evolution[M]. Princeton University Press, 2003.

[4] LEIBOLD M A. A graphical model of keystone predators in food webs trophic regulation of abundance, incidence, and diversity patterns in communities[J]. American Naturalist, 2011, 147: 784–812.

[5] ODLING-SMEE F J, ERWIN D H, PALKOVACS E P, et al. Niche Construction Theory: A Practical Guide for Ecologists[J]. The Quarterly Review of Biology, 2013, 88(1): 3–28.

[6] O’BRIEN M J, LALAND K N. Genes, culture, and agriculture: an example of human niche construction[J]. Current. Anthropology, 2012, 53: 434–470.

[7] MATTHEWS B, MEESTER L D, JONES C G, et al. Under niche construction: an operational bridge between ecology, evolution, and ecosystem science[J]. Ecological Monographs, 2014, 84(2): 245–263.

[8] PALKOVACS E P, KINNISON M T, CORREA C, et al. Fates beyond traits: ecological consequences of humaninduced trait change[J]. Evolutionary Applications, 2012, 5: 183–191.

[9] HUI Cang, LI Zizheng, YUE Dongxia. Metapopulation dynamics and distribution, and environmental heterogeneity induced by niche construction[J]. Ecological Modelling, 2004, 177: 107–188.

[10] HAN Xiaozhuo, LI Zizheng, HUI Cang, et al. Polymorphism maintenance in a spatially structured population: two-locus genetics model of niche construction[J]. Ecological Modelling, 2006, 192: 160–174.

[11] HAN Xiaozhuo, LI Zizheng, ZHANG Feng, et al. Effect of niche construction on metapopulation dynamics by modifying destroyed habitat[J]. Journal of Lanzhou University, 2006, 43(1): 56–64.

[12] KYLAFIS G, LOREAU M. Ecological and evolutionary consequences of niche construction for its agent[J]. Ecology Letters, 2008, 11: 1072–1081.

[13] KYLAFIS G, LOREAU M. Niche construction in the light of niche theory[J]. Ecology Letters, 2011, 14: 82–90.

[14] HAN Xiaozhuo, HUI Cang. Niche Construction on Environmental Gradients: The Formation of Fitness Valley and Stratified Genotypic Distributions [J]. PLoS One, 2014, 9(6): e99775.

[15] BUSER C C, NEWCOMBANNE R D, GASKETT A C, et al. Niche construction initiates the evolution of mutualistic interactions[J]. Ecology Letters, 2014, 17: 1257–1264

[16] BRODAW N, BUSING P T. Niche versus chance and tree diversity in forest gaps[J]. Trends in Ecology & Evolution, 2000, 15: 183–188.

[17] HUI Cang, ROURA-PASCUAL N, BROTONS L, et al. Flexible dispersal strategies in native and non-native ranges: environmental quality and the ‘good-stay, bad-disperse’rule[J]. Ecography, 2012, 35: 1024–32.

[18] HURFORD J R. Niche construction, co-evolution and domain specificity[J]. Behavioral and brain Sciences, 2008, 31(5): 526–526.

[19] IHARA Y, FELDMAN M W. Cultural niche construction and the evolution of small family size[J]. Theoretical Population Biology, 2004, 65: 105–111.

[20] BERTHOULY-SALAZAR C, VAN RENSBURG B J, LEROUX J J, et al. Spatial sorting drives morphological variation in the invasive bird, Acridotheres tristis[J]. PLoS ONE, 2012, 7: e38145.

[21] ERWIN D H. Macroevolution of ecosystem engineering, niche construction and diversity[J]. Trends in Ecology & Evolution, 2004, 23(6): 304–311.

[22] HORAK J, HUI Cang, ROURA-PASCUAL N, et al. Changing roles of propagule, climate, and land use during extralimital colonization of a rose chafer beetle[J]. Naturwissenschaften, 2013, 100: 327–36.

[23] HUI Cang, LI Zizheng, YUE Dongxia. Metapopulation dynamics and distribution, and environmental heterogeneity induced by niche construction[J]. Ecological Modelling, 2004, 177: 107–188.

[24] HAN Xiaozhuo, LI Zizheng, ZHANG Yanyu. Effects of time-lagged niche construction on metapopulation dynamics and environmental heterogeneity[J]. Applied Mathematics & Computation, 2009, 215: 449–458.

[25] PALMER M E, FELDMAN M W. Survivability is more fundamental than evolvability[J]. PLoS One, 2012, 7: e38025.

[26] HAN Xiaozhuo, LI Zizheng, HUI Cang, et al. Theory of niche construction and application of its two-locus population genetic model[J]. Acta Botanica Boreali-Occidentalia Sinica, 2004, 24(3): 558–562.

[27] HUI Cang. Simulated research on spatial Ecology: the ecological and evolutionary responses of metapopulation to environmental changes[D]. Lanzhou: Lanzhou University, 2004: 7–11.

[28] SILVERTOWN J W, DOUST J L. Introduction to plant population biology[M]. Blackwell, Oxford. 1993.

[29] 崔麗娟, 李偉, 趙欣勝, 等. 采砂跡地型濕地恢復(fù)過程中優(yōu)勢種群生態(tài)位研究[J]. 生態(tài)科學(xué), 2013, 32(1): 073–077.

[30] 袁叢軍, 喻理飛, 安明態(tài), 等. 巖生紅豆天然林主要樹種生態(tài)位研究[J]. 生態(tài)科學(xué), 2013, 32(4): 487–493.

[31] HASSELL M P, COMINS H N, MAY R M. Spatial structure and chaos in insect population dynamics[J]. Nature, 1991, 353: 255–258.

[32] TILMAN D. Competition and biodiversity in spatially structured habitats[J]. Ecology, 1994, 75: 2–16.

The effects of niche construction on prey-competition system

MO Changjian1,2, HAN Xiaozhuo1,*
1.School of Applied Mathematics,Guangdong University of Technology,Guangzhou510512,China
2.Foundation Department,Guangdong Songshan Polytechnic College,Shaoguan512000,China

Niche construction is the process whereby an organism improves its environment to enhance its growth and persistence. Based on a prey-competition system with two competitive consumers and one predator, we included direct and indirect ways of niche construction by consumers and also the factor of spatial diffusion. Through numerical simulation, the effect of niche construction on temporal dynamics and spatial distribution of different species in the system were explored. Our results indicated that the niche construction did not destroy the stability of the prey-competition system. Niche- constructing species, however, can enhance its own competitive strength and also maintain the system balance in adverse environments. In addition, the niche construction has the potential to improve a diffusivity of the niche-constructing species, which promotes them to immigrate for more suitable habitats and become the stable species in the system. Because of the limit resource, the trade-off relationship between the niche construction and the diffusion of species imposed restrictions on the distribution range, inhabiting species unlimited expansion.

niche construction; reaction-diffusion equations; trade-off; equilibrium state

10.14108/j.cnki.1008-8873.2016.05.011

Q141

A

1008-8873(2016)05-073-09

莫昌健, 韓曉卓. 生態(tài)位構(gòu)建作用對捕食‐競爭系統(tǒng)的影響[J]. 生態(tài)科學(xué), 2016, 35(5): 73-81.

MO Changjian, HAN Xiaozhuo. The effects of niche construction on prey-competition system[J]. Ecological Science, 2016, 35(5): 73-81.

2015-03-12;

2015-04-25

國家自然科學(xué)青年基金項目(31100308)

莫昌健(1986—), 男, 廣東韶關(guān)人, 碩士, 主要從事數(shù)學(xué)生態(tài)學(xué)研究

*通信作者: 韓曉卓(1978—), 女, 青海西寧人, 博士, 副教授, 主要從事數(shù)學(xué)生態(tài)學(xué)研究, E-mail: hanxzh@gdut.edu.cn