方應(yīng)波，張中旺，熊宏濤

?

近30年武漢市土地利用景觀格局動態(tài)變化

方應(yīng)波1，張中旺2，熊宏濤3

(1.廣州南洋理工職業(yè)學(xué)院，廣東廣州 510925；2.湖北文理學(xué)院管理學(xué)院，湖北襄陽 441053；3.湖北國土資源職業(yè)學(xué)院，湖北武漢 430090)

基于1980s、1990s、2005年Landsat TM遙感影像和2010年歐空局(ESA)Globcover分類數(shù)據(jù)，運用RS和GIS技術(shù)及景觀生態(tài)學(xué)方法，分析了1980s—2010年武漢市土地利用類型、土地利用動態(tài)度、土地利用景觀格局的動態(tài)變化. 結(jié)果表明：農(nóng)地是主要的土地利用類型，地塊連續(xù)性增強，成片分布；林地和草地1990s破壞嚴(yán)重，2010年區(qū)域植被得到一定程度的恢復(fù)；水域所占比重持續(xù)下降，減少速度之快，觸目驚心；受城市化進程影響，1990s—2010年建設(shè)用地增長最快，呈規(guī)?；?、成片發(fā)展；1980s—1990s未利用地得到大幅度開發(fā)，2010年未利用地面積明顯增加且趨于成片化發(fā)展，地塊連續(xù)性增強.

土地利用；景觀格局；武漢市

1995年由國際地圈-生物圈計劃(IGBP)和全球環(huán)境變化中的人文領(lǐng)域計劃(HDP)聯(lián)合提出“土地利用/土地覆蓋”(Land use and Land cover change, LUCC)研究計劃以來，引起國內(nèi)外學(xué)者的廣泛關(guān)注，此后LUCC成為全球環(huán)境變化研究的熱點[1-2]，其涵蓋的內(nèi)容包括以下3個方面：1)利用Landsat TM/ETM影像，經(jīng)人機交互遙感圖像解譯后，對土地利用變化進行監(jiān)測，評價與分析制圖. 如李謝輝等以河南省1980s和2005年兩期影像為數(shù)據(jù)源，對1980s—2005年土地利用／覆蓋主要類型的演變進行了現(xiàn)狀和動態(tài)分析[3]；邢容容等以1990年、2000年、2006年3期的Landsat TM影像和2011年衛(wèi)星影像，分析了1990—2011年青島市土地利用/覆被變化情況[4]. 2)進行土地利用/土地覆被變化動力學(xué)分析. 如李松等通過對新疆18個綠洲城市的1995年、2000年、2005年、2010年土地利用景觀演變發(fā)展過程及驅(qū)動力進行分析，認(rèn)為人文因素是綠洲城市土地利用景觀演變的主導(dǎo)動力[5]；潘竟虎等采用四期Landsat TM遙感影像對玉門市近30年土地利用的變化及其驅(qū)動力進行了分析[6]. 3)開展土地利用/土地覆被與生態(tài)環(huán)境，可持續(xù)發(fā)展研究. 艾建超等以2005年和2009年的土地利用遙感影像圖，通過構(gòu)建景觀格局指標(biāo)及生態(tài)風(fēng)險評價指標(biāo)，分析了吉林省鎮(zhèn)賚縣景觀格局變化和生態(tài)風(fēng)險[7]；韓黎陽等分析了三峽庫區(qū)蘭陵溪小流域土地利用景觀格局對水體氮、磷等輸出時空特征及對其產(chǎn)生的影響[8].

