基于ICESat/GLAS數(shù)據(jù)的可可西里地區(qū)湖泊水位變化研究*

李 龍，姚曉軍，李風(fēng)賢，李曉鋒

(1.蘭州資源環(huán)境職業(yè)技術(shù)學(xué)院，甘肅蘭州 730021； 2.西北師范大學(xué)地理與環(huán)境科學(xué)學(xué)院，甘肅蘭州 730070)

0 引言

自19世紀(jì)后期以來，全球平均地表溫度(GMST)持續(xù)上升。過去30年比任何有氣象資料以來的年份都暖，其中2000年以來的10年是有史以來最暖的10年； 相比氣溫變化而言，全球陸地表面降水總體變化并不大， 30°N～60°N降水增加較為明顯[1]。水是人類賴以生存的自然資源、經(jīng)濟(jì)資源和生態(tài)資源[2]。湖泊作為陸地水圈的重要組成部分，其水域變化是其所在流域水量平衡的綜合結(jié)果[3]，是氣候變化敏感的指示器[4]，湖泊面積變化會(huì)加劇了湖區(qū)洪澇災(zāi)害[5]。青藏高原是地球上海拔最高、數(shù)量最多、面積最大的高原湖群區(qū)，共有面積大于1.0km2湖泊1 055個(gè)，分別占我國(guó)湖泊總數(shù)量和總面積的39.2%和51.4%[6]。研究表明，青藏高原近50年年均氣溫以0.37℃/10年的速率上升，降水則以9.1mm/10年的速率增加[7]。

湖泊水位數(shù)據(jù)是湖泊水量平衡的重要指標(biāo)[8， 9]，也是湖泊潰決模擬研究的重要參數(shù)[10]。早期的湖泊水位監(jiān)測(cè)主要依賴于水位實(shí)測(cè)數(shù)據(jù)，但由于絕大多數(shù)湖泊都位于環(huán)境惡劣、人煙稀少的無人區(qū)，無法得到湖泊水位數(shù)據(jù)。隨著衛(wèi)星測(cè)高技術(shù)的不斷發(fā)展，衛(wèi)星測(cè)高數(shù)據(jù)已廣泛應(yīng)用于在全球范圍內(nèi)大型湖泊的水位動(dòng)態(tài)監(jiān)測(cè)。如Zhang等[11]利用ICESat/GLAS數(shù)據(jù)測(cè)得青藏高原部分鹽湖的水位變化，發(fā)現(xiàn)青藏高原大部分鹽湖水位上升明顯，才多茶卡水位上升速率最快，為0.80m/年； Wang等[12]利用ICESat/GLAS數(shù)據(jù)研究發(fā)現(xiàn)中國(guó)青藏高原湖泊水位上升，雅魯藏布江流域、內(nèi)蒙古北部地區(qū)和新疆地區(qū)湖泊水位下降，東部平原地區(qū)和云貴高原地區(qū)湖泊水位出現(xiàn)強(qiáng)烈波動(dòng)現(xiàn)象，但影響因素不盡相同； 同樣，Zhang等[13]基于ICESat/GLAS數(shù)據(jù)對(duì)我國(guó)10個(gè)最大湖泊進(jìn)行研究，發(fā)現(xiàn)青藏高原地區(qū)的青海湖、納木錯(cuò)、色林錯(cuò)和東北地區(qū)的興凱湖水位上升明顯，新疆地區(qū)的博斯騰湖和呼倫湖水位下降顯著，揚(yáng)子江流域的洞庭湖、鄱陽湖、太湖和洪澤湖水位無明顯變化； Krause等[14]利用分布式水文模型得到納木錯(cuò)水位近年來升高0.35m/年，Kropáek等[15]利用ICESat/GLAS數(shù)據(jù)測(cè)得納木錯(cuò)水位升高0.31m/年，對(duì)分布式水文模型精度做了很好驗(yàn)證。

