南滾河流域土地利用變化對亞洲象生境的影響

吳得卿,魏建華,樊 輝

1 云南大學(xué)國際河流與生態(tài)安全研究院,昆明 650091 2 六盤水師范學(xué)院旅游與歷史文化學(xué)院,六盤水 553001 3 中國市政工程華北設(shè)計(jì)研究總院有限公司 北京分公司,北京 100044

亞洲象(Elephasmaximus)是世界重要瀕危物種之一,為國家Ⅰ級保護(hù)野生動物,主要分布于云南西雙版納州、普洱市和臨滄市等地[1- 4],現(xiàn)存種群數(shù)量300余頭[4- 6]。亞洲象數(shù)量減少的主要原因是其生境的喪失和退化[1]。因土地利用活動加劇,尤其是近30年來全球橡膠需求驅(qū)動下橡膠林?jǐn)U張導(dǎo)致天然林面積減少[7- 8],嚴(yán)重威脅著區(qū)域生態(tài)環(huán)境及生物多樣性[9]。亞洲象生境與天然林分布高度疊合[2, 10],天然林減少導(dǎo)致的生境喪失使其種群生存難以為繼,人象沖突事件逐年攀升,保護(hù)亞洲象、緩解人象沖突已成為當(dāng)務(wù)之急。現(xiàn)有亞洲象保護(hù)研究主要致力于種群數(shù)量調(diào)查[5, 11-12]、生境評價(jià)與保護(hù)[1- 2, 13- 15]以及人象沖突等[16- 18]。但較少從土地利用變化的角度揭示亞洲象生境的長期變化過程。

基于Landsat時(shí)序影像的遙感變化檢測能診斷土地利用時(shí)序變化過程[19],具有長時(shí)序及較高空間分辨率優(yōu)勢[20-21],已逐漸成為一種主流的檢測方法[22]。生態(tài)位因子分析模型可根據(jù)亞洲象出現(xiàn)點(diǎn)的生態(tài)位特征預(yù)測亞洲象適宜生境范圍[23]。亞洲象近乎只生存在完整的天然林中[24],天然林的時(shí)序變化動態(tài)能較好地反映亞洲象生境變化的時(shí)空過程[25]。以中-緬邊境的南滾河跨境流域?yàn)槔?運(yùn)用基于LandTrendr算法的時(shí)序變化檢測法和生態(tài)位因子分析模型,揭示1988—2018年土地利用變化對亞洲象生境影響的時(shí)空過程,以期服務(wù)于跨境區(qū)域亞洲象生境保護(hù)、種群規(guī)模維持以及人象沖突緩解等。

1 研究區(qū)概況

南滾河流域位于東經(jīng)98°34′—99°12′,北緯22°58′—23°22′(圖1),為中、緬跨境流域,海拔413—2632 m,地形復(fù)雜多樣；流域總面積1666.67 km2,中、緬境內(nèi)面積分別為554.03 km2(占33.24%)和1112.64 km2(占66.76%)。該流域?qū)俚湫臀髂霞撅L(fēng)氣候區(qū),干濕季分明,雨季為5—10月,旱季為11月至次年4月,年均降水量1700—2500 mm,年均氣溫16—25 ℃；動植物熱帶性質(zhì)較為明顯[1, 26],生物多樣性豐富。流域內(nèi)主要土地利用類型包括六類：天然林、橡膠林、園地(茶園及核桃套種等)、耕地、建設(shè)用地以及水體(河流和水塘)。

南滾河流域中國境內(nèi)亞洲象現(xiàn)僅存12頭[27],主要棲息于海拔1000 m以下且坡度小于15°的河谷地區(qū)[1, 28]。20世紀(jì)80年代以前,亞洲象常遷徙于中、緬兩國之間覓食,但隨著土地利用變化加劇,中-緬邊境森林被破壞,導(dǎo)致亞洲象棲息地破碎化[26, 28],近30年已無跨境種群交流[27],形成極小種群,生存和繁衍受到嚴(yán)重威脅[29]。