2005年以來，伴隨著武漢城市圈的建立和發(fā)展，尤其是2007年國務(wù)院正式批準(zhǔn)武漢城市圈成為“兩型社會”建設(shè)綜合配套改革試驗區(qū)，武漢市如何發(fā)揮科學(xué)發(fā)展與社會和諧的示范和帶動作用，協(xié)調(diào)區(qū)域經(jīng)濟發(fā)展與環(huán)境保護的關(guān)系，不僅關(guān)系到武漢未來經(jīng)濟發(fā)展的方向，更關(guān)系到全國新型工業(yè)化、城市化發(fā)展成功的關(guān)鍵. 區(qū)域經(jīng)濟的發(fā)展勢必會改變土地利用現(xiàn)狀及景觀格局，繼而影響生態(tài)環(huán)境，而LUCC是全球環(huán)境變化的主要原因和組成部分之一，研究LUCC對了解區(qū)域生態(tài)環(huán)境變化、維持生態(tài)平衡、促進區(qū)域經(jīng)濟與環(huán)境的協(xié)調(diào)發(fā)展具有重要意義[9]. 因此，從長時間尺度研究武漢市土地利用結(jié)構(gòu)變化規(guī)律及現(xiàn)狀特征，是為“兩型社會”建設(shè)和武漢市區(qū)域經(jīng)濟發(fā)展摸清家底，提供決策基礎(chǔ). 本文基于GIS/RS技術(shù)，以土地利用景觀格局為切入點，對武漢市1980s—2010年土地利用類型、土地利用動態(tài)度、土地利用景觀格局進行定量分析.

1 數(shù)據(jù)來源及處理

1.1 數(shù)據(jù)來源

文章涉及的遙感影像數(shù)據(jù)為武漢市1980s(1984年和1985年多景遙感影像拼接裁剪)、1990s(1995年和1996年多景遙感影像拼接裁剪)、2005年TM遙感影像和2010年歐空局(ESA)Globcover分類數(shù)據(jù)(分辨率為300*300m)，根據(jù)研究需要將其合并為農(nóng)地、林地、草地、水域、建設(shè)用地和未利用地六種類型；武漢市矢量數(shù)據(jù)從全國1:400萬基礎(chǔ)地理數(shù)據(jù)庫中提取，包括shapfile格式面狀文件及邊界，保留基本地理屬性表信息.

1.2 數(shù)據(jù)處理

對武漢市的遙感影像進行幾何校正及波段處理，拼接裁剪并在此基礎(chǔ)上濾波后，再使用城市矢量面狀文件進行掩膜裁剪，得到研究區(qū)域完整的Landsat TM遙感影像數(shù)據(jù)，在此基礎(chǔ)上通過圖像增強處理，運用監(jiān)督分類的方法進行采樣，建立分類模板，圖像解譯，最后進行精度評價. 加上2010年歐空局(ESA)Globcover分類數(shù)據(jù)，生成武漢市四時相土地利用數(shù)據(jù)(農(nóng)地、林地、草地、水域、建設(shè)用地和未利用地)，然后將土地利用數(shù)據(jù)轉(zhuǎn)化為景觀格局?jǐn)?shù)據(jù)，利用Fragstats軟件計算武漢市景觀格局指數(shù)，分析景觀結(jié)構(gòu)特征.

2 結(jié)果與分析

2.1 土地利用類型動態(tài)變化

根據(jù)遙感影像解譯結(jié)果，在ARCGIS9.3中分別導(dǎo)出4個時期土地利用類型的統(tǒng)計數(shù)據(jù)，再將數(shù)據(jù)導(dǎo)入數(shù)據(jù)庫軟件ACCESS，然后利用數(shù)據(jù)庫的制表功能進行分類統(tǒng)計，最后獲得武漢市不同時期各土地利用類型的面積及比重，如表1所示. 由于不同的遙感影像分辨率，解譯誤差等主客觀因素限制，解譯結(jié)果必然存在相應(yīng)誤差，最后統(tǒng)計出4個時期的總面積分別為：859 698.70hm2、859 700.14hm2、859 698.70hm2、860 205.05hm2，其誤差相對于80多萬hm2的總面積來說可以忽略不計，解譯結(jié)果符合研究需要.

表1 1980s—2010年武漢市不同土地利用類型的面積及比重

