氣候變化和人類活動對灞河流域輸沙量影響的定量評估

夏 偉周維博何慶龍李文溢

(1.長安大學環境科學與工程學院,陜西西安 710054;

2.長安大學旱區地下水文與生態效應教育部重點實驗室,陜西西安 710054)

我國很多學者對我國各大流域的水沙特性及其影響因素也作了大量研究。劉成等[6]對我國主要入海河流水沙變化的分析表明,我國最主要的10條河流下游站1997—2006年輸沙量的平均值低于多年平均值的一半。姚文藝等[7]綜合利用水文法、水保法和數學模型法等多種方法對黃河流域1997—2006年的水沙變化情勢進行了研究,結果表明黃河實測輸沙量較1970年以前年均減小11郾80億t,其中人類活動的減沙作用與降雨的影響基本相當。史紅玲等[8]研究了淮河流域水沙變化趨勢及其成因,表明淮河流域來沙量呈現顯著減小的趨勢,發生突變的年份在1973—1975年之間。

目前,已有學者對灞河流域水沙特性進行了研究。趙建春等[9]對灞河流域水沙特性的研究表明,1985年后的來沙量基本保持在較小的水平;魏炳乾等[10]認為灞河流域的水沙量自20世紀70年代以后呈減小趨勢,其中20世紀90年代減小十分顯著,人類活動和氣候變化是灞河水沙量減小的主要因素。為研究各因素對灞河流域水沙量減小的影響程度,采用累積量斜率變化率比較法(SCRCQ)[11鄄15]對氣候變化和人類活動在灞河流域輸沙量變化中的相對貢獻率進行定量分析評估。

1 研究區概況

灞河系渭河的一級支流,發源于秦嶺北麓藍田縣灞源鄉麻家坡以北,位于東經 109毅00憶~109毅47憶,北緯 33毅50憶~34毅27憶,流經西安市灞橋區、未央區,于高陵區匯入渭河,全長104 km,流域總面積為2581km2,主要支流有浐河、清河、道溝峪等。灞河上游河道比降較大,平均為0郾9%;中游為平原彎曲型河道,平均比降為0郾235%;下游河道比較緩,平均比降為0郾158%。灞河流域的基本情況如圖1所示。

圖1 灞河流域

2 基礎數據與研究方法

2.1 基礎數據

馬渡王水文站設立于1952年6月,位于浐河匯入灞河的入匯口附近,距離渭河口30km,是灞河下游干流控制站,控制流域面積1601km2。馬渡王水文站觀測資料系統全面、時間序列較長,因此選用該水文站1960—2012年的實測水沙數據作為基礎數據。

2.2 研究方法

2.2.1 小波分析法

小波分析[16]屬于時頻分析的一種,它能夠清晰地揭示出隱藏在時間序列中不同的變化周期,充分反映系統在不同時間尺度上的變化趨勢,且能對系統在未來的變化趨勢作定性分析。其基本思想是:用一簇小波函數系來表示或逼近某一信號或函數,因此,小波函數是小波分析的關鍵,它是指具有震蕩性、能夠迅速衰減到零的一類函數,即小波函數追(t)沂L2(R)且滿足:

式中:追a,b(t)為子小波;a為尺度因子,反映小波的周期長度;b為平移因子,反映時間上的平移。

2.2.2 累積距平法

累積距平法[17鄄18]是一種由曲線直觀判斷離散數據點變化趨勢的統計方法,通過先計算每年輸沙量的距平值,然后按照時間序列逐年累加,從而得到累積距平值隨時間的變化趨勢。如果輸沙量累積距平曲線上升,則輸沙量為逐年增加的趨勢;反之,輸沙量為逐年減小的趨勢,趨勢的轉折點即為突變點。

2.2.3 累積量斜率變化率比較法

