李龍 姚曉軍 巨喜鋒 卜亞平 李曉鋒

摘?要：利用數(shù)字高程模型并結(jié)合中國流域劃分方案，得到中國內(nèi)流區(qū)一級流域分布范圍，基于2000—2014年Landsat TM、ETM+遙感影像，通過人工目視解譯提取中國內(nèi)流區(qū)面積大于100 km2的湖泊邊界矢量數(shù)據(jù)，并據(jù)此對中國內(nèi)流區(qū)主要湖泊的分布、面積變化及成因進行分析。結(jié)果表明：①在中國現(xiàn)有的7個一級內(nèi)流區(qū)中，面積大于100 km2的湖泊有88個，總面積32 018.46 km2，除伊犁河流域外均有分布;②2000—2014年，中國內(nèi)流區(qū)湖泊總面積以平均每年231.95 km2的速度持續(xù)增大，總計增加3 247.33 km2，其中柴達木盆地和青藏高原腹地湖泊面積持續(xù)增大，塔里木盆地和準噶爾盆地湖泊面積在2003年之后持續(xù)萎縮，陰山北部和祁連山北部湖泊在2003年以后面積平穩(wěn)增大;③在研究時段內(nèi)，絕大部分湖泊面積增大明顯，主要集中在可可西里地區(qū)和羌塘高原地區(qū)，面積減小區(qū)域零星散布于整個研究區(qū);④降水量增加、蒸發(fā)量減少、冰川和凍土加速消融是導致研究區(qū)湖泊面積整體增加的主要原因，局地氣候、各湖泊之間存在的水道聯(lián)系與人類活動的干擾造成了部分湖泊面積減小和波動。

關(guān)鍵詞：湖泊;面積變化;氣候變化;中國內(nèi)流區(qū)

中圖分類號：TV213.3文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2020.06.013

Variations of Lakes Area in Endorheic Basin of China During 2000-2014

LI Long1，YAO Xiaojun2， JU Xifeng3， BU Yaping1， LI Xiaofeng2

（1.Lanzhou Resources & Environment Voc-Tech College， Lanzhou 730721， China;

2.College of Geography and Environment Sciences， Northwest Normal University， Lanzhou 730070， China;

3.School of Resources & Environment Science， Xinjiang University， Urumqi 830046， China）

Abstract：The distribution range of first-level watershed in endorheic basin of China was gotten by combing digital elevation model and drainage basin partition plan of China. Besides， based on Landsat TM/ETM+/OLI remote sensing images in the period of 2000-2014， the vectorized boundary of 88 lakes with areas greater than 100 km2 were obtained by artificial visual interpretation technology and then the causes of lake variation were analyzed. The results show that： a） Among seven primary endorheic basins of China， there are 88 lakes that area greater than 100 km2 and the total area is 32，018.46 km2 except Ili River watershed. b） During the period of 2000-2014， the area of endorheic basin of China expends with the speed of 231.95 km2 per year， which totally increases to 2，347.33 km2. The lake area of Qaidam Basin and Qinghai-Tibet Plateau hinterland is continuously expending. Although the areas of Tarim Basin and Junggar Basin are shrinking after 2003， the lakes of the north of Yinshan and Qilian Mountains are increased smoothly. c） During researching， majority of lakes expend obviously， most of them concentrated on Hohxil and Qiangtang Plateau. The shrinking area scatters throughout the study area. d） The main reason why lakes expending are because of increasing precipitation， decreasing evaporation and melting of glaciers and frozen oil. Local climate， waterway of lakes and human activities also affect the reduction and fluctuation of some lakes.

Key words： lake; area change; climate change; endorheic basin of China

湖泊作為陸地水圈的組成部分，參與自然界的水分循環(huán)，其水域變化是其所在流域水量平衡的綜合結(jié)果[1]，對氣候的變化極為敏感[2-3]。據(jù)統(tǒng)計，全國共有面積大于1.0 km2的湖泊2 693個，總面積81 414.6 km2[4]。近年來，伴隨著湖區(qū)荒漠化、生態(tài)退化、自然災害頻發(fā)等問題日益嚴重，湖泊研究越來越受到人們的關(guān)注[5-6]。很多學者對湖泊進行了較為系統(tǒng)的研究，主要集中在長時間序列湖泊面積變化、湖泊水位變化和湖泊的水量平衡[7-9]，研究區(qū)域涉及青藏高原湖區(qū)、西北干旱湖區(qū)、云貴高原湖區(qū)、東北平原與山地湖區(qū)、東部平原湖區(qū)等典型區(qū)域[10-14]。

