張 杰， 張正棟， 萬露文， 楊傳訓(xùn)， 葉 晨

(華南師范大學(xué)地理科學(xué)學(xué)院, 廣州 510631)

氣候變化和人類活動(dòng)對(duì)汀江徑流量變化的貢獻(xiàn)

張 杰， 張正棟*， 萬露文， 楊傳訓(xùn)， 葉 晨

(華南師范大學(xué)地理科學(xué)學(xué)院, 廣州 510631)

基于汀江流域控制站(溪口站)1965—2012年日均徑流量實(shí)測(cè)數(shù)據(jù)及長汀、上杭、永定氣象站逐日降水量、逐日氣溫資料，采用年段均值比較、5a滑動(dòng)、趨勢(shì)線分析和累積距平等方法，分析流域近50年來的徑流量序列，揭示徑流量變化的趨勢(shì)，判別徑流量變化的突變年份，劃分徑流量變化的階段. 并應(yīng)用改進(jìn)后的累積量斜率變化率比較法定量估算了近50年來氣候變化和人類活動(dòng)對(duì)汀江徑流量變化的貢獻(xiàn)率. 結(jié)果表明：(1)徑流量總體呈現(xiàn)下降—上升—下降的波動(dòng)，判別突變年份為1972、2001年，劃分徑流變化階段為1965—1972年(A階段)、1973—2001年(B階段)、2002—2012年(C階段)；(2)與A階段相比，B、C階段氣候變化對(duì)徑流量變化的貢獻(xiàn)率分別為65.1%、50.5%，而人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率為34.9%、49.5%，可見人類活動(dòng)對(duì)徑流量變化貢獻(xiàn)日趨突出.

徑流量； 氣候變化； 人類活動(dòng)； 貢獻(xiàn)率； 汀江

氣候變化和人類活動(dòng)是影響徑流量變化的主要因素. 近幾十年來，隨著全球氣候不斷變化和人類活動(dòng)加劇，我國一些河流的徑流量明顯減少[1-3]，引起了水資源的嚴(yán)重短缺，激化了供需矛盾,而且對(duì)周邊及其下游、入海口產(chǎn)生一定程度的影響[4-6]. 因此，定量研究氣候變化和人類活動(dòng)對(duì)河流徑流量的影響，查明徑流量的變化特征和變化趨勢(shì)，分析其產(chǎn)生的原因，具有重要的實(shí)踐意義.

目前氣候變化和人類活動(dòng)對(duì)徑流量影響的定量研究方法主要是水文模型模擬法和數(shù)理統(tǒng)計(jì)分析法. 水文模型模擬法主要適用于小流域且具有較好的物理基礎(chǔ)，但參數(shù)的敏感性存在一定的不確定性，如果對(duì)模擬的結(jié)果不進(jìn)行驗(yàn)證，很有可能導(dǎo)致氣候變化對(duì)徑流量的影響偏大[7]. 數(shù)理統(tǒng)計(jì)分析法雖然比較簡(jiǎn)單，但需要較長的數(shù)據(jù)序列，且對(duì)水文氣象觀測(cè)資料要求較高，且長時(shí)間尺度的數(shù)據(jù)序列中的噪聲會(huì)對(duì)計(jì)算結(jié)果造成干擾. 2012年，王隨繼等[8]提出累積徑流量斜率變化率比較法，以消除噪聲對(duì)計(jì)算結(jié)果的干擾，定量研究氣候變化和人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率. 之后國內(nèi)大量學(xué)者在該方法的基礎(chǔ)上對(duì)長江流域、黃河流域和松花江流域徑流量變化進(jìn)行研究，發(fā)現(xiàn)這些流域徑流量在氣候變化和人類活動(dòng)的影響下基本呈現(xiàn)不同程度減小的趨勢(shì)[9-12]，并且人類活動(dòng)逐漸成為徑流量變化最重要的影響因素[13-15].

汀江哺育了閩西地區(qū)約70%(近200余萬)人口，是客家人的母親河，近年來隨著農(nóng)業(yè)、城市化進(jìn)程的加快和人口的不斷增多，流域內(nèi)人類活動(dòng)愈發(fā)活躍，使得水循環(huán)的時(shí)空分布特征以及生態(tài)水文都發(fā)生較為顯著變化，導(dǎo)致該流域出現(xiàn)水土流失、洪澇災(zāi)害頻發(fā)、極端降水等一系列的生態(tài)環(huán)境問題. 雖然對(duì)于汀江流域的徑流量已開展了一些研究[16-18]，但是很少有涉及氣候變化和人類活動(dòng)對(duì)徑流量變化貢獻(xiàn)的定量估算.

