黃土高原涇河流域氣候和土地利用變化對徑流產(chǎn)沙的影響

鄭培龍, 李云霞, 趙 陽, 張曉明, 寇馨月, 徐曉捂

(1.北京水保生態(tài)工程咨詢有限公司, 北京 100053; 2.中水北方勘測設(shè)計(jì)研究有限責(zé)任公司,北京 300222; 3.中國水利水電科學(xué)研究院, 北京 100048; 4.北京林業(yè)大學(xué), 北京 100083)

?

黃土高原涇河流域氣候和土地利用變化對徑流產(chǎn)沙的影響

鄭培龍1, 李云霞2, 趙 陽3, 張曉明3, 寇馨月4, 徐曉捂4

(1.北京水保生態(tài)工程咨詢有限公司, 北京 100053; 2.中水北方勘測設(shè)計(jì)研究有限責(zé)任公司,北京 300222; 3.中國水利水電科學(xué)研究院, 北京 100048; 4.北京林業(yè)大學(xué), 北京 100083)

研究黃土高原流域水沙對氣候和土地利用變化的響應(yīng)機(jī)理，對生態(tài)調(diào)控決策制定及區(qū)域水土保持措施優(yōu)化配置具有重要意義。以黃土高原涇河流域?yàn)檠芯繉ο?，采用Mann-Kendall檢驗(yàn)、累積距平和雙累積曲線等方法，定量分析了流域近40年降水、徑流和輸沙的演變規(guī)律及土地利用變化特征，探討了流域氣候和土地利用變化對流域水沙量的影響。結(jié)果表明：(1) 涇河流域徑流深和輸沙模數(shù)年際均呈顯著下降趨勢，在20世紀(jì)80，90年代尤為明顯；(2) 涇河流域土地利用結(jié)構(gòu)發(fā)生了重大變化，具體表現(xiàn)為坡耕地向梯田的轉(zhuǎn)換、耕地向林地草地的轉(zhuǎn)換和沙地裸地向草地灌叢的轉(zhuǎn)換；(3) 人類活動(dòng)影響下的土地利用變化對徑流泥沙貢獻(xiàn)率分別為82.90%，94.10%，氣候因素影響很小。

徑流泥沙; 土地利用變化; 突變; 涇河流域

近幾十年來，在氣候和土地利用變化的共同影響下，黃河流域的水沙量發(fā)生了明顯變化[1]。但徑流量、輸沙量和含沙量變化不一致，徑流量和輸沙量顯著降低，含沙量卻變化不明顯，汛期含沙量甚至還偏高[2]。目前，我國各流域水沙動(dòng)態(tài)變化對氣候和土地利用變化響應(yīng)的研究工作已成為環(huán)境、地理和生態(tài)等科學(xué)的研究熱點(diǎn)。王光謙等[3]認(rèn)為人類活動(dòng)是黃河流域土地利用和覆蓋變化的主要驅(qū)動(dòng)力，程磊等[4]對渭河流域20年間的土地利用時(shí)空變化進(jìn)行研究，同樣得出人類活動(dòng)是渭河流域土地利用變化的主要驅(qū)動(dòng)力的結(jié)論。但由于不同流域的研究尺度和研究方法存在相當(dāng)大的差別，從而造成了研究結(jié)果的不同。劉昌明等[5]研究認(rèn)為區(qū)域的降水變化是引起黃河上游流域徑流變化的主要原因，而黃河中游流域則人類活動(dòng)影響更為顯著，因此兩者對水沙動(dòng)態(tài)變化的影響程度是因地而異的。因此，由以上研究結(jié)果可知，土地覆被變化是一個(gè)復(fù)雜的過程，需要將人為和自然因素合并討論。

本文采用累積距平分析、Mann-Kendall趨勢檢驗(yàn)和雙累積曲線等方法，定量分析涇河流域近50年降水、徑流和輸沙的演變規(guī)律及研究區(qū)土地利用變化特征，探討流域氣候和土地利用變化對流域水沙動(dòng)態(tài)變化的影響。旨在揭示涇河流域水沙動(dòng)態(tài)變化對土地利用變化的響應(yīng)機(jī)理，為黃土高原地區(qū)土地利用生態(tài)調(diào)控、水土保持措施的優(yōu)化配置提供科學(xué)依據(jù)。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況

