王君波，彭 萍，馬慶峰，朱立平

(中國科學(xué)院青藏高原研究所，青藏高原環(huán)境變化與地表過程重點(diǎn)實(shí)驗(yàn)室，北京100101)

瑪旁雍錯(cuò)(30°34'～30°47'N，81°22'～81°27'E)和拉昂錯(cuò)(30°40'～30°51'N，81°06'～81°19'E)流域是藏南內(nèi)流水系中為數(shù)不多的面積較大的內(nèi)陸湖泊［1］，由于瑪旁雍錯(cuò)在宗教傳統(tǒng)上的特殊地位，在過去一個(gè)多世紀(jì)以來就為世人所廣泛關(guān)注.廣為流傳的是瑞典人斯文·赫定在1907年前后對(duì)瑪旁雍錯(cuò)的調(diào)查，據(jù)記載其測定的最大水深達(dá)81.8 m;1976年中國科學(xué)院青藏高原綜合科學(xué)考察隊(duì)對(duì)瑪旁雍錯(cuò)進(jìn)行了較為詳細(xì)的調(diào)查，在水溫垂直分布、湖水化學(xué)性質(zhì)等方面獲得了寶貴的資料［2］，而作為位于同一流域甚至以前與瑪旁雍錯(cuò)連為一體的姊妹湖——拉昂錯(cuò)則相對(duì)研究的很少，目前湖泊基礎(chǔ)資料幾乎仍為空白.

全球氣候變化特別是全球變暖對(duì)青藏高原湖泊造成了顯著的影響，由于大部分湖泊幾乎不受人類活動(dòng)的直接影響，因而成為研究全球變化影響及區(qū)域響應(yīng)的重要區(qū)域.近年來，研究者選擇位于青藏高原不同氣候區(qū)的不同類型湖泊進(jìn)行了大量的考察工作，在湖泊水深分布、水質(zhì)特征、水化學(xué)組成及影響因素等方面取得了很多研究成果［3-8］，這些考察和研究一方面豐富了青藏高原湖泊的基礎(chǔ)資料，另一方面為湖泊環(huán)境演變研究提供了基礎(chǔ).

在全球變暖導(dǎo)致冰川退縮、湖泊面積變化的背景下，瑪旁雍錯(cuò)和拉昂錯(cuò)流域因其內(nèi)部具有完整的現(xiàn)代冰川、河流、封閉湖泊等水文要素而成為良好的研究地點(diǎn).基于遙感和GIS的研究顯示，該流域近三十余年來冰川面積減少、湖泊退縮［9-10］;在近年來大部分湖泊面積擴(kuò)張的情況下，瑪旁雍錯(cuò)面積變化不大，甚至略有萎縮趨勢［11-12］;由于補(bǔ)給和蒸發(fā)的差值不同，瑪旁雍錯(cuò)和拉昂錯(cuò)湖內(nèi)鹽分有較大的差異，因而造成二者水色的差異［13］.以上研究都是基于遙感影像而獲取的，對(duì)湖泊本身的特征尚未涉及.基于瑪旁雍錯(cuò)和拉昂錯(cuò)在冰川-湖泊相互作用、水文過程、水循環(huán)、水文模型以及環(huán)境變化等方面較好的研究潛力，筆者于2009年9月和2010年7月分別對(duì)瑪旁雍錯(cuò)和拉昂錯(cuò)進(jìn)行了調(diào)查，本文簡要報(bào)道這2個(gè)湖泊的水深、水質(zhì)特征及現(xiàn)代沉積速率.

1 方法

瑪旁雍錯(cuò)湖泊測深利用德國Innomar Technologie GmbH生產(chǎn)的沉積物回聲測量儀(Sediment Echo Sounder SES-96)完成，拉昂錯(cuò)測深利用廣州中海達(dá)測繪儀器有限公司生產(chǎn)的HD-27數(shù)字單頻測深儀(外接同步HD8500型GPS接收器)進(jìn)行.湖泊水量計(jì)算根據(jù)水深數(shù)據(jù)和經(jīng)緯度坐標(biāo)點(diǎn)利用ARCMAP和Surfer軟件完成.

