破碎景觀建立保護(hù)區(qū)面積-數(shù)量模式界定

劉昊奇,呂光輝

1 新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,烏魯木齊 830046 2 新疆綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,烏魯木齊 830046 3 新疆大學(xué)干旱生態(tài)環(huán)境研究所, 烏魯木齊 830046

由于動(dòng)植物生存環(huán)境的改變,很多物種的生存變得日益困難[1],更嚴(yán)重的是大量珍稀野生動(dòng)植物瀕臨滅絕,如東北虎(Pantheratigrisaltaica)、大熊貓(Ailuropodamelanoleuca)和銀杉(Cathayaargyrophylla)等,因而保護(hù)瀕危物種成為當(dāng)前研究的熱點(diǎn)。自然保護(hù)區(qū)作為保護(hù)瀕危物種的一種有效工具[2-4],引起了研究者的廣泛關(guān)注。在設(shè)計(jì)自然保護(hù)區(qū)時(shí),有幾個(gè)問題特別引人關(guān)注,例如：應(yīng)該將可利用的生境建設(shè)成少量大保護(hù)區(qū)還是大量小保護(hù)區(qū)？是否應(yīng)在各保護(hù)區(qū)間建立廊道,將孤立的保護(hù)區(qū)連接成網(wǎng)絡(luò)[5-6]？第一個(gè)問題就是著名的FLOMS(few large or many small reserves, 少量大保護(hù)區(qū)還是大量小保護(hù)區(qū))問題,或稱之為SLOSS(single large or several small reserves, 單一大保護(hù)區(qū)還是多個(gè)小保護(hù)區(qū))問題,本文針對(duì)此問題進(jìn)行探究。

早在20世紀(jì)70年代中期,研究者就提出了6條自然保護(hù)區(qū)的設(shè)計(jì)準(zhǔn)則,其中之一：當(dāng)保護(hù)區(qū)面積固定時(shí),單一的大保護(hù)區(qū)模式更有助于物種生存[5,7]。但隨后這些結(jié)論受到了質(zhì)疑,FLOMS問題隨之成為研究的熱點(diǎn)[8-11]。部分研究表明少量大保護(hù)區(qū)對(duì)物種生存有利[12-15];部分研究則認(rèn)為大量小保護(hù)區(qū)更有利[10- 11,16-18];一些研究者認(rèn)為中等程度的保護(hù)區(qū)數(shù)量最有利[19];還有研究者認(rèn)為最優(yōu)保護(hù)區(qū)數(shù)量隨保護(hù)區(qū)環(huán)境承載量的增大而增加[20]。盡管已經(jīng)有了大量研究結(jié)果,但由于研究目標(biāo)和方法的不同,沒有得出統(tǒng)一的結(jié)論,綜合所有的研究結(jié)果只能得到：最優(yōu)保護(hù)區(qū)數(shù)量由物種自身特征和所處的環(huán)境決定,例如：環(huán)境承載量、個(gè)體擴(kuò)散和環(huán)境擾動(dòng)[11,19- 20]。然而,還有很多因子也可能對(duì)FLOMS問題有重要影響,但之前的研究卻未深入研究,例如：生境喪失及其與個(gè)體擴(kuò)散的共同作用。因而探討這些因子可幫助研究者了解物種自身特征和所處的環(huán)境如何影響FLOMS問題,并促進(jìn)此問題的最終解決。

生境喪失不僅是物種滅絕的首要因素[21],還造成現(xiàn)今很難找到一片未被破壞的景觀來建設(shè)保護(hù)區(qū),因而生境喪失已成為影響保護(hù)區(qū)建設(shè)的重要因素。研究表明：生境喪失可通過作用于個(gè)體擴(kuò)散、種群增長(zhǎng)和環(huán)境擾動(dòng)來影響FLOMS問題[20,22- 23]。個(gè)體擴(kuò)散也可能對(duì)FLOMS問題有重要影響。個(gè)體擴(kuò)散可通過“援救效應(yīng)”幫助物種擺脫滅絕[24],但個(gè)體擴(kuò)散也可能帶入疾病或寄生蟲,對(duì)物種造成負(fù)面影響[25]。關(guān)于個(gè)體擴(kuò)散,即使不考慮帶入疾病和寄生蟲,也存在著截然相反的觀點(diǎn)——大部分研究認(rèn)為擴(kuò)散對(duì)物種有利[26- 27],其他的則認(rèn)為不利[28]。出現(xiàn)兩個(gè)截然相反的結(jié)論的原因可能在于擴(kuò)散過程中的個(gè)體死亡。若死亡率低,擴(kuò)散可通過“援救效應(yīng)”幫助物種擺脫滅絕;否則,較高的死亡率會(huì)導(dǎo)致擴(kuò)散帶來的副面作用大于正面作用。

