石羊河下游民勤土地利用及景觀格局動態分析

陳 杰, 楊太保, 何 毅

(蘭州大學 資源環境學院 冰川與生態地理研究所, 蘭州 730000)

石羊河下游民勤土地利用及景觀格局動態分析

陳 杰, 楊太保, 何 毅

(蘭州大學 資源環境學院 冰川與生態地理研究所, 蘭州 730000)

基于遙感(RS)和地理信息系統(GIS)平臺，以石羊河下游民勤縣1991和2011年遙感影像為基礎數據源，分析了研究區20 a來各土地利用景觀類型的變化特征，并從水土保持角度出發，利用斑塊數量(NP)、斑塊密度(PD)、最大斑塊所占景觀面積的比例(LPI)、平均斑塊面積(MPS)、香農多樣性指數(SHDI)、香農均勻度指數(SHEI)及分離度(Fi)7種景觀指數對沙漠、重鹽堿地、植被和輕鹽堿地的斑塊變化進行分析評價。結果表明：(1) 該區域以沙漠化景觀為主，且20 a來鹽堿地和植被均有所減少，沙漠化擴張；(2) 1991—2011年期間，區域整體景觀格局的景觀破碎化程度減小，景觀的異質化程度呈下降趨勢，表明該區整體景觀格局將有利于向水土保持方向發展。(3) 各土地利用類型的景觀格局動態分析都顯示植被破碎化程度加劇，人類活動對區域植被的干擾程度較大。

石羊河下游; 景觀格局; 動態分析; 土地利用

近年來，在以GIS和RS為代表的新技術支持下，對景觀動態變化進行大尺度生態監測和建模研究成為景觀生態學的新熱點[1-2]。目前國內現有的景觀動態變化研究[3-6]大多采用RS與GIS相結合的手段，其優勢在于結合RS與GIS技術對景觀格局進行動態分析，可以快速地分析出景觀斑塊的轉移和變化，揭示其變化規律，并以科學的規劃與設計為調控手段，為景觀規劃和管理提供理論依據。景觀空間格局是大小和形狀各異的景觀要素在空間上的排列，它既是景觀異質性的重要表現，又是各種生態過程在不同尺度上相互作用的結果[7-9]。用景觀指數描述景觀格局及變化，建立格局與景觀過程之間的聯系，是景觀生態學最常用的定量化分析方法[10]。目前景觀格局的研究成果主要集中于兩個方面：一是靜態的格局分析，主要探討景觀的空間異質性問題[11]；二是在相對穩定的地域空間上，探討景觀格局隨時間的變化情況，即時間異質性問題。

石羊河下游民勤地區風蝕嚴重，沙化面積和荒漠草原枯死面積逐年擴大，生物多樣性等生態環境問題日益突出，逐漸引起人們的重視。為了加強石羊河下游民勤的生態環境保護，實現該地區的可持續發展，首先需要了解該地區的土地利用現狀，認識區域景觀格局特征及其動態變化規律，分析人類活動對景觀格局及生態過程的影響，在此基礎上努力探索該區人口、經濟與資源環境的相互協調關系，尋求可持續發展的途徑[12-13]。本文以石羊河下游民勤為研究對象，在RS與GIS技術的支持下，利用景觀格局指數定量分析該區域1991—2011年的景觀格局的動態變化，揭示其景觀時空動態變化的演變規律，為實現流域可持續發展提供科學依據和決策支持。

1 材料與方法

民勤位于甘肅省西北部，河西走廊東北部，石羊河流域下游(102°45′—103°55′E，38°20′—39°10′N)。東西北三面被騰格里和巴丹吉林兩大沙漠包圍，毗鄰武威、金昌和內蒙古左右旗。該區境內地勢平坦，總面積1.6萬 km2，全縣總人口27.43萬人，常住人口24.13萬人，其土地類型主要由沙漠、鹽堿地和綠地構成，土壤主要為灰鈣土、灰漠土、綠洲灌淤土、草甸土、荒漠林土、鹽土、風沙土等類型。多年平均降水量110 mm，蒸發量2 644 mm，光照豐富，熱量充足，全年日照時數3 208 h，晝夜溫差大，平均相對濕度45%，濕潤系數0.66，平均氣溫為7.8℃，最高氣溫38.1℃，風沙大，盛行西北風，平均風速2.8 m/s，最大風速31 m/s，是我國干旱、荒漠化最嚴重的地區之一，也是我國沙塵暴四大源區之一。沙漠以每年8～10 m的速度向綠洲推進，最大處每年超過120 m，大片綠洲淪為不毛之地，沙漠區內的喬、灌木及大片草甸萎縮、枯死，固沙能力減弱，覆蓋率降低，綠洲已由過去的阻沙天塹變為沙源，水干風起，沙逼人退。如今，綠洲在水資源短缺和沙害的影響下，生態環境日益惡劣，直接威脅著綠洲的存亡。

