陳瀟菡,張 磊 ,史 麗,楊 冉,黃燕壽
(1.中信銀行 昆明東風(fēng)東路支行,云南 昆明 650500;2.泰山學(xué)院 旅游學(xué)院,山東 泰安 271000;3.保山學(xué)院 經(jīng)濟(jì)管理學(xué)院,云南 保山 678000)
作為全球變化研究的關(guān)鍵問題,土地利用覆蓋變化(land use cover change,LUCC)已成為自然與人文科學(xué)兩大領(lǐng)域共同搭建的“橋梁工程”[1],特別自《土地利用/土地覆蓋變化科學(xué)研究計劃》發(fā)布以來,對LUCC 的研究已成為當(dāng)前學(xué)者關(guān)注的熱點(diǎn)之一。經(jīng)過30 余年的發(fā)展,圍繞土地利用變化的研究理論與方法已被廣泛應(yīng)用于區(qū)域與可持續(xù)發(fā)展、資源環(huán)境承載力及生態(tài)風(fēng)險評估等領(lǐng)域,且形成了相對完整的理論與方法體系。
目前,學(xué)者已對土地利用變化做了大量研究,隨著GIS 等技術(shù)的成熟,關(guān)于土地利用變化的分析正由單一評價向多元數(shù)據(jù)分析轉(zhuǎn)變。廖諶婳[2]和劉曉娜[3]對中老緬交界地區(qū)土地利用/覆被數(shù)據(jù)提取、分類與校正做了分析;林樹高等[4]、翟香等[5]和耿鑫等[6]分別基于“三生”空間、地形梯度及生態(tài)系統(tǒng)服務(wù)價值等對廣西邊境等地區(qū)的土地利用格局及演變進(jìn)行測度;張涵等[7]對鄉(xiāng)村轉(zhuǎn)型背景下山—壩系統(tǒng)土地利用景觀耦合與協(xié)調(diào)態(tài)勢進(jìn)行了剖析;付慧等[8]和王墾等[9]分別對區(qū)域協(xié)同發(fā)展與城鎮(zhèn)化等背景下土地利用變動機(jī)制及核心要素等進(jìn)行了識別;孫毅中等[10]和張月等[11]分別基于元胞自動機(jī)及CA-Markov 等模型對江陰市和艾比湖地區(qū)土地的利用未來格局及變化態(tài)勢進(jìn)行了模擬。可見,當(dāng)前學(xué)者對土地利用變化的研究主要集中于數(shù)據(jù)獲取與分類、格局與時空演變、景觀測度、變化機(jī)制以及模擬與預(yù)測等方面,基于時空雙維度對土地利用變化軌跡與漲落態(tài)勢進(jìn)行的復(fù)合研究不足,且多參照行政區(qū)域選定研究區(qū),對打破行政界線且同時具有地緣政治和經(jīng)濟(jì)特色的跨境地區(qū)研究較少,研究對象多側(cè)重山區(qū)、水域、沙漠和城鎮(zhèn)等生態(tài)涵養(yǎng)功能脆弱區(qū),對熱帶雨林廣布且自然生態(tài)本底正在發(fā)生顯著變遷區(qū)域的研究不足。
近年來,作為中國—中南半島經(jīng)濟(jì)走廊重要組分的“黃金四角”地區(qū)社會經(jīng)濟(jì)活動日趨頻繁,生態(tài)環(huán)境問題不斷凸顯,處理好資源開發(fā)與生態(tài)保護(hù)間的關(guān)系已成為實現(xiàn)該區(qū)域可持續(xù)發(fā)展急需解決的關(guān)鍵問題[12-14]。鑒于此,本研究在提取土地利用本底數(shù)據(jù)的基礎(chǔ)上,引入土地利用動態(tài)度、景觀質(zhì)心平移度及土地利用變化圖譜等模型,對中緬泰老“黃金四角”地區(qū)2000—2018 年土地利用的變化特征進(jìn)行評估,以期為摸清其土地資源本底與構(gòu)建生態(tài)保護(hù)路徑提供參考。
中緬泰老“黃金四角”(簡稱“黃金四角”)地區(qū)位于瀾滄江—湄公河中上游,地勢北高南低,位于N17°16'~24°50'、E98°01'~103°24'之間,總面積約20 萬km2,地理空間上主要包括中國云南西雙版納州和普洱市、老撾北部的上寮地區(qū)五省(瑯南塔省、烏多姆賽省、波喬省、豐沙里省和瑯勃拉邦省)、緬甸東北部的東撣邦地區(qū)四縣(景棟縣、孟別縣、大其力縣和孟薩縣)以及泰國北部的清邁和清萊兩府[15]。
研究所用數(shù)據(jù)主要由土地利用遙感影像及研究區(qū)基礎(chǔ)地理信息數(shù)據(jù)構(gòu)成。其中,遙感影像主要源自美國地質(zhì)調(diào)查局的Landsat5、7、8 TM/ETM+/OLI 數(shù)據(jù)(分辨率30 m),選取年份分別為2000年、2009 年和2018 年(Landsat 影像取景行列號分別為131/043-131/048,130/043-130/047,129/044-129/047,128/046)[15]。首先,對土地利用數(shù)據(jù)進(jìn)行分類提取,參考《土地利用現(xiàn)狀分類》(GB/T 21010—2007)并綜合考慮研究區(qū)自然本底識別,分為水體、耕地、裸土與未利用地、林地、灌叢、園草地及建筑用地等7 個地類;其次,對研究區(qū)土地利用數(shù)據(jù)進(jìn)行精度評價,在分類整合野外GPS 調(diào)查采樣數(shù)據(jù)及GoogleEarth 高分辨歷史影像數(shù)據(jù)的基礎(chǔ)上,隨機(jī)選取6 958 個真實地表ROI 數(shù)據(jù)對2000—2018 年“黃金四角”地區(qū)土地利用數(shù)據(jù)進(jìn)行精度評價,三期Kappa 系數(shù)分別為0.863、0.879 和0.902,均符合研究標(biāo)準(zhǔn)。
1.3.1 土地利用動態(tài)度分析
土地利用動態(tài)度模型可分析研究時段內(nèi)各土地利用類型的變動速度與幅度,依據(jù)分析對象差異可分為單一土地利用動態(tài)度與綜合土地利用動態(tài)度。
(1) 單一土地利用動態(tài)度
該指標(biāo)反映了某個時段內(nèi)區(qū)域某一土地利用類型的數(shù)量、空間及類型組合方式等變化,主要用于測度研究區(qū)不同土地利用類型變化的速度與幅度[16]。計算公式為:
式中:K為研究時段內(nèi)某一土地利用類型動態(tài)度;Ub和Ua分別為研究末期和初期某一土地利用類型面積,km2;T為研究時間,年。
(2) 綜合土地利用動態(tài)度
該指標(biāo)反映了在一定時段內(nèi)研究區(qū)整體土地利用類型的變化速率[16-17],不僅可從全域分析土地利用的空間變化過程,更能進(jìn)一步反映區(qū)域內(nèi)各因素對其土地利用變化的綜合影響[16-17]。計算公式為:
式中:S為綜合土地利用動態(tài)度;Si為研究開始時第i類土地利用類型面積;ΔSi-j為研究時段內(nèi)第i類土地利用類型轉(zhuǎn)為其他土地利用類型的面積總和[16];t為研究時間[17]。S值越大,表明研究區(qū)整體土地利用變動速率越快。
1.3.2 景觀質(zhì)心平移度分析
景觀質(zhì)心平移度模型是通過計算某一景觀類型斑塊的面積加權(quán)質(zhì)心變化來分析該景觀類型的空間分布變化趨勢及形成機(jī)制[18-19],可補(bǔ)充土地利用動態(tài)度模型在空間分析上的不足,已廣泛應(yīng)用于土地利用的研究[20]。計算公式為:
式中:Xt和Yt分別為第t年某用地類型質(zhì)心的經(jīng)、緯度坐標(biāo);n為 斑塊數(shù)量;xi和yi為 第i個斑塊的質(zhì)心經(jīng)、緯度;Cit為t年第i個斑塊的面積[18-20]。
土地利用景觀質(zhì)心移動距離測度公式為:
式中:(Xm,Ym)和(Xt,Yt)分別為第m和t年某用地類型的質(zhì)心坐標(biāo);E為常數(shù),取111.111 km[21]。
1.3.3 土地利用變化圖譜分析
土地利用變化圖譜可從時空雙維度對研究區(qū)土地利用變化態(tài)勢與軌跡進(jìn)行分析,參照地學(xué)信息圖譜理論并結(jié)合相關(guān)學(xué)者研究,本研究基于ArcGIS 10.3 采用整數(shù)地類代碼對研究區(qū)土地利用類型進(jìn)行重新編碼[22],并基于空間分析模塊對不同時段各土地利用類型進(jìn)行地圖代數(shù)運(yùn)算,通過測度各圖斑單元轉(zhuǎn)換過程的屬性編碼進(jìn)行最終的土地利用變化圖譜單元[22-23],以分析其土地利用轉(zhuǎn)移圖譜特征與規(guī)律[24-25]。計算公式為:
式中:T00、T09和T18分別為2000 年、2009 年和2018年的土地利用數(shù)據(jù);T(00-09-18)為土地利用的合成圖譜單元[22]。
參照已有學(xué)者的研究,進(jìn)一步對各合成圖譜單元的變化模式進(jìn)行分類:(1)前期變化型,僅在2000—2009 年發(fā)生用地類型轉(zhuǎn)變;(2)后期變化型,僅在2009—2018 年發(fā)生用地類型轉(zhuǎn)變;(3)穩(wěn)定型,在2000—2018 年均未發(fā)生用地轉(zhuǎn)變;(4)反復(fù)變化型,僅在2009 年發(fā)生用地轉(zhuǎn)變,2000 年和2018 年未發(fā)生變化;(5)持續(xù)變化型,3 個時期用地類型持續(xù)發(fā)生轉(zhuǎn)變[22]。
1.3.4 土地利用漲落勢圖譜分析
土地利用轉(zhuǎn)型包括轉(zhuǎn)入和轉(zhuǎn)出2 個動態(tài)過程,轉(zhuǎn)入是其他用地類型轉(zhuǎn)變?yōu)楸居玫仡愋偷倪^程,導(dǎo)致該地類面積上漲;轉(zhuǎn)出則是本用地類型轉(zhuǎn)為其他用地類型的過程,致使該地類面積萎縮[23,26]?;贏rcGIS 10.5 對研究區(qū)不同時段的土地利用數(shù)據(jù)進(jìn)行疊加處理,并以此為基礎(chǔ)計算某一地類的新增面積和萎縮面積,進(jìn)而生成漲勢圖譜和落勢圖譜[23]。
2.1.1 土地利用現(xiàn)狀分析
“黃金四角”地區(qū)土地利用以林地、耕地和灌叢為主,占土地利用總面積的95%以上。其中,林地為主導(dǎo)型用地類型(13.524 萬km2),耕地為次要用地類型(3.107 萬km2),灌叢面積(2.472 萬km2)略低于耕地;園草地和建筑用地分別占土地利用總面積的2.11%和1.50%,而水體及裸土與未利用地的面積則相對較小(圖1)。總體上,“黃金四角”地區(qū)形成了以林地、耕地和灌叢為基質(zhì),水體為廊道,園草地及建筑用地為斑塊的土地利用景觀格局。
圖1 “黃金四角”地區(qū)土地利用現(xiàn)狀及結(jié)構(gòu)特征Fig.1 Land use status and structural characteristics of the“Golden Four Corners” of China,Myanmar,Thailand and Laos
2.1.2 單一土地利用動態(tài)分析
不同時段間研究區(qū)各土地利用類型變動速率差異凸顯(圖2)。其中,裸土與未利用地及建筑用地擴(kuò)展速率不斷提升,分別由2000—2009年的9.49%和3.48%增至2009—2018 年的10.58%和5.80%;耕地和園草地擴(kuò)展速率不斷減緩,分別由1.95%和3.04%降至0.99%和1.77%;水體與灌叢用地動態(tài)度的變動存在時段差異(均由負(fù)轉(zhuǎn)正),而林地面積變動幅度較小。
圖2 2000—2018 年“黃金四角”地區(qū)單一土地利用動態(tài)度Fig.2 Dynamic degree of single land use in the “Golden Four Corners” of China,Myanmar,Thailand and Laos from 2000 to 2018
2.1.3 綜合土地利用動態(tài)分析
隨著時間演進(jìn),研究區(qū)土地利用變動速率整體減慢,用地結(jié)構(gòu)穩(wěn)定性大幅增強(qiáng)。2000—2009 年,綜合土地利用動態(tài)度為45.99%,土地利用變動活躍,地類間的轉(zhuǎn)換頻繁;2009—2018 年,綜合土地利用動態(tài)度降至10.44%,地類間的轉(zhuǎn)類規(guī)模與頻率均顯著下降,這與主導(dǎo)型地類(林地、灌叢和耕地)轉(zhuǎn)類規(guī)模下降以及空間分布穩(wěn)定性提升相關(guān)。
由圖3 可知:(1) 耕地質(zhì)心向東南偏移,位于東部和南部的泰、老兩國段耕地開發(fā)面積占比提升。研究前期,耕地質(zhì)心向正東偏移18.72 km,主要是受老撾段耕地開發(fā)與“替代種植”規(guī)模的提升以及泰國段第一產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化與水稻種植面積增加驅(qū)動;研究后期耕地質(zhì)心向西南偏移26.78 km。(2) 林地質(zhì)心整體向正東偏移但階段性差異顯著,老撾段林區(qū)占研究區(qū)林地總面積的比重提升;2000—2009 年,受老撾東南部瑯勃拉邦與烏多姆賽省經(jīng)濟(jì)林面積擴(kuò)大驅(qū)動,林地質(zhì)心向東南偏移3.24 km;2009—2018 年林地質(zhì)心向東北偏移3.39 km,其原因在于老撾北部豐沙里省和瑯南塔省林區(qū)恢復(fù)及中、老兩國交界區(qū)橡膠林面積的擴(kuò)大。(3) 隨著時間演進(jìn),灌叢質(zhì)心呈向西北偏移態(tài)勢,北部中國段和西部緬甸段的灌叢面積占比提升,而園草地質(zhì)心整體向東南偏移,其原因一方面在于中國段以茶園為代表的規(guī)?;瘓@地面積的擴(kuò)展,另一方面在于老撾段“替代種植”驅(qū)動下人工型園草地(臺地茶園等)面積的增長。(4) 隨著時間演進(jìn),建筑用地質(zhì)心持續(xù)偏向東南但幅度下降,研究前期以向南部偏移為主,研究后期則主要向東部偏移,一方面表明南部泰國段持續(xù)為建筑用地的首要分布區(qū),另一方面表明研究時段內(nèi)中、老兩國段建筑用地擴(kuò)展顯著并促進(jìn)了用地質(zhì)心的東移且主要集中在研究后期。(5) 水體質(zhì)心顯著向東北偏移,北部和東部的中、老兩國段水體面積占比提升;2000—2018年,水體質(zhì)心向東北偏移近70 km,表明中、老兩國段對水利設(shè)施及瀾滄江—湄公河流域治理的成效凸顯。(6) 裸土與未利用地質(zhì)心持續(xù)偏向東北但幅度降低,其原因一方面在于中國段城鄉(xiāng)建設(shè)導(dǎo)致人工型裸土與未利用地增加,另一方面受老撾段“刀耕火種”型耕作方式的時段性變動影響。
圖3 2000—2018 年“黃金四角”地區(qū)各土地利用類型空間分布質(zhì)心及變動軌跡Fig.3 The centroid of spatial distribution and change trajectory of various land use types in the “Golden Four Corners” of China,Myanmar,Thailand and Laos from 2000 to 2018
土地利用轉(zhuǎn)移圖譜特征及其規(guī)律分析結(jié)果(圖4、表1)表明:(1) 穩(wěn)定型圖譜單元用地面積占比最大(68.7%),且廣泛散布于整個研究區(qū)域;林地—林地—林地為面積最大的圖譜模式(113 714.24 km2),這體現(xiàn)了典型熱帶雨林區(qū)土地利用的分布及變動特征。(2) 前期變化型圖譜單元總面積為45 139.89 km2,最大面積圖譜模式為灌叢—林地—林地,占該圖譜類型總面積的29.9%,且在東部老撾段集聚態(tài)勢明顯。相較于人為干擾,自然本底因素在研究區(qū)前期土地利用變化中的作用更為明顯。(3) 后期變化型圖譜單元總面積為9 094.99 km2且集中于南北兩側(cè)的中、泰兩國段,林地—林地—耕地為最大面積圖譜變化模式(2 365.21 km2),這源于研究后期過度砍伐和不合理利用等造成的林地退化以及隨著土地利用程度加深林地與耕地間的高頻率轉(zhuǎn)換。此類圖譜變動的區(qū)域分異也體現(xiàn)出研究后期各國段土地利用開發(fā)強(qiáng)度差異已經(jīng)凸顯。(4) 持續(xù)變化型圖譜單元數(shù)量最多(211 個),但面積占比較小,空間分布破碎化程度高;最大面積圖譜模式為林地—耕地—灌叢(799.09 km2),表明林地、耕地及灌叢三大主導(dǎo)型地類間的轉(zhuǎn)化直接影響研究區(qū)用地圖譜的變動,同時,耕地的持續(xù)轉(zhuǎn)類也體現(xiàn)出該區(qū)域“刀耕火種”及棄耕等現(xiàn)象仍廣泛存在,土地利用效率的提升與綜合整治需加強(qiáng)。(5) 反復(fù)變化型圖譜單元用地面積占比最小(1.72%)且顯著集中在南部泰國段,耕地—灌叢—耕地和耕地—林地—耕地為其面積較大的圖譜模式,這與該區(qū)域獨(dú)特的休耕、輪耕等耕作制度相契合。
表1 2000—2018 年“黃金四角”地區(qū)土地利用變化圖譜分析Tab.1 Analysis of land use change atlas in the “Golden Four Corners” of China,Myanmar,Thailand and Laos from 2000 to 2018
圖4 2000—2018 年“黃金四角”地區(qū)土地利用變化合成圖譜Fig.4 Synthetic atlas of land use change in the “Golden Four Corners” of China,Myanmar,Thailand and Laos from 2000 to 2018
2.4.1 土地利用漲勢圖譜
由表2 可知:2000—2009 年研究區(qū)土地利用未變區(qū)域面積為14.6 萬km2,變化區(qū)域面積為5.44 萬km2;變化區(qū)域中,新增林地面積最大,占新增區(qū)域總面積的40.20%;其次為耕地和灌叢,面積分別為25.30%和24.10%;相較而言,新增水體、裸土與未利用地及建筑用地的面積則相對較小。對比發(fā)現(xiàn),研究區(qū)土地利用數(shù)量與結(jié)構(gòu)轉(zhuǎn)型存在一定國別差異,老撾、緬甸、泰國和中國段未變區(qū)域面積占比分別為32.38%、18.14%、15.47%和34.01%,表明中、老兩國段土地利用結(jié)構(gòu)較為穩(wěn)定,而泰、緬兩國段的土地利用變動較為活躍。
表2 “黃金四角”地區(qū)各時序單元土地利用漲勢圖譜結(jié)構(gòu)表(2000—2009 年)Tab.2 Structure table of land use rise atlas of each time series unit in the “Golden Four Corners” of China,Myanmar,Thailand and Laos (2000-2009) km2
由表3 可知:2009—2018 年,研究區(qū)土地利用變動活躍度有所下降,未變土地利用面積增為18.1 萬km2;新增耕地面積最大,占新增用地總面積的39.50%,其次為灌叢,占新增用地總面積的29.10%;相較而言,新增建筑用地的面積顯著提升,而新增園草地面積有所下降。中、緬、泰、老4 國段未變區(qū)域面積占比與前一時段無明顯差異,值得注意的是,2009—2018 年間,4 國段新增耕地面積占比均顯著提升。
表3 “黃金四角”地區(qū)各時序單元土地利用漲勢圖譜結(jié)構(gòu)表(2009—2018 年)Tab.3 Structure table of land use rise atlas of each time series unit in the “Golden Four Corners” of China,Myanmar,Thailand and Laos (2009-2018) km2
2.4.2 土地利用落勢圖譜
由表4 可知:2000—2009 年研究區(qū)灌叢面積萎縮最大,為20 571.04 km2,其次為林地,萎縮面積占變動區(qū)域總面積的36.5%;相較而言,水體和裸土與未利用地的萎縮面積較小。中、緬、泰、老4 國段主要萎縮地類存在較高的相似性,均呈以灌叢、林地及耕地為主的轉(zhuǎn)出特征。
表4 “黃金四角”地區(qū)各時序單元土地利用落勢圖譜結(jié)構(gòu)表(2000—2009 年)Tab.4 Structure table of land use fall atlas of each time series unit in the “Golden Four Corners” of China,Myanmar,Thailand and Laos (2000-2009) km2
由表5 可知:2009—2018 年研究區(qū)萎縮面積最大的地類為林地(45.80%),且規(guī)模顯著高于其他地類;相較于上一時期,該時段研究區(qū)土地利用萎縮總面積顯著減小且在4 個國段廣泛存在,其原因主要在于隨經(jīng)濟(jì)與社會結(jié)構(gòu)的穩(wěn)定,作為研究區(qū)土地利用基質(zhì)的林地、灌叢及耕地間轉(zhuǎn)類面積減小,土地利用結(jié)構(gòu)的穩(wěn)定性增強(qiáng)。
表5 “黃金四角”地區(qū)各時序單元土地利用落勢圖譜結(jié)構(gòu)表(2009—2018 年)Tab.5 Structure table of land use fall atlas of each time series unit in the “Golden Four Corners” of China,Myanmar,Thailand and Laos (2009-2018) km2
對“黃金四角”地區(qū)土地利用變化的分析不僅有助于為該區(qū)域資源合作開發(fā)、協(xié)同治理、共同保護(hù)的主體價值取向、協(xié)調(diào)機(jī)制與開發(fā)模式提供借鑒,也可為提升中國在該區(qū)域的輻射與示范效應(yīng)提供參考[27]。本研究在傳統(tǒng)土地利用動態(tài)與景觀質(zhì)心平移分析基礎(chǔ)上,基于整數(shù)地類代碼構(gòu)建土地利用變化圖譜,可從時空雙維度對土地利用變化軌跡與漲落態(tài)勢進(jìn)行復(fù)合研究[28],并為打破行政界線且同時具有地緣政治、經(jīng)濟(jì)特色的跨境地區(qū)研究提供本底支撐。本研究中的土地利用時空演變分析與劉曉娜等[13]的研究結(jié)果相似,補(bǔ)充了張磊等[20,29]對中老緬交界區(qū)及緬甸東撣邦地區(qū)土地利用景觀及生態(tài)狀況研究,并可為封順等[30]及趙姝嵐等[31]關(guān)于跨境安全復(fù)合體及替代種植的研究提供數(shù)據(jù)支持?!包S金四角”地區(qū)土地利用變化態(tài)勢表明該區(qū)形成了林地、耕地和灌叢為基質(zhì),水體為廊道,園草地及建筑用地為斑塊的土地利用景觀特征,這驗證了LIU 等[14]對老撾北部、緬甸東北部和泰國北部等地土地利用結(jié)構(gòu)的論述。基于土地利用變化圖譜的分析則發(fā)現(xiàn):社會、經(jīng)濟(jì)狀況與研究區(qū)土地利用效率與變化活躍度緊密相關(guān),土地利用漲落態(tài)勢更是直接受林地、灌叢及耕地間轉(zhuǎn)類的影響。隨著時間演進(jìn),“黃金四角”地區(qū)土地利用變化的國別分異逐步凸顯,中、泰兩國段土地利用的活躍度與合理性整體高于老、緬兩國段。
基于土地利用動態(tài)度、景觀質(zhì)心平移度及地圖代數(shù)運(yùn)算等模型對“黃金四角”地區(qū)土地利用變化的測度與評價不僅實現(xiàn)了對該區(qū)域土地資源利用子系統(tǒng)結(jié)構(gòu)組成、時空分布、變化規(guī)律及國別分異的量化測度,更有助于構(gòu)建跨境/跨流域合作與共生地區(qū)土地資源本底專題數(shù)據(jù)庫。強(qiáng)化對代表性區(qū)域與地類(瀾滄江—湄公河沿線、中緬泰老四國交界處和橡膠林等)的分析研究、多元剖析土地利用變動機(jī)制并構(gòu)建促進(jìn)土地資源合理利用的優(yōu)化路徑是課題組后續(xù)研究的重點(diǎn)。