可可西里地區(qū)是青藏高原湖泊集中分布區(qū)之一，目前有大于10km2的湖泊83個(gè)，并且湖泊之間存在合并現(xiàn)象[16]，湖水外溢事件格外受學(xué)者關(guān)注[17]。文章通過對(duì)可可西里地區(qū)主要湖泊文獻(xiàn)資料的整理，結(jié)合統(tǒng)計(jì)資料和地理信息技術(shù)獲取1993—2009年該地區(qū)主要湖泊水位變化數(shù)據(jù)，系統(tǒng)闡述湖泊水位變化特征并對(duì)其影響因素進(jìn)行探討。

1 數(shù)據(jù)來源與研究方法

1.1 研究區(qū)概況

可可西里地區(qū)深處青藏高原腹地，介于33°30′N～36°29′N， 81°56′E～94°06′E之間，面積約23.5萬km2(圖1)。在行政區(qū)劃上，可可西里地區(qū)橫跨我國(guó)的青海省和西藏自治區(qū)，包括格爾木、治多、班戈、尼瑪和改則等市(縣)部分區(qū)域。地理上的可可西里指的是以可可西里山為主體的鄰近山原湖盆地區(qū)，位于喀喇昆侖山以東，青藏公路以西，唐古拉山以北，昆侖山以南[14]，較人們一般認(rèn)為的可可西里地區(qū)范圍要大[18]。可可西里地區(qū)地勢(shì)高亢，平均海拔4 600m以上，最高峰海拔6 973m，為木孜塔格峰[19]。區(qū)內(nèi)中部較低緩，具有西部高而東部低的地勢(shì)特點(diǎn)?；镜孛差愋统媳边吘壣降貫榇?、中起伏的高山和極高山外，廣大地區(qū)主要為中小起伏的山高和高海拔丘陵、臺(tái)地和平原。山地起伏和緩，河谷盆地寬坦，是青藏高原上高原面保存最完整的地區(qū)[20]。

圖1 可可西里地區(qū)位置及其境內(nèi)湖泊分布

1.2 數(shù)據(jù)來源