盡管已經(jīng)有了如此多的研究,但回顧之前的研究可發(fā)現(xiàn)幾個(gè)缺點(diǎn)。其一,盡管被破壞生境的空間分布方式是非隨機(jī)的[29- 30],但之前的研究均假定其為隨機(jī)的[20,22- 23]。其二,個(gè)體擴(kuò)散可能與生境喪失共同作用于物種,從而影響FLOMS問題,但涉及到該問題的研究很少。

1 模型

模型描述根據(jù)Overview, Design concepts, Details protocol(ODD)[31- 32]基于代理模型(agent-based models, ABM)的描述方式展開。

1.1 目的

在一個(gè)已經(jīng)遭受過生境喪失的景觀上,選取一片正方形的區(qū)域,并調(diào)整此區(qū)域的面積以保證未被破壞的生境的面積為一個(gè)固定常數(shù),探討生境喪失和個(gè)體擴(kuò)散對(duì)FLOMS問題的影響。

1.2 載體和狀態(tài)變量

所選取的區(qū)域包含兩種類型的生境：被破壞的生境和未被破壞的生境,且所有未被破壞的生境將被建設(shè)成自然保護(hù)區(qū)。因此,模型共有3種載體—被破壞的生境,保護(hù)區(qū)和生存于保護(hù)區(qū)中的被保護(hù)物種。每個(gè)保護(hù)區(qū)對(duì)應(yīng)一個(gè)環(huán)境最大可持續(xù)承載量和空間位置;每個(gè)保護(hù)區(qū)都存在一個(gè)局部種群,局部種群的某些個(gè)體可在保護(hù)區(qū)間擴(kuò)散,這些擴(kuò)散的個(gè)體將各局部種群聯(lián)系起來形成一個(gè)集合種群;每個(gè)被破壞的生境也占據(jù)著一個(gè)空間位置。

1.3 過程綜述

一個(gè)離散時(shí)間和空間的模型被用于描述整個(gè)生態(tài)過程(圖1)。每一步：環(huán)境可持續(xù)承載量增長(zhǎng);局部種群大小增長(zhǎng);每個(gè)保護(hù)區(qū)都可能受到擾動(dòng)的影響,擾動(dòng)率為μ;局部種群中的某些個(gè)體可能擴(kuò)散到相鄰的保護(hù)區(qū)。如果在某一步,被保護(hù)物種的個(gè)體總數(shù)少于a,那么認(rèn)為此物種的個(gè)體數(shù)量太少,以至于已經(jīng)生態(tài)滅絕,并且終止整個(gè)過程。

圖1 過程,次序安排和偽代碼Fig.1 Processes, scheduling and pseudo-code

1.4 概念設(shè)計(jì)

因?yàn)榭杀唤ㄔO(shè)成保護(hù)區(qū)的生境的面積是恒定常數(shù),所以所有保護(hù)區(qū)的總面積是恒定常數(shù)。盡管在現(xiàn)實(shí)中,各保護(hù)區(qū)的面積、形狀和環(huán)境可持續(xù)承載量以及各局部種群的大小不會(huì)完全相同,忽略這些差異,還假定各被破壞的生境的面積和形狀與保護(hù)區(qū)相同。

依據(jù)被保護(hù)物種的滅絕概率決定建立大量小保護(hù)區(qū)還是少量大保護(hù)區(qū),最優(yōu)保護(hù)區(qū)數(shù)量為被保護(hù)物種的滅絕概率達(dá)到最低時(shí)對(duì)應(yīng)的保護(hù)區(qū)數(shù)量,此時(shí)還可得到單個(gè)保護(hù)區(qū)的面積,因?yàn)樗斜Ｗo(hù)區(qū)的總面積是固定的常數(shù),據(jù)此可以決定建設(shè)大量小保護(hù)區(qū)還是少量大保護(hù)區(qū)。將未被破壞的生境劃分成N2個(gè)面積相等的保護(hù)區(qū)(N2的取值范圍從1變化到很多),對(duì)于每個(gè)N2各進(jìn)行1000次獨(dú)立重復(fù)模擬實(shí)驗(yàn),并且從實(shí)驗(yàn)結(jié)果得到物種的滅絕概率。物種的滅絕概率定義為：在模擬過程中發(fā)生物種滅絕的那些模擬的次數(shù)與總模擬次數(shù)(1000)的比值。所以每個(gè)N2都對(duì)應(yīng)一個(gè)物種滅絕概率,據(jù)此可得知物種滅絕概率達(dá)到最低時(shí)N2的值,這個(gè)值就是最優(yōu)保護(hù)區(qū)數(shù)量。

