近30年對(duì)干旱區(qū)內(nèi)陸河流域土地利用/土地覆被轉(zhuǎn)換的監(jiān)測(cè)

摘要：基于1989、2000和2014年的3期Landsat衛(wèi)星圖像數(shù)據(jù)，采用Intensity analysis方法，從時(shí)間間隔、類(lèi)型和轉(zhuǎn)換3個(gè)水平上分析了新疆維吾爾自治區(qū)瑪納斯河流域2個(gè)連續(xù)的時(shí)間間隔中7個(gè)土地利用類(lèi)型之間的轉(zhuǎn)換特征。結(jié)果表明，從1989-2000年時(shí)間間隔到2000-2014年時(shí)間間隔內(nèi)，流域土地利用的變化整體上由活躍趨于平緩。從第一個(gè)時(shí)間間隔（1989-2000年）到第二個(gè)時(shí)間間隔（2000-2014年），未利用地、草地和耕地的年增加面積都較大，但增加速率減小了；草地減少的變化有所減緩，而增加的變化趨于活躍；未利用地減少的變化較為平緩，而增加的變化趨于活躍；水域和建設(shè)用地的增加和減少變化都很活躍。從第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔內(nèi)，草地的年減少面積主要是向未利用地轉(zhuǎn)換，其次是耕地；兩個(gè)時(shí)間間隔內(nèi)，未利用地的增量第一個(gè)時(shí)間間隔內(nèi)主要依靠草地和冰雪地的轉(zhuǎn)入，第二個(gè)時(shí)間間隔內(nèi)主要來(lái)自草地和耕地的轉(zhuǎn)入；兩個(gè)時(shí)間間隔內(nèi)，耕地的增量主要來(lái)自草地和未利用地的轉(zhuǎn)出，并且各土地利用類(lèi)型向耕地的轉(zhuǎn)出強(qiáng)度都減弱了。

關(guān)鍵詞：土地利用；土地覆被；時(shí)空變化特征；瑪納斯河流域

中圖分類(lèi)號(hào)：F301.24 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：0439-8114（2016）21-5500-07

DOI：10.14088/j.cnki.issn0439-8114.2016.21.016

Monitoring of Land Use/Land Cover Transformations in Inland River Basin of the Arid Area in the Recent 30 Years：Taking Manasi River Basin in Xinjiang as an Example

TANG Xiao-ying，TANG Xiang-ling

（Department of Geography， College of Science，Shihezi University，Shihezi 832000，Xinjiang，China）

Abstract： On the basis of landsat thematic mapper satellite imagery of 1989， 2000， and 2014， and by using intensity analysis method， the characters of transitions among seven land categories during two consecutive time intervals were analyzed from three levels of the time interval， category， and transition. The results showed that， the change of land use was slowing down in the basin from time interval of 1989 to 2000 to time interval of 2000 to 2014. The unused land， grassland and farmland had larger size of annual gains during both time intervals， but the increasing rate decreased； the change of grassland decreased gently， and the change of increase was more active； the unused land lost gently， but it gained actively. Water and building land gained and lost very actively. From the first time interval to second time interval， grassland lost first to unused land and second to farmland in terms of the area of the annual transition； the increase of unused land was mainly from the grassland， ice and snow in the first time interval and from the grassland and farmland in the second time interval； in the two time intervals， farmland gained mainly from grassland and unused land in terms of the area of the annual transition， and the intensity of the land use types to the farmland was reduced.

Key words：land use；land cover；temporal and spatial change characteristics；Manasi river basin

