桐柏大別山區土地利用時空變化圖譜分析

張榮華，劉霞，李歡，姚孝友，王友勝，張光燦

(1.山東農業大學林學院，山東省土壤侵蝕與生態修復重點實驗室，國家林業局泰山森林生態站，271018，山東泰安;2.北京師范大學地理學與遙感科學學院，100875，北京;3.水利部淮河水利委員會水土保持處，233001，安徽蚌埠)

桐柏大別山區土地利用時空變化圖譜分析

張榮華1，2，劉霞1?，李歡1，姚孝友3，王友勝1，張光燦1

(1.山東農業大學林學院，山東省土壤侵蝕與生態修復重點實驗室，國家林業局泰山森林生態站，271018，山東泰安;2.北京師范大學地理學與遙感科學學院，100875，北京;3.水利部淮河水利委員會水土保持處，233001，安徽蚌埠)

運用地學信息圖譜理論方法，基于地理信息系統時空分析模型，以桐柏大別山區20 a土地利用3期(1987、1997、2007)空間數據為基礎，通過建立桐柏大別山區土地利用變化圖譜和漲落勢圖譜，分析該區域土地利用變化的時空特征，預測未來土地利用變化的趨勢。結果表明：1)1987—2007年，研究區土地利用發生了很大變化，耕地面積持續性減少，林地和建設用地面積持續性增加，草地和未利用地面積先增加后減少，水域面積先減少后增加;2)空間格局上，耕地與林地相互轉化遍布整個研究區域，耕地和林地轉為建設用地主要分布在各縣市周邊地區，水域轉為耕地和林地主要位于水域與其他用地交接地帶;3)漲落勢圖譜中，林地漲勢最明顯，增加了5 904.06 km2，占新增圖譜單元61.31%，主要來源于耕地和水域轉出，耕地落勢最明顯，減少了6 394.71 km2，占萎縮圖譜單元66.41%，主要轉為林地和建設用地;4)未來10年，桐柏大別山區土地利用類型仍將以耕地、林地為主。

地學信息圖譜;時空變化;土地利用;桐柏大別山區

土地利用/土地覆被記錄了人類改變地球表面特征空間格局的活動，再現了地球表面景觀的時空動態過程，它既受自然因素的影響，又受社會、經濟、政治等人文因素的制約。進入20世紀90年代以來，土地利用/土地覆被變化(LUCC)受到越來越多的關注，成為全球環境變化研究的核心領域［1-3］。近年來，國內外學者探索將地學信息圖譜理論和方法應用于土地利用/覆被變化研究，以解決LUCC研究中時空復合分析的難題［4-9］。地學信息圖譜是以地球空間信息認知理論為基礎，以遙感、地理信息系統、虛擬現實、計算機制圖技術為支撐發展起來的一種時空復合分析方法。其圖譜單元既包含地理實體和地理現象的空間差異性信息，又包含地理過程的時序變化信息，將“表示空間單元特征的圖”與“表示時間發展之起點與過程的譜”合為一體，解決了空間與過程集成研究的問題［10-11］。筆者運用地學信息圖譜方法，通過對桐柏大別山區1987—2007以來的土地利用“空間格局”和“變化過程”進行定量化綜合研究，分析土地利用動態變化特征，預測未來土地利用變化的趨勢，以期為區域土地資源可持續利用和生態環境保護提供決策支持。

1 研究區概況

桐柏大別山區地處我國中部，屬淮河流域，E 113°16'～116°45'，N 30°57'～32°43'，跨皖、豫、鄂3省交界處，涉及16個縣(市、區)，總面積約2.72萬km2(圖1)。屬北亞熱帶過渡性氣候，年平均氣溫13.2～15.7℃，無霜期180～220 d，日照時間1 990～2 650 h，作物生長周期長，是我國重要的糧、油生產基地，在保證我國糧食安全和生態安全，促進國民經濟建設及可持續發展中具有極為重要的地位。境內地形復雜，地貌類型以山地、丘陵為主，占總面積的59%，主要分布在南部，母質大部分以花崗片麻巖為主，易風化，抗蝕性差，土壤以黃棕壤、水稻土和粗骨土為主。地帶性植被為落葉闊葉混交林型，年平均降水量1 200～1 500mm，6—8月超過全年60%，年際相差1～3倍，時空分布不均，易形成洪澇與干旱災害。

圖1 桐柏大別山區地理位置圖Fig．1 Geographical position of Tongbai-Dabiemountainous area

