基于地學(xué)信息圖譜的中緬泰老“黃金四角”地區(qū)土地利用變化分析*

陳瀟菡，張 磊 ，史 麗，楊 冉，黃燕壽

(1.中信銀行 昆明東風(fēng)東路支行，云南 昆明 650500；2.泰山學(xué)院 旅游學(xué)院，山東 泰安 271000；3.保山學(xué)院 經(jīng)濟(jì)管理學(xué)院，云南 保山 678000)

作為全球變化研究的關(guān)鍵問(wèn)題，土地利用覆蓋變化(land use cover change，LUCC)已成為自然與人文科學(xué)兩大領(lǐng)域共同搭建的“橋梁工程”[1]，特別自《土地利用/土地覆蓋變化科學(xué)研究計(jì)劃》發(fā)布以來(lái)，對(duì)LUCC 的研究已成為當(dāng)前學(xué)者關(guān)注的熱點(diǎn)之一。經(jīng)過(guò)30 余年的發(fā)展，圍繞土地利用變化的研究理論與方法已被廣泛應(yīng)用于區(qū)域與可持續(xù)發(fā)展、資源環(huán)境承載力及生態(tài)風(fēng)險(xiǎn)評(píng)估等領(lǐng)域，且形成了相對(duì)完整的理論與方法體系。

目前，學(xué)者已對(duì)土地利用變化做了大量研究，隨著GIS 等技術(shù)的成熟，關(guān)于土地利用變化的分析正由單一評(píng)價(jià)向多元數(shù)據(jù)分析轉(zhuǎn)變。廖諶婳[2]和劉曉娜[3]對(duì)中老緬交界地區(qū)土地利用/覆被數(shù)據(jù)提取、分類(lèi)與校正做了分析；林樹(shù)高等[4]、翟香等[5]和耿鑫等[6]分別基于“三生”空間、地形梯度及生態(tài)系統(tǒng)服務(wù)價(jià)值等對(duì)廣西邊境等地區(qū)的土地利用格局及演變進(jìn)行測(cè)度；張涵等[7]對(duì)鄉(xiāng)村轉(zhuǎn)型背景下山—壩系統(tǒng)土地利用景觀(guān)耦合與協(xié)調(diào)態(tài)勢(shì)進(jìn)行了剖析；付慧等[8]和王墾等[9]分別對(duì)區(qū)域協(xié)同發(fā)展與城鎮(zhèn)化等背景下土地利用變動(dòng)機(jī)制及核心要素等進(jìn)行了識(shí)別；孫毅中等[10]和張?jiān)碌萚11]分別基于元胞自動(dòng)機(jī)及CA-Markov 等模型對(duì)江陰市和艾比湖地區(qū)土地的利用未來(lái)格局及變化態(tài)勢(shì)進(jìn)行了模擬?？梢?jiàn)，當(dāng)前學(xué)者對(duì)土地利用變化的研究主要集中于數(shù)據(jù)獲取與分類(lèi)、格局與時(shí)空演變、景觀(guān)測(cè)度、變化機(jī)制以及模擬與預(yù)測(cè)等方面，基于時(shí)空雙維度對(duì)土地利用變化軌跡與漲落態(tài)勢(shì)進(jìn)行的復(fù)合研究不足，且多參照行政區(qū)域選定研究區(qū)，對(duì)打破行政界線(xiàn)且同時(shí)具有地緣政治和經(jīng)濟(jì)特色的跨境地區(qū)研究較少，研究對(duì)象多側(cè)重山區(qū)、水域、沙漠和城鎮(zhèn)等生態(tài)涵養(yǎng)功能脆弱區(qū)，對(duì)熱帶雨林廣布且自然生態(tài)本底正在發(fā)生顯著變遷區(qū)域的研究不足。

