張永輝
(北京國遙新天地信息技術有限公司,北京 100000)
水域形態特征的研究對了解區域水資源動態變化、生態環境變化、氣候變化以及可持續發展具有重要意義[1]。遙感技術的飛速發展為水域動態變化的監測提供了有利條件。楊日紅等[2]在對色林錯湖的遙感數據進行分析,得出可能是由于溫室效應使得湖水面積增大的結論;殷立瓊等[3]利用Landsat遙感數據,研究了近15年來太湖面積變化狀況,并分析了引起太湖水域面積減小的原因;謝元禮等[4]利用多期TM影像數據,提取渭河寶雞段的水域范圍,并通過總結近20年間水域面積的變化情況,分析了該區域水環境的變化情況和水工程治理的實際效果;李輝等[5]基于MODIS影像對鄱陽湖湖面積及與水位關系進行研究,常利等[6]基于資源三號衛星對丹江口水庫面積季節變化進行遙感監測,張行清等[7]基于HJ-1衛星數據對廣西水庫面積進行了監測;迪麗努·爾阿吉等[8]選用了不同時期的LANDSAT衛星遙感影像,對新疆典型湖泊水域面積的動態變化展開了監測和研究。然而這些研究仍然存在:①對水體檢測和邊界提取缺乏不同方法結果的對比分析,整體解譯精度有待提高;②缺乏足夠有效的資料支持和合理的研究方法,很難定量說明水域形態特征變化的主導原因等問題。
為此,本文以新疆阿勒泰西部地區河流湖泊為研究對象,研究該地區水域動態變化規律以及變化趨勢,充分利用遙感、氣象、水文等資料分析引起水域面積發生變化的主要因素,從而為合理開發和持續利用新疆湖泊資源以及制定湖泊資源保護政策提供有價值的參考,并對新疆地區應對湖泊的鹽漬化、修建水庫、抽水補給、工農業發展做出合理的建議。
阿勒泰地區下轄六縣一市[9],即阿勒泰市、布爾津縣、哈巴河縣、吉木乃縣、福海縣、富蘊縣、青河縣,本文以位于阿勒泰地區西部的前6個縣市作為研究區。阿勒泰地處歐亞大陸中心腹地,遠離海洋,屬于中溫帶大陸性氣候區。由于受路經西部山口、塔城和布克賽爾、阿爾泰山東部山谷不同性質氣流的影響,在加上地形復雜,高度懸殊,地面性質差異等原因,從而形成了春旱多風、夏短少炎、秋高氣爽、冬季漫長和多大風的氣候特點[10,11]。
阿勒泰地區是新疆豐水區之一。有額爾齊斯河、烏倫古河、吉木乃山溪三大水系,大小河流共56條,年總徑流量133.7億m3,其中,地表徑流量為122.5億m3,占自治區地表徑流量的14.7%。開發利用的地表徑流量僅28億m3,其余都流入哈薩克斯坦共和國境內。全區有各類水庫共48座(其中大型水庫2座、中型水庫11座、小型水庫35座),總庫容9.36億m3。
本文獲取的數據包括氣象、遙感等數據,具體如下所述。
(1)氣象數據:本研究選用的氣象數據為阿勒泰氣象觀測站2000~2013年各月溫度、降水數據及年平均溫度、降水數據,以上數據均由中國氣象局國家氣候中心網站(https://www.ncc-cma.net/cn/)免費提供。
(2)影像資料:藏北7、8、9月份為雨季,收集該時相的遙感影像具有很強的對比性且該時期有利于植被的提取。本研究所用的影像資料包括研究區的2000~2013年7、8、9月份的多期影像數據。影像范圍中包含有新疆第二大淡水湖烏倫古湖以及額爾齊斯河,數據來源于中國科學院計算機網絡信息中心國際科學數據鏡像網站(http://www.cnic.cas.cn/)。本文對阿勒泰地區多期遙感影像數據進行幾何校正、裁剪等預處理。
基于該區域的Landsat TM/ETM+/OLIS影像,本文采用“綜合分析法”逐期對影像中的水體等進行解譯,對最優的處理結果進行面積統計并分析變化趨勢。利用MAPGIS6.7軟件對耕地進行矢量化,并統計其面積。結合該地區近14年來的降水量與氣溫的數據信息,從水文氣象的角度分析引起阿勒泰地區水域面積變化的主要因素。
前人提出的水體提取的方法有很多,使用較為廣泛的有歸一化水體指數法[12]NDWI=(Green-NIR)/(Green+NIR)和改進的歸一化水體指數法[13]MNDWI=(Green-SWIR)/(Green+SWIR)。本文通過分析阿勒泰西部地區的光譜亮度值(圖1),嘗試用基于閾值的比值法提取水體,公式為:
(1)
式(1)中:SWIR2為中紅外波段,為TM/ETM+、OLI_TIRS影像中的第7波段;Blue為藍色波段,TM/ETM+的第1波段、OLI_TIRS的第2波段。
通過實驗對比NDWI、MNDWI和比值法的提取結果,發現將閾值取為[0.24,0.75]的比值法,提取效果優于NDWI和MNDWI,水域解譯精度可達90%左右,能較好地滿足實驗要求,因此選用比值法提取水域面積。

圖1 水體、耕地、裸地和建筑物的波譜特征
根據水體提取統計的數據結果(表1),2000~2013年水域面積總體增長126.93 km2,其中,2000~2005年水域面積減少了28.268 km2。統計結果表中,水域面積變化規律呈波動變化趨勢,2000~2005年呈持續遞減趨勢,2005~2013年水域面積持續增加,且2013年達到最大值。總體來看,研究區域的水域面積呈增長趨勢。

表1 水域面積統計
因為植被和耕地難以區分,本文采用目視解譯及與高分辨率影像(Googleearth)相結合的方法將耕地從植被中分離,對其進行矢量化,并統計其面積。對解譯出的矢量文件中選取檢驗點:計算經緯度,利用Google Earth查詢,檢驗點落在耕地區域內,視為該點正確。經計算耕地解譯精度在85%以上,能夠滿足分析要求。和通過遙感影像的對比發現,人為的耕地面積在不斷的擴大,且主要集中在額爾齊斯河流域北岸,以及烏倫古湖周圍。根據統計數據(表2)可知,2000~2013年耕地面積呈持續增長趨勢。2000年耕地面積為2090.886 km2,2013年的耕地面積增加到了3195.376 km2,間隔13年,面積增加了1104.49 km2。
從圖2可以看出,阿勒泰地區7、8、9月溫度和降雨變化趨勢與該時間段(表1)的水域面積變化趨勢相一致,并且全區的年均溫度持續升高,近年來該區水域面積也不斷增加。這反映了溫度的變化對該區的河流-湖泊面積的變化有著直接的關系。

表2 耕地面積統計
經分析可知:隨著全球氣候變暖,我國西北地區溫度呈持續升高[14~20],氣溫的升高,使得水的蒸發量增加與蒸發速度加快,促進了區域性的水汽循環,新疆地區的氣候由暖干型向暖濕型轉變,氣溫的升高,致使阿勒泰地區降水量呈上升趨勢,研究區水域面積則隨著氣象條件的變化做出一系列反應。同時,氣溫的升高會導致該區域的阿爾泰山冰雪融化速度加快,雪線上升,冰雪融水量增加,而這些現象可對研究區域的地表徑流、地下徑流產生直接的影響,促使徑流量增大、徑流速度加快,使得河流、湖泊水位上升,總體水域面積增加。

圖2 7~9月溫度和降雨變化
由研究區域的統計信息得出,人為的耕地面積在不斷擴大,近年來呈不斷的增長趨勢。隨著人口、耕地、畜牧業的大幅度增長,人們為了奪取農牧業豐收而不斷完善引水條件,不斷修建水庫、攔河大壩、河渠,由單純的引水發展到引儲兼備,大大改變了水分的空間分配規律。
由于人口的增多,農業、畜牧業經濟化的快速發展,人們需要不斷地擴大耕地、牧場規模,興修水庫、水電站、引水渠,這些水利土木的開發對河流徑流、湖泊注水量產生很大的影響。在研究的區域中,人為干擾對烏倫古湖的影響較大。烏倫古湖水位的升降有氣候方面的原因,更重要的卻是人類活動的影響。1970年以前,烏倫古河是烏倫古湖水量唯一補給來源,烏倫古湖水位升降及水質變化,由烏倫古河流入的水量所決定。1970年修建引額爾齊斯河水濟布倫托海工程,但引水量很小;1987年以后,“引額濟海”擴建工程竣工,使烏倫古湖水位變化從純天然狀態變為天然加人為干預的狀態。烏倫古湖流域不再是封閉流域。當烏倫古河水量減少乃至斷流時,湖泊補水主要依靠額爾齊斯河。近10多年烏倫古湖的湖岸向東北推進情況比較明顯,在湖的東北角形成了數個新的小型湖水域[21~24]。
隨著人們對生態環境和濕地重要性認識的進一步提高,以及國家“退耕還林”、“三北防護林工程”、“天然林保護工程”等政策的實施,2001年成立了阿勒泰科克蘇濕地自然保護區,2008年新疆阿勒泰山濕地成為歐盟資助項目,2009年阿勒泰市向科克蘇濕地自然保護區進行生態補水等一系列政策,使得額爾齊斯河流域濕地得到一定的修養和恢復,這些成為該區天然濕地面積擴大的政策驅動。此外,該區雖然水資源相對較豐富,但空間分布不均勻,額爾齊斯河上游和烏倫古河的徑流量年際變化較大,對農牧業生產穩定發展不利,為調節該區用水問題,修建了人工水庫坑塘等。這些人為修建水庫成為該區人工濕地面積增加的人文動因,同時人工水庫坑塘也為水澆地的開墾提供了客觀條件,使得水澆地面積持續增加,這些都成為該區人工濕地面積有所增加的人為影響。
通過對遙感數據與氣候數據的綜合分析,阿勒泰地區湖泊面積自2005年以來呈增長的趨勢。研究分析認為,阿勒泰地區水域面積變化的原因主要體現在以下幾個方面。
(1)隨著近年來溫度的升高,導致高山積雪融化加速,積雪量不斷減少,從而加大了河流湖泊的補給量。
(2)降水量的增多,直接引起河流補給量的增多,從而引起水域面積的增大。
(3)經調研發現,近年來阿勒泰地區耕地面積不斷擴大,對耕地的灌溉大多采用大水漫灌的方式,人為的將水從河流中上游引到下游,可直接使尾閭湖的水位抬高,面積增大。
(4)近年來,人類大量修建水庫、水電站,硬化地面,人為的調配水資源、改變水資源的分布格局,導致河流終端湖泊水位上升,面積增大。