巴丹吉林沙漠湖泊季節變化的遙感監測*

朱金峰，王乃昂＊＊，李卓侖，董春雨，陸 瑩，馬 寧

(1:蘭州大學資源環境學院，蘭州730000)

(2:蘭州大學干旱區水循環與水資源研究中心，蘭州730000)

巴丹吉林沙漠湖泊季節變化的遙感監測*

朱金峰1，2，王乃昂1，2＊＊，李卓侖1，2，董春雨1，2，陸 瑩1，2，馬 寧1，2

(1:蘭州大學資源環境學院，蘭州730000)

(2:蘭州大學干旱區水循環與水資源研究中心，蘭州730000)

利用2001年12月28日、2002年3月18日、7月24日、9月26日、12月15日五期ETM+遙感影像，提取了巴丹吉林沙漠湖泊年內季節變化信息，并分析了湖泊季節變化的統計特征與空間分布特征.結果顯示，巴丹吉林沙漠湖泊總面積和數量隨季節更替呈現出明顯的統計特征與空間分布特征.湖泊總面積和數量在當年春、夏、秋、冬季均依次減少，到第二年春季又逐漸恢復到前一年春季狀態;湖泊面積季節變化萎縮與增大區域呈集聚、連片分布態勢，部分時段湖泊面積變化表現出漸變分布特征.研究結果可為進一步深入研究巴丹吉林沙漠湖泊水循環、湖泊水量平衡、湖泊生態系統提供科學參考.

巴丹吉林沙漠;湖泊;遙感;季節變化;湖泊水循環

在巴丹吉林沙漠高大沙丘之間的低地分布有許多內陸小湖泊(海子)，總數約144個，其中常年有水的湖泊達74個，主要集中分布在沙漠的東南部［1］.1990s以來，關于巴丹吉林沙漠湖泊水源、湖泊水循環等科學問題一直是國內外學者研究的熱點.這些研究主要通過地球化學、地質外推等方法來解釋湖泊成因、湖泊水量在極度干旱的氣候環境下為何長久不衰，以及湖水的補給途徑和補給來源等問題，但得出的結論至今說法不一［2-12］.遙感技術能夠大范圍、及時快速地獲取地表環境信息，為湖泊水域動態變化監測和相關參數的反演研究帶來了便利，成為湖泊研究強有力的技術手段［13］.目前國內外學者在這些方面已經開展了大量的研究工作，并取得了很大進展［14-22］.但這些工作都圍繞著研究區湖泊近幾十年的動態變化及其原因展開的，而沒有關注湖泊季節變化特征.湖泊水域季節變化與降水、蒸發、地下水活動的季節性波動關系密切，特別是分布在西北干旱半干旱地區的內陸湖泊，通過遙感監測分析其季節變化特征，可為深入研究湖泊水循環、湖泊水量平衡、湖泊生態系統提供科學參考.本文基于以上思路，利用ETM+遙感影像，提取巴丹吉林沙漠湖泊年內季節變化信息，并分析湖泊季節變化的統計特征與空間分布特征.

1 研究區概況

巴丹吉林沙漠位于雅布賴鹽湖與雅布賴山之西北，宗乃山、省道S218公路之西，黑河正義峽出山口、弱水東岸至古日乃湖之東，橫跨5個經度，東西長約442km;合黎山、北大山、黑山頭之北，拐子湖、古居延澤之南，縱貫3個緯度，南北寬約354km，面積為5.2162×104km2，系我國僅次于塔克拉瑪干沙漠的第二大沙漠［23］.區內地勢總體呈南高北低、東高西低的趨勢，海拔1000-1700m，全年主要受W-NW大陸風系控制，多年平均降水量由東南向西北逐漸減少，東南部約為120mm，西北部不足40mm，多年平均蒸發量大于3000mm，年平均氣溫1-8℃，7月份最高達37-41℃，沙面則高達70-80℃，氣候極度干旱［5］.與世界上其他沙漠明顯不同的是，東南部的高大復合型沙山，相對高度一般達200-300m，最高的超過430m;沙山之間的洼地分布著許多大小不等的永久性湖泊，絕大多數屬礦化度很高的鹽堿湖，還有一定數量的淡水湖，這一奇特景觀類型為世界罕見.本文選取沙漠東南部常年有水的湖泊區域(39°30'-40°07'N，101°42'-102°40'E)(圖1)作為研究區.

