● 左天惠，陸飛揚(yáng)，黃 妤

（1.廣西自然資源調(diào)查監(jiān)測院，廣西 南寧 531002；2.梧州市自然資源局，廣西 梧州 543002）

隨著漓江流域內(nèi)經(jīng)濟(jì)活動規(guī)模的不斷擴(kuò)大，流域內(nèi)生態(tài)環(huán)境有逐步失衡的趨勢[1]。為全面了解漓江流域土地利用時(shí)空的變化，研究小組利用“3S”技術(shù)進(jìn)行動態(tài)監(jiān)測，通過提取1980—2020 年間共8 期各類土地利用時(shí)空變化數(shù)據(jù)后，運(yùn)用土地利用動態(tài)和轉(zhuǎn)移矩陣方法，分析漓江流域土地利用時(shí)空變化的特征及驅(qū)動力。研究結(jié)果表明：（1）流域最主要的土地利用類型為林地，多年平均占比為 62.21%，其次為耕地、草地、建設(shè)用地和水域；（2）1980—2020 年以來，流域各土地利用類型整體變化幅度不大，主要表現(xiàn)為耕地、林地和草地面積減小，建設(shè)用地和水域面積增加；（3）1980—2020 年期間流域內(nèi)各土地利用類型均發(fā)生不同程度變化，林地面積先減后增再減，整體變化不大。耕地面積先減少后增加，再急劇減少。水域面積先增加后減少再增加。草地面積先增加后減少。建設(shè)用地面積持續(xù)增加，且增加趨勢明顯。

1 流域概況及數(shù)據(jù)處理

1.1 流域概況

漓江流域位于廣西東北部、南嶺山系西南端，是珠江水系西江（珠江上游干流）的支流桂江上游河段的通稱，發(fā)源于廣西興安縣華江鄉(xiāng)（四百里越城嶺山脈的最高峰）貓兒山東北面海拔1.732 km（黃?；妫┑睦仙浇缒蟼?cè)，流域范圍主要包括象山、七星、疊彩、秀峰、雁山、臨桂等6 個(gè)區(qū)以及靈川縣、興安縣、陽朔縣、平樂縣部分區(qū)域，第七次全國人口普查數(shù)據(jù)及《2021 廣西統(tǒng)計(jì)年鑒》數(shù)據(jù)顯示，2020 年，該流域范圍內(nèi)總?cè)丝跒?06.97 萬，生產(chǎn)總值為1 454.78 億元，流域總面積為6 352.78 km2，干流全長214 km。整體地勢北高南低，以漓江為中軸線，呈南北向狹長帶狀分布，北部為碎屑巖中低山地貌，中南部主要為低山、丘陵、巖溶地貌，具有典型的喀斯特地貌特征。漓江流域?qū)儆谀蟻啛釒駶櫦撅L(fēng)氣候區(qū)，年平均氣溫在20℃左右，年降雨量超過1 000 mm，降雨空間分布特征為由北向南遞減[2]。

1.2 流域基礎(chǔ)數(shù)據(jù)來源

（1）流域的基礎(chǔ)數(shù)據(jù)主要來源于漓江流域SRTM_V4 DEM 90 m、ASTER GDEM 30 m、ALOS DEM 12.5 m 分辨率數(shù)字高程模型（來源自USGS、NASA 和JAXA），交通、水系、居民點(diǎn)、行政駐地、行政區(qū)劃界限和注記等。其中，交通數(shù)據(jù)包括鐵路、高速路、國道、縣道；水系包括主要的河流、水庫、濕地和湖泊；居民點(diǎn)包括農(nóng)村居民點(diǎn)和城鎮(zhèn)居民點(diǎn)；行政區(qū)劃分到鄉(xiāng)級并標(biāo)有行政駐地的位置和注記等。