由王隨繼等[11]提出的SCRCQ方法以客觀年份作為自變量,以各因子的累積量作為因變量,因累積量在一定程度上消除了實測數據年際變化的影響,使年份與累積量之間具有非常高的相關性,為進一步進行定量分析創造了條件,并避免了采用常用的多元回歸法進行定量分析時的人為因素的影響。其基本公式為式中追(t)為基小波函數。選用Morlet小波函數作為基小波函數,它可通過對基小波函數在尺度的伸縮和時間軸上的平移構成一族函數系:

式中:RS為累積輸沙量斜率變化率;KSa、KSb分別為累積輸沙量-年份線性關系式在拐點前后兩個時期的斜率,萬t/a;RP為累積降水量斜率變化率;KPb、KPa分別為累積降水量年份線性關系式在拐點前后兩個時期的斜率,mm/a;RE為累積蒸散發量斜率變化率;KEb、KEa分別為累積蒸散發量年份線性關系式在拐點前后兩個時期的斜率,mm/a。

根據SCRCQ方法,降水、蒸散發、人類活動對輸沙量的貢獻率CP、CE、CH分別為

式中:CT為氣溫對輸沙量的影響;CC為氣象條件的變化引起的天然植被的變化對輸沙量的影響。

由于氣溫主要是通過影響蒸散發進而影響輸沙量,所以不考慮氣溫對輸沙量的影響。氣候條件的變化引起的天然植被的變化會影響產沙量,進而影響輸沙量,而天然植被的變化主要是受人類活動的影響,可不考慮氣象條件的變化引起的天然植被的變化對輸沙量的影響。根據魏炳乾等[8]對灞河流域水沙特性的研究,蒸散發對輸沙量的影響很小,式(8)可以簡化為

3 研究結果與分析

3.1 年輸沙量、年降水量及汛期降水量變化特征

1960—2012年灞河流域年輸沙量、年降水量及汛期降水量的年際變化過程如圖2所示。由圖2可知,1960—2012年多年平均輸沙量為182郾86萬t,年輸沙量呈減小的趨勢,減小的速率為6郾150萬t/a;1960—2012年多年平均降水量為616郾31 mm,年降水量呈減小的趨勢,減小的速率為1郾761 mm/a;1960—2012年多年平均汛期(7—9月)降水量為292郾73 mm,占多年平均降水量的47郾5%,汛期降水量呈減小的趨勢,減小的速率為0郾221 2 mm/a。灞河流域年輸沙量小波分析時頻分布如圖3所示,可以看出,年輸沙量存在3~7 a、13~18 a和40~45 a 3類尺度的周期,其中第一主周期為15 a左右,在13~18a尺度上存在豐枯交替準3次震蕩,且該尺度的周期變化在整個分析階段非常穩定,具有全域性。

圖2 1960—2012年灞河流域年輸沙量、年降水量及汛期降水量的年際變化

圖3 灞河流域年輸沙量小波分析時頻分布

3.2 年輸沙量突變分析及階段劃分

采用累積距平法計算1960—2012年灞河流域年輸沙量突變年份,得灞河流域年輸沙量累積距平曲線如圖4所示。由圖4(a)可以看出,年輸沙量累積距平曲線在1960—1987年表現為上升趨勢,而在1988—2012年表現為下降趨勢,在1988年發生明顯轉折,可確定1988年為年輸沙量的突變年份。由圖4(b)可以看出,年輸沙量累積距平曲線表現出先下降后上升的趨勢,其突變年份為2002年。

圖4 灞河流域年輸沙量累積距平曲線

通過流域調查發現,灞河流域在20世紀80年代以前,受人類活動的影響相對較小,灞河流域年輸沙量變化主要是受氣候變化的影響,因此,可以將1960—1987年這一階段看作基準期;在20世紀90年代以來,人類活動對灞河流域年輸沙量變化的影響逐漸加大,灞河流域年輸沙量除了受氣候變化的影響,在一定程度上還受到人類活動的影響。

3.3 突變年份分割時期的年輸沙量及年降水量與年份之間的關系

