基于GF衛(wèi)星解譯巴丹吉林沙漠湖泊水量變化

曹樂 聶振龍 姜高磊

摘 要：遙感技術(shù)作為科學(xué)、快速、大面積調(diào)查監(jiān)測手段，在湖泊演化、動態(tài)變化研究中應(yīng)用廣泛。通過解譯國產(chǎn)GF-1、GF-2衛(wèi)星遙感影像信息，結(jié)合實地勘測，解譯了巴丹吉林沙漠110個湖泊地質(zhì)歷史時期與現(xiàn)代的湖面高程和邊界，總結(jié)了湖泊萎縮規(guī)律，估算了兩期水量變化。結(jié)果顯示：湖面高程平均降低9.76 m，面積共減少61.052 km2（占古湖總面積的75.49%），湖泊群水量共減少4.9億m3，說明湖泊萎縮程度高，沙漠干旱化趨勢明顯;研究區(qū)古、今地下水等水位線均表現(xiàn)出東南高、西北低的宏觀特征，反映了古、今沙漠區(qū)地下水相同的補給來源與徑流條件;基底凹陷區(qū)湖泊萎縮程度較低（如蘇木吉林湖），基底隆起區(qū)湖泊萎縮程度較高（如雅布賴山前湖泊），萎縮程度的差異性受區(qū)域構(gòu)造基底的控制與影響。

關(guān)鍵詞：遙感解譯;古湖高程;湖泊萎縮;水量變化;巴丹吉林沙漠

中圖分類號：TV213.3 ? 文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2020.07.009

Abstract：Remote sensing technology is widely used in the study of lake evolution and dynamic changes as a scientific， rapid and large-scale investigation and monitoring method. This study obtained the elevations and boundaries of 110 lakes in ancient and modern in the Badain Jaran Desert， based on the GF-1 and GF-2 satellite remote sensing images and regional survey. The results show that the average lakes elevation is decreased by 9.76 m and the area of lakes reduced totally 61.052 km2， which accounts for 75.49% of the ancient lakes area. The total water volume of lakes has been reduced 490 million m3， showing high degree of shrinkage of lakes and obvious drought trend. The counter lines of groundwater have the characters of higher in southeast and lower in northwest generally whether in ancient or at present and those demonstrate the source of groundwater recharge and the groundwater runoff condition remains unchanged. The lakes shrinking are obviously in the basement uplift areas （e.g. Yabulai piedmont）， while lightly in the basement depression areas （e.g. Lake Sumu Jaran）. The differences of lakes shrinkage are controlled by the regional tectonic basement.

Key words： remote sensing interpretation; ancient lakes elevation; lakes shrinking; status change; Badain Jaran Desert

巴丹吉林沙漠位于阿拉善高原西部，面積4.9萬km2，沙漠中有高大沙山以及100多個湖泊，為世界罕見[1-3]。關(guān)于湖泊成因、演變、補給來源等問題一直是第四紀(jì)地質(zhì)與水文地質(zhì)界研究的熱點，但尚無定論[4-13]。遙感技術(shù)能夠全方位、大區(qū)域地進(jìn)行湖泊季節(jié)性、年際動態(tài)變化監(jiān)測，是湖泊研究的重要技術(shù)手段[14]。朱金峰等[15]利用ETM+遙感影像，提取了湖泊年內(nèi)季節(jié)變化信息，認(rèn)為湖泊面積和數(shù)量在年內(nèi)隨季節(jié)依次減少，到第二年春季又恢復(fù)到前一年春季狀態(tài)。張振瑜等[16-17]利用Landsat影像分析了1973—2010年湖泊的時空變化特征，認(rèn)為地下水補給源水量變化和補給方式的不同控制著湖泊面積的時空變化。Zhang等[18]利用Landsat遙感影像分析了1967—2010年50個湖泊的面積變化特征，認(rèn)為目前當(dāng)?shù)亟邓畬囱a給有限，地下水補給源來自于古湖殘余水體。Jiao等[19]利用衛(wèi)星數(shù)據(jù)與重力恢復(fù)試驗證明了沙漠湖泊水位與地下水水量在緩慢減少。