內(nèi)流區(qū)是指地表徑流不與海洋相通的地區(qū)[15]，按照一級流域劃分方案，全國現(xiàn)有7個一級內(nèi)流區(qū)[16]，除伊犁河流域外均有面積大于100 km2的湖泊分布[2， 4]，主要包括青藏高原湖區(qū)和西北干旱區(qū)湖區(qū)[2， 16]，湖泊水資源一直是影響區(qū)域經(jīng)濟發(fā)展的關(guān)鍵因素[17]。在內(nèi)流區(qū)的大部分湖泊，只有入湖流量補給，而無流出，更能反映流域內(nèi)氣候變化，湖水損耗主要是湖面蒸發(fā)和人類活動[2]。由于研究目的不同，國內(nèi)學者尚未對我國內(nèi)流區(qū)湖泊動態(tài)變化進行詳細研究，因此缺少內(nèi)流區(qū)間的湖泊變化差異對比分析，亟需應(yīng)用遙感資料和地理信息技術(shù)建立我國內(nèi)流區(qū)2000年以來較長時間序列的湖泊空間數(shù)據(jù)庫，明晰我國不同內(nèi)流區(qū)湖泊變化差異及影響因素，從而為我國內(nèi)流區(qū)生態(tài)、經(jīng)濟建設(shè)乃至跨流域調(diào)水提供數(shù)據(jù)支持和理論依據(jù)。

1?研究區(qū)概況

我國內(nèi)、外流區(qū)大致沿大興安嶺西麓—內(nèi)蒙古高原南緣—陰山—賀蘭山—祁連山—日月山—巴顏喀拉山—念青唐古拉山—岡底斯山一線分布，此線以西，除額爾齊斯河流入北冰洋外，基本上屬于內(nèi)流區(qū)[4]。按照流域劃分體系，我國一級內(nèi)流區(qū)主要包括青藏高原腹地、柴達木盆地、準噶爾盆地、塔里木盆地、陰山北部、祁連山北部和伊犁河內(nèi)流區(qū)（見圖1）[16]。內(nèi)流區(qū)氣候干旱、降水稀少，區(qū)內(nèi)湖泊以鹽湖和咸水湖為主，淡水湖較少，有冰雪融水補給，來水較豐，因此湖泊具有較大的面積，深度也較大[2]?？紤]到湖泊受季節(jié)影響較為強烈[18]，湖泊面積越大越穩(wěn)定[17]，因此選取內(nèi)流區(qū)內(nèi)所有面積大于100 km2的88個湖泊作為研究對象，區(qū)內(nèi)除伊犁河流域外，其他流域均有此規(guī)模的湖泊分布，基本能夠代表我國內(nèi)流區(qū)湖泊變化特征。

2?數(shù)據(jù)和方法

提取中國一級內(nèi)流區(qū)范圍，采用地理空間數(shù)據(jù)云（http：//www.gscloud.cn）提供的空間分辨率為90 m的SRTM（Shuttle Radar Topography Mission航天飛機雷達地形測繪任務(wù)），參考《中國湖泊志》確定研究區(qū)湖泊位置及名稱[2]。湖泊矢量基礎(chǔ)數(shù)據(jù)參考“中國2005—2006年1∶25萬面積1.0 km2以上湖泊分布數(shù)據(jù)集”[19]、“2000—2011年可可西里地區(qū)主要湖泊邊界矢量數(shù)據(jù)”[7]、“2000—2014年中國西北地區(qū)湖泊邊界矢量數(shù)據(jù)”[18]。此外，于USGS網(wǎng)站（http：//earthexplorer.usgs.gov）下載317景Landsat TM/ETM+遙感影像，用以補充解譯缺失年份的湖泊矢量數(shù)據(jù)，從而形成“2000—2014年中國內(nèi)流區(qū)湖泊矢量數(shù)據(jù)集”。為保證補充解譯的湖泊矢量數(shù)據(jù)與參考數(shù)據(jù)的統(tǒng)一，用于補充解譯的影像時間集中于湖泊狀態(tài)較為穩(wěn)定的10月，僅有21景影像采集時段為9月和11月。

