氣候變化和人類活動(dòng)對(duì)水資源影響的初探——以陜西省為例

賀燕子,岳大鵬,達(dá) 興,程金文,張君茹

(陜西師范大學(xué) 旅游與環(huán)境學(xué)院，陜西 西安 710119)

?

氣候變化和人類活動(dòng)對(duì)水資源影響的初探
——以陜西省為例

賀燕子,岳大鵬*,達(dá) 興,程金文,張君茹

(陜西師范大學(xué) 旅游與環(huán)境學(xué)院，陜西 西安 710119)

針對(duì)陜西省14個(gè)小流域1961～2009年的徑流、降雨及潛在蒸散情況，運(yùn)用Mann-Kendall非參數(shù)檢驗(yàn)法、累積距平法分析了其變化趨勢(shì)及徑流突變點(diǎn)，使用敏感系數(shù)法分離了人類活動(dòng)和氣候變化對(duì)徑流的影響，并基于Budyko假設(shè)，采用張氏曲線研究了各流域ω、β以及γ的變化。結(jié)果表明:(1)研究區(qū)徑流呈現(xiàn)顯著下降的趨勢(shì)，降雨量呈現(xiàn)不顯著下降的趨勢(shì)，潛在蒸散發(fā)的變化各區(qū)域不一致并且不同流域徑流突變點(diǎn)出現(xiàn)的時(shí)間不一致，主要集中在1980年和1990年前后。(2)ω在陜北關(guān)中從北到南逐漸增大，在陜南比較??；β陜南普遍較高，γ陜南普遍較高。(3)陜北由于實(shí)施水土保持、退耕還林、黃土治理等措施，人類活動(dòng)對(duì)徑流變化的貢獻(xiàn)率高于氣候變化；關(guān)中地區(qū)由于工農(nóng)業(yè)用水以及生活用水，人類活動(dòng)的貢獻(xiàn)率有所下降；在陜南由于下墊面性質(zhì)，人類活動(dòng)的貢獻(xiàn)率亦占有很大比重。

徑流變化；Mann-Kendall；累積距平法；敏感系數(shù)法；陜西

0 引言

隨著全球氣候變化和人類土地利用方式的改變，水資源正發(fā)生著顯著的變化，主要表現(xiàn)為河道內(nèi)徑流量的不斷減少。隨著工農(nóng)業(yè)的發(fā)展，特別是進(jìn)入到20世紀(jì)80年代，水資源供需矛盾日漸突出，人類活動(dòng)對(duì)其影響日益顯著；徑流的變化特征已經(jīng)不單純地由氣候因素的時(shí)空變化規(guī)律所決定，人類活動(dòng)的影響使徑流的變化規(guī)律更趨向復(fù)雜化,氣候變化和人類活動(dòng)是引起流域水資源變化的主要原因，國(guó)內(nèi)外的學(xué)者均對(duì)徑流的變化及其受氣候和人類活動(dòng)的影響給予了很高的關(guān)注。Nemec等[1]指出若氣溫升高1 ℃，且降水減少10%，干旱和濕潤(rùn)地區(qū)徑流的變化幅度分別為25%和50%，這說(shuō)明了在徑流的響應(yīng)中氣候因子的變化會(huì)被放大。Pikounism等[2]利用分布式水文模型SWAT(Soil and Water Assessment Tool)研究了不同土地利用變化對(duì)水文循環(huán)的影響。李棟梁等[3]指出了青藏高原前期熱力作用及大氣環(huán)流特征對(duì)黃河流域夏秋季的徑流量預(yù)測(cè)具有一定的指示性。張建云等[4]對(duì)中國(guó)近50年來(lái)6大流域的年徑流變化趨勢(shì)進(jìn)行了研究，結(jié)果表明受人類活動(dòng)和氣候變化的共同影響，6大江河的實(shí)測(cè)徑流量整體呈下降趨勢(shì)，其中黃河中下游地區(qū)的徑流量減少十分明顯。雖然有關(guān)氣候變化和人類活動(dòng)對(duì)徑流的影響研究已有很多，然而以往的研究主要集中在流域[5-14]，針對(duì)具體行政管理單元的并不多見(jiàn)。