根據(jù)湖水化學(xué)特征推測，該地區(qū)湖泊演化了800～7 000 a[1， 20]，已有的高于現(xiàn)代湖水位的湖相沉積測年結(jié)果表明湖泊區(qū)在全新世中期處于泛湖期[21]。借助遙感技術(shù)，對巴丹吉林沙漠湖泊的季節(jié)性變化和近幾十年的年際變化特征已有了初步的認(rèn)識[15-18]，但對于地質(zhì)歷史時期湖泊面積、水量變化還未開展研究。筆者利用國產(chǎn)GF-1、GF-2衛(wèi)星PMS傳感器影像資料，結(jié)合Google Earth影像和實地調(diào)查對巴丹吉林沙漠110個湖泊進(jìn)行了演化解譯，確定了最大古湖的高程和邊界，以期為研究該區(qū)古地下水流場、湖泊萎縮特征等提供參考。

1 研究區(qū)概況

巴丹吉林沙漠位于雅布賴山與北大山以北、宗乃山以西、古日乃與弱水以東、拐子湖以南（見圖1）。

沙漠區(qū)整體地勢東南高、西北低，海拔900～1 600 m，地形起伏較大，東南部高大沙山居多，相對高差一般為200～300 m[5， 22]。高大沙山間分布有100多個湖泊，呈串珠狀排列，多為咸水湖，面積一般小于1 km2[21]。沙漠區(qū)降水稀少，多年平均降水量為30～120 mm，由東南向西北逐漸減少[23-24]，年平均水面蒸發(fā)量為1 450 mm[25]。沙漠區(qū)主要發(fā)育北東東向、北西向斷裂，控制著區(qū)內(nèi)基底結(jié)構(gòu)、隆起格架及中生代盆地的展布，使盆地形成凹凸相間的結(jié)構(gòu)特征，由南向北主要有雅布賴山隆起、蘇亥圖凹陷、宗乃山隆起、陶勒特凹陷、特羅西灘低隆起等構(gòu)造單元[26-27]。

2 研究方法

2.1 遙感解譯數(shù)據(jù)選取

湖泊在演化過程中處于各穩(wěn)定期次時，在地表殘留的湖相沉積物、各類植被及其微地貌形態(tài)所共同構(gòu)成的規(guī)則閉合邊界只有在高分辨率遙感影像中才能準(zhǔn)確解譯。為保證解譯準(zhǔn)確性，將用于湖泊演化期次解譯的國產(chǎn)GF-1衛(wèi)星16 m分辨率、GF-1衛(wèi)星2 m分辨率（融合后）、GF-2衛(wèi)星1 m分辨率（融合后）數(shù)據(jù)以及Google Earth提取的0.6 m分辨率數(shù)據(jù)進(jìn)行比對。選取音德爾圖作為比對湖泊，以此確定湖泊演化期次提取所需最優(yōu)空間分辨率的區(qū)間范圍（見圖2）。結(jié)果顯示，演化期次邊界在亞米級影像上可以準(zhǔn)確識別，最優(yōu)空間分辨率為0.6 ～1.0 m （分辨率高于0.6 m對解譯結(jié)果準(zhǔn)確性的提升影響不大）;而分辨率低于5 m的影像，難以識別湖泊演化期次邊界。故GF-1衛(wèi)星16 m分辨率影像不能滿足要求，GF-2衛(wèi)星及Google Earth所獲取的遙感數(shù)據(jù)較為適用。

2.2 湖泊邊界的確定

結(jié)合遙感影像與實地調(diào)查結(jié)果，將湖泊演化期次邊界分為古湖邊界、中間穩(wěn)定期次邊界、今湖邊界3類，各演化期次邊界及區(qū)域分布關(guān)系如圖3所示，古湖邊界以影像上可見面積最大的湖岸線為準(zhǔn)，今湖邊界以2015年3月影像中湖泊水體面積的極大值為準(zhǔn)。