根據灞河流域年輸沙量的兩個突變年份,將研究時間分為 1960—1987年、1988—2001年和2002—2012年3個時期,分別記為A、B和C。灞河流域累積輸沙量,累積降水量與年份之間的關系如圖5所示,圖中還給出了擬合關系式,各擬合關系式的相關系數R均大于0郾969,通過了顯著性水平為0郾01的Pearson顯著性檢驗。

圖5 灞河流域累積輸沙量、累積降水量與年份之間的關系

3.4 年降水量和人類活動對年輸沙量變化的貢獻率

根據SCRCQ方法,在基準期內,輸沙量主要受自然因子(降水)的影響,受人類活動的影響較小;在變異期內,輸沙量主要受自然因子(降水)和人類活動二者的綜合影響。

由圖5(a)擬合關系式可得,累積輸沙量-年份線性關系式在A、B和C時期的斜率分別為228萬t/a,76郾53萬t/a和101郾3萬t/a;B時期與A時期相比,斜率減小了151郾47萬t/a,減小率為66郾43%;C時期與A時期相比,斜率減小了126郾7萬t/a,減小率為55郾57%;C時期與 B時期相比,斜率增大了24郾77 萬 t/a,增大率為32郾37%。

由圖5(b)擬合關系式可得累積降水量-年份線性關系式在A、B和C時期的斜率分別為638郾2 mm/a,513郾9 mm/a和589郾8 mm/a;B時期與A時期相比,斜率減小124郾3 mm/a,減小率為19郾48%;C時期與A時期相比,斜率減小48郾4 mm/a,減小率為7郾58%;C時期與B時期相比,斜率增大75郾9 mm/a,增大率為 14郾77%。

根據式(6)~(9),可以計算得到B時期與A時期相比,降水量的減小對輸沙量減小的貢獻率為12郾22%,人類活動對輸沙量減小的貢獻率為87郾78%;C時期與A時期相比時,降水量的減小對輸沙量減小的貢獻率為6郾56%,人類活動對輸沙量減小的貢獻率為93郾44%;C時期與B時期相比,降水量的增大對輸沙量增大的貢獻率為0郾53%,人類活動對輸沙量增大的貢獻率為99郾47%,降水量對輸沙量的影響很小,可以忽略不計,即輸沙量增大完全受人類活動的影響。

3.5 年輸沙量減小的原因分析

灞河流域年輸沙量減小的原因有自然因素(降水量減小和降水強度減小)和非自然因素(水庫攔沙、水土保持、人工采砂和土地利用類型的變化)[18]。

a.自然因素。由于灞河流域降水量的減小,會減小雨水對地面的沖刷,減少土壤的流失,從而使得灞河流域年輸沙量在一定程度上有所減小;同時,由于灞河流域汛期降水強度的減小,也在一定程度上減少產沙量,從而使灞河流域年輸沙量有所減小。

b.非自然因素。自20世紀80年代以來,人類活動對灞河流域年輸沙量減小的影響越來越大,主要表現在:淤在灞河流域建了很多水庫(如紅旗水庫),流域內侵蝕的泥沙大量淤積在水庫內,使得下游年輸沙量減小[19];于灞河流域于1997年開始實施攔河造湖工程,在灞河上規劃建造10座橡膠壩,形成水面面積12萬m2,另外在灞河生態區上游布置大片生態水系涵養綠地,沿河道兩岸布置帶狀公園和防護綠地,在一定程度上遏制了水土流失,使灞河流域年輸沙量亦隨之減小[19];盂隨著西安市經濟的快速發展,以及大規模的基礎設施建設,在灞河流域存在大量的引水采砂作業,這也是灞河流域年輸沙量減小的一個重要原因;榆隨著大西安城市建設,灞河流域內修建了大量住宅和道路,土壤裸露面積減小,這使土壤流失有所減小,從而使灞河流域年輸沙量有所減小。

4 結 論