2 研究方法

2.1 數據來源與處理

以桐柏大別山區1987、1997、2007年4—5月份TM影像為數據源，空間分辨率為30m×30m。在遙感圖像處理軟件ERDAS 9.1中進行配準、鑲嵌、勻色及增強處理。利用1∶10萬地形圖選擇35個同名地面控制點(GCP)對影像進行校正，誤差控制在0.5個像元以內，校正中，采用北京1954坐標系，高斯-克呂格投影。

2.2 土地利用數據與圖譜合成

根據研究區土地利用現狀和土地資源特點，結合GB/T 21010—2007《土地利用現狀分類》，將研究區土地利用類型劃分為耕地、林地、草地、建設用地、水域和未利用地6類。針對遙感影像特征，結合野外實地調查，建立解譯判讀標志，采用監督分類和目視解譯相結合的方法對影像進行分類，提取土地利用信息，建立土地利用數據庫，并對分類結果運用Kappa指數進行精度評價，精度達88.36%。

將ERDAS解譯獲得的“.img”格式的土地利用數據，轉為“.grid”格式(柵格大小設為30m)，作為基本圖譜單元。以基本圖譜單元為基礎，采用Arc-GIS 9.2軟件中的空間分析模塊，利用地圖代數運算方法，對時間相鄰的土地利用數據進行合成，將不同采樣時刻上的地理單元的空間—屬性—過程一體化數據進行匹配、融合，建立土地利用變化圖譜。

2.3 土地利用漲落勢圖譜構建

土地利用變化包括2方面：一是轉入，即由其他用地類型轉移到本土地類型的部分;二是轉出，即原土地類型除不變區外，轉移到其他土地類型的部分。為了更直接把握土地利用變化的走向和趨勢，通過設定分類原則、重影射方法，弱化冗余信息、消除噪聲，將圖譜單元重新分類、綜合、組合、提取有效信息、重構圖譜，建立由轉入圖譜單元構建的漲勢系列圖譜，由轉出圖譜單元構建落勢系列圖譜。

2.4 土地利用動態度分析

土地利用動態度反映一定時間和范圍內某種土地利用類型的數量變化情況［12］，其數學表達式為

式中：Ri為研究時段內第i類土地利用類型面積的年變化，即土地利用動態度;Uai、Ubi分別為研究初期和末期第i類土地利用類型的面積，km2;T為研究時段，a。

2.5 馬爾柯夫預測模型

馬爾柯夫過程［5］是一個特殊的隨機過程，最大特點是t+1時刻的狀態只與t時刻的狀態有關，而與以前的狀態無關，這一特點決定了其適用于土地利用變化預測。關鍵在于確定各土地利用類型轉移概率矩陣P(式2)。根據柯爾可夫過程的性質和條件概率的定義，運用馬爾柯夫過程的基本方程(式3)，可模擬出任意一年的土地利用類型所占比例的變化情況。

式中：為第n期馬爾柯夫過程土地利用類型i轉化為土地利用類型j的轉移概率，轉移矩陣的每一項元素需滿足以下條件：0≤Pij≤1(i，j=1，2，…，N)，即各元素為非負值，且，即每行元素之和為1;k為中間變量;N為轉移概率矩陣的行列數。

3 結果與分析

3.1 土地利用時間序列變化圖譜分析

隨著人們對土地資源的開發利用，20年桐柏大別山區土地利用結構發生了很大變化(表1)。在面積減少的地類中，耕地減少最大，達4 175.22 km2，水域減少59.66 km2;同時，林地、草地、建設用地和未利用地不同程度增加，其中林地和建設用地增加最為顯著，分別增加了3 351.50和659.64 km2。從土地利用動態度來看，變化較大的是建設用地、耕地和林地，動態度分別為 +4.42%、-1.68%和+1.38%。對比2個時段(1987—1997年和1997—2007年)，耕地持續減少(后一時段動態度稍小)，林地持續增加(前一時段動態度較大)，這與該區20世紀90年代初實施林業二次創業、大力發展經濟林政策，90年代末國家頒布保護耕地資源的土地管理法規、退耕還林政策有關。建設用地持續增加，且后一時段增長較快，這與社會經濟發展、城市化、工業化進程加快密切相關。未利用地動態度最大，其原因可能為前一時段局部大規模進行坡地經濟林開發、毀林開荒、撂荒、棄耕等，導致土壤侵蝕、土地沙化加劇，進而導致未利用地增加，后一時段由于技術水平的提高和政府投入的增加，為未利用地的開墾提供了可能，土地開發利用率勢必提高。發展經濟林、退耕還林(草)等政策導向，使得研究區耕地面積大幅減少，林地和草地面積增加。