可可西里地區(qū)湖泊眾多，僅面積大于10km2的湖泊有83個(gè)，且多為咸水湖—半咸水湖，淡水湖和鹽湖分布較少[16]。將《中國(guó)湖泊志》統(tǒng)計(jì)的湖泊信息作為1993年的水位參考資料[21]，利用可可西里地區(qū)湖泊矢量數(shù)據(jù)[16]，通過美國(guó)國(guó)家冰雪數(shù)據(jù)中心(NSIDC，National Snow and Ice Data Center)發(fā)布的ICESat/GLAS陸地測(cè)高數(shù)據(jù)產(chǎn)品(GLA14)提取2003—2009年湖泊水位數(shù)據(jù)，獲得可可西里地區(qū)有ICESat/GLAS數(shù)據(jù)的40個(gè)面積大于10km2的湖泊水位數(shù)據(jù)，表1給出了可可西里地區(qū)湖泊水位數(shù)據(jù)獲得時(shí)間。此外，將0.5°插值年均氣溫、年降水量數(shù)據(jù)和距離可可西里地區(qū)最近的伍道梁、改則和獅泉河氣象臺(tái)站資料作為研究湖泊變化背景參考，該數(shù)據(jù)均由中國(guó)氣象數(shù)據(jù)網(wǎng)提供(http：//data.cma.cn)。

表1 可可西里地區(qū)主要湖泊ICESat/GLAS數(shù)據(jù)獲取時(shí)間序列

湖泊名稱ICESat/GLAS數(shù)據(jù)獲取時(shí)間(20??年??月??日)湖泊名稱ICESat/GLAS數(shù)據(jù)獲取時(shí)間(20??年??月??日)振泉湖03-10-06玉液湖07-10-12飲馬湖03-10-13浩波湖03-10-14永紅湖03-10-05黑石北湖03-10-24雙連湖04-11-07葫蘆池03-10-06太陽湖03-10-13美馬錯(cuò)03-11-18豌豆湖03-10-13阿魯錯(cuò)03-11-18可考湖04-05-23, 06-10-30, 07-03-16鹽湖04-10-28, 05-03-15, 07-04-05漣湖03-11-07, 04-10-26, 05-11-1206-11-16, 07-10-24, 08-12-06茍魯錯(cuò)03-10-29, 04-10-17, 05-11-0307-10-15, 08-10-17, 09-03-21銀波湖04-10-05, 05-05-22, 06-10-2607-10-04, 08-10-05, 09-10-02西金烏蘭湖03-11-11, 04-10-29, 05-11-1606-11-19, 07-10-28, 08-12-09月亮湖03-11-11, 04-10-29, 05-11-1606-06-18, 07-10-28, 08-12-09烏蘭烏拉湖03-11-07, 04-10-29, 05-11-1206-11-16, 07-10-28, 08-12-09向陽湖03-10-30, 04-03-02, 05-11-0406-03-08, 07-10-16雀莫錯(cuò)03-11-15, 04-11-03, 05-11-2006-11-24, 08-12-14, 09-04-07移山湖03-10-17, 04-05-18, 05-10-2207-10-03, 08-10-04, 09-10-01瑪章錯(cuò)欽04-10-12, 06-11-02, 07-10-1108-02-25, 09-10-09仙鶴湖04-10-05, 05-05-22, 06-10-2607-10-04, 08-10-05, 09-10-02庫賽湖03-03-02, 03-03-06, 03-03-1003-03-18, 03-09-30, 03-10-04江尼茶卡04-11-06, 05-03-24, 06-03-2707-04-14, 08-03-21, 09-04-11勒斜武擔(dān)湖03-11-03, 04-10-26, 05-11-1206-11-16, 07-10-24, 08-12-06節(jié)約湖03-11-07, 04-10-26, 05-11-1206-11-16, 07-10-24, 08-12-06庫水浣03-11-07, 04-10-26, 05-11-1206-11-16, 07-10-24, 08-12-06錯(cuò)尼03-11-12, 04-03-14, 05-11-1706-11-, 07-10-29, 08-12-10瑪爾蓋茶卡03-02-24, 03-03-12, 03-03-2003-03-28, 03-09-28朝陽湖03-11-08, 04-10-27, 05-11-1306-11-17, 07-10-25, 08-12-07, 09-03-31島湖03-02-25, 03-03-05, 03-03-1303-03-21, 03-09-29得雨湖03-11-08, 04-10-27, 05-11-1306-11-17, 07-10-25, 08-12-07, 09-03-31涌波湖03-10-02, 04-10-22, 05-11-0906-11-12, 07-10-21, 08-12-02, 09-03-27高臺(tái)湖03-11-15, 04-11-03, 05-11-2006-11-24, 07-04-10, 08-12-14, 09-04-07可可西里湖03-10-17, 04-10-04, 05-10-22,06-10-25, 07-10-03, 08-10-04, 09-10-01多格錯(cuò)仁03-10-26, 04-10-13, 05-10-2706-10-31, 07-10-08, 08-10-10, 09-10-10瑪日阿措04-05-23, 04-10-09, 05-10-2606-10-30, 07-10-07, 08-10-09, 09-03-13

1.3 研究方法

將ICESat/GLA14陸地地表測(cè)高數(shù)據(jù)進(jìn)行格式轉(zhuǎn)化，生成文本格式的GPS點(diǎn)位數(shù)據(jù)，并根據(jù)點(diǎn)位數(shù)據(jù)的經(jīng)緯度坐標(biāo)生成矢量點(diǎn)圖層。由于現(xiàn)有大部分水位數(shù)據(jù)基準(zhǔn)一般都采用平均海平面高程，所以該文將ICESat/GLAS的基于EMG96的大地水準(zhǔn)面高程轉(zhuǎn)換為基于平均海平面的高程，這樣以便以后與實(shí)測(cè)水位數(shù)據(jù)進(jìn)行對(duì)比驗(yàn)證。

獲得較高精度的湖泊水面高程信息，需要排除湖岸等因素的影響。為此，該文將2003—2009年可可西里地區(qū)湖泊邊界數(shù)據(jù)向湖心方向作一定距離的緩沖區(qū)(湖泊面積小于50km2緩沖區(qū)距離為200m，面積50～100km2緩沖區(qū)距離為500m，湖泊面積大于100km2緩沖區(qū)距離為1km)，取高程均值作為當(dāng)天的水位信息； 根據(jù)GLAS高程點(diǎn)的高程屬性信息，在ArcGIS 9.3軟件平臺(tái)下剔除被云、波浪等干擾的高程點(diǎn)，并根據(jù)單個(gè)觀測(cè)日內(nèi)最大/最小水位、標(biāo)準(zhǔn)方差，刪除當(dāng)日水位異常的高程點(diǎn)，使高程點(diǎn)位的方差小于0.1，按照觀測(cè)日期生成每個(gè)湖泊水位數(shù)據(jù)。從圖2可以看出，ICESat/GLAS衛(wèi)星測(cè)高點(diǎn)數(shù)據(jù)經(jīng)過美馬錯(cuò)和阿魯錯(cuò)，通過上述方法可以得出湖泊當(dāng)時(shí)水位數(shù)據(jù)。

圖2 ICESat衛(wèi)星過境時(shí)湖泊水位高程測(cè)量

2 結(jié)果與分析

2.1 可可西里地區(qū)主要湖泊水位統(tǒng)計(jì)

早在20世紀(jì)50年代，我國(guó)即建立了湖泊科學(xué)研究機(jī)構(gòu)，通過廣大湖泊科技工作者50多年來的努力工作，取得了豐碩的科研成果，出版了《中國(guó)湖泊志》[20]等多部科學(xué)專著，但都多缺乏系統(tǒng)、綜合性、動(dòng)態(tài)的水位研究工作。表2列出了ICESat/GLAS衛(wèi)星測(cè)高數(shù)據(jù)得到的可可西里地區(qū)主要湖泊水位數(shù)據(jù)，該數(shù)據(jù)可作為《中國(guó)湖泊志》水位數(shù)據(jù)驗(yàn)證數(shù)據(jù)。該研究的可可西里地區(qū)40個(gè)湖泊中，黑石北湖水位最高(5 050.83m)，鹽湖水位最低(4 440.05m)，可可西里地區(qū)湖泊水位與各湖泊所在湖盆海拔有關(guān)。

表2 ICESat數(shù)據(jù)獲得的各湖泊水位數(shù)據(jù)

湖泊名稱海拔(m)日期湖泊名稱海拔(m)日期烏蘭烏拉湖4 858.702008-12-09振泉湖4 796.352003-10-06多格錯(cuò)仁4 820.782009-10-10永紅湖4 773.752003-10-05西金烏蘭湖4 774.212008-12-09可考湖4 888.172007-03-16可可西里湖4 889.552009-10-01涌波湖4 886.522009-03-27庫賽湖4 476.182003-10-04得雨湖4 850.472009-03-31勒斜武擔(dān)湖4 873.332008-12-06島湖4 926.282003-09-29美馬錯(cuò)4 923.742003-11-18銀波湖4 884.842009-10-02瑪爾蓋茶卡4 798.452003-09-28仙鶴湖4 843.872009-10-02錯(cuò)尼4 926.682008-12-10漣湖4 910.362008-12-06玉液湖4 861.852007-10-12庫水浣5 005.302008-12-06飲馬湖4 911.012003-10-13雙連湖4 818.052004-11-07阿魯錯(cuò)4 937.262003-11-18月亮湖4 910.472008-12-09太陽湖4 881.612003-10-13移山湖4 845.292009-10-01向陽湖4 865.862007-10-16茍魯錯(cuò)4 665.712009-03-21黑石北湖5 050.832003-10-24浩波湖4 829.562003-10-14雀莫錯(cuò)4 924.632007-10-11豌豆湖4 858.152003-10-13葫蘆池4 798.9072003-10-06節(jié)約湖4 818.862008-12-06朝陽湖4 745.312009-03-31瑪日阿措4 648.782009-03-13瑪章錯(cuò)欽4 676.152009-10-09高臺(tái)湖5 015.712009-04-07鹽湖4 440.052007-04-05江尼茶卡4 789.582009-04-11

2.2 可可西里地區(qū)湖泊水位變化

從湖泊水位波動(dòng)幅度來看(圖3)，錯(cuò)尼湖水位波動(dòng)幅度最大，共計(jì)上升24.68m，在可可西里地區(qū)有23個(gè)湖泊水位波動(dòng)在0～5m， 9個(gè)湖泊水位波動(dòng)幅度在5～10m之間， 7個(gè)湖泊水位波動(dòng)幅度在10～15m之間。依據(jù)可可西里地區(qū)湖泊面積數(shù)據(jù)[16]， 1990—2009年，島湖面積擴(kuò)張接近60倍，瑪爾蓋茶卡面積擴(kuò)張6倍多，振泉湖和得雨湖面積擴(kuò)張接近5倍，錯(cuò)尼湖面積擴(kuò)張接近3倍，由此可以得出，可可西里地區(qū)湖泊水位波動(dòng)超過10m的湖泊面積都呈倍數(shù)級(jí)增加。從圖3可以看出，可可西里水位波動(dòng)較大的湖泊均集中于可可西里中部地區(qū)，東西部湖泊水位波動(dòng)并不大。

圖3 可可西里地區(qū)主要湖泊水位變化幅度

為動(dòng)態(tài)監(jiān)測(cè)可可西里地區(qū)湖泊，分別計(jì)算出1993—2009年區(qū)內(nèi)湖泊水位變化趨勢(shì)(圖4)， 21個(gè)湖泊水位有上升， 7個(gè)湖泊水位先降低后升高， 5個(gè)湖泊水位先升高后降低?？煽晌骼锏貐^(qū)湖泊雖然有很大一部分無數(shù)據(jù)或數(shù)據(jù)不全，且由于水位數(shù)據(jù)和面積數(shù)據(jù)時(shí)間段不一致，例如水位變化狀態(tài)為先升高后降低、降低的12個(gè)湖泊中， 3個(gè)湖泊水位采集時(shí)間為枯水期， 9個(gè)湖泊數(shù)據(jù)采集時(shí)間為豐水期，但僅有2003年一期數(shù)據(jù)，而在1970—2011年整體上湖泊面積是擴(kuò)張的[16]，兩者并不矛盾。由于數(shù)據(jù)限制，湖泊水位年內(nèi)變化尚不清楚，初步推斷每年3月份為湖泊水位最低點(diǎn)， 11月份為湖泊水位最高點(diǎn)。

圖4 可可西里地區(qū)主要湖泊水位變化趨勢(shì)

圖5 可可西里地區(qū)降水、氣溫變化

圖6 1970—2011年伍道梁、改則和獅泉河3個(gè)氣象臺(tái)站年潛在蒸散發(fā)變化

2.3 可可西里地區(qū)湖泊水位變化氣候因素分析

1970—2013年，可可西里地區(qū)氣候發(fā)生明顯變化，由0.5°×0.5°插值降水、氣溫?cái)?shù)據(jù)可知(圖5)，可可西里地區(qū)降水呈增加趨勢(shì)，氣候傾向率高達(dá)21.72mm/10年。參考其他學(xué)者對(duì)可可西里地區(qū)湖泊面積變化的研究，可可西里地區(qū)湖泊面積變化趨勢(shì)與0.5°×0.5°插值降水?dāng)?shù)據(jù)的降水變化具有很好的一致性， 1995年之后降水量迅速增加，而此時(shí)的湖泊面積規(guī)模也迅速擴(kuò)張[16]。尤其是在2002年，降水增加至歷史時(shí)期最高值，為343.68mm，這是引起可可西里地區(qū)湖泊水位上升的首要?dú)夂蛞蛩亍?/p>