1.1 景觀源數據

選取石羊河下游1991年和2011年Landsat TM影像數據，行列號為131/33,日期分別為1991-6-25和2011-6-24，波段組成為4,3,2 ( RGB) ，空間分辨率為30 m ，影像無云遮蓋，清晰易解譯。同期的遙感資料我們還選擇了1992-7-29和2011-7-26的遙感影像作為對比進行圖像解譯，以便解譯數據更準確。

1.2 數據處理

1.2.1 景觀類型劃分與遙感解譯 結合研究區的實際情況，參照國家《土地利用現狀調查技術規程》[14]中的分類標準，擬將研究區分為沙漠、重鹽堿地、植被和輕鹽堿地4種土地利用類型。

圖像經過幾何糾正、坐標變換和增強處理后，在ENVI中，根據解譯時建立的解譯標志和影像中地物的光譜特征來選取各類型的典型區域建立監督分類的樣本[15]，采用最大似然法進行監督分類[16]。由于在分類的過程中難免會產生一些小圖斑，需要對小圖斑進行剔除和重新分類，用聚類將周圍的小圖斑合并到大類中，用過濾剔除不符合的小圖斑；最后進行人工目視解譯。借助GIS對解譯結果進行數據合成，同時利用統計分析工具得到各種土地利用類型的定量信息，在GIS疊置分析功能的支持下，對四種土地利用類型進行變化檢測，已達到對研究區進行動態監測研究的目的，確立解譯標志和解譯精度。根據利用遙感技術開展土地資源分類調查及可能達到的規定精度，對低于應用陸地衛星TM圖像并配合采用細小地物成數抽樣方法量測達到預期精度的地類，均不能單獨列入分類系統。通過野外調繪、核實，本文遙感解譯精度為 96.7%。

1.2.2 景觀空間格局分析 石羊河下游景觀格局分析是在ArcGIS 10.0和 FRAGSTAT 3.3的聯合使用下完成的，在ArcGIS 10.0平臺支持下，首先將矢量格式的景觀類型分布圖轉換成柵格格式，然后再進行不同時段的景觀格局分析。

1.3 分析指標選取

目前已經有大量的景觀格局指數用于描述不同尺度和不同方面的景觀格局特征，但景觀格局指標的相關性和冗余性也受到生態學研究者的廣泛重視。根據前人[17]的研究經驗，本文從景觀水平和類型水平兩種尺度對石羊河下游的景觀結構及其變化進行分析。由于缺少判據斑塊核心面積和邊緣寬度的實測數據，因而不適合計算邊緣指標和核心面積指標及一部分鄰近指標；同時還排除了一部分相關性很高的指標；最后從FRAGSTATS軟件中選取了能夠確保計算精度的7個指標(表1)，景觀水平上選取斑塊數(NP)、平均斑塊面積(MPS)、香農多樣性指數(SHDI)和香農均勻度指數(SHEI)和分離度(Fi)5個指標；類型水平上選取斑塊數(NP)、斑塊密度(PD)、最大斑塊指數(LPI)、平均斑塊面積(MPS)4個指標。所選指標及其意義如表1所示。

表1 景觀格局指數

2 結果與分析

2.1 景觀要素類型面積變化

研究區總面積為940.92 km2，1991年和2011年均以沙漠和植被為主，兩種土地利用類型占總面積的近80%，鹽堿地所占比例較小(圖1)，由此反映出該區域沙漠分布和植被覆蓋率較高。

圖1 不同景觀類型面積分布及所占比例