土地利用/土地覆蓋變化（Land use and Land cover change，LUCC）是陸地生態(tài)系統(tǒng)變化的主要表現(xiàn)，是當(dāng)前全球變化的重要研究方向之一，也是全球環(huán)境變化研究領(lǐng)域的熱點(diǎn)和前沿課題之一[1-5]。土地利用是人類(lèi)根據(jù)一定的經(jīng)濟(jì)和社會(huì)目的進(jìn)行的土地開(kāi)發(fā)利用活動(dòng)。土地覆蓋是由于自然變化和人工建造活動(dòng)而形成的覆蓋狀況。土地利用變化通常會(huì)使其土地覆蓋發(fā)生變化，Turner等[6]將這種變化分成了轉(zhuǎn)換（Conversion）和漸變（Modification）兩種類(lèi)型，其中，轉(zhuǎn)換是指一種土地覆蓋類(lèi)型全部轉(zhuǎn)變成另一種類(lèi)型，土地屬性發(fā)生本質(zhì)變化，如森林轉(zhuǎn)變?yōu)楦兀粷u變則是同一種土地覆蓋類(lèi)型中的局部性變化，但土地屬性未發(fā)生本質(zhì)變化，如森林的間伐或疏伐等。此外，維護(hù)（Maintenance）是指讓土地覆蓋保持一定的狀態(tài)，也是人類(lèi)活動(dòng)影響土地覆蓋的一種形式[7]。土地利用/土地覆被變化是自然與人文過(guò)程交叉最為密切的問(wèn)題，體現(xiàn)了具體的人地關(guān)系[8-11]。土地利用變化可以使土地覆蓋發(fā)生變化，而土地覆蓋的變化將影響土地利用決策，進(jìn)而改變土地利用方式[12，13]。

在LUCC的研究中，通常采用轉(zhuǎn)換矩陣方法來(lái)定量分析同一區(qū)域、不同時(shí)間的土地利用類(lèi)型的相互轉(zhuǎn)換關(guān)系，常用二維表來(lái)表達(dá)，二維表可以呈現(xiàn)各個(gè)地類(lèi)間相互轉(zhuǎn)換的具體情況；也可以通過(guò)生成變化統(tǒng)計(jì)柵格圖（掩膜圖像）來(lái)呈現(xiàn)，該圖表示了前后兩幅土地分類(lèi)圖之間的地類(lèi)發(fā)生轉(zhuǎn)變的位置和類(lèi)別。但是這種傳統(tǒng)的分析方法并不能夠充分地對(duì)LUCC的特征進(jìn)行定量和系統(tǒng)的分析[14]。Aldwaik等[15]基于普拉姆1985、1991和1999年3年的土地利用圖，采用提出的Intensity analysis方法定量地分析了不同時(shí)間點(diǎn)、不同土地利用類(lèi)型的轉(zhuǎn)換過(guò)程，并與傳統(tǒng)的馬爾科夫方法進(jìn)行了比較，顯示該方法能夠提供更多的土地利用的變化信息。Mallinis等[16]基于Intensity analysis方法對(duì)雅典城郊山區(qū)的土地利用/土地覆被變化進(jìn)行了監(jiān)測(cè)，并分析了土地利用變化的驅(qū)動(dòng)因素。Huang等[14]和栗健等[17]也運(yùn)用了該方法分別對(duì)中國(guó)南方九龍河流域的土地利用變化和中國(guó)氣象干旱變化格局和過(guò)程進(jìn)行了分析，從時(shí)間間隔、干旱等級(jí)變化和轉(zhuǎn)換3個(gè)層面分析了各自的變化規(guī)律。

鑒于此，本研究基于1989、2000和2014年的Landssat衛(wèi)星影像數(shù)據(jù)，采用Aldwaik等[15]提出的Intensity analysis方法，對(duì)瑪納斯河流域1989-2014年的土地利用/土地覆蓋的時(shí)間和空間分布特征進(jìn)行了分析，以期能夠定量分析瑪納斯河流域土地利用的時(shí)空規(guī)律。

1 研究區(qū)概況與方法

1.1 研究區(qū)概況

瑪納斯河流域位于新疆天山北麓、準(zhǔn)噶爾盆地南緣，是一個(gè)典型的中國(guó)西北干旱區(qū)流域，是新疆最大的綠洲農(nóng)耕區(qū)，也是中國(guó)第四大灌溉農(nóng)業(yè)區(qū)。地理位置43°09′-46°55′N(xiāo)，84°35′-86°49′E，流域總面積約34 000 km2，地勢(shì)由東南向西北傾斜，海拔最高 5 259 m，最低168 m。由南向北依次分為山地、山前平原和沙漠三大地貌類(lèi)型區(qū)，流域?qū)賰?nèi)陸干旱區(qū)，夏季炎熱干燥，冬季寒冷多風(fēng)，年均氣溫6.8 ℃，年降水量為110～200 mm[18]，年蒸發(fā)量為1 500～2 000 mm。瑪納斯河起源于天山支脈依連哈比爾尕山的高山冰雪區(qū)，是流域內(nèi)最大的河流，向北經(jīng)紅山嘴流入綠洲平原，穿過(guò)古爾班通古特沙漠最后注入瑪納斯湖[19]。