圖1 研究區位置及湖泊分布Fig.1 The location of study area and distribution of desert lakes

2 數據與方法

2.1 數據源及其預處理

研究所用的數據為全球陸地覆蓋設施中心——地球科學數據交換站點(GLCF-ESDI)提供的五期Land-sat-7 ETM+L1T級遙感影像，軌道號(Path/Row)為132/32、132/33，時間分別為2001年12月28日、2002年3月18日、7月24日、9月26日、12月15日，影像質量良好，平均云量均小于1%，且已經過幾何校正［24］，不同月份數據間達到高精度配準.根據研究區范圍對影像進行拼接、裁剪.采用PCA方法對ETM+多光譜波段與全色波段數據進行融合，此方法的融合效果在影像保光譜特性、信息量和清晰度三方面效果最佳［25］，融合后數據的空間分辨率為15m.采用Albers等積投影對影像進行投影轉換.影像的拼接、裁剪、融合、投影轉換均借助ERDAS IMAGINE 9.2圖像處理軟件完成.

2.2 湖泊光譜特征分析

相對于其他地物而言，湖泊水體對入射光具有強吸收性，在大部分遙感傳感器的波長范圍內，呈現較弱的反射率.清澈水體的遙感信息模型根據其反射率可以近似表示為:藍光＞綠光＞紅光＞近紅外＞短波紅外［26-29］.從研究區湖泊、鹽堿地、陰影、植被、沙丘五種地物光譜曲線可以看出，湖泊在藍光波段(ETM+1波段)的DN值最大，在近紅外(ETM+4波段)、短波紅外(ETM+5、7波段)同時具有強烈的吸收性，而其他四種地物的DN值均高于湖泊，特別在ETM+4、5、7波段(圖2).利用湖泊這種典型特征，可建立提取湖泊信息的水體指數.

圖2 湖泊、鹽堿地、陰影、植被、沙丘的光譜曲線Fig.2 Spectral curves of lakes，saline and alkaline land，shadow，vegetation and sand dune from Landsat-ETM+images

2.3 湖泊信息提取

2.3.1 湖泊信息提取約定 為了真實反映湖泊數量、面積大小隨季節變化特征，在提取湖泊信息之前，約定如下:①湖泊指水域覆蓋地區，不包括湖泊內的島嶼、湖泊周圍沼澤化草甸、鹽生草甸等植被、現代干涸鹽堿湖盆等地區;②野外考察和影像初步解譯均得出在冬季(2001年12月28日、2002年12月15日)影像上，部分湖面邊緣水域出現結冰現象，在湖泊信息提取時，這些結冰區域應當包括在湖泊之內;③影像上大于等于4個像元(900m2)的湖泊均提取出來;④由于湖泊季節性萎縮，部分湖泊出現分裂現象，即一個大湖泊縮成幾個小湖泊，在湖泊信息后處理時，將這幾個小湖泊合并，其面積和數量統一記為一個湖泊.

2.3.2 湖泊信息提取方法 前人已經建立了多個水體指數提取水體信息，如歸一化水體指數NDWI(Normalized Difference Water Index)［26］、改進歸一化水體指數 MNDWI(Modified NDWI)［27］、增強型水體指數 EWI(Enhanced Water Index)［28］、新型水體指數 NWI(New Water Index)［29］.本研究根據沙漠湖泊典型的光譜特征，提出了一種用于提取沙漠湖泊信息的水體指數DLWI，其公式如下:

其中:bblue、bmir分別代表ETM+影像第1、5波段的DN值;C為常數，本研究取100，目的是將DLWI的數值區間進行拉伸，方便閾值的確定.

建立的DLWI指數可以輕易實現湖泊水體與其他地類的區分，這主要是由于:①湖泊在第1波段具有最大的DN值，而在第5波段的DN值很低(圖2)，由DLWI指數得到的水體信息值均在正值區間;② 鹽堿地、陰影、植被和沙丘在第5波段的DN值遠遠大于其各自在第1波段的DN值，由DLWI指數得到的這4種地物的信息值均在負值區間.因此用閾值法就可達到正確提取湖泊信息的目的.