（2）選取1980 年（1980 年末期）、1990 年、1995 年、2000 年、2005 年、2010 年、2015 年 和2020 年共8 期流域的土地利用數(shù)據(jù)，其中前7 期數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心，研究小組以Landsat TM/ETM 和Landsat 8 為主要信息源解譯以上數(shù)據(jù)，分辨率為30 m。該流域2020年的土地利用數(shù)據(jù)是通過Landsat 8 影像數(shù)據(jù)人工解譯獲得的。研究小組將該流域的1980—2020 年的土地利用分類系統(tǒng)分為2 級，一是5 個(gè)一級土地利用類型，分別為耕地、林地、草地、水域及建設(shè)用地，二是14 個(gè)二級土地利用類型。

（3）漓江流域內(nèi)的高分辨率衛(wèi)星遙感影像數(shù)據(jù)、DEM 數(shù)據(jù)、社會經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)以及相關(guān)的規(guī)劃資料等，主要有《漓江流域生態(tài)保護(hù)條例》、桂林市 1 ∶5 萬地形圖、2009 年桂林市1 ∶1 萬土地利用現(xiàn)狀圖、2020 年桂林市行政區(qū)劃圖、《桂林市志（1991—2005 年）》《桂林市社會經(jīng)濟(jì)統(tǒng)計(jì)年鑒（1973—2020 年）》《桂林市國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計(jì)公報(bào)（2013 年）》《桂林市土地利用總體規(guī)劃（2006—2020 年）》《桂林旅游發(fā)展總體規(guī)劃（2001—2020 年）》《桂林市城市總體規(guī)劃（2010—2020 年）》以及《桂林市重點(diǎn)旅游片區(qū)規(guī)劃（2014—2020 年）》等。

1.3 數(shù)據(jù)預(yù)處理

為保證漓江流域邊界及河網(wǎng)提取準(zhǔn)確性，研究小組利用90 m、30 m 和12.5 m 多種分辨率DEM數(shù)據(jù)，加入流域高分辨率影像數(shù)據(jù)自行矢量化的主河道數(shù)據(jù)，然后利用ArcSWAT 軟件中Burn-in功能、Pre-defined 功能實(shí)現(xiàn)對原始DEM 及流域邊界等進(jìn)行修正后提取[3]。最后，提取流域范圍對所有地理信息數(shù)據(jù)進(jìn)行裁剪、重投影、重分類等預(yù)處理，將投影統(tǒng)一定義為CSGS2000-49N。

1.4 2020 年漓江流域土地利用數(shù)據(jù)解譯關(guān)鍵技術(shù)

1.4.1 影像數(shù)據(jù)預(yù)處理

利用 ENVI 軟件，以流域高程數(shù)據(jù)和歷年影像數(shù)據(jù)為地理參考，在流域四周選擇易識別地物處20 個(gè)作為地面控制點(diǎn)，通過多項(xiàng)式變換模型，對2020 年遙感影像進(jìn)行幾何精校正，采用回歸分析去除大氣輻射后，利用正射校正降低山體陰影對影像的影響。首先，對校正后的影像進(jìn)行波段計(jì)算，使用Band Math 工具分別計(jì)算獲得NDVI、NDWI、WET、NDBSI、LST、NDBI、LCI 等 遙 感指數(shù)圖層；其次，對流域DEM 數(shù)據(jù)進(jìn)行坡度、坡向等常用地形參數(shù)計(jì)算；再次利用Gram-schmidt Pan Sharpening（GS）工具融合以上所有圖層；最后利用提取的流域矢量邊界進(jìn)行裁剪，獲得漓江流域遙感影像數(shù)據(jù)。

1.4.2 遙感影像解譯

