長江源區近32年徑流變化及影響因素分析

朱延龍，陳 進，陳廣才

長江源區近32年徑流變化及影響因素分析

朱延龍1，陳 進2a，陳廣才2b

（1.青海省水文局，西寧 810001；2.長江科學院，a.院長辦公室；b.水資源利用研究所，武漢 430010）

以長江源區控制站直門達1978－2009年實測徑流量序列為基礎，運用Mann-Kendall趨勢分析法和滑動T檢驗方法，對徑流量的趨勢變化和突變進行了詳細分析，并對其影響因素降水、氣溫和下墊面條件等進行了分析。結果表明：長江源區近32年來年徑流量總體呈現增加態勢，但趨勢不顯著，年徑流序列2004年前后發生突變；長江源區降水序列在2002年前后發生突變，氣溫序列呈顯著增加趨勢，長江源區降水增加和冰雪消融是長江源區徑流增加的主要原因。

長江源區；徑流變化；氣候變化；Mann-Kendall趨勢分析法；滑動T檢驗法

1 概 述

徑流是一定時期內氣候因素、下墊面自然因素、人類活動（下墊面人為因素）等綜合作用的產物。近年來，在氣候變化和人類活動的共同作用下，徑流序列形成的物理背景發生了較大變化，直接影響了流域水資源的合理配置、開發與利用，以及河流系統的物理、化學和生物過程［1］。長江源區是青藏高原上高原濕地主要分布地區之一，也是江河源區冰川分布最集中的地區，其冰川面積占整個三江源區的89%以上，冰川融水占長江源區徑流的25%以上［2］。長江源區是中國海拔最高的天然濕地和生物多樣性分布區以及生物物種形成、演化的區域之一，具有水源涵養與調節、生物多樣性保護、長江流域生態安全保障等生態功能［3］。然而，近幾十年來，在過度放牧等人為因素和全球氣候變化的共同影響下，長江源區氣候暖干化趨勢明顯，造成了一系列生態環境問題，如冰川退縮，凍土退化，土地荒漠化，濕地與草場退化，以及水文系統發生了變化［3］。氣候變化和人類活動對長江源區水資源和生態環境的影響已成為社會日益關注的焦點和研究的熱點問題。

本文將應用Mann-Kendall趨勢分析法、滑動T檢驗法等方法，分析長江源區徑流的變化趨勢，揭示近期長江源區的水資源變化規律及其制約因素，為長江源區水資源管理與生態環境保護提供支撐。

2 資料和方法

2.1 資料條件

長江源區范圍選擇以直門達水文站為界，大致范圍介于90°43′－96°45′E，32°30′－35°35′N之間，流域控制面積約13.78萬km2（見圖1）。地貌上以高原丘陵為主，河網水系發達，分布大小河流100余條，分屬沱沱河、當曲和楚瑪爾河為主的3個源區水系［4］。本文選取直門達1978－2009年共達32年的逐月徑流數據和長江源區沱沱河、五道梁、曲麻萊、玉樹、清水河共5個氣象站1978－2007年共30年的逐月降水、氣溫序列進行長江源區徑流變化規律及其影響因素的分析。

圖1 長江源區示意圖Fig.1 Study area of the headwaters of the Yangtze River

2.2 研究方法

2.2.1 Mann-Kendall趨勢分析法

利用Mann-Kendall非參數檢驗法分析長江源區水資源趨勢變化。Mann-Kendall非參數檢驗是Mann（1945）創建的一種時間序列趨勢分析方法，后經Kendall（1975）完善。該方法不需要樣本遵從一定的分布，也不受少數異常值的干擾，適用水文、氣象等非正態分布數據，計算簡便。它是世界氣象組織推薦的非參數檢驗方法之一，目前已被廣泛地用來分析降水、徑流和氣溫等要素時間序列的趨勢變化。

Mann-Kendall檢驗統計量可定義如下：

2.2.2 滑動 T檢驗法

利用滑動T檢驗法分析長江源區水資源要素的突變。

首先介紹傳統的T檢驗法。

設滑動點τ前后，2序列總體的分布函數各為F1（x）和 F2（x），今從總體 F1（x）和 F2（x）中分別抽取容量為n1和n2的2個樣本，要求檢驗原假設：F1（x）＝F2（x）。定義統計量為

式中：

T服從t（n1＋n2－2）分布，選擇顯著性水平 α，查 t分布表得到臨界值 tα／2，當｜T｜＞tα／2時，拒絕原假設，說明其存在顯著性差異，當｜T｜＜tα／2時，則接受原假設。