圖1 研究區(qū)位置及地形、野外調(diào)查樣線和亞洲象活動范圍Fig.1 The location and topography, and transect of field survey and Asian elephants′ home range國界線和省界線來源于全國地理信息資源目錄服務(wù)系統(tǒng) https://www.webmap.cn/main.do?method=index

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù)來源與處理

使用的數(shù)據(jù)包括：南滾河流域1988—2018年間Landsat逐年時(shí)序數(shù)據(jù)、地形數(shù)據(jù)、亞洲象出現(xiàn)點(diǎn)數(shù)據(jù)、1988年和2018年兩期土地利用/土地覆被遙感分類圖等。Landsat數(shù)據(jù)來源于https://www.usgs.gov/core-science-systems/nli/landsat,采用1988—2018年間云量小于10%的地表反射率TM、ETM+和OLI影像,并預(yù)處理,構(gòu)建每年旱季高質(zhì)量像元影像時(shí)序棧,其中無云像元比例達(dá)95%以上；地形數(shù)據(jù)來源于https://earthexplorer.usgs.gov/的航天雷達(dá)地形測繪任務(wù)的數(shù)字高程模型(SRTM DEM),空間分辨率為30 m；亞洲象出現(xiàn)點(diǎn)數(shù)據(jù)主要包括2018年的野外調(diào)查樣線和馮利民等[26]的樣線調(diào)查數(shù)據(jù)及亞洲象活動范圍圖,并參考郎學(xué)東等[30]提供的亞洲象活動痕跡樣點(diǎn)、Bohnett等[31]提供的紅外相機(jī)監(jiān)測數(shù)據(jù)、陳明勇團(tuán)隊(duì)[32]提供的無人機(jī)監(jiān)測數(shù)據(jù)(圖1)；1988年和2018年的土地利用土地/覆被遙感分類大致流程如下：根據(jù)南滾河流域中國境內(nèi)2018年土地利用現(xiàn)狀調(diào)查,建立土地利用分類系統(tǒng),結(jié)合Google EarthTM高空間分辨率影像繪制2018年南滾河流域土地利用分類地面參照樣本,并通過雙時(shí)相影像差異變化檢測,得到1988年南滾河流域土地利用分類地面參照樣本,基于1988年2月Landsat 5 TM 1—5、7和2018年2月Landsat 8 OLI 2—7波段影像,采用隨機(jī)森林分類器整合光譜、地形和紋理信息進(jìn)行監(jiān)督分類,分別得到1988和2018年該流域包括天然林、橡膠林、耕地、建設(shè)用地、園地和水體六類的土地利用/土地覆被遙感分類圖[33]。

2.2 研究方法

2.2.1基于ENFA模型的生境評價(jià)