研究小組基于大數(shù)據(jù)思維，結(jié)合流域?qū)嶋H情況以及已有的成果數(shù)據(jù)，利用自主研發(fā)的隨機(jī)森林影像分類軟件對流域影像進(jìn)行土地利用類型識別與分類，包括訓(xùn)練樣本制作、模型訓(xùn)練及參數(shù)調(diào)試、地物分類提取等步驟。研究小組制作訓(xùn)練樣本時(shí)主要參考了2015 年土地利用類型數(shù)據(jù)及流域高分辨率影像，分別針對5 類土地利用類型制作了100 個(gè)樣本；然后利用該樣本數(shù)據(jù)進(jìn)行模型訓(xùn)練和參數(shù)調(diào)試；最后利用優(yōu)化后的訓(xùn)練模型及其參數(shù)對2020 年漓江流域影像進(jìn)行土地利用數(shù)據(jù)。

1.4.3 分類后處理

研究小組利用聚類分析（Clump）和去除分析（Eliminate）工具將＜4 個(gè)像元的破碎小圖斑進(jìn)行地類歸并，然后結(jié)合土地利用現(xiàn)狀圖、高分辨率遙感影像等資料對分類后的影像進(jìn)行目視人工修正和檢驗(yàn)，應(yīng)注意調(diào)整山體陰影和水體的錯(cuò)分，并利用已有水系、交通矢量數(shù)據(jù)對分類結(jié)果進(jìn)行修正，以提高解譯結(jié)果的精度。

1.4.4 精度評價(jià)

研究小組采用隨機(jī)抽樣方式，在影像分類結(jié)果中隨機(jī)抽取約 300 個(gè)隨機(jī)點(diǎn)，結(jié)合土地利用現(xiàn)狀圖、流域高分辨率遙感影像數(shù)據(jù)等資料，分別把隨機(jī)點(diǎn)與其實(shí)際類型進(jìn)行比對，輸入各隨機(jī)點(diǎn)的實(shí)際類別，生成分類精度報(bào)告，經(jīng)精度評估的總體分類精度為89.73%，達(dá)到研究需求。

2 研究方法

2.1 土地利用變化特征

研究小組主要選取單一土地利用動態(tài)度（K）和綜合土地利用動態(tài)度（LC）2 個(gè)指標(biāo)，表示土地利用變化的速率與程度，從類型水平和整體來定量評估土地利用變化動態(tài)特征。單一土地利用動態(tài)度（K）用于表示某種土地利用類型在具體研究時(shí)段的變化速率，綜合土地利用動態(tài)度（LC）則用于表征某一時(shí)段內(nèi)研究區(qū)土地利用類型轉(zhuǎn)化的劇烈程度，LC的數(shù)值越大，說明土地利用變化越劇烈，反之則較弱，計(jì)算公式如下[4]：

式（1）中，Ua和Ub分別表示研究期初、期末某種土地利用類型的面積；LU表示土地利用類型的年際動態(tài)變化率；△LUi-j為研究時(shí)段T內(nèi)i類型土地轉(zhuǎn)為j類型土地的面積；T為監(jiān)測時(shí)長。

2.2 土地利用轉(zhuǎn)移矩陣分析

研究小組基于GDAL 庫進(jìn)行二次開發(fā)，開發(fā)了兩期土地利用轉(zhuǎn)移矩陣快速計(jì)算工具，對漓江流域1980—2020 年以來的8 個(gè)時(shí)段的土地利用數(shù)據(jù)進(jìn)行轉(zhuǎn)移矩陣分析，之后統(tǒng)計(jì)各土地利用類型轉(zhuǎn)移的趨向、面積及變化率，分析各研究階段土地利用轉(zhuǎn)移特征。

2.3 土地利用變化驅(qū)動力分析

研究小組采用定量分析與定性描述相結(jié)合的方法，分析不同時(shí)段土地利用變化的驅(qū)動力。通過對比不同時(shí)期土地利用轉(zhuǎn)移特征、三次產(chǎn)業(yè)比例、人口數(shù)、城市人口數(shù)等社會經(jīng)濟(jì)指標(biāo)變化特征，查閱相應(yīng)時(shí)段漓江流域內(nèi)出臺的政策措施等，分析不同時(shí)段土地利用變化的主要驅(qū)動力。

