氣候變化和人類活動對海河流域徑流變化的影響

王 磊

(河北省衡水水文水資源勘測局，河北 衡水 053000)

0 引 言

河流徑流量的變化既受氣候變化的影響，也受人類活動的影響[1]。氣候變化導致溫度上升和降水的空間分布格局的變化，其對區域徑流量的變化具有重要影響[2]。人類活動如森林砍伐、土地利用變化、灌溉和水庫大壩建設等也影響著徑流的變化[3-5]。在過去幾十年中，水量和水質的惡化已成為流域或區域范圍內水資源管理所面臨的日益嚴重的問題[6]。

海河流域是一個缺水地區，其在過去經歷過長期干旱。近年來，山區的徑流急劇減少。大多數河流的年徑流量都有減少的趨勢，比如灤河、潮白河、漳河。目前，大多數研究都應用水文模型來分析氣候變化和人類活動對徑流變化的影響，如VIC模型和GBHM模型等[7]。雖然這些模型是該方向研究的有力工具，但研究結果具有許多不確定性，這主要是由模型結構、參數校準和換算比例的缺陷造成的[8]。此外，水文模型需要大量的輸入數據[9]。作為替代以上的方法，統計分析已被證明是對徑流響應的各種影響因素檢驗的有效選擇[10]。長期記錄的水文數據間接展示了其受氣候和土地覆蓋變化造成的影響。從長期水文數據中分析這種變化不僅可以了解匯水區的徑流變化，還可以解讀氣候變化和人類活動對其造成的影響[11]。由于河北省內的海河流域人類活動強烈[12]，因此本文的研究目的為:①分析海河流域年徑流及季度徑流的突變及趨勢變化；②分析氣候變化(降水和蒸散發量)及人類活動(水庫水量調節和土地利用變化)對徑流變化影響的貢獻率。

1 研究區概況

本文以河北省境內的海河流域為研究區，N39°44′-N42°44′，E115°33′-E119°36′。該區域總面積為44 618 km2，占海河總流域面積的17.06%。地形以低山丘陵地貌為主，高程范圍為28～2 171 m。該地區屬于半濕潤的大陸性季風氣候，年平均溫度為9.5℃。降水量在時間上分布不均勻，6～9月份的降水量占全年的80%以上[13]。年平均降水量為527 mm。從多年月降水量的平均值(圖1)可以看出，7月份的降水量最大，其次為8月份，其降水量分別為166.2和144.6 mm。

圖1 多年平均月降水量

2 數據來源