表1反映了1980s—2010年武漢市六類土地利用類型面積及比重的變化情況，從中可以看出：按土地利用類型所占的比重排序依次是農(nóng)地、水域、林地、建設(shè)用地、草地、未利用地，其中農(nóng)地是主要的土地利用類型，占土地總面積的70%以上；從近30年土地利用的變化趨勢看，水域所占比重雖然偏高，但一直在持續(xù)下降，其比重從1980s的19.35%減少到2010年的11.86%，減少速度之快觸目驚心；農(nóng)地所占比重一直偏高但變化不明顯，1990s之前農(nóng)地開發(fā)速度較快，之后農(nóng)地開發(fā)速度逐漸放緩；林地在1990s之前破壞嚴(yán)重，其比重由1980s的8.25%下降到1990s的3.97%，1990s之后植被得到一定程度的恢復(fù)；草地所占比重一直偏低，且2005年以前持續(xù)緩慢下降，2005年以后草地得到一定程度的恢復(fù)，因前一時期比重過低，故2010年草地比重增幅達4.9倍；2005年之前建設(shè)用地所占比重持續(xù)增長，2010年以后增長速度放緩，其中2005年建設(shè)用地比重達近30年最高值；未利用地所占比重一直偏低，呈波動緩慢變化，1990s其比重達近30年最低值，說明1990s以前土地開發(fā)程度最強.

2.2 土地利用動態(tài)度測算

土地利用動態(tài)變化包括土地資源的數(shù)量、質(zhì)量隨時間的變化，也包括土地利用的空間變化及土地利用類型組合方式的變化. 本文主要利用遙感影像得到的土地數(shù)量數(shù)據(jù)，分析武漢市各類土地利用類型的數(shù)量變化. 土地利用數(shù)量變化可以用土地利用動態(tài)度來定量描述，分為單一土地利用動態(tài)度和綜合土地利用動態(tài)度[10-12]，文章選取單一土地利用動態(tài)度進行測算. 單一土地利用動態(tài)度以定量方式表達區(qū)域一定時間內(nèi)某一類用地面積的變化幅度與變化速度，以及區(qū)域土地利用變化中的類型差異，其公式[10-12]為：

式中，為研究時段內(nèi)某一種土地利用類型動態(tài)度，U、U分別為研究期初和期末某種土地利用類型的數(shù)量，為研究時段長. 當(dāng)?shù)臅r段設(shè)定為年時，值就是該研究區(qū)域內(nèi)某種土地利用類型的年變化率. 據(jù)此得出武漢市各土地利用類型的動態(tài)度變化情況如表2所示.

表2 1980s—2010年武漢市單一土地利用類型動態(tài)度 /%

表2反映了1980s—2010年武漢市單一土地動態(tài)度的變化情況，從中可以得出：1980s—1990s土地利用類型減少的有林地、草地、水域、未利用地，其中未利用地減少速度最快，其它兩類面積增加，其中建設(shè)用地增加最快；1990s—2005年土地利用類型減少的有農(nóng)地、草地、水域，林地得到一定程度的恢復(fù)、建設(shè)用地與未利用地快速增長；2005—2010年土地利用類型減少的有林地、水域、建設(shè)用地，農(nóng)地、草地、未利用地均有不同程度增長. 綜合近30年的整體情況，水域面積一直持續(xù)減少，2005—2010年速度有所放緩，其它土地利用面積交錯增減，其中建設(shè)用地變化幅度最大.

2.3土地利用景觀格局動態(tài)變化

在研究武漢市土地景觀格局時，文章選取拼塊類型面積(CA)、拼塊所占景觀面積的比例(%LAND)、斑塊個數(shù)(NP)、面積加權(quán)平均斑塊面積(AREA-MN)、最大斑塊面積百分比(Largest Patch Index，LPI)、周長面積分維指數(shù)(Perimeter-Area Fractal Dimension)幾項指數(shù)進行分析[13-15].

在ARCGIS 9.3 軟件環(huán)境下，將武漢市土地利用數(shù)據(jù)(*IMG格式)轉(zhuǎn)化為矢量數(shù)據(jù)(SHAPE格式)，進一步轉(zhuǎn)為柵格GRID格式，將轉(zhuǎn)化后的土地利用數(shù)據(jù)導(dǎo)入FRSGATATS3.5，選擇以上篩選的景觀指數(shù)進行計算，得到武漢市不同時期景觀指數(shù)如表3所示. 表3反映了武漢市1980s—2010年不同景觀在斑塊類型水平指數(shù)上的變化情況.