各類演化區(qū)域邊界鑒定特征見圖4：①網(wǎng)格狀沙丘區(qū)，由微型曲弧狀沙梁及下凹沙窩構(gòu)成，網(wǎng)格及蜂窩狀形態(tài)明顯，沙窩內(nèi)發(fā)育灌叢及草本植物，整體色調(diào)較淺，多為土黃色，與外側(cè)緩起伏沙地及沙山區(qū)別明顯，其邊緣指示古湖邊界;②草灌叢沙堆區(qū)、植被覆蓋區(qū)、鹽漬化區(qū)、潛水蒸發(fā)區(qū)影像上形態(tài)起伏明顯，表面平滑與粗糙相異，色彩深淺不同，兩區(qū)交界邊緣均指示某一湖泊穩(wěn)定期邊界;③季節(jié)性湖泊區(qū)，影像表面光滑，水體多呈不規(guī)則形態(tài)，水體間裸露少量土體，可見斑狀白色鹽堿，整體色調(diào)呈深藍(lán)色，水體季節(jié)性變化明顯，故其邊界年內(nèi)變化較大，指示今湖邊界;④常年性湖泊區(qū)，影像表面光滑，水體形態(tài)規(guī)則完整，色調(diào)為深藍(lán)色，水體面積季節(jié)性變化小，邊界年內(nèi)變化小，指示今湖邊界。

2.3 湖泊高程提取與驗證

對于研究古湖及今湖高程特征而言，其相對高程比絕對高程更為重要，為此從不同來源的數(shù)據(jù)獲取固定點對的高程差，判斷相對高程的誤差特征。比對野外移動GPS獲取的高程數(shù)據(jù)與室內(nèi)提取的DEM高程數(shù)據(jù)，采用RTK（Real-Time Kinematic）方法所采集的高精度高程值進(jìn)行驗證（控制點相對高程誤差控制在5 cm內(nèi)）。表1顯示，以RTK數(shù)據(jù)為基準(zhǔn)，DEM數(shù)據(jù)所得點對相對高程誤差小于2 m，無正負(fù)號偏差，在區(qū)域研究尺度上是可以接受的;而移動GPS的相對高程出現(xiàn)正負(fù)號偏差，且數(shù)值與DEM高程相差可達(dá)10 m以上。因此，本研究統(tǒng)一采用DEM高程數(shù)據(jù)。

古湖/今湖邊界的高程提取方法是利用已有的古湖/今湖矢量邊界，將其等分為多條（如100條）矢量線段后，獲取各矢量線段中點，分別提取各中點對應(yīng)的DEM高程數(shù)據(jù)并求其平均值，最終獲得的邊界高程均值為其高程值。高程提取主要步驟如圖5所示。

2.4 野外驗證

選擇代表性湖泊進(jìn)行演化期次邊界野外驗證，驗證古湖67個，占湖泊總數(shù)的61%;通過野外典型標(biāo)志的確定（見圖6），證明了解譯的準(zhǔn)確性。以昂斯格湖與巴丹湖區(qū)驗證為例，昂斯格湖影像（圖6（b））中箭頭①②為古湖邊界，實地調(diào)查時發(fā)現(xiàn)①處的古湖邊界處于沙山半坡上，明顯可見古湖邊界近水平出露于沙山表面（圖6（a））;②處指示湖泊東部古湖邊界，實地可見古湖沉積物出露地表，其高程與對岸古湖邊界一致（圖6（a））;③處指示網(wǎng)格狀沙丘區(qū)（圖6（a））。巴丹湖區(qū)東北部出露的古湖沉積物與遙感劃定古湖邊界位置相符（圖6（c）），據(jù)野外調(diào)查，巴丹湖及附近湖泊可能曾經(jīng)為同一個大湖，在后期演化萎縮過程中，古湖區(qū)被沙體掩埋覆蓋，最終形成現(xiàn)在的多個湖泊。

3 結(jié)果分析

3.1 湖泊高程與面積變化