3 流域土地利用變化情況分析

過去40 年中，伴隨著城鎮(zhèn)化進(jìn)程的加快、旅游活動的增加，特別是桂林國際旅游勝地建設(shè)的持續(xù)推進(jìn)，漓江流域的土地利用發(fā)生顯著變化。

3.1 流域土地利用總體變化

從土地利用類型構(gòu)成特征來看（見表1），1980—2020 年間漓江流域各土地利用類型均發(fā)生了不同程度的變化，到2020 年該流域的林地、耕地和草地面積減少，水域和建設(shè)用地面積增加，其中，耕地和建設(shè)用地面積變幅最大。林地一直為漓江流域的主要土地利用類型，占總面積的62%左右，主要分布在流域的北部貓兒山、越城嶺和東部的都龐嶺、海洋山等大型山脈，由1980 年前的3 969.10 km2減少到2020 年的3 937.00 km2，凈減32.1 km2。耕地作為流域第二大土地利用類型，主要分布在漓江河谷地勢平緩地帶，為城鎮(zhèn)化建設(shè)主要區(qū)域，因種植結(jié)構(gòu)調(diào)整以及建設(shè)用地的侵占致使面積持續(xù)減少，由1980年的1 520.71 km2減少到2020 年的 1 426.03 km2，凈減94.68 km2；草地面積減少了18.19 km2；水域面積增加了18.44 km2；建設(shè)用地占比不大，約3%～5%，主要分布在流域內(nèi)的桂林市區(qū)，興安、靈川、陽朔等縣城及周邊鄉(xiāng)鎮(zhèn)等，特別是臨桂新區(qū)和靈川縣等地的建設(shè)用地面積，呈快速增加之勢，面積由1980 年的 179.01 km2快速增加到 2020年的305.54 km2，凈增126.53 km2，增長幅度達(dá)到70.7%。這主要是由流域內(nèi)城鎮(zhèn)擴(kuò)張、工業(yè)用地以及旅游景區(qū)開發(fā)等因素驅(qū)動的。

表1 1980—2020 年漓江流域土地利用結(jié)構(gòu)特征表

3.2 流域土地利用動態(tài)變化分析

1980—2020 年間漓江流域各土地利用類型均發(fā)生不同程度變化，林地面積先減后增再減，幅度不大。耕地面積先減少再增加，之后急劇減少，2010—2015 年間減幅最大，凈減值為43.19 km2，約占1980—2020 年總減少量的45.6%，近年來桂林市持續(xù)進(jìn)行“旱改水”耕地提質(zhì)改造項(xiàng)目和“先補(bǔ)后占”“占優(yōu)補(bǔ)優(yōu)”的耕地占補(bǔ)平衡等措施使得耕地銳減趨勢得到有效控制。草地面積先增加后減少。水域面積先增加后減少，而后持續(xù)增加，且增幅較大。建設(shè)用地隨著靈川縣、臨桂新區(qū)建設(shè)和桂林市其他區(qū)、興安縣、陽朔縣、平樂縣等城鎮(zhèn)化建設(shè)的持續(xù)進(jìn)行而增加，2010—2015 年間增幅最大，5 年內(nèi)凈增69.7 km2，約占1980—2020 年總增加量的55.1%，增長趨勢較為明顯。

3.3 流域土地利用轉(zhuǎn)移特征分析

綜上所述，研究小組對1980—2020 年漓江流域的土地利用類型進(jìn)行分析，得到以下結(jié)論：

（1）漓江流域的耕地、林地、草地、水域和建設(shè)用地在1980—2020 年間的轉(zhuǎn)化率分別為：1.49、0.51、0.28、0.29、1.99，轉(zhuǎn)化率大小排列依次是建設(shè)用地＞耕地＞林地＞水域＞草地（見表1）。