表1 1987—2007年桐柏大別山區土地利用變化情況Tab．1 Land use changes of Tongbai-Dabiemountainous area during 1987—2007

3.2 土地利用空間序列變化圖譜分析

桐柏大別山區土地利用變化圖譜共計36類圖譜單元(圖2)。從表2矩陣可以看出，1987—1997年，73.60%的圖譜單元土地利用方式沒有發生變化，而變化的圖譜單元中耕地和林地轉移最明顯，水域和建設用地次之，草地和未利用地轉移較少。耕地是受人類活動影響最大的土地利用類型，其變化特征與發展動態已成為土地利用研究的重要內容。20世紀90年代農業結構調整、發展經濟林、比較利益驅動及城市發展對建設用地的需求，使得大量耕地被占用，其面積大幅減少。耕地中分別有3 889.27和283.75 km2轉為林地和建設用地，占耕地轉出面積85.79%、6.26%，主要分布在光山縣、潢川縣、固始縣北部，信陽市轄區北部，桐柏縣東北部，霍山縣東部，六安市轄區西北部。同時，由于該區作為我國重要的糧食生產基地，為保證一定的糧食產量和質量，加之人們比較利益的驅動，耕地中1 898.58 km2因其他土地利用類型轉入而增加，其增加來源主要為林地、水域，分別占轉入面積的79.63%、12.95%，主要分布在金寨縣北部、霍山縣西部、新縣南部及西部、羅山縣中部、桐柏縣西部。

表2 1987—2007年土地利用動態變化轉移矩陣Tab．2 Transitionmatrix of land use changes during 1987—2007 km2

1997—2007年土地利用方式沒有發生變化的圖譜單元占66.93%，低于前一時段，土地利用發生了較大規模的變化，但變化最大的仍是耕地與林地之間的轉化。由于土地管理政策法規的頒布實施、城鎮化水平和耕作技術水平的提高、退耕還林(草)的實行，使得耕地減少的趨勢得到了有效的遏制，減少了1 540.08 km2，主要轉為林地、建設用地，集中在桐柏縣中東部、信陽市轄區北部、羅山縣中北部、光山縣中部、潢川縣東北部、固始縣中南部、金寨縣東南部、霍山縣西北部、六安市轄區西北部。同時，當地政府面對人口增長、糧食產量、國內生產總值(GDP)、財政收入的壓力，以及人們比較利益的驅動，使得耕地中3 310.72 km2由其他土地利用類型轉變增加，集中在桐柏縣東北部、信陽市轄區西北部、羅山縣西部、光山縣南部、新縣南部、固始縣西北部、金寨縣中部、霍山縣西北部及東北部、六安市轄區東北部。耕地與林地之間的頻繁轉換，說明退耕還林與毀林開荒的現象同時存在，也說明土地利用結構不合理，人地關系緊張。

由表3排序表［7］可以看出：20 a耕地向林地轉化5 523.37 km2，林地向耕地轉化1 871.29 km2，二者之差即耕地向林地凈轉化3 652.08 km2，占總變化面積37.94%。耕地和林地轉為建設用地分別占5.21%和2.68%，轉為水域分別占3.01%和1.76%。水域轉為耕地和林地分別占2.60%和2.23%?？梢钥闯?，桐柏大別山區土地利用結構變化表現為3個特征：一是耕地與林地相互轉化是全局性的，遍及整個區域，相對而言，霍山、金寨、固始、潢川、光山、信陽市轄區、桐柏等更為明顯;二是耕地和林地轉為建設用地明顯，主要在各縣市區周邊以及交通線兩側，是經濟相對發達的地區，居民用地和工礦、交通用地快速擴展，大量蠶食耕地和林地;三是水域轉為耕地和林地較為明顯，大部分位于沿河或庫區周圍一些水域邊緣地帶，而其余土地利用類型的空間移動則成布朗運動態勢。

表3 1987—2007年桐柏大別山區主要土地利用類型變化排序表Tab．3 Order ofmain land use conversion during 1987—2007

3.3 土地利用漲落勢圖譜分析