氣候變化和人類活動(dòng)的加劇，使流域中的水文循環(huán)、水資源時(shí)空、數(shù)量分配產(chǎn)生了很大的變化并直接影響到流域水資源的蒸發(fā)、降雨以及徑流,在這種背景下,許多流域中的徑流量減少甚至出現(xiàn)了季節(jié)性斷流的現(xiàn)象,從而給流域水資源的開(kāi)發(fā)、利用、調(diào)配和管理帶來(lái)了很大的挑戰(zhàn)。所以對(duì)人類活動(dòng)和氣候變化對(duì)流域水資源的影響做出研究,能夠?yàn)榱饔蛩Y源的開(kāi)發(fā)利用、設(shè)計(jì)規(guī)劃等提供必要的依據(jù)。

陜西根據(jù)其地貌類型的不同，可以分為三部分：陜北黃土高原、關(guān)中平原、陜南秦巴山地。陜北黃土高原是生態(tài)環(huán)境脆弱區(qū)，是全國(guó)水土流失最為嚴(yán)重的地區(qū)之一。主要表現(xiàn)在2個(gè)方面:一是土地荒漠化，二是水土流失嚴(yán)重。關(guān)中平原是我國(guó)重要的工農(nóng)業(yè)集中區(qū)，是我國(guó)工農(nóng)業(yè)和文化發(fā)達(dá)地區(qū)之一，是全國(guó)重要麥、棉產(chǎn)區(qū)。陜南秦巴山地是我國(guó)南水北調(diào)中線水源地，對(duì)于我國(guó)水資源優(yōu)化配置、改善北方地區(qū)生態(tài)系統(tǒng)、保障經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展將發(fā)揮巨大作用。因此明晰不同區(qū)域的演變特征，對(duì)水資源管理和適應(yīng)氣候變化具有重要意義。

1 研究區(qū)概況

本次研究所選區(qū)域位于陜西省，根據(jù)其地貌類型的不同，陜西可以分為三部分：北部是在第三紀(jì)末起伏和緩的準(zhǔn)平原基礎(chǔ)上，歷經(jīng)第四紀(jì)以來(lái)多次黃土堆積和侵蝕作用的陜北黃土高原；中部是關(guān)中斷陷盆地，它南依秦嶺，北連黃土高原，為西狹東闊的新生代斷陷盆地，渭河橫貫其中，盆地兩側(cè)地形向渭河傾斜，由洪積傾斜平原、黃土臺(tái)塬、沖積平原組成，呈階梯狀地貌景觀；南部是由隴山余脈、秦嶺和巴山組成的陜南秦巴山地，是中生代末以來(lái)全面隆起的褶皺山地，以中山地貌為主體，高峰林立，斷陷盆地星散于群山之中，漢江谷地貫穿于秦嶺、巴山之間。由于其地貌類型、地形特征等地理因素的不同導(dǎo)致陜北、關(guān)中以及陜南的水資源使用特征有明顯的差異。陜北的水資源使用主要是在煤化工業(yè)，關(guān)中是我國(guó)重要的工農(nóng)業(yè)集中區(qū)，并且關(guān)中地區(qū)人口密集，生活用水量大，而陜南是我國(guó)南水北調(diào)工程的中線(圖1)。

圖1 各流域位置及水文站分布

2 數(shù)據(jù)和方法

2.1 數(shù)據(jù)來(lái)源

根據(jù)全國(guó)氣象站點(diǎn)分布，選用資料較完整的氣象站1961～2009年的年累積降雨量資料數(shù)據(jù)，氣象資料來(lái)源于中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)，采用FAO推薦的Penman-Monteith公式計(jì)算潛在蒸散發(fā)量[15]，流域徑流量選用由黃河及長(zhǎng)江水利委員會(huì)(水資源公報(bào))提供的各水文觀測(cè)站1961～2009年逐年實(shí)測(cè)徑流量資料。