湖水理化性質(zhì)測量利用美國哈希公司生產(chǎn)的Hydrolab DS5型多參數(shù)水質(zhì)監(jiān)測儀，該儀器可同時(shí)測量溫度、pH、電導(dǎo)率、溶解氧、光合有效輻射(PAR)以及水深等參數(shù).在瑪旁雍錯(cuò)共測量了13個(gè)水質(zhì)剖面，剖面的深度為0.3～58.1 m;在拉昂錯(cuò)不同湖區(qū)的6個(gè)點(diǎn)位進(jìn)行了水質(zhì)剖面測量，測點(diǎn)水深為10.3～43.0 m.此外，用2.5 L有機(jī)玻璃采水樣器在2個(gè)湖泊采集湖水樣，在離子色譜儀IC(Dionex-ICS900)下檢測主要陰陽離子的濃度，其中濃度根據(jù)離子平衡原理利用其他離子的濃度估算而得出.

利用重力采樣器在2個(gè)湖泊中采取了短巖芯，其中瑪旁雍錯(cuò)短巖芯(MPG09-1)長約43 cm，采樣點(diǎn)水深約為55 m，位于開闊湖區(qū)(圖1a);拉昂錯(cuò)短巖芯(LAG10-1)長約43.5 cm，采樣點(diǎn)水深約為45 m，位于主要湖區(qū)的中心開闊區(qū)(圖1b).短巖芯在野外按照0.5 cm間隔分樣帶回實(shí)驗(yàn)室，利用210Pb和137Cs方法測定沉積速率，具體測量儀器、方法及沉積速率計(jì)算方法參見文獻(xiàn)［14］.

2 結(jié)果

2.1 水深分布

瑪旁雍錯(cuò)湖盆形狀較為規(guī)則，從水下地形分布來看，也顯示了與流域湖盆較為相似的特征，四周湖岸坡度都較為平緩，中間深水區(qū)坡度更小，超過50 m水深的區(qū)域面積較大，實(shí)測最大水深72.60 m，位于中部深水開闊湖區(qū)的偏南部位置(圖1a);在北部淺水區(qū)，測深過程中發(fā)現(xiàn)水下地形起伏較大，可能存在湖泊低水位時(shí)期形成的沙丘.

拉昂錯(cuò)湖體形狀非常不規(guī)則，略呈湯勺狀，且湖中出露多個(gè)島嶼，北部湖區(qū)為深度10 m左右面積較小的開闊湖區(qū)，為河流沖積扇;中間部分是狹長的連接處，水深也不超過30 m;南部為較大的開闊湖區(qū)，中心部分平坦，水深超過45 m，最大實(shí)測水深為49.03 m，位于中間開闊湖區(qū)的偏南位置(圖1b);西部湖灣處未進(jìn)行實(shí)測，從地形、河流補(bǔ)給和島嶼分布情況推測，應(yīng)該為較淺的水下沖積扇區(qū)域，位于深水區(qū)的小島四周都非常陡峭.

2.2 湖水理化性質(zhì)

MC09-7(圖2a)是瑪旁雍錯(cuò)野外實(shí)測樣點(diǎn)中最深的水質(zhì)剖面(深度58 m，圖1a)，湖水溫度分層明顯，表層水溫在0～25 m之間穩(wěn)定在10.4℃左右，溫躍層介于25～35 m之間，溫度梯度為0.17℃/m，在35 m以下水溫變化趨緩，底層水溫穩(wěn)定在7.2℃左右.受溫度梯度變化的影響，pH值逐漸遞減，從表層至底層的湖水pH介于9.6～7.2之間，呈堿性至弱堿性的變化趨勢.瑪旁雍錯(cuò)是淡水湖，湖水鹽度為0.27 ppt，湖水電導(dǎo)率值介于532～546 μS/cm之間.水質(zhì)剖面MC09-7電導(dǎo)率從表層至底層呈現(xiàn)明顯的降低趨勢，尤其是介于25～35 m的溫躍層間.溶解氧在表層至30 m水深之間變化不大，約為10.6 mg/L，但是隨著深度進(jìn)一步增加，溶解氧呈現(xiàn)顯著階梯型降低趨勢，底層水的溶解氧為8.5 mg/L.光合有效輻射是表征湖水透光程度的參數(shù)，在剖面MC09-7中，光合有效輻射在表層至30 m水深之間，由1633 μmol/(m2·s)急劇降至30 μmol/(m2·s)，在50 m水深以下維持在1 μmol/(m2·s)，30 m水層內(nèi)的衰減系數(shù)為0.17 m－1.