從表2可以看出，近20 a間，研究區內各類景觀要素的面積發生了很大變化，沙漠面積增幅較大，所占總面積的比例由1991年的48.49%增加到2011年的54.60%，重鹽堿地大幅度的降低，由1991年的127.62 km2減少到2011年的89.53 km2,凈減38.09 km2,所占面積的比例由1991年的13.56%迅速下降到2011年的9.52%，植被和輕鹽堿地均有所減少。從變化數量看，在沙漠面積顯著增加的同時植被面積明顯減少，20 a間沙漠增加了57.49 km2，平均每年增加2.87 km2，同期植被減少了13.56 km2，平均每年減少0.68 km2；其原因是該區降水稀少，導致土地不能被充分的開發利用，植被減少則是由人類活動破壞造成的。鹽堿地的減少，意味著土地被開發利用或是進一步退化為沙漠；從土地利用景觀類型年變化率來看：重鹽堿地年變化率最大，平均每年減少1.49%，其次是沙漠，平均每年增加0.63%。

表2 1991-2011年不同景觀類型的變化

2.2 區域整體的景觀格局動態分析

在景觀水平上，斑塊密度是評價景觀完整性和破碎化的重要指標之一，其數量的大小可以反映景觀破碎化程度及空間異質性大小，數值越大說明破碎程度越大。平均斑塊面積也是評價景觀破碎化程度的指標之一，它與斑塊數量及斑塊的總面積有關，與破碎化呈反相關，平均斑塊面積越大，則破碎度越低。從研究區總體景觀格局來看(表3)，1991—2011年，斑塊數量(NP)、斑塊密度(PD)和最大斑塊面積百分比(LPI)分別由1991年的343，0.362 1，35.57變為2011年的252，0.271 4，31.40，逐漸減少，減少比率分別為26.53%，25.05%，11.72%；而平均斑塊面積(MPS)大幅度增加，增加比率為26.72%；香農多樣性指數和香農均勻度指數在研究時間段內分別減少了0.080 3和0.057 9，這些景觀指數的變化表明研究區域內景觀類型的豐富程度和均勻程度在下降，優勢景觀在上升，破碎度降低，對區域生態起調節作用的主要景觀類型對區域的生態調節作用在增強。20 a間研究區景觀破碎化程度逐漸減小，景觀的異質化程度呈下降趨勢，表明在經濟發展和城市化水平快速推進的背景下，人類對區域景觀的干擾程度下降。代表自然生態系統的土地覆被類型大幅度增加，從而使自然生態環境的凈化和維系能力上升，適于生物生存，物種的繁殖擴散遷移都得到了很好的保護。

表3 景觀級別指數動態變化

注：NP為斑塊個數；PD為斑塊密度；LPI為最大斑塊指數；MPS為平均斑塊面積；SHEI為香農均勻度指數；SHDI為香農多樣性指數。

2.3 景觀要素類型格局動態分析

各景觀要素的格局指數及其變化情況見圖2。沙漠是石羊河下游地區主要的景觀類型，其面積占到整個研究區面積的50%，在這20 a間，各類型的景觀格局發生了改變。

圖2 不同時期各景觀要素類型格局指數動態

(1) 斑塊數量(NP)分析。由圖2a可以看出，沙漠NP較小且有變小的趨勢，由1991年的14變為2011年的9，說明沙漠破碎化程度較低，表明受人類活動干擾較小，這與現實相符合，由于沙漠地區人類無法生活，很少受到人類的影響；植被和輕鹽堿地NP較大但也逐年減小,1991年分別為73，154，到2011年變為111，52；說明破碎化程度較高，受人類活動影響較大；從1991—2011年，植被NP增大，增率34.23%；反映出植被破碎化程度呈現增加趨勢，其主要是因為人類毀林開荒等活動的影響，導致植被斑塊破碎化程度加劇。沙漠、重鹽堿地和輕鹽堿地都是NP減少，表明其破碎化程度下降。

(2) 斑塊密度(PD)分析。由圖2b可以看出，重鹽堿地、植被和輕鹽堿地的平均斑塊密度較高，其值在1991年分別為0.108 4，0.077 6，0.163 7,都在0.1左右，2011年分別為0.084 8，0.117 7，0.055 1；表明這幾類景觀破碎化程度較大。研究時段內，植被斑塊密度指數增大了0.040 1，呈上升趨勢，植被破碎化程度加劇；重鹽堿地、沙漠和輕鹽堿地的平均斑塊密度都呈現下降趨勢，下降率為21.77%，36.24%，66.34%，其中輕鹽堿地變率最大；說明景觀破碎化程度減弱，表明了總體景觀組織活力較好，景觀與生物多樣性高，景觀格局復雜。