本文采用的年降水量和徑流量數據(36個雨量站及8個水文站)來自中國水文年鑒。這些數據目前尚未公布，但已打印出來供內部使用。12個氣象站數據來自國家氣象信息中心氣象數據共享服務網( http://data.cma.cn/)，1979、1990、1995和2000年4期土地利用數據來自國家地球系統科學數據共享服務平臺(http://www.geodata.cn/)。

3 研究方法

3.1 Mann-Kendall趨勢檢驗

Mann-Kendall非參數檢驗被廣泛用于檢測氣象要素時間序列的趨勢變化及是否發生顯著性變化，其計算方法如下[14]：

定義統計量S：

(1)

(2)

趨勢計算公式為：

(3)

(4)

Z值被用于判斷時間序列的變化趨勢，當Z大于0時，表明序列具有增加趨勢；當Z小于0時，表明序列具有減少趨勢。在顯著性水平α下，當|Z|≥Z1-α/2時，表明時間序列存在顯著性變化趨勢。

3.2 Mann-Kendall突變檢測

定義統計量Sk[14]：

(5)

(6)

定義統計量UF[14]：

(7)

(8)

(9)

將數據進行逆序排列，采用上述公式計算得到UB，當UF線與UB線的交點在±Uα之間時，且數據存在顯著性變化趨勢時候，說明交點為突變點。Uα為在顯著性水平α下，查正態分布表獲得的統計量值。

3.3 氣候變化和人類活動對徑流量的影響

氣候變化引起的年平均徑流量的變化為[15-17]：

ΔQclim=βΔP+γΔE0

(10)

(11)

(12)

(13)

式中：ΔQclim為由氣候變化引起的徑流變化；ΔP為突變點前后降水量變化；ΔE0為突變點前后蒸散發量變化；β、γ為徑流對降水及潛在蒸散量的敏感系數；x為干燥指數；ω為植物水分利用效率。

人類活動對徑流變化的影響為：

Qhum=ΔQ-Qclim

(14)

式中：ΔQ為突變點前后徑流量的變化；Qhum為由人類活動引起的徑流變化。

4 結果分析

4.1 年徑流的趨勢及突變分析

趨勢分析有助于了解水文和氣候變量的長期變化特點。本文采用M-K法對1956-2000年的年徑流數據進行了分析。Z統計值為2.96，在0.01顯著性水平下表現為顯著性下降趨勢。圖2給出M-K突變檢驗結果，曲線的交點出現在1979年，其年徑流量在0.05顯著水平上存在突變。根據1979年的變化點，徑流量可以分為兩部分，分別為基準期(1956-1979年)和變化期(1980-2000年)，見圖3。負的UF值表明年徑流量存在下降趨勢。在基準時期年平均徑流量為108.35 mm，而變化時期的值僅為51.63 mm，變化時期的徑流量大幅減少了56.72 mm(52.34%)。

圖2 年平均徑流量趨勢及突變點檢測結果

圖3 1956-2000年年平均徑流量變化

4.2 氣候變化和人類活動對年徑流變化的影響

表1給出基準期和變化期徑流及氣象因素的變化量。本文利用1960-2000年的氣象數據，采用Penman-Monteith模型[18-19]計算12個氣象站的潛在蒸散量。與基準時期相比，變化期年降水量減少了51.90 mm(8.65%)，而潛在蒸散量變化不大，僅為6.94 mm，僅增加了0.71%。

表1 不同時期年平均徑流量、降水及蒸散發的變化量

本文土地利用的覆蓋范圍通過海河流域的土地利用圖獲得，根據1979、1990、1995和2000年的土地利用圖，可獲得不同土地利用覆蓋范圍的比例，見表2。

表2 不同土地利用方式所占的比例

由表2可以看出，不同土地利用方式所占比例的大小分別為：森林>草地>農田>荒地>建設用地>水域，森林所占比例的平均值為41.02%，水域的均值為1.35%。根據式(10)計算可得，由氣候變化引起的徑流變化Qclim為23.19 mm，占40.89 %。而ΔQ為56.72 mm，因此由人類活動引起的徑流變化Qhum為33.53，占59.12%。

人類活動對河流徑流量的影響主要是由土地利用變化和水庫水量調節引起的。通過對比海河流域的基準流量和調整流量，可以判斷水庫流量調節的效果。調整后的流量資料從水文局獲得。比較年觀測流量和調整流量，結果表明在基準期(1956-1979年)，水庫對徑流影響較小，水庫流量減少僅為3.31 mm；而變化期影響明顯，減少流量為25.35 mm。因此，在總徑流變化ΔQ中，水庫的影響為22.04 mm，占38.86%。土地利用變化的影響為11.49 mm，占總徑流變化的20.26%。最終計算的氣候變化和人類活動對年徑流的影響見表3?？傮w來看，人類活動的影響大于氣候變化，在人類活動的影響因素中，水庫的影響大于土地利用變化，為其2倍左右。

表3 氣候變化和人類活動對徑流變化的影響

4.3 季節性徑流量的突變及趨勢變化

為了更好地分析氣候和人類活動對季節性徑流量變化的影響，根據其降水特征(圖1)，本文將年徑流量分為旱季(10月份至次年5月份)和雨季(6-9月份)。M-K突變點檢測的結果見圖4和圖5，在0.05顯著性水平下，季節性徑流量的突變點也發生在1979年。因此，季節性徑流量也可分為兩部分:基準周期(1956-1979年)和變化周期(1980-2000年)。表4給出兩個時期之間降水和徑流量的變化。由表4可以看出，兩個季節的徑流量均表現出明顯的減少；雨季減少了56.36%，旱季減少了40.43%。雨季降水量減少了12.19%，而旱季降水量增加了4.84%。

圖4 雨季徑流M-K突變檢驗結果

圖5 旱季徑流M-K突變檢驗結果

表4 基準期和變化期積極性降水量及徑流量的變化

4.4 氣候變化和人類活動對季節性徑流量的影響

為確定人類活動開始顯著影響徑流的時間，本文分析了1956-2000年的年徑流量和降水的雙累積曲線。雙累積曲線可以用來檢查降雨-徑流關系是否隨時間一致。曲線的顯著變化可表明降水以外的因素可能會對徑流變化造成影響徑流。在1956-1979年期間，本文的季度徑流量與降水之間呈一條直線變化。從1980年開始，干濕季節的徑流量均偏離了1980年以前延伸的直線(即延伸或預測線)，這表明自然產生的徑流比預測的要少。季節性累積徑流量與預測值之間的差異被稱為人類活動干擾引起的季節性流量偏差。這些徑流偏差是由于人類活動對海河流域季節性徑流的影響所造成的。計算結果表明，人類活動引起的徑流變化在雨季為-25.80 mm，在旱季為-11.95 mm。根據表5可知，徑流在雨季的整體變化為-45.69 mm，旱季為11.03 mm；氣候變化在旱季引起的徑流變化為-19.89 mm，在旱季為0.92 mm。對于人類活動對季節性流量的影響，水庫流量調節的影響也可通過實測流量和調節流量的比較來判斷，分析結果表明，水庫在雨季對徑流變化的影響為-18.04 mm，在旱季的影響為-5.24 mm。其次土地利用變化在雨季的影響為-7.76 mm，在旱季為-6.71 mm。氣候變化和人類活動對海河流域季節性流量的貢獻率見表5。結果表明，在雨季，氣候變化的影響占43.53%，人類活動占56.47%。人類活動的貢獻率中水庫調節占39.48%，而土地利用變化僅占16.99%。旱季河流流量的變化主要是由人類活動引起的，占92.85%，而氣候變化僅占7.15%。人類活動的貢獻率為92.85%，其中水庫調節占40.71%，土地利用變化占52.14%。

表5 在雨季和旱季氣候變化和人類活動對季節性徑流變化的影響

5 結 論

海河流域年徑流和季節性徑流的趨勢及突變分析結果表明，徑流有顯著的下降趨勢，其突變點出現在1979年；土地利用方式所占比例的大小分別為：森林>草地>農田>荒地>建設用地>水域，森林和水域所占比例的平均值分別為41.02%和1.35%；對于年徑流量的變化，氣候變化的貢獻率為40.89%，人類活動為59.11%，其中水庫調節占38.86%，土地利用變化占20.26%，人類活動的影響大于氣候變化，其中水庫的影響大于土地利用變化，為其2倍左右。雨季徑流變化中，氣候變化占43.53%，人類活動占56.47%，其中水庫調節占39.48%，土地利用變化占16.99%。旱季徑流變化主要是由人類活動引起的，占92.85%，而氣候變化僅占7.15%；在人類活動的貢獻中，水庫調節占40.71%，土地利用變化占52.14%?？傊?，對于海河流域，無論是年徑流還是季節性徑流，人類活動對徑流變化造成的影響均大于氣候變化。