近70 a 黃河輸沙量影響機制的定量評估

來源:用戶上傳      作者:王佳佳 王厚杰 吳曉 盧泰安

　　摘要：基于1951―2017年黃河中游實測降水量、徑流量以及輸沙量數據，采用克里金插值、Mann－Kendall 突變檢驗等方法定量評估人類活動和氣候變化對不同階段黃河輸沙量減少的貢獻率。研究發現，在第二階段（1969―1985年）、第三階段（1986―1999年）、第四階段（2000―2017年）人類活動對黃河輸沙量減少的貢獻率分別為76％、67％、97％。人類活動對黃河輸沙量減少的貢獻率較大，這表明流域內人類活動已成為控制沉積物向海洋輸送的主要力量。
　　關鍵詞：降水變化；人類活動；輸沙量；貢獻率；黃河
　　中圖分類號：P333.4；P338＋.5 文獻標志碼：A doi：10.3969／j.issn.1000－1379.2022.10.005
　　引用格式：王佳佳，王厚杰，吳曉，等.近70 a 黃河輸沙量影響機制的定量評估[J].人民黃河，2022，44（10）：25－30.
　　Quantitative Evaluation of the Mechanism on Sediment Transport in the Yellow River Basin over the Recent 70 Years
　　WANG Jiajia1，WANG Houjie1，2，WU Xiao1，LU Taian1
　?。?.College of Marine Geoscience ，Ocean University of China，Qingdao 266100，China；2.Key Laboratory of Submarine Geoscience and Prospecting Technology Ministry of Education，Qingdao 266100，China）
　　Abstract：Based on the measured precipitation，runoff and sediment transport data in the Middle Yellow River from 1951 to 2017，Kriging interpolation，Mann?Kendall catastrophe test and other methods were used to quantitatively evaluate the contribution rate of human activitiesand climate change to the reduction of sediment transport in different stages of the basin.The results indicate that the contributions of humanactivities on the Yellow River sediment reduction at different stages （1969－1985，1986－1999 and 2000－2017） was 76％，67％ and 97％re?spectively.The contribution rate of human activities to the decline of the sediment transport from the Yellow River into the sea is relatively larger，which suggests that human activity in the river basin has become a primary force that dominates the sediment delivery to the ocean.
　　Key words：precipitation change；human activity；sediment transport；contribution rate；Yellow River
　　河流作槁皆次鎦氏蠔Ｑ笫湓說鬧饕通道，在全球生物地球化學循環中發揮著關鍵作用[1－2]。黃河以其高含沙量著稱，每年向海洋輸送約10.8億 t 泥沙，但近幾十年來受氣候變化和人類活動的影響，輸沙量的時空分布發生顯著變化，科學評估其輸沙量的變化趨勢、厘清氣候變化和人類活動對其輸沙量的影響，不僅對認識全球變化背景下河流系統變化具有重要意義，還可以為實現黃河流域生態保護和高質量發展、河口三角洲生態環境可持續發展提供重要參考[3]。國內外眾多學者針對氣候變化和人類活動對黃河輸沙量的影響已開展大量研究，如：Wang 等[4]通過研究1950―2005年黃河入海泥沙變化過程，估算1950年以來氣候變化和人類活動對黃河入海泥沙量減少的貢獻率分別為30％和70％；劉通等[5]研究黃河流域內蒙古河段時發現，1990―2010年西柳溝區域分別有101.00 km2和158.10 km2的未利用地和沙地轉變為林地和草地，導致該時期內蒙古河段年輸沙量顯著下降。但上述研究中或選取時間跨度較短，或未對人類活動對年輸沙量的貢獻進行定量評估。本文系統收集1951―2017年黃河流域降水、徑流和泥沙資料，研究年代際尺度上降水變化和人類活動對輸沙量的影響，通過建立相關模型，定量估算不同階段降水變化和人類活動對黃河輸沙量減少的貢獻率。
　　1 研究區概況
　　黃河全長約5464 km，分為上游、中游及下游：上游河段從黃河起源地到河口鎮，河道全長3472 km，流域面積42.8萬 km2；中游河段從河口鎮到桃花峪，河道全長1206 km，流域面積34.4萬 km2，中游河段流經黃土高原，是黃河流域的主要產沙區；下游河段從桃花峪到入海口，河道全長786 km，流域面積2.3萬 km2[6]。

