高吉喜，王永財，侯 鵬，萬華偉，張文國

（生態(tài)環(huán)境部衛(wèi)星環(huán)境應(yīng)用中心，北京 100094）

1 研究背景

黃河是中華民族的母親河，也是我國第二長河。黃河流域在我國經(jīng)濟社會發(fā)展和生態(tài)安全方面具有十分重要的地位，2019年黃河流域涉及的地級行政區(qū)人口總量為3.24億，占全國人口的23.3%，2018年的國內(nèi)生產(chǎn)總值（GDP）為19.4 萬億元，占全國GDP的比重為21.55%[1]。作為我國主要的能源基地，黃河流域煤炭資源占全國的70%，石油儲量占50%。2019年9月18日黃河流域生態(tài)保護和高質(zhì)量發(fā)展座談會上，習(xí)近平總書記的講話中把黃河流域的生態(tài)保護和高質(zhì)量發(fā)展同京津冀協(xié)同發(fā)展、長江經(jīng)濟帶發(fā)展、粵港澳大灣區(qū)建設(shè)、長三角一體化發(fā)展一樣，上升為重大國家戰(zhàn)略[2]，進一步說明了黃河流域在國家安全中的重要地位。黃河流域以其占全國2.2%的徑流量，承擔(dān)著占全國15%的耕地和12%人口的供水任務(wù)[3-4]。在氣候變暖和人類活動的雙重影響下，黃河流域水資源的供需矛盾更加突出，已成為我國水資源極其短缺的地區(qū)之一[5-6]，并在未來仍將繼續(xù)面臨水資源嚴(yán)重短缺的嚴(yán)峻挑戰(zhàn)[7-8]。在此背景下，堅持黃河流域生態(tài)保護和高質(zhì)量發(fā)展，對黃河水資源綜合開發(fā)利用和保護提出了更高的要求。為此，快速、及時、準(zhǔn)確地監(jiān)測黃河流域陸表水域分布和時空變化，對黃河流域水域水資源管理、生態(tài)保護政策的制定等有著重要的意義。

遙感技術(shù)是目前水資源調(diào)查常用的技術(shù)手段之一，具有覆蓋面廣、時效性強等優(yōu)勢，為大范圍水體變化監(jiān)測提供了全新的技術(shù)手段，光學(xué)遙感、微波雷達等技術(shù)的不斷完善，遙感數(shù)據(jù)水體提取算法日趨完善[9-10]。同時，基于地理大數(shù)據(jù)平臺在水體中的應(yīng)用日益廣泛[11-12]，使得遙感在水體變化、洪水淹沒監(jiān)測等方面發(fā)揮越來越重要的作用[13-23]。本文利用遙感技術(shù)，分析2000—2019年黃河流域陸表水域分布變化，以期為黃河流域水資源保護與利用的決策提供依據(jù)。

2 研究區(qū)概況、資料與方法

2.1 研究區(qū)概況黃河流域包括青海、四川、甘肅、寧夏、內(nèi)蒙古、陜西、山西、河南、山東9個省份，全長5464 km，流域面積79.6萬km2，其中上游36.5萬km2，中游40.8萬km2，下游2.3萬km2。黃河上游是指從河源到內(nèi)蒙古自治區(qū)托克托縣河口鎮(zhèn)以上的黃河河段，主要位于青藏高原、寧蒙河套地區(qū)，但在提及上游陸地面積時，通常還包含鄂爾多斯地表內(nèi)流區(qū)。為此，黃河上游區(qū)域覆蓋青海省、甘肅省、寧夏回族自治區(qū)、內(nèi)蒙古自治區(qū)。黃河中游是指從河口鎮(zhèn)到河南鄭州花園口的河段，包括陜西省、山西省、內(nèi)蒙古自治區(qū)、甘肅省、寧夏回族自治區(qū)、河南省。黃河下游指黃河干流自花園口到進入渤海的區(qū)域，包括河南省、山東省，如圖1所示。黃河流域水資源二級分區(qū)總共8個，包括龍羊峽以上、龍羊峽至蘭州、蘭州至河口鎮(zhèn)、河口鎮(zhèn)至龍門、龍門至三門峽、三門峽至花園口、花園口以下、內(nèi)流區(qū)，如圖1所示。