ENFA模型根據(jù)收集的亞洲象出現(xiàn)點(diǎn)及其所處的生態(tài)環(huán)境,對比其生態(tài)位與環(huán)境因子間的異同,計(jì)算2018年亞洲象的適宜生境范圍[23],其建模與分析在Biomapper 4.0(https://www2.unil.ch/biomapper/)中完成?；谀蠞L河流域2018年的土地利用/土地覆被圖,獲取生境評價(jià)所需的植被因子和干擾要素,根據(jù)該流域自然條件[26, 28],結(jié)合已有研究篩選出的對亞洲象適宜生境貢獻(xiàn)率較高的因子[1, 14, 25],選取了12個(gè)生態(tài)地理變量(EGVs)作為模型輸入(表1),并計(jì)算其邊際值M、特異值S和耐受值T以描述亞洲象生態(tài)位特征[23, 33]。其中,M值為EGV在整個(gè)南滾河流域與該流域內(nèi)亞洲象分布區(qū)的差異,值越大,則亞洲象對EGV的選擇性越強(qiáng)；S為流域內(nèi)亞洲象對EGV的選擇范圍,值越大,其生態(tài)位的特化程度越高；T為亞洲象對EGV的耐受性,值越小,則亞洲象對EGV的耐受性越低、生存范圍越狹窄[23]。與主成分分析類似,ENFA模型將所選取的12個(gè)EGVs變換為包含主要信息的邊際因子和特異因子,并選取累計(jì)貢獻(xiàn)率超過90%的因子計(jì)算生境適宜性指數(shù)(HSI)。ENFA模型采用Boyce Index的交叉驗(yàn)證方法評價(jià)其準(zhǔn)確度,并根據(jù)Boyce Index曲線間斷點(diǎn)將HSI劃分為四個(gè)等級：最適生境、較適生境、邊際生境和非生境[14]。

由于1988年亞洲象出現(xiàn)點(diǎn)數(shù)據(jù)難以獲取,因此據(jù)ENFA模型評價(jià)得到的2018年各類EGVs對亞洲象適宜生境的貢獻(xiàn)率,假設(shè)亞洲象30年間對環(huán)境偏好不變,采用賦值累乘法,結(jié)合1988年土地利用/土地覆被圖為EGV賦予相應(yīng)權(quán)重,計(jì)算每個(gè)評價(jià)單元的亞洲象適宜性指數(shù)[15],得到1988年的亞洲象適宜生境及其空間分布范圍。

表1 生態(tài)地理變量(EGVs)及其處理方式

2.2.2森林?jǐn)_動時(shí)序變化檢測

亞洲象幾乎只生存在完整的天然林中[24],天然林狀態(tài)的變化則意味著其生境的適宜性發(fā)生改變。ARVI能準(zhǔn)確區(qū)分橡膠林與其它地類(含天然林)[34],減小大氣及地形對時(shí)序變化檢測的影響,公式為：

不過，風(fēng)影也有快樂的地方，終究不必再做苦行僧了，不用在青燈黃卷中念經(jīng)敲木魚了，不用在院落中一下一下地清掃那些永遠(yuǎn)也掃不完的落葉了。再說，能與依然青春逼人的紅琴朝夕相處，終究也不失為一樁人間樂事。人間四月天是美好的，盡管風(fēng)影喜歡秋天比喜歡春天要多一點(diǎn)點(diǎn)，秋天的風(fēng)帶著絲絲涼意，是從往年吹過來的，那種帶點(diǎn)淡淡的傷感與世事的滄桑的味道，似乎更適合他這個(gè)曾經(jīng)皈依佛門而現(xiàn)在已經(jīng)還了俗的人，更吻合他那種天涼好個(gè)秋的心境。

(1)

式中,ρBlue、ρRed、ρNIR為Landsat影像中藍(lán)、紅、近紅外波段的地表反射率?；?988年土地利用/土地覆被圖,將天然林設(shè)為時(shí)序變化檢測的森林初始分布,用LandTrendr算法檢測森林?jǐn)_動時(shí)序變化,擾動閾值取ΔARVI= 0.32[7]。若某一像元的ARVI時(shí)序系列中其值年間變化大于0.32且保持低值不連續(xù)大于3年,則將其標(biāo)記為擾動。

2.2.3精度評價(jià)

采用Kappa系數(shù)和基于混淆矩陣的傳統(tǒng)精度評價(jià)方法計(jì)算土地利用分類精度與森林?jǐn)_動時(shí)序變化檢測精度。Kappa系數(shù)能客觀評價(jià)檢測結(jié)果的準(zhǔn)確性[35]；基于混淆矩陣的傳統(tǒng)精度評價(jià)方法的指標(biāo)包括總體精度(OA)、制圖精度(PA)、用戶精度(UA)[36]。其中,土地利用分類精度評價(jià)中所采用的參考地類為根據(jù)實(shí)地調(diào)查結(jié)合Google EarthTM高空間分辨率影像繪制所得,通過分層隨機(jī)抽樣法從1988年和2018年地面參照樣本中分別抽取兩個(gè)年份的驗(yàn)證樣本；森林?jǐn)_動時(shí)序變化檢測的驗(yàn)證樣本采用TimeSync方法繪制[37],共獲取320余個(gè)擾動樣本,結(jié)合未擾動森林樣本可對森林?jǐn)_動時(shí)序變化檢測結(jié)果進(jìn)行逐年精度評價(jià)。