nlc202212201221

　　2 數據來源與研究方法
　　2.1 數據來源
　　由于黃河泥沙來源主要是黃土高原地區水土流失導致的輸沙入河，因此著重選取1951―2017年黃河中游及周邊14個重要氣象站點（分布在包頭、呼和浩特、榆林、太原、石家莊、延安、臨汾、長治、西峰鎮、西安、漢中、鄭州、南陽、鄖縣）的降水量數據，數據源自中國氣象局國家氣候中心。1951―2017年黃河上、中、下游主要水文站（三門峽、花園口、利津等）的泥沙數據源自國家地球系統科學數據共享服務平臺黃土高原科學數據中心以及水利部黃河水利委員會發布的《黃河泥沙公報》，其中：1951―2000年泥沙數據源自黃土高原科學數據中心，2001―2017年泥沙數據源自《黃河泥沙公報》。
　　2.2 研究方法
　　2.2.1 克里金插值法
　　克里金插值法是依據協方差函數對隨機過程／隨機場進行空間建模和預測（插值）的回歸算法，在特定的隨機過程中給出最優線性無偏估計[7]。假定插值數據的均值未知，待插入點 x0的插值數據 z （ x0）的計算公式為
　　式中：n 為實測點的總數；xi 為1～ n 范圍內的任一實測點；z （ xi ）為點 xi 的屬性值；λi為插值點與實測點 xi 的距離權重[8]。
　　2.2.2 泰森多邊形法
　　計算流域氣象數據平均值常采用泰森多邊形法，這是一種基于離散型氣象站點數據的加權平均方法。其思想是假設任何一個觀測站的測量值可以代表一定范圍內的降水量數據，這意味著任何區域的氣象數據都等于距其最近的多個觀測點的實測氣象值按照不同權重加和的結果[9]。
　　2.2.3 Mann－Kendall 突變檢驗
　　Mann－Kendall 突變檢驗能夠揭示時間序列數據的變化趨勢，該方法中數據無須遵循特定的分布規律，適用于分布異常的數據[10]，表達式[11－12]如下：
　　式中：S 為統計量；xj為時間序列的第j個數據值；xk為時間序列的第 k 個數據值；n 為數據樣本的長度； Var （S）為 S 的方差；Z 為標準化統計量，顯著性水平α＝0.05，即在95％的顯著性水平下檢驗 Z 值。若 Z 為負值，其絕對值大于顯著性水平，則趨勢為明顯遞減；若 Z 為正值，其絕對值大于顯著性水平，則趨勢為明顯遞增；若 Z 的絕對值小于顯著性水平，則不存在趨勢。
　　2.2.4 ENSO 事件強度判斷
　　已有研究表明，發生在赤道西太平洋的海表溫度出現周期性異常升溫或降溫的現象（稱為厄爾尼諾現象）與發生在太平洋、印度洋之間的反向氣壓振動現象（稱為南方濤動）實際上是相互聯系的，是同一種現象在海洋和大氣這兩種介質中的不同表現，因此將厄爾尼諾現象（ El Nino ）和南方濤動（ Southern Oscillation）合稱為 ENSO 事件[13]。ENSO 事件的發生與降水量的相關性顯著[14]，本文采用 Nino3.4區海溫距平指數作為判定 ENSO 事件的依據，美國海洋大氣局氣候預測中心的判定標準為海溫距平指數≥0.5℃且持續5個月以上定義為一次 ENSO 暖事件（厄爾尼諾事件），海溫距平指數≤－0.5℃且持續5個月以上定義為一次 ENSO 冷事件（拉尼娜事件）。另分別采用事件發生期間海溫距平指數極值（峰值或谷值）與這個事件持續時間內海溫距平指數的累加值作為 ENSO 事件強度等級的劃分指標，對2個指標進行標準化處理后將事件強度劃分為5個等級，見表1[15]。
　　3 黃河中游降水量變化特征
　　3.1 黃河中游降水量年際、月際變化特征
　　1951―2017年黃河中游多年平均降水量為493 mm，降水量年際、年代際變化趨勢見圖1。降水量呈現明顯的時間分布差異，自20世紀50年代以來，黃河中游年降水量變化幅度較大，最大年降水量可達750 mm，最小年降水量僅為321 mm；1970―1979年、1980―1989年、1990―1999年多年平均降水量明顯降低；2000―2009年、2010―2017年多年平均降水量明顯上升。
　　降水量月際變化趨勢見圖2，每年7―10月為汛期，汛期降水量占全年降水量的65％，且降水量的峰值主要出現在7月和8月。
　　3.2 極強 ENSO 事件影響下黃河中游的降水量變化特征
　　根據 ENSO 事件強度劃分標準確定自1950年以來，黃河流域共發生15次厄爾尼諾事件和20次拉尼娜事件（ La Nina），其中包括3次極強暖事件（1982―1983年、1997―1998年、2015―2016年）、2次極強冷事件（1954―1957年、1973―1974年）。1950―2017年黃河流域年降水量、年徑流量（利津水文站實測徑流量）的時間序列圖見圖3（圖中氣泡大小與 ENSO 事件強度相對應），可以看出，ENSO 事件發生年份對應某時段內年降水量極值的出現年份，即絕大多數厄爾尼諾事件發生年份對應年降水量峰值的出現年份，部分厄爾尼諾事件和拉尼娜事件發生年份對應年降水量谷值的出現年份，據此可推斷出ENSO事件確實對黃河流域的降水量變化起到一定作用。自 20世紀70年代之后，ENSO 事件出現的頻率明顯提高，分析降水量數據可知，1970―2017年黃河流域多年平均降水量只有481 mm，而1970年以前流域多年平均降水量約為522 mm，1970年以后減少了約8％。ENSO 事件與流域降水量的相關性表明，流域降水受全球氣候的強烈影響，海溫距平指數的變化是由赤道地區的東北信風引起的，信風將水汽輸送回大氣，改變了水汽的時空分布格局，最終改變了黃河流域的區域降水格局，從而影響了黃河流域的年徑流量。
　　為進一步研究強 ENSO 事件對黃河流域降水量的影響，選取1950年以來正常時間（無事件l生）、2次極強暖事件（ EL1和 EL2）發生時間和2次極強冷事件（ LA1和 LA2）發生時間為3個典型時段，分析黃河流域40 mm 月降水量等值線在3個時段內的變化情況（見圖4）?？梢钥闯?，強厄爾尼諾事件（極強暖事件）造成40 mm月降水量等值線明顯向東南方向移動，拉尼娜事件（極強冷事件）造成40 mm 月降水量等值線向東北方向移動。尤其是黃河中游地區，極強暖事件造成40 mm 月降水量等值線移至延河甚至北洛河以南，而極強冷事件造成40 mm月降水量等值線移至延河以北，使黃河中游大部分區域處于豐水帶。另對黃河中游月平均降水量數據進行分析發現，正常事件發生月份中游地區月平均降水量在25～45 mm 范圍內，極強暖事件發生月份中游地區月平均降水量在25～35 mm 范圍內，極強冷事件發生月份中游地區月平均降水量在35～55 mm 范圍內，極強暖事件造成黃河中游地區月平均降水量減少約10 mm，而極強冷事件造成黃河中游地區月平均降水量增加約10 mm 。綜上，極強暖事件造成黃河中游降水量的減少，而極強冷事件對黃河中游降水量具有增加作用。