在每一次模擬中,從初始時(shí)刻到模擬終止時(shí)刻之間,被保護(hù)物種個(gè)體數(shù)量的變動(dòng)以及所對(duì)應(yīng)的時(shí)間點(diǎn)都被記錄。每個(gè)保護(hù)區(qū)受到擾動(dòng)影響的概率都為μ,所以擾動(dòng)是一個(gè)隨機(jī)發(fā)生的事件,也是模型所包含的隨機(jī)因素。

1.5 初始化

使用元胞自動(dòng)機(jī)模擬本文的模型,元胞自動(dòng)機(jī)的每個(gè)格子對(duì)應(yīng)于一個(gè)保護(hù)區(qū)或被破壞的生境。為了保證保護(hù)區(qū)數(shù)量從1變化到很多,設(shè)定元胞自動(dòng)機(jī)有N1×N1(N1=2, 3, …, 31)個(gè)格子,所以保護(hù)區(qū)的數(shù)量為N2=N1×N1×(1-D)(D為被破壞生境的數(shù)量與格子數(shù)量的比值。因?yàn)榧俣ǜ鞅Ｗo(hù)區(qū)和各被破壞生境的面積都相同,所以D也為所有被破壞生境的總面積與所選擇區(qū)域的面積的比值,即生境喪失比例),被破壞生境的數(shù)量為N1×N1×D。因?yàn)槊總€(gè)保護(hù)區(qū)在初始時(shí)刻的環(huán)境可持續(xù)承載量都相同,所以單個(gè)保護(hù)區(qū)在初始時(shí)刻的環(huán)境可持續(xù)承載量為KTOT/N2,KTOT為所有保護(hù)區(qū)在初始時(shí)刻的環(huán)境可持續(xù)承載量。因?yàn)樯嬗诿總€(gè)保護(hù)區(qū)上的局部種群的大小在初始時(shí)刻也相同,所以假定初始時(shí)刻局部種群的大小為所在保護(hù)區(qū)環(huán)境承載量的一半——KTOT/(N2×2)。被破壞生境的位置由D和qE/E(從一個(gè)被破壞生境的鄰居中隨機(jī)挑一個(gè),其也為被破壞生境的概率,也稱之為聚集度[33])決定,保護(hù)區(qū)占據(jù)剩余的位置。

1.6 子模型

1.6.1 自然保護(hù)區(qū)模型

保護(hù)區(qū)共有N2=N1×N1×(1-D)(N1=2, 3, …, 31)個(gè),被破壞的生境共有N1×N1×D個(gè)。被破壞生境的位置由D和qE/E決定,保護(hù)區(qū)占據(jù)剩余的位置。若被破壞生境的空間分布方式是隨機(jī)的,則在任一被破壞生境的鄰居中隨機(jī)選一個(gè),其也為被破壞生境的概率(qE/E)為D,即D=qE/E;若D>qE/E或D

因?yàn)樗斜Ｗo(hù)區(qū)的總面積為固定常數(shù),每個(gè)保護(hù)區(qū)或被破壞生境的面積都相同,所以若保護(hù)區(qū)的總面積為S,單個(gè)保護(hù)區(qū)或被破壞生境的面積都為S/N2,且單個(gè)保護(hù)區(qū)的面積與保護(hù)區(qū)數(shù)量N2成反比。圖2為兩種類型的保護(hù)區(qū)模式：少量大保護(hù)區(qū)模式和大量小保護(hù)區(qū)模式。

圖2 兩種保護(hù)區(qū)模式Fig.2 Two types of nature reserve modelsa—e: 少量大保護(hù)區(qū)模式;f—j: 大量小保護(hù)區(qū)模式;D=0.5; 如1.4部分所述,每一幅圖中保護(hù)區(qū)(灰色格子)的總面積都相同且為一固定常數(shù),每幅圖中所有保護(hù)區(qū)環(huán)境承載量的總和也相同。圖中灰色格子為保護(hù)區(qū),白色格子為被破壞的生境,黑線為各保護(hù)區(qū)或各被破壞生境的分界線

圖3 保護(hù)區(qū)A和其他保護(hù)區(qū)間的個(gè)體擴(kuò)散 Fig.3 Individual diffusion between reserve A and other reserves圖3只給出保護(hù)區(qū)A和其他保護(hù)區(qū)間的個(gè)體擴(kuò)散,其他保護(hù)區(qū)之間的個(gè)體擴(kuò)散沒有給出?；疑褡訛楸Ｗo(hù)區(qū),黑色格子為被破壞的生境,黑色長(zhǎng)線段為保護(hù)區(qū)或被破壞生境的邊緣,黑色短線段為保護(hù)區(qū)A和其他保護(hù)區(qū)間的個(gè)體擴(kuò)散。保護(hù)區(qū)A共有8個(gè)鄰居,其中5個(gè)(N3)為保護(hù)區(qū),其余3個(gè)為被破壞的生境

