李慧婷，徐麗萍，郭 鵬，劉 琳

（石河子大學(xué) 理學(xué)院，新疆 石河子832003）

土地利用/覆被變化（LUCC）研究是全球變化的重要研究課題［1－2］。土地利用/覆被格局具有典型的空間異質(zhì)性，在空間上表現(xiàn)為不同土地利用類型斑塊的鑲嵌，反映了土地生態(tài)過程的作用結(jié)果。景觀格局指數(shù)可準(zhǔn)確定量反映土地利用方式變化的時空狀態(tài)［3－4］。因此，很多國際組織在土地利用/覆被變化研究項目中都結(jié)合景觀格局分析來研究土地利用/覆被變化的空間結(jié)構(gòu)及變化趨勢。關(guān)于流域土地利用/覆被變化及景觀格局變化，國內(nèi)外學(xué)者已經(jīng)進(jìn)行了大量研究，目前，國內(nèi)學(xué)者已對滇池流域、白洋淀流域、小江流域、石羊河流域、太湖流域、瀾滄江流域等開展了流域土地利用景觀格局研究［5－13］。不過這些研究多數(shù)采用景觀格局指數(shù)直接分析或者模型分析的方法，雖然能反映景觀格局的幾何特征，但往往難以揭示景觀格局演化的深層規(guī)律，而且模型方法也存在不確定性。而地形作為影響景觀分布最重要的原因之一，有研究表明，地形是影響景觀格局最重要的自然因子，是土地利用規(guī)劃和景觀規(guī)劃研究中必須考慮的自然因子，同時，地形也是景觀演變方向和速度的重要制約因素［14－18］。但在瑪河流域這樣的典型干旱地區(qū)鮮有基于地形對景觀格局演化規(guī)律的分析，本研究以1990年、2000年、2005年、2010年四個時期的土地利用數(shù)據(jù)以及60m分辨率的DEM數(shù)據(jù)為信息源，借助GIS和Fragstats軟件，通過利用空間數(shù)據(jù)疊加、景觀轉(zhuǎn)移矩陣揭示各景觀類型的時空轉(zhuǎn)化動態(tài)，并基于DEM數(shù)據(jù)的高程數(shù)據(jù)、坡度數(shù)據(jù)和坡向數(shù)據(jù)進(jìn)行空間疊加和數(shù)理統(tǒng)計，探索瑪河流域景觀格局變化的三維特征及影響因素，為實現(xiàn)流域土地資源合理開發(fā)利用、景觀格局的優(yōu)化及可持續(xù)發(fā)展提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)簡介

瑪納斯河流域地處新疆天山北麓準(zhǔn)噶爾盆地南緣，位于天山北坡經(jīng)濟帶的中心。由南向北依次為山地、山前平原、沙漠三大地貌單元（比例約為2.83∶1∶1.07）［19］。從東到西依次為塔西河、瑪納斯河、寧家河、金溝河、八音溝河及其相連的5片沖積洪積扇。各河流均發(fā)源于北天山中段，順天山北坡向北流入干旱的準(zhǔn)噶爾盆地，其中瑪納斯河全長400km，是準(zhǔn)噶爾盆地水量最大、流程最長的內(nèi)陸河流，注入準(zhǔn)噶爾盆地西北部的瑪納斯湖（已干枯）。流域行政上屬于瑪納斯縣、沙灣縣、石河子墾區(qū)和克拉瑪依部分地區(qū)。年平均降水量100～200mm，時空分布不均。年平均蒸發(fā)量1 500～2 100mm，年平均氣溫4.7～5.7℃，屬于典型的溫帶大陸性干旱氣候。

瑪納斯河流域在古代屬于純牧區(qū)，自清代以來成為新疆的重點開發(fā)區(qū)域之一，新中國成立以后，瑪納斯河流域成為新疆開墾的最大的人工綠洲和我國第四大灌溉農(nóng)業(yè)區(qū)［19］，社會經(jīng)濟取得巨大發(fā)展。但是近年來，由于自然環(huán)境的變化，人口的不斷增長，經(jīng)濟的發(fā)展，部分地區(qū)對水土資源的不合理利用等原因，使得流域土地利用發(fā)生劇烈的變化，出現(xiàn)了諸如尾閭湖干涸、草場退化、土地沙化、鹽漬化和生物多樣性降低等問題［19－20］。