2.2 研究方法

2.2.1 Mann—Kendall檢驗(yàn) 用Mann-Kendall非參數(shù)趨勢(shì)檢驗(yàn)法檢測(cè)要素序列變化趨勢(shì)[16-17]。該統(tǒng)計(jì)檢驗(yàn)方法的優(yōu)點(diǎn)是統(tǒng)計(jì)測(cè)試的樣本不需要服從一定的分布，也不受少數(shù)異常值的干擾，適用于非正態(tài)分布的水文氣象等數(shù)據(jù)。在Mann-Kendall檢驗(yàn)中，原假設(shè)H0為時(shí)間序列數(shù)據(jù)(x1,x2,…,xn),是n個(gè)獨(dú)立的，隨機(jī)分布的樣本；對(duì)于所有的k、j≤n,且k≠j，xk和xj的分布式是不相同的，檢驗(yàn)的統(tǒng)計(jì)變量S計(jì)算公式如下：

(1)

(2)

S為正態(tài)分布，其均值點(diǎn)為0，Var(S)=n(n-1)(2n+5)/18。當(dāng)n>10時(shí)，標(biāo)準(zhǔn)的正態(tài)統(tǒng)計(jì)變量通過(guò)下式計(jì)算：

(3)

在給定的α置信水平上，如果|Z|≥Z1-a/2，則原假設(shè)是不可接受的，即在α置信水平上，時(shí)間序列數(shù)據(jù)存在明顯的上升或下降趨勢(shì)。對(duì)于統(tǒng)計(jì)變量Z>0時(shí)，樣本是上升趨勢(shì)；而Z<0時(shí)，樣本則是下降趨勢(shì)。

累積距平是一種十分常用的、根據(jù)曲線直觀判斷變化趨勢(shì)的方法。對(duì)于某一序列x某一時(shí)刻t的累積距平表示為:

(4)

2.2.3 敏感系數(shù)法 采用多年平均水量平衡法，假設(shè)前后時(shí)段年均徑流的變化量可由下式計(jì)算：

△Q=Q2-Q1

(5)

作為一階近似值，同時(shí)前后時(shí)段年均徑流變化量也可由下式確定：

△Q==Qclim-QLUCC

(6)

Qclim通過(guò)Milly等[10]提出的方法計(jì)算：

Qclim=P+E0

(7)

式中：△Q為年均徑流變化量，單位為mm；Q1為前一時(shí)段年均徑流量，單位為mm；Q2為后一時(shí)段年均徑流量，單位為mm；Qclim為由于氣候變化引起的前后時(shí)段徑流變化量，單位為mm；QLUCC為由土地利用/覆被變化不同等人類活動(dòng)引起的徑流變化量，單位為mm；△P、△E0為前后時(shí)段降雨量(P)和潛在蒸發(fā)量(E0)的變化量，單位為mm；β和γ分別為徑流對(duì)降雨及蒸發(fā)的敏感系數(shù)，決定于區(qū)域氣候環(huán)境特征，可由Li等[19]提供的公式確定：

(8)

(9)

式中：x為干燥指數(shù)(E0/P)，其量綱為1；ω為反映流域植被變化的參數(shù)，其量綱為1。

2.2.4 估算地表系數(shù) Budyko認(rèn)為可用的能量和水是決定蒸散速度的主要因素，并且建立了一個(gè)基于干燥指數(shù)來(lái)確定蒸散量的模型[19]。張魯?shù)膹埵锨€模型與Budyko的假設(shè)類似，用于估算蒸散量[16,19]。該模型不僅符合先前的理論工作并與全世界的研究達(dá)到很高的一致[20-22]。

(10)

式中：E0是區(qū)域潛在蒸散發(fā)，ω是與植被類型有關(guān)的下墊面參數(shù)，E可由水量平衡推算，公式為：