1.6.2 物種模型

T時(shí)刻,保護(hù)區(qū)I上的局部種群為NI(T),KI(T)為保護(hù)區(qū)I的環(huán)境可持續(xù)承載量。初始時(shí)刻,所有保護(hù)區(qū)的環(huán)境可持續(xù)承載量為KTOT,K為單個(gè)保護(hù)區(qū)的環(huán)境可持續(xù)承載量。因?yàn)樵诔跏紩r(shí)刻,每個(gè)保護(hù)區(qū)的環(huán)境可持續(xù)承載量都相同,所以K=KTOT/N2,N2為保護(hù)區(qū)數(shù)量。

種群增長(zhǎng)：T時(shí)刻,保護(hù)區(qū)I上的局部種群為[34]

NI(T)=R1×NI(T-1)/(1+(R1-1)×NI(T-1)/KI(T-1))

(1)

R(R1≥1)為無競(jìng)爭(zhēng)情況下,局部種群的增長(zhǎng)率;而R1/(1+(R1-1)×NI(T-1)/KI(T-1))-1為局部種群的實(shí)際增長(zhǎng)率。從方程(1)可得到：局部種群越接近環(huán)境承載量,實(shí)際增長(zhǎng)率越低;當(dāng)局部種群達(dá)到環(huán)境承載量時(shí),實(shí)際增長(zhǎng)率為0;若其超過環(huán)境承載量,種群會(huì)逐步減小。

環(huán)境可持續(xù)承載量增長(zhǎng)：T時(shí)刻，保護(hù)區(qū)I的環(huán)境承載量為[34]

KI(T)=R2×KI(T-1)/(1+(R2-1)×KI(T-1)/K)

(2)

R2(R2>1)為環(huán)境承載量的增長(zhǎng)率。R2/(1+(R2-1)×KI(T-1)/K)-1為環(huán)境承載量的實(shí)際增長(zhǎng)率,其變動(dòng)方式與局部種群實(shí)際增長(zhǎng)率的變動(dòng)方式類似,但KI(T)≤K。

擾動(dòng)：每個(gè)保護(hù)區(qū)受到擾動(dòng)影響的概率都為μ,擾動(dòng)導(dǎo)致環(huán)境承載量降低和局部種群減小。T時(shí)刻,如果保護(hù)區(qū)I受到擾動(dòng)影響,環(huán)境承載量和局部種群變化如下：

NI(T)=NI(T-1)×(1-d1)

(3)

KI(T)=KI(T-1)×(1-d1)

(4)