涇河是渭河的第一大支流，水系發(fā)達(dá)，干流長約455 km，流經(jīng)面積45 421 km2。流域地處黃土高原中部，地勢西北高東南低，地貌類型主要有黃土高原溝壑、丘陵溝壑、丘陵林、階地和土石丘陵。屬大陸季風(fēng)性氣候，各區(qū)域四季溫差明顯，年平均氣溫9℃左右。流域是半濕潤半干旱過渡區(qū)，徑流主要來自降雨。年平均降水量550 mm，降水量由南向北遞減，年際變化大，夏季量大多暴雨，冬季干旱少雨有霜或降雪。主要土壤類型為黃土和灰褐土，土質(zhì)疏松，易水土流失。地植被為溫帶森林草原過渡類型，分布不均勻。流域西北部多為退化荒草地，草地質(zhì)量較低。東南部多為山地植被，類型多樣且覆蓋率較高。

1.2 數(shù)據(jù)處理

1.2.1 數(shù)據(jù)采集 涇河流域徑流泥沙數(shù)據(jù)(1971—2010年)來自涇河流域出口水文控制站——張家山水文站；降水?dāng)?shù)據(jù)(1971—2010年)來自于流域內(nèi)的11個(gè)水文站點(diǎn)：洪德站、賈橋站、慶陽站、雨落坪站、袁家庵站、洪河站、巴家嘴站、毛家河站、涇川站、蘆村河站和張家山站，考慮到流域水文站點(diǎn)較少，流域面降水量采用泰森多邊形法計(jì)算求得。涇河流域土地利用數(shù)據(jù)主要來源于區(qū)域?qū)ｎ}研究資料、流域所述行政區(qū)社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)年鑒和相關(guān)出版文獻(xiàn)資料[6-10]，研究時(shí)段主要包括：1986年、1995年和2000年3個(gè)時(shí)段。

1.2.2 研究方法

(1) Mann-Kendall檢驗(yàn)。Mann-Kendall突變檢驗(yàn)法[11]是檢驗(yàn)氣候突變，不易受少數(shù)異常值干擾的方法。具有n個(gè)樣本量的時(shí)間序列x，構(gòu)造一個(gè)秩序列Sk，它是第i時(shí)刻數(shù)值大于j時(shí)刻數(shù)值個(gè)數(shù)的累計(jì)數(shù)。

(1)

(2)

在時(shí)間序列隨機(jī)獨(dú)立的假設(shè)下，定義統(tǒng)計(jì)量：

(3)

式中：UF1=0；E(sk)，var(sk)——累計(jì)數(shù)sk的均值和方差，在x1，x2，…，xn相互獨(dú)立，且有相同連續(xù)分布時(shí)，通過下式計(jì)算：

(4)

UFk為標(biāo)準(zhǔn)正態(tài)分布，是按時(shí)間序列x順序計(jì)算出的統(tǒng)計(jì)量序列，按時(shí)間序列x逆序xn，xn-1，…，x1，再重復(fù)上述過程，同時(shí)使得UBk=-UFk(k=n，n-1，…，1)，UB1=0。繪制UFk和UBk序列曲線，若兩曲線相交且交點(diǎn)在置信區(qū)間之內(nèi)，則此交點(diǎn)可認(rèn)為是突變點(diǎn)。

(2) 累積距平曲線。累積距平法可顯現(xiàn)水文要素年際變化的階段性特征。該方法先計(jì)算每年的水文要素距平，然后按年序累加，得到距平累積序列。采用累積距平法判斷徑流量和輸沙量的序列突變年份。該方法是由曲線直觀判斷離散數(shù)據(jù)點(diǎn)變化趨勢的一種非線性統(tǒng)計(jì)方法[12]。

時(shí)間序列X(x1，x2，…，xn)，其在某一時(shí)刻t的累積距平為：

(5)

根據(jù)距平有正有負(fù)的特點(diǎn)，當(dāng)距平累積持續(xù)增大時(shí)，表明該時(shí)段內(nèi)水文要素值大于平均值；當(dāng)距平累積持續(xù)不變，表明該時(shí)段水文要素值等于平均值；當(dāng)距平累積持續(xù)減小時(shí)，表明時(shí)段內(nèi)水文要素值小于平均值。

(3) 雙累積曲線法。雙累積曲線法是進(jìn)行時(shí)間序列分析的一種常用方法[13]，它的基本原理是兩個(gè)變量按同一時(shí)間長度逐步累加，一個(gè)變量作為橫坐標(biāo)，另一個(gè)變量作為縱坐標(biāo)，降水與徑流的雙累積曲線是可以揭示徑流是否有趨勢性變化的一種常用方法。對于正常的水文年份，如果不受外界影響，徑流量累積值雖有波動(dòng)，但沒有系統(tǒng)偏離，降水與徑流的雙累積曲線應(yīng)該是一條直線；如果流域下墊面發(fā)生了改變，則降水與徑流雙累積曲線會(huì)發(fā)生明顯的轉(zhuǎn)折。

