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

（1．山東農(nóng)業(yè)大學(xué) 林學(xué)院 水土保持系，山東農(nóng)業(yè)大學(xué) 農(nóng)業(yè)生態(tài)與環(huán)境重點(diǎn)實(shí)驗(yàn)室，山東 泰安271018；2．水利部淮河水利委員會(huì)水土保持處，安徽 蚌埠233001）

景觀格局是景觀異質(zhì)性的具體表現(xiàn)，同時(shí)又是包括干擾在內(nèi)的各種生態(tài)過(guò)程在不同尺度上作用的結(jié)果［1－2］。景觀格局的形成及其變化是自然、社會(huì)和生物等多種要素相互作用的結(jié)果。景觀要素在空間分布的數(shù)量、位置、類(lèi)型、形狀、大小和方向構(gòu)成了景觀空間格局分析的主要內(nèi)容。景觀結(jié)構(gòu)和功能、格局與過(guò)程之間的聯(lián)系與反饋是景觀生態(tài)學(xué)的基本命題。用景觀指數(shù)描述景觀格局及變化，建立格局與景觀過(guò)程之間的聯(lián)系，是景觀生態(tài)學(xué)最常用的定量化研究方法［3－7］。本文采用1987年、2007年兩個(gè)典型時(shí)期4－5月的T M 影像為數(shù)據(jù)源，在Region Manager 5．5、Arc GIS和Fragstats 3．3軟件的支持下研究1987－2007年20 a間蒙陰縣土地利用景觀格局的特點(diǎn)及其動(dòng)態(tài)變化趨勢(shì)，分析各景觀要素類(lèi)型之間的空間轉(zhuǎn)化關(guān)系，反映景觀要素之間相互作用的性質(zhì)、強(qiáng)度和方式，進(jìn)而闡明景觀要素動(dòng)態(tài)演替的方向、過(guò)程，為實(shí)現(xiàn)流域可持續(xù)發(fā)展提供科學(xué)依據(jù)與決策支持。

1 研究區(qū)概況

蒙陰縣位于山東省中南部，117°45′－118°15′E，35°27′－36°02′N(xiāo)，地處蒙山北麓，東汶河上游，是沂蒙山區(qū)的腹地。東鄰沂水縣、沂南縣，西靠新泰市，南依蒙山與費(fèi)縣、平邑縣交界，北與沂源縣接壤，總面積1 610．85 k m2，為典型的北方土石山區(qū)，海拔一般為200～600 m，基巖主要由石灰?guī)r和頁(yè)巖組成，土壤類(lèi)型主要為棕壤類(lèi)。屬暖溫帶季風(fēng)型大陸氣候，年均日照時(shí)間2 257 h，年均氣溫12．8℃，無(wú)霜期200 d。年均降水量820 mm，降水分布不均，有70%主要集中分布在5－9月。地帶性植被類(lèi)型為暖溫帶落葉闊葉林，主要喬木樹(shù)種有側(cè)柏（Platycladus orientalis Linn．Franco）、刺槐（Robinia pseudoacacia L．）、油松（Pinus tabul aef or mis Carr．）、黃櫨（Cotinus coggyria Scop．）、遼東櫟（Quercus liaotungensis Koidz．）、山楂（Crataegus pinnatif ida Bge．）、柿樹(shù)（Diospyros kaki Linn F．）等；自然灌木與草本植物主要有黃荊（Vitex negundo L．）、胡枝子（Lespedeza bicolors）

2 研究方法

2．1 數(shù)據(jù)源

遙感數(shù)據(jù)源采用1987年、2007年兩個(gè)典型時(shí)期4－5月的T M影像，空間分辨率為30 m×30 m。蒙陰縣1996年土地利用圖，蒙陰縣1∶10萬(wàn)地形圖作為影像匹配基準(zhǔn)。

2．2 數(shù)據(jù)處理

（1）圖像處理。為消除遙感數(shù)據(jù)中存在的誤差，提高可判讀性，借助ERDAS軟件進(jìn)行配準(zhǔn)、鑲嵌、勻色及增強(qiáng)處理。利用地形圖選擇25～35個(gè)同名地面控制點(diǎn)（GCP），對(duì)影像進(jìn)行校正，誤差控制在0．5個(gè)像元以?xún)?nèi)。校正中，采用北京1954坐標(biāo)系，高斯克呂格投影。