1.2 數(shù)據(jù)來(lái)源與處理

選取的數(shù)據(jù)有瑪納斯河流域3個(gè)時(shí)期的Landsat遙感影像，分別為1989年的Landsat TM、2000年的Landsat ETM+和2014年的Landsat OLI數(shù)據(jù)（表1）；分辨率為30 m的ASTER GDEM V2版的DEM數(shù)據(jù)（數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心網(wǎng)站http：//www.gscloud.cn）。選取無(wú)云或少云的遙感影像和DEM，在ENVI5.1的支持下進(jìn)行了鑲嵌和裁剪，并將坐標(biāo)系統(tǒng)統(tǒng)一為Albers_Conic_Equal_Area投影，基準(zhǔn)面為WGS_1984，中央經(jīng)線為105°E。

本研究遙感影像分類(lèi)在ENVI5.1中完成。首先，數(shù)據(jù)預(yù)處理。分別進(jìn)行了1989、2000和2014年瑪納斯河流域遙感影像圖的NDVI和ISODATA非監(jiān)督分類(lèi)計(jì)算，在進(jìn)行ISODATA分類(lèi)時(shí)，最大分類(lèi)數(shù)量設(shè)置為10，迭代次數(shù)設(shè)置為10。將獲得的瑪納斯河流域的DEM、NDVI、ISODATA和近紅外、紅、綠波段進(jìn)行波段組合，得到一個(gè)6波段的文件。然后，選取訓(xùn)練樣本。根據(jù)瑪納斯河流域的實(shí)際土地利用情況和遙感影像的紋理信息和光譜特征，將研究區(qū)土地利用類(lèi)型劃分為冰雪地、建設(shè)用地、水域、耕地、林地、草地和未利用地7個(gè)類(lèi)型。在3期影像中，分別選取一定數(shù)量的7類(lèi)土地利用的訓(xùn)練樣本，并保證訓(xùn)練樣本的可分離度在1.8以上。最后，獲取規(guī)則和應(yīng)用規(guī)則分類(lèi)。根據(jù)選擇的訓(xùn)練樣本，基于CART算法獲取決策樹(shù)規(guī)則，進(jìn)行決策樹(shù)分類(lèi)，最終提取出瑪納斯河流域土地利用信息。使用混淆矩陣方法檢驗(yàn)遙感影像分類(lèi)的精度，1989、2000和2014年遙感影像分類(lèi)的Kappa指數(shù)分別為0.822、0.768和0.899。

1.3 Intensity analysis方法

研究主要采用了Aldwaik等[15]提出的Intensity analysis方法進(jìn)行分析。Intensity analysis方法是運(yùn)用轉(zhuǎn)移矩陣分析同一區(qū)域、不同時(shí)間點(diǎn)的土地利用和土地覆蓋的變化，分別在時(shí)間間隔、類(lèi)型和轉(zhuǎn)換3個(gè)水平上探究土地利用變化強(qiáng)度規(guī)律。根據(jù)文獻(xiàn)[15]，Intensity analysis的公式和含義見(jiàn)公式（1）-（8）。以下是這8個(gè)公式中出現(xiàn)字母的含義：J為類(lèi)型的數(shù)量；下標(biāo)i表示一個(gè)特定時(shí)間間隔的初始時(shí)間點(diǎn)的類(lèi)型；下標(biāo)j表示一個(gè)特定時(shí)間間隔的終止時(shí)間點(diǎn)的類(lèi)型；下標(biāo)m表示向其他土地利用類(lèi)型轉(zhuǎn)換的類(lèi)型；下標(biāo)n表示由其他土地利用類(lèi)型轉(zhuǎn)換而來(lái)的類(lèi)型；T為時(shí)間點(diǎn)的數(shù)量；Yt為時(shí)間點(diǎn)t的年份；下標(biāo)t表示時(shí)間間隔[Yt，Yt+1]的初始時(shí)間點(diǎn)，范圍是[1，T-1]；Ctij為從時(shí)間點(diǎn)Yt的i類(lèi)型轉(zhuǎn)換至?xí)r間點(diǎn)Yt+1的j類(lèi)型的象元個(gè)數(shù)。

