摘 要:若爾蓋濕地為黃河上游主要產(chǎn)流區(qū),為研究若爾蓋濕地水文情勢的變化,選取黃河干流吉邁站、瑪曲站和主要支流白河唐克站、黑河大水站的歷史逐年、日流量資料,應(yīng)用Mann?Kendall 法、距平累積法以及Pettitt 檢驗法分析4 站突變特征。采用水文指標(biāo)變動范圍法(IHA?RVA)和水文改變度法綜合評估各站對比時段和評價時段水文指標(biāo)改變程度。研究結(jié)果表明:4 站年徑流量無顯著突變趨勢;與對比時段相比,評價時段4 站整體水文改變度接近中度改變,其中唐克站達(dá)到中度改變;32 個水文改變指標(biāo)中,極端大流量事件均以減少為主,高、低流量脈沖頻率均以下降為主,改變度最大的是流量改變率及逆轉(zhuǎn)次數(shù),瑪曲、吉邁、大水站流量增大率和減小率均下降,逆轉(zhuǎn)次數(shù)減少改變度除唐克站外均達(dá)到中度甚至高度改變,這將給河流生物群落帶來一定的影響,也在一定程度上反映了濕地退化現(xiàn)象。
關(guān)鍵詞:水文情勢;IHA?RVA 法;水文改變度;若爾蓋濕地;黃河
中圖分類號:TV62;TV882.1 文獻(xiàn)標(biāo)志碼:Adoi:10.3969/ j.issn.1000-1379.2023.04.006
引用格式:許明一,蘇柳.基于IHA?RVA 法的若爾蓋濕地水文情勢研究[J].人民黃河,2023,45(4):35-39,101.
1 研究背景
黃河上游是由一系列河流、星羅棋布的大小湖泊、沼澤和草原以及位于高山之巔的冰川雪原和埋藏于地下的潛水組成的一個完整而復(fù)雜的水資源系統(tǒng)。瑪多以下進(jìn)入峽谷段,徑流補(bǔ)給以降水為主,其中吉邁至瑪曲段為多雨區(qū),至索宗寺和瑪曲處黃河與白河、黑河匯合,水量顯著增加。阿壩州黃河流域位于此河段,其多年平均(1956—2016 年)地表水資源量41.41 億m3,占黃河流域多年平均(1956—2016 年) 地表水資源量583.6 億m3的7.1%。
若爾蓋濕地位于阿壩州若爾蓋縣境內(nèi)青藏高原東北部的黃河上游,是世界上面積最大、保存最完好的高原高寒泥炭沼澤濕地,是黃河最重要的水源涵養(yǎng)區(qū)[1-2] 。濕地對調(diào)節(jié)氣候、改善生態(tài)環(huán)境和保護(hù)生物多樣性都有著極其重要的作用,被生態(tài)學(xué)家稱為“地球之腎”。同時,若爾蓋濕地是我國生物多樣性研究熱點(diǎn)地區(qū)之一。20 世紀(jì)60 年代以來,受人類活動和氣候變化的影響,若爾蓋濕地面臨一系列問題,比如濕地萎縮、功能退化、草甸含水量減小、水土流失、濕地荒漠化、生物多樣性下降等[3-4] 。
國內(nèi)外學(xué)者非常關(guān)注人類活動對河流水文情勢的影響[5-7] 。20 世紀(jì)80 年代,美國學(xué)者提出水資源研究不僅要考慮化學(xué)指標(biāo),更應(yīng)注重生態(tài)指標(biāo),90 年代研究成果增加,并發(fā)布了恢復(fù)生態(tài)完整性等多項指導(dǎo)原則,其他國家、地區(qū)、組織也出臺相關(guān)技術(shù)導(dǎo)則。我國在2015 年后出臺各種有關(guān)河湖生態(tài)保護(hù)的規(guī)范、標(biāo)準(zhǔn)、導(dǎo)則等,其中在水文情勢方面明確了各種水文要素時空變化的態(tài)勢和趨勢,常用流量、頻率、發(fā)生時機(jī)、延續(xù)時間和變化率等表示??紤]對河流生態(tài)系統(tǒng)有重要意義的這5 種水文情勢,國外學(xué)者建立了評估生態(tài)水文變化過程的IHA 指標(biāo)體系及對河流生態(tài)健康評價的RVA 方法。運(yùn)用IHA?RVA 法研究河湖生態(tài)流量,能夠有效評估河流生態(tài)水文情勢的變化特性,目前常用該方法評價河流水文情勢[8-9] 。
本文選取若爾蓋濕地黃河干流及主要支流斷面,通過Mann?Kendall 法、距平累積法、Pettitt 檢驗法3 種方法進(jìn)行突變點(diǎn)檢驗,基于IHA?RVA 法評價水文情勢指標(biāo)及變化特性,分析討論水文情勢變化原因,以期為若爾蓋濕地生態(tài)保護(hù)管理提供科學(xué)依據(jù)。
2 研究數(shù)據(jù)及方法
2.1 數(shù)據(jù)資料
選取若爾蓋濕地黃河干流瑪曲水文站及其上游吉邁水文站徑流系列,以及若爾蓋濕地重要黃河支流白河唐克水文站、黑河大水水文站徑流系列進(jìn)行分析。采用IHA?RVA 法分析2000—2019 年水文情勢相對于對比時段的變化?,斍?、吉邁水文站將建站至20世紀(jì)90 年代作為對比時段,即瑪曲站對比時段為1960—1989 年,吉邁站對比時段為1959—1989 年;大水、唐克水文站將20 世紀(jì)80 年代建站至1999 年作為對比時段,即大水站對比時段為1985—1999 年,唐克站對比時段為1981—1999 年。數(shù)據(jù)采用水文站實測日徑流量數(shù)據(jù)。
2.2 研究方法
2.2.1 突變檢驗方法
本研究選取若爾蓋濕地相關(guān)的水文站作為代表水文站,基于各水文站建站至2019 年的歷年年均流量資料,運(yùn)用Mann?Kendall 法[10](以下簡稱M-K 法)、距平累積法[11-12] 和Pettitt 檢驗法[13] 對年均流量進(jìn)行趨勢性和突變性分析。其中M-K 法、距平累積法比較常見,本文不作詳細(xì)介紹。
Pettitt 檢驗是一種與M-K 法相似的非參數(shù)檢驗方法,該方法不適用于突變點(diǎn)較多的情況,直接利用秩序列檢驗突變點(diǎn)[14] ,因而在生態(tài)水文學(xué)領(lǐng)域廣泛應(yīng)用。Pettitt 檢驗的原理是對于包含N 個樣本的時間序列數(shù)據(jù)X,計算其統(tǒng)計量U。
把Ut 中絕對值最大的點(diǎn)記為Kt ,設(shè)該點(diǎn)為最顯著的突變點(diǎn),并計算Kt所對應(yīng)的統(tǒng)計量P,如果P 小于給定的顯著水平,那么存在統(tǒng)計顯著的突變點(diǎn)。
2.2.2 IHA?RVA 方法
為定量分析河流自然水文情勢變化特征,Richter等[15] 建立了水文變化指標(biāo)體系IHA,其包括流量特征、頻率特征、時間特征、高低脈沖流量及其發(fā)生時機(jī)特征、河流水文流量日變化幅度特征,共5 組33 個評價指標(biāo),本次研究無流量為0 的情況,故有32 個評價指標(biāo)(見表1)[16] 。
變化范圍評價法(Range of Variability Approach,RVA)可分析不同時段IHA 指標(biāo)的變化程度即水文改變度,根據(jù)評價時段落在RVA 目標(biāo)區(qū)間內(nèi)的年數(shù)占預(yù)期年數(shù)的比例評價水文情勢變化程度[17] 。本文將對比時段各IHA 指標(biāo)值的25%至75%作為RVA 目標(biāo)區(qū)間。水文改變度定義為[16-20]
式中:Di為第i 個IHA 指標(biāo)的水文改變度;No為評價時段落到RVA 目標(biāo)區(qū)間內(nèi)的年數(shù);Ne 為評價時段落到RVA目標(biāo)區(qū)間內(nèi)的預(yù)期年數(shù),本文預(yù)期年數(shù)為評價時段總年數(shù)的一半,為10;i 為IHA 指標(biāo)編號,i =1,2,…,32。
Richter 建議,當(dāng)Di 落在[0%,33%)范圍內(nèi)時,水文改變度為無改變或低度改變;在[33%,67%)范圍內(nèi)時,為中度改變;在[67%,100%]范圍內(nèi)時,為高度改變[16] 。基于IHA?RVA 法,提出綜合指標(biāo)改變度,基于每個指標(biāo)的變異程度Di ,變異后的水文變異程度D 計算公式如下[16-17,21-22] :
式中:D 為流域整體的水文改變度。
3 結(jié) 果
3.1 突變性檢驗結(jié)果
采用M-K 法和Pettitt 檢驗法進(jìn)行趨勢檢驗,結(jié)果見表2,4 個站檢驗結(jié)果均為無顯著性改變,其中Pettitt 檢驗的K 值范圍均達(dá)不到相應(yīng)的顯著水平。
4 個站的距平累積變化過程見圖1?,斍?、吉邁站1990 年之后有來水減小的趨勢,吉邁站2004 年達(dá)到最低點(diǎn);唐克站、大水站2000 年之后有來水減小趨勢,2010 年來水開始增大。由此說明IHA?RVA 法中評價年份和對比年份選取較為合理。
3.2 IHA?RVA 評價結(jié)果
3.2.1 整體情況
IHA?RVA 評價結(jié)果見表3。與20 世紀(jì)90 年代前相比,2000—2019 年瑪曲、吉邁斷面河流整體屬于低度改變;與2000 年前相比,2000—2019 年大水?dāng)嗝婧恿髡w屬于低度改變,唐克斷面河流整體達(dá)到中度改變。對于32 個IHA 評價指標(biāo),4 個斷面均以低度改變?yōu)橹?,瑪曲斷面中度改變所占比例較大。唐克斷面河流整體屬中度改變,主要原因是基流量的增大。
3.2.2 月均流量
分析2000—2019 年4 個斷面各組IHA 指標(biāo)值相對于對比時段指標(biāo)值的偏離度。對于月均流量,2000—2019年瑪曲斷面1—12 月月均流量僅5 月、12月少量增大,其余均減小,改變度7 月份為40%,其他均在0%~20%之間,減小(增大)幅度9 月份最大為27.2%,其他在7.2% ~16.3%之間;吉邁斷面1—4 月、6—7 月、11—12 月月均流量增大,其余月份減小,改變度9 月份為50%、5 月份為40%,其他均在0% ~20%之間,增大(減?。?幅度1—4 月及12 月份平均為32.3%,其他在0.6%~19.9%之間;唐克斷面月均流量僅2—3 月、11 月增大,其余減小,改變度12 月份為40%,其他均在10%~30%之間,減?。ㄔ龃螅┓龋?月份為36.6%,2 月、3 月均在24% 以上,其他月份在1.05%~17.9%之間;大水?dāng)嗝嬖戮髁績H1—2 月、10月、12 月增大,其余減小,改變度3 月為70%,2 月、4月為40%,其他月份在0%~30%之間,減?。ㄔ龃螅┓茸兓草^大,其中1—2 月、4 月、6 月、9 月在30%以上,其他月份在0.7%~23.8%之間。
總體來說,汛期平均流量均不同程度減小,枯水期11 月至次年2 月則不同程度增大。各水文站流量序列整體水文改變度分級見表4。
3.2.3 極端流量及持續(xù)時間
最大1、3、7、30、90 d 平均流量,2000—2019 年4個水文站均小于對比時段,瑪曲站平均減?。保罚叮?,改變度平均為14%;吉邁站平均減?。保保梗?,改變度平均為16%;唐克站平均減?。保梗矗?,改變度平均為24%;大水站平均減小10%,改變度平均為16%。反映出年最大流量減小,說明大流量事件在一定范圍內(nèi)處于收緊狀態(tài)。
最?。?、3、7、30、90 d 平均流量,2000—2019 年瑪曲斷面小于對比時段,平均減?。保矗担?,改變度平均為42%,屬中度改變;吉邁斷面大于對比時段,平均增大37%,改變度平均為16%;唐克斷面大于對比時段,平均增大29%,屬中高度改變;大水?dāng)嗝娲笥趯Ρ葧r段,平均增大15.9%,屬中度改變。
4 個站基流變化情況:瑪曲減?。埃担?,改變度20%;吉邁增大26.5%,改變度高達(dá)90%;唐克增大幅度較大,為76.8%,改變度高達(dá)100%;大水增大59%,改變度為50%。
3.2.4 極端流量發(fā)生時間
瑪曲站年最大流量發(fā)生時間基本不變,最小流量推遲較多,由12—1 月推遲至1—2 月;吉邁站年最大、最小流量發(fā)生時間均基本不變;唐克站年最大流量發(fā)生時間基本不變,最小流量提前較多,由2—3 月提前至1—2 月;大水站年最大流量發(fā)生時間基本不變,最小流量不規(guī)律變化;改變度方面,大水站最小流量發(fā)生時間為中度改變,唐克站最小流量發(fā)生時間為高度改變,其余均為低度改變。
3.2.5 高、低流量脈沖頻率及歷時
瑪曲站高、低流量脈沖發(fā)生次數(shù)均減少,持續(xù)時間均延長,其中低流量脈沖發(fā)生次數(shù)及持續(xù)時間均為中度改變;吉邁站高、低流量脈沖發(fā)生次數(shù)均減少,持續(xù)時間高流量脈沖略有延長、低流量脈沖略微縮短,其中低流量脈沖發(fā)生次數(shù)為中度改變;唐克站高、低流量脈沖發(fā)生次數(shù)均減少,持續(xù)時間高流量脈沖縮短、低流量脈沖延長,均為低度改變;大水站高流量脈沖發(fā)生次數(shù)減少、持續(xù)時間延長,低流量脈沖發(fā)生次數(shù)增加、持續(xù)時間縮短,均為低度改變。
3.2.6 流量改變率及逆轉(zhuǎn)次數(shù)
瑪曲、吉邁、大水站流量增大率和減小率均下降,瑪曲站逆轉(zhuǎn)次數(shù)由87 次減少至68 次,改變度達(dá)50%;吉邁站逆轉(zhuǎn)次數(shù)由130 次減少至120 次,改變度高達(dá)90%;大水站逆轉(zhuǎn)次數(shù)由59 次減少至50 次,改變度達(dá)50%。唐克站流量增大率和減小率變化不大,逆轉(zhuǎn)次數(shù)由69 次減少至51 次,改變度為20%。圖2 為水文站對比時段與評價時段逆轉(zhuǎn)次數(shù)的變化。河流流量逆轉(zhuǎn)次數(shù)多少關(guān)系著河流的生態(tài)水文規(guī)律,河流生態(tài)系統(tǒng)穩(wěn)定將受到流量變化及頻率的影響。流量逆轉(zhuǎn)次數(shù)的改變,將對水生動植物的生物量和多樣性產(chǎn)生影響,危害水生動植物的生存環(huán)境[23] 。
4 討 論
河流水文情勢受到氣候變化、人類活動等多因素的復(fù)合影響?,斍鷶嗝嬉陨系狞S河干流及黑河、白河支流上游沒有控制性水利樞紐工程,基本處于天然狀態(tài)。流量過程變異程度可能與氣候、人類活動影響有關(guān),比如全球變暖、區(qū)域蒸發(fā)量增大、過度放牧、開挖溝渠等,濕地生態(tài)環(huán)境明顯惡化[24-26] 。在黃河源區(qū)的研究結(jié)果為氣溫升高,水面蒸發(fā)量增大,降水增加趨勢并不顯著,這可能是導(dǎo)致該區(qū)域徑流量減少的主要原因[27] 。人工開渠直接疏干成片沼澤濕地是沼澤濕地退化的重要原因[25] 。人類活動使鼠蟲天敵銳減,缺乏食物鏈約束的鼠蟲害肆虐,蠶食牧草,縱橫交錯的地下洞穴使土壤水分流失、牧草枯萎死亡,從而形成大面積的水土流失。
4 個站2000—2019 年多年平均天然徑流量比長系列徑流量有所減少,4 個站平均減少不到5%。在前人分析中,若爾蓋濕地徑流量呈下降趨勢,兩條主要支流白河以及黑河流域徑流量下降趨勢顯著,且各站大流量事件減少,瑪曲站小流量事件減少,其他站小流量事件增加,說明大流量減小,小流量可能會被消耗,在一定程度上引起濕地水位下降,沼澤溯源下切速率增大,加速了濕地萎縮與退化,也可能影響黃河上游的水資源量短缺情勢[28-29] 。
在后續(xù)研究中,需深入研究河流水文情勢與水生動植物之間動態(tài)關(guān)系,結(jié)合生態(tài)監(jiān)測資料,定量分析水文情勢改變對高原水生生物的影響程度[5,30] ,還需研究水文變化與濕地儲水的關(guān)系等,以期為若爾蓋濕地生態(tài)保護(hù)和恢復(fù)提供參考。
5 結(jié) 論
采用IHA?RVA 法分析了黃河干流吉邁、瑪曲,支流黑河大水、白河唐克4 個水文站長系列水文情勢演變過程,得到以下結(jié)論。
(1)對4 個斷面逐年的流量過程進(jìn)行突變性檢驗,都沒有通過顯著性檢驗,說明流量過程變異程度的變化均是漸進(jìn)的,沒有顯著的突變點(diǎn)。
(2)與對比時段相比,評價時段河流整體改變度為低度或中度,唐克站整體改變度較大;5—10 月月均流量以減小為主,其他月份除瑪曲站外均以增大為主;極端大流量事件均以減少為主,極端小流量事件的流量除瑪曲站外均以增大為主。高、低流量脈沖頻率均以下降為主。
(3)改變度最大的是流量改變率及逆轉(zhuǎn)次數(shù),瑪曲、吉邁、大水站流量增大率和減小率均下降,逆轉(zhuǎn)次數(shù)減少改變度除唐克站外均已達(dá)到中度甚至高度改變。
參考文獻(xiàn):
[1] 趙娜娜,王賀年,張貝貝,等.若爾蓋濕地流域徑流變化及其對氣候變化的響應(yīng)[J].水資源保護(hù),2019,35(5):40-47.
[2] XIANG S,GUO R Q,WU N,et al.Current Status and FutureProspects of Zoige Marsh in Eastern Qinghai?Tibet Plateau[J].Ecological Engineering,2009,35(4):553-562.
[3] 李志威,王兆印,張晨笛,等.若爾蓋沼澤濕地的萎縮機(jī)制[J].水科學(xué)進(jìn)展,2014,25(2):172-180.
[4] 閔泓翔.若爾蓋高原濕地退化現(xiàn)狀、成因及恢復(fù)對策研究[D].成都:四川農(nóng)業(yè)大學(xué),2012:19-21.
[5] 卓志宇,錢湛,姜恒,等.基于IHA?RVA 法的湘江流域水文情勢評價[J].中國農(nóng)村水利水電,2020(8):188-192.
[6] 郭文獻(xiàn),夏自強(qiáng),王乾.丹江口水庫對漢江中下游水文情勢的影響[J].河海大學(xué)學(xué)報(自然科學(xué)版),2008,36(6):733-737.
[7] YANG Z F,YAN Y,LIU Q.Assessment of the Flow RegimeAlterations in the Lower Yellow River,China[J].EcologicalInformatics,2012,10:56-64.
[8] 陳棟為,陳曉宏,李翀,等.基于RVA 法的水利工程對河流水文情勢改變的累積效應(yīng)研究:以東江流域為例[J].水文,2011,31(2):54-57.
[9] 班璇,姜劉志,曾小輝,等.三峽水庫蓄水后長江中游水沙時空變化的定量評估[J].水科學(xué)進(jìn)展,2014,25(5):650-657.
[10] 田小靖,趙廣舉,穆興民,等.水文序列突變點(diǎn)識別方法比較研究[J].泥沙研究,2019,44(2):33-40.
[11] 張舒羽,趙廣舉,穆興民,等.黃土高原典型流域河川徑流水文情勢變化及其歸因分析[J]. 水土保持通報,2021,41(4):1-8.
[12] 張建云,章四龍,王金星,等.近50 年來中國六大流域年際徑流變化趨勢研究[J].水科學(xué)進(jìn)展,2007,18(2):230-234.
[13] 高鵬.黃河中游水沙變化及其對人類活動的響應(yīng)[D].北京:中國科學(xué)院研究生院(教育部水土保持與生態(tài)環(huán)境研究中心),2010:21-25.
[14] 姚天次,盧宏瑋,于慶,等.近50 年來青藏高原及其周邊地區(qū)潛在蒸散發(fā)變化特征及其突變檢驗[J].地球科學(xué)進(jìn)展,2020,35(5):534-546.
[15] RICHTER B D,BAUMGARTNER J V,POWELL J,et al.AMethod for Assessing Hydrologic Alteration Within Ecosys?tems[J].Conservation Biology,1996,10(4):1163-1174.
[16] 尚文繡,暒齙隊蝊暒許明一,尚奕,等.龍羊峽水庫調(diào)度對徑流的影響及蓄補(bǔ)水規(guī)律[J].南水北調(diào)與水利科技(中英文),2022,20(3):451-458.
[17] 尚文繡,許明一,尚奕.黃河干流水文情勢長期演變趨勢研究[C] // 北京沃特咨詢有限公司.2021 第九屆中國水生態(tài)大會論文集.北京:北京沃特咨詢有限公司,2021:490-497.
[18] RICHTER B D,BAUMGARTNER J V,BRAUN D P,et al.A Spatial Assessment of Hydrologic Alteration Within aRiver Network[J].Regulated Rivers:Research & Manage?ment,1998,14(4):329-340.
[19] 李興拼,黃國如,江濤.RVA 法評估楓樹壩水庫對徑流的影響[J].水電能源科學(xué),2009,27(3):18-21.
[20] 李舒,齊青松.基于IHA?RVA 法的人類活動對窟野河流域水文情勢影響評估[J].水電能源科學(xué),2021,39(4):24-26,63.
[21] 羅志勇,楊炳良,管曉祥,等.不同水文模型在金溪流域的模擬與適應(yīng)性研究[J].華北水利水電大學(xué)學(xué)報(自然科學(xué)版),2018,39(3):6-12.
[22] 李紫妍,劉登峰,黃強(qiáng),等.基于多種水文學(xué)方法的漢江子午河生態(tài)流量研究[J].華北水利水電大學(xué)學(xué)報(自然科學(xué)版),2017,38(1):8-12.
[23] 曾金鳳,劉祖文,劉友存,等.基于IHA?RVA 法的東江源區(qū)生態(tài)水文情勢變化分析[J].水土保持通報,2021,41(6):157-164.
[24] 雍國瑋,石承蒼,邱鵬飛.川西北高原若爾蓋草地沙化及濕地萎縮動態(tài)遙感監(jiān)測[J].山地學(xué)報,2003,21(6):758-762.
[25] 劉佳,陳超,秦寧生,等.青藏高原若爾蓋生態(tài)區(qū)水資源對氣候變化的響應(yīng)[J].冰川凍土,2016,38(2):498-508.
[26] 李志威.三江源河床演變與濕地退化機(jī)制研究[D].北京:清華大學(xué),2013: 159-165.
[27] 張永勇,張士鋒,翟曉燕,等.三江源區(qū)徑流演變及其對氣候變化的響應(yīng)[J].地理學(xué)報,2012,67(1):71-82.
[28] 魯瀚友,李志威,胡旭躍,等.若爾蓋高原徑流量變化與儲水量計算[J].水資源與水工程學(xué)報,2019,30(6):12-19.
[29] 周澤江,覃光華,于春平,等.若爾蓋濕地黑河徑流分析及預(yù)測[J].水電與新能源,2013,27(3):18-22.
[30] ACREMAN M,ARTHINGTON A H,COLLOFF M J,et al.Environmental Flows for Natural,Hybrid,and Novel RiverineEcosystems in a Changing World[J].Frontiers in Ecologyand the Environment,2014,12(8):466-473.
【責(zé)任編輯 張 帥】