（2）解譯標(biāo)志的建立與圖像解譯。進(jìn)行野外考察并結(jié)合已有資料建立解譯標(biāo)志。本研究采用自動(dòng)分類(lèi)和人工判讀相結(jié)合的方法對(duì)研究區(qū)T M圖像進(jìn)行解譯，獲取景觀要素類(lèi)型等數(shù)據(jù)，運(yùn)用Arc GIS軟件建立數(shù)據(jù)庫(kù)。

（3）景觀指數(shù)的提取。利用 Arc GIS 9．2的空間分析功能及Fragstats 3．3軟件提取部分景觀格局指數(shù)，對(duì)蒙陰縣近20 a來(lái)的景觀格局動(dòng)態(tài)變化情況進(jìn)行分析，并分析造成其變化的原因。

2．3 景觀要素類(lèi)型劃分

本文主要依據(jù)全國(guó)土地利用現(xiàn)狀分類(lèi)系統(tǒng)，并結(jié)合研究目的和區(qū)域特征，將蒙陰縣整個(gè)區(qū)域內(nèi)的景觀類(lèi)型劃分為耕地、林地、草地、城鎮(zhèn)村及工礦交通用地、水域及水利設(shè)施用地和未利用土地6類(lèi)［8－10］。

2．4 景觀要素格局指數(shù)選擇與提取

通過(guò)景觀指數(shù)描述景觀格局，使數(shù)據(jù)獲得一定的統(tǒng)計(jì)性質(zhì)，可以比較、分析不同尺度上的格局，長(zhǎng)期以來(lái)一直倍受景觀生態(tài)學(xué)研究者的關(guān)注。對(duì)于景觀空間格局及其變化的考察，一套合理的指標(biāo)體系的選擇主要從景觀水平上和類(lèi)型水平上分析，盡可能地選擇相互獨(dú)立的景觀指數(shù)，且能較全面地描述土地利用景觀格局的數(shù)量和空間特征［11－13］。

Fragstats是一個(gè)功能強(qiáng)大的景觀指數(shù)計(jì)算軟件，可以計(jì)算出上百種景觀指數(shù)［14－16］。本文運(yùn)用Arc GIS的空間分析功能將文件轉(zhuǎn)化為grid格式的文件，然后利用Fragstats軟件提取景觀、類(lèi)型水平上的部分景觀指數(shù)，從而對(duì)研究區(qū)的景觀格局動(dòng)態(tài)進(jìn)行分析，本文選取的指標(biāo)及其計(jì)算公式與生態(tài)含義見(jiàn)表1。

3 結(jié)果與分析

3．1 景觀要素類(lèi)型面積的變化

研究區(qū)的總面積1 610．85 k m2，1987年和2007年兩期蒙陰縣的景觀類(lèi)型均以耕地、林地、草地為主，三種土地利用類(lèi)型占總面積的90%；其次是水域及水利設(shè)施用地，占4%左右；未利用土地較少，僅占1%左右。由此可知該區(qū)域植被覆蓋率高。1987年耕地面積幾乎占據(jù)整個(gè)蒙陰縣面積的一半，可見(jiàn)研究區(qū)農(nóng)業(yè)發(fā)展較快；2007年耕地、林地面積均為研究區(qū)的30%以上，水資源缺乏，制約著蒙陰縣的經(jīng)濟(jì)發(fā)展。