基于SRTM地形數(shù)據(jù)，利用ArcGIS 9.3水文模型提取中國一級內(nèi)流區(qū)范圍，結(jié)合Google Earth軟件以及Ma等[4]和張國平等[16]的研究成果，對提取結(jié)果進行驗證，確保提取的流域與實際湖泊水系特征相吻合。湖泊邊界矢量通過人工目視解譯獲取，解譯標準參照中國國家科技基礎(chǔ)性工作專項“中國湖泊水質(zhì)、水量與生物資源調(diào)查”制定的原則[4]，解譯精度控制在一個像元以內(nèi)。

3?結(jié)果與討論

3.1?中國內(nèi)流區(qū)湖泊現(xiàn)狀

根據(jù)流域劃分方案[4，16]，中國共有29個一級流域，其中內(nèi)流區(qū)2014年面積大于100 km2的湖泊共88個（見表1），總面積32 018.46 ?km2。青藏高原腹地湖泊數(shù)量最多、面積最大，分別占內(nèi)流區(qū)湖泊總數(shù)的77%、總面積的69%;其次是柴達木盆地，分別占湖泊總數(shù)的14%、總面積的20%;塔里木盆地和祁連山北部僅有博斯騰湖和東居延海面積大于100 km2，伊犁河流域無大于100 km2的湖泊。

3.2?中國內(nèi)流區(qū)湖泊面積變化特征

2000—2014年中國內(nèi)流區(qū)湖泊面積持續(xù)增大（見圖2）。2000年內(nèi)流區(qū)88個面積大于100 km2的湖泊總面積為28 771.13 km2，2014年為32 018.46 km2，14 a增加了3 247.33 km2，平均每年增加231.95 km2。2000—2002年與2009—2012年內(nèi)流區(qū)湖泊總面積增加最為明顯，均超過了1 000 km2。

3.3?中國內(nèi)流區(qū)各流域湖泊變化特征

由圖3可知：2000—2014年青藏高原腹地和柴達木盆地湖泊面積呈增大趨勢，2014年比2000年分別增加3 336.19、311.28 km2，增速分別為222.41、20.75 km2/a;塔里木盆地和準噶爾盆地湖泊面積在2002年、2003年最大，之后呈萎縮趨勢，2014年比2000年分別減小216.32、228.74 km2;陰山北部和祁連山北部湖泊面積在2003年最小，之后呈增大趨勢，2014年陰山北部湖泊面積比2000年減小2.79 km2，而祁連山北部湖泊面積增大47.71 km2。

3.4?中國內(nèi)流區(qū)湖泊面積變化空間差異

2000—2014中國內(nèi)流區(qū)主要湖泊面積變化及分布情況見圖4，可以看出，各湖泊面積變化存在一定的差異性。把湖泊面積變化分為增大、減小、先增大后減小、先減小后增大4種情況，其中：面積增大的湖泊最多，有61個，全部位于青藏高原腹地和柴達木盆地;其次為面積減小和面積先減小后增大兩種湖泊，均有11個，散布在整個內(nèi)流區(qū);面積先增大后減小的湖泊有5個，分別為青藏高原腹地的拜惹布錯和班戈錯、準噶爾盆地的艾比湖和賽里木湖、柴達木盆地的卓乃湖。

為比較湖泊面積變化的差異性，分別計算2014年與2000年各湖泊面積變化百分比（見圖5）。其中：2014年72個湖泊面積比2000年有所增大，占內(nèi)流區(qū)湖泊數(shù)量的82%。卓乃湖和艾比湖面積減幅接近40%，博斯騰湖等14個湖泊面積減幅小于30%，青海湖和納木錯等38個湖泊面積增幅小于20%，28個湖泊面積增幅為20%～80%，6個湖泊面積成倍增大。從空間分布來看，2014年比2000年面積明顯增大的湖泊主要集中在可可西里地區(qū)和羌塘高原地區(qū)，面積減小的湖泊則散布于整個內(nèi)流區(qū)。