(4) 以氣候和土地利用變化對泥沙影響分析為例，本研究以降雨—輸沙量雙累積曲線為依據(jù)，以曲線的轉(zhuǎn)折點(diǎn)為劃分的基準(zhǔn)期和措施期，以此來評估人類活動(dòng)下的土地利用變化對輸沙量的影響。

經(jīng)回歸分析求得累積年降雨量與累積年輸沙量的線性相關(guān)方程，將流域治理后不同時(shí)段的累積年降雨量值分別代入方程中，求得累積年輸沙量，然后與同期實(shí)測值比較，即得減沙效應(yīng)。計(jì)算公式如下:

(6)

式中：R1——利用累積年降雨量—累積年輸沙量的線性相關(guān)方程計(jì)算出的年輸沙量；R2——流域治理后的實(shí)測年輸沙量。

通過建立基準(zhǔn)期降雨—輸沙關(guān)系，用以計(jì)算措施期預(yù)期自然輸沙量，將其與基準(zhǔn)期平均輸沙量相比得到降雨變化對輸沙的影響量，將其與實(shí)際觀測輸沙量對比，得到人類活動(dòng)影響量[14]：

ΔS=QS2-QS1=ΔS1+ΔS2

(7)

(8)

(9)

2 結(jié)果與分析

2.1 流域降水及水沙過程年際變化規(guī)律

由圖1可知，涇河流域降水呈現(xiàn)出下降趨勢。從5年滑動(dòng)平均趨勢線可看出涇河流域20世紀(jì)90年代前降水量無明顯的趨勢性變化，但1991—2001年降水量下降顯著，2003年之后開始回升。涇河平均年徑流深較小，其中最大和最小徑流深分別為51.03，13.96 mm。從線性趨勢線的走勢來看，涇河年徑流深逐年下降。涇河的最大輸沙模數(shù)發(fā)生在1973年，最小輸沙模數(shù)發(fā)生在1972年，均為早期的研究時(shí)期，且輸沙模數(shù)的年際變化很大。從5年滑動(dòng)平均趨勢線可看出，雖然年際間存在波動(dòng)變化，但從70年代開始輸沙模數(shù)有了明顯下降。泥沙量的減少主要是由于流域徑流近些年急劇下降造成的。涇河的水沙動(dòng)態(tài)變化基本同步。

2.2 流域水沙過程突變分析

由圖2可知，徑流深的年際變化特征十分明顯，呈現(xiàn)先增大后減小的趨勢。涇河流域突變點(diǎn)發(fā)生在1985年，發(fā)生第一次突變前累積距平值呈增加趨勢，此階段徑流數(shù)據(jù)大于平均值。之后在1996年出現(xiàn)突變拐點(diǎn)，80—90年代中期累積距平值由上升變?yōu)榧眲∠陆底兓f明進(jìn)入徑流偏小的階段，這與90年代降水的減少有明顯關(guān)系。輸沙模數(shù)累積矩平變化與徑流的變化基本相一致。涇河流域輸沙模數(shù)在1996年發(fā)生第一次突變，其前后累積距平值為顯著上升和顯著下降的趨勢，說明輸沙模數(shù)由偏高階段轉(zhuǎn)為偏低階段，這與徑流的變化有明顯相關(guān)性。

圖1 涇河水文要素變化趨勢

圖2 涇河流域的徑流量和輸沙量的突變特征

2.3 流域土地利用變化

從表1可以看出，涇河流域1986—1995年林地、草地和水域3種土地利用類型面積略有減少，相應(yīng)的耕地、城鎮(zhèn)用地和未利用地面積略有增加。然而，1995—2000年林地、草地和水域逐漸恢復(fù)接近原來的面積，耕地和未利用地減少，城鎮(zhèn)用地繼續(xù)增加?？梢?，涇河流域土地利用類型在兩個(gè)時(shí)期內(nèi)發(fā)生著由生態(tài)型向人工型轉(zhuǎn)化，再由人工型恢復(fù)到生態(tài)型的循環(huán)。

表1 涇河流域土地利用面積變化

2.4 土地利用變化對水沙過程的定量影響

