龍川江上游徑流量變化及其對氣候變化的響應

李 瀾, 丁文榮

(云南師范大學 旅游與地理科學學院, 昆明 650092)

?

龍川江上游徑流量變化及其對氣候變化的響應

李 瀾, 丁文榮

(云南師范大學 旅游與地理科學學院, 昆明 650092)

基于1970—2012年龍川江上游小河口水文站的徑流資料和楚雄的氣象資料，運用M-K法、相關系數法、累積距平法等方法，分析了龍川江上游徑流量、降水量和氣溫的年內、年際和年代際的變化以及徑流量對氣候變化的響應。結果表明：龍川江上游徑流量年內分配極不均，在1994年后漸趨均勻，年徑流量呈不顯著的波動下降趨勢，而且徑流量在70，80年代和2010—2012年偏枯，90和00年代偏豐；年平均氣溫和年降水量都呈上升趨勢，但后者的趨勢不明顯；年降水量對年徑流量的影響很明顯，而年平均氣溫對徑流量的影響較小，四個季節中春季和秋季的徑流量受降雨量變化的影響較大，降水量和氣溫的變化對夏季徑流的影響都較為明顯，而冬季相反。總之，不同季節的徑流量對不同氣候因素的響應是不同的。

龍川江上游; 徑流量; 氣候變化; 響應

水循環是由多個環節組成的自然過程，其中降水、蒸發和徑流是水循環過程中最主要的三個環節，這三者構成的水循環決定著全球的水量平衡，也決定著一個地區的水資源總量。水循環具有調節地球各圈層之間能量的作用，是氣候變化的重要影響因素，同時氣候變化也對水循環的各個環節具有重要的影響，尤其是對降水、蒸發和徑流等環節。

氣候的變化主要體現在氣溫和降水這兩大因子，而氣候變化對徑流的影響也主要體現在這兩個方面，而且氣溫和降水都在某種程度上影響徑流量的變化[1-3]。在以冰川冰雪融水補給為主的地區氣溫的影響是非常顯著的[4-5]，例如在青藏高原地區、新疆地區和中亞地區等。但很多研究表明降水的變化對徑流的影響要比氣溫的影響要大，降水對徑流量起主導作用[6-8]，其中夏季降水及其變化的影響突出[9-10]，而研究表明在西南的大部分地區徑流對氣溫的響應相對不顯著[11]。龍川江為金沙江南岸一級支流，發源于楚雄彝族自治州的南華縣天子廟坡東側，由西向東流，在元謀北部的江邊鄉匯入金沙江。龍川江流域不僅是楚雄州經濟、文化、政治中心，也是滇中經濟區重要組成部分，還是長江上游重要的產水區[12]。目前對龍川江流域已經進行了多方面的研究，但研究的范圍基本上都是整個區域[13-15]，對于龍川江上游的研究還很少，目前只有文安邦、張信寶等對上游泥沙輸移比進行研究[16]。由于龍川江上游是該流域主要的降水區，同時也是主要的產流區，本研究對于龍川江流域水資源的合理使用具有重要意義。

1　研究區概況

龍川江上游指的是小河口水文站以上地區，流經沙橋鎮、南華縣和楚雄市，流域面積為1 788 km2，上游河段總長96 km，河床比降3‰[16]。龍川江上游地區地形起伏較大，山高坡陡，河床切割較深。屬于低緯度高原季風氣候，氣候類型為亞熱帶季風氣候，深受南亞季風的影響，年平均氣溫為16.17℃，年平均降水861 mm。

2　資料與方法

本文所用資料為1970—2012年龍川江流域內的楚雄氣象站的逐月氣溫和降水量數據，以及小河口水文站的逐月徑流量數據，均來自于云南省水資源局。運用Mann-Kendall法分析徑流量、降水量和氣溫的變化趨勢[3，7]。運用累積距平法[17]和相關系數法分析徑流量與氣候因子關系及其響應。具體的方法：

(1) 非參數Man-Kendall趨勢檢驗。

定義檢驗統計量S：

(1)

式中：sgn()為符號函數，當xk-xi小于、等于或大于0時，sgn(xk-xi)分別為-1，0或1；M-K統計量公式S大于、等于或小于0時分別為：

(2)

Z為正值表示增加趨勢，負值表示減少趨勢，Z的絕對值在大于等于1.96時，表示通過了信度99.5%的顯著性檢驗。

(2) 非參數Mann-Kendall法突變檢測。

對于具有n個樣本量的時間序列，構造一秩序列：