傳統的T檢驗法只能對已知變異點進行驗證而無法找出變異點，滑動T檢驗法是利用傳統的T檢驗法對序列逐點進行檢驗，對于滿足｜T｜＞tα／2所有可能的點τ，選擇使T統計量達到極大值的那一點作為所求的最可能變異點τ0。

圖2 直門達水文站年徑流量序列滑動T檢驗圖Fig.2 Slide T test of annual runoff series at Zhimenda Station

3 徑流變化特征

3.1 徑流年內分布特征

地表徑流的年內分配，主要取決于河流的補給類型［5］。長江源區徑流主要以降水和融雪補給為主，因而徑流年內分配主要受降水和氣溫影響。長江源區徑流年內分配不均，6－9月份徑流總量約占全年徑流量的72.5%。

3.2 變化趨勢分析

取顯著性水平α＝0.05，對直門達水文站逐月平均流量序列和年徑流序列進行趨勢分析，計算結果見表1。

表1 長江源區徑流序列趨勢分析匯總表Table 1 Statistical analysis on the trend of runoff series

由表1可以看出，直門達水文站徑流序列5月和7月呈遞減態勢，其它月徑流序列呈增加態勢，8－10月徑流序列增加最為明顯，年徑流序列也呈增加態勢，但趨勢都不顯著。

3.3 突變分析

取顯著性水平α＝0.05對直門達水文站逐月平均流量序列和年徑流序列進行突變分析（見圖3），計算結果見表2。

圖3 直門達水文站年徑流序列突變圖Fig.3 The change points of annual runoff at Zhimenda Station

表2 長江源區徑流序列突變分析匯總表Table 2 Statistical analysis on the sudden change point of runoff series

由表2可以看出，1－3月、8－10月徑流序列、年徑流序列發生顯著突變，變異時間多在2004年前后。

4 徑流變化的影響因素分析

4.1 氣候因素

4.1.1 降 水

取顯著性水平 α＝0.05，采用 Mann-Kendall法和滑動T檢驗法對長江源區典型站1978－2007年的年降水序列進行趨勢分析和突變，結果見3。分析發現，長江源區不同區域降水呈現不同的態勢（如圖4所示），玉樹站呈現減少態勢，其他站則呈現增加態勢，不過趨勢均不顯著。突變分析表明，長江源區降水序列發生了顯著的變異，五道梁和玉樹站在80年代發生突變，其他站突變點多在2002年前后（見圖4，表3）。

圖4 長江源區典型站年降水趨勢圖Fig.4 Trends of annual precipitation change at weather stations in the headwaters of Yangtze River

4.1.2 氣 溫

在顯著性水平 α＝0.05下，采用 Mann-Kendall趨勢檢驗對長江源區典型站1978－2007年的年平均氣溫序列分析發現，長江源區各站年平均氣溫均呈顯著增加趨勢（圖5）。

表3 長江源區年降水序列變化分析結果匯總表Table 3 Analysis on annual precipitation variation

圖5 長江源區典型站年平均氣溫序列趨勢圖Fig.5 Trends of average annual air temperature change at weather stations in the headwaters of Yangtze River

氣候變暖使得長江源區冰川融化。據調查［6］，2009年長江源區冰川總面積為1 051 km2，比《長江源志》記載的1 247 km2減少了196 km2，冰川年消融量達9.89億m3。長江源區冰川主要分布在北部的昆侖山、南部的唐古拉山及西部的祖爾肯烏拉山；在冰川持續大規模消融中，昆侖山玉珠峰冰川與1971年相比，冰舌退縮了1 500 m，平均每年退縮達42.86 m；唐古拉山口東側冰川側向最大退縮量為125 m，正面退縮265 m，與1970年相比，正面每年退縮量為7.57 m，退縮速度驚人。

4.2 下墊面因素

下墊面因素直接影響植被截流、蒸騰、下滲等水文過程，對徑流產生一定的影響。長江源區以草地為主，其面積占源區面積的58.9%，未利用土地面積占源區面積的32.3%，水域、林地、耕地、城鄉工礦居民用地，分別占源區面積的7.2%，1.6%，0.01%，0.004%。氣候變暖使高寒草甸發生逆向演替趨勢：高寒草甸→輕度退化高寒草地→中度退化高寒草地→重度退化高寒草甸→荒漠化地區［2］。自20世紀80年代以來，長江源區不斷加劇的高寒草甸草地退化使流域徑流系數減小、流域常遇徑流減少、稀遇洪水徑流發生頻率增加，這與嚴重退化極易形成暴雨洪水有關［4］。