圖7為所有古湖與今湖的高程與面積，可以看出，所有今湖湖面高程下降，平均降幅為9.76 m，面積萎縮。湖面高程下降0～10 m的有61個湖泊，下降10～20 m的有43個湖泊，下降20～25 m的有6個湖泊。從解譯的古湖面積來看，最大的為昂斯格，面積3.412 6 km2，最小的為葫蘆斯臺（東湖），面積僅0.007 2 km2;從解譯的今湖面積來看，最大的為東諾爾圖，面積1.417 5 km2，最小的為葫蘆斯臺（中湖），面積僅為0.000 4 km2。作為地表水體，地質(zhì)歷史時期最大古湖面積總計80.877 km2，今湖現(xiàn)存水體面積總計19.825 km2，減小61.052 km2，占古湖總面積的75.49%。表2顯示，110個湖泊的萎縮程度（（1-今湖面積/古湖面積）×100%）均大于25%，大于90%的湖泊有48個，大于95%的有33個。不同湖泊的萎縮程度各不相同，萎縮程度最大的湖泊是達(dá)布蘇圖巴潤敖格欽，今湖水體面積僅為古湖面積的0.12%;萎縮程度最小的湖泊是扎拉特，今湖水體面積占古湖面積的74.96%。

3.2 湖泊水量平衡

圖8為110個湖泊的水量平衡計算結(jié)果（湖泊編號按湖泊面積從大到小排列），可以看出，所有湖泊體積減小，面積越大的古湖水量減少越多。計算得到110個湖泊水量總共減少4.9億m3，其中：蘇木吉林與巴潤蘇木吉林作為一個古湖萎縮分裂成的兩個湖泊，減少的水量最大，為0.32億m3;未分裂的單個湖泊中，伊和吉格德減小量最大，為0.27億m3。

3.3 古、今湖泊等水位線特征

地表水位的分布特征一定程度上受地下水系統(tǒng)影響，因此可利用湖面高程對地下水系統(tǒng)特征進(jìn)行研究。圖9顯示，利用古、今湖邊界平均高程數(shù)據(jù)得出的等水位線分布特征基本一致，二者均顯示出地下水位東南高、西北低的宏觀特征;東南部地下水水力梯度較大，西北部相對較小，說明地下水流向未有明顯變化。在圖9中選取3點P1、P2、P3，分段計算水力梯度，結(jié)果見表3，顯示今湖的地下水水力梯度均大于古湖的，P1—P2段的水力梯度較大，為0.707%;從東南部水力梯度較大區(qū)到西北部，水力梯度逐漸減小，P2—P3段今湖的平均水力梯度為0.070%，P1—P3段今湖的地下水平均水力梯度為0.219%。

0.070%～0.707%的水力梯度變幅與張競等[13， 29]的研究結(jié)果0.08%～0.79%很相近，說明古湖與今湖地下水由東南向西北的主流向一致，反映地下水補給來源與補給模式未發(fā)生改變。地下水是歷史時期東南部雅布賴山區(qū)降水徑流及山前古河道滲漏補給而來，馬金珠等[30]、黃天明等[31]利用環(huán)境示蹤劑研究認(rèn)為是更新世晚期至全新世早期周邊山區(qū)降水徑流補給的古水。

3.4 湖泊萎縮特征

圖10顯示了湖泊萎縮程度區(qū)的分布特征，為湖泊萎縮程度差分模擬結(jié)果。區(qū)內(nèi)東南部、東北部、中北部及西部湖泊面積萎縮程度較大，而中部、中東部、南部及北部湖泊面積萎縮程度較小。結(jié)合區(qū)域地質(zhì)構(gòu)造特征（見圖1（a）），凹陷區(qū)基底較深，地下水匯集于此，水量豐富，當(dāng)氣候干旱時，有地下水補給的湖泊萎縮較慢;相反，隆起區(qū)湖泊對氣候干旱、水量減少反應(yīng)更為敏感。圖10中Ⅰ區(qū)、Ⅱ區(qū)、Ⅲ區(qū)分別位于雅布賴山隆起、因格井南凸起及宗乃山隆起帶，整體萎縮程度較大。Ⅳ區(qū)、Ⅴ區(qū)、Ⅵ區(qū)屬于蘇亥圖凹陷區(qū)，其中：Ⅴ區(qū)為沙漠大湖分布區(qū)，核心區(qū)為蘇木吉林湖，位于蘇亥圖凹陷區(qū)中部，萎縮程度最為緩慢;Ⅳ區(qū)、Ⅵ區(qū)萎縮程度略大于Ⅴ區(qū)，表現(xiàn)出凹陷區(qū)中心至邊界基底深度的漸變（Ⅵ區(qū)）與局部變化（Ⅳ區(qū)），造成湖泊萎縮程度有差異。Ⅶ區(qū)位于因格井凹陷南端，相對較深的基底保證了更厚的含水層，湖泊萎縮程度相對較小，地下水接受其東側(cè)的宗乃山山區(qū)降水徑流補給。