本文基于1965—2012年汀江流域下游控制站(溪口站)實(shí)測(cè)水文數(shù)據(jù)和流域內(nèi)上中下游3個(gè)國家氣象站(長汀站、上杭站和永定站)逐日降水量、逐日氣溫?cái)?shù)據(jù)，采用年段均值比較、5a滑動(dòng)、趨勢(shì)線分析、累積距平和改進(jìn)后的累積徑流量斜率變化率比較法，對(duì)近50年來氣候變化和人類活動(dòng)對(duì)汀江徑流量變化進(jìn)行定量分析，為汀江流域水資源管理和可持續(xù)利用提供有意義的參考.

1 研究區(qū)概況

本次研究區(qū)選擇的是汀江流域(圖1). 汀江是閩西最大河流，位于東經(jīng)115°59′～117°10′、北緯24°28′～26°02′，位于福建省西部，發(fā)源于武夷山南段東南一側(cè)的寧化縣治平鄉(xiāng)境內(nèi)木馬山北坡，植被覆蓋度比較好，流經(jīng)長汀縣、武平縣、上杭縣和永定縣，在永定縣峰市鎮(zhèn)進(jìn)入廣東省，至大埔縣三河壩與梅江匯合后稱韓江. 流域面積為9 022 km2,干流長度約285 km;支流眾多，流域面積大于500 km2以上的支流有：濯田河、桃瀾溪、舊縣河、黃潭河、永定河和金豐溪. 汀江地處亞熱帶季風(fēng)區(qū)，氣候溫和，降水充足，多年平均降水量為1 500～2 000 mm，雨量由北向南遞減，5—7月雨量集中，約占全年總雨量的60%. 汀江水哺育閩西約占70%的近200余萬人口，其下游的韓江水則是沿岸1 717.21余萬粵東百姓的生命之源.

圖1 汀江流域圖

2 數(shù)據(jù)來源及研究方法

2.1 研究數(shù)據(jù)

水文數(shù)據(jù)收集了汀江流域控制站(溪口站)1965—2012年日徑流量的實(shí)測(cè)數(shù)據(jù)，溪口站控制的流域面積為8 657 km2，占總流域面積的96%，實(shí)測(cè)數(shù)據(jù)較為完整，且研究時(shí)間序列較長，可以滿足研究的需要.

氣象數(shù)據(jù)選用了汀江流域長汀站、上杭站和永定站等3個(gè)國家氣象站(圖1)1965—2012年逐日降水量、逐日氣溫?cái)?shù)據(jù)，數(shù)據(jù)主要來源于中國氣象數(shù)據(jù)網(wǎng)(https://data.cma.cn/). 流域的面降水量采用3個(gè)站點(diǎn)的算術(shù)平均值表示. 雖然汀江流域內(nèi)只有3個(gè)國家站數(shù)據(jù)，氣象站點(diǎn)較少，但長汀站位于汀江上游，上杭站位于汀江中游，永定站位于汀江下游支流，站點(diǎn)分布較均勻，具有一定的代表性，所以選取這3個(gè)站可以滿足研究的需要，且各氣象站點(diǎn)氣象數(shù)據(jù)無明顯的突變點(diǎn)和隨機(jī)變化，數(shù)據(jù)變化較為一致，數(shù)據(jù)可靠，可以代表區(qū)域氣候狀況.

2.2 研究方法

2.2.1 數(shù)理統(tǒng)計(jì)分析方法 采用年段均值比較、5a滑動(dòng)、趨勢(shì)線分析和累積距平的方法，分析1953—2012年汀江流域徑流量、降水量、氣溫的變化特征和趨勢(shì). 采用累積距平法不僅可以分析徑流量、降水量、氣溫長時(shí)間序列變化特征，而且識(shí)別的徑流量、氣溫和降水量長時(shí)間序列突變年份清晰，突變前后雙累積量的線性關(guān)系擬合性較高，可以避免近年來被廣泛采用的降水量-徑流量雙累積曲線方法在判斷突變年份上存在的缺陷[19].

2.2.2 改進(jìn)的累積量斜率變化率比較法 本研究改進(jìn)了累積量斜率變化率比較法[8]，使用改進(jìn)后的累積量斜率變化率比較法分析氣候變化和人類活動(dòng)對(duì)徑流量變化的影響程度.

參考的累積量斜率變化率比較法過程如下：

假設(shè)累積徑流量-年份線性關(guān)系式的斜率(km3/a)在拐點(diǎn)前后2個(gè)時(shí)期分為SRb和SRa，累積降水量-年份線性關(guān)系式的斜率(mm/a)在拐點(diǎn)前后2個(gè)時(shí)期分為SPb和SPa，則累積徑流量-年份線性關(guān)系式的斜率變化率RSR(%)為:

RSR=100×(SRa-SRb)/SRb=(SRa/SRb-1)×100，

(1)

累積降水量-年份線性關(guān)系式的斜率變化率RSp(%)為:

RSP=100×(SPa-SPb)/SPb=(SPa/SPb-1)×100，

(2)

其中，RSR、RSP為正數(shù)時(shí)表示斜率增大，為負(fù)數(shù)時(shí)表示斜率減小. 降水量變化對(duì)徑流量變化的貢獻(xiàn)率CP(%)可以表示為：

CP=100×(RSP/RSR)=100×(SPa/SPb-1)/(SRa/SRb-1).

(3)

氣溫變化導(dǎo)致蒸發(fā)量變化而引起徑流量變化，如果用CET表示流域內(nèi)累積蒸散總量一年之間線性關(guān)系的斜率變化率(%)，則人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率CH(%)可以表示為:

CH=100-CP-CET.

(4)

徑流量變化是氣候變化和人類活動(dòng)共同作用的結(jié)果，所以，可將對(duì)徑流量的影響只分為氣候變化和人類活動(dòng)兩大類；而氣候變化對(duì)徑流量的影響主要體現(xiàn)為水、熱方面，且長時(shí)間序列持續(xù)觀測(cè)的蒸散量數(shù)據(jù)較難獲取，因此，可將各類氣候因子變化歸結(jié)為兩大因素：水和熱，體現(xiàn)為降水量和氣溫. 為了滿足氣候?qū)搅髁坑绊懹?jì)算的合理性，體現(xiàn)水熱綜合變化對(duì)徑流量變化影響，本文采用累積降水量-累積平均氣溫線性關(guān)系式的斜率變化率來代替單因子-年份線性關(guān)系式的斜率變化率.