Q=P-E

(11)

3 結(jié)果與分析

3.1 徑流降水以及潛在蒸散的變化趨勢(shì)分析

根據(jù)全國(guó)氣象站點(diǎn)分布，在ArcGIS 10.0中，將獲取到的站點(diǎn)尺度的降雨數(shù)據(jù)以及計(jì)算的站點(diǎn)尺度的潛在蒸散發(fā)數(shù)據(jù)，使用克里金(kring)插值的方法生成研究區(qū)面狀降雨數(shù)據(jù)及潛在蒸散發(fā)數(shù)據(jù)，并統(tǒng)計(jì)得到流域尺度降雨以及潛在蒸散發(fā)的平均值。借助Matlab執(zhí)行Mann-Kendall趨勢(shì)分析方法，各流域徑流、降雨以及潛在蒸散發(fā)的變化趨勢(shì)(圖2)。

從圖2以及表1的徑流變化趨勢(shì)中可以看出：除鎮(zhèn)安外，各流域的徑流呈現(xiàn)顯著下降的趨勢(shì)。但不同區(qū)域徑流變化的斜率不等，其中在陜北趙石窯、關(guān)中林家村、陜南石泉斜率最小，徑流下降最明顯。從圖2以及表1的降雨變化趨勢(shì)中可以看出：各流域降雨均呈現(xiàn)不顯著下降的趨勢(shì)。但不同區(qū)域徑流變化的斜率不等，其中陜北橫山、關(guān)中林家村、陜南石泉斜率最小，降雨下降最明顯。從圖2以及表1的潛在蒸散發(fā)變化趨勢(shì)中可以看出：陜北除了丁家溝外，其他流域的潛在蒸散發(fā)均呈現(xiàn)不顯著下降的趨勢(shì)，在關(guān)中林家村是不顯著下降，其余是不顯著上升，在陜南除竹林關(guān)外，均為不顯著下降。并且不同區(qū)域徑流變化的斜率不等，其中陜北趙石窯下降最明顯，關(guān)中甘谷驛上升最明顯，陜南白河下降最明顯。

圖2 徑流、降水量及潛在蒸散發(fā)的變化趨勢(shì)

3.2 各小流域徑流突變及要素變化分析

根據(jù)公式(4)，累積距平的方法對(duì)各小流域的徑流進(jìn)行分析，得到各小流域在研究時(shí)段徑流發(fā)生突變的位置(表1)。

表1 秦嶺南北14個(gè)小流域水文—?dú)夂?/p>

注：“***”表示顯著變化；“ns”表示不顯著變化。

從表1各小流域徑流變化的突變點(diǎn)中可以看出：不同的研究區(qū)徑流突變的年份存在差異，對(duì)于研究時(shí)間從1961年開(kāi)始的區(qū)域，突變點(diǎn)主要集中在20世紀(jì)70年代末80年代初，狀頭和甘谷驛出現(xiàn)了異常，其突變點(diǎn)出現(xiàn)在20世紀(jì)90年代中期；對(duì)于研究時(shí)間從1980年開(kāi)始的流域，突變點(diǎn)主要集中在20世紀(jì)80年代末90年代初期的1989年和1990年。

此外，不同區(qū)域各要素的變化情況及其變化程度也是有差異的，從圖3-a中可以發(fā)現(xiàn)：屬于陜南的5個(gè)小流域它們的徑流與基準(zhǔn)期相比較下降最明顯，徑流減少最大的流域?yàn)槭跔铑^、甘谷驛以及綏德徑流減少量比較小。從圖3-b降雨變化圖中可以看出，降雨的變化情況和變化程度與徑流的變化比較一致，屬于陜南的5個(gè)小流域降雨減少量最大，降雨減少最多的為石泉，而在甘谷驛、綏德、丁家溝降雨的減少量比較小。從圖3-c潛在蒸散發(fā)變化圖中可以看出，潛在蒸散發(fā)的變化情況與徑流、降雨的變化存在差異。石泉、白河以及橫山的潛在蒸散發(fā)減少最明顯，而溫家川、竹林關(guān)以及荊紫關(guān)的變化比較小。