d1(0

擴(kuò)散：假設(shè)任一保護(hù)區(qū)中的個(gè)體只能擴(kuò)散到其周圍的保護(hù)區(qū)中,且不能穿過被破壞的生境。任一保護(hù)區(qū)中的個(gè)體只能擴(kuò)散到其周圍八個(gè)鄰居(摩爾鄰居關(guān)系),且其鄰居必須也是保護(hù)區(qū)(圖3)。MIJ(T) 為T時(shí)刻從保護(hù)區(qū)I擴(kuò)散到保護(hù)區(qū)J中的個(gè)體。

(5)

(6)

2 結(jié)果

2.1 生境喪失對(duì)FLOMS問題的影響

圖4展示隨機(jī)或非隨機(jī)的生境喪失對(duì)FLOMS問題的影響。表明：隨機(jī)的生境喪失(D=qE/E)下,生境喪失比例為0.2時(shí),最優(yōu)保護(hù)區(qū)數(shù)量超過170;若生境喪失比例上升,最優(yōu)保護(hù)區(qū)數(shù)量將大幅減少,當(dāng)生境喪失比例上升到0.9時(shí),最優(yōu)保護(hù)區(qū)數(shù)量已不足20。所以隨機(jī)的生境喪失下,生境喪失比例越高,少量大保護(hù)區(qū)模式的優(yōu)勢(shì)越明顯(保護(hù)區(qū)數(shù)量與各保護(hù)區(qū)的面積成反比)。

圖4b表明：若生境喪失比例恒定,當(dāng)聚集度為0.1時(shí),最優(yōu)保護(hù)區(qū)數(shù)量不足20;聚集度提高,最優(yōu)保護(hù)區(qū)數(shù)量將大幅增加,當(dāng)聚集度提高到0.9時(shí),最優(yōu)保護(hù)區(qū)數(shù)量已超過120。所以即使生境喪失比例保持不變,被破壞生境的空間分布方式也會(huì)有重要影響 - 被破壞生境的空間聚集程度越高,大量小保護(hù)區(qū)模式的優(yōu)勢(shì)越明顯。

圖4 生境喪失比例和被破壞生境的聚集程度對(duì)最優(yōu)保護(hù)區(qū)數(shù)量的影響Fig.4 The effect of proportion of lost habitats and the degree of clustering of lost habitats on the optimal reserve number參數(shù)：R1=R2=1.013, kTOT=2000, μ=0.01, ?=0.1, d1=0.5, d2=0, a=50, T=1100;聚集度: D=0.5

2.2 擴(kuò)散(擴(kuò)散率和擴(kuò)散死亡率)對(duì)FLOMS問題的影響

通過圖5,探究生境喪失比例(D)恒定,但被破壞生境的空間分布(用qE/E衡量)變化時(shí),擴(kuò)散率和擴(kuò)散死亡率對(duì)FLOMS問題的影響。擴(kuò)散率為0時(shí),最優(yōu)保護(hù)區(qū)數(shù)量非常少(不足20);當(dāng)擴(kuò)散率增加時(shí),最優(yōu)保護(hù)區(qū)數(shù)量也增加。被破壞生境的聚集程度較高時(shí)(聚集度不小于0.7),最優(yōu)保護(hù)區(qū)數(shù)量的增幅更加明顯;即使被破壞生境的聚集程度非常低時(shí)(聚集度為0.2),最優(yōu)保護(hù)區(qū)數(shù)量也會(huì)增加,但增幅較小。所以增加擴(kuò)散率可導(dǎo)致從少量大保護(hù)區(qū)更有利于物種到大量小保護(hù)區(qū)的轉(zhuǎn)變,且被破壞生境的聚集程度越高,轉(zhuǎn)變的程度越高。

圖5還表明：當(dāng)擴(kuò)散死亡率較高時(shí)(0.9),最優(yōu)保護(hù)區(qū)數(shù)量很少;若擴(kuò)散死亡率降低,最優(yōu)保護(hù)區(qū)數(shù)量增加,且被破壞生境的聚集程度較高時(shí)(聚集度不小于0.7),最優(yōu)保護(hù)區(qū)的增幅更加明顯。所以降低擴(kuò)散死亡可導(dǎo)致從少量大保護(hù)區(qū)更有利于物種到大量小保護(hù)區(qū)的轉(zhuǎn)變,且被破壞生境的聚集程度越高,轉(zhuǎn)變的程度越高。

圖5 個(gè)體擴(kuò)散率和個(gè)體的擴(kuò)散死亡率對(duì)最優(yōu)保護(hù)區(qū)數(shù)量的影響Fig.5 The optimal number of reserves changes with diffusion rate and diffusion mortality rate參數(shù)：擴(kuò)散率：D=0.5, R1=R2=1.012, kTOT=2000, μ=0.016, ?=0.2, d1=0.5, a=50, T=1200;擴(kuò)散死亡率：D=0.5, R1=R2=1.013, kTOT=2000, μ=0.01, d1=0.5, d2=0, a=50, T=1200

3 實(shí)例

作為“活化石”和“中國(guó)國(guó)寶”的大熊貓?jiān)?jīng)廣泛分布在華北、西北、華東、西南、華南以至越南和緬甸北部。但自近代以來,由于人口急劇增加,導(dǎo)致所需的土地越來越多,很多大熊貓?jiān)?jīng)的生境已轉(zhuǎn)變?yōu)槿祟惿畹某鞘泻袜l(xiāng)村,剩余的生境也因?yàn)槿祟悓?duì)于森林的過度砍伐而遭到嚴(yán)重破壞。因此,大熊貓的數(shù)量急劇減小到不足1000只,成為瀕危物種。為保護(hù)大熊貓,自臥龍大熊貓自然保護(hù)區(qū)創(chuàng)建以來,中國(guó)已陸續(xù)建立起67個(gè)大熊貓保護(hù)區(qū)。保護(hù)措施取得了明顯的成就,一系列大范圍國(guó)家調(diào)查顯示,此前的大熊貓數(shù)量下滑趨勢(shì)已被扭轉(zhuǎn)為穩(wěn)定上升,目前野外大熊貓的數(shù)量已上升到2000只左右。但隨著保護(hù)區(qū)周圍地區(qū)經(jīng)濟(jì)的發(fā)展,城市和鄉(xiāng)村不斷擴(kuò)大,也由于交通的發(fā)展,公路和鐵路密度越來越大,這些都造成了大熊貓生境的進(jìn)一步喪失和空間聚集程度進(jìn)一步降低。為了應(yīng)對(duì)這一變化,國(guó)家決定將所有67個(gè)大熊貓保護(hù)區(qū)整合為由岷山、邛崍山——大相嶺、秦嶺和白水江這4個(gè)片區(qū)組成的大熊貓公園,以加強(qiáng)以大熊貓為核心的生物多樣性保護(hù)。大熊貓生境生境喪失程度的增加和空間聚集程度的降低,導(dǎo)致現(xiàn)今近70個(gè)保護(hù)區(qū)將被整合為4個(gè),這成為證明“生境喪失比例增加或被破壞生境的空間聚集程度降低時(shí),少量大保護(hù)區(qū)更優(yōu)”的現(xiàn)實(shí)例子。