以CC表示氣候變化對(duì)徑流量變化的貢獻(xiàn)率(計(jì)算方法同式(3))，則改進(jìn)后人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率可表示為:

CH=100-CC.

(5)

3 結(jié)果與分析

3.1 徑流量分析

由圖2～圖4及表1可知：(1)汀江流域的年徑流量變化波動(dòng)較大，尤其在1971—1980、2001—2012年：1971—1980年與1965—1970年相比，年代平均徑流量增加1.55 km3，即增加了21.08%；2001—2012年與1991—2000年相比，年代平均徑流量減少1.08 km2，即減少了11.9%. 年徑流量近50 a總體呈現(xiàn)微弱上升趨勢(shì)，但在近10 a徑流量明顯減少. (2)近50 a降水量總體呈現(xiàn)微弱上升的趨勢(shì)，在1971—1980、2001—2012年波動(dòng)較大：1971—1980年與1965—1970年相比，年代降水量增加101 mm，即增加了6.55%；2001—2012年與1991—2000年相比，年代降水量減少82 mm，即減少了4.82%. 多年平均氣溫總體呈現(xiàn)上升的趨勢(shì).

圖2 汀江流域徑流量年際變化和總體趨勢(shì)(1965—2012年)

Figure 2 Annual variations and general tendency of runoff of the Tingjiang Basin from 1965 to 2012

圖3 汀江流域降水量年際變化和總體趨勢(shì)(1965—2012年)

Figure 3 Annual variations and general tendency of precipitation of the Tingjiang Basin from 1965 to 2012

圖4 汀江流域年平均氣溫年際變化和總體趨勢(shì)(1965—2012年)

Figure 4 Annual variations and general tendency of air temperature of the Tingjiang Basin from 1965 to 2012

表1 汀江流域徑流量、降水量、平均氣溫年代變化及相對(duì)變化率(1965—2012年)Table 1 Decade changes of runoff，precipitation，air temperature and relative variation ratios of the Tingjiang Basin from 1965 to 2012

由此可知，在人類活動(dòng)影響較小的情況下，徑流量隨著氣候的變化呈現(xiàn)有規(guī)律的變化，但隨著人類利用和改造自然能力的不斷加強(qiáng)，日趨活躍的人類活動(dòng)也會(huì)對(duì)徑流量變化幅度產(chǎn)生影響，所以根據(jù)人類活動(dòng)的強(qiáng)弱，將徑流量變化的影響因素分為氣候變化作用、氣候變化和人類活動(dòng)共同作用，接下來根據(jù)不同影響因素劃分不同階段.

3.2 徑流量、降水量及平均氣溫突變分析及階段劃分

由圖5可知，汀江流域徑流量在1972 年前后呈現(xiàn)從減少到增加的趨勢(shì)，2001年前后呈現(xiàn)從增加到減少的趨勢(shì)，顯然徑流量發(fā)生突變年份為1972、2001年. 汀江流域的年降水量變化大致也呈現(xiàn)減少—增加—減少的趨勢(shì)，降水量突變的年份大致也為1972、2001年，年均溫大致呈現(xiàn)減少—增加的趨勢(shì)，1972年呈現(xiàn)減少趨勢(shì)，2001年呈現(xiàn)上升趨勢(shì)，因此可以根據(jù)累積距平值確定汀江流域的突變年份為1972、2001年.