近年來(lái)，作為中國(guó)—中南半島經(jīng)濟(jì)走廊重要組分的“黃金四角”地區(qū)社會(huì)經(jīng)濟(jì)活動(dòng)日趨頻繁，生態(tài)環(huán)境問(wèn)題不斷凸顯，處理好資源開(kāi)發(fā)與生態(tài)保護(hù)間的關(guān)系已成為實(shí)現(xiàn)該區(qū)域可持續(xù)發(fā)展急需解決的關(guān)鍵問(wèn)題[12-14]。鑒于此，本研究在提取土地利用本底數(shù)據(jù)的基礎(chǔ)上，引入土地利用動(dòng)態(tài)度、景觀(guān)質(zhì)心平移度及土地利用變化圖譜等模型，對(duì)中緬泰老“黃金四角”地區(qū)2000—2018 年土地利用的變化特征進(jìn)行評(píng)估，以期為摸清其土地資源本底與構(gòu)建生態(tài)保護(hù)路徑提供參考。

1 材料與方法

1.1 研究區(qū)概況

中緬泰老“黃金四角”(簡(jiǎn)稱(chēng)“黃金四角”)地區(qū)位于瀾滄江—湄公河中上游，地勢(shì)北高南低，位于N17°16＇～24°50＇、E98°01＇～103°24＇之間，總面積約20 萬(wàn)km2，地理空間上主要包括中國(guó)云南西雙版納州和普洱市、老撾北部的上寮地區(qū)五省(瑯南塔省、烏多姆賽省、波喬省、豐沙里省和瑯勃拉邦省)、緬甸東北部的東撣邦地區(qū)四縣(景棟縣、孟別縣、大其力縣和孟薩縣)以及泰國(guó)北部的清邁和清萊兩府[15]。

1.2 數(shù)據(jù)來(lái)源與處理

研究所用數(shù)據(jù)主要由土地利用遙感影像及研究區(qū)基礎(chǔ)地理信息數(shù)據(jù)構(gòu)成。其中，遙感影像主要源自美國(guó)地質(zhì)調(diào)查局的Landsat5、7、8 TM/ETM+/OLI 數(shù)據(jù)(分辨率30 m)，選取年份分別為2000年、2009 年和2018 年(Landsat 影像取景行列號(hào)分別為131/043-131/048，130/043-130/047，129/044-129/047，128/046)[15]。首先，對(duì)土地利用數(shù)據(jù)進(jìn)行分類(lèi)提取，參考《土地利用現(xiàn)狀分類(lèi)》(GB/T 21010—2007)并綜合考慮研究區(qū)自然本底識(shí)別，分為水體、耕地、裸土與未利用地、林地、灌叢、園草地及建筑用地等7 個(gè)地類(lèi)；其次，對(duì)研究區(qū)土地利用數(shù)據(jù)進(jìn)行精度評(píng)價(jià)，在分類(lèi)整合野外GPS 調(diào)查采樣數(shù)據(jù)及GoogleEarth 高分辨歷史影像數(shù)據(jù)的基礎(chǔ)上，隨機(jī)選取6 958 個(gè)真實(shí)地表ROI 數(shù)據(jù)對(duì)2000—2018 年“黃金四角”地區(qū)土地利用數(shù)據(jù)進(jìn)行精度評(píng)價(jià)，三期Kappa 系數(shù)分別為0.863、0.879 和0.902，均符合研究標(biāo)準(zhǔn)。

1.3 研究方法與指標(biāo)選取

1.3.1 土地利用動(dòng)態(tài)度分析

土地利用動(dòng)態(tài)度模型可分析研究時(shí)段內(nèi)各土地利用類(lèi)型的變動(dòng)速度與幅度，依據(jù)分析對(duì)象差異可分為單一土地利用動(dòng)態(tài)度與綜合土地利用動(dòng)態(tài)度。

(1) 單一土地利用動(dòng)態(tài)度

該指標(biāo)反映了某個(gè)時(shí)段內(nèi)區(qū)域某一土地利用類(lèi)型的數(shù)量、空間及類(lèi)型組合方式等變化，主要用于測(cè)度研究區(qū)不同土地利用類(lèi)型變化的速度與幅度[16]。計(jì)算公式為：