2.2.4生境景觀格局分析

采用斑塊數(shù)(NP)、斑塊密度(PD)、最大斑塊指數(shù)(MPI)、平均斑塊面積(AREA_MN)、斑塊分維數(shù)(FRAC)、斑塊連通度指數(shù)(CONNECT)等6個(gè)常用景觀指數(shù)來分析亞洲象生境景觀格局特征[38],這些指標(biāo)由FRAGSTATS 3.3(https://www.umass.edu/landeco/research/fragstats/)計(jì)算得出。其中,NP、MPI、PD、AREA_MN值越小,則表示亞洲象生境退化越嚴(yán)重；FRAC取值范圍為[1,2],值越接近1,表示人類活動對生境適宜性影響越大,反之則??；CONNECT取值范圍[0,1],0表示生境間不存在生態(tài)連接、生境破碎化程度高,1代表生境間處于完全連接狀態(tài)。對比1988年到2018年亞洲象的生境景觀格局指數(shù),可得到30年間亞洲象的生境面積、受人類活動影響程度以及生境連通性等的變化情況。

3 研究結(jié)果

3.1 南滾河流域土地利用變化特征

南滾河流域1988年和2018年土地利用土地/覆被分類圖的OA、UA和PA均在90%以上(除1988年耕地的UA為87.27%)。從圖2可知,1988年到2018年,南滾河流域天然林面積急劇減少,主要向橡膠林轉(zhuǎn)化,其次是轉(zhuǎn)化為耕地、建設(shè)用地和園地。30年間,天然林從1495.45 km2減少到1291.91 km2,年均減少6.73 km2；橡膠林面積從48.07 km2擴(kuò)張到194.09 km2,年均增長4.87 km2。

圖2 南滾河流域1988年和2018年土地利用土地/覆被分類圖及其土地利用土地/覆被類型轉(zhuǎn)移圖Fig.2 Land use/land cover classification map of 1988 and 2018, and land use/land cover conversion from 1988 to 2018 in Nangun River Basin

3.2 南滾河流域亞洲象生境格局

ENFA模型中,前4個(gè)因子共解釋了91%的信息,可用于計(jì)算亞洲象生境適宜度,其中邊際因子解釋了100%的邊際性信息和18.97%的特異性信息[33](表2)。表2表明亞洲象生境主要分布于海拔低、坡度小、水源與食物充足、遠(yuǎn)離經(jīng)濟(jì)林和居民點(diǎn)的河谷地區(qū)。亞洲象生態(tài)位邊際值M=1.691>1,特異值S=2.865,耐受值T=0.349,表明其對生境的選擇非隨機(jī)、生境特化程度高,對環(huán)境有一定耐受性,但生境范圍較窄,受環(huán)境制約較大。據(jù)交叉檢驗(yàn)所得Boyce Index=0.825±0.2021,表明ENFA對亞洲象生境預(yù)測能力較強(qiáng)。根據(jù)HSI的自然裂點(diǎn)分級,亞洲象生境可分為最適生境(100≥HSI>62)、較適生境(62≥HSI>37)、邊際生境(37≥HSI>12)和非生境(12 ≥HSI>0)4個(gè)等級(圖3)。

表2 生態(tài)地理變量得分矩陣

圖3表明,1988年和2018年亞洲象生境主要分布于河流兩岸的河谷區(qū),最適和較適生境面積較少,邊際生境面積較大；最適生境相對集中,但其斑塊面積較小且斑塊間連通性低,生態(tài)位較窄。相較于1988年,2018年亞洲象生境范圍變窄,上、下游流域段最適生境被完全分隔；較適生境和邊際生境面積均大幅減少,分別減少32.38 km2和44.12 km2。

圖3 南滾河流域1988年和2018年亞洲象生境分布Fig.3 The habitat distribution of Asian elephants in Nangun River Basin in 1988 and 2018

圖4 1988年和2018年亞洲象生境景觀指數(shù)Fig.4 Landscape indices of Asian elephants′ habitat in 1988 and 2018

3.3 南滾河流域亞洲象生境時(shí)空變化過程

3.3.1森林?jǐn)_動變化檢測及精度