根據(jù)上述突變年份把汀江流域徑流量變化劃分為3個(gè)階段：1965—1972年為基準(zhǔn)期(A階段)，徑流量主要受到氣候變化的影響；1973—2001年為第2個(gè)階段(B階段)，2002—2012為第3個(gè)階段(C階段)，B、C階段徑流量主要受到氣候變化和人類活動(dòng)的雙重影響，氣候變化和人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率是需要定量分析的問題.

3.3 氣候變化和人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)

利用徑流量的突變年份為分界點(diǎn)，分別對(duì)不同時(shí)期年份與累積徑流量、降水量與平均氣溫累積變化進(jìn)行線性回歸分析(圖6、圖7). 兩類關(guān)系式的決定系數(shù)R2均大于0.99，表明相關(guān)性很強(qiáng)，均達(dá)到0.01的顯著水平.

圖5 汀江流域徑流量、降水量、平均氣溫累積距平值(1965— 2012年)

Figure 5 The variations of cumulative anomaly of runoff，precipitation，and air temperature of the Tingjiang Basin from 1965 to 2012

圖6 汀江流域累積徑流量-年份線性關(guān)系(1965—2012年)

Figure 6 Linear relationships between year and cumulative runoff on the Tingjiang River Basin from 1965 to 2012

圖7 汀江流域累積降水量-累積平均氣溫線性關(guān)系(1965—2012年)

Figure 7 Linear relationships between cumulative air temperature and cumulative precipitation on the Tingjiang River Basin from 1965 to 2012

由表2、表3可知， B階段與A階段相比，累積徑流量-年份線性關(guān)系式的斜率變化量為1.836 5 km3/a，增加率為25.8%，累積降水量-累積平均氣溫線性關(guān)系式的斜率變化量為9.78 mm/℃ ，增加率為16.8%，由此根據(jù)式(3)、(4)可以得到氣候變化和人類活動(dòng)對(duì)汀江徑流量增加的貢獻(xiàn)率分別為65.1%、34.9%(表4)；C階段與A階段相比，累積徑流量-年份線性關(guān)系式的斜率變化量為1.437 2 km3/a，增加率為20.2%，累積降水量-累積平均氣溫線性關(guān)系式的斜率變化量為5.97 mm/℃，增加率為10.2%，由式(3)、(5)可以得到氣候變化和人類活動(dòng)對(duì)汀江徑流量增加的貢獻(xiàn)率分別為50.5%、49.5%(表4).

表2汀江流域累積徑流量-年份的線性關(guān)系式斜率及其變化率(1965—2012年)

Table 2 The slope of the linear relation ship between cumulative runoff and year and its change rate in Tingjiang Basin from 1965 to 2012

階段累積徑流量-年份線性關(guān)系式斜率/(km3·a-1)B、C階段與A階段比較C階段與B階段比較關(guān)系式斜率變化量/(km3·a-1)相對(duì)變化率/%關(guān)系式斜率變化量/(km3·a-1)相對(duì)變化率/%A7．1302————B8．96671．836525．8——C8．56741．437220．2-3．99-4．5

表3汀江流域累積降水量-累積平均氣溫的線性關(guān)系式斜率及其變化率(1965—2012年)

Table 3 The slope of the linear relationship between cumulative precipitation and accumulated mean temperature and its change rate in Tingjiang Basin from 1965 to 2012

階段累積降水量-累積平均氣溫線性關(guān)系式斜率/(mm·℃-1)B、C階段與A階段比較C階段與B階段比較關(guān)系式斜率變化量/(mm·℃-1)相對(duì)變化率/%關(guān)系式斜率變化量/(mm·℃-1)相對(duì)變化率/%A58．23————B68．019．7816．8——C64．205．9710．2-3．81-5．6

表4氣候變化和人類活動(dòng)對(duì)汀江流域徑流量變化的貢獻(xiàn)率(1965—2012年)

Table 4 The contribution of precipitation and human activities to runoff change in Tingjiang Basin from 1965 to 2012

階段B、C階段與A階段比較C階段與B階段比較氣候變化CC/%人類活動(dòng)CH/%氣候變化CC/%人類活動(dòng)CH/%A————B65．134．9——C50．549．5124-24