1.3.1 時(shí)間間隔水平 時(shí)間間隔水平的分析，可以說(shuō)明各個(gè)時(shí)間間隔內(nèi)，各土地利用類(lèi)型整體上的年變化速率是相對(duì)快還是慢。公式（1）中，定義了研究區(qū)景觀從上一個(gè)時(shí)間間隔到下一個(gè)時(shí)間間隔內(nèi)年平均變化率，U是Intensity analysis均值線的值。公式（2）中，St定義為某個(gè)時(shí)間間隔內(nèi)景觀的年變化率，是時(shí)間段[Yt，Yt+1]的年變化強(qiáng)度。

U=×100%=×100% （1）

St=×100%=×100% （2）

1.3.2 類(lèi)型水平 類(lèi)型水平上的分析結(jié)果可以說(shuō)明各時(shí)間間隔內(nèi)各土地利用類(lèi)型的變化是相對(duì)平緩還是活躍。公式（3）中，Gtj為[Yt，Yt+1]間隔中j類(lèi)型的年增加強(qiáng)度。公式（4）中，Lti為[Yt，Yt+1]間隔中i類(lèi)型的年減少?gòu)?qiáng)度。

Gtj=×100%=×100% （3）

Lti=×100%=×100% （4）

1.3.3 轉(zhuǎn)換水平 表示各個(gè)土地利用類(lèi)型在各個(gè)時(shí)間間隔內(nèi)與其他土地利用類(lèi)型相互轉(zhuǎn)換的強(qiáng)度。公式（5）中，Wtn為[Yt，Yt+1]間隔中從時(shí)間點(diǎn)Yt的非n土地利用類(lèi)型轉(zhuǎn)換至n類(lèi)型的年平均強(qiáng)度。公式（6）中，Rtin為[Yt，Yt+1]間隔中從i土地利用類(lèi)型轉(zhuǎn)換至n（n≠i）類(lèi)型的年轉(zhuǎn)換強(qiáng)度。公式（7）中，Vtm為[Yt，Yt+1]間隔中從末時(shí)間點(diǎn)Yt+1的m土地利用類(lèi)型轉(zhuǎn)換至其他所有類(lèi)型的年轉(zhuǎn)換的平均強(qiáng)度。公式（8）中，Qtmj為[Yt，Yt+1]間隔中從m土地利用類(lèi)型轉(zhuǎn)換至j（j≠m）類(lèi)型的年轉(zhuǎn)換強(qiáng)度。

Wtn=×100%=×100% （5）

Rtin=×100%=×100% （6）

Vtm=×100%=×100% （7）

Qtmj=×100%=×100% （8）

1.4 Intensity analysis的數(shù)據(jù)處理

為了進(jìn)行Intensity analysis，對(duì)研究區(qū)1989-2000年和2000-2014年這兩個(gè)時(shí)間間隔研究區(qū)的面積（像元數(shù)目）計(jì)算了轉(zhuǎn)移矩陣，具體見(jiàn)表2和表3。

2 結(jié)果與分析

2.1 瑪納斯河流域年際土地利用變化

圖1顯示了瑪納斯河流域1989、2000和2014年3個(gè)年份的土地利用情況。1989-2014年，瑪納斯河流域土地利用狀況發(fā)生了較大的變化，整體上未利用地、草地、冰雪地、林地和水域這5類(lèi)土地利用類(lèi)型呈減少趨勢(shì)，而耕地和建設(shè)用地表現(xiàn)出較明顯的增加趨勢(shì)。

研究區(qū)的總面積約為34 140.7 km2，1989-2000年，草地、冰雪地和林地的面積都呈減少的趨勢(shì)，而未利用地、耕地、水域和建設(shè)用地的面積呈現(xiàn)增加的趨勢(shì)。期間，草地面積從11 192.5 km2減少到9 992.0 km2，占流域面積的比例從32.8%降低到29.3%；冰雪地減少的幅度較大，12年間面積減少了925.6 km2，所占總面積比例從5.8%降低為3.1%；林地面積由1 048.3 km2降為479.2 km2，面積所占比例從3.1%降為1.4%。未利用地在此期間增幅較大，面積從14 885.3 km2增加到16 646.6 km2，占研究區(qū)總面積的比例由43.6%增加到48.8%；水域面積也呈增加的趨勢(shì)，面積增加了184.3 km2；耕地面積增加了717.6 km2，占流域總面積的比例從14.1%升至16.2%；建設(shè)用地增幅顯著，面積由21.4 km2增至53.5 km2。