南滾河流域森林?jǐn)_動由天然林邊緣逐漸向森林內(nèi)部延伸(圖5)；1989—2018年間總擾動面積230.48 km2,占1988年森林總面積的15.61%,各年受擾動面積各異。據(jù)混淆矩陣評價(jià),森林?jǐn)_動檢測的OA為91.25%,未擾動區(qū)域的PA與UA分別為91.71%和92.40%,而擾動區(qū)域的PA與UA為92.42%和91.74%, Kappa系數(shù)0.88。

圖5 南滾河流域森林?jǐn)_動及亞洲象生境變化Fig.5 The map of forest disturbance and Asian elephants′ habitat change in Nangun River Basin

3.3.2森林?jǐn)_動時(shí)序變化特征

1988—2018年間,南滾河流域以橡膠林?jǐn)U張為主的森林?jǐn)_動明顯(圖6)。1989—1993年間擾動面積波動較大,1989年和1992年分別為24.69 km2和26.35 km2,遠(yuǎn)高于其他年份；1993—2003年擾動面積總體較小；2004—2008年擾動面積快速遞增,2008年達(dá)25.71 km2；此后(2012年除外)擾動面積回落至較低水平。但1989—2018年間該流域的中國境內(nèi)、外區(qū)域森林?jǐn)_動面積存在較大差異(圖6),其總面積分別為58.48 km2和171.89 km2,占1988年初期森林總面積的3.96%和11.64%；從時(shí)序上看,境外區(qū)域(緬甸境內(nèi))擾動面積大于中國境內(nèi)區(qū)域,以2001—2013年間最突出,其擾動面積在2004—2008年快速遞增,且在2009—2013年仍維持較大面積,而在中國境內(nèi)區(qū)域擾動面積在2009—2018年出現(xiàn)大幅減小。

圖6 1989—2018年南滾河流域森林?jǐn)_動面積Fig.6 The areas of disturbed forest in Nangun River Basin from 1989 to 2018

3.3.3融合森林?jǐn)_動的南滾河流域亞洲象生境變化時(shí)空動態(tài)

1988—2018年間,亞洲象適宜生境范圍內(nèi)森林?jǐn)_動總面積32.2 km2,占1988年其適宜生境的13.84%。1989年、1992年、2007年、2008年和2012年擾動面積較大,主要集中于流域中段(圖5)。30年間亞洲象生境面積逐漸縮小,破碎化加重,生境間連通性降低,各類生境均由高等級向次級、更次級乃至非生境轉(zhuǎn)出。最適生境減少了35.76%,較適生境約50%轉(zhuǎn)化為非生境；邊際生境主要轉(zhuǎn)化為非生境；適宜生境總面積從232.61 km2減至141.42 km2,減少量為91.19 km2,占1988年適宜生境總面積的39.20%(圖7)。

圖7 1988—2018年南滾河流域生境面積轉(zhuǎn)化和生境內(nèi)各時(shí)期森林?jǐn)_動面積Fig.7 The map of transferred habitat area from 1988 to 2018 and the area and area percentage of disturbed habitat in five-year intervals in Nangun River Basin