究其原因，人類活動(dòng)的影響增大導(dǎo)致汀江流域徑流量也處于增加的態(tài)勢(shì)，主要是由于人類活動(dòng)使流域內(nèi)植被減少、水土流失、城鎮(zhèn)化等，導(dǎo)致下墊面條件發(fā)生變化，地面不透水面積增加，流域內(nèi)透水性能和滯水性能變差，下滲減少，從而引起流域內(nèi)徑流量的增加[20]：(1)20世紀(jì)七八十年代，汀江流域城市化進(jìn)程加快，經(jīng)濟(jì)發(fā)展，人口增加，住宅和糧食問題日益突出，以及缺乏科學(xué)的指導(dǎo)，出現(xiàn)不合理的土地利用方式，1980—2010年流域土地利用方式發(fā)生較大變化，建筑用地增加238.03 km2，草地面積減少34.38 km2，林地面積減少34.38 km2，水域面積增多57.25 km2(數(shù)據(jù)來自地理空間數(shù)據(jù)云https://www.gscloud.cn/). 流域內(nèi)植被覆蓋度下降，不透水面積增加，水土流失嚴(yán)重，使流域徑流量增加. (2)汀江流域多山地，人們不合理的耕作方式(粗放坡地耕種)加劇了水土流失，導(dǎo)致匯入河流徑流量增多. 如農(nóng)民在25°以上的陡坡開荒，且順坡做成旱田，時(shí)種時(shí)荒. 雖然近年來，區(qū)域內(nèi)大于25°坡耕地基本退耕還林，但由于整個(gè)區(qū)域內(nèi)屬于山多地少、人多地缺矛盾較為突出的地區(qū)，且原有的耕作習(xí)慣短時(shí)間內(nèi)無法改變，仍存在較多的坡耕地種植旱糧等經(jīng)濟(jì)作物[21]. (3)汀江流域有5 個(gè)主要縣區(qū)存在不同程度的崩崗侵蝕,共有崩崗2 677個(gè),占福建省崩崗個(gè)數(shù)的40%,平均密度達(dá)0.297 個(gè)/km2,是全省平均值的5.5倍,汀江流域土地面積占福建省土地面積的7%,而崩崗總面積達(dá)2 086.86 hm2,占福建省崩崗總面積的79%，其中以長汀、上杭兩縣最為嚴(yán)重[22].

C階段與B階段相比，累積徑流量-年份線性關(guān)系式的斜率變化量為-0.399 a/km3，減少率為4.5%，累積降水量-累積平均氣溫線性關(guān)系式的斜率變化量為-3.81 mm/℃，減少率為5.6%，累積降水量-累積平均氣溫線性關(guān)系式的斜率減少率大于累積徑流量-年份線性關(guān)系式的斜率減少率，氣候變化和人類活動(dòng)對(duì)于徑流量變化的貢獻(xiàn)率為124%、-24% .

C階段與B階段相比，是假設(shè)B階段為基準(zhǔn)期，該時(shí)段徑流量變化只受到氣候變化的影響， C階段受到氣候變化和人類活動(dòng)的雙重影響. 通過式(3)、(5)計(jì)算可知，氣候變化的貢獻(xiàn)率超過了100%，所以說這個(gè)假設(shè)是不成立的，累積徑流量-年份線性關(guān)系式的斜率變化率下降的幅度小于累積降水量-累積平均氣溫線性關(guān)系式的斜率變化率的下降幅度. 究其原因，2002年以后一些大型水利工程的投入使用對(duì)徑流起到一定影響[23]. 如2002年永定縣境內(nèi)汀江干流棉花灘峽谷河段中部棉花灘水電站通過檢查驗(yàn)收投入使用，壩址以上控制流域面積7 907 km2，水庫正常蓄水位173.0 m，死水位146.0 m，調(diào)節(jié)庫容11.22億m3，校核洪水位177.8 m，相應(yīng)的總庫容20.35億m3[24]. 水利工程的建設(shè)對(duì)徑流量起到了調(diào)蓄作用，使累積徑流量-年份線性關(guān)系式的斜率變化率下降的幅度小于累積降水量-累積平均氣溫線性關(guān)系式的斜率變化率的下降幅度.

4 結(jié)論與討論

本文基于汀江流域控制站(溪口站)1965—2012年日均徑流量實(shí)測(cè)數(shù)據(jù)及長汀、上杭、永定氣象站逐日降水量、逐日氣溫?cái)?shù)據(jù)，采用數(shù)理統(tǒng)計(jì)方法和改進(jìn)后的累積量斜率變化率比較法揭示了汀江流域近半個(gè)世紀(jì)以來徑流量變化的趨勢(shì)及主要影響因素的貢獻(xiàn)率，主要結(jié)論如下：

(1)1965—2012年汀江流域徑流量變化波動(dòng)較大，呈現(xiàn)下降—上升—下降的趨勢(shì)，其中1971—1980年徑流量相對(duì)于60年代增加21.8%，2001—2012年徑流量相對(duì)于90年代減少了11.9%. 降水量總體呈現(xiàn)微弱上升的趨勢(shì)，在1971—1980、2001—2012年波動(dòng)較大：1971—1980年相對(duì)于1965—1970年增加6.55%，2001—2012年相對(duì)于1991—2000年減少4.82%. 平均氣溫總體呈現(xiàn)上升的趨勢(shì).