圖3和表4漲落勢圖譜中，耕地、水域以落勢為主，林地、建設用地以漲勢為主，草地和未利用地相對穩定。漲勢系列圖譜記錄了各土地利用類型的增長信息(圖3(a))：新增圖譜單元最多的是林地，3個時段分別占新增土地面積的56.82%、50.67%、61.31%，以耕地和水域減少為依托;其次為耕地，分別占31.55%、36.79%、23.05%，主要由砍伐森林和侵占水域導致;水域、建設用地次之，主要來源于對耕地、林地的侵占;其他用地新增較少。落勢系列圖譜記錄了各土地利用類型的萎縮信息(圖3(b))：耕地萎縮最明顯，3個時段分別占萎縮土地面積的60.24%、53.90%、66.41%，主要轉為林地和城鎮建設用地;其次為林地，分別占30.24%、38.33%、24.43%，主要轉為耕地和水域;水域萎縮次之，主要是被耕地、林地侵占;由于人類活動的強烈干預，建設用地萎縮較小;其他用地萎縮較少。

圖2 桐柏大別山區土地利用變化圖譜Fig．2 Land use changes atlas in Tongbai-Dabiemountainous area

圖3 桐柏大別山區土地利用漲落勢圖譜Fig．3 Arising/declining atlas in Tongbai-Dabiemountainous area

表4 1987—2007年土地利用漲落勢圖譜結構Tab．4 Arising/declining atlas analysis of land use during 1987—2007 km2

4 土地利用結構預測

以1987年各土地利用類型所占面積比例表示初始狀態矩陣A0。以1987—2007年2期土地利用類型轉移概率矩陣為基礎，求出多年加權平均轉移概率矩陣，則該矩陣即為初始狀態的轉移概率矩陣P0。從初始狀態經過n=20、21、22步轉移到2007、2008、2009年轉移概率矩陣，預測各土地利用類型所占的面積比例。在SPSS中對預測值和實測值進行方差分析，顯著性(Sig.=1)大于0.05，表明用馬爾柯夫預測數據與實測情況吻合較好，也說明利用馬爾柯夫模型來預測土地利用變化是可行的。

表5為2010—2019年的預測結果。結果表明：在保持人為干擾強度不變的前提下，桐柏大別山區未來土地利用變化的趨勢是耕地面積繼續減少，至7 370 km2趨于穩定值;林地、建設用地逐年增加，這些變化都是比較劇烈的過程，草地、水域面積在不斷減少，但變化速度相對緩慢。表明未來10年，桐柏大別山區仍將以耕地、林地為主要土地利用類型，其他地類基本沒有明顯變化。

表5 2010—2019年各土地利用類型面積預測值Tab．5 Predictions of land use types from 2010 to 2019

5 結論與討論

1)20年間耕地持續性減少，林地和建設用地持續性增加。

2)耕地與林地之間的相互轉化遍布整個區域，耕地、林地轉為建設用地主要發生在各縣市區周邊經濟相對發達的地區，水域轉為耕地、林地大部分位于水域和其他用地交接地帶。

3)漲勢圖譜單元最多的是林地，占新增面積61.31%，主要來源于耕地和水域;落勢圖譜單元最大的是耕地，占萎縮面積66.41%，主要轉變為林地和建設用地。

4)未來10年，桐柏大別山區仍將以耕地、林地為主要土地利用類型。

在政策導向、經濟發展、比較利益、人口增長等主要驅動因素的影響下，各種土地利用類型發生了不同程度的變化。隨著社會經濟發展，城市化進程加快，城市周邊的優質耕地被擴建為居民區、開發區、交通運輸等建設用地，部分山區應對人口壓力和糧食短缺，通過毀林開荒來增加耕地面積。作為我國重要的糧、油生產基地，遏制耕地減少、保障糧食安全、生態安全、促進區域生態經濟可持續發展是桐柏大別山區的首要任務。因此，結合區域實際情況，首要任務是改進農業生產技術，提高農業生產效率，同時加強耕地資源管理;其次，制定相應的法規政策，合理分配有限的土地資源，著重對當前已利用土地的深度和廣度進行二次開發，正確引導土地利用發展方向;再次，統籌規劃城市的發展，嚴格控制建設用地規模，加強土地資源保護和合理利用，達到土地利用和社會經濟的協調發展。

6 參考文獻

［1］張樹文.東北地區土地利用/覆被時空特征分析［M］.北京：科學出版社，2006