3.3 氣候?qū)搅鞯拿舾邢禂?shù)變化

根據(jù)公式(10)和(11)可以準(zhǔn)確計(jì)算出各個(gè)流域的下墊面參數(shù)ω、降雨對(duì)徑流變化的敏感系數(shù)β以及潛在蒸散發(fā)對(duì)徑流變化的敏感系數(shù)γ。不同的研究區(qū)由于地形、氣候、水熱條件等因素的不同導(dǎo)致不同的ω、β和γ，如圖4所示。

從圖4-a中可以發(fā)現(xiàn)：陜北以及關(guān)中的9個(gè)小流域,除了突變點(diǎn)狀頭和趙石窯外其他下墊面參數(shù)ω從北到南在增大，其中北部區(qū)域的ω值小于南部區(qū)域，北部的ω值最低為0.09，出現(xiàn)在高家川，而偏南部的ω值最高可達(dá)14.09，出現(xiàn)在狀頭。在陜南各流域的ω值普遍偏小。從圖3-b降雨對(duì)徑流變化的敏感系數(shù)β可以看出，降雨對(duì)徑流變化的β各流域存在差異，在陜北從北到南，β以不等速率減小，在關(guān)中地區(qū)增大，在陜北和關(guān)中β最大值為0.45出現(xiàn)在高家川，β最大值為0.05出現(xiàn)在趙石窯。在陜南β普遍較大，最大的β為0.75，出現(xiàn)在石泉。從圖3-c潛在蒸散發(fā)對(duì)徑流變化的敏感系數(shù)γ可以看出，部分區(qū)域的γ的變化情況和變化程度與下墊面參數(shù)ω的變化接近一致，其中在陜北關(guān)中，最大值出現(xiàn)在趙石窯和狀頭，為-0.01,在陜南γ的取值比較接近最大值出現(xiàn)在白河，為-0.19，最小值出現(xiàn)在石泉，為-0.28。

圖3 徑流、降雨以及潛在蒸散發(fā)的變化

圖4 下墊面參數(shù)，降雨以及潛在蒸散對(duì)徑流變化的敏感系數(shù)

4 結(jié)論及討論

4.1 討論

4.1.1 氣候變化和人類活動(dòng)的貢獻(xiàn)率 徑流的變化主要?dú)w因于2個(gè)方面：氣候變化和人類活動(dòng)。在敏感性分析中，通過(guò)公式(8)和(9)求得徑流對(duì)降雨及蒸發(fā)的敏感系數(shù)，計(jì)算由氣候因素變化導(dǎo)致的徑流的變化量，從而分離了氣候變化和人類活動(dòng)對(duì)徑流的影響。研究區(qū)14個(gè)流域氣候變化與人類活動(dòng)對(duì)徑流影響分離的結(jié)果如表2所示。

從表2中可以看出：在陜北人類活動(dòng)對(duì)徑流變化的貢獻(xiàn)率明顯高于氣候變化，平均人類活動(dòng)貢獻(xiàn)率為77.47%。除過(guò)狀頭這個(gè)突變區(qū)域外，人類活動(dòng)對(duì)徑流變化的最低貢獻(xiàn)率也達(dá)到54.6%(甘谷驛)，超過(guò)一半，各個(gè)小流域徑流的變化主要受人類活動(dòng)控制，越是偏北，人類活動(dòng)在徑流變化中占的地位越高，人類活動(dòng)對(duì)徑流變化的影響越大；并且觀察陜北關(guān)中的9個(gè)小流域的研究結(jié)果發(fā)現(xiàn)，根據(jù)地理位置從北到南，人類活動(dòng)對(duì)徑流變化的貢獻(xiàn)率呈現(xiàn)逐漸降低趨勢(shì)，偏北部流域人類活動(dòng)對(duì)徑流變化的貢獻(xiàn)率普遍高于偏南部區(qū)域。也就是說(shuō)，從北到南，人類活動(dòng)對(duì)徑流變化的影響逐漸減弱，氣候因素占的份量變大；關(guān)中地區(qū)人類活動(dòng)所占比份有所下降，平均人類活動(dòng)貢獻(xiàn)率為50.34%；在表2中，陜南的5個(gè)小流域發(fā)現(xiàn)人類活動(dòng)對(duì)徑流變化的貢獻(xiàn)率除竹林關(guān)的42.69%外，其余均達(dá)到50%以上，平均人類活動(dòng)貢獻(xiàn)率為58.3%，也就是說(shuō)，陜南的這些小流域，人類活動(dòng)對(duì)徑流變化的影響亦占有很重要的份量。