式中：K為研究時(shí)段內(nèi)某一土地利用類(lèi)型動(dòng)態(tài)度；Ub和Ua分別為研究末期和初期某一土地利用類(lèi)型面積，km2；T為研究時(shí)間，年。

(2) 綜合土地利用動(dòng)態(tài)度

該指標(biāo)反映了在一定時(shí)段內(nèi)研究區(qū)整體土地利用類(lèi)型的變化速率[16-17]，不僅可從全域分析土地利用的空間變化過(guò)程，更能進(jìn)一步反映區(qū)域內(nèi)各因素對(duì)其土地利用變化的綜合影響[16-17]。計(jì)算公式為：

式中：S為綜合土地利用動(dòng)態(tài)度；Si為研究開(kāi)始時(shí)第i類(lèi)土地利用類(lèi)型面積；ΔSi-j為研究時(shí)段內(nèi)第i類(lèi)土地利用類(lèi)型轉(zhuǎn)為其他土地利用類(lèi)型的面積總和[16]；t為研究時(shí)間[17]。S值越大，表明研究區(qū)整體土地利用變動(dòng)速率越快。

1.3.2 景觀(guān)質(zhì)心平移度分析

景觀(guān)質(zhì)心平移度模型是通過(guò)計(jì)算某一景觀(guān)類(lèi)型斑塊的面積加權(quán)質(zhì)心變化來(lái)分析該景觀(guān)類(lèi)型的空間分布變化趨勢(shì)及形成機(jī)制[18-19]，可補(bǔ)充土地利用動(dòng)態(tài)度模型在空間分析上的不足，已廣泛應(yīng)用于土地利用的研究[20]。計(jì)算公式為：

式中：Xt和Yt分別為第t年某用地類(lèi)型質(zhì)心的經(jīng)、緯度坐標(biāo)；n為 斑塊數(shù)量；xi和yi為 第i個(gè)斑塊的質(zhì)心經(jīng)、緯度；Cit為t年第i個(gè)斑塊的面積[18-20]。

土地利用景觀(guān)質(zhì)心移動(dòng)距離測(cè)度公式為：

式中：(Xm，Ym)和(Xt，Yt)分別為第m和t年某用地類(lèi)型的質(zhì)心坐標(biāo)；E為常數(shù)，取111.111 km[21]。

1.3.3 土地利用變化圖譜分析

土地利用變化圖譜可從時(shí)空雙維度對(duì)研究區(qū)土地利用變化態(tài)勢(shì)與軌跡進(jìn)行分析，參照地學(xué)信息圖譜理論并結(jié)合相關(guān)學(xué)者研究，本研究基于ArcGIS 10.3 采用整數(shù)地類(lèi)代碼對(duì)研究區(qū)土地利用類(lèi)型進(jìn)行重新編碼[22]，并基于空間分析模塊對(duì)不同時(shí)段各土地利用類(lèi)型進(jìn)行地圖代數(shù)運(yùn)算，通過(guò)測(cè)度各圖斑單元轉(zhuǎn)換過(guò)程的屬性編碼進(jìn)行最終的土地利用變化圖譜單元[22-23]，以分析其土地利用轉(zhuǎn)移圖譜特征與規(guī)律[24-25]。計(jì)算公式為：

式中：T00、T09和T18分別為2000 年、2009 年和2018年的土地利用數(shù)據(jù)；T(00-09-18)為土地利用的合成圖譜單元[22]。

參照已有學(xué)者的研究，進(jìn)一步對(duì)各合成圖譜單元的變化模式進(jìn)行分類(lèi)：(1)前期變化型，僅在2000—2009 年發(fā)生用地類(lèi)型轉(zhuǎn)變；(2)后期變化型，僅在2009—2018 年發(fā)生用地類(lèi)型轉(zhuǎn)變；(3)穩(wěn)定型，在2000—2018 年均未發(fā)生用地轉(zhuǎn)變；(4)反復(fù)變化型，僅在2009 年發(fā)生用地轉(zhuǎn)變，2000 年和2018 年未發(fā)生變化；(5)持續(xù)變化型，3 個(gè)時(shí)期用地類(lèi)型持續(xù)發(fā)生轉(zhuǎn)變[22]。