降雨—徑流量、降雨量—輸沙量的雙累積曲線的斜率代表單位降水所產(chǎn)生的徑流量和輸沙量。在自然狀態(tài)下，降雨與徑流輸沙量的關(guān)系應(yīng)該是一定的。如果流域下墊面發(fā)生了改變，則其雙累積曲線斜率會(huì)發(fā)生明顯的轉(zhuǎn)折。當(dāng)曲線發(fā)生轉(zhuǎn)折時(shí)，說明下墊面受人類活動(dòng)影響改變，從而影響了徑流輸沙。從圖3—4可以看出，在各流域降水—徑流、降雨量—輸沙量的雙累積曲線圖中存在轉(zhuǎn)折點(diǎn)。由此可見流域徑流對降水變化響應(yīng)劃分為基準(zhǔn)期和變化期。第一次突變前的時(shí)期，相對于突變后的時(shí)期可看作為徑流輸沙的自然基準(zhǔn)期，流域下墊面產(chǎn)流產(chǎn)沙受降雨外的其他條件影響極弱。我們可以發(fā)現(xiàn)各流域的變化期降水—徑流、降雨量—輸沙量的雙累積曲線圖的斜率均較上一時(shí)段減小，說明其為顯著下降階段，下墊面受人類活動(dòng)影響顯著。并且，各流域受影響程度有差異，各階段擬合方程的判定系數(shù)大部分都大于0.95，擬合度高。

圖3 不同時(shí)段降雨量-徑流量雙累積曲線

根據(jù)研究流域降水—徑流、降水—輸沙的雙累積曲線在研究系列的變化趨勢及突變點(diǎn)，對流域徑流輸沙研究時(shí)間序列進(jìn)行了階段劃分。并分別計(jì)算了研究流域氣候變化和土地利用變化對流域徑流輸沙的定量影響。

由表2—3可知，涇河流域人類活動(dòng)因素對流域徑流影響貢獻(xiàn)率為82.90%，是氣候貢獻(xiàn)率的4.85倍；涇河流域人類活動(dòng)因素對流域輸沙量影響貢獻(xiàn)率為94.10%，氣候貢獻(xiàn)率僅為5.90%。人類活動(dòng)因素是研究流域徑流發(fā)生變異的主要驅(qū)動(dòng)因素，土地利用變化是主要的人類活動(dòng)表現(xiàn)。雖然氣候?qū)α饔驈搅髯兓a(chǎn)生影響，但影響貢獻(xiàn)率均較小。

圖4 不同時(shí)間段降雨量-輸沙量雙累積曲線

表2 線性回歸方程擬合

3 結(jié)論與討論

(1) 涇河流域降水、徑流輸沙均呈減少趨勢，且徑流輸沙減小趨勢更為明顯。徑流突變時(shí)間為1985年和1996年，與降雨變化有關(guān)；輸沙突變時(shí)間為1996年。

(2) 涇河流域主要土地利用類型為林地、草地和耕地。主要變化是坡耕地向梯田的轉(zhuǎn)換、耕地向林地、草地的轉(zhuǎn)換和沙地、裸地向草地、灌叢的轉(zhuǎn)換。

(3) 土地利用變化是涇河流域徑流輸沙發(fā)生變異的主要驅(qū)動(dòng)因素，貢獻(xiàn)率均大于80%；而降水影響貢獻(xiàn)率則較小。但較之人類活動(dòng)對徑流影響而言，人類活動(dòng)對流域輸沙量的影響較大。

表3 涇河流域徑流輸沙量對氣候和土地利用變化的響應(yīng)

引起了區(qū)域氣候變化的因子主要有降水、溫度、輻射、濕度、風(fēng)速等，它們綜合作用于氣候影響下的環(huán)境。但由于數(shù)據(jù)采集獲取等有限，時(shí)間序列長度的溫度影響不夠強(qiáng)烈，本研究氣候變化主要指降水變化，這對研究結(jié)果的準(zhǔn)確性有一定影響，條件允許的情況下應(yīng)當(dāng)加強(qiáng)影響生態(tài)水文過程的其他氣候因子的分析。

[1] Liu C, Sui J, Wang Z Y. Changes in runoff and sediment yield along the Yellow River during the period from 1950 to 2006[J]. Journal of Environmental Informatics,2008,12(2):129-139.

[2] 張學(xué)成,王玲,高貴成.黃河泥沙動(dòng)態(tài)分析[J].泥沙研究,2002(1):43-52.

[3] 王光謙,王思遠(yuǎn),陳志祥.黃河流域的土地利用和土地覆蓋變化[J].清華大學(xué)學(xué)報(bào):自然科學(xué)版,2004,44(9):1218-1222.