1.2 數(shù)據(jù)獲取與處理

采用1990年、2000年、2005年與2010年的Landsat TM（分辨率30m×30m）遙感影像數(shù)據(jù)，輔助資料為瑪納斯河流域行政區(qū)劃圖、土地利用圖、Gooleearth等。數(shù)據(jù)處理平臺為Envi 4.8，ArcGIS 10.0和Fragstats 3.4。首先利用Envi 4.8對得到的遙感影像進(jìn)行預(yù)處理，之后利用監(jiān)督分類法結(jié)合DEM，Gooleearth等輔助對圖像進(jìn)行分類?；趧⒓o(jì)遠(yuǎn)等土地利用分類標(biāo)準(zhǔn)并結(jié)合研究區(qū)自身特點，將原標(biāo)準(zhǔn)進(jìn)行修訂最終將二級地類合并為耕地、林地、草地、水域、鹽堿地、荒漠、冰川/永久積雪和城鄉(xiāng)建設(shè)用地八個一級地類，并將柵格數(shù)據(jù)分辨率重采樣為60m，最后在Fragstats 3.4中計算景觀格局指數(shù)。

1.3 景觀格局指標(biāo)選擇

本研究根據(jù)數(shù)據(jù)的易取性和計算精度兩方面選取斑塊形狀指數(shù)，包括斑塊邊界長度（TE）、邊界密度（ED）、面積加權(quán)平均形狀指數(shù)（AWMS1）、面積加權(quán)平均斑塊分維數(shù)（AWMPFD）；景觀面積（TA）、景觀類型面積比（PLAND）等進(jìn)行景觀格局分析。以上這些指數(shù)的計算方法以及生態(tài)學(xué)意義詳見參考文獻(xiàn)［4，21］，本文不再一一詳述。

2 結(jié)果與分析

2.1 瑪河流域景觀格局現(xiàn)狀及時間變化

瑪納斯河流域2010年景觀類型總面積為3 404 998 hm2，把荒漠、鹽堿地等未利用地作為主導(dǎo)景觀類型，草地和耕地其次，草地和荒漠兩者占景觀總面積的百分比分別為35.8%和31.8%，控制著該地區(qū)的景觀；其次為耕地，面積百分比為21.85%，鹽堿地、水域、林地、冰川/永久積雪以及城鄉(xiāng)建設(shè)用地五種景觀類型所占面積較少，五類之和為11.6%。從面積比重看，草地應(yīng)是研究區(qū)景觀的基質(zhì)類型（見表1）。

1990—2010年20a間，耕地以及城鄉(xiāng)建設(shè)用地的增加幅度最大，耕地總面積共增加了283 761hm2，增長率為61.66%；城鄉(xiāng)建設(shè)用地雖然增加的面積不大，但增長率卻非常高，達(dá)到39.1%。草地以及冰川/永久積雪減少的幅度最大，草地由1990—2010年總面積減少了184 045hm2，冰川/永久積雪的減少率更是高達(dá)29.85%。另外，荒漠、鹽堿地也有不同程度的減少，林地與水域的變化幅度不大，總之，耕地和城鄉(xiāng)建設(shè)用地面積驟增的同時伴隨著草地、荒漠等其他地類面積的減少（見圖1）。

2.2 瑪河流域景觀斑塊形狀的時間變化

描述景觀斑塊形狀的指數(shù)有很多，包括斑塊邊界總長度TE、邊界密度ED、面積加權(quán)平均形狀指數(shù)AWMSI以及面積加權(quán)平均斑塊分維數(shù)AWMPFD。

表1 瑪納斯河流域景觀單元類型面積及百分比

圖1 瑪納斯河流域景觀單元類型面積百分比

表2 瑪納斯河流域景觀斑塊形狀指數(shù)