表3 1980s—2010年武漢市不同景觀斑塊類型的水平指數(shù)

農(nóng)地景觀變化 自1980s—2010年農(nóng)地景觀拼塊類型面積CA和拼塊所占景觀面積的比例%LAND均較大，變化趨勢為先增加后下降且增加的幅度超過下降的幅度，說明農(nóng)地一直是優(yōu)勢景觀，其變化趨勢是農(nóng)業(yè)在產(chǎn)業(yè)結(jié)構(gòu)調(diào)整中地位的反映；農(nóng)地斑塊個數(shù)NP持續(xù)下降，面積加權(quán)平均斑塊面積AREA-MN持續(xù)上升，說明農(nóng)地景觀破碎程度越來越低，地塊連續(xù)性增強，成片分布；農(nóng)地最大斑塊面積百分比LPI總趨勢上升，2010年略有下降，周長面積分維指數(shù)Perimeter持續(xù)下降，說明人類對農(nóng)地景觀的影響在逐漸增強，由于人為平整土地，農(nóng)地斑塊形狀更有規(guī)律，斑塊的幾何形狀趨于簡單.

林地景觀變化 自1980s—2010年林地景觀拼塊類型面積CA和拼塊所占景觀面積比例%LAND均較小，說明林地不是區(qū)域優(yōu)勢景觀，其變化趨勢反映1990s植被破壞嚴(yán)重，2010年區(qū)域植被得到一定程度的恢復(fù)；林地斑塊個數(shù)NP均較大且持續(xù)下降，面積加權(quán)平均斑塊面積AREA-MN先降后升，表明林地景觀整體分散零碎，其中1990s因植被破壞，林地最為破碎；林地最大斑塊面積百分比LPI均較小且持續(xù)下降，周長面積分維指數(shù)Perimeter均較大，表明林地受人類影響，零碎分布，斑塊形狀極不規(guī)則.

草地景觀變化 自1980s—2010年草地景觀拼塊類型面積CA和拼塊所占景觀面積的比例%LAND均先減小后快速增加，1990sCA與%LAND均達最小值，2010年CA與%LAND快速增長，其變化趨勢表明草地于1990s遭受破壞，零碎分布，2010年后草地恢復(fù)速度較快；草地面積加權(quán)平均斑塊面積AREA-MN與最大斑塊面積百分比LPI均先降后升，其變化趨勢同樣可以反映出草地在1990s遭受破壞，2010年后草地迅速恢復(fù)且呈連片分布；草地周長面積分維指數(shù)Perimeter持續(xù)下降，說明受人類干擾程度越來越強.

水域景觀變化 自1980s—2010年水域景觀拼塊類型面積CA、拼塊所占景觀面積的比例%LAND、斑塊個數(shù)NP、最大斑塊面積百分比LPI、周長面積分維指數(shù)Perimeter五項指標(biāo)均持續(xù)下降，表明水域受人類影響最強烈，水域面積和水域景觀個數(shù)減少嚴(yán)重；水域拼塊所占景觀面積的比例由1980s的19.42%下降到2010年的11.89%，“百湖之市”的稱謂已成歷史；水域面積加權(quán)平均斑塊面積AREA-MN持續(xù)上升，1990s—2010年AREA-MN由226.55上升到444.49，說明小湖泊消失較多，1990s—2010年水域面積萎縮最快.

建設(shè)用地景觀變化 自1980s—2010年建設(shè)用地景觀拼塊類型面積CA、拼塊所占景觀面積的比例%LAND和最大斑塊面積百分比LPI均持續(xù)上升，反映城市化進程的步伐，其中1990s—2010年建設(shè)用地快速增加；建設(shè)用地斑塊個數(shù)NP、面積加權(quán)平均斑塊面積AREA-MN波動變化，其中1990s—2010年變化最為明顯，表明此階段建設(shè)用地呈規(guī)?；?、成片發(fā)展；建設(shè)用地周長面積分維指數(shù)Perimeter整體偏小且持續(xù)下降，2010年該指標(biāo)達六項景觀類型的最小值，表明建設(shè)用地景觀形狀最有規(guī)律，受人類干擾最大.