［2］李秀彬.全球環境變化研究的核心領域：土地利用/覆蓋變化的國際研究動向［J］.地理學報，1996，51(6)：553-558

［3］邵懷勇，仙巍，楊武年，等.三峽庫區近50年間土地利用/覆被變化［J］.應用生態學報，2008，19(2)：453-458

［4］Jin S D，Ke W，Yang H，et al.Spatio-temporal dynamics and evolution of land use change and landscape pattern in response to rapid urbanization［J］.Landscape and Urban Planning，2009，92：187-198

［5］Pinkim，Jane S.Protection vs.commercialmanagement：Spatial and temporal analysis of land cover changes in the tropical forests of Central India［J］.Forest Ecology andmanagement，2010，259：1009-1017

［6］張國坤，鄧偉，張洪巖，等.新開河流域土地利用格局變化圖譜分析［J］.地理學報，2010，65(9)：1111-1120

［7］鮑文東，侯志華，吳泉源.基于地學信息圖譜的土地利用動態變化研究：以山東省龍口市為例［J］.地域研究與開發，2007，26(3)：80-84

［8］葉慶華，劉高煥，田國良，等.黃河三角洲土地利用時空復合變化圖譜分析［J］.中國科學：D輯：地球科學，2004，34(5)：461-474

［9］葉慶華，劉高煥，姚一鳴，等.黃河三角洲新生濕地土地利用變化圖譜［J］.地理科學進展，2003，22(2)：141-148

［10］葉慶華，劉高煥，陸洲，等.基于GIS的時空復合體：土地利用變化圖譜模型研究方法［J］.地理科學進展，2002，21(4)：349-357

［11］陳述彭.地學信息圖譜探索研究［M］.北京：商務印書館，2001

［12］閆正龍，黃強，暢建霞，等.基于3S技術的塔里木河干流土地利用動態監測［J］.水土保持學報，2008，22(1)：190-193

Atlas analysis on spatial-temporal changes of land use in Tongbai-Dabiemountainous area

Zhang Ronghua1，2，Liu Xia1，Li Huan1，Yao Xiaoyou3，Wang Yousheng1，Zhang Guangcan1
(1.Forestry College of Shandong Agricultural University，Key Laboratory of Soil Erosion and Ecological Restoration，Taishan Forest Ecosystem Research Station，271018，Tai’an，Shandong;2.School of Geography，Beijing Normal University，100875，Beijing;3.Soil and Water Conservation Division of Huaihe River Commission of Water Resourcesministry，233001，Bengbu，Anhui：China)

Based on GIS spatial and temporal analysismodel，geo-information atlas theory were used in this paper.Land use data were extracted from spatial data of three periods(1987，1997，2007)in Tongbai-Dabiemountainous area.Land use change atlas and arising/declining atlas from 1987 to 2007 were established.Land use changes characteristic were analyzed，and land use changes were predicted.Results indicate that land use changes significantly during the past 20 years.The area of farmland decreased，and forest land and construction land increased.Grassland and unused land increased firstly，then decreased，vice versewater land decreased firstly，then increased.From the spatial series，mutual transformation between farmland and forest land was obvious in the whole research area.Farmland and forest land transfered to construction land in the nearby of town.Waters transfered to farmland and forest land in the nearby of waters.In arising/declining atlas，the largest expanding atlas was forest land，with an increase of 5904.06 km2，accounting for 61.31%in expanding atlas，mainly from farmland and water land.And the largest shrinking atlas was farmland，with a decrese of 6 394.71 km2，accounting for 66.41%in shrinking atlas，mainly transforming into forest land and construction land.In the future，farmland and forest land will continue to be chief land use type.

geo-information atlas;spatio-temporal changes;land use;Tongbai-Dabiemountainous area

2011-02-18

2011-05-12

水利部淮河水利委員會“淮河流域生態修復機理和評價指標體系研究”(hwsbc-2007-001)、“淮河流域水土保持區劃及分區防治模式研究”(SBJ-2010-002);林業公益性行業科研專項“生境脅迫立地植被恢復與重建技術研究”(201104002-6)

張榮華(1984—)，女，博士研究生。主要研究方向：水土保持信息技術。E-mail：zrhua5766@163.com

?責任作者簡介：劉霞(1971—)，女，教授，博士。主要研究方向：土壤侵蝕與生態修復。E-mail：liuxia@sdau.edu.cn

(責任編輯：宋如華)