表2 秦嶺南北14個(gè)小流域水文-氣候要素變化及貢獻(xiàn)率分析

4.1.2 陜西水資源變化特征及對(duì)策 (1)在陜北人類活動(dòng)對(duì)徑流變化占主導(dǎo)作用，是徑流變化的主要原因，人類活動(dòng)在此處是多樣的，主要是為了控制泥沙修建了大量的水土保持工程。隨著國(guó)家黃土高原治理規(guī)劃(2013～2020年)、陜西省水保工程的實(shí)施以及林草覆蓋的加強(qiáng)，導(dǎo)致徑流可能將進(jìn)一步銳減，人類活動(dòng)對(duì)徑流的影響可能會(huì)逐步加強(qiáng)，在這樣的背景下，必須依靠調(diào)水解決區(qū)域的水資源問(wèn)題。

(2)渭河林家村徑流的減少，人類活動(dòng)占主導(dǎo)地位。由于本區(qū)域是陜西省的糧食主產(chǎn)區(qū)和經(jīng)濟(jì)社會(huì)中心，所以需要加強(qiáng)對(duì)渭河主要徑流區(qū)的甘肅省天水市提供生態(tài)補(bǔ)償，從而確保流入陜西的水量及水質(zhì)。

(3)陜南的徑流變化主要是氣候變化和人類活動(dòng)共同帶來(lái)的。其中下墊面特征所導(dǎo)致的氣候變化所占比重同樣較大，因此需要通過(guò)調(diào)整人類活動(dòng)來(lái)更好地適應(yīng)氣候變化。具體的人類活動(dòng)包括增加林草覆蓋、保護(hù)環(huán)境降低溫室效應(yīng)、減少耗水工業(yè)以及高污染的產(chǎn)業(yè)等。

由于研究資料等原因，未對(duì)區(qū)域人類活動(dòng)進(jìn)一步細(xì)化并且研究區(qū)各流域研究時(shí)間不一致；此外，本文的研究將氣候變化與人類活動(dòng)分離，事實(shí)上，兩者是互相聯(lián)系的，如何衡量2個(gè)要素的相互聯(lián)系以及這種聯(lián)系對(duì)研究結(jié)果的影響均需要進(jìn)一步深入研究。

4.2 結(jié)論

(1)各流域降雨和徑流都呈現(xiàn)不同程度的下降趨勢(shì)，其中徑流呈現(xiàn)顯著下降的趨勢(shì)，降雨量呈現(xiàn)不顯著下降的趨勢(shì)。潛在蒸散發(fā)在丁家溝、甘谷驛呈現(xiàn)不顯著上升的趨勢(shì)，其他區(qū)域均呈現(xiàn)不顯著下降的趨勢(shì)。不同流域徑流突變點(diǎn)出現(xiàn)的時(shí)間不一致，主要集中在1980年和1990年前后。

(2)在陜北它們的下墊面參數(shù)ω從北到南以不等的速率增大，偏北部區(qū)域的ω小于偏南部區(qū)域。在關(guān)中，下墊面參數(shù)達(dá)到最大；在陜南各個(gè)流域的下墊面參數(shù)ω普遍偏小。在陜北β以不等速率減小之后在關(guān)中又增大,在陜南β普遍較大。γ從北到南逐漸增大，陜南γ普遍較大。