未利用地景觀變化 自1980s—2010年未利用地景觀拼塊類型面積CA、拼塊所占景觀面積的比例%LAND波動變化，其中1980s—1990s變化最快，表明該時段未利用地得到大幅度的開發(fā)，2010年未利用地面積明顯增加；未利用地斑塊個數(shù)NP和面積加權(quán)平均斑塊面積AREA-MN變化表明未利用地趨于成片化發(fā)展，地塊連續(xù)性增強；最大斑塊面積百分比LPI、周長面積分維指數(shù)Perimeter變化幅度不大，變化趨勢同樣反映出2010年地塊呈片狀發(fā)展.

3 結(jié)語

基于GIS/RS技術(shù)，通過定量分析1980s—2010年武漢市土地利用類型、土地利用動態(tài)度、土地利用景觀格局三個方面的變化情況，主要結(jié)論如下：

1)土地利用類型變化. 從土地利用類型所占的比重看，農(nóng)地是主要的土地利用類型，占土地總面積的70%以上；從近30年土地利用的變化趨勢看，水域所占比重近30年一直持續(xù)下降，減少速度之快觸目驚心；農(nóng)地在1990s之前開發(fā)速度較快，之后開發(fā)速度逐漸放緩；林地在1990s之前破壞嚴(yán)重，1990s之后植被得到一定程度的恢復(fù)；草地所占比重一直偏低，2005年以前持續(xù)緩慢下降，2005年以后草地得到恢復(fù)；建設(shè)用地2005年之前持續(xù)增長，其中2005年建設(shè)用地比重達近30年最高值；未利用地所占比重一直偏低，呈波動緩慢變化.

2)土地利用動態(tài)變化度. 水域面積一直持續(xù)減少，2005—2010年變化速度有所放緩，其它土地利用動態(tài)變化度交錯增減，其中建設(shè)用地變化幅度最大.

3)土地景觀格局變化. 農(nóng)地景觀破碎程度越來越低，地塊連續(xù)性增強，成片分布，同時由于人為平整土地，農(nóng)地斑塊形狀更有規(guī)律，斑塊的幾何形狀趨近于簡單；林地1990s破壞嚴(yán)重，2010年區(qū)域植被得到一定程度的恢復(fù)，林地景觀整體分散零碎，斑塊形狀極不規(guī)則；草地于1990s遭受破壞，零碎分布，2010年后草地恢復(fù)速度較快且呈連片分布；水域受人類影響最強烈，其中1990s—2010年水域面積萎縮最快，水域面積和水域景觀個數(shù)減少嚴(yán)重，小湖泊消失較多，“百湖之市”的稱謂已成歷史；由于受城市化進程影響，1990s—2010年建設(shè)用地增長較快，此階段建設(shè)用地呈規(guī)模化、成片發(fā)展，其景觀形狀趨于簡單；1980s—1990s未利用地得到大幅度開發(fā)，2010年未利用地面積明顯增加且趨于成片化發(fā)展，地塊連續(xù)性增強.

[1] TUNER II B L, SKOLE D, SANDERSON S, et al. Land-use and land-cover Change Science /Research Plan[R]. Stockholm: Royal Swedish Academy of Sciences, 1995.

[2] 陳佑啟, 楊 鵬. 國際上土地利用/土地覆蓋變化研究的新進展[J]. 經(jīng)濟地理, 2001, 21(1): 95-100.

[3] 李謝輝, 王 磊. 河南省土地利用/覆蓋景觀格局變化分析[J]. 水土保持研究, 2012, 19(5): 83-89.

[4] 邢容容, 劉樂軍, 馬安青. 青島市1990—2011年土地利用/覆被變化研究[J]. 水土保持通報, 2013, 33(3): 220-225.