就整個瑪納斯河流域而言，1990—2010年，研究區(qū)的斑塊邊界總長度TE、邊界密度ED、面積加權(quán)平均形狀指數(shù)AWMSI以及面積加權(quán)平均斑塊分維數(shù)AWMPFD總體都呈下降趨勢，尤其是斑塊邊界總長度下降趨勢最為明顯（見表2），這說明瑪納斯河流域斑塊形狀總體趨向規(guī)則與簡單的發(fā)展趨勢。這可能是由近年來流域的不斷開發(fā)所引起的，隨著人類社會與科學(xué)技術(shù)的快速發(fā)展，大量的天然用地被開發(fā)為人工用地，使得原來小面積的人工用地發(fā)展成大面積的、連片集中的景觀類型。

就瑪納斯河流域各景觀類型來看，耕地的邊界密度ED、面積加權(quán)平均形狀指數(shù)AWMSI以及面積加權(quán)平均斑塊分維數(shù)AWMPFD在各個時期都是最大的，其次為草地、荒漠，林地、水域和城鄉(xiāng)建設(shè)用地的各項指數(shù)均較?。ㄒ姳?）。這說明在研究期內(nèi)，耕地、草地與荒漠的斑塊形狀指數(shù)較為不規(guī)則和復(fù)雜，這可能是由于草地與荒漠在研究區(qū)所占面積較大，與其他景觀類型鑲嵌分布，因此其形狀不規(guī)則，而耕地則是由于分布較分散，且單個斑塊面積較小，造成斑塊體形狀復(fù)雜。林地、水域和城鄉(xiāng)建設(shè)用地的斑塊形狀最為簡單，這主要是由于這幾種景觀類型斑塊形狀都比較穩(wěn)定與規(guī)則。

表3 瑪納斯河流域景觀類型斑塊形狀指數(shù)

1990—2010年，不同景觀類型的面積加權(quán)平均形狀指數(shù)AWMSI以及面積加權(quán)平均斑塊分維數(shù)AWMPFD變化趨勢是不同的，荒漠、草地、林地與冰川＼永久積雪整體呈下降趨勢（見表3）。城鄉(xiāng)建設(shè)用地與耕地有明顯的增加趨勢，這兩種景觀類型受人類活動影響最顯著，尤其是近十年來，面積不斷增加，在擴張的過程中導(dǎo)致其斑塊形狀趨于復(fù)雜。鹽堿地、水域變化趨勢不明顯。

2.3 景觀單元類型轉(zhuǎn)換的地形梯度特征

通過對研究區(qū)1990年與2000年的土地利用數(shù)據(jù)進(jìn)行疊加分析，得到景觀演變面積最大的十種轉(zhuǎn)化類型，分別為：草地轉(zhuǎn)為耕地、荒漠轉(zhuǎn)為耕地、冰川/永久積雪轉(zhuǎn)為荒漠、冰川/永久積雪轉(zhuǎn)為草地、荒漠轉(zhuǎn)為草地、鹽堿地轉(zhuǎn)為荒漠、耕地轉(zhuǎn)為城鄉(xiāng)建設(shè)用地、耕地轉(zhuǎn)為草地、草地轉(zhuǎn)為城鄉(xiāng)建設(shè)用地、水域轉(zhuǎn)為草地。

基于研究區(qū)數(shù)字高程模型，利用ArcMap軟件獲取高程圖、坡度圖與坡向圖，并將其分級，然后與主要景觀類型演變數(shù)據(jù)進(jìn)行疊加分析，統(tǒng)計得出各景觀演變類型在不同高程、坡向、坡度的分布情況。