1970—2013年，可可西里地區(qū)氣溫上升明顯(圖5)，氣候傾向率為0.45℃/10年，但氣溫增加并未引起湖泊蒸發(fā)量增加。鑒于可可西里地區(qū)湖泊多為內(nèi)陸湖，又由于該區(qū)域缺少相關(guān)實(shí)測(cè)資料，在不計(jì)其他因素的條件下，湖泊水量流失主要取決于湖面蒸發(fā)，且湖面蒸發(fā)與由氣象臺(tái)站資料得到的潛在蒸發(fā)具有很好的線性關(guān)系[22]，因此該文利用Penman-Monteith的方法分別計(jì)算伍道梁、獅泉河和改則3個(gè)氣象臺(tái)站的年潛在蒸散發(fā)量及其變化情況(圖6)。顯然，1971—2011年，改則、獅泉河和伍道梁3個(gè)氣象臺(tái)站的年潛在蒸散發(fā)總量都呈減少趨勢(shì)，減少速率分別為-25mm/10年、-18mm/10年和-11mm/10年。進(jìn)一步分析發(fā)現(xiàn)，這3個(gè)氣象臺(tái)站2000—2011年平均潛在蒸散發(fā)量均小于研究時(shí)段內(nèi)多年平均值，表明2000—2011年研究區(qū)潛在蒸散發(fā)整體上下降了一個(gè)臺(tái)階，這是可可西里地區(qū)湖泊水位上升的另一氣候因素。