(3)為了控制泥沙修建了大量的水土保持工程，國(guó)家黃土高原治理規(guī)劃(2013～2020年)、陜西省水保工程的實(shí)施以及林草覆蓋的加強(qiáng)使得在陜北人類活動(dòng)對(duì)徑流變化的貢獻(xiàn)率明顯高于氣候變化，各小流域徑流的變化主要受人類活動(dòng)控制，越是偏北，人類活動(dòng)在徑流變化中占的地位越高，人類活動(dòng)對(duì)徑流變化的影響越大，平均貢獻(xiàn)率為77.47%；關(guān)中地區(qū)由于主要受工農(nóng)業(yè)以及生活用水等人類活動(dòng)的影響，其人類活動(dòng)所占比分有所下降，平均貢獻(xiàn)率為50.34%；在陜南發(fā)現(xiàn)由于其下墊面的特征人類活動(dòng)對(duì)徑流變化的影響亦占有很重要的分量，平均貢獻(xiàn)率為58.3%。

[1] Nemec J, Schaake J. Sensitivity of water resource systems to climate variation[J]. Hydrologic Science Journal, 1982, 27(2): 327-343.

[2] Pikounism, Varanoue, Baltase, et al. Application of the SWAT model in the pinios river basin under different land-use scenarios[J]. Global Nest: The Int J, 2003, 5(2): 71-79.

[3] 李棟梁,張佳麗,全建瑞，等.黃河上游徑流量演變特征及成因研究[J].水科學(xué)進(jìn)展,1998,9(1):22-28.

[4] 張建云,章四龍,王金星，等.近50年來(lái)中國(guó)六大流域年際徑流變化趨勢(shì)研究[J].水科學(xué)進(jìn)展，2007，18(2)：230-234.

[5] 張強(qiáng),張存杰,白虎志，等.西北地區(qū)氣候變化新動(dòng)態(tài)及對(duì)干旱環(huán)境的影響——總體暖干化,局部出現(xiàn)暖濕跡象[J].干旱氣象，2009，28(1):1-7.

[6] Chen J F, Li X B, Zhang M. Simulating the impacts of climate variation and land-cover changes on basin hydrology:A case study of the Suomo basin[J]. Science in China (Series D: Earth Sciences), 2005, 48(9): 1501-1508.

[7] 李志,劉文兆,鄭粉莉，等.黃土塬區(qū)氣候變化和人類活動(dòng)對(duì)徑流的影響[J].生態(tài)學(xué)報(bào),2010,30(9):2379-2386.

[8] Budyko M I. Climate and Life[M]. New York: Academic Press, 1974: 508.

[9] Dooge J C I, Bruen M, Parmentier B. A simple model for estimating the sensitivity of runoff to long-term changes in precipitation without a change in vegetation[J]. Advance in Water Resources, 1999, 23(2): 153-163.

[10] Milly P C D, Dunne K A. Macroscale water fluxes 2. Water and energy supply control of their inter annual variability[J]. Water Resource Research, 2002, 38(10): 241-249.

[11] 胡珊珊,鄭紅星,劉昌明，等.氣候變化和人類活動(dòng)對(duì)白洋淀上游水源區(qū)徑流的影響[J].地理學(xué)報(bào),2012,67(1)：62-70.

[12] 曹小雪，黃建武，揭毅.近52年來(lái)武漢市氣候變化特征分析[J].江西農(nóng)業(yè)學(xué)報(bào)2014，26(9)：80-85.

[13] 許炯心,孫季.嘉陵江流域年徑流量的變化及其原因[J].山地學(xué)報(bào),2007,25(2):153-159.

[14] 王雁,丁永建,葉柏生，等.黃河與長(zhǎng)江流域水資源變化原因[J].中國(guó)科學(xué):地球科學(xué),2013(7)：14.