2.3.3 湖泊信息后處理 對DLWI指數法提取的湖泊信息，依據“湖泊信息提取約定”進行人工處理.對于春、夏、秋三季對應的影像，目視解譯對湖泊邊界進行修正，并去除小于4個像元的湖泊信息.對于冬季對應的影像，存在兩個問題并對其處理如下:①由于成像時太陽高度角為一年中最小時段，高大沙山陰影與湖泊水體出現“異物同譜”現象，對提取后的信息目視去除所有沙山陰影;②由于冰面與湖泊水體光譜特征存在差異，去除沙山陰影后的湖泊信息沒有包含冰面，故目視補充解譯冰面信息.

根據地形圖信息為提取的湖泊空間數據賦屬性值，最終形成研究區湖泊的春、夏、秋、冬四季空間數據庫.總體研究技術路線如圖3所示.

圖3 研究技術路線Fig.3 Completing routines for research

3 結果與分析

3.1 湖泊季節變化統計特征

利用ETM+影像提取的研究區湖泊各季節信息及變化情況可以看出，巴丹吉林沙漠湖泊總數量和總面積隨季節更替均有所變化，春季湖泊總數量和總面積達到年內最大值，夏、秋、冬季湖泊總數量和總面積均逐漸減小，到了冬季湖泊總數量和總面積減小到年內最小值.從2001年冬季到2002年春季湖泊總數量和總面積又開始增大，重新達到年內最大值(表1).

表1 湖泊信息提取結果統計Tab.1 The statistical results of desert lakes interpreted from Landsat-ETM+images

依據研究區湖泊面積大小變化范圍，將湖泊分為大于 1、0.5-1、0.5km2以下三種級別，分別稱為“大”、“中”、“小”湖泊，其分級信息如表2所示;并將湖泊季節變化分級的轉化信息列出如表3所示，進一步分析湖泊季節變化的分級效應發現，“大”湖泊隨季節更替的變化最明顯，湖泊總面積由春至冬依次遞減，春、夏、秋、冬季分別為6.84、4.73、3.64、3.61km2，其中，春夏之交2個“大”湖泊萎縮成“中”湖泊，夏秋之交1個“大”湖泊萎縮成“中”湖泊，而由冬至春，有3個“中”湖泊增大成“大”湖泊，補給了湖泊面積;“中”湖泊隨季節更替的變化也比較明顯，湖泊面積春、夏、秋、冬季分別為1.69、3.16、4.03、3.85km2，呈先增大后減小的變化趨勢，這主要是由于在春夏之交和夏秋之交個別“大”湖泊萎縮成“中”湖泊形成對“中”湖泊面積的補給，由冬至春3個“中”湖泊又增大成“大”湖泊，致使“中”湖泊面積有所減小;“小”湖泊隨季節更替的面積變化相對不太明顯，但湖泊數量變化由春至冬呈明顯的減少特征，春、夏、秋、冬季分別為100、87、79、76個，排除春夏之交1個“中”湖泊萎縮成“小”湖泊對“小”湖泊數量的補給后可看出，春夏之交、夏秋之交、秋冬之交干涸的湖泊數量分別為14、7、3個，而由冬至春25個“小”湖泊再現水域，補給了湖泊面積.

由此可知，湖泊總面積和數量隨季節更替均有較大變化，當年春、夏、秋、冬季湖泊總面積和數量均依次減少，到第二年春季，湖泊總面積和數量又逐漸恢復到前一年狀態;由于湖泊面積的季節變化性，“大”、“中”、“小”湖泊之間存在相互轉化特點，部分“小”湖泊表現出“干涸”與“再現水域”的變化現象，這些變化均表現在湖泊總面積與數量季節總體變化特征的過程中.

表2 湖泊分級信息Tab.2 The classified information of desert lakes

表3 湖泊季節變化分級轉化信息Tab.3 The transformation of classified information of desert lakes

3.2 湖泊季節變化空間分布特征