3 討論

該研究以可可西里地區(qū)為研究區(qū)，利用ICESat/GLAS衛(wèi)星測(cè)高數(shù)據(jù)對(duì)其境內(nèi)有衛(wèi)星過境的40個(gè)湖泊的水位進(jìn)行測(cè)算，并對(duì)各個(gè)湖泊水位變化幅度和變化規(guī)律進(jìn)行詳細(xì)闡述。該研究將《中國(guó)湖泊志》中的湖泊水位與ICESat/GLAS衛(wèi)星測(cè)高數(shù)據(jù)放在同一的參考系統(tǒng)內(nèi)，并采用了緩沖區(qū)算法和水位均值算法，消除了ICESat/GLAS衛(wèi)星測(cè)高數(shù)據(jù)中湖泊岸堤影響和浪涌影響。詳細(xì)列出了可可西里地區(qū)ICESat/GLAS衛(wèi)星過境的40個(gè)湖泊水位數(shù)據(jù)，其中有7個(gè)湖泊2003年水位較《中國(guó)湖泊志》水位上升超過10m，錯(cuò)尼湖水位上升超過20m，變化超出湖泊一般變化規(guī)律，依據(jù)姚曉軍[16]等的研究成果，可可西里地區(qū)湖泊面積在1990年達(dá)到歷史最低值， 2000—2002年3年期間湖泊面積陡然上升，湖泊擴(kuò)張總面積超過20%，表明該研究數(shù)據(jù)可靠，并可為《中國(guó)湖泊志》數(shù)據(jù)進(jìn)行精度驗(yàn)證。近年來，可可西里地區(qū)主要湖泊水位上升明顯，這與該地區(qū)湖泊面積擴(kuò)張現(xiàn)象所一致[16]，也符合青藏高原湖泊面積擴(kuò)張的整體趨勢(shì)[23]。