4 討論

生境喪失不僅造成很多珍稀的動(dòng)植物瀕臨滅絕,還造成其所在的生境破碎化[35- 36]。雖然可以借助生態(tài)技術(shù)幫助動(dòng)植物恢復(fù)部分被破壞的生境[37- 38],但要恢復(fù)所有生境,并任意選擇所需的生境建立自然保護(hù)區(qū),已沒有可能。因而在建立保護(hù)區(qū)時(shí),必須處理生境喪失所帶來的問題。本文選擇其中的部分問題進(jìn)行探討,在一個(gè)已經(jīng)遭受過生境喪失的景觀上,選取一片區(qū)域,并保證此區(qū)域中未被破壞生境的面積為一個(gè)固定常數(shù),然后探討應(yīng)該將未被破壞生境設(shè)成少量大保護(hù)區(qū)還是大量小保護(hù)區(qū)。結(jié)果表明：(1)生境喪失比例和被破壞生境的空間分布形式共同影響FLOMS問題——若被破壞生境的空間分布形式保持不變,生境喪失比例越高,少量大保護(hù)區(qū)模式的優(yōu)勢(shì)越明顯;若生境喪失比例恒定,被破壞生境的空間聚集程度越高,大量小保護(hù)區(qū)模式的優(yōu)勢(shì)越明顯。(2)生境喪失還對(duì)FLOMS問題存在間接影響,擴(kuò)散率增加或擴(kuò)散死亡率降低可導(dǎo)致從少量大保護(hù)區(qū)更有利于物種到大量小保護(hù)區(qū)最優(yōu)的轉(zhuǎn)變,且被破壞生境的聚集程度越高,轉(zhuǎn)變程度越劇烈。

之前的研究表明生境喪失會(huì)通過對(duì)其他因子(種群增長(zhǎng)率、擾動(dòng)或環(huán)境承載量)產(chǎn)生作用來影響FLOMS問題[20,22]。例如：當(dāng)被破壞的生境隨機(jī)分布在景觀上時(shí),Reed[23]的研究表明,種群增長(zhǎng)率越高,越適宜采用大量小保護(hù)區(qū)模式。本文的部分研究也考慮生境喪失通過對(duì)其他因子(個(gè)體擴(kuò)散)產(chǎn)生作用來影響FLOMS問題,但之前的研究均假設(shè)被破壞的生境隨機(jī)分布在景觀上,而相關(guān)研究表明被破壞生境的分布方式是非隨機(jī)的[29-30],因而這些假設(shè)與實(shí)際不符。此外,這些研究也僅涉及到生境喪失對(duì)FLOMS問題的間接作用,而沒有探討直接作用。本文解決了以上問題,結(jié)果表明：不但生境喪失比例直接對(duì)FLOMS問題造成影響,而且即使生境喪失比例保持不變,被破壞生境的空間分布特征對(duì)FLOMS問題也造成影響。所以,在探討生境喪失時(shí),必須同時(shí)考慮生境喪失比例和被破壞生境的空間分布特征(即：生境喪失的空間特征),尤其是后者,因?yàn)橐郧暗难芯客匾暻罢?而忽視后者[20,22- 23]。本文的研究結(jié)果還表明：被破壞生境的聚集程度還會(huì)影響個(gè)體擴(kuò)散(擴(kuò)散率和擴(kuò)散死亡率)對(duì)FLOMS問題的作用——被破壞生境的聚集度較高時(shí),保護(hù)區(qū)最優(yōu)數(shù)量對(duì)個(gè)體擴(kuò)散的反應(yīng)更劇烈。