[15] Allen R G, Pereira L S, Raes D, et al. Crop evapotranspiration-Guidelines for computing crop water requirements[M]. Rome: FAO，1998: 4-7.

[16] Liu C M, Zheng H. Changes in components of the hydrological cycle in the Yellow River basin during the second half of the 20th century[J]. Hydrological Process, 2004, 18(12): 2337-2345.

[17] Yu P S, Yang T C, Wu C K. Impact of climate change on water resource in southern Taiwan[J]. Journal of Hydrology, 2011, 22(5): 161-175.

[18] 劉昌明,鄭紅星.黃河流域水循環(huán)要素變化趨勢(shì)分析[J].自然資源學(xué)報(bào),2003,18(2):129-135.

[19] Li L J, Zhang L, Wang H, et al. Assessing the impact of climate variability and human activities on streamflow from the Wuding River basin in China[J]. Hydrological Process, 2007, 21(25): 3485-3491.

[20] Liu Q, Yang Z, Cui B, et al. Temporal trends of hydro-climatic variables and runoff response to climatic variability and vegetation changes in the Yiluo River basin, China[J]. Hydrological Processes, 2009, 23(21): 3030-3039.

[21] Shao Q, Traylen A, Zhang L. Nonparametric method for estimating the effects of climatic and catchment characteristics on mean annual evapotranspiration[J]. Water Resources Research, 2012, 48(3): 949-960.

[22] Zhang L, Dawes W R, Walker G R. Response of mean annual ET to vegetation changes at catchment scale[J]. Water Resources Research, 2001, 37(3): 701-708.

(責(zé)任編輯：管珊紅)

A Preliminary Study on Impact of Climatic Change and Human Activity on Water Resources: Taking Shaanxi Province as An Example

HE Yan-zi, YUE Da-peng*, DA Xing, CHENG Jin-wen, ZHANG Jun-ru

(College of Tourism and Environmental Sciences, Shaanxi Normal University, Xi’an 710119, China)

According to the status of runoff, rainfall and potential evaporation in 14 small watersheds of Shaanxi province during 1961～2009, the author analyzed the variation trend and abrupt-change points of the runoff by using the methods of Mann-Kendall nonparametric test and cumulative departure, separated the impacts of human activity and climatic change on the runoff by using the sensitive coefficient method, and studied the changes in omega, beta, and gamma values of various watersheds based on the Budyko assumption and by adopting Zhang’s curve. In the studied area, the runoff showed a significant reduction tendency, rainfall revealed an un-significant reduction trend; the changes in potential evaporation in different regions were inconsistent. The abrupt-change points of runoff in different watersheds were different, and mainly concentrated in 1980 and 1990 or so. In Northern Shaanxi and Guanzhong areas, the omega value increased gradually from the north to the south. The omega value was lower in Southern Shaanxi. Both the beta value and the gamma value were generally higher in Southern Shaanxi. In Northern Shaanxi, due to the implementation of such measures as water conservation, returning farmland to forestry, and loess management, human activity had a higher contribution rate to the runoff changes than climatic change. In Guanzhong area, due to the water use for industrial and agricultural sectors, and the domestic water, the contribution rate of human activity to the runoff changes was lower. In Southern Shaanxi, the contribution rate of human activity to the runoff changes also occupied a large proportion due to the nature of the underlying surface.

Shaanxi; Runoff change; Mann-Kendall; Accumulative departure method; Sensitive coefficient method

2015-07-09

國(guó)家自然科學(xué)基金項(xiàng)目(41071012)；陜西省水土流失補(bǔ)償費(fèi)資助項(xiàng)目“陜西省能源開(kāi)發(fā)項(xiàng)目水土保持監(jiān)測(cè)技術(shù)研究及信息系統(tǒng)研究”。

賀燕子(1992—)，女，甘肅寧縣人，碩士研究生，研究方向：水土資源評(píng)價(jià)與規(guī)劃。*通訊作者：岳大鵬。

K928.42

A

1001-8581(2016)12-0087-07