蒙陰縣在1987－2007年20 a間，隨著對(duì)各種土地資源的開(kāi)發(fā)，各類(lèi)景觀要素的面積發(fā)生了很大變化（表2）。林地、城鎮(zhèn)村及工礦交通用地，水域及水利設(shè)施用地面積增加，耕地等其他用地面積減少。從變化數(shù)量看，最為顯著的是林地的增加和草地、耕地的減少。20 a來(lái)林地面積凈增加133．35 k m2，平均每年增加6．67 k m2，林地面積的大量增加是大面積開(kāi)墾草地與未利用土地的結(jié)果，尤其是開(kāi)發(fā)條件較差的未利用土地與草地，隨著農(nóng)業(yè)機(jī)械化程度與交通條件的便利，使許多開(kāi)發(fā)初期不便開(kāi)發(fā)的宜林荒地得到了進(jìn)一步開(kāi)發(fā)。草地凈減少43．64 k m2，平均每年減少2．18 k m2，草地不斷減少主要是由于大面積林地種植，尤其是自然條件良好、交通相對(duì)方便、人口密度較大的區(qū)域，大量適宜開(kāi)墾的草地被開(kāi)發(fā)為林地。耕地凈減少107．49 k m2，平均每年減少5．37 k m2。為提高土地生產(chǎn)力，保證農(nóng)業(yè)的穩(wěn)定發(fā)展，蒙陰縣采取了一系列措施，使農(nóng)業(yè)技術(shù)裝備水平顯著提高，由傳統(tǒng)農(nóng)業(yè)逐步向現(xiàn)代化農(nóng)業(yè)轉(zhuǎn)變。使人們開(kāi)始尋找提高收入的新的土地經(jīng)營(yíng)方式，農(nóng)業(yè)的集約化程度越來(lái)越高，使得不需要大面積的耕地即可供養(yǎng)更多的人口，所需要的耕地面積減少。從土地利用景觀類(lèi)型年變化率來(lái)看，未利用土地的年變化率最大，平均每年減少2．83%，其次是城鎮(zhèn)村及工礦交通用地，平均每年增加2．52%，其面積膨脹較大，是以耕地的減少為依托的。林地年平均增加1．47%，草地年平均減少0．75%。水域及水利設(shè)施用地年平均增加0．74%，變化相對(duì)較小。耕地年平均減少0．74%，雖然耕地的年變化率低于城鎮(zhèn)村及工礦交通用地，但由于其面積比重大，因此變化的總面積僅低于林地。

表1 景觀格局指數(shù)

表2 1987－2007年期間不同時(shí)期蒙陰縣不同景觀類(lèi)型的變化

3．2 景觀格局動(dòng)態(tài)分析

3．2．1 景觀水平上格局動(dòng)態(tài)分析 從對(duì)研究區(qū)總體景觀格局來(lái)看（表3），1987年斑塊總量為792塊，2007年為3 318塊，斑塊數(shù)量呈現(xiàn)顯著增加。斑塊密度PD大幅度增加，最大斑塊面積百分比LPI明顯減小，說(shuō)明1987－2007年研究區(qū)景觀破碎化程度逐漸加大。研究區(qū)破碎化程度的加大，表明在經(jīng)濟(jì)發(fā)展和城市化水平快速推進(jìn)的背景下，人類(lèi)干擾活動(dòng)日益增多和深入，代表自然生態(tài)系統(tǒng)的土地覆被類(lèi)型大幅度減少，景觀的異質(zhì)化程度呈上升趨勢(shì)，從而使自然生態(tài)環(huán)境的凈化和維系能力下降，適于生物生存的環(huán)境減少，直接影響到物種的繁殖、擴(kuò)散、遷移和保護(hù)。研究區(qū)域1987年分維數(shù)為1．289 3，2007年為1．473 8，增長(zhǎng)了0．184 5。分維數(shù)越大，表明斑塊形狀越復(fù)雜。研究時(shí)間段內(nèi)，分維數(shù)變大，說(shuō)明從總體而言，斑塊形狀趨于復(fù)雜和不規(guī)則，這主要是由于景觀要素類(lèi)型受到人類(lèi)活動(dòng)的影響較大，耕地、林地、草地等景觀類(lèi)型，受到建設(shè)用地的侵入和分割，使景觀的整體形狀變得較為復(fù)雜。1987－2007年蔓延度指數(shù)呈下降趨勢(shì)，表明1987年以來(lái)，耕地由優(yōu)勢(shì)景觀、連接性好的景觀逐漸向連接性較差的人工景觀格局演變，非農(nóng)化和人工化情況嚴(yán)重，充分說(shuō)明耕地景觀類(lèi)型的不完整性，呈現(xiàn)出區(qū)域占優(yōu)勢(shì)的耕地在空間上的距離加大，影響到耕地地塊之間物質(zhì)、能量和信息的流通不暢，進(jìn)而影響到區(qū)域的整個(gè)生態(tài)系統(tǒng)。香農(nóng)多樣性指數(shù)和香農(nóng)均勻度指數(shù)在研究時(shí)間段內(nèi)分別增長(zhǎng)了0．011 5和0．006 4，說(shuō)明區(qū)域景觀類(lèi)型的豐富程度和均勻程度在上升，優(yōu)勢(shì)景觀在下降，破碎度增加，說(shuō)明對(duì)區(qū)域生態(tài)起調(diào)節(jié)作用的主要景觀類(lèi)型對(duì)區(qū)域的生態(tài)調(diào)節(jié)作用減弱。