a.自1960年以來,灞河流域年輸沙量一直呈減小趨勢,減小速率為6郾150萬t/a;年輸沙量存在3~7 a、13~18 a和40~45 a 3類尺度的周期,其中第一主周期為15 a左右;灞河流域年降水量也一直呈減小趨勢,減小速率為1郾761 mm/a。

b.采用累積距平方法確定了灞河流域年輸沙量在1988年和2002年前后發生突變。3個時期累積輸沙量和累積降水量與年份之間的關系較好,相關系數均大于 0郾969,且通過了 0郾01的顯著性檢驗。

c.采用累積量斜率變化率比較方法分析了降水量和人類活動在灞河流域輸沙量減小中的相對貢獻率,B時期與A時期相比時,降水量和人類活動在灞河流域輸沙量減小中的相對貢獻率分別為12郾22%和87郾78%;C時期與A時期相比時,降水量和人類活動在灞河流域輸沙量減小中的相對貢獻率分別為6郾56%和93郾44%;C時期與B時期相比時,降水量和人類活動在灞河流域輸沙量增加中的相對貢獻率分別為0郾53%和99郾47%。由此可見,人類活動是影響灞河流域輸沙量變化的主要因素。

[1]ZHANG Q,XU C Y,BECKER S,et al.Sediment and runoff changes in the Yangtze River Basin during last 50 years[J].Journal of Hydrology,2006,331(3/4):511鄄523.

[2]葉許春,劉健,李相虎,等.氣候和人類活動對贛江徑流變化的作用分析[J].河海大學學報(自然科學版),2013,41(3):196鄄203.(YE Xuchun,LIU Jian,LI Xianghu,et al.Effects of climate variability and human activities on runoff variation of Ganjiang River Basin[J].Journal of Hohai University(Natural Sciences),2013,41(3):196鄄203.(in Chinese))

[3]吳麗,張愛靜.氣候變化和人類活動對大凌河上游流域徑流的影響[J].水利水電科技進展,2016,36(2):10鄄15.(WU Li,ZHANG Aijing.Effects of climate change and human activities on runoff in upper Daling River Basin[J].Advances in Science and Technology of Water Resources,2016,36(2):10鄄15.(in Chinese))

[4]WALLING D E,FANG D.Recent trends in the suspended sediment loads of the world爺s rivers[J].Global&Planetary Change,2003,39(1):111鄄126.

[5]WALLING D E.Human impact on land鄄ocean sediment transfer by the world爺s rivers[J].Geomorphology,2006,79(3/4):192鄄216.

[6]劉成,王兆印,隋覺義.我國主要入海河流水沙變化分析[J].水利學報,2007,38(12):1444鄄1452.(LIU Cheng,WANG Zhaoyin,SUI Jueyi.Analysis on variation of seagoing water and sediment load in main rivers of China[J].Journal of Hydraulic Engineering,2007,38(12):1444鄄1452.(in Chinese))

[7]姚文藝,冉大川,陳江南.黃河流域近期水沙變化及其趨勢預測[J].水科學進展,2013,24(5):607鄄616.(YAO Wenyi,RAN Dachuan,CHEN Jiangnan.Recent changes in runoffand sedimentregimesand future projections in the Yellow River Basin[J].Advances in Water Science,2013,24(5):607鄄616.(in Chinese))

[8]史紅玲,胡春宏,王延貴,等.淮河流域水沙變化趨勢及其成因分析[J].水利學報,2012,43(5):571鄄579.(SHI Hongling,HU Chunhong,WANG Yangui,et al.Analysis of trends and reasons of variations in runoff and sediment load of the Huaihe River[J].Advances in Water Science,2012,43(5):571鄄579.(in Chinese))

[9]趙建春,王新宏,李九發,等.灞河流域水沙特性研究[J].人民黃河,2008,30(4):20鄄22.(ZHAO Jianchun,WANG Xinhong,LI Jiufa,et al.Research on characteristics of flow and sediment of Bahe Basin[J].Yellow River,2008,30(4):20鄄22.(in Chinese))