3.5?中國內(nèi)流區(qū)湖泊變化驅(qū)動因素分析

影響中國內(nèi)流區(qū)湖泊面積變化的主要因素有地表徑流、湖面降水、冰川融水、湖面蒸發(fā)和地下水交換、凍土、沼澤等[9]。2000—2014年，全球氣溫升高明顯加快[20]，中國西北地區(qū)由暖干逐漸向暖濕轉(zhuǎn)化，冰川退化嚴重[21-22]。青藏高原腹地近年來降水增加明顯，蒸發(fā)顯著減少，冰川融水增加，湖泊面積增大主要因素為有利的氣候條件[23-26]。柴達木盆地近年來處于豐水期，蒸發(fā)減少，降水量增加，冰雪融水增多，導致該區(qū)2001年以來湖泊面積顯著增大[27-28]。內(nèi)蒙古地區(qū)氣溫升高，蒸發(fā)量顯著增大，降水量明顯減小，這種不利的氣候條件導致區(qū)內(nèi)湖泊面積萎縮[29]。2000—2014年，祁連山北部地區(qū)和陰山北部地區(qū)降水量先增大后減小再增大[30-31]，與該區(qū)湖泊面積變化規(guī)律一致。祁連山北部地區(qū)2000年以來氣溫降低，蒸發(fā)量減少，湖泊面積略有增大[30];陰山北部自2000年后氣溫升高明顯、蒸發(fā)量呈先增大再減小后增大的趨勢[31-32]，湖泊面積略有減小。1961—2013年，準噶爾盆地降水量呈增大趨勢，蒸發(fā)量顯著減小，湖泊面積增大主要因素為有利的氣候條件[26， 33]。塔里木盆地降水量略有增大，蒸發(fā)量顯著減小，區(qū)內(nèi)除博斯騰湖外其余湖泊面積均呈增大趨勢[18， 34]。

2014年中國內(nèi)流區(qū)共有16個湖泊面積較2000年減小，原因與各湖泊之間存在水道聯(lián)系、局地氣候和人類活動的干擾有關(guān)，如色林錯流域的班戈錯、錯鄂、吳如錯、格仁錯與色林錯存在水道聯(lián)系[35]，西昆侖地區(qū)阿魯錯與美馬錯亦存在水道聯(lián)系[7]，卓乃湖、海丁諾爾、鹽湖和庫賽湖的水道聯(lián)系引起湖水串聯(lián)性外溢事件[6]，郭扎錯、太陽湖、飲馬湖流域局地氣候變化引起冰川補給變化[7， 36]，柴達木盆地茶卡鹽湖、準噶爾盆地艾比湖、塔里木盆地博斯騰湖、青藏高原腹地瑪旁雍錯、陰山北部達來諾爾地區(qū)人類用水供需矛盾加劇[37-39]。

4?結(jié)?論

（1）中國現(xiàn)有7個一級內(nèi)流區(qū)，2014年區(qū)內(nèi)大于100 km2的湖泊有88個，總面積32 018.46 ?km2，除伊犁河流域外均有分布，其中青藏高原腹地湖泊數(shù)量最多、面積最大，分別占內(nèi)流區(qū)湖泊總數(shù)的77%、總面積的69%。

（2）2000—2014年，中國內(nèi)流區(qū)湖泊面積呈增大趨勢，共增加3 247.33 km2，平均每年增加231.95 km2。柴達木盆地和青藏高原腹地湖泊面積持續(xù)增大，塔里木盆地和準噶爾盆地湖泊面積2003年之后持續(xù)萎縮，陰山北部和祁連山北部湖泊面積2003年之后平穩(wěn)增大。湖泊面積增大明顯區(qū)域主要集中在可可西里地區(qū)和羌塘高原地區(qū)，面積減小區(qū)域散布于整個內(nèi)流區(qū)。