nlc202212201221

　　4 黃河流域主要人類活動
　　4.1 水利樞紐工程建設
　　水庫作為重要水利設施，通過汛期攔蓄洪水、非汛期泄放興利流量，調平年內流量過程，實現黃河水資源的有效分配。同時水庫作為重要的攔調泥沙工程，可有效緩解下游河道淤積，還可通過人造洪峰、水沙調節等方式沖刷河床、輸沙入海，從而長期發揮減淤作用，因此水庫建設一直是黃河流域防洪減淤的重要舉措。據《黃河水資源公報》統計，黃河流域共有大中小型水庫3147座，其中：大型水庫34座，中型水庫185座，大中型總庫容可達574億 m3，相當于1950―2000年黃河流域的年平均天然徑流量。所有水庫中對黃河流域影響較大的4座水庫分別是三門峽、劉家峽、龍羊峽以及小浪底水庫，這4座水庫的總庫容占黃河流域總庫容的97％[16]。
　　4.2 水土保持
　　為加快恢復黃土高原生態環境、保護水土資源，國家啟動實施了針對黃土高原水土流失治理的多個項目工程。目前，水土保持的措施有水土保持農業技術措施、水土保持林草措施和水土保持工程措施三大類[17]。度陽[18]關于黃土高原退耕還林還草工程效果的研究表明，2000―2015年耕地比重凈減少2.45％，林地比重凈增加0.7％，草地比重凈增加0.18％。水土保持工程措施是在小流域內建造的水土攔截工程，包括山坡防護工程、山溝治理工程以及小型蓄水引水工程。根據水利部黃河水利委員會2010年發布的《黃河流域水土保持公報》，截至2010年，黃河流域已累計治理水土流失面積22.56萬 km2，建設淤地壩9.1萬座，各類小型水土保持工程措施實施區域184萬處，入河泥沙減少量超過3億 t，水土保持工程效果顯著。
　　4.3 引水灌溉
　　黃河作為北方工農業生產發展供水的主要源區，引水灌溉工程的建設一直受到國家高度重視。據水利部黃河水利委員會2017年統計，黃河流域及其下游引黃灌區有大型灌區（規模大于0.7萬 hm2）87處、中型灌區663處，灌溉面積達0.08億 hm2，比新中國建立初期增長近10倍。
　　5 輸沙量影響機制的定量評估
　　5.1 突變性分析
　　1951―2017年黃河流域主要水文站輸沙量的 M－K突變檢驗結果見表2，主要水文站大都在1999年左右發生輸沙量的顯著突變，各水文站輸沙量在突變點后呈現明顯的減少趨勢，減少比例均在70％以上，其中花園口水文站輸沙量減少比例最大（為92％），蘭州水文站輸沙量減少比例最小（為72％）。
　　5.2 利津水文站輸沙量年際變化特征分析
　　1951―2017年利津水文站年降水量、年徑流量、年輸沙量的階段性變化見圖5，可以明顯看出，1968年、1985年、1999年為3個典型的分界點，具體表現為3個分界點之后的年徑流量和年輸沙量的峰值未曾高過分界點之前的年徑流量和年輸沙量的峰值。結合水利樞紐工程建成時間，這3個年份與劉家峽水庫（1968年）、龍羊峽水庫（1986年）以及小浪底水庫（1999年）的投入運用年份基本一致[19]。這表明以水利工程的建成時間為分界點，通過分析汛期（7―10月）黃河的降水量和輸沙量來研究不同時段氣候變化和人類活動對黃河輸沙量的貢獻率是合理的。
　　5.3 階段性定量評估
　　以劉家峽、龍羊峽和小浪底水庫的運用年份為分界點，把黃河輸沙量變化過程劃分成4個階段，分別為第一階段（1951―1968年）、第二階段（1969―1985年）、第三階段（1986―1999年）、第四階段（2000―2017年）。以第一階段為基準期，由于1960―1961年利津水文站部分月份輸沙量缺失，因此第一階段去掉這兩年數據缺失的異常值。建立降水量―輸沙量的統計模型：y＝0046.x－117（相關系數 R.2＝059）.，表明在人類活動影響較小的基準期內，降水量和輸沙量具有顯著相關性，在其他自然條件保持一定或變化較小時，降水量和輸沙量呈線性關系。定義侵蝕比 E 即可消除降水量對輸沙量的影響[20]，公式為
　　式中：S 為輸沙量；S′為除降水量之外的其他因素產生的輸沙量；P 為降水量。以基準期的多年平均輸沙量 S0（2.8億 t）、降水量 P0（86 mm）、侵蝕比 E0（0.046）為基準，定量評估人類活動和氣候變化對黃河輸沙量的貢獻率，其余各階段輸沙量與基準期輸沙量的差值ΔS 計算公式為
　　式中：Si 為不同階段的實測輸沙量；Ei為不同階段的侵蝕比；Pi 為不同階段的年平均降水量；ΔE 為不同階段侵蝕比與基準期侵蝕比的差值；P 為基準期與不同階段的平均降水量；E 為基準期與不同階段的平均侵蝕比。
　　當ΔE＝0時（侵蝕比保持不變），ΔS＝EΔP 代表降水量對輸沙量的貢獻量；當ΔP＝0時（降水量保持不變），ΔS＝PΔE 代表人類活動對輸沙量的貢獻量，降水量和人類活動對輸沙量的貢獻率分別為×100％、×100％。人類活動和氣候變化對黃河輸沙量的貢獻率見表3。
　　第二階段：汛期總降水量為5523 mm ，占全階段總降水量的67％，月平均降水量為81 mm；實測汛期總輸沙量為122億 t，占全階段總輸沙量的85％，月平均輸沙量為1.75億 t 。計算可得侵蝕比為0.038，人類活動和氣候變化對汛期輸沙量減少的貢獻量分別為0.668億、0.210億 t，人類活雍推候變化對輸沙量的貢獻率分別為76％、24％。
　　第三階段：汛期總降水量為3911 mm ，占全階段總降水量的62％，月平均降水量為70 mm；實測汛期總輸沙量為49億 t，占全階段總輸沙量的88％，月平均輸沙量為0.88億 t 。計算可得侵蝕比為0.030，人類活動和氣候變化對汛期輸沙量減少的貢獻量分別為1.248億、0.608億 t，人類活動和氣候變化對輸沙量的貢獻率分別為67％、33％。
　　第四階段：汛期總降水量為5916 mm ，占全階段總降水量的66％，月平均降水量為82 mm；實測汛期總輸沙量為14.56億 t，占全階段總輸沙量的71％，月平均輸沙量為0.20億 t 。計算可得侵蝕比為0.006，人類活動和氣候變化對汛期輸沙量減少的貢獻量分別為3.360億、0.104億 t，人類活動對輸沙量的貢獻率增大至97％，氣候變化的貢獻率僅為3％，表明在第四階段黃河輸沙量基本由人類活動控制，氣候變化的影響已經微乎其微。