為了分析湖泊季節變化空間分布特征，利用ArcGIS 9.2軟件對各季節湖泊空間數據進行疊加分析，得到各季節湖泊面積變化量，處理如下:①根據空間矢量數據分級符號化方法［30］，按照Natural Breaks方案對湖泊面積季節變化量進行分級符號化處理，得到其點層上的空間分布格局;② 基于GIS的空間插值方法，在采用交叉驗證法對反距離加權法(Inverse Distance Weighting，IDW)、普通克里格法(Ordinary Kriging，OK)、泛克里格法(Universal Kriging，UK)三種插值效果驗證的基礎上(表4)［31］，選擇IDW和OK法作為精度和效果最優的方法，分別對夏-秋、秋-冬和冬-春、春-夏湖泊面積季節變化量進行插值模擬，得到其面層上的空間分布格局(圖4).現結合圖1和圖4，對湖泊季節變化空間分布特征解釋如下:

(1)由春至夏，研究區西北部格力圖、達布蘇圖和中北部包爾準圖、準吉格德周圍形成兩大湖泊萎縮區域，梧桐圖、諾爾圖到音德爾圖一帶也是湖泊萎縮較大的區域，萎縮的湖泊共計66個，其中位于兩大湖泊萎縮區域中干涸的季節性湖泊達到7個;與此同時，東南部寶日陶勒蓋、準敖格欽一帶的淡水湖形成湖泊增大區域，其他地區亦有個別增大湖泊呈零星分布，增大的湖泊共計43個.在總體空間趨勢上，從研究區西北部到東南部，湖泊面積變化表現出萎縮-基本不變-增大的漸變分布特征.

表4 三種空間插值方法的交叉驗證結果Tab.4 Cross-validation results of the three spatial interpolation methods

圖4 湖泊面積季節變化量空間分布格局Fig.4 Spatial pattern of the lake area seasonal variation

(2)由夏至秋，研究區西北部木日圖和哈拉木格圖形成兩大湖泊萎縮區域，音德爾圖、浩尼吉林一帶也是湖泊萎縮較大的區域，萎縮的湖泊共計92個，其中位于兩大湖泊萎縮區域中干涸的季節性湖泊有2個;期間，增大的湖泊只有3個，分別是位于中部的布爾特和東南部的哈拉特日圖、巴丹湖.與春夏之交相比，此期間萎縮的湖泊數量大大增加.在總體上湖泊面積變化趨勢雖不明顯，但仍可以看出從西北到東南湖泊萎縮量逐漸減少的變化特征.

(3)由秋至冬，研究區中東部音德爾圖、呼隆特和西南部陶來圖、昭爾格圖一帶形成兩大湖泊萎縮區域，北部腦滾諾爾、陶蘇吉林周圍也是湖泊萎縮較大的區域，萎縮的湖泊共計68個，其中干涸的季節性湖泊共計15個;西北部達布蘇圖、烏蘭吉林和東南部希勒克、通古圖一帶形成湖泊增大區域，增大的湖泊共計31個，其中再現水域的季節性湖泊共計12個.與春夏之交、夏秋之交相比，此期間湖泊面積變化總體空間趨勢很明顯，增大的湖泊在西北-東南方向上形成一帶狀區域，而在此帶狀區域的中部兩邊則為湖泊萎縮區域.

(4)由冬至春，研究區中東部音德爾圖、諾爾圖和西南部布爾特、伊和吉格德周圍形成兩大湖泊增大區域，北部的腦滾諾爾、南部的杜貴吉林一帶也是湖泊增大較大的區域，增大的湖泊共計96個，其中再現水域的季節性湖泊達到25個;期間亦有個別面積減小的湖泊，如北部的蘇勒圖、南部的哈拉特日圖等.此期間是湖泊總面積與數量補給的時間段，湖泊面積變化表現出總體增大的空間分布格局，只在南部淡水湖區域出現較小萎縮區域.

從以上分析看出，湖泊面積季節變化在空間分布上特征明顯，湖泊萎縮與增大的區域呈集聚、連片分布態勢，部分時段湖泊面積變化表現出漸變分布特征.

4 結論