（2）漓江流域林地主要轉(zhuǎn)化為至建設(shè)用地。早期主要轉(zhuǎn)化為耕地，主要由開荒種糧、種菜等引起；后期則主要轉(zhuǎn)化為建設(shè)用地，主要由流域內(nèi)城鎮(zhèn)化建設(shè)而導(dǎo)致的；而后，退耕還林和石漠化治理等措施使林地面積得到一定補(bǔ)充，減緩林地流出趨勢。

（3）耕地的變更范圍比較廣泛，與其他各類要素間均有轉(zhuǎn)換發(fā)生，40 年間轉(zhuǎn)出變化率最高，達(dá)到1.49%（見表1），主要轉(zhuǎn)出流向建設(shè)用地、林地和水域，耕地轉(zhuǎn)出建設(shè)用地過程也印證了流域內(nèi)2 次城鎮(zhèn)化建設(shè)高潮，即1980—1990 年改革開放后城市快速發(fā)展和2010—2015 年間桂林主城區(qū)發(fā)展及靈川縣、臨桂新區(qū)建設(shè)，同時(shí)也反映出流域內(nèi)土地利用結(jié)構(gòu)發(fā)生了改變，從農(nóng)業(yè)逐步調(diào)整為工業(yè)、旅游產(chǎn)業(yè)等。

（4）流域內(nèi)水域面積增加主要來源于耕地和林地。早期主要來源于耕地，是因?yàn)轸~類的養(yǎng)殖效益高，導(dǎo)致大量耕地被用于挖塘養(yǎng)魚；2000 年后水域面積逐年增加，也主要來源于耕地和林地，且增加趨勢較為明顯。2011 年后漓江流域持續(xù)進(jìn)行的生態(tài)環(huán)境修復(fù)，建成斧子口、小榕江、川江等多座漓江補(bǔ)水樞紐工程以及建設(shè)漓江蓄水壩等的建設(shè)，增加流域的水域面積。

（5）草地的變更較為單一，變化幅度也最小。1980—1995 年間草地面積增加主要來源于林地，而后2010—2020 年間因經(jīng)濟(jì)發(fā)展，大量農(nóng)民進(jìn)城務(wù)工，開荒種糧種果樹減少，且部分耕作條件差的農(nóng)田被撂荒變?yōu)椴莸兀圆莸孛娣e大量增加；草地面積減少主要是由于轉(zhuǎn)為林地、建設(shè)用地和耕地等，反映出當(dāng)?shù)赝烁€林和生態(tài)修復(fù)的效果，有效改善了流域內(nèi)生態(tài)環(huán)境，同時(shí)也保障了經(jīng)濟(jì)發(fā)展用地需求。

4 結(jié) 語

研究小組以漓江流域?yàn)檠芯繉ο螅?980—2020 年以來8 期土地利用數(shù)據(jù)，運(yùn)用“3S”技術(shù)和相關(guān)分析方法，綜合分析1980—2020 年間漓江流域土地利用時(shí)空變化及驅(qū)動力，得到以下結(jié)論：

（1）1980—2020 年漓江流域最主要土地利用類型為林地，多年平均占比為 62.20%，其次為耕地、草地、建設(shè)用地和水域，多年平均占比分別為23.29%、9.45%、3.58%、1.48%。林地分布廣泛，耕地和草地分布較為分散，建設(shè)用地主要沿河分布，集中在流域中部河谷盆地。

（2）1980—2020 年以來，漓江流域各土地利用類型整體變化幅度不大，主要表現(xiàn)為耕地、林地和草地面積減小，建設(shè)用地和水域面積增加。在此期間，經(jīng)濟(jì)進(jìn)一步發(fā)展加快城鎮(zhèn)化、工業(yè)化和旅游產(chǎn)業(yè)化步伐，使得社會經(jīng)濟(jì)用地大量侵占生態(tài)環(huán)境用地，社會經(jīng)濟(jì)用地驟增。