圖1 黃河流域主要水系及河段分布

2.2 研究資料與方法為獲取研究區(qū)近20年水體信息，本文使用了研究區(qū)2000—2018年長時間序列的Landsat5、Landsat7、Landsat8的表面反射率Tier 1產(chǎn)品，該產(chǎn)品可以在谷歌地球引擎（Google Earth Engine，GEE）上進行查詢、可視化和分析。通過使用Landsat 5、Landsat7和Landsat8 表面反射代碼（LarSRC）校正了大氣散射和氣溶膠吸收，包括云、陰影和每個像素的飽和度，獲得了黃河流域2000—2019年共45 424景影像。水體提取具體過程如下。

首先，通過SRTM高程數(shù)據(jù)計算流域相對高程，具體計算方法可參考Donchyts等[24]的計算方法。基于Landsat表面反射率數(shù)據(jù)計算植被指數(shù)、水體指數(shù)、干裸指數(shù)，具體各個指數(shù)的計算方法如下：

式中：NDVI為歸一化植被指數(shù)；MNDWI為改進的歸一化差異水體指數(shù)；AWEI為自動水體提取指數(shù)；DBSI為干裸指數(shù)；ρblue為藍波段；ρgreen為綠波段；ρNIR為近紅外波；ρred為紅波段；ρSWIR為短波紅外波段；ρSWIR1為短波紅外波段1；ρSWIR2為短波紅外波段2。

其次，基于遙感水體指數(shù)、流域相對高程Hand根據(jù)閾值提取水體，并結(jié)合植被指數(shù)、干裸指數(shù)過濾非水體。

最后，利用每年每期水體數(shù)據(jù)構(gòu)建水淹頻度（一年中某一地點被水淹沒的次數(shù)占有效觀測次數(shù)的比例），為了分析不同年份水域面積動態(tài)變化，利用62.0%的閾值提取水域分布，相當(dāng)于一年62.0%時間淹水的范圍。經(jīng)過抽樣和人工目視與每年多期遙感影像比對，提取水體分布與影像中水體分布一致。

3 結(jié)果與分析

3.1 流域陸表水域面積空間分布及變化特征流域陸表水域空間面積分布分析表明，2014—2019年期間，黃河流域的陸表水域面積平均為6705.84 km2。從不同河段來看，河口鎮(zhèn)以上的上游地區(qū)陸表水域面積平均為4046.03 km2，約占整個流域陸表水域面積的比例為60.3%；龍羊峽以上黃河源區(qū)的陸表水域面積平均為2572.20 km2，約占整個流域陸表水域面積的比例為38.4%；河口鎮(zhèn)至花園口的中游地區(qū)陸表水域面積平均為1965.38 km2，約占整個流域陸表水域面積的比例為29.3%；花園口以下的下游地區(qū)陸表水域面積平均為694.44 km2，約占整個流域陸表水域面積的比例為10.4%。從不同省份來看，青海省和內(nèi)蒙古自治區(qū)的比例最大，四川省的比例最??；青海省的陸表水域面積平均為2701.63 km2，占流域陸表水域面積的比例為40.3%；內(nèi)蒙古自治區(qū)的陸表水域面積平均為1137.47 km2，占整個流域陸表水域面積的比例為17.0%；四川省的陸表水域面積平均為52 km2，占流域陸表水域面積的比例為0.8%。2000—2019年間，黃河流域的陸表水域面積總體呈增加趨勢、約為103 km2/a，2015年之后，陸表水域面積增加趨勢更加明顯，如圖2所示。