4 結(jié) 論

（1）巴丹吉林沙漠現(xiàn)代湖泊湖面高程平均低于古湖面9.76 m，湖泊總面積萎縮減小61.052 km2，占古湖總面積的75.49%，湖泊水量共減少4.9億m3，說明湖泊萎縮程度較高，沙漠干旱化趨勢明顯。

（2）湖泊區(qū)古、今等水位線均表現(xiàn)出東南高、西北低的宏觀特征，古、今水力梯度變化不明顯，反映了長時間尺度地下水有相同的補給來源，可能是周邊山區(qū)特別是東南部雅布賴山區(qū)古降水徑流及山前古河道下滲補給。

（3）湖泊萎縮的差異性受區(qū)域內(nèi)構(gòu)造基底的控制與影響，凹陷區(qū)萎縮程度相對較低，是大湖的主要分布區(qū)。

參考文獻(xiàn)：

[1] YANG X P， WILLIAMS M A J. The Ion Chemistry of Lakes and Late Holocene Desiccation in the Badain Jaran Desert， Inner Mongolia， China[J]. Catena， 2003，51（1）：45-60.

[2] ANDREOTTI B， FOURRIERE A， OULD-KADDOU F， et al. Giant Aeolian Dune Size Determined by the Average Depth of the Atmospheric Boundary Layer[J]. Nature， 2009，457：1120-1123.

[3] DONG Z B， QIAN G Q， LYU P， et al. Investigation of the Sand Sea with the Tallest Dunes on Earth： Chinas Badain Jaran Sand Sea[J]. Earth-Science Reviews， 2013，120：20-39.

[4] CHEN J S， LI L， WANG J Y， et al. Groundwater Maintains Dune Landscape[J]. Nature，2004，432：459-460.

[5] 陳建生， 趙霞， 盛雪芬， 等.巴丹吉林沙漠湖泊群與沙山形成機理研究[J].科學(xué)通報，2006，51（23）：2789-2796.

[6] YANG X P， SCUDERI L A. Hydrological and Climatic Changes in Deserts of China Since the Late Pleistocene[J]. Quaternary Research，2010，73（1）：1-9.

[7] GATES J B， EDMUNDS W M， DARLING W G， et al. Conceptual Model of Recharge to Southeastern Badain Jaran Desert Groundwater and Lakes from Environmental Tracers[J]. Applied Geochemistry， 2008，23（12）：3519-3534.

[8] GATES J B， EDMUNDS W M， MA J Z， et al. Estimating Groundwater Recharge in a Cold Desert Environment in Northern China Using Chloride[J]. Hydrogeology Journal， 2008，16（5）：893-910.

[9] MA J Z， EDMUNDS W M. Groundwater and Lake Evolution in the Badain Jaran Desert Ecosystem， Inner Mongolia[J]. Hydrogeology Journal，2006，14（7）：1231-1243.

[10] MA J Z， DING Z， GATES J B， et al. Chloride and the Environmental Isotopes as the Indicators of the Groundwater Recharge in the Gobi Desert， Northwest China[J]. Environmental Geology，2008，55（7）：1407-1419.

[11] WANG Z T， CHEN T Y， LIU S W， et al. Aeolian Origin of Inter Dune Lakes in the Badain Jaran Desert， China[J]. Arabian Journal of Geosciences，2016，9（3）：190.

[12] 丁宏偉，郭瑞，藍(lán)永超，等.再論巴丹吉林沙漠湖泊水的補給來源、補給模式與高大沙山的形成機理[J].冰川凍土，2015，37（3）：783-792.

[13] 張競，王旭升，胡曉農(nóng)，等.巴丹吉林沙漠地下水流場的宏觀特征[J].中國沙漠，2015，35（3）：774-782.