利用ETM+遙感影像提取了巴丹吉林沙漠湖泊年內季節變化信息，從統計結果和空間分布兩方面分析了湖泊季節變化特征，得出以下結論:(1)巴丹吉林沙漠湖泊總面積和數量在當年春、夏、秋、冬季均依次減少，到第二年春季湖泊總面積和數量又逐漸恢復到前一年狀態，在此變化過程中伴有“大”、“中”、“小”湖泊之間的相互轉化以及部分“小”湖泊“干涸”和“再現水域”的現象特征.(2)在空間分布上，巴丹吉林沙漠湖泊面積季節變化特征明顯，湖泊萎縮與增大區域呈集聚、連片分布態勢，部分時段湖泊面積變化表現出漸變分布特征.

由于缺少2002年6月ETM+數據，本文在遙感影像獲取時間與季節變化對應問題上尚存缺陷.本研究組已于2009年末建成巴丹吉林沙漠野外觀測實驗站.利用全年12月份遙感數據逐月監測沙漠湖泊變化，將其結果與沙漠實驗站觀測資料對比分析，進一步解釋巴丹吉林沙漠湖泊年內季節變化的原因、湖泊水循環、湖泊水量平衡等科學問題，是下一步研究的重點.

致謝:本文在寫作和修改過程中得到蘭州大學地球系統科學研究所張建明副教授、黃銀洲博士、李育博士等的幫助，兩位審稿專家對本文修改也提出了重要建議，在此一并感謝.

［1］朱震達，吳 正，劉 恕等.中國沙漠概論.北京:科學出版社，1980:73-76.

［2］Yang XP，Martin AJW.The ion chemistry of lakes and late Holocene desiccation in the Badain Jaran Desert，Inner Mongolia，China.Catena，2003，51(1):45-60.

［3］Yang XP.Chemistry and late Quaternary evolution of ground and surface waters in the area of Yabulai Mountains，western Inner Mongolia，China.Catena，2006，66(1-2):135-144.

［4］楊小平.巴丹吉林沙漠腹地湖泊的水化學特征及其全新世以來的演變.第四紀研究，2002，22(2):97-104.

［5］馬妮娜，楊小平.巴丹吉林沙漠及其東南邊緣地區水化學和環境同位素特征及其水文學意義.第四紀研究，2008，28(4):702-711.

［6］馬金珠，陳發虎，趙 華.1000年以來巴丹吉林沙漠地下水補給與氣候變化的包氣帶地球化學記錄.科學通報，2004，49(1):22-26.

［7］馬金珠，黃天明，丁貞玉等.同位素指示的巴丹吉林沙漠南緣地下水補給來源.地球科學進展，2007，22(9):922-930.

［8］Gates JB，Edmunds WM，Darling WG et al.Conceptual model of recharge to southeastern Badain Jaran Desert groundwater and lakes from environmental tracers.Applied Geochemistry，2008，23(12):3519-3534.

［9］陳建生，趙 霞，盛雪芬等.巴丹吉林沙漠湖泊群與沙山形成機理研究.科學通報，2006，51(23):2789-2796.

［10］陳建生，凡哲超，汪集旸等.巴丹吉林沙漠湖泊及其下游地下水同位素分析.地球學報，2003，24(6):497-504.

［11］趙 霞，陳建生.相似優先比法研究巴丹吉林沙漠及周邊地區地下水補給.湖泊科學，2006，18(4):407-413.

［12］Chen JS，Li L，Wang JY et al.Water resources:Groundwater maintains dune landscape.Nature，2004，432:459-460.

［13］王海波，馬明國.基于遙感的湖泊水域動態變化監測研究進展.遙感技術與應用，2009，24(5):674-684.

［14］Zhao SQ，Fang JY，Miao SL et al.The 7-decade degradation of a large freshwater lake in central Yangtze river，China.Environmental Science and Technology，2005，39:431-436.

［15］殷立瓊，江 南，楊英寶.基于遙感技術的太湖近15年面積動態變化.湖泊科學，2005，17(2):139-142.

［16］魯安新，王麗紅，姚檀棟.青藏高原湖泊現代變化遙感方法研究.遙感技術與應用，2006，21(3):173-177.