表3 研究區(qū)域景觀格局動(dòng)態(tài)變化

3．2．2 類(lèi)型水平上格局動(dòng)態(tài)分析

（1）景觀分維數(shù)分析。分維數(shù)值界于1～2之間，趨于l，斑塊形狀越有規(guī)律，表明受干擾的程度越大。由圖1可以看出，城鎮(zhèn)村及工礦交通用地、水域及水利設(shè)施用地的平均分維數(shù)較高，在1．5以上，表明這兩類(lèi)斑塊形狀較為復(fù)雜。其中，水域及水利設(shè)施用地的平均分維數(shù)最大（1．715 9），反映其形狀最復(fù)雜、最不規(guī)則，表明它是受人類(lèi)活動(dòng)影響最小的地類(lèi)。耕地、林地、草地平均分維數(shù)分別為1．317 6，1．398 1，1．368 9，形狀相對(duì)規(guī)則，表明受人類(lèi)活動(dòng)影響較大。分維數(shù)變化最大的景觀類(lèi)型是人為影響林地，從1987年的1．271 9增加到2007年的1．524 3，由此可以看出人為影響林地斑塊形狀逐漸趨向復(fù)雜化，受人為活動(dòng)的影響較多，存在一定的環(huán)境問(wèn)題，景觀類(lèi)型不穩(wěn)定，并且比較容易造成水土流失。耕地、草地的斑塊分維數(shù)也有較大范圍的變化，但是都沒(méi)有人為影響林地變化的幅度大，表明耕地、草地的平均斑塊形狀也在逐漸發(fā)生變化。城鎮(zhèn)村及工礦交通用地是人類(lèi)活動(dòng)的集中地，受人類(lèi)活動(dòng)干擾最強(qiáng)，其形狀趨于規(guī)則的矩形，大致沿交通線(xiàn)兩側(cè)平行于交通線(xiàn)延伸，其平均分維數(shù)較大。

圖1 不同時(shí)期各景觀要素類(lèi)型格局指數(shù)動(dòng)態(tài)

（2）景觀破碎化指數(shù)分析。由圖1可以看出，耕地、林地、草地的平均斑塊密度較高，在0．3以上，表明這幾類(lèi)景觀破碎化程度較大。研究時(shí)間段內(nèi)，各景觀類(lèi)型斑塊密度指數(shù)均呈上升趨勢(shì)，研究區(qū)景觀破碎化程度加劇，表明了總體景觀組織活力較差，景觀與生物多樣性低，景觀格局簡(jiǎn)單。1987年耕地斑塊總量為225塊，2007年為935塊，斑塊數(shù)量呈現(xiàn)總體增加的趨勢(shì)，耕地、草地斑塊密度增加速度較快，最大斑塊指數(shù)大幅度下降，表明耕地和草地在向建設(shè)用地轉(zhuǎn)化過(guò)程中，原有斑塊被分割，破碎度增大。林地斑塊數(shù)量增加745，斑塊密度增大幅度最大，最大斑塊指數(shù)顯著下降，這主要是因?yàn)槿祟?lèi)毀林開(kāi)荒等活動(dòng)的影響，導(dǎo)致斑塊破碎化程度加大。隨著城市化和工業(yè)化的推進(jìn)，城鎮(zhèn)村及工礦交通用地的斑塊密度呈上升趨勢(shì)。水域、未利用土地占流域面積小，且研究時(shí)段內(nèi)變化小，均為單一的斑塊數(shù)目。所以水域、未利用土地斑塊密度發(fā)生變化較小，破碎化程度變化也較小。