5 結 語

通過對長江源區徑流序列和降水、氣溫序列的分析發現：

（1）長江源區近32年來年徑流量總體呈現增加態勢，但趨勢不顯著，且滑動T檢驗表明年徑流序列2004年前后發生突變。

（2）長江源區近32年來年降水序列和年平均氣溫序列呈現顯著增加趨勢，區域暖化趨勢沒有發生顯著變化，但近期干化趨勢有所緩解。降水序列在2002年前后發生突變，長江源區降水增加和冰雪消融是長江源區徑流增加的主要原因。

［1］ 王兆禮，陳曉宏，楊 濤.近50年東江流域徑流變化及影響因素分析［J］.自然資源學報，2010，25（8）：1365－1374.（WANG Zhao-li，CHEN Xiao-hong，YANG Tao.Runoff Variation and Its Impacting Factors in the Dongjiang River Basin During 1956－2005［J］.Journal of Natural Resources，2010，25（8）：1365－1374.（in Chinese））

［2］ 逯軍峰，董治寶，胡光印，等.長江源區土地利用／覆蓋現狀及成因分析［J］.中國沙漠，2009，29（6）：1043－1048.（LU Jun-feng，DONG Zhi-bao，HU Guang-yin，et al.Land Use／Land Cover Status and Forming Causes in Source Regions of Yangtze River［J］.Journal of Desert Research，2009，29（6）：1043－1048.（in Chinese））

［3］ 陳 芳，馬英芳，申紅艷，等.長江源區近44年氣候變化的若干統計分析［J］.氣象科技，2007，35（3）：340－344.（CHEN Fang，MA Ying-fang，SHEN Hong-yan，et al.Statistical Analysis of Climate Changes for Recent 44 Years over Source Region of Yangtze River［J］.Meteorological Science and Technology，2007，35（3）：340－344.（in Chinese））

［4］ 徐 明，馬超德.長江流域氣候變化脆弱性與適應性研究［M］.北京：中國水利水電出版社，2009.（XU Ming，MA Chao-de.Vulnerability and Adaptation of Climate Change in Yangtze River Basin［M］.Beijing：China Water Power Press，2009.（in Chinese））

［5］ 王菊英，丘玉俐.長江源年徑流量變化趨勢分析［J］.水利水電快報，2008，29（S1）：62－64，67.（WANG Juying，QIU Yu-li.Analysis on Annual Runoff Trend in the Headwaters of the Yangtze River［J］.Express Water Resources＆Hydropower Information，2008，29（S1）：62－64，67.（in Chinese））

［6］ 《節能環保》編輯部.長江源區冰川急劇退縮消融退縮速度令人堪憂［J］.節能環保，2009，（3）：6.（Editorial Office of Energy Conservation and Environmental Protection.Rapid Glacier Melting in the Headwaters of Yangtze River［J］.Energy Conservation and Environmental Protection，2009，（3）：6.（in Chinese） ）

Runoff Variation and Its Impacting Factors in the Headwaters of the Yangtze River in Recent 32 Years

ZHU Yan-long1，CHEN Jin2，CHEN Guang-cai2
（1.Qinghai Hydrology Bureau，Xining 810001，China；2.Yangtze River Scientific Research Institute，Wuhan 430010，China）

The runoff data measured from 1978 to 2009 at Zhimenda hydrological station in the headwaters of Yangtze River are analyzed by Mann-Kendall test and slide T test to study the trend of runoff variation.Impacting factors including precipitation，temperature，and underlying surface condition are also elaborated in this paper.The results show that the overall annual runoff increased slightly during 1978-2009；while there was an abrupt change around 2004 according to slide T test.With regard to the impacting factors，the annual precipitation witnessed slight changes in general and a sudden change around 2002，and the annual air temperature series within the basin increased remarkably.It is concluded that the runoff increase during 1978-2009 is mainly caused by the increase of precipitation and ice melting.

headwaters of the Yangtze River；runoff variation；climate change

S157

A

1001－5485（2011）06－0001－04

2010-02-24

水利部公益性行業科研專項長江流域用水總量控制管理與示范研究（201001005）

朱延龍（1962-），男，湖北襄陽人，高級工程師，主要從事水文測量及預報工作，（電話）13519703118（電子信箱）zhuyl208＠126.com。

（編輯：趙衛兵）