黃土高原典型流域水沙變化歸因?qū)Ρ确治?/h1>
2020-04-10 07:25蔣凱鑫于坤霞曹文洪張曉明
農(nóng)業(yè)工程學(xué)報(bào) 2020年4期
蔣凱鑫，于坤霞，曹文洪，張曉明，李 鵬，孫 倩，劉 昱
黃土高原典型流域水沙變化歸因?qū)Ρ确治?/p>
蔣凱鑫1,2，于坤霞1,2※，曹文洪3，張曉明3，李鵬1,2，孫倩1,2，劉昱1,2
（1. 西安理工大學(xué)省部共建西北旱區(qū)生態(tài)水利國家重點(diǎn)實(shí)驗(yàn)室，西安 710048；2. 西安理工大學(xué)旱區(qū)生態(tài)水文與災(zāi)害防治國家林業(yè)局重點(diǎn)實(shí)驗(yàn)室，西安 710048；3. 中國水利水電科學(xué)研究院泥沙所流域水循環(huán)模擬與調(diào)控國家重點(diǎn)實(shí)驗(yàn)室，北京 100048）
近幾十年來，黃河中游黃土高原地區(qū)水沙變化顯著。眾多學(xué)者以該地區(qū)典型流域無定河開展了大量研究，但局限于方法的突破而忽視了方法間的對比分析。為此，該研究采用雙累積曲線法、水文法、彈性系數(shù)法和水保法對無定河流域1957—2010年水沙變化進(jìn)行歸因?qū)Ρ确治觥＝Y(jié)果表明：無定河流域年降水量未發(fā)生顯著變化，年潛在蒸散發(fā)、年徑流量和年輸沙量呈顯著減少趨勢且突變年份依次為1982、1971、1971年。徑流歸因?qū)Ρ确治霰砻?，水文法與彈性系數(shù)法結(jié)果較為接近，水保法人類活動貢獻(xiàn)率偏??；輸沙歸因?qū)Ρ确治霰砻鳎７ㄅc其他3種方法結(jié)果一致。歸因分析表明，人類活動是流域水沙變化的主控因素且逐年增強(qiáng)，人類活動中對徑流變化起主導(dǎo)作用的是灌溉引水和造林；對輸沙變化起主導(dǎo)作用的是淤地壩和造林。淤地壩的淤滿或失效對流域水保措施減水減沙效益影響較大。因此，未來無定河流域的治理應(yīng)繼續(xù)落實(shí)退耕還林還草政策，并加強(qiáng)淤地壩等工程措施的布局優(yōu)化和后期維護(hù)。
徑流；流域；輸沙；黃土高原；無定河；歸因?qū)Ρ确治?/p>
0 引 言
黃土高原是中國水土流失的重災(zāi)區(qū)，也是中國生態(tài)建設(shè)治理的重點(diǎn)區(qū)域[1-2]。近幾十年來，氣候變化、日益頻繁發(fā)生的極端降雨事件[3]及愈加強(qiáng)烈的人類活動[4]，使得黃土高原地區(qū)徑流、輸沙發(fā)生了顯著變化[5]。黃土高原作為黃河流域的主要組成部分[6]，貢獻(xiàn)了入黃泥沙的近90%[7]，自20世紀(jì)70年代以來治理后，減沙幅度最大的區(qū)域主要分布在該地區(qū)以無定河為中心的極強(qiáng)烈侵蝕區(qū)和汾河流域的大部分地區(qū)[8]。因此諸多學(xué)者以無定河流域?yàn)槔_展水沙變化研究。張鑫等[9]應(yīng)用重標(biāo)度極差分析法對無定河流域分析發(fā)現(xiàn)趙石窯站和白家川站的徑流量與輸沙量具有繼續(xù)下降的趨勢，人類活動是引起這種變化的主要原因。Liang等[10]基于Budyko的彈性系數(shù)法和分解法量化了無定河流域徑流對氣候變化和人類活動的響應(yīng)，結(jié)果表明徑流變異性對人類活動最敏感，人類活動貢獻(xiàn)率為65%。孫倩等[11]對黃河中游多沙粗沙區(qū)典型流域徑流與輸沙量分析表明，徑流量和輸沙量均有顯著減少趨勢，人類活動對無定河流域徑流減少、輸沙減少的貢獻(xiàn)率分別為75.5%、79.8%。Han等[12]通過收集無定河流域水文氣象數(shù)據(jù)，將其劃分為5個時期，結(jié)果表明，水保措施對降雨-徑流-輸沙關(guān)系的變化起著重要的作用，且不同時期具體發(fā)揮作用的水保措施存在差異。
既有研究的水文要素資料序列及來源站點(diǎn)不一致等造成劃分時段不同、分析結(jié)果存在差異。同時關(guān)于方法間計(jì)算結(jié)果對比及各種方法適用性的研究較少。為此本研究統(tǒng)一劃分時期，采用不同的貢獻(xiàn)率計(jì)算方法，分析對比不同影響因素對流域徑流輸沙變化的貢獻(xiàn)率，討論不同貢獻(xiàn)率計(jì)算方法的適用性，以期為黃河中游黃土高原地區(qū)急劇變化的水沙關(guān)系提供一定的理論支持。
1 材料與方法
1.1 研究區(qū)概況
無定河是黃河中游一級支流，是陜北榆林地區(qū)最大的河流[13]。無定河發(fā)源于陜北定邊、吳起、靖邊三縣交界的白于山北麓，流經(jīng)靖邊、橫山、米脂、綏德、清澗，最終匯入黃河。無定河流域干流全長491.1 km，總面積30 261 km2，河道平均比降1.97‰，把口站白家川水文站集水面積為29 662 km2，占流域面積的98.02%。由于流域內(nèi)惡劣的自然條件和人類活動，其土壤侵蝕速率是黃土高原地區(qū)最高的[14]。無定河流域位置見圖1。
圖1 無定河流域位置圖
1.2 數(shù)據(jù)來源
1957—2010年白家川水文站年徑流量、年輸沙量及流域內(nèi)雨量站的點(diǎn)年降水量資料均來自于《黃河流域水文年鑒》。無定河流域面降水量通過泰森多邊形法求得。流域年潛在蒸發(fā)數(shù)據(jù)來自于陜西省水文局。流域梯田、造林、種草、壩地及封禁治理等面積來自《榆林市水土保持統(tǒng)計(jì)年鑒》。
1.3 研究方法
1.3.1 趨勢分析
本研究采用Mann-Kendall趨勢檢驗(yàn)、線性趨勢檢驗(yàn)、滑動平均法及Pettitt突變點(diǎn)檢驗(yàn)對無定河流域1957—2010年的年降水量、年潛在蒸發(fā)量、年徑流量及年輸沙量進(jìn)行趨勢分析和突變點(diǎn)檢驗(yàn)。
1.3.2 歸因分析
1）雙累積曲線法
雙累積曲線法是將突變年份前的時段作為研究流域的基準(zhǔn)期，即認(rèn)為此時段人類活動可忽略；突變年份之后作為研究期，即認(rèn)為此時段流域的徑流、輸沙受到了氣候變化和人類活動的共同作用。將研究期降雨數(shù)據(jù)代入基準(zhǔn)期建立的降雨-徑流、降雨-輸沙關(guān)系，進(jìn)而獲得研究期無人類活動下的徑流、輸沙模擬值，其與研究期徑流、輸沙實(shí)測值的差值作為人類活動對水沙變化的影響。
2）水文法
水文法是通過對流域基準(zhǔn)期的降雨-徑流、降雨-輸沙關(guān)系進(jìn)行回歸分析優(yōu)選，參數(shù)率定，模型確定，然后計(jì)算出流域大規(guī)模開展水土流失治理時期的徑流輸沙模擬值，扣除實(shí)測值，獲得流域水利水保措施的減水減沙效益。無定河流域常見水文法模型[15-18]見表 1所示。
為了優(yōu)選出相對最佳的模型，本研究選取納什效率系數(shù)（nash-sutcliffe efficiency，NSE）和相對誤差（relative error，RE）作為評價(jià)指標(biāo)進(jìn)行模型評價(jià)。
表1 無定河常見水文法模型
注：指徑流模擬值；指輸沙模擬值；P指年降水量；P指汛期（6—9月）降雨量；7+8指7—8月降雨量；6+9指6月和9月降雨量；、、為模型參數(shù)；下同。
Note:refers to simulated value of runoff;refers to simulated value of sediment transport; Prefers to the annual precipitation; Prefers to precipitation during flood season;7+8refers to precipitation from July to August;6+9refers to precipitation in June and September;,andrefer to the model parameters; the same below.
3）彈性系數(shù)法
彈性系數(shù)定義為因變量變化率與自變量變化率的比值，它被認(rèn)為是徑流對其他變量（如降水、潛在蒸散發(fā)）敏感性的有效指標(biāo)[19]。假設(shè)研究流域地形和土壤在研究時段內(nèi)不發(fā)生變化，則在長時段內(nèi)，徑流（或輸沙）的變化量可表示為：
式中D表示徑流（或輸沙）的總變化量；DQ和DQ分別表示氣候變化和人類活動導(dǎo)致的徑流（或輸沙）變化量。
氣候變化對徑流（或輸沙）變化量的影響量為：
式中和分別為年降水量和年潛在蒸散發(fā)量對年徑流量的彈性系數(shù)；D和D0分別為降水和潛在蒸散發(fā)的變化。
① 計(jì)算徑流貢獻(xiàn)率的彈性系數(shù)
本研究徑流貢獻(xiàn)率的降水彈性系數(shù)采用Sankarasubramanian等[20]通過實(shí)測水文氣象序列得到的。()函數(shù)采用Zhang等[21]的公式。
式中表示干燥系數(shù)，即年潛在蒸散發(fā)與年降水的比值；表示一個地區(qū)植被覆被條件的參數(shù)，本研究通過繪制0/&/相關(guān)圖，與Zhang等[21]繪制的0/&/圖進(jìn)行比較，求得=2.0。
② 計(jì)算輸沙貢獻(xiàn)率的彈性系數(shù)
為了消除流域面積對徑流和輸沙關(guān)系的影響，本研究對流域徑流、輸沙均除以流域面積，即擬合徑流深和輸沙模數(shù)的關(guān)系：
輸沙公式推導(dǎo)具體如下：
4）水保法
本研究采用的水保法為指標(biāo)法，即核實(shí)流域內(nèi)各水利水保措施面積等，確定不同地貌類型區(qū)不同條件、時期下各措施減水減沙指標(biāo)，分項(xiàng)計(jì)算，逐項(xiàng)累加，得到流域水利水保措施的減水減沙量。參考水利部第一期黃河水沙變化基金項(xiàng)目研究成果[18,22-24]，確定流域各項(xiàng)水利水保措施減水減沙指標(biāo)，見表2。采用無定河流域2011年水利普查的骨干壩數(shù)據(jù)、2009年淤地壩安全大檢查數(shù)據(jù)等[25]分析計(jì)算無定河流域截至2010年底流域內(nèi)淤地壩實(shí)際攔沙總量，按不同年份淤地壩控制面積、產(chǎn)沙降雨指標(biāo)權(quán)系數(shù)分配的方法，反推不同時期淤地壩攔沙量。流域灌溉定額參考《陜西省水資源利用》取5 940 m3/hm2，灌溉水回歸系數(shù)取15%，灌溉引水平均含沙量據(jù)調(diào)查取10 kg/m3。
表2 無定河流域水保措施減水減沙效益
2 結(jié)果與分析
2.1 趨勢與變點(diǎn)分析
對無定河流域年降水量、年潛在蒸散發(fā)、年徑流量和年輸沙量進(jìn)行趨勢分析，見圖2。結(jié)果表明無定河流域的年降水量（圖2a）變化趨緩，其5a滑動平均值在系列均值附近波動；而流域年潛在蒸散發(fā)（圖2b）、年徑流量（圖2c）和年輸沙量（圖2d）均呈明顯減少趨勢。
本研究進(jìn)一步采用M-K檢驗(yàn)、Pettitt檢驗(yàn)分別進(jìn)行系列一致性分析和突變分析，結(jié)果見表3。結(jié)果表明：年降水量呈不顯著增加趨勢，年潛在蒸散發(fā)、年徑流量及年輸沙量均呈極顯著減少趨勢；年潛在蒸散發(fā)突變年份為1982年，年徑流量和年輸沙量的突變年份均為1971年。上述分析結(jié)果與牟夏等[26]、任宗萍等[27]分析結(jié)果一致。
2.2 水文法模型優(yōu)選
不同水文法模型擬合結(jié)果見表4。徑流方面，除模型1以外的3個水文法模型納什效率系數(shù)NSE和相對誤差RE基本一致；輸沙方面，模型1、2、3的NSE和RE基本一致，而模型4表現(xiàn)相對較好。綜合來看，模型4的NSE更接近于1且RE更小，即模擬效果相對更好，故本研究采用模型4代表水文法對無定河流域水沙變化進(jìn)行歸因分析。
圖2 無定河流域水文要素年際變化曲線
表3 無定河流域水文要素M-K檢驗(yàn)與Pettitt檢驗(yàn)計(jì)算結(jié)果
注：、分別表示M-K檢驗(yàn)、Pettitt檢驗(yàn)統(tǒng)計(jì)量；***為99.9%置信水平。
Note:andvalue represent Mann-Kendall test and Pettitt test statistics respectively. *** is 99.9% confidence level.
表4 無定河水文法模型參數(shù)率定及優(yōu)選
2.3 歸因?qū)Ρ确治?/h3>
在無定河流域內(nèi)，分別采用彈性系數(shù)法、雙累積曲線法、水文法和水保法計(jì)算氣候變化和人類活動對流域徑流輸沙變化的貢獻(xiàn)率，其中基準(zhǔn)期為1957—1971年。結(jié)果見表5。總體來看，在無定河流域，人類活動為徑流輸沙變化的主導(dǎo)因素，其對輸沙減少的作用較對徑流減少的作用大。人類活動中灌溉引水和造林對流域徑流變化起主導(dǎo)作用；淤地壩和造林對流域輸沙變化起主導(dǎo)作用。淤地壩對流域徑流減少的效益較小，但對流域輸沙量的減少占據(jù)了主導(dǎo)作用，貢獻(xiàn)率達(dá)43.3%。
徑流方面，雙累積曲線法、水文法、彈性系數(shù)法的計(jì)算結(jié)果較為一致，水保法對人類活動的計(jì)算結(jié)果相對偏?。惠斏撤矫?，4種方法的結(jié)果均表現(xiàn)出人類活動占據(jù)主導(dǎo)作用，具體為雙累積曲線法與彈性系數(shù)法的計(jì)算結(jié)果較為一致，水文法和水保法的計(jì)算結(jié)果較為一致。
表5 無定河流域年徑流輸沙變化歸因分析
進(jìn)一步對無定河流域內(nèi)各時期氣候變化和人類活動對水沙變化的影響進(jìn)行分析，見圖3。人類活動對流域徑流、輸沙變化的貢獻(xiàn)率（圖3b、圖3d）明顯較氣候變化對徑流、輸沙變化的貢獻(xiàn)率（圖3a、圖3c）大，且總體來看人類活動對徑流、輸沙變化的貢獻(xiàn)率呈增加的趨勢。人類活動對徑流、輸沙的貢獻(xiàn)率在1990年左右出現(xiàn)一低谷，同時該低谷對徑流的影響（圖3b）較對輸沙的影響（圖3d）大，即與對輸沙的影響相比，人類活動對徑流的影響進(jìn)一步減弱。
圖3 無定河流域水沙變化不同方法對比分析
3 討 論
本研究雙累積曲線法僅考慮年降水量單一自變量，并將其與徑流輸沙間的關(guān)系概化為線性，這在水保措施少、治理標(biāo)準(zhǔn)低、大雨產(chǎn)大沙特征明顯的基準(zhǔn)期有較好的適用性[28]，但并未考慮黃土高原場次降雨中產(chǎn)流產(chǎn)沙的侵蝕性降雨[29]、降雨年內(nèi)分布與基準(zhǔn)相比發(fā)生的變化[30]及坡面、溝道治理措施的持續(xù)開展和減水減沙作用的滯后顯現(xiàn)，即存在年降水量之外的因素對徑流、輸沙產(chǎn)生著不可忽視的作用，或造成計(jì)算結(jié)果中人類活動對水沙變化的貢獻(xiàn)率計(jì)算偏小。為提高突變后研究期模擬值的準(zhǔn)確發(fā)生，一種是在降雨-徑流、降雨-輸沙關(guān)系中引入一個反映年內(nèi)降雨分配的參數(shù)，使關(guān)系式的物理意義更加清晰，如本研究中采用的水文法；另一種是從流域內(nèi)長時段水量平衡原理出發(fā)，考慮各水沙變化潛在影響因素并計(jì)算相關(guān)貢獻(xiàn)率，如本研究中采用的彈性系數(shù)法。在降雨-徑流歸因分析結(jié)果中，見表5，可發(fā)現(xiàn)水文法和彈性系數(shù)法的結(jié)果基本一致，與雙累積曲線法相比結(jié)果更為合理。水保法尚無法考慮影響流域徑流輸沙變化的所有因素，如不同措施之間的耦合作用、淤地壩淤積比及其異地減沙作用、流域內(nèi)極端事件（干旱或暴雨）對水利水保措施的影響、流域內(nèi)社會經(jīng)濟(jì)發(fā)展帶來的道路建設(shè)及取用河沙等。
人類活動對徑流、輸沙貢獻(xiàn)率在1990年左右出現(xiàn)低谷的原因在于流域內(nèi)20世紀(jì)70、80年代修建完成的大部分中小型淤地壩因攔蓄壽命多小于20 a，至1990年左右時已淤滿或已失效[31]，而退耕還林還草政策、淤地壩作為之一的水利部“亮點(diǎn)”工程等尚未出臺實(shí)施。而低谷對徑流影響較輸沙大的原因在于淤地壩淤滿之后雖失去了攔蓄泥沙的庫容，但因攔蓄泥沙而抬高的壩地有效增大了溝道過水?dāng)嗝?，使攜帶的泥沙沿程沉降，而徑流則可以越過淤滿的壩頂進(jìn)入河道。流域內(nèi)輸沙減少幅度較徑流的大，水利水保措施對下墊面的改變有效減緩或阻止了產(chǎn)匯流的形成，而下滲的徑流又促進(jìn)了林草措施進(jìn)一步發(fā)揮作用。
4 結(jié) 論
本研究選取黃河中游黃土高原地區(qū)典型流域無定河流域，以其年降水量、年潛在蒸散發(fā)、年輸沙量和年徑流量及水利水保措施為研究對象，通過趨勢突變分析和歸因?qū)Ρ确治?，得出不同方法下氣候變化和人類活動對徑流輸沙的貢獻(xiàn)率，結(jié)果表明：
1）通過M-K檢驗(yàn)和Pettitt檢驗(yàn)得出，無定河流域1957—2010年年降水量呈不顯著增加趨勢；年徑流量和年輸沙量呈顯著減少趨勢，其突變年份均為1971年。
2）不同方法分析結(jié)果表明，徑流方面，水文法和彈性系數(shù)法物理意義較雙累積曲線法更為清晰，三者計(jì)算結(jié)果較水保法更為合理；輸沙方面，4種方法分析結(jié)果均為人類活動占據(jù)主導(dǎo)作用，區(qū)別在于對人類活動的水沙變化貢獻(xiàn)率具體大小方面，即水保法和水文法結(jié)果較為一致，彈性系數(shù)法和雙累積曲線法結(jié)果較為一致。
3）不同徑流輸沙影響因素貢獻(xiàn)率計(jì)算方法的結(jié)果表明，人類活動占據(jù)主導(dǎo)作用且其作用逐漸增強(qiáng)。水利水保措施對輸沙變化的作用較對徑流變化的強(qiáng)。人類活動中，徑流方面，灌溉引水是一項(xiàng)重要因素；輸沙方面，淤地壩是短期內(nèi)減少流域輸沙的有效措施，而20世紀(jì)集中修建的淤地壩的淤滿或失效在短期內(nèi)將會削弱水利水保措施的減水減沙效益。長遠(yuǎn)來看，流域內(nèi)持續(xù)開展的退耕還林還草等生態(tài)建設(shè)對徑流、輸沙的影響將逐漸發(fā)揮主導(dǎo)作用。
[1]劉國彬，上官周平，姚文藝，等. 黃土高原生態(tài)工程的生態(tài)成效[J]. 中國科學(xué)院院刊，2017，32(1)：11－19. Liu Guobin, Shangguan Zhouping, Yao Wenyi, et al. Ecological effects of soil conservation in Loess Plateau[J]. Bulletin of Chinese Academy of Sciences, 2017, 32(1): 11－19. (in Chinese with English abstract)
[2]Wang F, Mu X, Li R, et al. Co-evolution of soil and water conservation policy and human–environment linkages in the Yellow River Basin since 1949[J]. Science of The Total Environment, 2015, 508: 166－177.
[3]Zhao Y, Cao W, Hu C, et al. Analysis of changes in characteristics of flood and sediment yield in typical basins of the Yellow River under extreme rainfall events[J]. Catena, 2019, 177: 31－40.
[4]Kong D, Miao C, Wu J, et al. Impact assessment of climate change and human activities on net runoff in the Yellow River Basin from 1951 to 2012[J]. Ecological Engineering, 2016, 91: 566－573.
[5]Zheng H, Miao C, Wu J, et al. Temporal and spatial variations in water discharge and sediment load on the Loess Plateau, China: A high-density study[J]. Science of the Total Environment, 2019, 666: 875－886.
[6]高海東，李占斌，李鵬，等. 基于土壤侵蝕控制度的黃土高原水土流失治理潛力研究[J]. 地理學(xué)報(bào)，2015，70(9)：1503－1515. Gao Haidong, Li Zhanbin, Li Peng, et al. The capacity of soil loss control in the Loess Plateau based on soil erosion control degree[J]. Acta Geographica Sinica, 2015, 70(9): 1503－1515. (in Chinese with English abstract)
[7]Wang S, Fu B, Piao S, et al. Reduced sediment transport in the Yellow River due to anthropogenic changes[J]. Nature Geoscience, 2015, 9: 38－41.
[8]王萬忠，焦菊英. 黃土高原侵蝕產(chǎn)沙強(qiáng)度的時空變化特征[J]. 地理學(xué)報(bào)，2002，57(2)：210－217. Wang Wanzhong, Jiao Juying. Temporal and spatial variation characteristics of erosion and sediment yield in the Loess Plateau[J]. Acta Geographica Sinica, 2002, 57(2): 210－217. (in Chinese with English abstract)
[9]張鑫，蔡煥杰，尹曉楠. 應(yīng)用重標(biāo)度極差分析法（R/S）分析無定河流域水沙變化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2010，26(增刊2)：212－217. Zhang Xin, Cai Huanjie, Yin Xiaonan. Variation trends analysis of runoff and sediment time series based on R/S method in Wuding River Basin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(Supp.2): 212－217. (in Chinese with English abstract)
[10]Liang W, Bai D, Jin Z, et al. A study on the streamflow change and its relationship with climate change and ecological restoration measures in a sediment concentrated region in the Loess Plateau, China[J]. Water Resources Management, 2015, 29(11): 4045－4060.
[11]孫倩，于坤霞，李占斌，等. 黃河中游多沙粗沙區(qū)水沙變化趨勢及其主控因素的貢獻(xiàn)率[J]. 地理學(xué)報(bào)，2018，73(5)：945－956. Sun Qian, Yu Kunxia, Li Zhanbin, et al. The trends of streamflow and sediment and their driving factors in the middle reaches of the Yellow River[J]. Acta Geographica Sinica, 2018, 73(5): 945－956. (in Chinese with English abstract)
[12]Han J, Gao J, Luo H. Changes and implications of the relationship between rainfall, runoff and sediment load in the Wuding River basin on the Chinese Loess Plateau[J]. Catena, 2019, 175: 228－235.
[13]陳思. VIC大尺度陸面水文模型在無定河流域的應(yīng)用研究[D]. 北京：中國地質(zhì)大學(xué)(北京)，2018. Chen Si. The Application of Variable Infiltration Capacity (VIC) Macroscale Hydrologic Model in Wuding River Basin[D]. Beijing: China University of Geosciences(Beijing), 2018. (in Chinese with English abstract)
[14]Shi C, Zhou Y, Fan X, et al. A study on the annual runoff change and its relationship with water and soil conservation practices and climate change in the middle Yellow River basin[J]. Catena, 2013, 100: 31－41.
[15]賈紹鳳，梁季陽. 黃土高原降雨、徑流、產(chǎn)沙相互關(guān)系的研究[J]. 水土保持學(xué)報(bào)，1992(3)：42－47. Jia Shaofeng, Liang Jiyang. The correlation among rainfall, runoff and sediment generation in the Loess Plateau[J]. Journal of Soil and Water Conservation, 1992(3): 42－47. (in Chinese with English abstract)
[16]Li L, 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 Processes, 2007, 21: 3485－3491.
[17]湯立群，陳國祥. 水土保持減水減沙效益計(jì)算方法研究[J]. 河海大學(xué)學(xué)報(bào)：自然科學(xué)版，1999(1)：82－87. Tang Liqun, Chen Guoxiang. Study on calculation method of water and soil conservation water and sediment reduction benefit[J]. Journal of Hohai University(Natural Sciences), 1999(1): 82－87. (in Chinese with English abstract)
[18]張經(jīng)濟(jì)，冀文慧，馮曉東. 無定河流域水沙變化現(xiàn)狀、成因及發(fā)展趨勢的研究[M]//汪崗，范昭. 黃河水沙變化研究（第二卷）. 鄭州：黃河水利出版社，2002. 393－429.
[19]Li Y Z, Liu C M, Zhang D, et al. Reduced runoff due to anthropogenic intervention in the Loess Plateau, China[J]. Water, 2016, 8(458): 1－16.
[20]Sankarasubramanian A, Vogel R M. Hydroclimatology of the continental United States[J]. Geophysical Research Letters, 2003, 30(7):314－326.
[21]Zhang L, Dawes W R, Walker G R. Response of mean annual evapotranspiration to vegetation changes at catchment scale[J]. Water Resources Research, 2001, 37(3): 701－708.
[22]惠養(yǎng)瑜，冀文慧，同信奇，等. 無定河流域水沙變化及其發(fā)展趨勢預(yù)測研究[M]//汪崗，范昭. 黃河水沙變化研究（第一卷：下冊）. 鄭州：黃河水利出版社，2002. 689－714.
[23]左仲國. 清水河流域水沙變化及發(fā)展趨勢預(yù)測[M]//汪崗，范昭. 黃河水沙變化研究（第一卷：上冊）. 鄭州：黃河水利出版社，2002. 439－448.
[24]冉大川，張棟，焦鵬，等. 西柳溝流域近期水沙變化歸因分析[J]. 干旱區(qū)資源與環(huán)境，2016，30(5)：143－149. Ran Dachuan, Zhang Dong, Jiao Peng, et al. Analysis on the attribution of water and sediment changes in Xiliugou basin[J]. Journal of Arid Land Resources and Environment, 2016, 30(5): 143－149. (in Chinese with English abstract)
[25]韓向楠，謝世友，高云飛. 近年無定河流域淤地壩攔沙作用研究[J]. 人民黃河，2018，40(11)：9－12，41. Han Xiangnan, Xie Shiyou, Gao Yunfei. Study on sediment retaining effect of check dams in Wuding River Basin in recent years[J]. Yellow River, 2018, 40(11): 9－12, 41. (in Chinese with English abstract)
[26]牟夏，劉強(qiáng)，平凡，等. 1960—2007年無定河流域徑流情勢變化及其歸因分析[J]. 北京師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2017，53(4)：465－471. Mou Xia, Liu Qiang, Ping Fan, et al. Analysis of changes in annual streamflow and its influence factors in Wuding River basin in the last 50 years (1960—2007)[J]. Journal of Beijing Normal University(Natural Sciences), 2017, 53(4): 465－471. (in Chinese with English abstract)
[27]任宗萍，馬勇勇，王友勝，等. 生態(tài)建設(shè)條件下無定河不同地貌區(qū)徑流變化歸因分析[J]. 生態(tài)學(xué)報(bào)，2019，39(7)：1－9. Ren Zongping, Ma Yongyong, Wang Yousheng, et al. Runoff changes and attribution analysis in tributaries of different geomorphic regions in Wuding River under ecological construction[J]. Acta Ecologica Sinica, 2019, 39(7): 1－9. (in Chinese with English abstract)
[28]李勉，李平，楊二，等. 黃土丘陵區(qū)淤地壩建設(shè)后小流域泥沙攔蓄與輸移特征[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(18)：80－86. Li Mian, Li Ping, Yang Er, et al. Characteristics of sediment retention and transport in small watershed after construction of check dams in Loess Hilly Area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(18): 80－86. (in Chinese with English abstract)
[29]張攀，姚文藝，唐洪武，等. 黃土坡面細(xì)溝形態(tài)變化及對侵蝕產(chǎn)沙過程的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34(5)：114－119. Zhang Pan, Yao Wenyi, Tang Hongwu, et al. Rill morphology change and its effect on erosion and sediment yield on loess slope[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(5): 114－119. (in Chinese with English abstract)
[30]Tang X, Miao C, Xi Y, et al. Analysis of precipitation characteristics on the loess plateau between 1965 and 2014, based on high-density gauge observations[J]. Atmospheric Research, 2018, 213: 264－274.
[31]高云飛，郭玉濤，劉曉燕，等. 陜北黃河中游淤地壩攔沙功能失效的判斷標(biāo)準(zhǔn)[J]. 地理學(xué)報(bào)，2014，69(1)：73－79. Gao Yunfei, Guo Yutao, Liu Xiaoyan, et al. Failure criteria of the warping dams on sediment interception in the Middle Yellow River in northern Shaanxi[J]. Acta Geographica Sinica, 2014, 69(1): 73－79. (in Chinese with English abstract)
Attributional comparative analysis of runoff and sediment change in typical basin of Loess Plateau
Jiang Kaixin1,2, Yu Kunxia1,2※, Cao Wenhong3, Zhang Xiaoming3, Li Peng1,2, Sun Qian1,2, Liu Yu1,2
(1.-,,710048,; 2.,,710048,; 3.,,100048,)
In recent decades, the runoff and sediment of the Yellow River middle reaches located in the Loess Plateau have changed significantly. Many scholars have carried out a lot of researches on the Wuding River, a typical basin of the Loess Plateau, but these researches are focused on the breakthrough of methods and ignore the comparative analysis among methods. This study compared the difference of four methods, i.e., double cumulative curve method, hydrological method, elastic coefficient method and soil and water conservation method, to conduct attributional comparative analysis of runoff and sediment change in the Wuding River Basin from 1957 to 2010. Both double cumulative curve method and hydrological method separate the influence of climate change and human activity by substituting the precipitation data in the research period into the precipitation-runoff and precipitation-sediment transport relationships established in the reference period, and the hydrological method takes into account the annual distribution of precipitation. The elastic coefficient method generates runoff and sediment elasticities to quantify the effects of climate change and human activities. The soil and water conservation method refers to determining the runoff and sediment reduction indicators of each measure under different conditions and periods in different geomorphic areas. The results of trend and change point analysis showed that the annual precipitation in the Wuding River Basin has not changed significantly. The annual potential evaporation, annual runoff and annual sediment transport decreased significantly and the year of change point was 1982, 1971 and 1971, respectively. The attributional comparative analysis of runoff showed that the results of hydrological method were close to those of elastic coefficient method, and the contribution rate of human activities of soil and water conservation method was relatively small, and the attributional comparative analysis of sediment transport showed that the results of soil and water conservation method were consistent with those of the other three methods. Attribution analysis results indicated human activities were the dominant factor of runoff and sediment change in the basin, and their effects on runoff and sediment change were increasing year by year. Among different human activities, irrigation water diversion and afforestation played a leading role in runoff change, and check dam and afforestation played a leading role in sediment transport change. The full up or failure of check dam built in the 1970s and 1980s would largely weaken the runoff and sediment reduction benefits of water conservancy and soil conservation measures in the short term. In the long term, the impact of ecological construction, such as returning cultivated land to forestland or grassland, on runoff and sediment transport will gradually play a major role. Therefore, in the future, the governance of the Wuding River Basin should continue to implement the policy of conversion of cropland to forest and grassland, and strengthen the layout optimization and post-maintenance of engineering measures such as check dam.
runoff; watershed; sediment transport; Loess Plateau; Wuding River basin; attributional comparative analysis
蔣凱鑫，于坤霞，曹文洪，張曉明，李鵬，孫倩，劉昱. 黃土高原典型流域水沙變化歸因?qū)Ρ确治鯷J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(4)：143－149. doi：10.11975/j.issn.1002-6819.2020.04.017 http://www.tcsae.org
Jiang Kaixin, Yu Kunxia, Cao Wenhong, Zhang Xiaoming, Li Peng, Sun Qian, Liu Yu. Attributional comparative analysis of runoff and sediment change in typical basin of Loess Plateau[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(4): 143－149. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.04.017 http://www.tcsae.org
2019-07-08
2019-10-16
國家自然科學(xué)基金面上項(xiàng)目（51879281，41877077）；國家重點(diǎn)研發(fā)計(jì)劃（2016YFC0402407）
蔣凱鑫，博士生，主要從事流域侵蝕動力學(xué)研究。Email：jiangkaizn@qq.com
于坤霞，博士，副教授，主要從事水文水資源研究。Email：yukunxia@126.com
10.11975/j.issn.1002-6819.2020.04.017
S157.1
A
1002-6819(2020)-04-0143-07
猜你喜歡
格陵蘭島積雪區(qū)地表徑流增加研究水利水電快報(bào)(2022年8期)2022-11-23 黃丘區(qū)小流域暴雨徑流動力特征對淤地壩配置的響應(yīng)農(nóng)業(yè)工程學(xué)報(bào)(2022年10期)2022-08-22 基于SWAT模型的布爾哈通河流域徑流模擬研究黑龍江大學(xué)自然科學(xué)學(xué)報(bào)(2022年1期)2022-03-29 關(guān)于新時期淤地壩建設(shè)管理的思考現(xiàn)代農(nóng)村科技(2021年9期)2021-12-06 山西省淤地壩建設(shè)管理專項(xiàng)提升三年（2021-2023年）行動方案山西水土保持科技(2021年4期)2021-12-05 黃土高原地區(qū)淤地壩存在問題分析農(nóng)業(yè)科技與信息(2020年15期)2020-12-18 高中生學(xué)習(xí)成績與成就歸因方式的相關(guān)研究魅力中國(2020年39期)2020-12-08 雅魯藏布江河川徑流變化的季節(jié)性規(guī)律探索品牌研究(2020年32期)2020-08-09 近40年來蒲河流域徑流變化及影響因素分析人民珠江(2019年4期)2019-04-20 為何越來越愛愛你(2018年22期)2018-11-14
蔣凱鑫，于坤霞，曹文洪，張曉明，李 鵬，孫 倩，劉 昱
黃土高原典型流域水沙變化歸因?qū)Ρ确治?/p>
蔣凱鑫1,2，于坤霞1,2※，曹文洪3，張曉明3，李鵬1,2，孫倩1,2，劉昱1,2
（1. 西安理工大學(xué)省部共建西北旱區(qū)生態(tài)水利國家重點(diǎn)實(shí)驗(yàn)室，西安 710048；2. 西安理工大學(xué)旱區(qū)生態(tài)水文與災(zāi)害防治國家林業(yè)局重點(diǎn)實(shí)驗(yàn)室，西安 710048；3. 中國水利水電科學(xué)研究院泥沙所流域水循環(huán)模擬與調(diào)控國家重點(diǎn)實(shí)驗(yàn)室，北京 100048）
近幾十年來，黃河中游黃土高原地區(qū)水沙變化顯著。眾多學(xué)者以該地區(qū)典型流域無定河開展了大量研究，但局限于方法的突破而忽視了方法間的對比分析。為此，該研究采用雙累積曲線法、水文法、彈性系數(shù)法和水保法對無定河流域1957—2010年水沙變化進(jìn)行歸因?qū)Ρ确治觥＝Y(jié)果表明：無定河流域年降水量未發(fā)生顯著變化，年潛在蒸散發(fā)、年徑流量和年輸沙量呈顯著減少趨勢且突變年份依次為1982、1971、1971年。徑流歸因?qū)Ρ确治霰砻?，水文法與彈性系數(shù)法結(jié)果較為接近，水保法人類活動貢獻(xiàn)率偏??；輸沙歸因?qū)Ρ确治霰砻鳎７ㄅc其他3種方法結(jié)果一致。歸因分析表明，人類活動是流域水沙變化的主控因素且逐年增強(qiáng)，人類活動中對徑流變化起主導(dǎo)作用的是灌溉引水和造林；對輸沙變化起主導(dǎo)作用的是淤地壩和造林。淤地壩的淤滿或失效對流域水保措施減水減沙效益影響較大。因此，未來無定河流域的治理應(yīng)繼續(xù)落實(shí)退耕還林還草政策，并加強(qiáng)淤地壩等工程措施的布局優(yōu)化和后期維護(hù)。
徑流；流域；輸沙；黃土高原；無定河；歸因?qū)Ρ确治?/p>
0 引 言
黃土高原是中國水土流失的重災(zāi)區(qū)，也是中國生態(tài)建設(shè)治理的重點(diǎn)區(qū)域[1-2]。近幾十年來，氣候變化、日益頻繁發(fā)生的極端降雨事件[3]及愈加強(qiáng)烈的人類活動[4]，使得黃土高原地區(qū)徑流、輸沙發(fā)生了顯著變化[5]。黃土高原作為黃河流域的主要組成部分[6]，貢獻(xiàn)了入黃泥沙的近90%[7]，自20世紀(jì)70年代以來治理后，減沙幅度最大的區(qū)域主要分布在該地區(qū)以無定河為中心的極強(qiáng)烈侵蝕區(qū)和汾河流域的大部分地區(qū)[8]。因此諸多學(xué)者以無定河流域?yàn)槔_展水沙變化研究。張鑫等[9]應(yīng)用重標(biāo)度極差分析法對無定河流域分析發(fā)現(xiàn)趙石窯站和白家川站的徑流量與輸沙量具有繼續(xù)下降的趨勢，人類活動是引起這種變化的主要原因。Liang等[10]基于Budyko的彈性系數(shù)法和分解法量化了無定河流域徑流對氣候變化和人類活動的響應(yīng)，結(jié)果表明徑流變異性對人類活動最敏感，人類活動貢獻(xiàn)率為65%。孫倩等[11]對黃河中游多沙粗沙區(qū)典型流域徑流與輸沙量分析表明，徑流量和輸沙量均有顯著減少趨勢，人類活動對無定河流域徑流減少、輸沙減少的貢獻(xiàn)率分別為75.5%、79.8%。Han等[12]通過收集無定河流域水文氣象數(shù)據(jù)，將其劃分為5個時期，結(jié)果表明，水保措施對降雨-徑流-輸沙關(guān)系的變化起著重要的作用，且不同時期具體發(fā)揮作用的水保措施存在差異。
既有研究的水文要素資料序列及來源站點(diǎn)不一致等造成劃分時段不同、分析結(jié)果存在差異。同時關(guān)于方法間計(jì)算結(jié)果對比及各種方法適用性的研究較少。為此本研究統(tǒng)一劃分時期，采用不同的貢獻(xiàn)率計(jì)算方法，分析對比不同影響因素對流域徑流輸沙變化的貢獻(xiàn)率，討論不同貢獻(xiàn)率計(jì)算方法的適用性，以期為黃河中游黃土高原地區(qū)急劇變化的水沙關(guān)系提供一定的理論支持。
1 材料與方法
1.1 研究區(qū)概況
無定河是黃河中游一級支流，是陜北榆林地區(qū)最大的河流[13]。無定河發(fā)源于陜北定邊、吳起、靖邊三縣交界的白于山北麓，流經(jīng)靖邊、橫山、米脂、綏德、清澗，最終匯入黃河。無定河流域干流全長491.1 km，總面積30 261 km2，河道平均比降1.97‰，把口站白家川水文站集水面積為29 662 km2，占流域面積的98.02%。由于流域內(nèi)惡劣的自然條件和人類活動，其土壤侵蝕速率是黃土高原地區(qū)最高的[14]。無定河流域位置見圖1。
圖1 無定河流域位置圖
1.2 數(shù)據(jù)來源
1957—2010年白家川水文站年徑流量、年輸沙量及流域內(nèi)雨量站的點(diǎn)年降水量資料均來自于《黃河流域水文年鑒》。無定河流域面降水量通過泰森多邊形法求得。流域年潛在蒸發(fā)數(shù)據(jù)來自于陜西省水文局。流域梯田、造林、種草、壩地及封禁治理等面積來自《榆林市水土保持統(tǒng)計(jì)年鑒》。
1.3 研究方法
1.3.1 趨勢分析
本研究采用Mann-Kendall趨勢檢驗(yàn)、線性趨勢檢驗(yàn)、滑動平均法及Pettitt突變點(diǎn)檢驗(yàn)對無定河流域1957—2010年的年降水量、年潛在蒸發(fā)量、年徑流量及年輸沙量進(jìn)行趨勢分析和突變點(diǎn)檢驗(yàn)。
1.3.2 歸因分析
1）雙累積曲線法
雙累積曲線法是將突變年份前的時段作為研究流域的基準(zhǔn)期，即認(rèn)為此時段人類活動可忽略；突變年份之后作為研究期，即認(rèn)為此時段流域的徑流、輸沙受到了氣候變化和人類活動的共同作用。將研究期降雨數(shù)據(jù)代入基準(zhǔn)期建立的降雨-徑流、降雨-輸沙關(guān)系，進(jìn)而獲得研究期無人類活動下的徑流、輸沙模擬值，其與研究期徑流、輸沙實(shí)測值的差值作為人類活動對水沙變化的影響。
2）水文法
水文法是通過對流域基準(zhǔn)期的降雨-徑流、降雨-輸沙關(guān)系進(jìn)行回歸分析優(yōu)選，參數(shù)率定，模型確定，然后計(jì)算出流域大規(guī)模開展水土流失治理時期的徑流輸沙模擬值，扣除實(shí)測值，獲得流域水利水保措施的減水減沙效益。無定河流域常見水文法模型[15-18]見表 1所示。
為了優(yōu)選出相對最佳的模型，本研究選取納什效率系數(shù)（nash-sutcliffe efficiency，NSE）和相對誤差（relative error，RE）作為評價(jià)指標(biāo)進(jìn)行模型評價(jià)。
表1 無定河常見水文法模型
注：指徑流模擬值；指輸沙模擬值；P指年降水量；P指汛期（6—9月）降雨量；7+8指7—8月降雨量；6+9指6月和9月降雨量；、、為模型參數(shù)；下同。
Note:refers to simulated value of runoff;refers to simulated value of sediment transport; Prefers to the annual precipitation; Prefers to precipitation during flood season;7+8refers to precipitation from July to August;6+9refers to precipitation in June and September;,andrefer to the model parameters; the same below.
3）彈性系數(shù)法
彈性系數(shù)定義為因變量變化率與自變量變化率的比值，它被認(rèn)為是徑流對其他變量（如降水、潛在蒸散發(fā)）敏感性的有效指標(biāo)[19]。假設(shè)研究流域地形和土壤在研究時段內(nèi)不發(fā)生變化，則在長時段內(nèi)，徑流（或輸沙）的變化量可表示為：
式中D表示徑流（或輸沙）的總變化量；DQ和DQ分別表示氣候變化和人類活動導(dǎo)致的徑流（或輸沙）變化量。
氣候變化對徑流（或輸沙）變化量的影響量為：
式中和分別為年降水量和年潛在蒸散發(fā)量對年徑流量的彈性系數(shù)；D和D0分別為降水和潛在蒸散發(fā)的變化。
① 計(jì)算徑流貢獻(xiàn)率的彈性系數(shù)
本研究徑流貢獻(xiàn)率的降水彈性系數(shù)采用Sankarasubramanian等[20]通過實(shí)測水文氣象序列得到的。()函數(shù)采用Zhang等[21]的公式。
式中表示干燥系數(shù)，即年潛在蒸散發(fā)與年降水的比值；表示一個地區(qū)植被覆被條件的參數(shù)，本研究通過繪制0/&/相關(guān)圖，與Zhang等[21]繪制的0/&/圖進(jìn)行比較，求得=2.0。
② 計(jì)算輸沙貢獻(xiàn)率的彈性系數(shù)
為了消除流域面積對徑流和輸沙關(guān)系的影響，本研究對流域徑流、輸沙均除以流域面積，即擬合徑流深和輸沙模數(shù)的關(guān)系：
輸沙公式推導(dǎo)具體如下：
4）水保法
本研究采用的水保法為指標(biāo)法，即核實(shí)流域內(nèi)各水利水保措施面積等，確定不同地貌類型區(qū)不同條件、時期下各措施減水減沙指標(biāo)，分項(xiàng)計(jì)算，逐項(xiàng)累加，得到流域水利水保措施的減水減沙量。參考水利部第一期黃河水沙變化基金項(xiàng)目研究成果[18,22-24]，確定流域各項(xiàng)水利水保措施減水減沙指標(biāo)，見表2。采用無定河流域2011年水利普查的骨干壩數(shù)據(jù)、2009年淤地壩安全大檢查數(shù)據(jù)等[25]分析計(jì)算無定河流域截至2010年底流域內(nèi)淤地壩實(shí)際攔沙總量，按不同年份淤地壩控制面積、產(chǎn)沙降雨指標(biāo)權(quán)系數(shù)分配的方法，反推不同時期淤地壩攔沙量。流域灌溉定額參考《陜西省水資源利用》取5 940 m3/hm2，灌溉水回歸系數(shù)取15%，灌溉引水平均含沙量據(jù)調(diào)查取10 kg/m3。
表2 無定河流域水保措施減水減沙效益
2 結(jié)果與分析
2.1 趨勢與變點(diǎn)分析
對無定河流域年降水量、年潛在蒸散發(fā)、年徑流量和年輸沙量進(jìn)行趨勢分析，見圖2。結(jié)果表明無定河流域的年降水量（圖2a）變化趨緩，其5a滑動平均值在系列均值附近波動；而流域年潛在蒸散發(fā)（圖2b）、年徑流量（圖2c）和年輸沙量（圖2d）均呈明顯減少趨勢。
本研究進(jìn)一步采用M-K檢驗(yàn)、Pettitt檢驗(yàn)分別進(jìn)行系列一致性分析和突變分析，結(jié)果見表3。結(jié)果表明：年降水量呈不顯著增加趨勢，年潛在蒸散發(fā)、年徑流量及年輸沙量均呈極顯著減少趨勢；年潛在蒸散發(fā)突變年份為1982年，年徑流量和年輸沙量的突變年份均為1971年。上述分析結(jié)果與牟夏等[26]、任宗萍等[27]分析結(jié)果一致。
2.2 水文法模型優(yōu)選
不同水文法模型擬合結(jié)果見表4。徑流方面，除模型1以外的3個水文法模型納什效率系數(shù)NSE和相對誤差RE基本一致；輸沙方面，模型1、2、3的NSE和RE基本一致，而模型4表現(xiàn)相對較好。綜合來看，模型4的NSE更接近于1且RE更小，即模擬效果相對更好，故本研究采用模型4代表水文法對無定河流域水沙變化進(jìn)行歸因分析。
圖2 無定河流域水文要素年際變化曲線
表3 無定河流域水文要素M-K檢驗(yàn)與Pettitt檢驗(yàn)計(jì)算結(jié)果
注：、分別表示M-K檢驗(yàn)、Pettitt檢驗(yàn)統(tǒng)計(jì)量；***為99.9%置信水平。
Note:andvalue represent Mann-Kendall test and Pettitt test statistics respectively. *** is 99.9% confidence level.
表4 無定河水文法模型參數(shù)率定及優(yōu)選
2.3 歸因?qū)Ρ确治?/h3>
在無定河流域內(nèi)，分別采用彈性系數(shù)法、雙累積曲線法、水文法和水保法計(jì)算氣候變化和人類活動對流域徑流輸沙變化的貢獻(xiàn)率，其中基準(zhǔn)期為1957—1971年。結(jié)果見表5。總體來看，在無定河流域，人類活動為徑流輸沙變化的主導(dǎo)因素，其對輸沙減少的作用較對徑流減少的作用大。人類活動中灌溉引水和造林對流域徑流變化起主導(dǎo)作用；淤地壩和造林對流域輸沙變化起主導(dǎo)作用。淤地壩對流域徑流減少的效益較小，但對流域輸沙量的減少占據(jù)了主導(dǎo)作用，貢獻(xiàn)率達(dá)43.3%。
徑流方面，雙累積曲線法、水文法、彈性系數(shù)法的計(jì)算結(jié)果較為一致，水保法對人類活動的計(jì)算結(jié)果相對偏?。惠斏撤矫?，4種方法的結(jié)果均表現(xiàn)出人類活動占據(jù)主導(dǎo)作用，具體為雙累積曲線法與彈性系數(shù)法的計(jì)算結(jié)果較為一致，水文法和水保法的計(jì)算結(jié)果較為一致。
表5 無定河流域年徑流輸沙變化歸因分析
進(jìn)一步對無定河流域內(nèi)各時期氣候變化和人類活動對水沙變化的影響進(jìn)行分析，見圖3。人類活動對流域徑流、輸沙變化的貢獻(xiàn)率（圖3b、圖3d）明顯較氣候變化對徑流、輸沙變化的貢獻(xiàn)率（圖3a、圖3c）大，且總體來看人類活動對徑流、輸沙變化的貢獻(xiàn)率呈增加的趨勢。人類活動對徑流、輸沙的貢獻(xiàn)率在1990年左右出現(xiàn)一低谷，同時該低谷對徑流的影響（圖3b）較對輸沙的影響（圖3d）大，即與對輸沙的影響相比，人類活動對徑流的影響進(jìn)一步減弱。
圖3 無定河流域水沙變化不同方法對比分析
3 討 論
本研究雙累積曲線法僅考慮年降水量單一自變量，并將其與徑流輸沙間的關(guān)系概化為線性，這在水保措施少、治理標(biāo)準(zhǔn)低、大雨產(chǎn)大沙特征明顯的基準(zhǔn)期有較好的適用性[28]，但并未考慮黃土高原場次降雨中產(chǎn)流產(chǎn)沙的侵蝕性降雨[29]、降雨年內(nèi)分布與基準(zhǔn)相比發(fā)生的變化[30]及坡面、溝道治理措施的持續(xù)開展和減水減沙作用的滯后顯現(xiàn)，即存在年降水量之外的因素對徑流、輸沙產(chǎn)生著不可忽視的作用，或造成計(jì)算結(jié)果中人類活動對水沙變化的貢獻(xiàn)率計(jì)算偏小。為提高突變后研究期模擬值的準(zhǔn)確發(fā)生，一種是在降雨-徑流、降雨-輸沙關(guān)系中引入一個反映年內(nèi)降雨分配的參數(shù)，使關(guān)系式的物理意義更加清晰，如本研究中采用的水文法；另一種是從流域內(nèi)長時段水量平衡原理出發(fā)，考慮各水沙變化潛在影響因素并計(jì)算相關(guān)貢獻(xiàn)率，如本研究中采用的彈性系數(shù)法。在降雨-徑流歸因分析結(jié)果中，見表5，可發(fā)現(xiàn)水文法和彈性系數(shù)法的結(jié)果基本一致，與雙累積曲線法相比結(jié)果更為合理。水保法尚無法考慮影響流域徑流輸沙變化的所有因素，如不同措施之間的耦合作用、淤地壩淤積比及其異地減沙作用、流域內(nèi)極端事件（干旱或暴雨）對水利水保措施的影響、流域內(nèi)社會經(jīng)濟(jì)發(fā)展帶來的道路建設(shè)及取用河沙等。
人類活動對徑流、輸沙貢獻(xiàn)率在1990年左右出現(xiàn)低谷的原因在于流域內(nèi)20世紀(jì)70、80年代修建完成的大部分中小型淤地壩因攔蓄壽命多小于20 a，至1990年左右時已淤滿或已失效[31]，而退耕還林還草政策、淤地壩作為之一的水利部“亮點(diǎn)”工程等尚未出臺實(shí)施。而低谷對徑流影響較輸沙大的原因在于淤地壩淤滿之后雖失去了攔蓄泥沙的庫容，但因攔蓄泥沙而抬高的壩地有效增大了溝道過水?dāng)嗝?，使攜帶的泥沙沿程沉降，而徑流則可以越過淤滿的壩頂進(jìn)入河道。流域內(nèi)輸沙減少幅度較徑流的大，水利水保措施對下墊面的改變有效減緩或阻止了產(chǎn)匯流的形成，而下滲的徑流又促進(jìn)了林草措施進(jìn)一步發(fā)揮作用。
4 結(jié) 論
本研究選取黃河中游黃土高原地區(qū)典型流域無定河流域，以其年降水量、年潛在蒸散發(fā)、年輸沙量和年徑流量及水利水保措施為研究對象，通過趨勢突變分析和歸因?qū)Ρ确治?，得出不同方法下氣候變化和人類活動對徑流輸沙的貢獻(xiàn)率，結(jié)果表明：
1）通過M-K檢驗(yàn)和Pettitt檢驗(yàn)得出，無定河流域1957—2010年年降水量呈不顯著增加趨勢；年徑流量和年輸沙量呈顯著減少趨勢，其突變年份均為1971年。
2）不同方法分析結(jié)果表明，徑流方面，水文法和彈性系數(shù)法物理意義較雙累積曲線法更為清晰，三者計(jì)算結(jié)果較水保法更為合理；輸沙方面，4種方法分析結(jié)果均為人類活動占據(jù)主導(dǎo)作用，區(qū)別在于對人類活動的水沙變化貢獻(xiàn)率具體大小方面，即水保法和水文法結(jié)果較為一致，彈性系數(shù)法和雙累積曲線法結(jié)果較為一致。
3）不同徑流輸沙影響因素貢獻(xiàn)率計(jì)算方法的結(jié)果表明，人類活動占據(jù)主導(dǎo)作用且其作用逐漸增強(qiáng)。水利水保措施對輸沙變化的作用較對徑流變化的強(qiáng)。人類活動中，徑流方面，灌溉引水是一項(xiàng)重要因素；輸沙方面，淤地壩是短期內(nèi)減少流域輸沙的有效措施，而20世紀(jì)集中修建的淤地壩的淤滿或失效在短期內(nèi)將會削弱水利水保措施的減水減沙效益。長遠(yuǎn)來看，流域內(nèi)持續(xù)開展的退耕還林還草等生態(tài)建設(shè)對徑流、輸沙的影響將逐漸發(fā)揮主導(dǎo)作用。
[1]劉國彬，上官周平，姚文藝，等. 黃土高原生態(tài)工程的生態(tài)成效[J]. 中國科學(xué)院院刊，2017，32(1)：11－19. Liu Guobin, Shangguan Zhouping, Yao Wenyi, et al. Ecological effects of soil conservation in Loess Plateau[J]. Bulletin of Chinese Academy of Sciences, 2017, 32(1): 11－19. (in Chinese with English abstract)
[2]Wang F, Mu X, Li R, et al. Co-evolution of soil and water conservation policy and human–environment linkages in the Yellow River Basin since 1949[J]. Science of The Total Environment, 2015, 508: 166－177.
[3]Zhao Y, Cao W, Hu C, et al. Analysis of changes in characteristics of flood and sediment yield in typical basins of the Yellow River under extreme rainfall events[J]. Catena, 2019, 177: 31－40.
[4]Kong D, Miao C, Wu J, et al. Impact assessment of climate change and human activities on net runoff in the Yellow River Basin from 1951 to 2012[J]. Ecological Engineering, 2016, 91: 566－573.
[5]Zheng H, Miao C, Wu J, et al. Temporal and spatial variations in water discharge and sediment load on the Loess Plateau, China: A high-density study[J]. Science of the Total Environment, 2019, 666: 875－886.
[6]高海東，李占斌，李鵬，等. 基于土壤侵蝕控制度的黃土高原水土流失治理潛力研究[J]. 地理學(xué)報(bào)，2015，70(9)：1503－1515. Gao Haidong, Li Zhanbin, Li Peng, et al. The capacity of soil loss control in the Loess Plateau based on soil erosion control degree[J]. Acta Geographica Sinica, 2015, 70(9): 1503－1515. (in Chinese with English abstract)
[7]Wang S, Fu B, Piao S, et al. Reduced sediment transport in the Yellow River due to anthropogenic changes[J]. Nature Geoscience, 2015, 9: 38－41.
[8]王萬忠，焦菊英. 黃土高原侵蝕產(chǎn)沙強(qiáng)度的時空變化特征[J]. 地理學(xué)報(bào)，2002，57(2)：210－217. Wang Wanzhong, Jiao Juying. Temporal and spatial variation characteristics of erosion and sediment yield in the Loess Plateau[J]. Acta Geographica Sinica, 2002, 57(2): 210－217. (in Chinese with English abstract)
[9]張鑫，蔡煥杰，尹曉楠. 應(yīng)用重標(biāo)度極差分析法（R/S）分析無定河流域水沙變化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2010，26(增刊2)：212－217. Zhang Xin, Cai Huanjie, Yin Xiaonan. Variation trends analysis of runoff and sediment time series based on R/S method in Wuding River Basin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(Supp.2): 212－217. (in Chinese with English abstract)
[10]Liang W, Bai D, Jin Z, et al. A study on the streamflow change and its relationship with climate change and ecological restoration measures in a sediment concentrated region in the Loess Plateau, China[J]. Water Resources Management, 2015, 29(11): 4045－4060.
[11]孫倩，于坤霞，李占斌，等. 黃河中游多沙粗沙區(qū)水沙變化趨勢及其主控因素的貢獻(xiàn)率[J]. 地理學(xué)報(bào)，2018，73(5)：945－956. Sun Qian, Yu Kunxia, Li Zhanbin, et al. The trends of streamflow and sediment and their driving factors in the middle reaches of the Yellow River[J]. Acta Geographica Sinica, 2018, 73(5): 945－956. (in Chinese with English abstract)
[12]Han J, Gao J, Luo H. Changes and implications of the relationship between rainfall, runoff and sediment load in the Wuding River basin on the Chinese Loess Plateau[J]. Catena, 2019, 175: 228－235.
[13]陳思. VIC大尺度陸面水文模型在無定河流域的應(yīng)用研究[D]. 北京：中國地質(zhì)大學(xué)(北京)，2018. Chen Si. The Application of Variable Infiltration Capacity (VIC) Macroscale Hydrologic Model in Wuding River Basin[D]. Beijing: China University of Geosciences(Beijing), 2018. (in Chinese with English abstract)
[14]Shi C, Zhou Y, Fan X, et al. A study on the annual runoff change and its relationship with water and soil conservation practices and climate change in the middle Yellow River basin[J]. Catena, 2013, 100: 31－41.
[15]賈紹鳳，梁季陽. 黃土高原降雨、徑流、產(chǎn)沙相互關(guān)系的研究[J]. 水土保持學(xué)報(bào)，1992(3)：42－47. Jia Shaofeng, Liang Jiyang. The correlation among rainfall, runoff and sediment generation in the Loess Plateau[J]. Journal of Soil and Water Conservation, 1992(3): 42－47. (in Chinese with English abstract)
[16]Li L, 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 Processes, 2007, 21: 3485－3491.
[17]湯立群，陳國祥. 水土保持減水減沙效益計(jì)算方法研究[J]. 河海大學(xué)學(xué)報(bào)：自然科學(xué)版，1999(1)：82－87. Tang Liqun, Chen Guoxiang. Study on calculation method of water and soil conservation water and sediment reduction benefit[J]. Journal of Hohai University(Natural Sciences), 1999(1): 82－87. (in Chinese with English abstract)
[18]張經(jīng)濟(jì)，冀文慧，馮曉東. 無定河流域水沙變化現(xiàn)狀、成因及發(fā)展趨勢的研究[M]//汪崗，范昭. 黃河水沙變化研究（第二卷）. 鄭州：黃河水利出版社，2002. 393－429.
[19]Li Y Z, Liu C M, Zhang D, et al. Reduced runoff due to anthropogenic intervention in the Loess Plateau, China[J]. Water, 2016, 8(458): 1－16.
[20]Sankarasubramanian A, Vogel R M. Hydroclimatology of the continental United States[J]. Geophysical Research Letters, 2003, 30(7):314－326.
[21]Zhang L, Dawes W R, Walker G R. Response of mean annual evapotranspiration to vegetation changes at catchment scale[J]. Water Resources Research, 2001, 37(3): 701－708.
[22]惠養(yǎng)瑜，冀文慧，同信奇，等. 無定河流域水沙變化及其發(fā)展趨勢預(yù)測研究[M]//汪崗，范昭. 黃河水沙變化研究（第一卷：下冊）. 鄭州：黃河水利出版社，2002. 689－714.
[23]左仲國. 清水河流域水沙變化及發(fā)展趨勢預(yù)測[M]//汪崗，范昭. 黃河水沙變化研究（第一卷：上冊）. 鄭州：黃河水利出版社，2002. 439－448.
[24]冉大川，張棟，焦鵬，等. 西柳溝流域近期水沙變化歸因分析[J]. 干旱區(qū)資源與環(huán)境，2016，30(5)：143－149. Ran Dachuan, Zhang Dong, Jiao Peng, et al. Analysis on the attribution of water and sediment changes in Xiliugou basin[J]. Journal of Arid Land Resources and Environment, 2016, 30(5): 143－149. (in Chinese with English abstract)
[25]韓向楠，謝世友，高云飛. 近年無定河流域淤地壩攔沙作用研究[J]. 人民黃河，2018，40(11)：9－12，41. Han Xiangnan, Xie Shiyou, Gao Yunfei. Study on sediment retaining effect of check dams in Wuding River Basin in recent years[J]. Yellow River, 2018, 40(11): 9－12, 41. (in Chinese with English abstract)
[26]牟夏，劉強(qiáng)，平凡，等. 1960—2007年無定河流域徑流情勢變化及其歸因分析[J]. 北京師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2017，53(4)：465－471. Mou Xia, Liu Qiang, Ping Fan, et al. Analysis of changes in annual streamflow and its influence factors in Wuding River basin in the last 50 years (1960—2007)[J]. Journal of Beijing Normal University(Natural Sciences), 2017, 53(4): 465－471. (in Chinese with English abstract)
[27]任宗萍，馬勇勇，王友勝，等. 生態(tài)建設(shè)條件下無定河不同地貌區(qū)徑流變化歸因分析[J]. 生態(tài)學(xué)報(bào)，2019，39(7)：1－9. Ren Zongping, Ma Yongyong, Wang Yousheng, et al. Runoff changes and attribution analysis in tributaries of different geomorphic regions in Wuding River under ecological construction[J]. Acta Ecologica Sinica, 2019, 39(7): 1－9. (in Chinese with English abstract)
[28]李勉，李平，楊二，等. 黃土丘陵區(qū)淤地壩建設(shè)后小流域泥沙攔蓄與輸移特征[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(18)：80－86. Li Mian, Li Ping, Yang Er, et al. Characteristics of sediment retention and transport in small watershed after construction of check dams in Loess Hilly Area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(18): 80－86. (in Chinese with English abstract)
[29]張攀，姚文藝，唐洪武，等. 黃土坡面細(xì)溝形態(tài)變化及對侵蝕產(chǎn)沙過程的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34(5)：114－119. Zhang Pan, Yao Wenyi, Tang Hongwu, et al. Rill morphology change and its effect on erosion and sediment yield on loess slope[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(5): 114－119. (in Chinese with English abstract)
[30]Tang X, Miao C, Xi Y, et al. Analysis of precipitation characteristics on the loess plateau between 1965 and 2014, based on high-density gauge observations[J]. Atmospheric Research, 2018, 213: 264－274.
[31]高云飛，郭玉濤，劉曉燕，等. 陜北黃河中游淤地壩攔沙功能失效的判斷標(biāo)準(zhǔn)[J]. 地理學(xué)報(bào)，2014，69(1)：73－79. Gao Yunfei, Guo Yutao, Liu Xiaoyan, et al. Failure criteria of the warping dams on sediment interception in the Middle Yellow River in northern Shaanxi[J]. Acta Geographica Sinica, 2014, 69(1): 73－79. (in Chinese with English abstract)
Attributional comparative analysis of runoff and sediment change in typical basin of Loess Plateau
Jiang Kaixin1,2, Yu Kunxia1,2※, Cao Wenhong3, Zhang Xiaoming3, Li Peng1,2, Sun Qian1,2, Liu Yu1,2
(1.-,,710048,; 2.,,710048,; 3.,,100048,)
In recent decades, the runoff and sediment of the Yellow River middle reaches located in the Loess Plateau have changed significantly. Many scholars have carried out a lot of researches on the Wuding River, a typical basin of the Loess Plateau, but these researches are focused on the breakthrough of methods and ignore the comparative analysis among methods. This study compared the difference of four methods, i.e., double cumulative curve method, hydrological method, elastic coefficient method and soil and water conservation method, to conduct attributional comparative analysis of runoff and sediment change in the Wuding River Basin from 1957 to 2010. Both double cumulative curve method and hydrological method separate the influence of climate change and human activity by substituting the precipitation data in the research period into the precipitation-runoff and precipitation-sediment transport relationships established in the reference period, and the hydrological method takes into account the annual distribution of precipitation. The elastic coefficient method generates runoff and sediment elasticities to quantify the effects of climate change and human activities. The soil and water conservation method refers to determining the runoff and sediment reduction indicators of each measure under different conditions and periods in different geomorphic areas. The results of trend and change point analysis showed that the annual precipitation in the Wuding River Basin has not changed significantly. The annual potential evaporation, annual runoff and annual sediment transport decreased significantly and the year of change point was 1982, 1971 and 1971, respectively. The attributional comparative analysis of runoff showed that the results of hydrological method were close to those of elastic coefficient method, and the contribution rate of human activities of soil and water conservation method was relatively small, and the attributional comparative analysis of sediment transport showed that the results of soil and water conservation method were consistent with those of the other three methods. Attribution analysis results indicated human activities were the dominant factor of runoff and sediment change in the basin, and their effects on runoff and sediment change were increasing year by year. Among different human activities, irrigation water diversion and afforestation played a leading role in runoff change, and check dam and afforestation played a leading role in sediment transport change. The full up or failure of check dam built in the 1970s and 1980s would largely weaken the runoff and sediment reduction benefits of water conservancy and soil conservation measures in the short term. In the long term, the impact of ecological construction, such as returning cultivated land to forestland or grassland, on runoff and sediment transport will gradually play a major role. Therefore, in the future, the governance of the Wuding River Basin should continue to implement the policy of conversion of cropland to forest and grassland, and strengthen the layout optimization and post-maintenance of engineering measures such as check dam.
runoff; watershed; sediment transport; Loess Plateau; Wuding River basin; attributional comparative analysis
蔣凱鑫，于坤霞，曹文洪，張曉明，李鵬，孫倩，劉昱. 黃土高原典型流域水沙變化歸因?qū)Ρ确治鯷J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(4)：143－149. doi：10.11975/j.issn.1002-6819.2020.04.017 http://www.tcsae.org
Jiang Kaixin, Yu Kunxia, Cao Wenhong, Zhang Xiaoming, Li Peng, Sun Qian, Liu Yu. Attributional comparative analysis of runoff and sediment change in typical basin of Loess Plateau[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(4): 143－149. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.04.017 http://www.tcsae.org
2019-07-08
2019-10-16
國家自然科學(xué)基金面上項(xiàng)目（51879281，41877077）；國家重點(diǎn)研發(fā)計(jì)劃（2016YFC0402407）
蔣凱鑫，博士生，主要從事流域侵蝕動力學(xué)研究。Email：jiangkaizn@qq.com
于坤霞，博士，副教授，主要從事水文水資源研究。Email：yukunxia@126.com
10.11975/j.issn.1002-6819.2020.04.017
S157.1
A
1002-6819(2020)-04-0143-07
猜你喜歡
水利水電快報(bào)(2022年8期)2022-11-23
農(nóng)業(yè)工程學(xué)報(bào)(2022年10期)2022-08-22
黑龍江大學(xué)自然科學(xué)學(xué)報(bào)(2022年1期)2022-03-29
現(xiàn)代農(nóng)村科技(2021年9期)2021-12-06
山西水土保持科技(2021年4期)2021-12-05
農(nóng)業(yè)科技與信息(2020年15期)2020-12-18
魅力中國(2020年39期)2020-12-08
品牌研究(2020年32期)2020-08-09
人民珠江(2019年4期)2019-04-20
愛你(2018年22期)2018-11-14
