新疆葉爾羌河流域胡楊林時空格局特征

李華林,白林燕,馮建中,高華端,冉啟云,于 濤,高旺旺

1 中國科學院遙感與數字地球研究所數字地球重點實驗室,北京 100094 2 中國農業科學院農業信息研究所,北京 100081 3 貴州大學林學院,貴陽 550025

葉爾羌河是塔里木盆地西部最大的河流,孕育著新疆最大的綠洲之一—葉爾羌綠洲,是保障新疆塔里木河生態系統健康運行的重要源流之一,對維護流域乃至整個新疆生態環境安全、促進新疆經濟發展具有重要的作用[1- 2]。葉爾羌河流域地處塔里木河上游地區,屬于典型極端干旱區[3],生態系統脆弱,極易受到氣候變化和人類活動的影響。近年來,在自然和強烈的人為干擾下,流域內河道斷流現象頻繁發生[4],以胡楊林為代表的荒漠河岸林日益衰敗,減弱了其防風固沙、調節氣候、維持綠洲生態可持續發展等作用[5-6]。因此,準確掌握葉爾羌河胡楊林時空格局變化特征,對保護和恢復葉爾羌河流域生態環境、保障“一帶一路”生態文明建設順利實施有著十分重要的意義。

胡楊(PopuluseuphraticaOliv)是楊柳科楊屬中最古老、最原始的木本植物,它具有抗寒、抗熱、抗干旱、抗風沙、耐鹽堿等優良特性,是新疆干旱荒漠區天然林的建群樹種和珍貴的抗逆種質資源,有著重要的生態作用[7]。研究表明,水分條件仍然是制約胡楊林生長發育的主要因素,決定其群落結構和功能及植物個體形態和生理特征[8]。胡楊根系是吸收土壤水分和養分的關鍵,根系的分布決定著水資源的利用效率[9]；水分分布和鹽分條件影響著胡楊母樹和種子育苗的繁殖模式,對胡楊林的更新產生深刻的影響[10]；生態輸水能有效地改善胡楊林的局部生態環境[11-12]。目前,對新疆胡楊林流域尺度的研究主要集中在塔里木河干流。例如：崔彥軍等[13]、胡新振等[14]、金文斌等[15]、萬紅梅等[16]利用遙感方法提取塔里木河干流胡楊林信息特征；黃青等[17]、王亮[18]、趙虎等[19]通過遙感分析方法對塔里木河干流土地利用景觀格局進行評價與分析；朱緒超等[20]、陳海燕等[21]、張繪芳等[22]利用遙感方法研究塔里木河干流植被群落結構特征及變化趨勢；王家強等[23]采用高光譜技術研究塔里木河上游胡楊和灰葉胡楊葉片氮素含量。而對塔里木河上游葉爾羌河流域胡楊林研究較少,只涉及到胡楊林的碳、氮含量特征、胡楊林演變規律以及胡楊林生態保護對策等方面,對利用遙感方法監測葉爾羌河流域胡楊林時空格局的研究極少。如：王晶等[24]研究葉爾羌河流域荒漠河岸林胡楊碳、氮含量特征,結果表明不同季節胡楊葉片C、N含量不同；李靖等[25]采用系統動力學和隸屬度結合的方法建立了流域生態承載力模型,確定了葉爾羌河流域的最優發展方案；陳超群[26]探討了葉爾羌河流域胡楊林的演變規律及生態保護政策；程仲雷等[27]、胡林金等[28]、凌紅波等[29]、孫本國等[30]利用水文資料對葉爾羌河徑流量與降水量變化進行了分析研究。鑒于此,本文基于1995—2015年Landsat TM/OLI長時間序列土地利用遙感監測數據,通過土地利用程度變化來揭示區域土地利用狀況和變化趨勢；利用景觀指數分析胡楊林時空格局動態變化特征,旨在揭示葉爾羌河流域胡楊林時空格局動態變化規律及其對氣候變化和人類活動的響應,為保護葉爾羌綠洲和恢復葉爾羌河流域生態系統提供科學依據。