2000-2014年，未利用地和水域面積有一定減少，其他5類(lèi)土地利用面積都呈增加趨勢(shì)。耕地增量最大，面積所占總面積的比例由16.2%升至22.8%，增加了2 281.6 km2；建設(shè)用地和冰雪地也有一定的增幅，面積分別增加了51.6 km2和630.5 km2；未利用地面積降幅顯著，這15年間面積減少了3 274.6 km2，所占總面積的比例由48.8%降為39.2%；水域也呈現(xiàn)較為明顯的下降趨勢(shì)，面積減少了222.0 km2。

2.2 瑪納斯河流域土地利用轉(zhuǎn)換時(shí)間強(qiáng)度分析

圖2為時(shí)間間隔水平的強(qiáng)度分析結(jié)果，時(shí)間間隔水平上的分析能夠說(shuō)明不同時(shí)間間隔的變化規(guī)模和速度。圖2中，垂直軸表示時(shí)間間隔，坐標(biāo)原點(diǎn)的左側(cè)水平橫條表示每個(gè)時(shí)間間隔的土地利用類(lèi)型轉(zhuǎn)換的總面積占總面積的比例，坐標(biāo)原點(diǎn)的右側(cè)水平橫條表示每個(gè)時(shí)間間隔內(nèi)年的轉(zhuǎn)換強(qiáng)度，用年變化面積所占總面積的百分比表示。水平的虛線即均值線，表示年平均變化強(qiáng)度。如果水平橫條超過(guò)均值線，表示該時(shí)間間隔內(nèi)土地利用類(lèi)型轉(zhuǎn)換相對(duì)較快，反之則相對(duì)較為緩慢。

間隔水平上，圖2左側(cè)呈現(xiàn)了第一個(gè)時(shí)間間隔內(nèi)（1989-2000年）的土地利用轉(zhuǎn)換規(guī)模大于第二個(gè)時(shí)間間隔（2000-2014年）的規(guī)模，而且第一個(gè)時(shí)間間隔的持續(xù)期比第二個(gè)時(shí)間間隔少3年，圖2右側(cè)顯示了第一個(gè)時(shí)間間隔內(nèi)的年變化速率實(shí)際上比第二個(gè)時(shí)間間隔快0.73個(gè)百分點(diǎn)，這就解釋了研究區(qū)土地利用轉(zhuǎn)換規(guī)模第一個(gè)時(shí)間間隔比第二個(gè)時(shí)間間隔大的現(xiàn)象。該水平上的分析僅表明了兩個(gè)時(shí)間間隔內(nèi)，瑪納斯河流域土地利用的總體轉(zhuǎn)換速率減緩了。

2.3 瑪納斯河流域土地利用類(lèi)型轉(zhuǎn)換強(qiáng)度

圖3顯示的是土地利用類(lèi)型水平上的強(qiáng)度分析結(jié)果，左邊的橫條表示瑪納斯河流域某一種土地利用類(lèi)型年增加或減少的面積，右邊的橫條表示的是某一種土地利用類(lèi)型年增加或減少的強(qiáng)度。

對(duì)圖3a和圖3b進(jìn)行對(duì)比可以看出，整體上第一個(gè)時(shí)間間隔內(nèi)各土地利用類(lèi)型的年減少面積和增加面積都較第二個(gè)時(shí)間間隔更大，這說(shuō)明研究區(qū)各土地利用類(lèi)型的減少和增加變化從第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔內(nèi)明顯地減速了。由圖3還可以看出，兩個(gè)時(shí)間間隔中，未利用地年增量都最大，并且年增量從第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔速度減緩了；草地和耕地的增量分別排第二位和第三位，與未利用地有著相似的變化規(guī)律。就減少變化具體而言，第一個(gè)時(shí)間間隔內(nèi)，冰雪地、未利用地、草地、耕地和林地面積都有較大規(guī)模的減少；而第二個(gè)時(shí)間間隔內(nèi)，各土地利用類(lèi)型的年減少面積都很小。從1989-2014年，草地、耕地和未利用地的年變化面積位列前三位并依次遞減，而在1989、2000和2014年這3個(gè)時(shí)間點(diǎn)上，這3種土地利用類(lèi)型占總面積的比例也是最大的三類(lèi)，未利用地在各個(gè)時(shí)間點(diǎn)上所占面積最大，草地次之，耕地居第三位，在某種程度上表明，在土地利用類(lèi)型轉(zhuǎn)換過(guò)程中，大面積的土地利用類(lèi)型會(huì)有大面積的轉(zhuǎn)換量，也可能是由土地利用類(lèi)型的活躍程度所致。