（4）2000—2014年，中國內(nèi)流區(qū)降水量增大、蒸發(fā)量減小、冰川加速退縮，導致區(qū)內(nèi)湖泊面積整體增大，個別湖泊面積減小或波動與局地氣候、湖泊之間存在水道聯(lián)系、人類活動干擾有關(guān)。

參考文獻：

[1]?丁永建，劉時銀，葉柏生，等. 近50 a中國寒區(qū)與旱區(qū)湖泊變化的氣候因素分析[J]. 冰川凍土，2006，28（5）：623-632.

[2]?王蘇民，竇鴻身.中國湖泊志[M].北京：科學出版社，1998：5-77.

[3]?魯安新，姚檀棟，王麗紅，等.青藏高原典型冰川和湖泊變化遙感研究[J].冰川凍土，2005，27（6）：783-792.

[4]?MA R，YANG G，DUAN H，et al. Chinas Lakes at Present：Number，Area and Spatial Distribution[J]. Science China Earth Sciences，2011，54（2）：283-289.

[5]?張振克，楊達源.中國西北干旱區(qū)湖泊水資源-環(huán)境問題與對策[J].干旱區(qū)資源與環(huán)境，2001，15（2）：7-10.

[6]?姚曉軍，孫美平，宮鵬，等.可可西里鹽湖湖水外溢可能性初探[J].地理學報，2016，71（9）：1520-1527.

[7]?姚曉軍，劉時銀，李龍，等.近40年可可西里地區(qū)湖泊時空變化特征[J].地理學報，2013，68（7）：886-896.

[8]?葛莉，習曉環(huán)，王成，等.ICESat-1/GLAS數(shù)據(jù)湖泊水位監(jiān)測研究進展[J].遙感技術(shù)與應(yīng)用，2017，32（1）：14-19.

[9]?LI L，LI J，YAO X J，et al. Changes of the Three Holy Lakes in Recent Years and Quantitative Analysis of the Influencing Factors[J]. Quaternary International，2014，349：339-345.

[10]?張國慶.青藏高原湖泊變化遙感監(jiān)測及其對氣候變化的響應(yīng)研究進展[J].地理科學進展，2018，37（2）：214-223.

[11]?FANG Leqi，TAO Shengli，ZHU Jiangling，et al. Impacts of Climate Change and Irrigation on Lakes in Arid Northwest China[J]. Journal of Arid Environments，2018，154：34-39.

[12]?肖茜，楊昆，洪亮.近30 a云貴高原湖泊表面水體面積變化遙感監(jiān)測與時空分析[J].湖泊科學，2018，30（4）：1083-1096.

[13]?李寧，劉吉平，王宗明.2000—2010年東北地區(qū)湖泊動態(tài)變化及驅(qū)動力分析[J].湖泊科學，2014，26（4）：545-551.

[14]?劉新，何隆華，周馳.長江中下游近30年來湖泊的水域面積變化研究[J].華東師范大學學報（自然科學版），2008（4）：124-129.

[15]?侯光才，趙振宏，王曉勇，等.黃河中游鄂爾多斯高原內(nèi)流區(qū)與閉流區(qū)的形成機理：基于水循環(huán)的分析[J].地質(zhì)通報，2008（8）：1107-1114.

[16]?張國平，趙琳娜，許鳳雯，等.基于流域結(jié)構(gòu)分析的中國流域劃分方案[J].北京師范大學學報（自然科學版），2010，46（3）：417-423.

[17]?饒恩明，肖燚，歐陽志云，等.中國湖泊水量調(diào)節(jié)能力及其動態(tài)變化[J].生態(tài)學報，2014，34（21）：6225-6231.

[18]?李曉鋒，姚曉軍，孫美平，等.2000—2014年我國西北地區(qū)湖泊面積的時空變化[J].生態(tài)學報，2018，38（1）：96-104.

[19]?國家地球系統(tǒng)科學數(shù)據(jù)中心.國家地球系統(tǒng)科學數(shù)據(jù)共享服務(wù)平臺湖泊-流域科學數(shù)據(jù)庫[EB/OL].（2017-11-22）[2019-06-18].http：//lake.geodata.cn.