1 研究區概況

葉爾羌河流域位于新疆維吾爾族自治區西南部,塔里木盆地西部(圖1),地理坐標為74°28′—80°54′E,34°50′—40°31′N[5]。葉爾羌河發源于喀喇昆侖山,全長1097 km,多年平均徑流量為66.3億m3,且年際和年內流量分布不均,6—8月來水量較大,徑流主要靠冰川融雪和降雨補給,最終匯入塔里木河。在帕米爾高原、喀喇昆侖山和南天山的阻礙作用下,使大西洋和印度洋暖、濕氣流無法入境,導致區域常年干旱少雨,空氣干燥,多年平均降雨量為30—60 mm,屬于典型干旱大陸性氣候。流域內光熱資源豐富,日照時間長,蒸發強烈,年均氣溫為11.9—12.1℃,最高溫為39.4℃,最低氣溫為-25℃[31]。土壤主要以沙壤土、板結鹽堿土為主,pH值為7.5—8.5[32]。

圖1 研究區域位置Fig.1 The study area

2 材料與方法

2.1 數據來源

本文數據主要包括土地利用數據、社會經濟數據和氣象數據。土地利用數據為北京宇世藍圖信息技術有限公司基于Landsat TM/OLI 遙感影像數據,通過人工目視解譯方法生產的1995年、2005年、2015年3期柵格產品數據,其空間分辨率為30 m,土地利用類型劃分為耕地、林地、草地、水域、居民用地和未利用地6個一級利用類型。本研究采用葉爾羌河流域1995年、2005年Landsat5 TM(30 m)和2015 年Landsat5 OLI(30 m)7、8、9月的遙感影像數據(http://www.gscloud.cn/),基于RULEGEN決策樹分類方法[33]并結合人工目視解譯提取天然胡楊林,并更新至已有的土地利用產品數據中,最終得到包括耕地、胡楊林、草地、其他林地、水域、居民地、未利用地7個類型的土地利用數據。每期數據選擇3個固定驗證樣區進行精度驗證,1995年使用Landsat5 TM(30 m)數據,2005年使用Terra ASTER(15 m)數據(https://glovis.usgs.gov/),2015年使用GF1 PMS1(8 m)數據(http://210.72.27.21:8790/SNFFWeb/WebUI/HomePage.jsp)通過人工目視解譯方法得到驗證樣區土地利用分類數據,樣區土地利用數據與基于RULEGEN的決策樹分類的土地利用數據構建混淆矩陣進行精度驗證評價[34](表1)。社會經濟數據來源于《新疆統計年鑒》(1995—2016年)各地、州、市、縣(市)的耕地面積、人口數、地區生產總值以及自然資源狀況等數據。氣象數據來自于中國氣象數據網(http//cdc.cma.gov.com)“中國地面氣候資料年值數據集”,數據集包括氣壓、氣溫、水汽壓、相對濕度、風向風速、降水、蒸發量和日照等8個要素的歷年年值數據,本文選取研究區內氣象站點1995—2015年21期年平均氣溫和年降水量數據。

表1 土地利用分類精度驗證

2.2 研究方法

圖2 土地利用轉移過程Fig.2 Land use transfer processT：初期,The beginning of land use types；T+1：末期,The end of land use types；A、B、C、D、E、F：不同土地利用類型,Different land use types

本文基于1995—2015年Landsat TM/OLI長時間序列土地利用遙感監測數據,利用土地利用程度變化指數、轉移矩陣以及景觀指數對新疆葉爾羌河流域胡楊林時空變化進行監測與分析。(1)野外樣方調查：采用樣帶調查方法對葉爾羌河流域天然胡楊林進行野外調查,在葉爾羌河中、下游天然胡楊林分布區,選擇具有代表性的6個調查地點(圖1),在調查地點垂直河岸不同距離(250 m,500 m,1000 m,1500m)設立20 m×20 m調查樣方進行群落調查。(2)轉移矩陣(圖2)：土地利用轉移矩陣反映某一區域某一時段初期(T)和期末(T+1)各土地利用類型之間相互轉化的動態過程信息[35],(3)土地利用指數(表2)：本文利用單一土地利用動態度、單一土地利用空間動態度、單一土地利用變化狀態指數評價土地利用的變化速率和現狀。單一土地利用動態度表征某土地利用類型數量變化；單一土地利用空間動態度表征某類型在空間變化及其變化頻率；單一土地利用變化狀態指數表征某類型的趨勢和狀態。(4)景觀指數(表2)：以充分表征景觀組成、結構和反映景觀全局及各類型變化為原則,結合研究目標和尺度選擇具有相同生態學意義的景觀指數,最終選取斑塊密度(Patch density)、最大斑塊指數(Largest patch index)、散布與并列指數(Interspersion and juxtaposition index)、結合度指數(Patch cohesion index)共4個類型水平景觀指數；選取景觀形狀指數(Landscape shape index)、蔓延度指數(Contagion index)、分離度指數(Splitting index)、Shannon多樣性指數(Shannon′s diversity index)、Shannon均勻度指數(Shannon′s evenness index)共5個景觀水平指數[36]。

表2 土地利用與景觀分析

3 結果與分析

3.1 區域土地利用時—空變化

3.1.1土地利用變化狀況

1995—2015年間,葉爾羌河流域土地利用類型面積發生了顯著變化(圖3),流域內草地、水域、未利用地占主要類型,其總面積在三期土地利用中分別占91.03%、88.90%和92.60%。耕地和居民地面積呈上升趨勢,比重由1995年的5.73%增加到2015年的9.72%,面積增加3810.14 km2；胡楊林面積呈下降趨勢,比重由1995年的2.01%下降至2015年的1.73%,面積減少263.90 km2。1995—2005年,其他林地(除胡楊林外)和水域面積分別增加306.57 km2和2435.09 km2,而2005—2015年,其他林地和水域面積分別減少622.85 km2和4370.36 km2。1995—2005年,草地面積減少1942.60 km2,而2005—2015年,草地面積增加1268.95 km2。

圖3 1995—2015年葉爾羌河流域土地利用類型Fig.3 Land use types of the Yarkant River Basin from 1995 to 2015

根據野外調查和圖3可知：胡楊林沿主河道、支流以走廊式分布在葉爾羌河中、下游地區,而上游地區分布較少。胡楊林主要分布在四十八團至三河口河道兩側、低階地、河漫灘以及沖積平原等區域。四十八團至永安壩段受小海子、永安壩水庫和葉爾羌河地下水的補給,其土壤濕潤,養分充足,胡楊林長勢旺盛,植物群落復雜,具有喬、灌、草三層結構；永安壩至夏河林場段胡楊林受攔截洪水漫灌,胡楊林長勢較好,幼成林較多,植物群落單一；夏河林場至三河口段地下水較深,胡楊林長勢較差,林地稀疏,無幼成林,呈散生狀態,且胡楊林向河道逐漸縮減,受沙漠侵襲較嚴重。總體上,隨著遠離河道主軸線,植物種類逐漸減少和植被結構趨于簡單,胡楊林個體長勢逐漸下降和衰敗。葉爾羌河下游水資源缺乏,胡楊林逐漸衰退,植被生態系統活力、結構自我調節能力和恢復能力趨向退化,加上人為干預破壞系統承載力,導致胡楊林大面積死亡。

3.1.2土地利用動態變化

根據葉爾羌河流域1995—2015年土地利用動態變化情況(圖4)可知,近20年耕地和居民地土地利用動態度和變化狀態指數較高,表明該土地利用類型變化幅度較大,增速較快,面積有持續增加的趨勢。1995—2005年時段內耕地土地利用動態度和變化狀態指數均高于2005—2015年時段,表明耕地面積增幅降低和增速有減緩的趨勢,而居民地面積增速加快。其他林地土地利用動態度由從1995—2005年的2.61%降低到2005—2015年的-4.20%,表明其他林地在1995—2005年內增速和增幅較大,而在2005—2015年內減少,在研究時間序列內變化劇烈。從土地利用空間動態度來看,近20年其他林地和居民地土地利用空間動態度較大,表明該時期其他林地和居民地轉移面積較多,空間轉移頻繁；而水域、未利用地及草地空間動態度較少,表明此三種土地利用類型空間變化頻率小。

胡楊林土地利用動態度從1995—2005年的1.04%降低到2005—2015年的0.38%；土地利用變化狀態指數從1995—2005年的-0.10升高到2005—2015年的-0.06,表明胡楊林在研究時段內面積持續減少,并幅度減少和速度有減緩的趨勢。從空間上來看,胡楊林土地利用空間動態度從1995—2005年的10.35%降低到2005—2015年的6.77%,表明胡楊林在1995—2005年內土地利用類型轉換頻繁,空間變化頻率較高,而2005—2015年有減緩的趨勢。1995—2005年胡楊林的減少速度較快,而且一直呈現持續減少的趨勢,這與人口對耕地的需求逐漸增加、人們對資源的合理開發意識較低以及對胡楊林的保護意識較差有關,最終導致胡楊林面積減少。而2005—2015年胡楊林減少的速度變緩,減少的趨勢在降低,這與葉爾羌河流域逐步完善生態用水配套體系、強化水資源節約保護意識、新疆實施天然林保護措施以及強有力的水資源管理制度有顯著的關系。但總體上葉爾羌河下游地區水資源匱乏,胡楊林呈減少的趨勢,與王芳[37]、徐培軍等[38]的研究一致。

圖4 1995—2015年葉爾羌河流域土地利用變化Fig.4 The land use types change in the Yarkant River Basin from 1995 to 2015

3.2 胡楊林轉移時-空變化

圖5 1995—2015葉爾羌河流域胡楊林變化Fig.5 The change of Populus euphratica forests in the Yarkant River Basin from 1995 to 2015

為了揭示胡楊林空間格局的變化情況,通過統計、計算和制圖得到研究區胡楊林與其他土地利用類型間的占比和空間轉移分布圖(圖5)。1995—2005胡楊林轉移較為劇烈,其總轉移面積占胡楊林總面積的70.60%。胡楊林轉為其他土地利用類型面積較多,占胡楊林總面積的38.84%,主要向未利用地和草地轉移；而由其他土地利用類型轉為胡楊林的面積占胡楊林總面積的31.76%,主要由草地和其他林地轉入。2005—2015年胡楊林轉為其他土地利用類型面積較多,占胡楊林總面積的27.10%,主要向未利用地和草地轉移；而由其他土地利用類型轉為胡楊林的面積占胡楊林總面積的24.23%,主要由草地和其他林地轉入。1995—2005年,永安壩至夏河林場段受人類活動的影響劇烈,胡楊林和耕地過渡區域有大片的胡楊林被迫轉為耕地,耕地向葉爾羌河道擴張；而夏河林場至三河口段,在遠離葉爾羌河道的區域胡楊林面積逐漸減少,轉為草地和未利用地。胡楊林增加主要分布在河道兩側和水庫周邊,胡楊林更新較快,存在天然幼齡林。2005—2015年,隨著葉爾羌河下游水資源的逐漸減少,夏河林場至三河口段遠離河道的胡楊林逐漸衰敗,向未利用地轉移。在人類生態環境的重視以及產業結構轉型的條件下,胡楊林亂砍亂伐以及毀林開荒的現象減小,胡楊林增加區域主要分布在河漫灘與沖積平原地區,該區域土壤濕度較大,土壤含水量較高,土壤生物較多,能為胡楊林提供充足的水分和養分滿足胡楊生長需求。

圖6 1995—2015年葉爾羌河流域胡楊林時空間轉移Fig.6 Spatio-temporal transfer of Populus euphratica forests in the Yarkant River Basin from 1995 to 2015UTP：未利用地向胡楊林轉移,Unused land transfer to Populus euphratica forests；RTP：居民地向胡楊林轉移,Residential area transfer to Populus euphratica forests；WTP：水域向胡楊林轉移,Water area land transfer to Populus euphratica forest forests；GTP：草地向胡楊林轉移,Grassland land transfer to Populus euphratica forests；OTP：其他林地向胡楊林轉移,Other forestland transfer to Populus euphratica forests；FTP：耕地向胡楊林轉移,Cultivated land transfer to Populus euphratica forests；PTU：胡楊林向未利用地轉移,Populus euphratica forests transfer to unused land；PTR：胡楊林向居民地轉移,Populus euphratica forests transfer to residential area；PTW：胡楊林向水域轉移,Populus euphratica forests transfer to water area；PTG：胡楊林向草地轉移,Populus euphratica forests transfer to grassland land；PTO：胡楊林向其他林地轉移,Populus euphratica forests transfer to other forestland；PTF：胡楊林向耕地轉移,Populus euphratica forests transfer to cultivated land

為了定量揭示胡楊林與其他土地利用類型的相互轉移情況,通過土地利用類型轉移矩陣反映各土地利用類型之間相互轉移動態過程信息(圖6)。1995—2015胡楊林主要向草地、其他林地、未利用地和耕地轉移,同時部分未利用地、草地和其他林地轉為胡楊林。胡楊林向居民地和水域轉移面積較少,同時水域、耕地及居民地向胡楊林轉移較少。1995—2005年胡楊林向草地、其他林地和耕地面積分別轉移667.3 km2、252.66 km2、116.13 km2；草地、其他林地和未利用地分別向胡楊林轉移547.13 km2、199.85 km2、135.09 km2。2005—2015年胡楊林主要向未利用地、草地和其他林地面積分別轉移310.35 km2、221.07 km2、55.76 km2,而草地、其他林地和未利用地分別向胡楊林轉移296.46 km2、173.57 km2、67.25 km2。總體上,近20年內葉爾羌河流域天然胡楊林與草地、未利用地和其他林地轉移頻繁,而與居民地和水域相互轉移較少。

3.3 景觀格局動態變化

3.3.1基于景觀類型水平動態變化特征

草地和胡楊林斑塊密度(PD)較大,耕地、水域以及居民地斑塊密度較小(表3),表明草地和胡楊林分布較集中,而水域和居民地分布比較分散；胡楊林PD從1995年的0.19增加到2005年的0.56,然后再減少到2015年的0.35,總體呈增加趨勢。草地和未利用地最大斑塊指數(LPI)最大,而其他地類值較小,說明草地和未利用地占主要類型；胡楊林LPI從1995年的0.35增加到2015年0.58,說明胡楊林在研究時序內空間分布上逐漸連片。居民地散布與并列指數(IJI)最低,說明居民地分布較分散,鄰近度低,而其他地類相對集中。1995—2015年胡楊林IJI指數一直增加,從50.11增加到60.23,說明胡楊林分布越來越集中,彼此為鄰,分離度逐漸降低。其他林地和居民地斑塊結合度指數(COHESION)較低,而其他的五類土地利用類型斑塊結合度指數較高接近100,說明居民地和其他林地的連通性較差,分布相對分散,空間連接性較低。胡楊林斑塊結合度指數變化幅度較小,但呈增加趨勢,說明其分布逐漸集中,空間連接性越來越高。胡楊林景觀格局變化劇烈,胡楊林斑塊密度、最大斑塊指數和散布與并列指數呈增加趨勢,斑塊結合度指數相對較高。1995—2015年,永安壩至夏河林場段胡楊林逐漸向葉爾羌河道縮減；而夏河林場至三河口段因水資源減少和地下水位下降,在遠離葉爾羌河區域的胡楊林逐漸趨于死亡,導致胡楊林分布逐漸集中、景觀破碎度降低、連通性呈上升趨勢。

表3 葉爾羌河流域景觀類型水平指數變化

3.3.2基于景觀水平動態變化特征

圖7 葉爾羌河流域胡楊林景觀指數變化 Fig.7 Changes of landscape indices of Populus euphratica forests in the Yarkant River BasinSHDI：Shannon多樣性指數,Shannon′s diversity index；SHEI：Shannon均勻度指數,Shannon′s evenness index；CONTAG：蔓延度指數,Contagion index；SPLIT：分離度指數,Splitting index；LSI：景觀形狀指數,Landscape shape index

整體景觀形狀指數(LSI)呈先增后減的趨勢(圖7),分別由1995年的84.67增加到2005年的93.44,再下降至2015年的90.56,整體上增加5.89。說明研究區在20年內景觀形狀越來越不規則,越來越復雜多樣。Shannon多樣性指數(SHDI)和Shannon均勻度指數(SHEI)呈現相同的變化趨勢,1995—2005年,呈現增加的趨勢,分別由1995年的1.35和0.69增加至2005年的1.41和0.72,表明景觀內土地利用越來越豐富,破碎化程度較高,景觀異質性、多樣性及均勻度增加,景觀類型向復雜化發展。2005—2015年,呈現減少趨勢,分別下降至2015年的1.38和0.71,表明景觀內土地利用越來越單一,優勢度明顯,景觀異質性、多樣性及均勻度減弱,景觀類型向單一發展。蔓延度指數(CONTAG)呈現先減后增的趨勢,分離度指數(SPLIT)呈現先增后減的趨勢,表明1995—2005年,景觀要素連接性逐漸變弱,聚集程度降低,離散和分離程度增加,景觀趨于復雜化；2005—2015年,景觀要素連接性逐漸變強,聚集程度增加,離散和分離程度減弱,優勢度斑塊類型比例逐漸增加,景觀趨于穩定。總體上,1995—2015年間葉爾羌河流域LSI、SPLIT、SHEI、SHDI指數增加,CONTAG指數減少。景觀異質性、多樣性及均勻度增加、破碎和分離程度較高、景觀類型向復雜化發展。

3.4 歸因分析

3.4.1氣候變化

氣候變化是影響土地利用和景觀格局變化的主要因素,降水和氣溫變化又是氣候要素中相對活躍的影響因子；葉爾羌河流域是典型的冰雪融水補給河流[39],其出山口卡群水文站水文資料對研究葉爾羌河流域水資源具有重要的價值。因此,本文以降水、氣溫和卡群站徑流量三個氣象因子分析氣候變化對葉爾羌河流域胡楊林時空格局的影響。由圖8可見,1995—2005年,平均年降水量為67.29 mm,其中1996年與2002年降水量較大,分別為153.70 mm和152.70 mm,年際降水波動較大,平均降水量低；年平均氣溫為12.54℃,其中,1996年平均氣溫最低11.30℃,而2004年平均氣溫最高13.30℃,氣溫趨于平穩,徑流量變化較小。2005—2015年,平均降水量為60.84 mm,降水量趨于平穩,年際降水波動較小；年平均氣溫為12.61℃,其中,2012年平均氣溫最低為11.60℃,氣溫波動幅度較大,平均氣溫有所升高,徑流量變化較大。整體上,降水年際變化明顯,降水量呈減少趨勢；年際平均氣溫年際變化明顯,平均氣溫呈上升趨勢；徑流量呈增加趨勢；總體呈現出“冷濕-暖干”變化的趨勢。溫度升高加快冰雪融化,徑流量有增加的趨勢,而降水直接對冰川和徑流進行補給,降水增加使冰川區積累量增加,調節河道徑流量[40]。

葉爾羌河流域屬于典型的干旱地區,降水量、溫度和徑流量是影響胡楊林生長的主要因素。四十八團至永安壩段水資源豐富,胡楊林受氣候變化的影響較少,胡楊林有增加的趨勢。而永安壩至夏河林場段,胡楊林受氣候變化較為明顯,溫度升高使蒸發加快,而降水減少使該區域更加干旱,使胡楊林逐漸衰敗。夏河林場至三河口段自然條件更為惡劣,在氣候向“暖干”變化的趨勢下,胡楊林向河道逐漸縮減,逐漸轉為未利用地。總體上：溫度、徑流量、降水三者相互影響,氣溫升高,降水減少使胡楊林逐漸衰敗；徑流量有增加的趨勢,但在上游新修水利工程和人類活動的影響下,使中、下游地區徑流量減少,地下水位逐漸變深[41],導致下游胡楊林逐漸轉化為未利用地,夏河林場至三河口段尤為明顯。可見氣候變化對土地利用變化、胡楊林生長情況、草地等景觀格局影響較大。

圖8 1995—2015年葉爾羌河流域氣象因子統計Fig.8 Statistical chart of meteorological factors in the Yarkant River Basin from 1995 to 2015

3.4.2人類活動

人類活動是人類為了生存發展和提升生活水平而進行的活動,而這些活動是通過人口、地區生產總值等經濟社會指標表征。人口是人類社會經濟中最主要的因素；國內生產總值(GDP)是地區農業、工業和服務業等發展狀況的綜合反映,也反映了一個地區的經濟發展水平；第一產業對胡楊林、林地、耕地等景觀結構變化影響極大。因此,本文以人口、第一生產比重、GDP三個社會因子分析人類活動對葉爾羌河流域胡楊林時空格局的影響。從圖9可知：1995—2005年,人口增長速率較快,增加21.29%；GDP增長平穩,增長幅度小；第一產業比重先升高后降低。2005—2015年,人口增長較前10年較緩,增加18.02%；GDP增長迅速,增長幅度巨大,是2005年的4倍；第一產業比重持續下降。總體上,人口持續增加,從1995年的175.45萬人增加到20015年的251.12萬人,增加43.13%。人口增長使城鎮快速擴張,增加居民地面積,而為保障農民生活需發展農業經濟,因此破壞草地、其他林地以及胡楊林地,共使140.59 km2胡楊林轉移到耕地。GDP持續增長,呈現先緩后迅速的增長趨勢；第一產業比重呈現先增加后持續降低的趨勢,表明隨著經濟的發展,產業越來越多樣化,經濟發展迅速,第一產業比重先增加后持續下降,而使耕地面積增加,隨著第一產業比重的下降,草地和未利用地面積增加。

圖9 1995—2015年葉爾羌河流域人口、GDP及第一產業比重變化趨勢Fig.9 Change trends of population, GDP, and proportion of primary industry in Yarkant River Basin from 1995 to 2015

4 討論

葉爾羌河流域是塔里木河流域重要的生態系統,同時處于塔克拉瑪干沙漠北部,是葉爾羌綠洲的屏障,其地理位置極其特殊與重要。本文基于長時間序列遙感數據,通過土地利用指數、轉移矩陣以及景觀指數對葉爾羌河流域胡楊林時空格局變化進行動態監測與分析,為葉爾羌河流域生態系統恢復及胡楊林保護提供參考。本次研究時間序列適中,但研究時間間隔較長,由于胡楊林在沙漠化地區植被覆蓋度較低,同時與其他地類差異較小,通過遙感影像準確提取胡楊林難度較大。因此,后續研究可以根據遙感影像其顏色特征、紋理特征、光譜特征,結合計算機深度學習建立胡楊林特征數據庫,實現全自動提取胡楊林面積分布,更高效,更準確,更快速地完成分類,增加本研究時間序列和減少研究間隔。

本研究表明,耕地、居民地增加,而水域、草地、胡楊林減少,與前人研究一致。例如：張廣朋等[42]探討綜合治理項目實施前后葉爾羌河流域生態與環境結果表明：流域土地利用變化強烈,且2000—2010年比1990—2000年土地利用變化速度比慢,草地面積呈減少趨勢；許艾文等[43]通過遙感數據研究了1978—2015年喀喇昆侖山克勒青河流域(葉爾羌河支流)冰川變化情況,結果表明：近二十年來水域面積呈減少的趨勢,流域冰川近40年面積減少145.78 km2；王芳[37]研究表明耕地和居民地面積一直增加,區域土地利用變化是自然和人類活動共同作用的結果,而人口快速增長占主要因素。對葉爾羌河流域胡楊林研究表明[44- 46]：葉爾羌河流域氣候變化和人類活動共同影響胡楊林生境,而不同區域影響胡楊林景觀格局的主導因子不同。四十八團至夏河林場段,水資源相對豐富,盲目開荒、過度樵采、不合理采伐等人類活動是胡楊林減少的主要因素,胡楊林主要轉為耕地；而在夏河林場至三和口段,遠離城市,胡楊林主要受氣候變化影響,同時受葉爾羌河上、中游河道攔河筑壩、攔截洪水的影響,導致下游徑流量減少,地下水位下降,使胡楊向河道縮減逐漸變為沙漠。干旱區域水資源是影響區域景觀格局的關鍵因子,應如何恢復胡楊林生態系統,保障流域生態安全？李靖等[25]研究表明通過節流、治污、開源、合理調控水資源、提高水資源利用率和加強流域生態建設是提高該流域生態承載力和恢復植被生態系統的具體措施[25]。綜上所述,在氣候變化和人類活動的干預下,葉爾羌河流域胡楊林景觀格局發生了顯著變化,胡楊林逐漸衰敗。因此,構建生態用水法律保障,嚴格執行流域限額用水制度,統籌協調生活、生產和生態用水,調整生產結構,控制人口增長速度以及實施胡楊林保護措施,有望實現葉爾羌河流域水資源可持續利用和生態系統恢復的目標。

5 結論

(1)胡楊林面積持續減少,從1995年的1916.15 km2減少到2015年的1652.25 km2,面積減少了263.90 km2。胡楊林土地利用動態度和空間動態度分別從2005年之前的1.04%和10.35%降低到2005年之后的0.38%和6.77%；而變化狀態指數卻從1995—2005年的-0.10升高到2005—2015年的-0.06,表明胡楊林面積持續減少,減少速度漸緩,空間轉移劇烈。胡楊林主要分布在葉爾羌河下游四十八團至三河口河道兩側,且隨著遠離河道主軸線,植物種類逐漸減少,結構趨于簡單,胡楊林逐漸衰退。

(2)1995—2015年胡楊林地轉移較劇烈,主要向草地、耕地以及未利用地轉移。永安壩至夏河林場段,胡楊林主要向耕地轉移；在遠離葉爾羌河道的腹地(夏河林場至三河口段),胡楊林逐漸趨于死亡,轉為未利用地,而胡楊林主要增加在河漫灘、沖積平原地區、河道兩岸以及水庫周邊。

(3)草地和未利用地占主要的景觀類型,居民地分布比較分散；胡楊林斑塊密度、最大斑塊指數和散布與并列指數分別從1995年的0.19、0.35、50.11增加到2015年的0.35、0.58、60.23,胡楊林分布逐漸集中、景觀破碎度降低、連通性呈上升趨勢。整體上,景觀異質性、多樣性及均勻度增加,破碎化程度較高,聚集程度和優勢度斑塊類型比例降低,景觀類型向復雜化發展。

(4)氣候總體表現出“冷濕-暖干”的變化趨勢。1995—2005年在氣候變化和人類活動的強烈干預下,胡楊林面積減少較快,而2005—2015年隨著水資源節約保護意識的增強、胡楊林保護措施的實施、水資源管理制度的執行、第一產業比重的下降以及生態輸水等舉措的實施,胡楊林衰退減緩,但對胡楊林仍有巨大的威脅。