根據(jù)主要景觀單元演變類型分級統(tǒng)計結(jié)果，研究區(qū)景觀類型演變存在明顯的地域分異規(guī)律，在低海拔地區(qū)以草地向耕地的演變?yōu)橹饕葑兎绞?，尤其是在海拔?50～450m的地區(qū)，草地轉(zhuǎn)為耕地的比例高達(dá)86%；在0～350m的地區(qū)，主要的景觀演變類型包括荒漠轉(zhuǎn)為耕地、草地轉(zhuǎn)為耕地以及鹽堿地轉(zhuǎn)為荒漠；3 000～5 000m的區(qū)域，以冰川/永久積雪轉(zhuǎn)為荒漠以及冰川/永久積雪轉(zhuǎn)為草地為主（見表4）。這說明近20a來，隨著人類活動的加劇，在人類的活動范圍內(nèi)，即低海拔地區(qū)，大量的天然用地演變?yōu)槿斯び玫?，包括耕地以及城鄉(xiāng)建設(shè)用地。在高海拔地區(qū)主要表現(xiàn)為冰川/永久積雪轉(zhuǎn)為荒漠以及冰川/永久積雪轉(zhuǎn)為草地兩種景觀單元演變類型，這與高海拔地區(qū)冰川/永久積雪的大面積分布有關(guān)。

表4 主要景觀單元演變類型高程分級統(tǒng)計

景觀演變類型與坡向的關(guān)系不太密切，每一種景觀演變類型在各坡向的分布差異不大，只有草地轉(zhuǎn)為耕地類型在北坡的分布面積明顯高于南坡，冰川/永久積雪轉(zhuǎn)為荒漠以及冰川/永久積雪轉(zhuǎn)為草地兩種景觀單元演變類型主要分布在南坡，前者主要是由瑪納斯河流域的南高北低的地形特征決定的，后者主要是因為南坡溫度高于北坡，造成南坡冰雪消融。整體看來，草地轉(zhuǎn)為耕地在各個坡向的分布均較廣，其次是荒漠轉(zhuǎn)為耕地，水域轉(zhuǎn)為草地的區(qū)域在各個坡向均最少，此分布規(guī)律與景觀演變類型的面積排序一致，說明坡向?qū)坝^演變類型的分布影響不大（見圖2）。

坡度對景觀類型演變的影響比較明顯，草地向耕地的演變主要位于坡度介于0°～10°的區(qū)域，荒漠向耕地的演變主要集中在0°～3°的區(qū)域，這都與耕地的限制因素相關(guān)；冰川/永久積雪向其他景觀類型的演變主要位于大于10°的坡度較高的地區(qū)，在0°～3°的區(qū)域沒有分布；其他幾種主要的景觀演變類型多分布在坡度小于10°的區(qū)域，尤其是鹽堿地向荒漠的演變以及耕地向城鄉(xiāng)建設(shè)用地的演變（見表5）。此外，值得關(guān)注的是：在10°～20°的較陡坡區(qū)域，仍有大面積的草地轉(zhuǎn)化為耕地，甚至在40°～90°的陡坡也有部分草地轉(zhuǎn)為耕地，這對生態(tài)環(huán)境的保護是極其不利的。

圖2 主要景觀單元演變類型坡向分級統(tǒng)計

總體分析主要景觀單元演變類型在各高程梯度、坡向梯度以及坡度梯度分布的數(shù)據(jù)發(fā)現(xiàn)：景觀類型演變存在明顯的地域分異規(guī)律，尤其與高程和坡度因子密切相關(guān)。

表5 主要景觀單元演變類型坡度分級統(tǒng)計

3 結(jié)論與建議

瑪納斯河流域1990—2010年土地利用類型變化顯著，荒漠、草地、水域以及林地面積都出現(xiàn)了不同程度的減少，耕地以及城鄉(xiāng)建設(shè)用地面積不斷增加，各景觀類型面積差異不斷減少。1990—2010年，瑪納斯河流域斑塊形狀總體趨向規(guī)則與簡單的發(fā)展趨勢，城鄉(xiāng)建設(shè)用地與耕地斑塊形狀趨于復(fù)雜，尤其是近十年來，面積不斷增加，在擴張的過程中導(dǎo)致其斑塊形狀趨于復(fù)雜。鹽堿地、水域變化趨勢不明顯。景觀類型演變存在明顯的地域分異規(guī)律，尤其與高程和坡度因子密切相關(guān)。