(3)

(4)

(5)

(6)

定義統計量：

(7)

(3) 累積距平法。

(8)

式中:CAi為第一年至第j年的累積距平；x為平均值。CAi的正負極值點可能出現突變，對序列提取信號最強的幾個極點，且不從首尾各四年處提取。

3　結果與分析

3.1龍川江上游徑流量的變化特征

3.1.1月平均徑流量由于龍川江上游徑流的補給主要來自降水，而其氣候又受南亞季風的影響，干濕季分明，導致龍川江上游的徑流量年內分配極不均勻。徑流量主要集中在每年的6至10月(圖1)，此時的徑流量占全年的81.06%。而每年的1—5月和11—12月的徑流量只占全年的18.94%。因此龍川江上游徑流量集中分布在夏半年，每年6至10月為汛期。此外，龍川江上游的月徑流量在分別在1984年和2002年前后發生轉變，在1984年以前月徑流量不斷下降，之后一直處于上升狀態，而在2002年以后又轉為下降。而且在1984年前后年內徑流量分布曲線，由以雙峰曲線為主轉變為以單峰曲線為主，但在2010—2012年波峰不明顯。

總的來說，汛期徑流量占全年總徑流量的百分比是波動起伏的，而且在1994年前后汛期徑流量占全年總徑流量的百分比趨勢有所改變(圖2)。在1994年以前汛期徑流量占全年的百分比處于上升趨勢，汛期徑流量曾占到全年的93%，年內分配越來越不均勻。1994年以后汛期徑流量占全年的百分比呈現出明顯的下降趨勢(Z=-2.01，通過了0.05的顯著性檢驗)，即1994年以后龍川江上游徑流量年內分配漸趨均勻。

圖1　月平均徑流量分布

圖2　汛期徑流量占年徑流量的百分比

從年內季節分布來說，徑流量主要集中在夏秋兩季，分別占全年徑流量的52.33%和34.52%，而冬季徑流量占全年的百分比最小。而且在各年代的徑流量變化中，夏季徑流量變化最大，占總變化量的百分比也最大，其次是秋季，冬季最小(表1)。因此夏季徑流量是年徑流量變化的主要部分，冬季徑流量對年徑流的影響最小。

表1　各季節徑流量變化量占徑流總變化量的百分比　%

3.1.2年徑流量在1970—2012年期間，龍川江上游徑流量年際波動較大，徑流量變差系數為0.62。龍川江上游徑流量在近40多年里經歷了3個階段(圖3)：1970—1984年期間年徑流量下降；1984—2002年期間年徑流量波動上升；2002—2012年下降。通過非參數M-K單調趨勢檢驗計算得到的結果(表2)也與年徑流量5 a滑動曲線所呈現的階段趨勢一致，但均未通過顯著性檢驗。

通過非參數M-K突變檢測計算表明，在近40 a里，年徑流量只在1980—1984年和2002年顯著變化(前者下降，后者上升)(圖4)。并且年徑流量分別在1975年、1985年和2009年發生突變。非參數M-K單調趨勢檢驗計算得Z=-1.84，但未通過0.05的顯著性檢驗，即近40多年來龍川江上游年徑流量呈下降趨勢，但趨勢不明顯。

圖3　年徑流分布圖表2　龍川江上游年徑流量M-K單調趨勢檢驗

M-K單調趨勢檢驗1970—20121970—19841984—20022002—2012Z-1.8-1.051.26-0.54

注:│Z│≥1.96表示通過α=0.05顯著性水平，*表示通過顯著性檢驗。

圖4　年徑流量M-K突變檢測

3.1.3年代際徑流量龍川江上游各年代徑流量均有不同，差別較大(表3)。70年代徑流量偏枯，80年代徑流量最小，90年代徑流量上升幅度較大，是徑流量最大。00年代的徑流量有所減少，而2010—2012年三年的平均年徑流量比多年平均年徑流量小很多。總的來說，龍川江上游70，80年代和2010—2012年徑流量偏枯，90和00年代徑流量偏豐。

表3　龍川江上游各年代徑流量

3.2龍川江上游氣候變化特征

3.2.1氣溫龍川江上游地區地處亞熱帶季風氣候區，全年溫和，年溫差較小。在1970—2012年期間龍川江上游多年平均氣溫16.17℃，年內最大溫差為3.33℃，最高溫出現在6月，最低溫出現在1月。而且各年代各月平均氣溫具有上升的趨勢，以冬季變化最為明顯(圖5)。

圖5　年內月平均氣溫分布