(2)汀江流域突變年份為1972 年和2001年，1965—1972年這個(gè)階段看做基準(zhǔn)期(A階段)，徑流量主要受到氣候變化的影響. 1973—2001年為第2個(gè)階段(B階段)，2002—2012為第3個(gè)階段(C階段)，B、C階段徑流量主要受到氣候變化和人類活動(dòng)的雙重影響. B、C階段相對(duì)于A階段徑流量變化率分別為25.8%、20.2%，降水-平均氣溫變化率分別為16.8%、10.2%，氣候變化對(duì)徑流量變化的貢獻(xiàn)率分別為65.1%、50.5%，人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率分別為34.9%、49.5%,這與林凱榮等[25]采用改進(jìn)的SCS模型將氣候變化和人類活動(dòng)對(duì)東江流域(濕熱地區(qū))徑流量影響的貢獻(xiàn)進(jìn)行分解研究得到各子流域中氣候變化和人類活動(dòng)對(duì)徑流影響分量相當(dāng)?shù)慕Y(jié)論相一致. 研究結(jié)果表明，人類活動(dòng)對(duì)徑流量變化的影響在不斷增強(qiáng)，即將超過氣候變化對(duì)徑流量的影響程度，因此需高度重視該流域的水資源利用.

本文對(duì)氣候變化和人類活動(dòng)對(duì)徑流量變化的影響主要采用了改進(jìn)后的累積量斜率變化率比較法，但只進(jìn)行了局部改進(jìn)，在研究氣候變化對(duì)徑流量變化的貢獻(xiàn)率時(shí)，氣候變化主要考慮了降水量和平均氣溫的變化(水熱變化)，并沒有考慮蒸發(fā)對(duì)徑流量變化的影響，如何將降水量、平均氣溫、蒸發(fā)量對(duì)徑流量影響結(jié)合起來，是今后開展定量研究需要進(jìn)一步攻克的地方. 目前徑流量變化的定量研究方法主要分為統(tǒng)計(jì)分析法和水文模型模擬法，所以在之后的研究中采用水文模型模擬法[26-28]從形成機(jī)理上對(duì)徑流量變化進(jìn)行定量研究，對(duì)2種方法的研究結(jié)果進(jìn)行對(duì)比是檢驗(yàn)其合理性的方法之一.

[1] 盧曉寧,洪佳,王玲玲. 近50年來黃河入海水沙的多時(shí)間尺度特性分析[J]. 水土保持研究,2015,22(1):115-122.

LU X N,HONG J,WANG L L. Multiple time-scale characteristics variations of runoff and sedimentloading of yellow river to sea in recent 50 years[J]. Research on Soil and Water Conservation,2015,22(1):115-122.

[2] 劉宇峰,孫虎,原志華. 基于小波分析的汾河河津站徑流與輸沙的多時(shí)間尺度特征[J]. 地理科學(xué),2012，32(6):764-770.

LIU Y F,SUN H,YUAN Z H. Multitime scale features of runoff and sediment discharge on wavelet analysis at Hejin Station of Fenhe River Basin[J]. Scientia Geographica Sinica,2012,32(6):764-770.

[3] 陳亞寧,李稚,范煜婷,等. 西北干旱區(qū)氣候變化對(duì)水文水資源影響研究進(jìn)展[J]. 地理學(xué)報(bào),2014,69(9):1295-1304.

CHEN Y N,LI Z,FANG Y T,et al. Research progress on the impact of climate change on water resources in the arid region of Northwest China[J]. Journal of Geographical Sciences,2014,69(9):1295-1304.

[4] 吳創(chuàng)收,楊世倫,黃世昌,等. 1954-2011年間珠江入海水沙通量變化的多尺度分析[J]. 地理學(xué)報(bào),2014,69(3):422-432.

WU C S,YANG S L,HUANG S C,et al. Multi-scale variability of water discharge and sediment load in the Pearl River during 1954-2011[J]. Journal of Geographical Sciences,2014,69(3):422-432.

[5] 任惠茹,李國勝,崔林林,等. 近60年來黃河入海水沙通量變化的階段性與多尺度特征[J]. 地理學(xué)報(bào),2014,69(5):619-631.

REN H R,LI S G,CUI L L,et al. Phases and periodic changes of water discharge and sediment load from the Yellow River to the Bohai Sea during 1950-2011 [J]. Journal of Geographical Sciences,2014,69(5):619-631.

[6] 胡光偉,毛德華,李正最,等. 荊江三口60a來入湖水沙變化規(guī)律及其驅(qū)動(dòng)力分析[J]. 自然資源學(xué)報(bào),2014,29(1):129-142.

HU G W,MAO D H,LI Z Z,et al. Characteristics and driving factors of runoff and sediment changes fluxes into the Dongting Lake from three outlets in Jingjiang River during past 60 years[J]. Journal of Natural Resources,2014,29(1):129-142.

[7] LEGESSE D，VALLET-COULOMB C，GASSE F. Hydrological response of a catchment to climate and landuse changes in Tropical Africa case study South Central Ethiopia[J]. Journal of Hydrology,2003,275(1):67-85.