[5] 李 松, 張小雷, 李壽山, 等. 新疆綠洲城市土地利用景觀格局時空演化及驅(qū)動力研究[J]. 經(jīng)濟地理, 2013, 33(12): 161-168.

[6] 潘竟虎, 蘇有才, 黃永生, 等. 近30年玉門市土地利用與景觀格局變化及其驅(qū)動力[J]. 地理研究, 2012, 31(9): 1631-1639.

[7] 艾建超, 趙 力, 程龍飛, 等. 土地整理進程中鎮(zhèn)賚項目區(qū)景觀格局變化及生態(tài)風(fēng)險[J]. 東北師大學(xué)報: 自然科學(xué)版: 2013, 45(2): 146-152.

[8] 韓黎陽, 黃志霖, 肖文發(fā), 等. 三峽庫區(qū)蘭陵溪小流域土地利用及景觀格局對氮磷輸出的影響[J]. 環(huán)境科學(xué), 2014, 35(3): 1091-1097.

[9] 于興修, 楊桂山. 中國土地利用/覆被變化研究的現(xiàn)狀與問題[J]. 地理科學(xué)進展, 2002, 21(1): 51-57.

[10] 李仁東, 程學(xué)軍, 隋曉麗. 江漢平原土地利用的時空變化及其驅(qū)動因素分析[J]. 地理研究, 2003, 22(4): 423-431.

[11] 王秀蘭, 包玉海. 土地利用動態(tài)變化研究方法探討[J]. 地理科學(xué)進展, 1999, 18(1): 81-87.

[12] 羅湘華, 倪晉仁. 土地利用/土地覆蓋變化研究進展[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報, 2000, 8(3): 262-272.

[13] 何 鵬, 張會儒. 常用景觀指數(shù)的因子分析和篩選方法研究[J]. 林業(yè)科學(xué)研究, 2009, 22(4): 470-474.

[14] 彭 建, 王仰麟, 張 源, 等. 土地利用分類對景觀格局指數(shù)的影響[J]. 地理學(xué)報, 2006, 61(2): 157-168.

[15] 周青山, 彭軍還, 余德清, 等. 基于遙感信息的區(qū)域景觀指數(shù)自動提取[J]. 西部探礦工程, 2008, 20(3): 67-71.

Dynamic Changes of Land Use and Landscape Pattern in Wuhan City in the Past 30 Years

FANG Yingbo1,ZHANG Zhongwang2,XIONG Hongtao3

(1.Guangzhou Nanyang College, Guangzhou 510925, China;2. School of Management, Hubei University of Arts and Science, Xiangyang 441053, China;3. Hubei Land Resources Vocational College, Wuhan 430090, China)

Combined with the integrated technology of RS/GIS technology and landscape ecological method, based on the Landsat TM image of 1980s, 1990s, 2005 and ESA’ Globcover data of 2010, this paper analyzed the changes of land use/cover and landscape pattern of Wu’han City form 1980s to 2010. The results show that agricultural land, which is continuous, is the main land use type. Forestland and grassland, which were partly restored in 2010, were damaged seriously in1990s. The proportion of waters decreased continuously from 1980s to 2010. With the promotion of urbanization process, construction land increased most rapidly from1990s to 2010. Unutilized land, which increased gradually large-scalely and continuously in 2010, developed greatly from 1980s to 1990s.

Land use; Landscape pattern; Wuhan City

(責(zé)任編輯：陳 丹)

F291.1

A

2095-4476(2015)02-0034-05

2014-10-13;

2015-01-05

湖北省科技支撐計劃項目(2014BDF044); 湖北省社會科學(xué)基金項目(2012088)

方應(yīng)波(1980— ), 男, 湖北大悟人, 廣州南洋理工職業(yè)學(xué)院講師.

張中旺(1966— ), 男, 湖北孝感人, 湖北文理學(xué)院管理學(xué)院教授, 博士, 主要研究方向：區(qū)域水安全與可持續(xù)發(fā)展.