[10]魏炳乾,楊川,石忠科,等.灞河水沙特性變化分析[J].水資源與水工程學報,2008,19(2):38鄄40.(WEI Bingqian,YANG Chuan,SHI Zhongke,et al.Analysis on the changing characteristics of sediment and water in the Bahe River[J].Journal of Water Resources and Water Engineering,2008,19(2):38鄄40.(in Chinese))

[11]王隨繼,閆云霞,顏明,等.皇甫川流域降水和人類活動對徑流量變化的貢獻率分析:累積量斜率變化率比較方法的提出及應用[J].地理學報,2012,67(3):388鄄397.(WANG Suiji,YAN Yunxia,YAN Ming,et al.Contributions of precipitation and human activities to the runoff change of the Huangfuchuan Drainage Basin:application of comparative method of the slope changing ratio of cumulative quantity[J].Acta Geographica Sinica,2012,67(3):388鄄397.(in Chinese))

[12]WANG S,YAN M,YAN Y,et al.Contributions of climate change and human activities to the changes in runoff increment in different sections of the Yellow River[J].Quaternary International,2012,282(60):66鄄77.

[13]王隨繼,李玲,顏明.氣候和人類活動對黃河中游區間產流量變化的貢獻率[J].地理研究,2013,32(3):395鄄402.(WANG Suiji,LI Ling,YAN Ming.The contributions of climate change and human activities to the runoff yield changesin the middle Yellow River Basin[J].Geographical Research,2013,32(3):395鄄402.(in Chinese))

[14]張調風,朱西德,王永劍,等.氣候變化和人類活動對湟水河流域徑流量影響的定量評估[J].資源科學,2014,36(11):2256鄄2262.(ZHANG Tiaofeng,ZHU Xide,WANG Yongjian,et al.The impact of climate variability and human activity on runoff changes in the Huangshui River Basin[J].Resources Science,2014,36(11):2256鄄2262.(in Chinese))

[15]王彥君,王隨繼,蘇騰.降水和人類活動對松花江徑流量變化的貢獻率[J].自然資源學報,2015,30(2):304鄄314.(WANG Yanjun,WANG Suiji,SU Teng.Contributions of precipitation and human activities to runoff change in the Songhua River Basin[J].Journal of Natural Resouces,2015,30(2):304鄄314.(in Chinese))

[16]李春暉,鄭小康,楊志峰,等.黃河天然徑流量變化趨勢及其影響分析[J].北京師范大學學報(自然科學版),2009,45(1):80鄄85.(LI Chunhui,ZHENG Xiaokang,YANG Zhifeng,et al.Trends of annual natural runoff in the Yellow River Basin[J].Journal of Beijing Normal University(Natural Science),2009,45(1):80鄄85.(in Chinese))

[17]RAN L,WANG S,FAN X.Channel change at Toudaoguai Station and its responses to the operation of upstream reservoirs in the upper Yellow River[J].Journal of Geographical Sciences,2010,20(2):231鄄247.

[18]劉惠英,程冬兵,張平倉,等.三峽庫區龍河流域降雨和人類活動對水沙關系變化貢獻率分析[J].長江科學院院 報,2015,32(3):20鄄26.(LIU Huiying,CHENG Dongbing,ZHANG Pingcang,etal.Contributionsof precipitation and human activities to the runoff and sediment variation in Longhe River Watershed of Three Gorges Reservoir Area[J].Journal of Yangtze River Scientific Research Institute,2015,32(3):20鄄26.(in Chinese))

[19]李洋洋,白潔芳,周維博,等.50年來灞河流域水沙變化特征分析[J].泥沙研究,2017,42(3):20鄄25.(LI Yangyang, BAI Jiefang, ZHOU Weibo, et al.Characteristics of flow and sediment of Bahe Basin in recent 50 years[J].Journal of Sediment Research,2017,42(3):20鄄25.(in Chinese))