[8] 王隨繼,閆云霞,顏明,等. 皇甫川流域降水和人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率分析——累積量斜率變化率比較方法的提出及應(yīng)用[J]. 地理學(xué)報(bào),2012,67(3):388-397.

WANG S J,YAN Y X,YAN M,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.

[9] 王隨繼,李玲,顏明. 氣候和人類活動(dòng)對(duì)黃河中游區(qū)間產(chǎn)流量變化的貢獻(xiàn)率[J]. 地理研究,2013,32(3):395-402.

WANG S J,LI L,YAN M. The contributions of climate change and human activities to the runoff yield changes in the middle Yellow River Basin[J]. Geographical Research,2013,32(3):395-402.

[10] 王彥君,王隨繼,蘇騰. 1955-2010年松花江流域不同區(qū)段徑流量變化影響因素定量評(píng)估[J]. 地理科學(xué)進(jìn)展,2014,33(1):65-75.

WANG Y J,WANG S J,SU T. Quantification of impacts of precipitation and human activities on runoff changes in different sections of Songhua River during 1955-2010[J]. Progress in Geography,2014,33(1):65-75.

[11] 張曉婭,楊世倫. 流域氣候變化和人類活動(dòng)對(duì)長江徑流量影響的辨識(shí)(1956～2011)[J]. 長江流域資源與環(huán)境,2014,23(12):1729-1739.

ZHANG X Y,YANG S L. Climatic and anthropogenic impacts on wathr discharge in the Yangtze River over the last 56 years(1956—2011)[J]. Resources and Environment in the Yangtze River Basin,2014,23(12):1729-1739.

[12] 王彥君,王隨繼,蘇騰. 降水和人類活動(dòng)對(duì)松花江徑流量變化的貢獻(xiàn)率[J]. 自然資源學(xué)報(bào),2015,30(2):304-314.

WANG Y J,WANG S J,SU T,et al. Contributions of precipitation and human activities to runoff change in the Songhua River Basin [J]. Journal of Natural Resources,2015,30(2):304-314.

[13] 馬龍,劉廷璽,馬麗,等. 氣候變化和人類活動(dòng)對(duì)遼河中上游徑流變化的貢獻(xiàn)[J]. 冰川凍土,2015,37(2):470-479.

MA L,LIU T X,MA L,et al. The effect of climate change and human activities on the runoff in the upperand middle reaches of the Liaohe River Inner Mongoli [J]. Journal of Glaciology and Geocryology,2015,37(2):470-479.

[14] 邱臨靜,鄭粉莉,尹潤生,等. 降水變化和人類活動(dòng)對(duì)延河流域徑流影響的定量評(píng)估[J]. 氣候變化研究進(jìn)展,2011,7(5):357-362.

QIU L J,ZHENG F L,YIN R S,et al. Quantification of impacts of precipitation change and human activities on streamflow in the Yanhe River Basin[J]. Advances in Climate Change Research,2011,7(5):357-362.

[15] 夏智宏,劉敏,王苗,等. 1990s以來氣候變化和人類活動(dòng)對(duì)洪湖流域徑流影響的定量辨識(shí)[J]. 湖泊科學(xué),2014,26(4):515-521.

XIA Z H,LIU M,WANG M,et al. Quantitative identification of the impact of climate change and human activity on runoff in Lake Honghu Basin since 1990s[J]. Journal of Lake Sciences,2014,26(4):515-521.

[16] 邱宇. 汀江流域水環(huán)境安全評(píng)估[J]. 環(huán)境科學(xué)研究,2013,26(2):152-159.

QIU Y． The evaluation of the Tingjiang River Basin water environment safety [J]. Research of Environmental Sciences,2013，26(2):152-159.

[17] 周素萍,張興奇. 汀江上游水沙變化特征分析[J]. 水文,2010,30(2):89-92.

ZHOU S P,ZHANG X Q． Characteristics analysis of runoff and sediment in upper reach of Tingjiang River [J]. Journal of China Hydrology,2010,30(2):89-92.

[18] 王海龍. 基于小波分析的汀江水沙通量變化規(guī)律研究[J]. 泥沙研究,2012(2):61-67.

WANG H L． Characteristics of water discharge and sediment load variability of Tingjiang River by wavelet analysis [J]. Journal of Sediment Research,2012(2):61-67.

[19] RAN L S,WANG S J,FAN X L. Channel change at tou-daoguai 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.

[20] 李艷,陳曉宏,王兆禮. 人類活動(dòng)對(duì)北江流域徑流系列變化的影響初探[J]. 自然資源學(xué)報(bào),2006,21(6):910-915.

LI Y,CHENG X H,WANG Z L. A tentative discussion on the impact of human activities on the variability of runoff series of the Beijing River Basin[J]. Journal of Natural Resources,2006,21(6):910 -915.