（3）景觀結(jié)合度指數(shù)分析。從景觀整體來(lái)看，斑塊結(jié)合度指數(shù)均比較高，而且差異不是很大，斑塊結(jié)合度指數(shù)以耕地、林地、草地較高，接近99%，說(shuō)明它們的連通性較好，并且在景觀中占有較大比重，類(lèi)型間的物質(zhì)和能量遷移比較通暢。但是居民點(diǎn)的斑塊結(jié)合度相對(duì)較小，說(shuō)明居民點(diǎn)由于人類(lèi)活動(dòng)的阻塞，使得斑塊之間相對(duì)分散，連通性明顯低于其它景觀要素類(lèi)型。在研究時(shí)間段內(nèi)，各景觀類(lèi)型斑塊結(jié)合度指數(shù)不斷下降，說(shuō)明景觀斑塊邊緣復(fù)雜程度持續(xù)加大，斑塊分布分散且各景觀之間的相互變化更加頻繁，影響到各類(lèi)地塊之間物質(zhì)、能量和信息的流通不暢，進(jìn)而影響到區(qū)域的整個(gè)生態(tài)系統(tǒng)。1987年、2007年蒙陰縣未利用土地斑塊平均結(jié)合度指數(shù)最小，說(shuō)明該類(lèi)型相對(duì)分散，空間連接性相對(duì)較小，且發(fā)生較大變化。在大型建成區(qū)特別是城鎮(zhèn)斑塊不斷膨脹和交通道路擴(kuò)展過(guò)程中，吸收和合并的斑塊數(shù)量遠(yuǎn)遠(yuǎn)大于新增獨(dú)立小斑塊的數(shù)量，從而導(dǎo)致城鎮(zhèn)、農(nóng)村居民點(diǎn)及工礦建設(shè)用地斑塊的破碎度下降，平均斑塊面積大幅度增長(zhǎng)。

4 結(jié) 論

（1）在時(shí)間序列上，20 a間研究區(qū)景觀結(jié)構(gòu)發(fā)生了明顯變化，林地、城鎮(zhèn)村及工礦交通用地面積增加，耕地等其他用地面積減少；從變化數(shù)量看，最為顯著的是林地的增加和草地、耕地的減少，從土地利用年變化率來(lái)看，未利用土地的年變化率最大，其次是城鎮(zhèn)村及工礦交通用地。在空間序列上，各類(lèi)土地利用變化在空間上呈現(xiàn)明顯的區(qū)域分異。這種變化以林地、城鎮(zhèn)村及工礦交通用地的大面積增加和耕地等其他用地面積的減少為特征，但各種地類(lèi)之間有著復(fù)雜的相互轉(zhuǎn)換關(guān)系。

（2）通過(guò)景觀多樣性指數(shù)、結(jié)合度指數(shù)等評(píng)價(jià)指標(biāo)表明該區(qū)景觀類(lèi)型相對(duì)來(lái)說(shuō)穩(wěn)定，為該地區(qū)生物類(lèi)型的穩(wěn)定性創(chuàng)造了良好條件；景觀的分維數(shù)在1．5左右，且呈上升趨勢(shì)，表明斑塊形狀逐漸趨向復(fù)雜化，破碎化程度加劇。研究區(qū)耕地、林地和草地這幾種主要的景觀類(lèi)型之間的破碎化程度增大，結(jié)合度降低，在區(qū)域景觀要素中的優(yōu)勢(shì)不明顯，對(duì)區(qū)域整個(gè)生態(tài)系統(tǒng)的調(diào)節(jié)作用較差，說(shuō)明區(qū)域的生態(tài)環(huán)境在區(qū)域整體尺度上不利于生態(tài)系統(tǒng)效益的發(fā)揮，區(qū)域的生態(tài)系統(tǒng)脆弱性在加大，并且向惡性方向發(fā)展。