從6個(gè)不同時(shí)期來看(圖7),南滾河流域亞洲象適宜生境森林?jǐn)_動面積1989—1993年間較大(尤其在上游流域區(qū)),而1994—1998年、1999—2003年兩個(gè)時(shí)段相對較小,且均主要發(fā)生在較適及邊際生境內(nèi),但2001年以后最適生境受森林?jǐn)_動影響較大；2004—2008年間擾動面積最大,達(dá)15.97 km2,擾動比為49.59%,且主要集中在流域中段的適宜生境區(qū)；2009—2013年間擾動面積仍較大,主要發(fā)生在中下游流域區(qū)；但2014年以后森林?jǐn)_動面積極少。

綜合圖5和圖7可知,因森林?jǐn)_動導(dǎo)致亞洲象生境退化表現(xiàn)為：2001年以前(1989年和1992年除外)擾動面積較小,外圍生境退化；2001—2013年間擾動面積增大,最適生境受擾明顯,適宜生境大面積減少,特別是邊境地區(qū)適宜生境受擾后被分隔成上、下游兩個(gè)孤立區(qū)域；2013年以后適宜生境變動較少,趨于穩(wěn)定。

4 討論與結(jié)論

4.1 討論

(1)土地利用變化與森林?jǐn)_動的跨境差異：1988—2018年間,南滾河流域橡膠林種植擴(kuò)張導(dǎo)致毀林,這與以前眾多學(xué)者的研究結(jié)果一致[7- 8, 24, 39- 41]。因亞洲象幾乎只生存在完整的天然林中[1, 24],天然林的時(shí)空變化在較大程度上反映其生境時(shí)空演變特征。南滾河流域天然林的時(shí)序變化檢測結(jié)果顯示,近30年間以橡膠林?jǐn)U張導(dǎo)致毀林為主的森林?jǐn)_動面積中、緬跨境差異明顯。2001年以前中、緬境內(nèi)森林?jǐn)_動面積大致相當(dāng),波動性均較大,主要發(fā)生在邊際生境和較適生境區(qū)域；2001—2013年中-緬邊境天然林面積銳減,導(dǎo)致南滾河中游流域區(qū)亞洲象適宜生境喪失,被分割為上、下游兩個(gè)孤立區(qū)域,從地理上阻斷了中、緬亞洲象種群交流的機(jī)會[25, 27]；該時(shí)期緬甸境內(nèi)森林?jǐn)_動面積明顯大于中國境內(nèi),這與兩國發(fā)展水平差異以及相關(guān)社會經(jīng)濟(jì)政策實(shí)施有關(guān)。20世紀(jì)90年代后期,緬甸開展罌粟替代種植,從而促進(jìn)了橡膠林種植的迅速擴(kuò)張；而中國實(shí)施了“集體林權(quán)制度改革”和“天然林保護(hù)計(jì)劃”等措施,抑制了橡膠種植擴(kuò)張導(dǎo)致的毀林,表明政策是影響該流域土地利用變化的重要因素；2013年之后中、緬境內(nèi)森林?jǐn)_動面積均減少,亞洲象生境趨于穩(wěn)定,這與該流域內(nèi)橡膠適宜種植區(qū)趨于飽和,以及全球橡膠價(jià)格長期低迷有關(guān)[42]。

(2)土地利用變化對亞洲象生境的影響：流域內(nèi)土地利用變化導(dǎo)致亞洲象生境面積嚴(yán)重減少,破碎化程度增加。受天然林減少的持續(xù)影響,2018年南滾河流域亞洲象最適生境總面積僅剩26.34 km2,遠(yuǎn)低于其生存最小家域面積(105 km2)[43],也低于其在旱季的最小活動范圍需求35.67 km2[10]。大面積森林?jǐn)_動導(dǎo)致的適宜生境破碎化,阻斷了亞洲象往返于上、下游的遷徙通道,致使中、緬亞洲象種群交流中斷了近30年[25, 27]。橡膠林等熱作種植的擴(kuò)張擠占了亞洲象的生存空間,削減了天然食物來源,迫使亞洲象進(jìn)入人類生產(chǎn)生活區(qū)覓食,誘發(fā)人象沖突。最新調(diào)查顯示中國境內(nèi)的南滾河流域區(qū)目前僅存12頭亞洲象[27],適宜生境喪失與破碎化是制約亞洲象種群增長的重要原因[6, 25],而極小種群效應(yīng)是目前該流域亞洲象保護(hù)所面臨的最大挑戰(zhàn)[29]。此外,研究也表明南滾河國家級自然保護(hù)區(qū)是亞洲象的最適生境區(qū)域,且受森林?jǐn)_動影響極微,表明保護(hù)區(qū)的設(shè)立和建設(shè)對亞洲象生境保護(hù)起到了明顯成效。但上、下游生境地理隔離,長期孤島式保護(hù)不利于種群繁衍[25]。因此,亟需突破國界線的人為限制,謀劃建立流域尺度的跨境亞洲象保護(hù)區(qū),有計(jì)劃恢復(fù)或修復(fù)亞洲象遷移廊道,改善和恢復(fù)適宜生境連通性,從而促進(jìn)境內(nèi)外亞洲象種群間的交流。