1.3.4 土地利用漲落勢(shì)圖譜分析

土地利用轉(zhuǎn)型包括轉(zhuǎn)入和轉(zhuǎn)出2 個(gè)動(dòng)態(tài)過(guò)程，轉(zhuǎn)入是其他用地類(lèi)型轉(zhuǎn)變?yōu)楸居玫仡?lèi)型的過(guò)程，導(dǎo)致該地類(lèi)面積上漲；轉(zhuǎn)出則是本用地類(lèi)型轉(zhuǎn)為其他用地類(lèi)型的過(guò)程，致使該地類(lèi)面積萎縮[23,26]?；贏rcGIS 10.5 對(duì)研究區(qū)不同時(shí)段的土地利用數(shù)據(jù)進(jìn)行疊加處理，并以此為基礎(chǔ)計(jì)算某一地類(lèi)的新增面積和萎縮面積，進(jìn)而生成漲勢(shì)圖譜和落勢(shì)圖譜[23]。

2 結(jié)果與分析

2.1 土地利用時(shí)序變化分析

2.1.1 土地利用現(xiàn)狀分析

“黃金四角”地區(qū)土地利用以林地、耕地和灌叢為主，占土地利用總面積的95%以上。其中，林地為主導(dǎo)型用地類(lèi)型(13.524 萬(wàn)km2)，耕地為次要用地類(lèi)型(3.107 萬(wàn)km2)，灌叢面積(2.472 萬(wàn)km2)略低于耕地；園草地和建筑用地分別占土地利用總面積的2.11%和1.50%，而水體及裸土與未利用地的面積則相對(duì)較小(圖1)。總體上，“黃金四角”地區(qū)形成了以林地、耕地和灌叢為基質(zhì)，水體為廊道，園草地及建筑用地為斑塊的土地利用景觀(guān)格局。

圖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 單一土地利用動(dòng)態(tài)分析

不同時(shí)段間研究區(qū)各土地利用類(lèi)型變動(dòng)速率差異凸顯(圖2)。其中，裸土與未利用地及建筑用地?cái)U(kuò)展速率不斷提升，分別由2000—2009年的9.49%和3.48%增至2009—2018 年的10.58%和5.80%；耕地和園草地?cái)U(kuò)展速率不斷減緩，分別由1.95%和3.04%降至0.99%和1.77%；水體與灌叢用地動(dòng)態(tài)度的變動(dòng)存在時(shí)段差異(均由負(fù)轉(zhuǎn)正)，而林地面積變動(dòng)幅度較小。