nlc202212201221

　　綜上，與氣候變化相比，人類活動對黃河輸沙量減少的貢獻率較大，表明人類活動已成為控制沉積物向海洋輸送的主要力量。
　　6 結論
　　（1）黃河中游降水量呈明顯年際、年代際分布差異，起伏波動較大；強 ENSO 事件對降水量具有明顯的影響，極強暖事件造成黃河中游降水量減少，而極強冷事件則對黃河中游降水量具有增加作用。
　?。?）突變性分析和徑流量、輸沙量階段性分析結果表明，1968年、1985年、1999年為徑流泥沙變化的3個典型分界點，這3個年份與劉家峽（1968年）、龍羊峽（1986年）及小浪底水庫（1999年）的投入運用年份基本一致。
　　（3）1968年以來，人類活動對黃河輸沙量減少的貢獻始終占據較大比重，第二階段（1969―1985年）、第三階段（1986―1999年）、第四階段（2000―2017年）貢獻率分別為76％、67％、97％。
　　參考文獻：
　　[1] MILLIMAN J D ，MEADE R H.World?Wide Delivery of RiverSediment to the Oceans [J].The Journal of Geology，1983，91（1）：1－21.
　　[2] WALLING D E ，FANG D.Recent Trends in the SuspendedSediment Loads of the World，s Rivers [J].Global and Plane? tary Change ，2003，39（1）：111－126.
　　[3] WU X，WANG H J ，BI N S ，et al.Climate and Human Battlefor Dominance Over the Yellow River’s Sediment Discharge：From the Mid?Holocene to the Anthropocene [J].Marine Ge ? ology，2020，425：106188.
　　[4] WANG H J ，YANG Z S ，SAITO Y ，et al.StepwiseDecreasesof the Huanghe （ Yellow River ） Sediment Load （1950―2005）：Impacts of Climate Change and Human Activities[J].Global and Planetary Change ，2007，57（3）：331－354.
　　[5] 劉通，黃河清，邵明安，等.氣候變化與人類活動對鄂爾多斯地區西柳溝流域入黃水沙過程的影響[J].水土保持學報，2015，29（2）：17－22.
　　[6] WANG H J，YANG Z S，LI Y H，et al.Dispersal Pattern ofSuspended Sediment in the Shear Frontal Zone off theHuanghe （ Yellow River ） Mouth [J].Continental Shelf Research，2006，27（6）：854－871.
　　[7] LEILA N ，MAZDA K Z ，REZA S A ，et al.GroundwaterDepthand Elevation Interpolation by Kriging Methods in MohrBasin of Fars Province in Iran [J].Environmental Monitoring and Assessment，2010，166（1－4）：387－407.
　　[8] 罟α鰨張桂敏，李士心.泛克里金插值法在地磁圖中的應用[J].中國慣性技術學報，2008，16（2）：162－166.
　　[9] 柴元方，孫雪嵐，王超，等.基于泰森多邊形的多點法測 U形明渠流量[J].人民黃河，2016，38（5）：139－142.
　　[10] KISI O ，AY M.Comparison of Mann?Kendall and InnovativeTrend Method for Water Quality Parameters of the KizilirmakRiver，Turkey[J].Journal of Hydrology，2014，513：362－375.
　　[11] RIBEIRO L，KRETSCHMER N，NASCIMENTO J，et al.Evalua?ting Piezometric Trends Using the Mann?Kendall Test on the Alluvial Aquifers of the Elqui River Basin，Chile [J].Hydro? logical Sciences Journal，2015，60（10）：1840－1852.
　　[12] MANN H B.NonparametricTests Against Trend[J].Econo?metrica：Journal of the Econometric Society，1945，13（3）：245－259.
　　[13] WANG H J ，YANG Z S ，SAITO Y，et al.InterannualandSeasonal Variation of the Huanghe （ Yellow River ） WaterDischarge Over the Past 50 years：Connections to Impactsfrom ENSO Events and Dams [J].Global and Planetary Change ，2006（50）：212－225.
　　[14] 黃婷婷，林青霞，吳志勇，等.黃河流域干旱時空特征及其與 ENSO 的關聯性分析[J].人民黃河，2021，43（11）：52－58.
　　[15] 許武成，王文，馬勁松，等.1951―2007年的 ENSO 事件及其特征值[J].自然災害學報，2009，18（4）：18－24.
　　[16] 張佳.黃河中游主要支流輸沙量變化及其對入海泥沙通量的影響[D].青島：中國海洋大學，2011：20－21.
　　[17] WANG H J，WU X，BI N S，et al.Impacts of the Dam?OrientatedWater?Sediment Regulation Scheme on the Lower Reaches and Delta of the Yellow River （ Huanghe ）：A Review [J]. Global and Planetary Change ，2017，157：93－113.
　　[18] 度陽.黃土高原地區退耕還林（草）工程效果研究[D].咸陽：西北農林科技大學，2018：25－26.
　　[19] 柳莎莎.氣候變化和人類活動對現代黃河輸沙量影響的甄別[D].青島：中國海洋大學，2013：37－43.
　　[20] 牛最榮，崔亮.氣候變化和人類活動對洮河輸沙量變化的影響分析[J].水利水電技術，2016，47（5）：80－83，88.
　　【責任編輯 栗銘】

nlc202212201221

轉載注明來源:http://www.hailuomaifang.com/1/view-15443539.htm