圖2 黃河流域2000—2019年水域面積變化

從流域陸表水域面積變化來看，2019年，陸表水域面積為7532.06 km2，比2000年增加了約2475.68 km2，增加的面積相當(dāng)于半個青海湖面積，增幅約為49.0%。從不同河段來看，上游地區(qū)的陸表水域面積每年增加約54 km2，近20年間陸表水域面積共增加約1181.50 km2、增加幅度為34.7%，對流域的陸表水域面積增加貢獻度為47.7%；中游地區(qū)的陸表水域面積每年增加約44 km2，近20年間陸表水域面積共增加1058.50 km2，對流域的陸表水域面積增加貢獻度為42.7%；下游地區(qū)的陸表水域面積每年增加約7.5 km2，近20年間陸表水域面積共增加約235.68 km2，對流域的陸表水域面積增加貢獻度為9.6%。從不同省份來看，陸表水域面積增加最大的省份是青海省和內(nèi)蒙古自治區(qū)，增幅最大的是山西??；青海省的陸表水域面積增加約619.80 km2，對流域的陸表水域面積增加貢獻度為25.7%；內(nèi)蒙古自治區(qū)的陸表水域面積增加580.63 km2，對流域的陸表水域面積增加貢獻度為24.1%；山西省的陸表水域面積增加幅度為108.28%，對流域的陸表水域面積增加貢獻度為11.3%。如表1所示。

表1 黃河流域陸表水域面積變化及對流域水域面積變化的貢獻度

圖3 2019年相比2000年湖庫水域面積變化

圖4 黃河源區(qū)水域面積與年降水變化關(guān)系

3.2 流域內(nèi)主要湖庫水域面積變化特征與2000年相比，2019年水域面積大于1 km2的湖庫中有217個發(fā)生了顯著變化，空間分布如圖3所示。其中，水域面積擴大的有115個，位于上游的66個、中游的40個、下游的9個；水域面積縮小的有45個，位于上游的26個、中游的16個、下游的3個；新增加的湖庫有57個，位于上游的24個、中游的31個、下游的2個。

上游地區(qū)的湖庫多數(shù)是高原湖泊，特點是面積小、數(shù)量多，兩個面積較大湖泊為扎陵湖和鄂陵湖。龍羊峽以上的黃河源區(qū)超過1 km2的湖庫中，面積發(fā)生變化的有46個，其中面積擴大的湖庫有45個，這種變化與降水變化趨勢基本一致，如圖4所示。龍羊峽-龍門鎮(zhèn)的干旱半干旱區(qū)河段的面積大于1 km2湖庫中，水域面積縮小的有33個；中游內(nèi)流區(qū)是鄂爾多斯高原鹽湖的聚集區(qū)，面積大于1 km2的32個湖庫中，面積縮小的有8個。通過對內(nèi)流區(qū)重點53個湖泊進行分析，發(fā)現(xiàn)有35個存在不同程度的開發(fā)活動，無序的開發(fā)活動和水開采可能加劇區(qū)域水資源緊張程度。2000—2019年間，下游地區(qū)的湖庫水面積增加以人工水庫、水塘為主，大于1 km2的新增加和面積擴大的人工水域有11個。

作為源頭區(qū)的重要湖泊，近20年間扎陵湖和鄂陵湖水域面積增加80 km2，入湖徑流量增加是扎陵湖、鄂陵湖水面擴大的主要原因。紅堿淖濕地作為重要的沙漠淡水湖和遺鷗繁殖棲息地，近20年間濕地面積總體在下降，2000—2015年間減少約34%。同時鄂爾多斯?jié)竦孛娣e也在減少，減少約85%。初步分析表明，導(dǎo)致濕地面積下降的主要原因有氣候持續(xù)暖干、周邊水庫建設(shè)和煤礦開采共同導(dǎo)致的補給河流被截斷，地下水層遭到破壞[25-28]，2015年采取補水措施之后紅堿淖濕地和鄂爾多斯?jié)竦厮蛎娣e逐年增加。作為國家南水北調(diào)東線工程重要樞紐，東平湖水域面積在2014年以前先增后減，2015年之后水域面積持續(xù)增加，2019年的水域面積增加了2.4 km2。如圖5所示。