拉昂錯(cuò)水質(zhì)剖面LAH 10-5(圖2b)位于開闊湖區(qū)偏北位置(深度40.5 m，圖1b)，在5～10 m水深之間水溫變化明顯，溫度梯度為0.16℃/m，在10～30 m之間水溫緩慢降低到8.3℃并維持至底層.pH在表層至10 m水深之間變化顯著，介于8.9～6.0之間，由表層的弱堿性變?yōu)槿跛嵝?隨著深度的增加，pH值有緩慢回升趨勢，但仍然維持在7.0以下，底層湖水環(huán)境仍為弱酸條件.湖水剖面電導(dǎo)率變化介于1452～1460 μS/cm之間，在水溫迅速降低的5～10 m水深間有顯著降低，隨后又緩慢升高.溶解氧從表層開始有所升高，但總體變化不大，介于10.9～11.1 mg/L之間，30 m水深以下呈明顯降低趨勢.光合有效輻射在表層急劇降低，但從整個(gè)剖面來看透光性較好，30 m水層內(nèi)的衰減系數(shù)為0.15 m－1.

對(duì)瑪旁雍錯(cuò)湖水剖面MCW 09檢測7個(gè)10 m間隔的水樣離子濃度，結(jié)果顯示主要陽離子中平均濃度大小順序?yàn)?Na+＞Mg2+＞Ca2+＞ K+，其平均濃度分別為 49.09、28.73、26.86 和 5.97 mg/L;而主要陰離子順序?yàn)椋維O24－＞Cl－＞F－，其平均濃度分別為 300.71、29.37、13.83 和1.17 mg/L(表1)，可見和Na+是占優(yōu)勢的離子，據(jù)此計(jì)算湖水平均礦化度為455.7 mg/L.

在拉昂錯(cuò)則檢測3個(gè)水樣點(diǎn)(LAW 10-1、LAW 10-2、LAW 10-3)的表層和底層水樣離子濃度，結(jié)果顯示主要陽離子中平均濃度大小順序?yàn)?Na+＞Mg2+＞K+＞Ca2+，其平均濃度分別為176.02、107.49、20.20和12.95 mg/L;而主要陰離子平均濃度大小順序?yàn)椋荆?Cl－＞F－，其平均濃度分別為 813.61、107.20、68.78 和2.98 mg/L(表1)，可見拉昂錯(cuò)湖水中和Na+也是占優(yōu)勢的離子，據(jù)此計(jì)算湖水平均礦化度為 1309.3 mg/L.

表1 瑪旁雍錯(cuò)和拉昂錯(cuò)湖水主要離子組成特征(mg/L)Tab.1 Major ions composition of lake water from Mapam Yumco and La'ang Co

2.3 現(xiàn)代沉積速率

瑪旁雍錯(cuò)和拉昂錯(cuò)短巖芯的210Pb和137Cs的測試結(jié)果如圖3、圖4所示，基于過剩210Pb活度分別用CRS(恒定補(bǔ)給速率)和CIC(常量初始濃度)方法計(jì)算沉積速率［14］，并與137Cs的峰值進(jìn)行了對(duì)比，結(jié)果顯示利用CIC方法得出的平均沉積速率與137Cs的時(shí)標(biāo)較為接近，且從計(jì)算結(jié)果來看，其擬合曲線的方差解釋量(R2)分別為0.7904和0.7558，說明擬合效果較好，也指示這2個(gè)湖泊的近代沉積速率都較為穩(wěn)定，因而本文采用此方法的結(jié)果.

圖3 瑪旁雍錯(cuò)PMG 09-1巖芯210Pb和137Cs隨深度的變化及年代-深度關(guān)系Fig.3210Pb and137Cs activities and age depth model of gravity core PMG 09-1 from Mapam Yumco