圖2 2000—2018 年“黃金四角”地區(qū)單一土地利用動(dòng)態(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 綜合土地利用動(dòng)態(tài)分析

隨著時(shí)間演進(jìn)，研究區(qū)土地利用變動(dòng)速率整體減慢，用地結(jié)構(gòu)穩(wěn)定性大幅增強(qiáng)。2000—2009 年，綜合土地利用動(dòng)態(tài)度為45.99%，土地利用變動(dòng)活躍，地類(lèi)間的轉(zhuǎn)換頻繁；2009—2018 年，綜合土地利用動(dòng)態(tài)度降至10.44%，地類(lèi)間的轉(zhuǎn)類(lèi)規(guī)模與頻率均顯著下降，這與主導(dǎo)型地類(lèi)(林地、灌叢和耕地)轉(zhuǎn)類(lèi)規(guī)模下降以及空間分布穩(wěn)定性提升相關(guān)。

2.2 土地利用空間演變分析

由圖3 可知：(1) 耕地質(zhì)心向東南偏移，位于東部和南部的泰、老兩國(guó)段耕地開(kāi)發(fā)面積占比提升。研究前期，耕地質(zhì)心向正東偏移18.72 km，主要是受老撾段耕地開(kāi)發(fā)與“替代種植”規(guī)模的提升以及泰國(guó)段第一產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化與水稻種植面積增加驅(qū)動(dòng)；研究后期耕地質(zhì)心向西南偏移26.78 km。(2) 林地質(zhì)心整體向正東偏移但階段性差異顯著，老撾段林區(qū)占研究區(qū)林地總面積的比重提升；2000—2009 年，受老撾東南部瑯勃拉邦與烏多姆賽省經(jīng)濟(jì)林面積擴(kuò)大驅(qū)動(dòng)，林地質(zhì)心向東南偏移3.24 km；2009—2018 年林地質(zhì)心向東北偏移3.39 km，其原因在于老撾北部豐沙里省和瑯南塔省林區(qū)恢復(fù)及中、老兩國(guó)交界區(qū)橡膠林面積的擴(kuò)大。(3) 隨著時(shí)間演進(jìn)，灌叢質(zhì)心呈向西北偏移態(tài)勢(shì)，北部中國(guó)段和西部緬甸段的灌叢面積占比提升，而園草地質(zhì)心整體向東南偏移，其原因一方面在于中國(guó)段以茶園為代表的規(guī)?；瘓@地面積的擴(kuò)展，另一方面在于老撾段“替代種植”驅(qū)動(dòng)下人工型園草地(臺(tái)地茶園等)面積的增長(zhǎng)。(4) 隨著時(shí)間演進(jìn)，建筑用地質(zhì)心持續(xù)偏向東南但幅度下降，研究前期以向南部偏移為主，研究后期則主要向東部偏移，一方面表明南部泰國(guó)段持續(xù)為建筑用地的首要分布區(qū)，另一方面表明研究時(shí)段內(nèi)中、老兩國(guó)段建筑用地?cái)U(kuò)展顯著并促進(jìn)了用地質(zhì)心的東移且主要集中在研究后期。(5) 水體質(zhì)心顯著向東北偏移，北部和東部的中、老兩國(guó)段水體面積占比提升；2000—2018年，水體質(zhì)心向東北偏移近70 km，表明中、老兩國(guó)段對(duì)水利設(shè)施及瀾滄江—湄公河流域治理的成效凸顯。(6) 裸土與未利用地質(zhì)心持續(xù)偏向東北但幅度降低，其原因一方面在于中國(guó)段城鄉(xiāng)建設(shè)導(dǎo)致人工型裸土與未利用地增加，另一方面受老撾段“刀耕火種”型耕作方式的時(shí)段性變動(dòng)影響。

圖3 2000—2018 年“黃金四角”地區(qū)各土地利用類(lèi)型空間分布質(zhì)心及變動(dòng)軌跡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

2.3 土地利用轉(zhuǎn)移圖譜分析