針對以上研究結(jié)論對土地利用規(guī)劃及景觀規(guī)劃提出以下幾個方面的建議：（1）控制耕地與城鄉(xiāng)建設(shè)用地的面積，尤其是耕地的面積，響應(yīng)國家退耕還林與退耕還草的政策。為保證糧食產(chǎn)量，應(yīng)該在保護生態(tài)環(huán)境的基礎(chǔ)上，提高耕地質(zhì)量。與此同時，應(yīng)采用一定的措施防治土壤鹽漬化與土地沙漠化。（2）針對城鄉(xiāng)建設(shè)用地與耕地斑塊形狀趨于復(fù)雜的問題，應(yīng)對城鄉(xiāng)建設(shè)用地進(jìn)行合理規(guī)劃，優(yōu)化城鄉(xiāng)景觀格局，提高城市生產(chǎn)效率。（3）根據(jù)不同地形梯度上景觀格局演變類型分布結(jié)果，在低海拔、坡度小的區(qū)域，是人類主要聚集地，目前也是環(huán)境污染最嚴(yán)重的區(qū)域，因此在發(fā)展經(jīng)濟的同時，要嚴(yán)格控制并管理重污染型工廠、防治城市點源污染與耕地面源污染；在中低海拔、坡度較緩的區(qū)域，開發(fā)耕地會加大水土流失發(fā)生的概率，但種植果樹有助于防風(fēng)固沙，并且能增加經(jīng)濟效益；在海拔較高、坡度較大的地區(qū)，更要減少耕地的分布，把已有的耕地轉(zhuǎn)為林地、草地，保持生態(tài)環(huán)境的平衡；在坡度陡的高海拔區(qū)域，植被分布主要是高山草原、草甸，景觀類型主要是冰川/永久積雪以及裸地，因此該區(qū)域生態(tài)環(huán)境及其脆弱，根據(jù)生態(tài)系統(tǒng)的穩(wěn)定性原理，此種生態(tài)環(huán)境一旦遭到破壞便很難恢復(fù)，所以更應(yīng)該加強生態(tài)環(huán)境的保護，具體做到規(guī)范并控制旅游開發(fā)，嚴(yán)禁破壞。

4 討 論

瑪納斯河流域作為新疆開墾的最大的人工綠洲和我國第四大灌溉農(nóng)業(yè)區(qū)［19］，社會經(jīng)濟取得巨大發(fā)展。近年來，由于自然環(huán)境的變化，人口的不斷增長，經(jīng)濟的發(fā)展，部分地區(qū)對水土資源的不合理利用等原因，瑪納斯河流域的景觀格局發(fā)生了變化。其中，人口增加是瑪納斯河流域景觀格局變化的主要驅(qū)動因素，進(jìn)入21世紀(jì)以來，人口的增加導(dǎo)致研究區(qū)耕地與城鄉(xiāng)建設(shè)用地的面積迅速擴大，與此同時，草地與荒漠的面積逐漸減少，進(jìn)而造成各景觀類型面積所占比例差異減少，這一點，與以往研究得出的瑪納斯河流域景觀演變類型主要以草地與未利用土地轉(zhuǎn)化為耕地與建設(shè)用地的結(jié)論一致［22］。研究中發(fā)現(xiàn)，隨著人口的增加，對耕地的需求量不斷上升，在坡度較陡的區(qū)域，仍有大面積的草地轉(zhuǎn)化為耕地，這對生態(tài)環(huán)境的保護是極其不利的，會大大加大水土流失的發(fā)生頻率。要實現(xiàn)區(qū)域景觀格局的優(yōu)化和生態(tài)功能的良性循環(huán)，必須重視土地利用格局優(yōu)化，維護景觀生態(tài)過程與格局的連續(xù)性。

本文采用景觀單元類型轉(zhuǎn)移矩陣、主要景觀演化類型的地形梯度分析和空間數(shù)據(jù)疊加等方法，分別從高程、坡向、坡度三個方面對景觀格局演變的分布規(guī)律進(jìn)行分析，結(jié)果表明：景觀類型演變存在明顯的地域分異規(guī)律，尤其與高程和坡度因子密切相關(guān)。但不同高程、坡向以及坡度帶上景觀演變類型對生態(tài)環(huán)境的影響及效應(yīng)還未明確，在今后的研究中會對此問題進(jìn)行進(jìn)一步探討。