[14] 王海波，馬明國.基于遙感的湖泊水域動態(tài)變化監(jiān)測研究進(jìn)展[J].遙感技術(shù)與應(yīng)用，2009，24（5）：674-684.

[15] 朱金峰，王乃昂，李卓侖，等.巴丹吉林沙漠湖泊季節(jié)變化的遙感監(jiān)測[J].湖泊科學(xué)，2011，23（4）：657-664.

[16] 張振瑜，王乃昂，馬寧，等.近40 a巴丹吉林沙漠腹地湖泊面積變化及其影響因素[J].中國沙漠，2012，32（6）：1743-1750.

[17] 張振瑜，王乃昂，吳月，等.1973—2010年巴丹吉林沙漠腹地湖泊面積空間變化的遙感分析[J].湖泊科學(xué)，2013，25（4）：514-520.

[18] ZHANG Z， DONG Z B， YAN C Z， et al. Change of Lake Area in the Southeastern Part of Chinas Badain Jaran Sand Sea and Its Implications for Recharge Sources[J]. Journal of Arid Land，2015，7（1）：1-9.

[19] JIAO J J， ZHANG X T， WANG X S. Satellite-Based Estimates of Groundwater Depletion in the Badain Jaran Desert， China[J]. Scientific Reports，2015，5（1）：8960.

[20] GONG Y P， WANG X S， HU B X， et al. Groundwater Contributions in Water-Salt Balances of the Lakes in the Badain Jaran Desert， China[J]. Journal of Arid Land，2016，8（5）：694-706.

[21] WANG N A， NING K， LI Z L， et al. Holocene High Lake-Levels and Pan-Lake Period on Badain Jaran Desert[J]. Science China：Earth Sciences，2016，59（8）：1633-1641.

[22] 朱金峰，王乃昂，陳紅寶，等. 基于遙感的巴丹吉林沙漠范圍與面積分析[J].地理科學(xué)進(jìn)展，2010，29（9）：1087-1094.

[23] 王乃昂，馬寧，陳紅寶，等.巴丹吉林沙漠腹地降水特征的初步分析[J].水科學(xué)進(jìn)展，2013，24（2）：153-160.

[24] 張克存，姚正毅，安志山，等.巴丹吉林沙漠及其毗鄰地區(qū)降水特征及風(fēng)沙環(huán)境分析[J].中國沙漠，2012，32（6）：1507-1511.

[25] 胡文峰，王乃昂，趙力強，等.巴丹吉林沙漠典型湖泊湖氣界面水-熱交換特征[J].地理科學(xué)進(jìn)展，2015，34（8）：1061-1071.

[26] 劉建利，申安斌，陳小龍.大地電磁測深方法在內(nèi)蒙古西部銀根-額濟納旗盆地石炭系—二疊系油氣地質(zhì)調(diào)查中的應(yīng)用[J].地質(zhì)通報，2011，30（6）：993-1000.

[27] 嚴(yán)云奎，袁炳強，楊高印，等.內(nèi)蒙古西部銀根-額濟納旗盆地重力場與斷裂構(gòu)造的特征[J].地質(zhì)通報，2011，30（12）：1962-1968.

[28] 張國慶，XIE Hongjie，姚檀棟，等.基于ICESat和Landsat的中國十大湖泊水量平衡估算[J].科學(xué)通報，2013，58（26）：2664-2678.

[29] 張競，王旭升，賈鳳超，等.對內(nèi)蒙古阿拉善西部地下水流向問題的新認(rèn)識[J].現(xiàn)代地質(zhì)，2015，29（1）：213-219.

[30] 馬金珠，黃天明，丁貞玉，等.同位素指示的巴丹吉林沙漠南緣地下水補給來源[J].地球科學(xué)進(jìn)展，2007，22（9）：922-930.

[31] 黃天明，龐忠和.應(yīng)用環(huán)境示蹤劑探討巴丹吉林沙漠及古日乃綠洲地下水補給[J].現(xiàn)代地質(zhì)，2007，21（4）：624-631.

【責(zé)任編輯 張華興】