瑪旁雍錯(cuò)MPG09-1巖芯在6 cm深度處(圖中取樣品中間位置深度為5.75 cm，下同)過剩210Pb活度接近于零，據(jù)此利用CIC方法計(jì)算出的平均沉積速率(ASR)為0.31 mm/a(圖3A);137Cs活度非常低，3 cm深度處開始檢測到137Cs的活度，峰值出現(xiàn)在1.5 cm處，活度值為11.73 Bq/kg(圖3B);利用平均沉積速率獲得的年代-深度關(guān)系與137Cs峰值處對(duì)應(yīng)于1963年進(jìn)行對(duì)比，可以發(fā)現(xiàn)二者吻合很好(圖3C)，說明本研究中獲得的平均沉積速率具有較好的可信度.

拉昂錯(cuò)LAG 10-1巖芯的過剩210Pb和137Cs的活度值與瑪旁雍錯(cuò)的較為相似，說明這2個(gè)湖區(qū)接受的核素沉降值非常接近.LAG10-1巖芯在9 cm深度處過剩210Pb活度接近于零，利用CIC方法計(jì)算得出平均沉積速率(ASR)為0.65 mm/a(圖4A);137Cs在5 cm深度處開始出現(xiàn)，峰值為3.5 cm，活度值為11.67 Bq/kg(圖4B)，2.5和1.5 cm深度處各出現(xiàn)一個(gè)較高值，反映了該湖表層沉積物可能擾動(dòng)較為劇烈的特征，此處取3.5 cm處的峰值作為1963年的年齡時(shí)標(biāo);圖4C顯示了利用CIC 方法計(jì)算出的年代-深度關(guān)系與1963年的137Cs峰值對(duì)應(yīng)很好.由此可見，拉昂錯(cuò)的近代平均沉積速率約是瑪旁雍錯(cuò)的2倍.

3 討論與結(jié)論

本次對(duì)瑪旁雍錯(cuò)的水深測定較為詳細(xì)地繪制了該湖的水下地形分布，并確定瑪旁雍錯(cuò)目前的最大水深為72.6 m，由于調(diào)查測線基本上覆蓋了所有湖區(qū)，特別是中間開闊湖區(qū)測線分布較為密集，因而可以確定水深數(shù)據(jù)不會(huì)有大的遺漏.關(guān)于瑞典人斯文·赫定1907年對(duì)瑪旁雍錯(cuò)的水深測定，目前有兩種版本流傳，一種是其最大實(shí)測水深是81.8 m，根據(jù)1976年中國科學(xué)院組織的科學(xué)考察而出版的《西藏河流與湖泊》中關(guān)于瑪旁雍錯(cuò)沒有詳細(xì)的水深數(shù)據(jù)，僅提到了斯文·赫定實(shí)測最大水深是81.8 m，并據(jù)此估算了瑪旁雍錯(cuò)的貯水量達(dá)200×108m3，使其成為地球上高海拔地區(qū)淡水最多的湖泊之一［2］.之后1998年出版的較為權(quán)威的《中國湖泊志》中瑪旁雍錯(cuò)的資料也使用了這一數(shù)據(jù)及結(jié)論［1］.2008年姚檀棟等［15］的《青藏高原及毗鄰地區(qū)冰川湖泊圖》上也將瑪旁雍錯(cuò)的最大水深標(biāo)注為81.8 m，顯然這些說法都來自于相同的文獻(xiàn)，經(jīng)查應(yīng)該是1917年在瑞典出版的基于斯文·赫定第三次探險(xiǎn)考察活動(dòng)的《Southern Tibet》，然而筆者未能查證此數(shù)據(jù)的可靠性.

圖4 拉昂錯(cuò)LAG 10-1巖芯210Pb和137Cs隨深度的變化及年代-深度關(guān)系Fig.4210Pb and137Cs activities and age depth model of gravity core LAG 10-1 from La'ang Co

第二種版本則可見于中文出版物中，新疆人民出版社亞洲探險(xiǎn)之旅叢書中《失蹤雪域750天》是斯文·赫定考察日記式的著作，該書詳細(xì)記錄了斯文·赫定在1907年7月27日夜晚和8月7日白天對(duì)瑪旁雍錯(cuò)的兩次考察，期間測得的最大水深是77 m［16］.雖無法確定斯文·赫定當(dāng)時(shí)測量時(shí)的精確路線及最大水深的位置，但根據(jù)其描述可判斷其測線主要在南部湖區(qū)，即瑪旁雍錯(cuò)的深水區(qū)，且繩測的結(jié)果應(yīng)該較為準(zhǔn)確，因而基本上可以把斯文·赫定測得的77 m作為瑪旁雍錯(cuò)當(dāng)時(shí)的最大水深.