從年平均氣溫的非參數M-K單調趨勢檢驗Z=5.369，通過了0.05的顯著性檢驗，即近40多年來龍川江上游氣溫顯著上升。年代際平均氣溫也呈現出上升趨勢，并且升高1.61℃。對龍川江上游氣溫的非參數M-K突變檢測得出，氣溫在1986年出現突變，并且在1986年之后氣溫顯著升高(圖6)。綜上所述，在1970—2012年期間龍川江上游年際氣溫顯著升高。

圖6　氣溫M-K突變檢測曲線

3.2.2降水由于受亞熱帶季風氣候的影響龍川江上游地區干濕季分明，降水年內分配極不均勻，降水量最大值出現在7月。降水量主要分布在每年6—10月(圖7)，即汛期，此時的降水量占全年降水量的88.3%。而且1984年以后汛期降水量占全年降水量的百分比幾乎都在85%以上，年內降水量更加分布不均。

圖7　年內月平均降水量分布

從年際來看，龍川江上游年降水量波動很大，最大相差857 mm。從百分比幾乎都在85%以上，年內降水量更加分布不均。趨勢線可以看出，年降水量呈現幅度很小的上升(圖8A)。但非參數M-K單調趨勢檢驗得出Z=0，即年降水量變化趨勢不明顯。年降水量M-K突變檢測表明年降水量在1972年，1977年和1984年左右發生突變，而且只在2002年發生明顯變化(圖8B)。另外，龍川江上游降水量在80年代和2010—2012年降水較少，90，00年代降水量豐沛(表4)。綜上所述，1970—2012年龍川江上游年降水量呈上升趨勢，且很不明顯。

圖8　1970-2012年龍川江上游年降水量變化曲線表4　1970-2012年龍川江上游各年代降水量分布

項目70s80s90s00s2010—2012年降水量/mm8514.37723.994249355.72020

3.3龍川江上游徑流量與氣候的關系

3.3.1徑流與降雨的關系年內各月降水量與徑流量的相關系數為0.928，為非常顯著正相關關系。徑流量與降水量都集中于6—10月，汛期與雨季相重合，其他月份降水量和徑流量都相對很小。降水量和徑流量的距平的最大值分別出現在7月和8月，即徑流量的波峰出現的時間要比降水量的推遲一個月，這是因為降水需要經過下滲、填洼等過程才能轉化為徑流。

從年際降水量與徑流量的關系來講，除了個別年份年降水量與年徑流量的變化趨勢相反，大多數年份的年降水量與年徑流量都呈顯著地正相關(圖9)，即降水量增加，徑流量增加；反之亦然，而且兩者的相關系數為0.868。從年代際來看，年代際的降水量與徑流量都在80年代和2010—2012年出現最小值，90年代出現最大值。從汛期與非汛期來講，徑流量與年降水量、汛期降水量、非汛期降水量都為正相關關系，而且通過了顯著性檢驗(表5)。所以綜上所述，降水量和徑流量的正相關關系顯著。

圖9　年際氣溫、降水量與徑流量距平

圖10　月氣溫徑流量的年內分布表5　龍川江上游徑流量與降水量、氣溫的相關系數

項目汛期徑流量非汛期徑流量年徑流量年降水量0.839**0.67**0.868**汛期降水量0.87**0.59**0.881**非汛期降水量0.421**0.556**0.475**年平均氣溫-0.0330.034-0.014汛期平均氣溫-0.126-0.028-0.108非汛期平均氣溫0.0100.0570.270

注:相關系數絕對值≥0.38表示通過α=0.01的顯著性水平，**表示通過顯著性檢驗。

3.3.2徑流與氣溫的關系月平均氣溫與徑流量的相關系數為0.69，表現出較為顯著地正相關性。而且兩者的年內分布都是單峰曲線，趨勢基本一致，但徑流量最大值比氣溫的推遲了兩個月(圖10)。從年平均氣溫與徑流量來看，兩者的相關系數只有-0.014 3，未能通過α=0.05的相關性檢驗，而且兩者的變化趨勢相差很大，因此年平均氣溫與年徑流量的相關性不明顯。年代際平均氣溫與徑流量都在80—90年代上升，90—00年代略微下降，而且它們的相關系數為-0.486 9，所以年代際平均氣溫與徑流量具有較為顯著的負相關。從汛期與非汛期來說，汛期平均氣溫與徑流量的相關系數都為負值，而非汛期氣溫相反，但均未通過顯著性檢驗(表5)。綜上所述，年內和年代際的平均氣溫與徑流量的相關性較為顯著，但年際的相關性不明顯。