（3）從整體變化的趨勢(shì)來(lái)看，研究區(qū)景觀類(lèi)型的變化趨勢(shì)不利于區(qū)域的生態(tài)環(huán)境平衡。由于蒙陰縣所在的區(qū)域本身就屬于生態(tài)環(huán)境的脆弱地區(qū)，環(huán)境容量小，加上人為活動(dòng)的加劇，研究區(qū)的生態(tài)環(huán)境尤其難以維持其穩(wěn)定態(tài)。加大景觀格局的動(dòng)態(tài)監(jiān)測(cè)可以進(jìn)一步指導(dǎo)區(qū)域的土地利用規(guī)劃，為其生態(tài)環(huán)境保護(hù)提供參考，同時(shí)也是城市可持續(xù)發(fā)展的重要技術(shù)支撐。在發(fā)展經(jīng)濟(jì)的同時(shí)，管理和決策部門(mén)應(yīng)該引起高度重視并制定科學(xué)、合理的發(fā)展規(guī)劃，注意保護(hù)生態(tài)環(huán)境，走人地和諧、可持續(xù)發(fā)展的道路。

［1］ 王繼夏．南水北調(diào)中線(xiàn)工程水源地板凳河流域景觀格局特征變化研究［D］．西安：陜西師范大學(xué)，2007．

［2］ 傅伯杰，陳利頂，馬克明，等．景觀生態(tài)學(xué)原理及應(yīng)用［M］．北京：科學(xué)出版社，2001．

［3］ 申懷飛．基于3S的豫西黃河流域景觀格局變化研究［D］．河南開(kāi)封：河南大學(xué)，2007．

［4］ 展鵬飛．基于3S的伏牛山自然保護(hù)區(qū)景觀格局變化研究［D］．河南開(kāi)封：河南大學(xué)，2007．

［5］ 趙東娟．縣域土地利用景觀空間格局動(dòng)態(tài)分析及生態(tài)分區(qū)優(yōu)化調(diào)控研究［D］．山東泰安：山東農(nóng)業(yè)大學(xué)，2007．

［6］ 鄔建國(guó)．景觀生態(tài)學(xué)：格局、過(guò)程、尺度與等級(jí)［M］．北京：高等教育出版社，2007．

［7］ 何原榮，周青山．基于SPOT影像與Fragstats軟件的區(qū)域景觀指數(shù)提取與分析［J］．海洋測(cè)繪，2008，28（1）：86－97．

［8］ 肖鵬峰，劉順喜，馮學(xué)智，等．基于中分辨率遙感圖像的土地利用與覆被分類(lèi)系統(tǒng)構(gòu)建［J］．中國(guó)土地科學(xué)，2006，20（2）：33－38．

［9］ 李偉，賈寶全，王成，等．北京市景觀格局動(dòng)態(tài)變化［J］．東北林業(yè)大學(xué)學(xué)報(bào)，2010，38（4）：37－41．

［10］ 鄭毅，任友理，魏翔．劍川縣土地利用現(xiàn)狀多因子分析［J］．中國(guó)西部科技，2010，9（2）：6－7．

［11］ 鄭建蕊，蔣衛(wèi)國(guó)，周廷剛，等．洞庭湖區(qū)濕地景觀指數(shù)選取與格局分析［J］．長(zhǎng)江流域資源與環(huán)境，2010，19（3）：305－310．

［12］ 馮湘蘭，李際平，陳端呂．景觀空間格局指數(shù)的相關(guān)性研究［J］．中南林業(yè)調(diào)查規(guī)劃，2010，29（2）：52－55．

［13］ 張林艷，夏既勝，葉萬(wàn)輝．景觀格局分析指數(shù)選取芻論［J］．云南地理環(huán)境研究，2008，20（5）：38－43．

［14］ 高義，蘇奮振，孫曉宇，等．珠江口濱海濕地景觀格局變化分析［J］．熱帶地理，2010，30（3）：215－220．

［15］ 呂亞軍，魯建偉．基于Fragstats的徐州市景觀格局變化分析［J］．北京工業(yè)職業(yè)技術(shù)學(xué)院學(xué)報(bào)，2009，8（2）：1－4．

［16］ 陳文波，肖篤寧，李秀珍．景觀指數(shù)分類(lèi)、應(yīng)用及構(gòu)建研究［J］．應(yīng)用生態(tài)學(xué)報(bào)，2002，13（1）：121－125．