[20]?STOCKER T F， QIN D， PLATTNER G K，et al. IPCC Climate Change 2013： the Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[J]. Computational Geometry，2013，18（2）：95-123.

[21]?施雅風，沈永平，李棟梁，等.中國西北氣候由暖干向暖濕轉(zhuǎn)型的特征和趨勢探討[J].第四紀研究，2003（2）：152-164.

[22]?施雅風，劉時銀，上官冬輝，等.近30 a青藏高原氣候與冰川變化中的兩種特殊現(xiàn)象[J].氣候變化研究進展，2006（4）：154-160.

[23]?張寧瑾，肖天貴，假拉.1979—2016年青藏高原降水時空特征[J].干旱氣象，2018，36（3）：373-382.

[24]?賀添，邵全琴.基于MOD16產(chǎn)品的我國2001—2010年蒸散發(fā)時空格局變化分析[J].地球信息科學學報，2014，16（6）：979-988.

[25]?楊珂含，姚方方，董迪，等.青藏高原湖泊面積動態(tài)監(jiān)測[J].地球信息科學學報，2017，19（7）：972-982.

[26]?ZHANG Guoqing，YAO Tandong，CHEN Wenfeng，et al. Regional Differences of Lake Evolution Across China During 1960-2015 and Its Natural and Anthropogenic Causes[J]. Remote Sensing of Environment，2019，221：386-404.

[27]?李燕，金永明，孫永壽.柴達木盆地近55 a降水量變化趨勢的時空特征[J].人民黃河，2013，35（9）：48-50.

[28]?徐浩杰，楊太保.1981—2010年柴達木盆地氣候要素變化特征及湖泊和植被響應(yīng)[J].地理科學進展，2013，32（6）：868-879.

[29]?HU Qi，PAN Feifei，PAN Xuebiao，et al. Spatial Analysis of Climate Change in Inner Mongolia During 1961-2012，China[J]. Applied Geography，2015，60：254-260.

[30]?李虹雨，馬龍，劉廷璽，等.1951—2014年內(nèi)蒙古地區(qū)氣溫、降水變化及其關(guān)系[J].冰川凍土，2017，39（5）：1098-1112.

[31]?劉海月，沈渭壽，李海東，等.1957—2013年內(nèi)蒙古達里諾爾自然保護區(qū)氣候變化特征[J].云南大學學報（自然科學版），2016，38（4）：586-594.

[32]?佟斯琴，劉桂香，武娜.1961—2010年錫林郭勒盟氣溫和降水時空變化特征[J].水土保持通報，2016，36（5）：340-345，351.

[33]?康麗娟，巴特爾·巴克，羅那那，等.1961—2013年新疆氣溫和降水的時空變化特征分析[J].新疆農(nóng)業(yè)科學，2018，55（1）：123-133.

[34]?YAO Junqiang，HU Wenfeng，CHEN Yaning，et al. Hydro-Climatic Changes and Their Impacts on Vegetation in Xinjiang，Central Asia[J]. Science of the Total Environment，2019，660：724-732.

[35]?旺堆杰布，丹增尼瑪，白馬多吉，等.近40年西藏高原北部4個內(nèi)陸湖泊面積變化及氣候要素分析[J].高原科學研究，2018，2（2）：65-74.

[36]?紀鵬，郭華東，張露.近20年西昆侖地區(qū)冰川動態(tài)變化遙感研究[J].國土資源遙感，2013，25（1）：93-98.

[37]?于升松，劉興起，譚紅兵，等.茶卡鹽湖水文、水化學及資源開發(fā)研究[J].鹽湖研究，2005（3）：10-16.

[38]?朱長明，李均力，張新，等.近40 a來博斯騰湖水資源遙感動態(tài)監(jiān)測與特征分析[J].自然資源學報，2015，30（1）：106-114.

[39]?馬超，劉瑋瑋，趙鵬飛，等.1962—2016年達里諾爾流域水、草退化及氣候響應(yīng)[J].地理研究，2017，36（9）：1755-1772.

【責任編輯?張華興】