(3) 最大斑塊所占景觀面積的比例(LPI)分析。由圖2c可以看出，1991—2011年期間，植被的最大斑塊所占景觀面積的比例大幅度增加，從1991年的10.48變為2 011的25.96，說明受到人類活動的影響最大；重鹽堿地和輕鹽堿地都有不同程度的減少，沙漠小幅度上升，在20 a間增加了6.22。

(4) 平均斑塊面積(MPS)分析。從整體上看圖2d，除沙漠外，其它類型平均斑塊面積較低，在2.5左右；20 a間，除植被平均斑塊面積有1991年的3.54減小到2011年的2.21，減少了1.33以外，其他景觀要素都是增加，也表明植被的破壞化程度較高。

(5) 分離度(Fi)分析。其描述景觀類型中不同斑塊個體分布的分離程度。經公式計算可知，在1991年，研究區域內斑塊類型的分離度較大的是鹽堿地，其值分別為輕鹽堿地5.963 0，重鹽堿地3.295 7，分離度最小的斑塊是沙漠 0.180 6?？梢姡芯繀^域內各景觀的分離程度差異較大，鹽堿地在整個地域空間上常呈孤立的散點狀, 被其他類型斑塊所隔離，故其值較大。沙漠面積較大，且分布較集中連片，分離度較小。從變化幅度來看，植被的分離度從0.967 4增加到1.293 4，從而反映出植被的破壞程度加劇。

綜上所述，從各景觀要素的格局指數及其變化分析可知，植被破碎化程度加劇，表明在當今經濟發展和城市化水平快速推進的背景下，人類活動對區域植被的干擾程度逐漸增加，植被遭到大量破壞。

3 討論與結論

景觀格局、功能和變化是自然因素和人為因素共同作用的結果，自然因素常常是在較大時空尺度上作用于景觀。研究區域內各景觀類型面積轉移的主要去向是沙漠。民勤處于風沙沿線的前沿地帶。據氣象資料，民勤地區20世紀90 a代以來8級以上大風月平均29日，最大風力11級，年均風沙日數達139日。另據《民勤縣生態環境綜合治理情況匯報》，綠洲風沙危害性近10 a來日益加劇，流沙每年以3～4 m的速度向綠洲腹部逼近，個別地段風沙前移速度高達8～10 m。從理論上講，由于風沙化的存在，綠洲中與沙地相鄰的耕地斑塊存在著一定的風蝕潛在危險性。其相鄰斑塊邊界的長度越長，數量越多，土壤風蝕的潛在危險性越大。從遙感影像上看到1991—2011年期間綠洲邊緣地帶大量開荒造成了一些零星分布于半固定沙丘中的，沒有防護林體系保護耕地的存在，面臨沙害的危險性很大。經遙感影像解譯所得到的景觀變化結果可知：由于景觀斑塊類型的空間變化，使得各灌區植被斑塊體與沙地斑塊體相鄰邊界長度與條數，呈現出逐年增加的發展趨勢。加之在地表水逐年減少和人類不合理的土地利用雙重作用下，最終導致持續超采地下水，破壞地下水的采補平衡、引起地下水埋深大幅度下降，進而破壞植被和鹽堿地的防護體系，從而引發了植被和鹽堿地整體水平的逐年下降。

通過GIS和RS平臺，利用遙感數據結合景觀指數分析了石羊河下游民勤20 a景觀格局動態變化，結果表明：

(1) 石羊河下游地區的組成不合理，沙漠景觀所占面積的比率嚴重偏高；并且20 a來鹽堿地和植被面積均有所減少，而沙漠面積增加，沙漠化加劇。

(2) 1991—2011年，區域整體的景觀格局景觀破碎化程度逐漸減小，景觀的異質化程度呈下降趨勢，表明研究區的整體景觀格局將有利于向水土保持方向發展。

(3) 各景觀要素類型格局動態分析都顯示植被破碎化程度加劇，由于人為活動的影響，研究區的生態環境難以維持其穩定性，因此，加大景觀格局的動態監測可以進一步指導區域的土地利用規劃，并為其生態環境保護提供參考。

[1] 劉惠明,尹愛國.3S技術及其在生態學研究中的應用[J].生態科學,2002,21(1):82-85.