本文基于元胞自動(dòng)機(jī)模型,研究生境喪失和個(gè)體擴(kuò)散的共同作用對(duì)自然保護(hù)區(qū)設(shè)計(jì)中關(guān)鍵問題——FLOMS問題——的影響。研究結(jié)果可為生物多樣性保護(hù)和保護(hù)區(qū)設(shè)計(jì)提供理論基礎(chǔ),有助于人們加深對(duì)物種如何響應(yīng)生境喪失的了解,增強(qiáng)人們對(duì)生境破壞的重視程度,提高人們保護(hù)自然環(huán)境的意識(shí)。但個(gè)體不僅能短距離擴(kuò)散到鄰近地區(qū),還可以借助廊道,長(zhǎng)距離擴(kuò)散到相距較遠(yuǎn)的區(qū)域;生境喪失可能與短距離擴(kuò)散和長(zhǎng)距離擴(kuò)散共同作用于種群動(dòng)態(tài)[39],從而影響影響FLOMS問題。但這些本文都沒涉及,此外本文還存在需要改進(jìn)的問題,雖然很多自然保護(hù)區(qū)都有其主要保護(hù)的物種,如中國(guó)的大熊貓自然保護(hù)區(qū),但自然保護(hù)區(qū)內(nèi)一般會(huì)有多個(gè)需要保護(hù)的物種,因此在設(shè)計(jì)自然保護(hù)區(qū)時(shí),除了考慮其主要保護(hù)的對(duì)象外,還必須兼顧所有需要保護(hù)的物種,對(duì)生境喪失如何影響FLOMS問題的研究還需深入討論。

參考文獻(xiàn)(References):

[1] 陳成忠, 葛緒廣, 孫琳, 邵冬, 柯文莉. 物種急劇喪失·生態(tài)嚴(yán)重超載·跨越“地球邊界”·區(qū)域公平失衡·“一個(gè)地球”生活——《地球生命力報(bào)告2014》解讀. 生態(tài)學(xué)報(bào), 2016, 36(9): 2779- 2785.

[2] Soule M E. Conservation: tactics for a constant crisis. Science, 1991, 253(5021): 744- 750.

[3] Van Cuong C, Dart P, Dudley N, Hockings M. Factors influencing successful implementation of biosphere reserves in vietnam: challenges, opportunities and lessons learnt. Environmental Science & Policy, 2017, 67: 16- 26.

[4] Gong M H, Fan Z Y, Wang J Y, Liu G, Lin C. Delineating the ecological conservation redline based on the persistence of key species: giant pandas (Ailuropodamelanoleuca) inhabiting the Qinling mountains. Ecological Modelling, 2017, 345: 56- 62.

[5] Diamond J M. The island dilemma: lessons of modern biogeographic studies for the design of natural reserves. Biological Conservation, 1975, 7(2): 129- 146.

[6] Wilson E O, Willis E O. Applied biogeography//Cody M L, Diamond J M, eds. Ecology and Evolution of Communities. Cambridge, Massachusetts: Belknap Press, 1975: 522- 534.

[7] May R M. Islands biogeography and the design of wildlife preserves. Nature, 1975, 254(5497): 177- 178.

[8] Gilpin M E, Diamond J M. Subdivision of nature reserves and the maintenance of species diversity. Nature, 1980, 285(5766): 567- 568.

[9] Groeneveld R. Economic considerations in the optimal size and number of reserve sites. Ecological Economics, 2005, 52(2): 219- 228.

[10] Baker T P, Jordan G J, Dalton P J, Baker S C. Impact of distance to mature forest on the recolonisation of bryophytes in a regenerating Tasmanian wet eucalypt forest. Australian Journal of Botany, 2013, 61(8): 633- 642.

[11] Lindenmayer D B, Wood J, McBurney L, Blair D, Banks S C. Single large versus several small: the SLOSS debate in the context of bird responses to a variable retention logging experiment. Forest Ecology and Management, 2015, 339: 1- 10.

[12] Etienne R S, Heesterbeek J A P. On optimal size and number of reserves for metapopulation persistence. Journal of Theoretical Biology, 2000, 203(1): 33- 50.

[13] Burkey T V. Metapopulation extinction in fragmented landscapes: using bacteria and protozoa communities as model ecosystems. The American Naturalist, 1997, 150(5): 568- 591.

[14] Tj?rve E. How to resolve the SLOSS debate: lessons from species-diversity models. Journal of Theoretical Biology, 2010, 264(2): 604- 612.

[15] McCarthy M A, Thompson C J, Moore A L, Possingham H P. Designing nature reserves in the face of uncertainty. Ecology Letters, 2011, 14(5): 470- 475.

[16] Boecklen W J. Nestedness, biogeographic theory, and the design of nature reserves. Oecologia, 1997, 112(1): 123- 142.