［17］馬國明，宋 怡，王雪梅.1973-2006年新疆若羌湖泊群遙感動態監測研究.冰川凍土，2008，30(2):189-195.

［18］牛沂芳，李才興，習曉環等.衛星遙感檢測高原湖泊水面變化及與氣候變化分析.干旱區地理，2008，31(2):284-290.

［19］王麗紅，趙 杰，魯安新等.基于遙感的西藏S301公路沿線湖泊變化研究.遙感技術與應用，2008，23(6):658-661.

［20］賈珅玥，肖鵬峰.基于多時相圖譜的青藏高原湖泊變化檢測研究.國土資源遙感，2009，82:78-85.

［21］李晶晶，賈建華，郝景研.基于RS的松嫩平原大安湖泊群面積提取與動態變化分析.遙感信息，2009，(3):44-48.

［22］熊 波，陳學華，宋孟強等.基于RS和GIS的沙漠湖泊動態變化研究——以巴丹吉林沙漠為例.干旱區資源與環境，2009，23(8):91-98.

［23］朱金峰，王乃昂，陳紅寶等.基于遙感的巴丹吉林沙漠范圍與面積分析.地理科學進展，2010，29(9):1087-1094.

［24］Williams D.Landsat7 Science Data Users Handbook［EB/OL］.http://landsathandbook.gsfc.nasa.gov/handbook.html.

［25］楊麗萍，夏敦勝，陳發虎.Landsat7ETM+全色與多光譜數據融合算法的比較.蘭州大學學報(自然科學版)，2007，43(4):7-11.

［26］McFeeters SK.The use of normalized difference water index(NDWI)in the delineation of open water features.Internation-al Journal of Remote Sensing，1996，17(7):1425-1432.

［27］徐涵秋.利用改進的歸一化差異水體指數(MNDWI)提取水體信息的研究.遙感學報，2005，9(5):589-595.

［28］閆 霈，張友靜，張 元.利用增強型水體指數(EWI)和GIS去噪音技術提取半干旱地區水體信息的研究.遙感信息，2007，6:62-67.

［29］丁 鳳.一種基于遙感數據快速提取水體信息的新方法.遙感技術與應用，2009，24(2):167-171.

［30］湯國安，楊 昕.ArcGIS地理信息系統空間分析實驗教程.北京:科學出版社，2006:125-134.

［31］林忠輝，莫興國，李宏軒等.中國陸地區域氣象要素的空間插值.地理學報，2002，57(1):47-54.

RS-based monitoring seasonal changes of lake in Badain Jaran Desert

ZHU Jinfeng1，2，WANG Naiang1，2，LI Zhuolun1，2，DONG Chunyu1，2，LU Ying1，2＆ MA Ning1，2
(1:College of Earth and Environmental Science，Lanzhou University，Lanzhou 730000，P.R.China)
(2:Center for Hydrologic Cycle and Water Resources in Arid Region，Lanzhou Univesity，Lanzhou 730000，P.R.China)

Using ETM+remote sensing images acquired on Dec.28 in 2001，March 18，July 24，Sept.26，Dec.15 in 2002，respectively，this study extracted the information of seasonal changes of the number and areas of lakes and analyzed the statistical characteristics and spatial distribution of them.The following conclusions were drawn.Firstly，the number and total area of desert lakes showed a decreasing trend from spring，summer，autumn to winter in the first year and regained the original condition gradually in the next spring.Secondly，the seasonal variation of lake area displayed an obvious characteristic that the atrophic and increased area of lake present a distribution of gathers and continual pieces and gradual changes in the spatial distribution showed in some time section.The automatic weather observation station of Badain Jaran Desert had been completed in winter in 2009.Our next focus is to study the scientific issues of lake water cycle further，water balance both using RS imagery data and meteorological observation data.

Badain Jaran Desert;lake;remote sensing;seasonal changes;lake water cycle

＊ 國家基礎科學人才培養基金項目(J0730536)和高等學校博士學科點專項科研基金項目(20090211110025)聯合資助.2010-07-06收稿;2010-09-02收修改稿.朱金峰，男，1985年生，碩士研究生;E-mail:zhujf08@lzu.cn.

＊＊ 通訊作者;E-mail:wangna@lzu.edu.cn.