3.4龍川江上游徑流量對氣候變化的響應

徑流量與氣候之間具有相當密切的關系，兩者相互聯系相互影響，氣候的變化影響徑流量的變化，在某種程度上來說徑流量的變化又是氣候變化的一種表現，即徑流量的變化是氣候變化的一種響應。降水和氣溫是氣候最主要的兩個因子，降水量和氣溫的變化都會直接或間接的影響徑流量的變化。

通過對比可見，年徑流量和降水量的變化過程幾乎一致，發生突變的年份也相同(圖11)，而且兩者的相關系數為0.868，年降水量的變化對年徑流的影響很大。而年徑流量與年平均氣溫的變化過程大概一致，但在某些階段的趨勢相差較大，例如在1984—1993年和2002—2012年。

圖11　龍川江上游年徑流量、年降水量和年平均氣溫累積距平

從四季來講，除了冬季以外，其他季節的徑流量和降水量的變化過程基本一致(圖12A和12B)，而且徑流量與降水量的相關系數都通過了0.01的顯著性檢驗，所以降水量的變化對春夏秋三個季節的徑流量影響很明顯，其中對夏季的影響最明顯。雖然夏秋冬三個季節的徑流量和氣溫的變化過程大概一致，但比徑流量與降水量的一致性差，一致性最差的是春季(圖12C和12D)。所以氣溫的變化對夏秋冬三季的徑流量的影響較小，對春季的影響最小。

從季節徑流與同期或者前期的降水量、氣溫的相關關系來看，春夏秋三個季節的徑流量與同期的降水量正相關顯著，而與同期氣溫為負相關(表6)，說明氣溫升高會使蒸發增大，導致徑流量減少，而且夏季氣溫與同期徑流量為顯著負相關。但春秋兩季的前期氣溫與徑流量為正相關，這可能是氣溫升高促進了水文循環，使得降水量增加的原因[1]。同時除春季以外，其他三個季節的徑流量都與前期降水量呈正相關，而且夏秋為顯著正相關，這是因為春季的降水量集中在后期，夏季和秋季又是雨季，而且降水形成徑流需要經過一定的過程與時間。而四季中冬季徑流量與同期和前期的降水、氣溫的相關性系數都沒有通過顯著性檢驗。因此，春季和秋季的徑流量受降雨量變化的影響較大，降水量和氣溫的變化對夏季徑流的影響都較為明顯，而冬季相反。總之，不同季節的徑流量對不同的氣候因素的響應是不同的。

圖12　龍川江上游四季徑流量、降水量和氣溫累積距平表6　龍川江上游季節徑流量與同期、前期年降水量、年平均氣溫的相關系數

徑流量同期降水量同期平均氣溫前期降水量前期平均氣溫春季0.648**-0.174-0.2650.025夏季0.868**-0.326*0.598**-0.256秋季0.618**-0.0080.593**0.036冬季-0.0110.0740.227-0.071

注:**和*分別表示通過0.01，0.05顯著性檢驗。

4　結 論

通過對龍川江上游43 a的徑流量、降水量和氣溫數據的分析表明：由于受亞熱帶季風氣候的影響，龍川江上游徑流的年內分配極其不均，但是在1994年以后有所變化，而且夏季徑流量對年徑流的影響最大；年徑流量呈不顯著的下降趨勢，并分為3個階段：1970—1984年(下降)，1984—2002年(波動上升)和2002—2012年(下降)；氣溫升高趨勢很明顯，而且在冬季特別顯著；雖然降水量也呈上升趨勢，但很不顯著，而且年內降水量分布極不均勻；降水量與徑流量具有顯著正相關關系，雖然年內和年代際平均氣溫與徑流量具有顯著正相關，但年均溫與年徑流量之間不存在顯著關系；相應的，徑流量對降水量的響應較為明顯，而對氣溫的響應則不明顯，其中夏季徑流量對降水量和氣溫的響應都很顯著，而冬季相反。

在龍川江上游，氣溫和降水量都在一定程度上影響著徑流量的變化，降水量的影響較為明顯，而氣溫則相反。這個結果與孫悅、陳玲飛等人的研究結論是一致的，但徑流量對氣候的響應率還有待進一步研究。而且現在人類活動對自然環境的影響越來越大，徑流量也不可避免的受其影響，在以后的研究中要注重人類活動的影響，并把其納入到相應的研究中。