[4] 程磊,徐宗學(xué),羅睿,等.渭河流域1980—2000年LUCC時(shí)空變化特征及其驅(qū)動(dòng)力分析[J].水土保持研究,2009,16(5):1-6,30.

[5] 劉昌明,張學(xué)成.黃河干流實(shí)際來水量不斷減少的成因分析[J].地理學(xué)報(bào),2004,59(3):323-330.

[6] 王小軍,蔡煥杰,張鑫,等.皇甫川流域水沙變化特點(diǎn)及其趨勢分析[J].水土保持研究,2009，16(1):222-226.

[7] 李慶云,呂素冰,劉蕾,等.黃土高原藉河流域近35 a來土地利用/覆被變化過程[J].華北水利水電學(xué)院學(xué)報(bào),2013,34(6):32-35.

[8] 吳琳娜,楊勝天,劉曉燕,等.1976年以來北洛河流域土地利用變化對人類活動(dòng)程度的響應(yīng)[J].地理學(xué)報(bào),2014,69(1):54-63.

[9] 楊麗,謝高地,甄霖,等.涇河流域土地利用格局的時(shí)空變化分析[J].資源科學(xué),2005,27(4):26-32.

[10] 張程,李小兵,張立,等.皇甫川流域土地利用/覆被變化對生態(tài)服務(wù)價(jià)值的影響[J].北京師范大學(xué)學(xué)報(bào):自然科學(xué)版,2009,45(4):399-403.

[11] 符淙斌,王強(qiáng).氣候突變的定義和監(jiān)測方法[J].大氣科學(xué),1992,12(4):482-493.

[12] 趙廣舉,穆興民,溫仲明,等.皇甫川流域降水和人類活動(dòng)對水沙變化的定量分析[J].中國水土保持科學(xué),2013,11(4):1-8.

[13] 穆興民,巴桑赤烈,Zhang Lu,等.黃河河口鎮(zhèn)至龍門區(qū)間來水來沙變化及其對水利水保措施的響應(yīng)[J].泥沙研究, 2007,52(2):36-41.

[14] 王國慶.氣候變化對黃河中游水文水資源影響的關(guān)鍵問題研究[D].南京:河海大學(xué),2006.

Effects of Climate Variation and Land Use Change on Runoff in Jinghe Basin of the Loess Plateau

ZHENG Peilong1, LI Yunxia2, ZHAO Yang3, ZHANG Xiaoming3, KOU Xinyue4, XU Xiaowu4

(1.BeijingWaterConservationEcologicalEngineeringConsultingCo.,Ltd.,Beijing100053,China; 2.NorthernInvestigationofChinaWaterResources,DesignandResearchCo.,Ltd.,TianJin300222,China; 3.ChinaInstituteofWaterResourcesandHydropowerResearch,Beijing100048,China; 4.BeijingForestryUniversity,Beijing100083,China)

In the Loess Plateau, study on the response mechanism of runoff and sediment to land use change is of great importance for ecological regulation decision-making and the optimal allocation of soil and water conservation measures. Therefore, based on Jinghe Basin in the Loess Plateau as the research example, the Mann-Kendall test, accumulative-annual anomaly and double mass analysis were used to quantitatively examine the precipitation, runoff and sediment change and the characteristics of land use change in nearly 40 years, then the influences of the climate and land use change on runoff and sediment were discussed. Results showed that: (1) trends of annual precipitation of Jinghe Basin were different and not obvious; annual runoff depth and sediment transport modulus of the five watersheds showed the significant decline, especially in the 1980s—1990s; (2) it turned out that land use structure greatly changed in recent decades; major changes were the conversions of slope land to terrace, cultivated land to forestland and grassland, sandy bare land to grassland and thickets; (3) the influence of human activities on runoff and sediment contribution rates were 82.90% and 94.10%, respectively, the impacts of climate factors were very slight.

runoff and sediment; land use change; extreme characteristic; Jinghe Basin

2014-09-03

2014-10-08

國家自然科學(xué)基金(51379008，41501041);流域水循環(huán)模擬與調(diào)控國家重點(diǎn)實(shí)驗(yàn)室自主研究課題(泥ZY1304,2014QN04)

鄭培龍(1976—),男,山西盂縣人,碩士,工程師,研究方向?yàn)樗帘３?。E-mail:zhengpeilong1122@163.com

趙陽(1986—),男,河北棗強(qiáng)人,博士,工程師,主要從事水土保持研究。E-mail:zhaoyang1224@163.com

S157; P333

1005-3409(2015)05-0020-05