圖3的右側(cè)顯示了水域在兩個(gè)時(shí)間間隔內(nèi)增加和減少的橫條都超過(guò)了均值線，這表明水域在1989-2014年增加和減少的變化都是非常活躍的，由于水域占研究區(qū)總面積的比例較小，因此在兩個(gè)時(shí)間間隔內(nèi)水域增加和減少的面積都很小。建設(shè)用地的增加強(qiáng)度僅次于水域，但占景觀總面積的比例很小，所以變化的規(guī)模也很小。冰雪地有著相似的變化規(guī)律。未利用地減少的橫條都未達(dá)到均值線，表明在兩個(gè)時(shí)間間隔中，未利用地減少的強(qiáng)度小于景觀水平。而未利用地增加的橫條第一個(gè)時(shí)間間隔內(nèi)接近于均值線，而在第二個(gè)時(shí)間間隔中超出了均值線，這表明未利用地從第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔中增加的變化趨于活躍。草地減少的變化從第一個(gè)時(shí)間間隔到第二時(shí)間間隔有所減弱，而草地增加的橫條兩個(gè)時(shí)間間隔內(nèi)都延伸至均值線右側(cè)，這表明草地比其他景觀增加的強(qiáng)度更強(qiáng)。耕地的減少變化都相對(duì)平緩，增加變化趨于活躍，由于耕地在總面積的比例較大，所以耕地在兩個(gè)時(shí)間間隔中都有著大的變化規(guī)模。

2.4 瑪納斯河流域土地利用轉(zhuǎn)換強(qiáng)度分析

圖4顯示了草地、未利用地和耕地增加或者減少的轉(zhuǎn)換強(qiáng)度分析結(jié)果。圖4左側(cè)的水平橫條顯示了某種土地利用類(lèi)型和其余類(lèi)型之間轉(zhuǎn)換的年度變化面積，右側(cè)的水平橫條顯示了某種土地利用類(lèi)型和其余類(lèi)型之間的轉(zhuǎn)換強(qiáng)度，垂直的虛線表示年變化強(qiáng)度的平均值。圖4a和圖4b給出了草地減少的轉(zhuǎn)換強(qiáng)度分析，可以看出，草地減少時(shí)，轉(zhuǎn)出為未利用地是最大的轉(zhuǎn)換，其次是耕地。再分析圖4a和圖4b的右側(cè)，未利用地的右側(cè)橫條第一個(gè)時(shí)間間隔內(nèi)極接近于均值線左側(cè)，第二個(gè)時(shí)間間隔內(nèi)略微超出均值線右側(cè)，這意味著第一個(gè)時(shí)間間隔中，草地向未利用地的轉(zhuǎn)換較為平緩，第二個(gè)時(shí)間間隔中，草地向未利用地的轉(zhuǎn)換又有所加強(qiáng)；而耕地的右側(cè)橫條第一個(gè)時(shí)間間隔內(nèi)超出均值線較多，但第二個(gè)時(shí)間間隔內(nèi)距離均值線左側(cè)還有一段距離，這表明第一個(gè)時(shí)間間隔中，草地主要向耕地轉(zhuǎn)出，而到第二個(gè)時(shí)間間隔中，草地向耕地的轉(zhuǎn)出的強(qiáng)度又平緩了。由此可知，草地在1989-2014年整體上一直呈現(xiàn)轉(zhuǎn)出的趨勢(shì)，但其間轉(zhuǎn)換的方向與速度卻是在變化的。