[21] 李守中,余麗凡,鄒怡靜,等. 長汀侵蝕退化地區(qū)生態(tài)恢復(fù)過程中植被多樣性與土壤肥力的變化研究[J]. 安徽農(nóng)學(xué)通報(bào),2008,14(19):49-51;81.

LI S Z,YU L F,ZOU Y J,et al. Changes of vegetation diversity and soil fertility during the ecological restoration process in an eroded area of Changting County[J]. Anhui Agricultural Science Bulletin,2008,14(19):49-51;81.

[22] 朱頌茜. 汀江流域水土流失現(xiàn)狀及特點(diǎn)分析[J]. 亞熱帶水土保持,2005,17(1):44-47.

[23] 李宇. 梅江、汀江調(diào)峰影響河段設(shè)計(jì)水位的推算及航道整治研究[D]. 南京：河海大學(xué),2007.

LI Y. Study on the influence of Meijiang and Dingjiang reckoning and waterway regulation reach peak water level design [D]. Nanjing:Hohai University,2007.

[24] 夏讓欣. 福建山區(qū)水利建設(shè)的調(diào)查與思考[J]. 中國水利,2002(9):110-111.

[25] 林凱榮,何艷虎,陳曉宏. 氣候變化及人類活動(dòng)對(duì)東江流域徑流影響的貢獻(xiàn)分解研究[J]. 水利學(xué)報(bào),2012,43(11):1312-1321.

LIN K R,HE Y H,CHEN X H. Identifying the quantitative effect of climate change and human activity on runoff in the Dongjiang River Basin [J]. Journal of Hydraulic Engineering,2012,43(11):1312-1321.

[26] 趙陽,余新曉,鄭江坤,等. 氣候和土地利用變化對(duì)潮白河流域徑流變化的定量影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(22):252-260.

ZHAO Y,YU X X,ZHENG J K,et al. Quantitative effects of climate variations and land-use changes on annual streamflow in Chaobai River Basin[J]. Transactions of the Chinese Society of Agricultural Engineering,2012,28(22):252-260.

[27] 袁宇志,張正棟,蒙金華. 基于SWAT模型的流溪河流域土地利用與氣候變化對(duì)徑流的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào),2015,26(4):989-998.

YUAN Y Z,ZHANG Z D,MENG J H. Impact of changes in land use and climate on the runoff in Liuxihe Watershed based on SWAT model[J]. Journal of Applied Ecology,2015,26(4):989-998.

[28] 鄧曉宇,張強(qiáng),孫鵬,等. 氣候變化和人類活動(dòng)對(duì)信江流域徑流影響模擬[J]. 熱帶地理,2014,34(3):293-301.

DENG X Y,ZHANG Q,SUN P,et al. Assessment of impacts of climate change and human activities on runoff with HSPF for the Xinjiang River Basin[J]. Tropical Geo-graphy,2014,34(3):293-301.

The Contributions of Climate Change and Human Activities on Runoff of the Tingjiang River

ZHANG Jie， ZHANG Zhengdong*， WAN Luwen， YANG Chuanxun， YE Chen

(School of Geography Science, South China Normal University, Guangzhou 510631, China)

Based on 1965-2012 average daily runoff measured data of control hydrometric station (Xi Kou) in the Tingjiang Basin，and daily precipitation and daily temperature of three meteorological stations in the Tingjiang Basin (Chang Ting, Shang Hang,Yong Ding). Paragraph mean comparison, 5a sliding trend line analysis and cumulative anomaly methods are used to analyze runoff series in recent 50 years. The trend of runoff variation is revealed. The mutation years of runoff variation are identified. The stages of runoff variation are divided. The improved cumulative slope rate comparison method is used to quantitatively analyze the contribution rate of the climate change and human activities to runoff variation of Tingjiang in nearly 50 years. The results show that：(1)The runoff shows a down-up-down fluctuation in general, and the mutation year is 1972 and 2001, the variation stages of runoff are 1965-1972 (stage A),1973-2001 (stage B),2002-2012 (stage C). (2)Comparing with stage A, the contribution rates of climate change to runoff variation in stage B and C are 65.1% and 50.5%, respectively, while the contribution rate of human activities to runoff variation is 34.9% and 49.5%, respectively. This shows that the contribution of human activities to the change of runoff is increasing prominent.

runoff； climate change； human activity； contribution rate； Tingjiang River

2016-05-03 《華南師范大學(xué)學(xué)報(bào)(自然科學(xué)版)》網(wǎng)址：http://journal.scnu.edu.cn/n

國家自然科學(xué)基金項(xiàng)目(41471147)

*通訊作者:張正棟，教授，Email:zhangzdedu@163.com.

P333.6

A

1000-5463(2017)06-0084-08

【中文責(zé)編：莊曉瓊 英文審校：葉頎】