（3）1980—1990 年間流域內(nèi)各土地利用類型轉(zhuǎn)移的主要去向?yàn)榻ㄔO(shè)用地和草地，此階段由于改革開放發(fā)展，城鎮(zhèn)化建設(shè)迅速起步，建設(shè)用地大量增加；1990—2000 年間各土地利用類型轉(zhuǎn)移的主要去向?yàn)榱值?、耕地和水域，此階段流域內(nèi)的主要活動為農(nóng)業(yè)生產(chǎn)，各類型土地間的轉(zhuǎn)移主要受農(nóng)林種植結(jié)構(gòu)調(diào)整、挖塘養(yǎng)魚和占用草地用于農(nóng)林開發(fā)等因素驅(qū)動；2000—2010 年后，耕地、建設(shè)用地、林地和水域成為流入的主要土地利用類型，耕地和林地則是土地利用轉(zhuǎn)出的主要類型，造成這一時(shí)期土地利用變化的主要驅(qū)動力是城鎮(zhèn)化進(jìn)程加速、產(chǎn)業(yè)結(jié)構(gòu)調(diào)整、果業(yè)種植等因素；2010—2020 年間林地、耕地和建設(shè)用地是土地利用的主要流入類型，城鎮(zhèn)快速擴(kuò)張、工業(yè)化、旅游業(yè)快速發(fā)展、農(nóng)業(yè)種植結(jié)構(gòu)調(diào)整以及林業(yè)、生態(tài)建設(shè)是此階段土地利用變化的主要驅(qū)動因素。

研究小組分析了漓江流域 1980—2020 年間8期土地利用變化的時(shí)空特征以及驅(qū)動因素，受遙感影像數(shù)據(jù)和分辨率等因素限制，所提取的土地利用類型僅為一級大類，包括耕地、林地、草地、水域和建設(shè)用地，僅能宏觀分析流域的土地利用時(shí)空變化情況及趨勢，很難滿足流域內(nèi)土地利用精細(xì)化的監(jiān)管要求。利用“空—天—地”立體化自然資源調(diào)查監(jiān)測手段可提高流域內(nèi)影像數(shù)據(jù)獲取頻次和分辨率，實(shí)現(xiàn)流域自然資源利用精細(xì)化時(shí)空變化監(jiān)測及模擬預(yù)測。

綜上所述，研究小組為漓江流域的經(jīng)濟(jì)與生態(tài)環(huán)境保護(hù)的協(xié)調(diào)發(fā)展提出以下3 點(diǎn)建議：（1）嚴(yán)格落實(shí)耕地保護(hù)政策，持續(xù)進(jìn)行“旱改水”耕地提質(zhì)改造項(xiàng)目和“先補(bǔ)后占”“占優(yōu)補(bǔ)優(yōu)”的耕地占補(bǔ)平衡等措施，有效解決流域耕地減少問題；（2）合理確定規(guī)模布局，強(qiáng)化節(jié)約集約用地，嚴(yán)格控制建設(shè)用地中交通、商業(yè)、工業(yè)用地的規(guī)模，并按照總體規(guī)劃的布局，引導(dǎo)向控制協(xié)調(diào)區(qū)集中。（3）持續(xù)開展流域生態(tài)環(huán)境修復(fù)，增加流域林地及水域面積。特別是在流域源頭核心區(qū)域，要按照生態(tài)系統(tǒng)演替的規(guī)律和方向，結(jié)合生態(tài)環(huán)境保護(hù)政策，處理好開發(fā)與保護(hù)的關(guān)系、人與自然的關(guān)系，按自然資源保護(hù)規(guī)劃“分級保”“分區(qū)護(hù)”的要求，合理、科學(xué)地進(jìn)行退耕還林、退園還林，加強(qiáng)封山育林，合理利用現(xiàn)有草地規(guī)模，轉(zhuǎn)為灌木林地、喬木林地等林地，增大林地的范圍，提高森林覆蓋率。