圖5 黃河流域典型湖庫水域面積變化

3.3 流域內(nèi)主要子流域的變化特征作為黃河生態(tài)保護與高質(zhì)量發(fā)展的關(guān)注區(qū)域，渭河和汾河子流域的水文特征和水環(huán)境特征在近20年間均發(fā)生了很大變化。在渭河流域，2000—2019年，渭河流域陸表水域面積擴大了113.17 km2，增幅約為139.9%，約有62.3%區(qū)域內(nèi)植被覆蓋度呈明顯增加趨勢。根據(jù)水文統(tǒng)計年鑒資料，渭河流域華縣水文監(jiān)測站的含沙量連續(xù)多年下降，下降幅度達到67.1%。在汾河流域，2000—2019年，汾河流域陸表水域面積擴大了55.75 km2、增幅約152.4%，約有69.7%區(qū)域植被覆蓋度呈明顯增加趨勢。根據(jù)水文統(tǒng)計年鑒，流域內(nèi)河津站水文監(jiān)測的含沙量連續(xù)多年下降、下降幅度達到65.0%。參見圖6、圖7和圖8。

圖6 2000—2019年渭河流域水域面積和含沙量變化

圖7 2000—2019年汾河流域水域面積和含沙量變化

圖8 2000—2019年渭河和汾河流域植被覆蓋度變化

4 主要結(jié)論

本文基于衛(wèi)星遙感數(shù)據(jù)提取了2000—2019年黃河流域陸表水域分布，開展了陸表水域面積的動態(tài)變化監(jiān)測，以期為黃河流域水資源保護與利用的決策提供依據(jù)，對于了解全球變化和人類活動對水資源的影響具有重要意義。

（1）流域陸表水域空間分布不均衡，空間分布主要集中在上游。近5年，黃河流域陸表水域平均面積為6705.84 km2。作為黃河流域水源涵養(yǎng)區(qū)和干流源頭區(qū)，上游陸表水域面積約為4046.03 km2，占黃河流域陸表水域總面積的60.3%；作為黃河流域重要水土保持區(qū)，中游地區(qū)陸表水域面積約為1965.38 km2，占黃河流域陸表水域總面積的29.3%。下游地區(qū)陸表水域面積約為694.44 km2，占黃河流域陸表水域總面積的10.4%。從不同省份來看，黃河流域陸表水域面積的比例，青海省最大、內(nèi)蒙古次之、四川省最小。

（2）近20年流域陸表水域面積呈動態(tài)增加趨勢。2019年流域陸表水域總面積為7532.06 km2，比2000年增加了2475.68 km2，增幅為49.0%。其中，上游黃河源區(qū)湖庫面積整體在擴大，中上游干旱半干旱區(qū)部分湖泊面積在萎縮，下游人工庫塘在擴大。2000—2019年期間，上游地區(qū)對流域陸表水域總面積增加貢獻度達到47.7%，中游地區(qū)貢獻度為42.7%，下游地區(qū)貢獻度為9.6%。從省域來看，青海省的流域陸表水域面積增加最大，增加約619.80 km2，對流域陸表水域總面積增加的貢獻度為25.7%。

（3）2000—2019年期間，渭河和汾河流域的陸表水域面積分別擴大了113.17 km2和55.75 km2，增幅約為139.9%和152.4%。近年來，渭河和汾河流域有60.0%以上地區(qū)的植被覆蓋度呈現(xiàn)出明顯的改善，導(dǎo)致含沙量持續(xù)下降，下降幅度分別為67.1%和65.0%。