［1］ Turner Ii B L，Skole D，Sanderson S，et al.Land－use and land－cover change.Science/Research plan［R］.IHDP Report，1995.

［2］ 馮仕超，高小紅，亢健，等.西寧市30多年來土地利用/土地覆被變化及城市擴展研究［J］.干旱區(qū)研究，2012，29（1）：129－136.

［3］ 傅伯杰，邱揚，陳利頂.景觀生態(tài)學(xué)的原理及應(yīng)用［M］.北京：科學(xué)出版社，2001.

［4］ 鄔建國.景觀生態(tài)學(xué)：格局、過程、尺度與等級［M］.北京：高等教育出版社，2000.

［5］ 程剛，張祖陸，呂建樹.基于CA－Markov模型的三川流域景觀格局分析及動態(tài)預(yù)測［J］.生態(tài)學(xué)雜志，2013，32（4）：999－1005.

［6］ 劉鐵冬.雜谷腦河流域近30年景觀格局時空演變研究［J］.水土保持研究，2014，21（3）：98－102.

［7］ 宮繼萍，石培基，楊雪梅，等.干旱區(qū)內(nèi)陸河流域景觀格局及其與環(huán)境因子的關(guān)系研究①［J］.土壤，2013，45（1）：143－151.

［8］ 張飛，丁建麗，王伯超，等.干旱區(qū)綠洲土地利用/覆被及景觀格局變化特征：以新疆精河縣為例［J］.生態(tài)學(xué)報，2009，29（3）：1251－1263.

［9］ 魏偉，石培基，周俊菊，等.近20多年來石羊河流域景觀格局演變特征［J］.干旱區(qū)資源與環(huán)境，2013，27（2）：156－161.

［10］ 張國坤，鄧偉，呂憲國，等.新開河流域濕地景觀格局動態(tài)變化過程研究［J］.自然資源學(xué)報，2007，22（2）：204－210.

［11］ 么欣欣，韓春蘭，劉洪彬，等.基于RS與GIS的沈陽市土地利用及景觀格局變化［J］.水土保持研究，2014，21（2）：158－161.

［12］ 張明陽，王克林，何萍，等.白洋淀流域景觀空間格局變化研究［J］.資源科學(xué)，2005，27（2）：134－140.

［13］ Tischendorf L.Can landscape indices predict ecological processes consistently ［J］.Landscape Ecology，2001，16（3）：235－254.

［14］ 李丹，劉丹丹，趙金祥.基于DEM的山區(qū)土地利用變化分析［J］.水土保持研究，2014，21（1）：66－70.

［15］ 陳丹，周啟剛，何昌華，等.重慶山地都市區(qū)1985—2010年土地利用變化地形特征分異研究［J］.水土保持研究，2013，20（5）：210－220.

［16］ 高彥凈，鞏杰，賈珍珍，等.甘肅白龍江流域土地利用在地形梯度上的空間分布［J］.蘭州大學(xué)學(xué)報：自然科學(xué)版，2014，50（5）：680－686.

［17］ 衛(wèi)國，文倩，劉艷艷，等.基于DEM的醴陵市土地利用空間格局分析［J］.資源科學(xué)，2008，30（2）：228－234.

［18］ 卜心國，王仰麟，沈春竹，等.深圳市地形對土地利用動態(tài)的影響［J］.地理研究，2009，28（4）：1011－1021.

［19］ 賴先齊.綠洲鹽漬化棄耕地生態(tài)重建研究［M］.北京：中國農(nóng)業(yè)出版社，2007.

［20］ 張宏鋒，歐陽志云，鄭華，等.新疆瑪納斯河流域景觀格局變化及其生態(tài)效應(yīng)［J］.應(yīng)用生態(tài)學(xué)報.2009.20（6）：1408－1414.

［21］ McGarigal K，Cushman S A，Neel M C，et al.FRAGSTATS：spatial pattern analysis program for categorical maps［EB/OL］.www.umass.edu/landeco/research/fragstats/fragstats.html，2002.

［22］ 姜亮亮，劉海隆，包安明.瑪納斯河流域景觀格局演變特征與驅(qū)動機制分析［J］.水土保持研究，2014，21（4）：256－262.