就湖泊變化的驅(qū)動(dòng)因素而言，利用0.5°插值氣溫、降水?dāng)?shù)據(jù)并結(jié)合研究區(qū)周邊3個(gè)氣象臺(tái)站數(shù)據(jù)，能夠很好分析出可可西里地區(qū)湖泊水位整體上升的主要影響因素為降水增加和蒸發(fā)減少。除降水補(bǔ)給之外，冰川融水和凍土中的水分釋放也是青藏高原湖泊補(bǔ)給的主要形式[24， 25]。根據(jù)中國(guó)冰川目錄[26]和最新一期冰川目錄資料[27]可知，在20世紀(jì)70年代初期，可可西里地區(qū)冰川面積為2 423.60km2，到2000年以后減少為2 243.35km2，共減少180.25km2，表明在可可西里地區(qū)，冰川整體上呈退縮狀態(tài)。然而在整個(gè)研究區(qū)，太陽湖距馬蘭冰帽距離較近(約2.13km)，但其面積變化不大，其余湖泊距冰川均較遠(yuǎn)，若考慮冰川融水沿途水分蒸發(fā)與下滲，以及部分湖泊并無冰川融水注入，冰川融水對(duì)湖泊的補(bǔ)給量應(yīng)該十分有限。就凍土中水分釋放而言，伍道梁地區(qū)凍土變化自1980年以來凍土厚度變薄，活動(dòng)層增厚，凍結(jié)持續(xù)天數(shù)縮短[28]，但凍土層中地下冰融化很難轉(zhuǎn)入可流動(dòng)的地下水，且比例很小[29]。