(3)森林?jǐn)_動時(shí)序變化檢測方法在物種生境研究中的推廣：LandTrendr算法用于南滾河流域1989—2018年森林?jǐn)_動檢測,為揭示亞洲象生境時(shí)空變化過程提供了準(zhǔn)確度較高、可靠性較強(qiáng)的數(shù)據(jù)支撐。時(shí)序變化檢測方法能精準(zhǔn)發(fā)現(xiàn)亞洲象生境變化的關(guān)鍵時(shí)間點(diǎn)與空間位置,利于生境變化歸因,也可為生境修復(fù)的先后次序劃定提供依據(jù),從而提高修復(fù)的精準(zhǔn)性和效率。該方法同樣可應(yīng)用于其他物種生境時(shí)空演變過程的診斷,為相應(yīng)物種的生境保護(hù)與恢復(fù)等提供科學(xué)依據(jù)。

4.2 結(jié)論

為探明土地利用變化對亞洲象生境影響的時(shí)空過程及機(jī)制,以南滾河流域?yàn)槔?采用基于Landsat影像的土地利用/土地覆被遙感制圖和基于LandTrendr算法的時(shí)序變化檢測等方法,分析了1988—2018年間流域土地利用變化及其導(dǎo)致的森林?jǐn)_動過程,結(jié)合ENFA模型評定的亞洲象適宜生境格局,探究以天然林?jǐn)_動為主的土地利用變化對亞洲象生境的影響,得出以下結(jié)論：

(1)1988—2018年間,南滾河流域土地利用變化主要表現(xiàn)為天然林向其他土地利用類型轉(zhuǎn)化,其中橡膠林?jǐn)U張導(dǎo)致毀林最為明顯；森林?jǐn)_動強(qiáng)度在中、緬境內(nèi)存在明顯差異,在2001年之前兩國森林?jǐn)_動面積大體相當(dāng),但其后緬甸境內(nèi)流域森林?jǐn)_動面積高于中國境內(nèi)流域。中、緬兩國經(jīng)濟(jì)發(fā)展水平和政策差異及全球橡膠價(jià)格變動是驅(qū)動該流域土地利用變化的重要因素。

(2)南滾河流域以橡膠林種植擴(kuò)張導(dǎo)致毀林為主的土地利用變化,導(dǎo)致亞洲象生境喪失,適宜生境斑塊間連通性降低,破碎化程度較高。

(3)該流域亞洲象生境喪失過程表現(xiàn)為由生境邊緣逐漸向內(nèi)部侵蝕,從而導(dǎo)致上、下游生境斑塊完全隔離。2001年以前,亞洲象最適生境區(qū)受土地利用變化影響較少,生境喪失主要表現(xiàn)為較適生境和邊際生境的減少；2001—2013年間中、緬邊境地區(qū)劇烈的土地利用變化導(dǎo)致流域中游段亞洲象適宜生境的喪失,從而使上、下游適宜生境區(qū)完全被隔離,進(jìn)而影響亞洲象的跨境遷徙和種群交流。