根據(jù)以上數(shù)據(jù)，則可推斷與1907年相比，現(xiàn)在的瑪旁雍錯(cuò)湖面至少下降了4.4 m，若取第一種說法的81.8 m作為當(dāng)時(shí)的最大水深，則湖面至少下降了9.2 m.瑪旁雍錯(cuò)現(xiàn)在湖水水量約為146×108m3，也比當(dāng)時(shí)估算的約200×108m3減少了很多.瑪旁雍錯(cuò)主要依靠湖面降水補(bǔ)給，現(xiàn)在湖水變淺說明與1907年相比湖泊補(bǔ)給減少了，反推20世紀(jì)初期瑪旁雍錯(cuò)流域降水量較大.青海湖近百年來(1880s-1980s)水位下降了約12 m，20世紀(jì)初期同樣維持了高湖面直到1930s湖面開始迅速下降［17］;同樣位于藏南地區(qū)的羊卓雍錯(cuò)近百年(1912-1992年)水位變化在4～5 m之間，且水位呈平緩下降趨勢［18］.說明這些湖泊近百年來水位變化具有相似的特征，而20世紀(jì)初期是一個(gè)降水較多的濕潤時(shí)期.

瑪旁雍錯(cuò)和拉昂錯(cuò)同屬藏南湖區(qū)，地質(zhì)及氣候條件相似［2］，但是相比之下，瑪旁雍錯(cuò)周邊河流比拉昂錯(cuò)多，接納了更多來自北部及東部低地的冰川融水，在相同的氣候條件下湖水鹽度比拉昂錯(cuò)更低［13］，因而2個(gè)湖泊在湖水電導(dǎo)率上存在著近3倍差異.從2個(gè)湖泊湖水的主要離子組成來看，和Na+都是占主要優(yōu)勢的離子，且2個(gè)湖泊主要陰離子濃度大小順序也一致，主要陽離子中除Ca2+的濃度外，其余離子濃度也都顯示了相似的濃度大小順序.拉昂錯(cuò)Ca2+的濃度約是瑪旁雍錯(cuò)的1/2，而其余離子濃度平均為瑪旁雍錯(cuò)離子濃度的3倍左右(表1).拉昂錯(cuò)湖水中較低的Ca2+濃度與湖水中CaCO3的沉淀有關(guān)，這與納木錯(cuò)湖水具有相似性［19］.

2個(gè)湖泊水質(zhì)剖面測試于不同年份幾乎相同的時(shí)間段內(nèi)，剖面表層水與底層水的溫差相近，但是剖面水溫變化趨勢卻存在很大差異.瑪旁雍錯(cuò)湖水呈現(xiàn)穩(wěn)定的正溫層分布，與青藏高原其他湖泊相比有著相似的熱學(xué)特征［8］，而拉昂錯(cuò)湖水溫度從表層至底層呈現(xiàn)急劇遞減趨勢.考慮到瑪旁雍錯(cuò)補(bǔ)給來源遠(yuǎn)多于拉昂錯(cuò)，冷的冰川融水大量匯入湖中，上、下水層密度存在差異引起湖水混合運(yùn)動(dòng)，而且瑪旁雍錯(cuò)較深的湖盆環(huán)境可能更有利于湖水的均勻混合.但是，拉昂錯(cuò)的補(bǔ)給主要來自于北部冰川融水補(bǔ)給，南部湖盆也相對(duì)較寬淺，湖水溫度變化相對(duì)單調(diào)，未達(dá)到穩(wěn)定分層，水溫呈現(xiàn)緩慢降低趨勢.因此，推測2個(gè)湖泊剖面水溫的差異可能是由于補(bǔ)給不同疊加湖泊深度差異造成的.瑪旁雍錯(cuò)湖水剖面酸堿度相對(duì)穩(wěn)定，從表層至底層保持堿性-弱堿性變化，說明湖水混合相對(duì)均勻.但是拉昂錯(cuò)10 m深度以下呈現(xiàn)出明顯的弱酸性變化，對(duì)比溶解氧變化來看，pH與溶解氧呈現(xiàn)反相變化趨勢，可能溶解氧的豐富程度與水體酸堿度間存在著某種聯(lián)系.而野外采樣中發(fā)現(xiàn)拉昂錯(cuò)表層沉積物呈罕見的泥黑色且夾雜大量植物碎屑，推測該湖水生植物相對(duì)茂盛，生物呼吸作用占優(yōu)勢，釋放更多CO2到水體中，導(dǎo)致下層水體pH呈弱酸性.對(duì)比湖水光學(xué)性質(zhì)來看，瑪旁雍錯(cuò)和拉昂錯(cuò)湖水30 m 水層內(nèi)的光合有效輻射衰減系數(shù)分別為 0.17 m－1和 0.15 m－1，與納木錯(cuò)(0.12 m－1)［3-4］和普莫雍錯(cuò)(0.15 m－1)［5］、扎日南木錯(cuò)(0.14 m－1)和當(dāng)惹雍錯(cuò)(0.16 m－1)［8］有很好的可比性.