土地利用轉(zhuǎn)移圖譜特征及其規(guī)律分析結(jié)果(圖4、表1)表明：(1) 穩(wěn)定型圖譜單元用地面積占比最大(68.7%)，且廣泛散布于整個(gè)研究區(qū)域；林地—林地—林地為面積最大的圖譜模式(113 714.24 km2)，這體現(xiàn)了典型熱帶雨林區(qū)土地利用的分布及變動(dòng)特征。(2) 前期變化型圖譜單元總面積為45 139.89 km2，最大面積圖譜模式為灌叢—林地—林地，占該圖譜類(lèi)型總面積的29.9%，且在東部老撾段集聚態(tài)勢(shì)明顯。相較于人為干擾，自然本底因素在研究區(qū)前期土地利用變化中的作用更為明顯。(3) 后期變化型圖譜單元總面積為9 094.99 km2且集中于南北兩側(cè)的中、泰兩國(guó)段，林地—林地—耕地為最大面積圖譜變化模式(2 365.21 km2)，這源于研究后期過(guò)度砍伐和不合理利用等造成的林地退化以及隨著土地利用程度加深林地與耕地間的高頻率轉(zhuǎn)換。此類(lèi)圖譜變動(dòng)的區(qū)域分異也體現(xiàn)出研究后期各國(guó)段土地利用開(kāi)發(fā)強(qiáng)度差異已經(jīng)凸顯。(4) 持續(xù)變化型圖譜單元數(shù)量最多(211 個(gè))，但面積占比較小，空間分布破碎化程度高；最大面積圖譜模式為林地—耕地—灌叢(799.09 km2)，表明林地、耕地及灌叢三大主導(dǎo)型地類(lèi)間的轉(zhuǎn)化直接影響研究區(qū)用地圖譜的變動(dòng)，同時(shí)，耕地的持續(xù)轉(zhuǎn)類(lèi)也體現(xiàn)出該區(qū)域“刀耕火種”及棄耕等現(xiàn)象仍廣泛存在，土地利用效率的提升與綜合整治需加強(qiáng)。(5) 反復(fù)變化型圖譜單元用地面積占比最小(1.72%)且顯著集中在南部泰國(guó)段，耕地—灌叢—耕地和耕地—林地—耕地為其面積較大的圖譜模式，這與該區(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 土地利用漲落勢(shì)圖譜分析

2.4.1 土地利用漲勢(shì)圖譜

由表2 可知：2000—2009 年研究區(qū)土地利用未變區(qū)域面積為14.6 萬(wàn)km2，變化區(qū)域面積為5.44 萬(wàn)km2；變化區(qū)域中，新增林地面積最大，占新增區(qū)域總面積的40.20%；其次為耕地和灌叢，面積分別為25.30%和24.10%；相較而言，新增水體、裸土與未利用地及建筑用地的面積則相對(duì)較小。對(duì)比發(fā)現(xiàn)，研究區(qū)土地利用數(shù)量與結(jié)構(gòu)轉(zhuǎn)型存在一定國(guó)別差異，老撾、緬甸、泰國(guó)和中國(guó)段未變區(qū)域面積占比分別為32.38%、18.14%、15.47%和34.01%，表明中、老兩國(guó)段土地利用結(jié)構(gòu)較為穩(wěn)定，而泰、緬兩國(guó)段的土地利用變動(dòng)較為活躍。

表2 “黃金四角”地區(qū)各時(shí)序單元土地利用漲勢(shì)圖譜結(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ū)土地利用變動(dòng)活躍度有所下降，未變土地利用面積增為18.1 萬(wàn)km2；新增耕地面積最大，占新增用地總面積的39.50%，其次為灌叢，占新增用地總面積的29.10%；相較而言，新增建筑用地的面積顯著提升，而新增園草地面積有所下降。中、緬、泰、老4 國(guó)段未變區(qū)域面積占比與前一時(shí)段無(wú)明顯差異，值得注意的是，2009—2018 年間，4 國(guó)段新增耕地面積占比均顯著提升。

表3 “黃金四角”地區(qū)各時(shí)序單元土地利用漲勢(shì)圖譜結(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 土地利用落勢(shì)圖譜

由表4 可知：2000—2009 年研究區(qū)灌叢面積萎縮最大，為20 571.04 km2，其次為林地，萎縮面積占變動(dòng)區(qū)域總面積的36.5%；相較而言，水體和裸土與未利用地的萎縮面積較小。中、緬、泰、老4 國(guó)段主要萎縮地類(lèi)存在較高的相似性，均呈以灌叢、林地及耕地為主的轉(zhuǎn)出特征。

表4 “黃金四角”地區(qū)各時(shí)序單元土地利用落勢(shì)圖譜結(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ū)萎縮面積最大的地類(lèi)為林地(45.80%)，且規(guī)模顯著高于其他地類(lèi)；相較于上一時(shí)期，該時(shí)段研究區(qū)土地利用萎縮總面積顯著減小且在4 個(gè)國(guó)段廣泛存在，其原因主要在于隨經(jīng)濟(jì)與社會(huì)結(jié)構(gòu)的穩(wěn)定，作為研究區(qū)土地利用基質(zhì)的林地、灌叢及耕地間轉(zhuǎn)類(lèi)面積減小，土地利用結(jié)構(gòu)的穩(wěn)定性增強(qiáng)。

表5 “黃金四角”地區(qū)各時(shí)序單元土地利用落勢(shì)圖譜結(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

3 討論

對(duì)“黃金四角”地區(qū)土地利用變化的分析不僅有助于為該區(qū)域資源合作開(kāi)發(fā)、協(xié)同治理、共同保護(hù)的主體價(jià)值取向、協(xié)調(diào)機(jī)制與開(kāi)發(fā)模式提供借鑒，也可為提升中國(guó)在該區(qū)域的輻射與示范效應(yīng)提供參考[27]。本研究在傳統(tǒng)土地利用動(dòng)態(tài)與景觀(guān)質(zhì)心平移分析基礎(chǔ)上，基于整數(shù)地類(lèi)代碼構(gòu)建土地利用變化圖譜，可從時(shí)空雙維度對(duì)土地利用變化軌跡與漲落態(tài)勢(shì)進(jìn)行復(fù)合研究[28]，并為打破行政界線(xiàn)且同時(shí)具有地緣政治、經(jīng)濟(jì)特色的跨境地區(qū)研究提供本底支撐。本研究中的土地利用時(shí)空演變分析與劉曉娜等[13]的研究結(jié)果相似，補(bǔ)充了張磊等[20,29]對(duì)中老緬交界區(qū)及緬甸東撣邦地區(qū)土地利用景觀(guān)及生態(tài)狀況研究，并可為封順等[30]及趙姝嵐等[31]關(guān)于跨境安全復(fù)合體及替代種植的研究提供數(shù)據(jù)支持。“黃金四角”地區(qū)土地利用變化態(tài)勢(shì)表明該區(qū)形成了林地、耕地和灌叢為基質(zhì)，水體為廊道，園草地及建筑用地為斑塊的土地利用景觀(guān)特征，這驗(yàn)證了LIU 等[14]對(duì)老撾北部、緬甸東北部和泰國(guó)北部等地土地利用結(jié)構(gòu)的論述?；谕恋乩米兓瘓D譜的分析則發(fā)現(xiàn)：社會(huì)、經(jīng)濟(jì)狀況與研究區(qū)土地利用效率與變化活躍度緊密相關(guān)，土地利用漲落態(tài)勢(shì)更是直接受林地、灌叢及耕地間轉(zhuǎn)類(lèi)的影響。隨著時(shí)間演進(jìn)，“黃金四角”地區(qū)土地利用變化的國(guó)別分異逐步凸顯，中、泰兩國(guó)段土地利用的活躍度與合理性整體高于老、緬兩國(guó)段。

4 結(jié)論

基于土地利用動(dòng)態(tài)度、景觀(guān)質(zhì)心平移度及地圖代數(shù)運(yùn)算等模型對(duì)“黃金四角”地區(qū)土地利用變化的測(cè)度與評(píng)價(jià)不僅實(shí)現(xiàn)了對(duì)該區(qū)域土地資源利用子系統(tǒng)結(jié)構(gòu)組成、時(shí)空分布、變化規(guī)律及國(guó)別分異的量化測(cè)度，更有助于構(gòu)建跨境/跨流域合作與共生地區(qū)土地資源本底專(zhuān)題數(shù)據(jù)庫(kù)。強(qiáng)化對(duì)代表性區(qū)域與地類(lèi)(瀾滄江—湄公河沿線(xiàn)、中緬泰老四國(guó)交界處和橡膠林等)的分析研究、多元剖析土地利用變動(dòng)機(jī)制并構(gòu)建促進(jìn)土地資源合理利用的優(yōu)化路徑是課題組后續(xù)研究的重點(diǎn)。