圖4c和圖4d給出了1989-2000年和2000-2014年兩個(gè)時(shí)間間隔內(nèi)未利用地轉(zhuǎn)換水平上的增加情況。圖4c和圖4d左側(cè)顯示，未利用地的增量主要來(lái)自草地、耕地和冰雪地這3種土地利用類(lèi)型的轉(zhuǎn)入，第一個(gè)時(shí)間間隔內(nèi)未利用地的增量主要來(lái)自草地和冰雪地，而第二個(gè)時(shí)間間隔內(nèi)未利用地的增量主要來(lái)自草地和耕地的轉(zhuǎn)入。再分析圖4c和圖4d的右側(cè)，兩個(gè)時(shí)間間隔內(nèi)冰雪地轉(zhuǎn)出到未利用地的強(qiáng)度是最大的，草地轉(zhuǎn)出到未利用地的強(qiáng)度由第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔有明顯的減弱，而耕地和水域的轉(zhuǎn)出強(qiáng)度都有明顯的增強(qiáng)。但由于冰雪地和水域在總面積中所占的百分比很小，所以在研究期間它們的變化面積也較小。

圖4e、圖4f分別代表1989-2000年和2000-2014年時(shí)間間隔各土地利用類(lèi)型轉(zhuǎn)入耕地的情況。圖4e、圖4f左側(cè)表明，耕地的增量主要來(lái)自草地和未利用地的轉(zhuǎn)出。圖4e、圖4f右側(cè)表明，從第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔各土地利用類(lèi)型轉(zhuǎn)出到耕地的轉(zhuǎn)換強(qiáng)度整體上減弱了，草地、水域和建設(shè)用地轉(zhuǎn)出到耕地的過(guò)程較活躍，而未利用地等轉(zhuǎn)出到耕地的過(guò)程比較緩慢。

3 結(jié)論

基于遙感影像、ENVI軟件、GIS軟件和Intensity analysis方法，分析了新疆瑪納斯河流域1989、2000和2014年3個(gè)時(shí)間點(diǎn)間的LUCC，得到了以下結(jié)論。

1）1989-2014年，瑪納斯河流域土地利用狀況發(fā)生了較大的變化，整體上未利用地、草地、冰雪地、林地和水域這5類(lèi)土地利用類(lèi)型呈減少趨勢(shì)，而耕地和建設(shè)用地表現(xiàn)出較明顯的增加趨勢(shì)。從1989-2000年時(shí)間間隔到2000-2014年時(shí)間間隔內(nèi)，瑪納斯河流域土地利用的變化由活躍趨于平緩，總的變化面積也由26.36%下降至23.40%。

2）從第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔內(nèi)，所有土地利用類(lèi)型的增加和減少變化都明顯減慢。1989-2000年和2000-2014年兩個(gè)時(shí)間間隔內(nèi)，未利用地、草地和耕地的年增加面積都較大，而且從第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔內(nèi)，增加速率減小了。水域和建設(shè)用地的增加和減少變化都較為活躍，由于在景觀中所占面積比例小，面積增量都不大。未利用地的減少變化較為平緩，而增加變化是趨于活躍的。草地的減少變化比較平緩，而增加變化兩個(gè)時(shí)間間隔內(nèi)都很活躍。耕地的減少變化兩個(gè)時(shí)間間隔中都較平緩，而增加變化從第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔內(nèi)由活躍轉(zhuǎn)向于平緩。

3）轉(zhuǎn)換層次的分析表明，從第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔內(nèi)，草地減少的主要方向是未利用地，其次是耕地。就轉(zhuǎn)換強(qiáng)度的分析表明，草地的減少第一個(gè)時(shí)間間隔內(nèi)傾向于耕地，而第二個(gè)時(shí)間間隔內(nèi)轉(zhuǎn)向于未利用地。未利用地的增量第一個(gè)時(shí)間間隔內(nèi)主要來(lái)自草地和冰雪地，而第二個(gè)時(shí)間間隔內(nèi)主要來(lái)自草地和耕地的轉(zhuǎn)入。兩個(gè)時(shí)間間隔內(nèi)，未利用地的增加都傾向于來(lái)自冰雪地的轉(zhuǎn)入，但由于冰雪地的面積占景觀的比例小，所以并沒(méi)有大規(guī)模的向未利用地轉(zhuǎn)入。未利用地的增加從第一個(gè)時(shí)間間隔到第二個(gè)時(shí)間間隔內(nèi)已經(jīng)從傾向于草地轉(zhuǎn)向避免草地。研究期間，耕地的增量主要來(lái)自草地和未利用地的轉(zhuǎn)出，并且各土地利用類(lèi)型向耕地的轉(zhuǎn)出強(qiáng)度都減弱了。