此外，可可西里地區(qū)存在地下水補(bǔ)給，且補(bǔ)給量較大[20]，如Yao等[30]通過對(duì)湖泊冰情的研究發(fā)現(xiàn)可可西里地區(qū)存在很多地下溫泉，溫泉在夏季不易發(fā)現(xiàn)但在冬季影響湖泊結(jié)冰且較為明顯。從圖3可以看出，錯(cuò)尼湖泊水位波動(dòng)幅度最大，共計(jì)上升24.68，研究發(fā)現(xiàn)錯(cuò)尼水溫垂直分布呈“S”型，水深0～15m水溫逐漸下降， 15m以下逐漸升高， 30～35m達(dá)16.3～17.3°C，并據(jù)此認(rèn)為該湖湖底存在溫泉補(bǔ)給[20]。因此，兩湖連接處偏東的部位應(yīng)是該溫泉所處位置，同時(shí)，這也表明區(qū)域泉水活動(dòng)對(duì)于湖泊的空間演化也起著重要作用[30]。綜上所述，在研究時(shí)段內(nèi)，降水增加導(dǎo)致湖泊輸入增多和蒸發(fā)減少導(dǎo)致湖泊輸出減少應(yīng)該是可可西里地區(qū)湖泊水位發(fā)生變化的主要原因，冰川融水和凍土中的水分釋放以及區(qū)內(nèi)地下水補(bǔ)給形式可能是引起可可西里地區(qū)湖泊水位上升的原因，但應(yīng)非決定性因素。

由于可可西里地區(qū)為無人區(qū)，缺少相應(yīng)實(shí)測(cè)資料，該研究未能從水量平衡角度分析湖泊水位變化，未能給出降水、蒸發(fā)、下滲、地下水等補(bǔ)給參數(shù)的量值； 此外，受數(shù)據(jù)限制，ICESat/GLAS衛(wèi)星測(cè)高數(shù)據(jù)年限僅為2003—2009年，并且部分湖泊該時(shí)間序列也并不完整，該研究未對(duì)個(gè)湖泊水位變化做年內(nèi)分析，今后將利用多源數(shù)據(jù)進(jìn)行補(bǔ)充。

4 結(jié)論

(1)1993—2009年，可可西里地區(qū)錯(cuò)尼湖泊水位上升最高，為24.68m，銀波湖下降最為明顯，下降6.16m。在可可西里地區(qū)中部湖泊水位波動(dòng)較大，整個(gè)區(qū)域大部分湖泊水位波動(dòng)在0～5m之間，有7個(gè)湖泊水位波動(dòng)幅度在10～15m之間，應(yīng)該引起人們重點(diǎn)關(guān)注。

(2)1993—2009年，可可西里地區(qū)湖泊水位整體上呈上升趨勢(shì)，但呈現(xiàn)出一定的區(qū)域差異。湖泊水位上升的湖泊最多，約占湖泊總數(shù)量的60%以上，分布于整個(gè)可可西里地區(qū)。水位下降的湖泊數(shù)量較少，零星分布在整個(gè)區(qū)域。

(3)受資料限制，該文尚無法從水量平衡角度給出可可西里地區(qū)湖泊水量變化及各輸入輸出要素變化量值，初步分析認(rèn)為可可西里地區(qū)湖泊動(dòng)態(tài)變化主要與氣候變化有關(guān)，其中降水增加(或減少)和湖面蒸發(fā)減少(或增加)是導(dǎo)致湖泊水位上升(或下降)的決定性因素，而氣候變暖引起的冰川融水增加、凍土水分釋放以及區(qū)內(nèi)地下水補(bǔ)給形式可能是次要原因。