[17] Skaggs R W, Boecklen W J. Extinctions of montane mammals reconsidered: putting a global-warming scenario on ice. Biodiversity & Conservation, 1996, 5(6): 759- 778.

[18] Le Roux D S, Ikin K, Lindenmayer D B, Manning A D, Gibbons P. Single large or several small? Applying biogeographic principles to tree-level conservation and biodiversity offsets. Biological Conservation, 2015, 191: 558- 566.

[19] Ovaskainen O. Long-term persistence of species and the SLOSS problem. Journal of Theoretical Biology, 2002, 218(4): 419- 433.

[20] Robert A. The effects of spatially correlated perturbations and habitat configuration on metapopulation persistence. Oikos, 2009, 118(10): 1590- 1600.

[21] Saunders D A, Hobbs R J, Margules C R. Biological consequences of ecosystem fragmentation: a review. Conservation Biology, 1991, 5(1): 18- 32.

[22] Pelletier J D. Model assessments of the optimal design of nature reserves for maximizing species longevity. Journal of Theoretical Biology, 2000, 202(1): 25- 32.

[23] Reed D H. Extinction risk in fragmented habitats. Animal Conservation, 2004, 7(2): 181- 191.

[24] Richards C M. Inbreeding depression and genetic rescue in a plant metapopulation. The American Naturalist, 2000, 155(3): 383- 394.

[25] Harding K C, McNamara J M. A unifying framework for metapopulation dynamics. The American Naturalist, 2002, 160(2): 173- 185.

[26] Gonzalez A, Lawton J H, Gilbert F S, Blackburn T M, Evans-Freke I. Metapopulation dynamics, abundance, and distribution in a microecosystem. Science, 1998, 281(5385): 2045- 2047.

[27] Debinski D M, Holt R D. A survey and overview of habitat fragmentation experiments. Conservation Biology, 2000, 14(2): 342- 355.

[28] Quinn J F, Hastings A. Extinction in subdivided habitats. Conservation Biology, 1987, 1(3): 198- 208.

[29] Boakes E H, Mace G M, McGowan P J K, Fuller R A. Extreme contagion in global habitat clearance. Proceedings of the Royal Society B: Biological Sciences, 2010, 277(1684): 1081- 1085.

[30] Xie Y J, Ng C N. Exploring spatio-temporal variations of habitat loss and its causal factors in the Shenzhen river cross-border watershed. Applied Geography, 2013, 39: 140- 150.

[31] Grimm V, Berger U, Bastiansen F, Eliassen S, Ginot V, Giske J, Goss-Custard J, Grand T, Heinz S K, Huse G, Huth A, Jepsen J U, J?rgensen C, Mooij W M, Müller B, Pe′er G, Piou C, Railsback S F, Robbins A M, Robbins M M, Rossmanith E, Rüger N, Strand E, Souissi S, Stillman R A, Vab? R, Visser U, DeAngelis D L. A standard protocol for describing individual-based and agent-based models. Ecological Modelling, 2006, 198(1/2): 115- 126.

[32] Grimm V, Berger U, DeAngelis D L, Polhill J G, Giskee J, Railsback S F. The ODD protocol: a review and first update. Ecological Modelling, 2010, 221(23): 2760- 2768.

[33] Hiebeler D. Populations on fragmented landscapes with spatially structured heterogeneities: landscape generation and local dispersal. Ecology, 2000, 81(6): 1629- 1641.

[34] Johst K, Drechsler M. Are spatially correlated or uncorrelated disturbance regimes better for the survival of species? Oikos, 2003, 103(3): 449- 456.

[35] De Montis A, Martín B, Ortega E, Ledda A, Serra V. Landscape fragmentation in Mediterranean Europe: a comparative approach. Land Use Policy, 2017, 64: 83- 94.

[36] 陳利頂, 劉雪華, 傅伯杰. 臥龍自然保護(hù)區(qū)大熊貓生境破碎化研究. 生態(tài)學(xué)報(bào), 1999, 19(3): 291- 297.

[37] 劉春霞, 韓烈保. 高速公路邊坡植被恢復(fù)研究進(jìn)展. 生態(tài)學(xué)報(bào), 2007, 27(5): 2090- 2098.

[38] 劉二佳, 張曉萍, 謝名禮, 陳妮, 張亭亭, 郭敏杰, 張建軍. 生態(tài)恢復(fù)對(duì)流域水沙演變趨勢(shì)的影響——以北洛河上游為例. 生態(tài)學(xué)報(bào), 2015,35(3): 622- 629.

[39] Hiebeler D. Competition between near and far dispersers in spatially structured habitats. Theoretical Population Biology, 2004, 66(3): 205- 218.