[1]孫悅,李棟梁,朱擁軍.渭河徑流變化及其對氣候變化與人類活動的響應研究進展[J].干旱氣象,2013,31(2):396-405.

[2]魯鳳,徐建華,陳亞寧，等.葉爾羌河源流區近50年來年徑流變化及其對氣候變化的響應[J].第四紀研究,2010(1):152-158.

[3]許海麗,潘云,宮輝力，等.1959—2000年媯水河流域氣候變化與水文響應分析[J].水土保持研究,2012,19(2):43-47.

[4]Chikita K A, Wada T, Kudo I, et al. The Intra-Annual Variability of Discharge, Sediment Load and Chemical Flux from the Monitoring:The Yukon River, Alaska[J]. Journal of Water Resource and Protection, 2012,4(4):173-179.

[5]鞏同梁,劉昌明,劉景時.拉薩河冬季徑流對氣候變暖和凍土退化的響應[J].地理學報,2006,61(5):519-526.

[6]劉士余,章俊霞,羅志軍，等.近50 a贛西北大坑小流域徑流對降雨的響應[J].水土保持研究,2012,19(1):19-22.

[7]卓嘎,建軍,邊巴次仁.1960—2004年金沙江徑流量特征及其對氣候變化的響應[J].冰川凍土,2011,33(2):405-415.

[8]劉明春,李玲萍,史志娟，等.石羊河流域徑流量分布特征及對氣候變化的響應[J].干旱地區農業研究,2013,31(1):193-198.

[9]張一馳,李寶林,程維明，等.開都河流域徑流對氣候變化的響應研究[J].資源科學,2004,26(6):69-76.

[10]李卓侖,王乃昂,李育，等.近50 a來黑河出山徑流對氣候變化的響應[J].水土保持通報，2012,32(2):7-16.

[11]陳玲飛,王紅亞.中國小流域徑流對氣候變化的敏感性分析[J].資源科學，2004,26(6):62-67.

[12]何進花,丁文榮.龍川江流域近50 a氣溫、降水及徑流的變化趨勢分析[J].大理學院學報,2011,10(4):54-58.

[13]馬平森,顧世祥,卯昌書，等.龍川江流域徑流量變化趨勢及水資源合理配置研究[J].中國農村水利水電,2011(5):6-10.

[14]許志敏,臧慶春.龍川江流域年徑流特性的地區綜合法分析[J].人民長江,2011,42(10):73-76.

[15]丁文榮.龍川江水資源對氣候變化的敏感性和適應性分析[J].人民長江,2011,42(9):15-18.

[16]文安邦,張信寶,王玉寬，等.云貴高原區龍川江上游泥沙輸移比研究[J].水土保持學報,2003,17(4):139-141.

[17]李運剛,何大明,葉長青.云南紅河流域徑流的時空分布變化規律[J].地理學報,2008,63(1):41-49.

Runoff Variation and Its Response to Climate Change in the Upstream of Longchuanjiang River

LI Lan, DING Wenrong

(CollegeofTourismandGeographicScience,YunnanNormalUniversity,Kunming650092,China)

Based on the hydrologic data of Xiaohekou, the upstream of Longchuanjiang River, and the meteorological data of Chuxiong from 1970 to 2012, we used the methods of accumulative anomalies, Mann-Kendall and correlation coefficient to analyze the variation of inner-annual, inter-annual and inter-decadal for the runoff of the upstream of Longchuanjiang River, precipitation, the runoff and its response to climatic change. The results indicated that: (1) the annual runoff distributed unevenly, and it gradually become even, and annual runoff showed a wavelike decrease, but it wasn′t obvious, and runoff volume was less than normal years in 1970s, 1980s and the the period from 2010 ot 2012. The runoff volume in 1990s and 2000s was more than normal years; (2) temperature and annual runoff had showed a wavelike increase, but the latter was indistinctive; (3) annual precipitation effects on annual runoff was obvious, but the annual average temperature had less effect on runoff, precipitation of spring and autumn had larger effect on runoff, influences of precipitation, and temperature on runoff were obvious in summer, but were just opposite in the winter. In summary, responses of runoff in different seasons to different climate factors are different.

upstream of Longchuangjiang River; runoff; climatic change; response

2015-06-30

2015-08-07

國家自然科學基金項目(41101099、41261044);云南省社會發展科技計劃(2012CA024)

李瀾(1990—),女,山東臨沂人,碩士研究生,主要從事自然地理學研究工作。Email:lan900726@126.com

P339

A

1005-3409(2016)04-0083-06