與前人在1976年7月份進(jìn)行的測量相比［2］，瑪旁雍錯(cuò)的水溫分布趨勢非常一致，其溫躍層深度及湖底溫度非常接近，說明30余年來該湖水體的溫度結(jié)構(gòu)變化不大，只是前人測量的湖水表面最大水溫約為12℃，應(yīng)該與測量季節(jié)的不同有關(guān).1976年考察時(shí)湖水礦化度約為405.8 mg/L［2］，而本次考察根據(jù)湖水主要離子含量計(jì)算的礦化度約為455.7 mg/L，反映瑪旁雍錯(cuò)湖水正在變咸，水面下降和水量的減少都是湖水變咸的直接原因.1976年不同測點(diǎn)的湖水表層pH為8.0～8.4，而現(xiàn)在該湖上層水的pH已高達(dá)9.6，若排除不同儀器所帶來的差異，則瑪旁雍錯(cuò)湖水的pH值在過去30多年來具有明顯的增加，也反映了湖水的咸化過程.瑪旁雍錯(cuò)湖水pH值和礦化度的升高與藏南羊卓雍錯(cuò)流域若干個(gè)湖泊的情況一致，即該流域內(nèi)湖水的pH和礦化度從1979年以來總體上也呈增加趨勢［20］.

已有的研究顯示，青藏高原上的湖泊現(xiàn)代沉積速率大多比較低，基本上都是每年毫米的沉積量，甚至低于此數(shù)值.如納木錯(cuò)多支巖芯指示的現(xiàn)代沉積速率為0.43～0.98 mm/a，與青藏高原第一大湖——青海湖較為接近［13］.由于較低的沉積速率，以及核素沉降后可能發(fā)生的沉積后過程特別是137Cs的垂直遷移，致使在一些湖泊中利用210Pb和137Cs方法進(jìn)行現(xiàn)代沉積速率的研究較為困難，尤其是在一些湖泊中137Cs并無峰值出現(xiàn)，因而結(jié)果準(zhǔn)確性較差.本研究中的瑪旁雍錯(cuò)和拉昂錯(cuò)近代平均沉積速率分別為0.31 mm/a和0.65 mm/a，且利用210Pb計(jì)算的平均沉積速率得到了137Cs年齡時(shí)標(biāo)較好的驗(yàn)證，因而具有較好的可信度.

作為藏南內(nèi)流區(qū)2個(gè)重要的湖泊，除了瑪旁雍錯(cuò)仍將以其第一神湖的聲譽(yù)繼續(xù)受到世人矚目外，它也會(huì)和拉昂錯(cuò)一起在全球變化、水循環(huán)和人類活動(dòng)的影響與適應(yīng)等重要科學(xué)研究領(lǐng)域受到越來越多的關(guān)注.本文基于初步考察所報(bào)道的瑪旁雍錯(cuò)和拉昂錯(cuò)的水深、水質(zhì)特征以及現(xiàn)代沉積速率將會(huì)為該地區(qū)的水量平衡、水文模型、湖泊流域管理和過去環(huán)境變化研究等提供基礎(chǔ)資料和重要參考.