參考文獻(xiàn)：

[1] 唐華俊，陳佑啟，邱建軍，等.中國(guó)土地利用/土地覆蓋變化研究[M].北京：中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社，2004.

[2] 王秀蘭，包玉海.土地利用動(dòng)態(tài)變化研究方法探討[J].地理科學(xué)進(jìn)展，1999，18（1）：81-87.

[3] 史培軍，宮 鵬，李曉兵，等.土地利用/覆蓋變化研究的方法與實(shí)踐[M].北京：科學(xué)出版社，2000.

[4] 張曉明，余新曉，武思宏，等.黃土丘陵溝壑區(qū)典型流域土地利用/土地覆被變化水文動(dòng)態(tài)響應(yīng)[J].生態(tài)學(xué)報(bào)，2007，27（2）：414-423.

[5] 婁和震，楊勝天，周秋文，等.延河流域2000-2010年土地利用/覆蓋變化及驅(qū)動(dòng)力分析[J].干旱區(qū)資源與環(huán)境，2014，28（4）：15-21.

[6] TURNER B L， MEYER W B， SKOLE D L. Global land-use/land-cover change：Towards an integrated study[J].Ambio，1994， 23（1）：91-95.

[7] 擺萬(wàn)奇，柏書(shū)琴.土地利用和覆蓋變化在全球變化研究中的地位與作用[J].地域研究與開(kāi)發(fā)，1999，18（4）：13-16.

[8] 馬 禮，蘇海霞.河北沽源縣土地利用動(dòng)態(tài)變化研究[J].干旱區(qū)資源與環(huán)境，2010，24（1）：131-136.

[9] 馬松增，史明昌，楊貴森，等.基于GIS的土地利用時(shí)空動(dòng)態(tài)變化分析——以塔里木盆地農(nóng)墾區(qū)為例[J].水土保持研究，2013， 20（1）：177-181.

[10] 李屹峰，羅躍初，劉 綱，等.土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)功能的影響——以密云水庫(kù)流域?yàn)槔齕J].生態(tài)學(xué)報(bào)，2013，33（3）：726-736.

[11] 王淵剛，羅格平，趙樹(shù)斌，等.新疆耕地變化對(duì)區(qū)域碳平衡的影響[J].地理學(xué)報(bào)，2014，69（1）：110-120.

[12] 陳佑啟.試論城鄉(xiāng)交錯(cuò)帶土地利用的形成演變機(jī)制[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃，2000，21（5）：22-25.

[13] 邱炳文，陳崇成.基于多目標(biāo)決策和CA模型的土地利用變化預(yù)測(cè)模型及其應(yīng)用[J].地理學(xué)報(bào)，2008，63（2）：165-174.

[14] HUANG J L，PONTIUS JR R G，LI Q S，et al. Use of intensity analysis to link patterns with processes of land change from 1986 to 2007 in a coastal watershed of southeast China[J].Applied Geography，2012，34：371-384.

[15] ALDWAIK S Z，PONTIUS JR R G. Intensity analysis to unify measurements of size and stationarity of land changes by interval，category，and transition[J].Landscape and Urban Planning，2012，106（1）：103-114.

[16] MALLINISA G，KOUTSIASB N，ARIANOUTSOUC M. Monitoring land use/land cover transformations from 1945 to 2007 in two peri-urban mountainous areas of Athens metropolitan area，Greece[J].The Science of the Total Environment，2014， 490：262-278.

[17] 栗 健，岳耀杰，潘紅梅，等.中國(guó)1961-2010年氣象干旱的時(shí)空規(guī)律——基于SPEI和Intensity analysis方法的研究[J].災(zāi)害學(xué)，2014，29（4）：176-182.

[18] 馮異星，羅格平，周德成，等.近50a土地利用變化對(duì)干旱區(qū)典型流域景觀格局的影響——以新疆瑪納斯河流域?yàn)槔齕J].生態(tài)學(xué)報(bào)，2010，30（16）：4295-4305.

[19] 劉金巍，靳甜甜，劉國(guó)華，等.新疆瑪納斯河流域2000-2010年土地利用/覆蓋變化及影響因素[J].生態(tài)學(xué)報(bào)，2014，34（12）：3211-3223.