[2] 薛麗霞,王佐成,李永樹.基于遙感的區域景觀特征及其動態變化研究：以重慶市北部新區為例[J].國土資源遙感,2005(2):64-68.

[3] 劉常富,孫冉,李小馬,等.基于RS與GIS的沈陽城市森林景觀格局動態變化[J].東北林業大學學報,2009,37(4):13-15.

[4] 胡金龍,周志翔,金葉,等.1986—2006年桂林市區景觀格局演變分析[J].水土保持研究,2013,20(2):48-53.

[5] 陸元昌,陳敬忠,洪玲霞,等.遙感影像分類技術在森林景觀分類評價中的應用研究[J].林業科學研究,2005,18(1):31-35.

[6] 陳振雄,李際平,劉素青.基于RS和GIS的天鵝山林場森林景觀格局研究[J].林業調查規劃,2006,31(3):6-11.

[7] 陳頂利,傅伯杰.黃河三角洲地區人類活動對景觀結構的影響分析：以山東省東營市為例[J].生態學報,1996,16(4):337-344.

[8] 高小紅,王一謀,楊國靖.基于RS與GIS的榆林地區景觀格局動態變化研究[J].水土保持學報,2004,18(1):168-171.

[9] Fu B J, Lu Y H. The progress and perspectives of landscape ecology in China[J]. Progress in Physical Geography, 2006(36):319-344.

[10] 李歡,劉霞,姚孝友,等.蒙陰縣土地利用景觀格局動態分析[J].水土保持研究,2011,18(5):43-47.

[11] 王婷婷,侯淑濤,唐軍利,等.七臺河市土地利用景觀格局動態變化分析[J].水土保持研究,2012,19(3):90-93.

[12] Forman R T T, Gordon M. Landscape Ecology [M].New York: John Wiley and Sons,1986.

[13] Richards J A. Thematic mapping from multitemporal image data using the principal components transformation[J]. Remote Sensing of Environment,1984,16(1):35-46.

[14] 全國農業區劃委員會.土地利用現狀調查技術規程[S].1984.

[15] 梅安新,彭望祿,秦其明,等.遙感導論[M].北京:高等教育出版社,2001.

[16] 王靜.土地資源遙感監測與評價方法[M].北京:科學出版社,2006:60-79.

[17] 郭濼,杜世宏,趙松婷,等.黔東南地區景觀格局的動態變化及驅動力[J].山地學報.2011,29(5):543-550.

AnalysisonDynamicCharacteristicsofLandscapePatternsinMinqinofShiyangRiverDownstream

CHEN Jie, YANG Tai-bao, HE Yi

(InstituteofGlaciologyandEcogeography,CollegeofEarthandEnvironmentalSciences,LanzhouUniversity,Lanzhou730000,China)

The spatiotemporal changes of land use and cover of Minqin county in the lower reaches of Shiyang River during 1991 and 2001 were analyzed quantitatively by taking Landsat TM images as data sources and based on GIS and RS. At the same time, NP， PD，LPI ，MPS，SHDI, SHEI and Fi were used to analyze and evaluate the patch changes of desert，heavy saline soil，vegetation and saline-alkali soil. Analyses indicated that :(1) desertification landscapes counted for the largest proportion in the region, the proportion of saline-alkali soil and vegetation reduced in recent 20 years, while desertification expanded; (2) the overall regional landscape pattern and landscape fragmentation decreased and the degree of landscape heterogeneity decreased during 1991—2011, which indicated that the overall landscape pattern would be conducive to soil and water conservation; (3) various types of pattern analysis on landscape elements showed that levels of fragmentation of vegetation had been increasing, indicating that impacts of human activities on regional vegetation degree of interference were great.

Shiyang River downstream; landscape; dynamic analysis; land use

2014-01-15

：2014-03-05

國家自然科學項目“黃河上游晚更新世河流階地和冰川演化對千年尺度氣候變化的響應”(41271024)

陳杰(1990—),男,甘肅張掖人,碩士,研究方向為GIS和RS應用。E-mail:chenjie12@lzu.cn

楊太保(1962—),男,山西平陸人,博士,教授,博導,主要從事自然地理環境變化研究。E-mail：yangtb@lzu.edu.cn

F321.1P901;X171.1

：A

：1005-3409(2014)06-0251-05