致謝:德國耶拿大學(xué)地理系Dr.Gerhard D參與了部分野外工作，中國科學(xué)院青藏高原研究所高少鵬高級(jí)工程師完成了210Pb和137Cs測試，楊瑞敏在湖泊水量計(jì)算中給予了幫助，野外考察過程中得到了當(dāng)?shù)卣嚓P(guān)部門的協(xié)助和支持，在此一并表示感謝.

［1］王蘇民，竇鴻身.中國湖泊志.北京:科學(xué)出版社，1998:1-580.

［2］關(guān)志華，陳傳友，區(qū)裕雄等.西藏河流與湖泊.北京:科學(xué)出版社，1984:1-238.

［3］王君波，朱立平，Gerhard D等.西藏納木錯(cuò)水深分布及現(xiàn)代湖沼學(xué)特征初步分析.湖泊科學(xué)，2009，21(1):129-135.

［4］Wang JB，Zhu LP，Daut G et al.Investigation of bathymetry and water quality of Lake Nam Co，the largest lake on the central Tibetan Plateau，China.Limnology，2009，10:149-158.

［5］Mitamura O，Seike Y，Kondo K et al.First investigation of ultraoligotrophic alpine Lake Puma Yumco in the pre-Himalaya，China.Limnology，2003，4:167-175.

［6］Murakami T，Terai H，Yoshiyama Y et al.The second investigation of Lake Puma Yum Co located in the Southern Tibetan Plateau，China.Limnology，2007，8:331-335.

［7］Zhu LP，Ju JT，Wang JB et al.Further discussion about the features of Lake Puma Yum Co，South Tibet，China.Limnology，2010，11:281-287.

［8］王君波，彭 萍，馬慶峰等.西藏當(dāng)惹雍錯(cuò)和扎日南木錯(cuò)現(xiàn)代湖泊基本特征.湖泊科學(xué)，2010，22(4):630-634.

［9］郭柳平，葉慶華，姚檀棟等.基于GIS的瑪旁雍錯(cuò)流域冰川地貌及現(xiàn)代冰川湖泊變化研究.冰川凍土，2007，29(4):517-524.

［10］葉慶華，姚檀棟，鄭紅星等.西藏瑪旁雍錯(cuò)流域冰川與湖泊變化及其對(duì)氣候變化的響應(yīng).地理研究，2008，27(5):1178-1190.

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

［12］拉 巴，邊 多，次 珍等.西藏瑪旁雍錯(cuò)流域湖泊面積變化及成因分析.干旱區(qū)研究，2012，29(6):992-996.

［13］萬 瑋，馮學(xué)智，肖鵬峰等.青藏高原瑪旁雍錯(cuò)和拉昂錯(cuò)水色差異遙感分析.遙感技術(shù)與應(yīng)用，2008，23(6):667-671.

［14］王君波，朱立平，汪 勇等.西藏納木錯(cuò)現(xiàn)代沉積速率的空間分布特征及近60年來的變化研究.第四紀(jì)研究，2011，31(3):535-543.

［15］姚檀棟，任賈文，徐柏青.青藏高原及毗鄰地區(qū)冰川湖泊圖:中文版.西安:西安地圖出版社，2008.

［16］斯文·赫定著.包菁萍譯.失蹤雪域750天.烏魯木齊:新疆人民出版社，2000:4.

［17］馮 松，湯懋蒼，周陸生.青海湖近600年的水位變化.湖泊科學(xué)，2000，12(3):205-210.

［18］劉天仇.西藏羊卓雍錯(cuò)水位動(dòng)態(tài)研究.地理科學(xué)，1995，15(1):55-62.

［19］Wang JB，Zhu LP，Wang Y et al.Comparisons between the chemical compositions of lake water，inflowing river water，and the lake sediment in Nam Co，central Tibetan Plateau，China and their controlling mechanisms.Journal of Great Lakes Research，2010，36:587-595.

［20］張雪芹，孫 瑞，朱立平.藏南羊卓雍錯(cuò)流域主要湖泊水質(zhì)狀況及其評(píng)價(jià).冰川凍土，2012，34(4):950-958.