李海防，衛(wèi) 偉，陳利頂，郭二輝，黃 勇

（1.中國科學(xué)院 生態(tài)環(huán)境研究中心 城市與區(qū)域生態(tài)國家重點(diǎn)實(shí)驗(yàn)室，北京100086；2.桂林理工大學(xué) 旅游學(xué)院，廣西 桂林541004；3.環(huán)境保護(hù)部 環(huán)境工程評估中心，北京100012）

在黃土高原，水資源是影響當(dāng)?shù)厣鷳B(tài)環(huán)境改善的主要制約因子。黃土高原土層深厚，結(jié)構(gòu)疏松，持水能力強(qiáng)，是巨大的天然蓄水庫［1］。土壤儲水既有收入，也有支出，處于動態(tài)平衡狀態(tài)。收入項(xiàng)主要以大氣降雨為主，部分地區(qū)和地段還包括少量地下水補(bǔ)給和凝結(jié)水；支出項(xiàng)包括植被蒸散耗水、冠層截留、水分滲漏、地表和地下徑流等；水量收入和支出的差值反映了土壤儲水量的變化［2］。目前，黃土高原林草地總面積約8.78萬km2［3－4］，近十幾年來，通過“退耕還林還草”工程措施的實(shí)施，黃土高原林草地面積有一定程度的提高，但由于不同植被類型和植被條件下，土壤性質(zhì)及土壤水分的利用狀況不同，從而影響土壤水分循環(huán)，進(jìn)而又對地上植被的生長產(chǎn)生重要影響［5－6］。大面積恢復(fù)林草地改變了土壤的水量平衡，導(dǎo)致土壤干層擴(kuò)張［3，7］及“小老樹”現(xiàn)象的產(chǎn)生［8］。因而，研究黃土高原林草地土壤水量平衡，對指導(dǎo)地方林業(yè)生態(tài)建設(shè)，改善地區(qū)生態(tài)環(huán)境，具有重要的理論和實(shí)踐意義［2］。

1 黃土高原土壤水分入滲

1.1 降雨時空變化

降雨是黃土高原土壤水分的唯一來源，其分布具有明顯的時空分異性。從大的空間范圍看，黃土高原降雨量從東南向西北方向逐漸減少，年平均降雨量從東南部的550mm遞減到西北部的不足200mm［4］。黃土高原大部分降雨集中于夏季7—9月份，約占全年總降雨量的50%～80%。同時，黃土高原年均降雨量在最近幾年呈現(xiàn)出逐年減少的趨勢［9－10］。

1.2 林草地土壤水分入滲過程

黃土高原大部分地區(qū)屬壤質(zhì)土壤，土層深厚，除因雨強(qiáng)、坡度等因素造成超滲產(chǎn)流外，大部分降雨可以滲入土中，成為土壤水［3］。有關(guān)黃土高原土壤水分入滲過程，國內(nèi)外學(xué)者在不同立地條件下進(jìn)行了大量研究。劉賢趙等［11］研究黃土高原溝壑區(qū)土壤入滲過程認(rèn)為，土壤積水入滲過程中，土壤含水量一般經(jīng)過“穩(wěn)定不變、緩慢上升、急劇上升和再穩(wěn)定”4個階段；不同土地利用類型各階段經(jīng)歷的時間長短不同，積水越深，土壤含水率、入滲變化越明顯，濕潤鋒推移的速度越快。一般認(rèn)為，土壤初始含水率直接影響著雨后土壤水分入滲。Bodman和Colinan［12］認(rèn)為，在入滲初期，隨著土壤初始含水率的增加，土壤入滲速率逐漸減少；隨著時間的延續(xù)，土壤水分趨于飽和，含水率對入滲的影響也逐漸變小，最終可以忽略。解文艷和樊貴盛［13］也認(rèn)為，土壤入滲量隨土壤含水量的增加而減?。挥绊懲寥罎B透性的因子首先是非毛管孔隙度，其次是坡度和水穩(wěn)性團(tuán)聚體，而與土壤容重和總孔隙度相關(guān)性很小［6，14］。在降雨動力推動下，一部分雨水可以下滲到一定土層范圍內(nèi)［15］，另一部分雨水盡管有重力推動，但往往由于上層土壤很難達(dá)到水分飽和狀態(tài)，雨水很難進(jìn)一步向深層土壤入滲［3］。李玉山等［16］的研究也指出，降雨對黃土高原土壤水分的影響一般可以達(dá)到200cm以下，在大雨條件下可以達(dá)到500cm。

1.3 地上植被與土壤水分入滲的關(guān)系

關(guān)于地上植被與水分入滲的關(guān)系，余新曉等［17］對不同蓋度的油松林進(jìn)行研究發(fā)現(xiàn)，林分蓋度越高，土壤入滲性能越好，水土保持效益越顯著；地表植被覆蓋越低，地表徑流越大；減少草地增加林地，能有效降低地表徑流。李斌兵等［18］、潘成忠和上官周平［19］研究了坡耕地、草地和林地的降雨入滲率發(fā)現(xiàn)，不同植被的入滲率排序一般為：林地＞草地＞坡耕地。韓冰等［20］也發(fā)現(xiàn)，油松林的初滲率、平均入滲速率和穩(wěn)滲率均顯著高于荒坡地，具有較強(qiáng)的入滲能力。劉賢趙等［21］也得出類似的結(jié)論，林地土壤穩(wěn)滲率是農(nóng)地的2.7倍；初滲率是撂荒地的3.3倍。林草地對土壤物理性質(zhì)的影響主要表現(xiàn)在植被能有效地改善0—20cm層的土壤結(jié)構(gòu)。與荒坡地相比，林草地土壤表層容重較小，有效地改善了表層土壤的結(jié)構(gòu)［22］。林地地表積累大量枯枝落葉，這些凋落物的分解，增加了土壤有機(jī)質(zhì)的含量，使土壤容重降低、孔隙度增大，增加了土壤的入滲能力。農(nóng)田與荒草地的土壤特性相似，但苜蓿地情況較為特殊，由于苜蓿具有強(qiáng)大的根系，其地表土壤特性與林地相似，就整個土壤剖面而言，土壤容重大小次序?yàn)椋很俎５兀玖值兀巨r(nóng)田＞荒草地［23］。人工林可以通過增加土壤入滲減少地表徑流量，張曉明等［24］研究認(rèn)為，通過擴(kuò)大林草地，治理流域比未治理流域在豐水年、平水年和枯水年的徑流系數(shù)分別減少50%、85%和90%。劉昌明等［25］指出，黃土高原林區(qū)的年徑流深顯著低于非林區(qū)，林區(qū)的徑流系數(shù)較非林區(qū)小40%～60%，非林區(qū)的年徑流量為林區(qū)的1.7～3.0倍。但Chen等［26］發(fā)現(xiàn)，由于地表堅(jiān)實(shí)度較高和缺少地表覆蓋物，相對農(nóng)田和草地，林地的地表徑流量較大。草地也能夠明顯改善土壤的滲透性能，在干旱半干旱地區(qū)，大力種植草本植物，能起到增加入滲改善地表徑流的作用，如在黃土半干旱區(qū)陡坡地通過林草混交，可提高植被覆蓋度和林下植被生物量，大大降低地表徑流［17］。就混交林和純林而言，陳杰等［27］認(rèn)為，混交林的土壤容重均略小于純林，這主要是由于混交林密度較大，活地被物及根系較多，對表層土壤起到一定的疏松作用，土壤容重降低，孔隙度增大，但其改良作用十分有限?？傊?，水分進(jìn)入土壤的過程是一個復(fù)雜的水文過程，它與枯落物覆蓋、根系分布、表土結(jié)構(gòu)、土壤容重、濕度、坡向和坡位等多種因素有關(guān)［20］。例如，李斌兵和鄭粉莉［18］研究認(rèn)為，林地的入滲率對土壤飽和導(dǎo)水率變化比較敏感，草地的入滲率對土壤初始含水率變化比較敏感。

2 黃土高原土壤水分蒸散發(fā)

2.1 裸地土壤水分的蒸發(fā)

黃土高原土壤水分蒸散發(fā)是指一定時間內(nèi)植物蒸騰量與土壤蒸發(fā)量之和，是土壤水分的主要支出項(xiàng)。尤其是當(dāng)土壤側(cè)滲和地表徑流很少發(fā)生時，蒸散發(fā)就成為土壤水分最主要的支出方式［14］。Chen等［15］研究認(rèn)為，地表蒸散可以消耗掉土壤入滲的大部分水分，即使在降雨較多的年份也可以消耗掉土壤超過80%的水分。當(dāng)?shù)乇頉]有植被而處于裸露狀態(tài)時，土壤水分的主要支出項(xiàng)則為土壤蒸發(fā)。黃土高原降水量少，地表蒸發(fā)量大，年蒸發(fā)量可達(dá)800～1 000 mm［2］。目前，有關(guān)黃土高原裸地土壤水分蒸發(fā)的研究較少。一般認(rèn)為，裸地土壤水分蒸發(fā)較大，土壤水分有效性比有植被覆蓋的低［28］。黃土高原土壤持水性與含中、輕和沙壤土的質(zhì)地有關(guān)，但不受重黏土的影響［29］。

2.2 林草地土壤水分的蒸騰

蒸騰作用是植物消耗水分的主要途徑，由于植物種類的不同以及地區(qū)間的差異，不同植被土壤水分的虧缺量不同，但總體上林地和草地土壤水分的蒸騰量較高［14］。胡夢珺等［2］對幼齡刺槐林的研究表明，在降水量為416.2mm的情況下，長在陰坡的刺槐林年總耗水量為685.6mm，長在陽坡的刺槐林年總耗水量達(dá)752mm。楊海軍等［30］研究晉西水土保持林地土壤的水量平衡認(rèn)為，刺槐、油松等喬木林的蒸散量超過同期降水量，易發(fā)生土壤干旱；沙棘、虎榛子等灌木林發(fā)生干旱的危險性次之，草地與裸地的危險性最小。王進(jìn)鑫等［31］對黃土高原丘陵區(qū)刺槐、側(cè)柏和油松林地土壤水分的研究結(jié)果表明，由于植被的蒸騰作用，0—40cm層的土壤水分虧缺最大，刺槐、側(cè)柏和油松分別為72.22～75.36mm、62.70～65.82mm和69.56～71.08mm。張建軍等［32］研究人工林和次生林土壤蓄水量，結(jié)果表明：0—150cm土層中次生林地蓄水量是人工林地的114倍，刺槐與油松人工林地的耗水量明顯大于次生林。

關(guān)于草地土壤水分的蒸騰，李玲芬等［33］在研究草地的土壤水分動態(tài)時發(fā)現(xiàn)，土壤水分含量與植物耗水密切相關(guān)；多年生植被生長期長，根系分布深，年蒸發(fā)蒸騰量大于一年生植被。Zhao等［34］研究了渭北人工草地的水量平衡，結(jié)果表明，紫花苜蓿人工草地的蒸散量很大，即使在豐水年也能消耗深層的土壤儲水，處在生長發(fā)育期的紫花苜蓿年耗水量高達(dá)800～1 000mm，為年降水量的1.4～1.7倍；3a生大葉苜蓿在年降水量為255.8mm的情況下，其年蒸騰量達(dá)312.6mm。在同樣的氣候、地形與土壤條件下，多年生自然牧草地的蒸發(fā)量和蒸騰量在植物旺盛生長時期，也明顯大于年降水量［35］。與林地內(nèi)的土壤水分相比，草地土壤含水量要明顯偏高，且隨深度增加呈增加趨勢，而林地則正好相反［33］。

2.3 地上植被和土壤水分蒸散的關(guān)系

在黃土高原地區(qū)，植被對土壤水分含量有雙重影響。一方面，林草植被的攔蓄降水、截留地表徑流等功能以及植被對地表的直接或間接覆蓋，能減小裸地蒸發(fā)；林草地的水文功能主要表現(xiàn)在樹冠截留、樹干滯流、林下植被及枯枝落葉層滯流和增加土壤入滲［36］。另一方面，植被的蒸騰作用會加強(qiáng)土壤水分向空氣中散失。大量資料表明，黃土高原林草地土壤水量支出往往大于收入，土壤水分嚴(yán)重虧缺，致使后期植被生長緩慢乃至衰敗枯死，土壤水分是影響地上植被的最重要的環(huán)境因子［37］。因而，黃土高原植被和土壤水分蒸散的關(guān)系研究，關(guān)系到生態(tài)環(huán)境和植被群落的可持續(xù)發(fā)展［1］。

植被耗水量與樹種構(gòu)成、年齡大小及生物量密切相關(guān)，不同的植被類型，高生長量的人工牧草和木本植物，根系發(fā)達(dá)，耗水量大。不同林型的植被耗水量一般為喬木林＞灌木林＞草地＞裸地［38－40］。不同的土地利用方式土壤水分含量一般為：農(nóng)田＞荒草地＞林地＞苜蓿地［23］，其土壤水分的季節(jié)變化趨勢一致，但苜蓿地和林地土壤水分含量的季節(jié)變化較農(nóng)田和荒草地更為劇烈。就林地而言，不同植被類型的生理耗水量不同，土壤水分虧缺程度不同，虧缺次序?yàn)椋簷帡l＞刺槐＞苜蓿＞側(cè)柏＞楊樹＞油松＞荒坡＞杏＞棗＞農(nóng)地［41］。李軍等［42］對黃土高原53類林地、草地和農(nóng)田深層土壤濕度的觀測結(jié)果表明：林地、草地和農(nóng)田土壤水分過耗量平均值分別為521，491，30 mm，林地、草地和農(nóng)田土壤干燥化指數(shù)平均值分別為39%，42%和96%。從草地到林地，土壤可補(bǔ)給水量減少了25%～45%［43－44］，根系越深，土壤水分蒸散越多［45－46］。Wang等［47］研究發(fā)現(xiàn)，苜蓿、檸條和油松分別可以吸收地下15.5，22.4，21.5m 深的水分。農(nóng)田土壤失水層可達(dá)地下500cm，而果園和草地分別為940cm和800～1 000cm［48］。樊軍等［49］研究發(fā)現(xiàn)，陜西省長武縣15a生的蘋果園和15a生的紫花苜蓿地，土壤干層比10a生的蘋果園和7a生的刺槐林地要深。邵薇薇等［50］以黃土高原51個小流域?yàn)檠芯康攸c(diǎn)，分析植被指數(shù)（NDVI）與水分蒸散的關(guān)系。結(jié)果表明：在降水量較多的地區(qū)，植被密度越大，實(shí)際蒸散發(fā)量也越多。但有研究發(fā)現(xiàn)：林地的生物量與總耗水量相關(guān)性不顯著，而與林分蒸騰耗水顯著相關(guān)。例如，紫花苜宿由于其根系強(qiáng)大，具有很強(qiáng)的水分蒸騰能力［51］。也有研究證實(shí)，不同林地類型總耗水量接近，而蒸騰耗水量差別較大；相同地點(diǎn)、同種植被條件，不同生物量的林分蒸騰耗水量差別為69～365mm，而蒸散耗水為1～400mm；凡是蒸騰耗水量少的林分，其蒸散耗水量必定較大［52－53］。葉面積指數(shù)是影響黃土高原植被蒸騰的主要因子［54－56］。25a果齡的果園100—400cm土層土壤含水量明顯低于15a果齡果園［57］。植被由低級階段向高級階段演替過程是“土壤水庫”的耗水過程［58］。此外，土壤水分的蒸散還受降雨量、氣溫、濕度等因素的影響，耗水量也有很大差別。

3 黃土高原土壤干層

3.1 黃土高原土壤干層的形成

土壤干層是指位于降雨入滲補(bǔ)給深度以下，因植被蒸散導(dǎo)致土壤水分負(fù)平衡而形成的干土層。它具有以下主要特征：（1）位于土體某一深度范圍內(nèi)，一般在降水入滲深度以下；（2）具有相對持久性，不會因?yàn)橥恋乩妙愋偷母淖兒徒涤甑娜霛B而消失；（3）具有一定的濕度范圍，下限為凋萎濕度，上限為毛管斷裂濕度或土壤穩(wěn)定濕度［47］；（4）干層深度與植物根系分布相對應(yīng)，根系越長，干層越深；（5）干燥強(qiáng)度因植物種類和生長年限而定，并與降水量和蒸散量的比值相關(guān)［59］。一般認(rèn)為，“低降水、高蒸發(fā)”的氣候環(huán)境是黃土高原土壤干層形成的主要原因［60］。近50a來，黃土高原地區(qū)氣候變暖和降雨減少可能是土壤干層形成的直接原因，而植被類型選擇不當(dāng)、群落密度過大和生產(chǎn)力過高則會加劇深層土壤的干燥化進(jìn)程［61］，有植被的土壤干層要比裸地嚴(yán)重的多［7，15］。

3.2 地上植被特征與土壤干層的關(guān)系

大量研究表明，黃土高原現(xiàn)有的天然林、灌木林和果林普遍存在著土壤干層［62］。研究表明，刺槐、白榆、油松是黃土高原干旱、半干旱區(qū)的主要人工林類型，林分密度普遍偏高，容易出現(xiàn)土壤干化現(xiàn)象［63］。也有研究認(rèn)為，干化層的形成主要受區(qū)域本身的土壤、地理、氣候等因素影響，而非林木生長的結(jié)果［64］。但綜合前人的研究，地上植被與土壤干層的形成是顯著相關(guān)的。杜鵑等［65］研究發(fā)現(xiàn)，楊樹林、法國梧桐林和中國梧桐林下1.5～4m土層發(fā)育了弱的土壤干化層，而麥地和草地下0～6m水分狀況良好，未出現(xiàn)干層現(xiàn)象。萬素梅等［66］發(fā)現(xiàn)生長年限低于8a的苜蓿草地出現(xiàn)輕度干層，生長年限超過8a，出現(xiàn)中度干層，苜蓿生長18a后，深層土壤通體干化。黃土高原典型植被土壤干層厚度表現(xiàn)為：油松、遼東櫟、狼牙刺和檸條林地＞刺槐林地＞蘋果和沙棘林地［42，67］。Wang等研究確定苜蓿、檸條和油松林植被下的耗水深度分別達(dá)到15.5，22.4，21.5m。人工林地和草地由于耗水量大，且具有多年連續(xù)性，往往是年耗水量大于年降水量；由于沒有足夠的降雨補(bǔ)充，深層土壤長時間水分虧缺，即使在豐水年也不可能得到完全恢復(fù)，長此累積則會形成永久性干層，且隨年限增加，干層厚度逐漸加大［68］，并在短時間內(nèi)難以改善［46］。此外，土壤干層并不完全是人工植被的產(chǎn)物，天然林地和天然草地同樣有土壤干層現(xiàn)象發(fā)生［60］。天然荒草地土壤水分循環(huán)強(qiáng)度大于裸地，表現(xiàn)為荒草地蒸騰作用層深度較大，二者分別為200cm和180cm。在相同的降雨量條件下，荒草地土壤水分循環(huán)強(qiáng)度高于裸地，表現(xiàn)為荒草地降雨入滲補(bǔ)給深度較?。?9］。此外，土壤干層的發(fā)育還受地形、坡向、坡度和輪作方式等因素的影響。一般認(rèn)為，陽坡形成的干層較陰坡嚴(yán)重；坡度愈大，土壤干化愈劇烈［70］。與連作苜蓿相比，輪作苜蓿不會大量消耗土壤深層水分而形成深厚的土壤干層，有利于土壤水分的可持續(xù)利用［71］。針闊葉植被配置模式土壤水分狀況要稍好于闊葉純林的配置模式［67］。

4 黃土高原水量平衡模擬

研究黃土高原植物與土壤水分之間的動態(tài)平衡關(guān)系，對于建立穩(wěn)定的植被生態(tài)系統(tǒng)，實(shí)現(xiàn)生態(tài)系統(tǒng)的高效可持續(xù)運(yùn)行，具有重要作用［72］。為解釋土壤干層的形成機(jī)理，圍繞黃土高原水循環(huán)動態(tài)過程，國內(nèi)外研究者開展了大量的建模嘗試和分析。其中，SWAP、CoupModel、WAVES等水量平衡模型應(yīng)用較多。但這些模型中，有些模型需要大量的參數(shù)，應(yīng)用起來比較復(fù)雜；有些模型則比較簡單，對水量平衡動態(tài)的描述又有欠缺［73］。黃土高原較為簡單的水量平衡模型為：

P＝E＋T＋D＋R＋ΔS＋Ei

式中：P——降雨量；E——地表蒸發(fā)；T——植被蒸騰；R——地表徑流；D——水分滲漏；ΔS——土壤儲水量；Ei——林冠穿透雨蒸發(fā)。一般情況下，地表徑流和林冠穿透雨蒸發(fā)量較少，水量平衡模型可簡寫為：P－ΔS＝E＋T＋D［74］。在黃土高原丘陵溝壑區(qū)，陳洪松等［69］（2005）提出簡化的坡面土壤水平衡模型為：

式中：ΔW——觀測期始末土壤儲水量的差值（mm）；Pα——坡面承雨量（mm）；Ps——地表徑流量（mm）；I——冠層截留量（mm）；E——蒸散（發(fā)）量（mm）；P——降落在平坦地面上的雨量（mm）；α——地面坡度（°）；β——雨點(diǎn)著地軌跡線與水平面的夾角（°）。佘冬立等［35］利用土壤—水—大氣—植物耦合模型（Soil－Water－Atmosphere－Plant，SWAP）模擬黃土高原水蝕風(fēng)蝕交錯區(qū)林草地SPAC系統(tǒng)水分傳輸過程，揭示了坡面典型植被SPAC系統(tǒng)水量平衡狀況。SWAP模型則主要用于模擬土壤—植物—大氣傳輸系統(tǒng)中水分運(yùn)動、溶質(zhì)運(yùn)移、熱量傳輸及作物生長過程［75］；CoupModel模型也經(jīng)常被應(yīng)用于水量平衡模擬［76－77］，該模型是一維模型，遵循質(zhì)量守恒定律，由水勢梯度和溫度梯度作為土壤水運(yùn)動的驅(qū)動力，在此基礎(chǔ)上模擬了飽和與非飽和條件下水分的傳導(dǎo)機(jī)制，以及SPAC系統(tǒng)包括地表徑流、水分蒸散、土壤表面的輻射吸收平衡、水分入滲、碳氮循環(huán)等生態(tài)過程［78］；該模型同時融合了SOIL［79］和SOIL—N［80］兩個模型。WAVES模型用四大模塊來模擬能量、降雨入滲、徑流、植物根系水分吸收、水分再分配、地下水和碳循環(huán)等生態(tài)過程，冠層蒸散量和土壤蒸發(fā)用Penman－Monteith公式計(jì)算，動力阻力根據(jù) Monm－Obukhov的相似理論來確定，冠層阻力利用Ball和Leuning模型來確定。而 Huang等［73］則用改進(jìn)的WAVES模型對黃土高原不同降雨條件下降雨、植被、水量平衡之間的聯(lián)系進(jìn)行動態(tài)模擬。

5 研究展望

如何解決黃土高原干旱缺水這一難題，是實(shí)現(xiàn)當(dāng)?shù)厣鷳B(tài)健康和區(qū)域人地關(guān)系協(xié)調(diào)發(fā)展的根本。自1999年國家實(shí)施“退耕還林還草”工程以來，黃土高原進(jìn)行了大面積的植樹種草，生態(tài)環(huán)境建設(shè)取得了一定成效。但同時，大面積植樹種草不可避免地對土壤水量平衡產(chǎn)生了重要影響，引起科學(xué)界的普遍關(guān)注。本文系統(tǒng)回顧了國內(nèi)外主要的相關(guān)研究與進(jìn)展，認(rèn)為未來有關(guān)黃土高原水量平衡方面的研究應(yīng)該在以下兩個方面進(jìn)行深入研究。

第一，加大黃土高原林草地水源涵養(yǎng)功能研究力度。目前，很多有關(guān)水量平衡的研究多集中于考慮植被的蒸騰蒸散對土壤水量平衡的影響，而沒有考慮林草地自身的水源涵養(yǎng)功能，特別是成熟林草地的水源涵養(yǎng)功能對土壤水分蒸散的抵消作用。

第二，林草地生態(tài)水文過程的尺度擴(kuò)展。從已有的研究來看，林草地生態(tài)正向效應(yīng)大于其生態(tài)負(fù)向效應(yīng)，但這個結(jié)論是在研究小區(qū)域生態(tài)水文效應(yīng)基礎(chǔ)上得出的，因此，有必要從大時空尺度研究地表植被對土壤干層的影響，為研究黃土高原土壤水量平衡提供新的內(nèi)容。

［1］ 孫寶勝，楊開寶，拓文俊.黃土高原丘陵溝壑區(qū)土壤水資源平衡利用與生態(tài)植被可持續(xù)發(fā)展［J］.西北農(nóng)業(yè)學(xué)報，2005，14（4）：92－96.

［2］ 胡夢珺，劉文兆，趙姚陽.黃土高原農(nóng)、林、草地水量平衡異同比較分析［J］.干旱地區(qū)農(nóng)業(yè)研究，2003，21（4）：113－116.

［3］ Shangguan Z P，Zheng S X.Ecological properties of soil water and effects on forest vegetation in the Loess Plateau［J］.International Journal of Sustainable Development and World Ecology，2006，13（4）：307－314.

［4］ 何永濤，李文華，郎海鷗.黃土高原降水資源特征與林木適宜度研究［J］.干旱區(qū)研究，2009，26（3）：407－412.

［5］ Priha O，Lehto T，Smolander A.Mycorrhizas and C and N transformations in the rhizospheres of Pinus sylvestris，Picea abies and Betula pendula seedlings［J］.Plant and Soil，1999，206（2）：191－204.

［6］ Cao S X，Chen L，Xu C G，et al.Impact of three soil types on afforestation in China′s Loess Plateau：growth and survival of six tree species and their effects on soil properties［J］.Landscape and Urban Planning，2007，83（2／3）：208－217.

［7］ Chen H S，Shao M A，Li Y Y.Soil desiccation in the Loess Plateau of China［J］.Geoderma，2008，143（1／2）：91－100.

［8］ 侯慶春，韓蕊蓮，韓仕鋒.黃土高原人工林草地“土壤干層”問題初探［J］.中國水土保持，1999（5）：11－14.

［9］ Tao F，Yokozawa M，Xu Y，et al.Climate changes and trends in phenology and yields of field crops in China，1981—2000［J］.Agricultural and Forest Meteorology，2006，138（1／4）：82－92.

［10］ Yang Y，F(xiàn)eng Z，Huang H Q，et al.Climate－induced changes in crop water balance during 1960—2001in Northwest China［J］.Agriculture，Ecosystems and Environment，2008，127（1／2）：107－118.

［11］ 劉賢趙，康紹忠.黃土高原溝壑區(qū)小流域土壤入滲分布規(guī)律的研究［J］.吉林林學(xué)院學(xué)報，1997，13（4）：203－208.

［12］ Bodman G B，Colinan E A.Moisture and energy condition during down ward entry of water into soil［J］.Soil Science，1994，8（2）：166－182.

［13］ 解文艷，樊貴盛.土壤含水量對土壤入滲能力的影響［J］.太原理工大學(xué)學(xué)報，2004，35（3）：272－275.

［14］ 李慶云，余新曉，信忠保，等.黃土高原典型流域不同土地利用類型土壤物理性質(zhì)分析［J］.水土保持研究，2010，17（6）：106－114.

［15］ Chen H S，Shao M A，Li Y Y.The characteristics of soil water cycle and water balance on steep grassland under natural and simulated rainfall conditions in the Loess Plateau of China［J］.Journal of Hydrology，2008，360（1／4）：242－251.

［16］ 李玉山.苜蓿生產(chǎn)力動態(tài)及其水分生態(tài)環(huán)境效應(yīng)［J］.土壤學(xué)報，2002，39（3）：404－411.

［17］ 余新曉，張曉明，武思宏，等.黃土區(qū)林草植被與降水對坡面徑流和侵蝕產(chǎn)沙的影響［J］.山地學(xué)報，2006，24（1）：19－26.

［18］ 李斌兵，鄭粉莉.黃土坡面不同土地利用下的降雨入滲模擬與數(shù)值計(jì)算［J］.干旱地區(qū)農(nóng)業(yè)研究，2008，26（5）：118－223.

［19］ 潘成忠，上官周平.黃土區(qū)次降雨條件下林地徑流和侵蝕產(chǎn)沙形成機(jī)制：以人工油松林和次生山楊林為例［J］.應(yīng)用生態(tài)學(xué)報，2005，16（9）：1597－1602.

［20］ 韓冰，吳欽孝，李秧秧，等.黃土丘陵區(qū)人工油松林地土壤入滲特征的研究［J］.防護(hù)林科技，2004，62（5）：1－3.

［21］ 劉賢趙，黃明斌.黃土丘陵溝壑區(qū)森林土壤水文行為及其對河川徑流的影響［J］.干旱地區(qū)農(nóng)業(yè)研究，2003，21（2）：72－75.

［22］ 張思祖，李凱榮，易亮，等.黃土高原溝壑區(qū)杏樹林地土壤水分時空動態(tài)變化及適宜性研究［J］.水土保持研究，2010，17（5）：21－25.

［23］ 劉春利，邵明安.黃土高原六道溝流域不同土地利用方式下土壤水力特性及其對土壤水分的影響［J］.應(yīng)用生態(tài)學(xué)報，2008，19（11）：2400－2407.

［24］ 張曉明，余新曉，武思宏，等.黃土丘陵溝壑區(qū)典型流域土地利用／土地覆被變化水文動態(tài)響應(yīng)［J］.生態(tài)學(xué)報，2007，27（2）：414－423.

［25］ 劉昌明，曾燕.植被變化對產(chǎn)水量影響的研究［J］.中國水利，2002（10）：112－117.

［26］ Chen L D，Huang Z L，Gong J.The effect of land cover／vegetation on soil water dynamic in the hilly area of the loess plateau，China［J］.Catena，2007，70（2）：200－208.

［27］ 陳杰，劉文兆，張勛昌，等.黃土高原溝壑區(qū)不同樹種的水土保持效益及其適應(yīng)性評價［J］.西北農(nóng)林科技大學(xué)學(xué)報：自然科學(xué)版，2008，36（6）：97－104.

［28］ Latta J，O′Leary G J.Long－term comparison of rotation and fallow tillage systems of wheat in Australia［J］.Field Crops Research，2003，83（2）：173－190.

［29］ Kimura R，F(xiàn)an J，Zhang X C，et al.Evapotranspiration over the Grassland Field in the Liudaogou Basin of the Loess Plateau，China［J］.Acta Oecologica，2006，29（1）：45－53.

［30］ 楊海軍，孫立達(dá)，余新曉.晉西黃土區(qū)水土保持林水量平衡的研究［J］.北京林業(yè)大學(xué)學(xué)報，1993，15（3）：42－50.

［31］ 王進(jìn)鑫，羅偉祥，劉廣全，等.黃土高原人工林根區(qū)土壤水分虧缺狀況與空間分布［J］.西北林學(xué)院學(xué)報，2004，19（4）：1－4.

［32］ 張建軍，張巖，張波.晉西黃土區(qū)水土保持林地的土壤水分［J］.林業(yè)科學(xué)，2009，45（11）：63－69.

［33］ 李玲芬，延軍平，劉冬梅，等.干旱／半干旱地區(qū)不同植被條件下土壤含水量變化及植被建設(shè)途徑分析［J］.水土保持通報，2009，29（1）：18－22.

［34］ Zhao C Y，F(xiàn)eng Z D，Chen G D.Soil water balance simulation of alfalfa（Medicago sativa L.）in the semiarid Chinese Loess Plateau［J］.Agricultural Water Management，2004，69（2）：101－114.

［35］ 佘冬立，邵明安，俞雙恩.黃土高原典型植被覆蓋下SPAC系統(tǒng)水量平衡模擬［J］.農(nóng)業(yè)機(jī)械學(xué)報，2011，42（5）：73－78.

［36］ 李海防，楊章旗，韋理電，等.廣西華山林場5種典型人工林水文功能評價［J］.安徽農(nóng)業(yè)大學(xué)學(xué)報，2011，38（2）：146－150.

［37］ Kimura R，Liu Y，Takayama N，et al.Heat and water balances of the bare soil surface and the potential distribution of vegetation in the Loess Plateau，China［J］.Journal of Arid Environments，2005，63（2）：439－457.

［38］ 尹忠東，朱清科，畢華興，等.黃土高原植被耗水特征研究進(jìn)展［J］.人民黃河，2005，27（6）：35－37.

［39］ Wang Y Q，Shao M A，Zhu Y J，et al.Impacts of land use and plant characteristics on dried soil layers in different climatic regions on the Loess Plateau of China［J］.Agricultural and Forest Meteorology，2011，151（4）：437－448.

［40］ 盧愛剛，董雯，趙景波.陜西長安和曲江人工林地土壤含水量對比研究［J］.水土保持研究，2009，16（3）：155－158.

［41］ 王延平.陜北黃土區(qū)陡坡地土壤水分植被承載力研究［D］.陜西楊凌：西北農(nóng)林科技大學(xué)，2009.

［42］ 李軍，陳兵，李小芳，等.黃土高原不同降水類型區(qū)林地：草地與農(nóng)田土壤干燥化效應(yīng)比較［J］.土壤學(xué)報，2008，45（1）：40－48.

［43］ Bastrup－Birk A，Gundersen P.Water quality improvements from afforestation in an agricultural catchment in Denmark illustrated with the INCA model［J］.Hydrology Earth System Sciences，2004，8（4）：764－777.

［44］ Farley K A，Jobbágy E G，Jackson R B.Effects of afforestation on water yield：aglobal synthesis with implications for policy［J］.Global Change Biology，2005，11（10）：1565－1576.

［45］ Huang M B，Zhang L，Gallichand J.Runoff responses to afforestation in watershed of the Loess Plateau，China［J］.Hydrological Processes，2003，17（13）：2599－2609.

［46］ 穆興民，徐學(xué)選，王文龍，等.黃土高原人工林對區(qū)域深層土壤水環(huán)境的影響［J］.土壤學(xué)報，2003，40（2）：210－217.

［47］ Wang Z Q，Liu B Y，Liu G，et al.Soil water depletion depth by planted vegetation on the Loess Plateau［J］.Science in China Series D：Earth Sciences，2009，52（6）：835－842.

［48］ 李玉山.黃土高原森林植被對陸地水循環(huán)影響的研究［J］.自然資源學(xué)報，2001，16（5）：427－432.

［49］ 樊軍，郝明德，邵明安.黃土旱塬農(nóng)業(yè)生態(tài)系統(tǒng)土壤深層水分消耗與水分生態(tài)環(huán)境效應(yīng)［J］.農(nóng)業(yè)工程學(xué)報，2004，20（1）：61－64.

［50］ 邵薇薇，楊大文，孫福寶，等.黃土高原地區(qū)植被與水循環(huán)的關(guān)系［J］.清華大學(xué)學(xué)報，2009，49（12）：1958－1962.

［51］ Crawford M C，Macfarlane M R.Lucerne reduces soil moisture and increases livestock production in an area of high groundwater recharge potential［J］.Australian Journal of Experimental Agriculture，1995，35（2）：171－180.

［52］ 李代瓊，梁一民，侯喜祿，等.黃土高原沙棘建造植被的生態(tài)功能及效益試驗(yàn)研究［J］.沙棘，2003，16（3）：16－21.

［53］ 程積民，萬惠娥，王靜.黃土丘陵區(qū)山桃灌木林地土壤水分過耗與調(diào)控恢復(fù)［J］.土壤學(xué)報，2003，40（5）：691－694.

［54］ Wang L，Wei S P，Horton R，et al.Effects of vegetation and slope aspect on water budget in the hill and gully region of the Loess Plateau of China［J］.Catena，2011，87（1）：90－100.

［55］ Granier A，Biron P，Lemoine D.Water balance，transpiration and canopy conductance in two beech stands［J］.Agricultural and Forest Meteorology，2000，100（4）：291－308.

［56］ Roberts J.The influence of physical and physiological characteristics of vegetation on their hydrological response［J］.Hydrological Processes，2000，14（16／17）：2885－2901.

［57］ 武阿鋒，劉文兆.長武塬區(qū)不同土地利用條件下土壤水分深剖面分布特征比較［J］.干旱地區(qū)農(nóng)業(yè)研究，2009，27（5）：133－136.

［58］ 龐敏.黃土高原植被建設(shè)的生態(tài)水文效應(yīng)研究［D］.西安：西安理工大學(xué)，2010.

［59］ 陳洪松，王克林，邵明.黃土區(qū)人工林草植被深層土壤干燥化研究進(jìn)展［J］.林業(yè)科學(xué)，2005，41（4）：155－161.

［60］ 陳寶群，趙景波，李艷花.黃土高原土壤干層形成原因分析［J］.地理與地理信息科學(xué)，2009，25（3）：85－89.

［61］ 陳洪松，邵明安，王克林.黃土區(qū)荒草地土壤水平衡的數(shù)值模擬［J］.土壤學(xué)報，2005，42（3）：353－359.

［62］ 王曉賢，張學(xué)培.晉西黃土高原土壤水分垂直變化規(guī)律及生態(tài)影響［J］.水土保持研究，2009，16（2）：200－204.

［63］ 劉晨峰，王正寧，賀康寧，等.黃土高原半干旱區(qū)幾種人工林的土壤水分、光照變化及其對林分的影響［J］.西部林業(yè)科學(xué)，2004，33（3）：34－41.

［64］ 許喜明，陳海濱，原煥英，等.黃土高原半干旱區(qū)人工林地土壤水分環(huán)境的研究［J］.西北林學(xué)院學(xué)報，2006，21（5）：60－64.

［65］ 杜娟，趙景波.黃土高原南部人工植被作用下的土壤水分研究［J］.土壤，2010，42（2）：262－267.

［66］ 萬素梅，賈志寬，韓清芳，等.黃土高原半濕潤區(qū)苜蓿草地土壤干層形成及水分恢復(fù)［J］.生態(tài)學(xué)報，2008，28（3）：1045－1051.

［67］ 楊磊，衛(wèi)偉，陳利頂，等.半干旱黃土丘陵區(qū)人工植被深層土壤干化效應(yīng)［J］.地理研究，2012，31（1）：71－80.

［68］ 張北贏，徐學(xué)選，李貴玉，等.土壤水分基礎(chǔ)理論及其應(yīng)用研究進(jìn)展［J］.中國水土保持科學(xué)，2007，5（2）：122－129.

［69］ 陳洪松，邵明安，王克林，等.黃土區(qū)深層土壤干燥化與土壤水分循環(huán)特征［J］.生態(tài)學(xué)報，2005，25（10）：2491－2498.

［70］ 王力，邵明安，侯慶春，等.黃土高原土壤干層初步研究［J］.西北農(nóng)林科技大學(xué)學(xué)報：自然科學(xué)版，2001，29（4）：34－38.

［71］ 李慧成，郝明德，楊曉，等.黃土高原苜蓿草地在不同種植方式下的土壤水分變化［J］.西北農(nóng)業(yè)學(xué)報，2009，18（3）：141－146.

［72］ 游松財，邸蘇闖.黃土高原地區(qū)1961—2000年間土壤水分變化模擬與分析［J］.自然資源學(xué)報，2010，25（7）：1206－1216.

［73］ Huang M B，Shao M A，Li Y S.Comparison of a modified statistical－dynamic water balance model with the numerical model WAVES and field measurements［J］.Agricultural Water Management，2001，48（1）：21－35.

［74］ Zhang S L，Simelton E，L?vdahl L，et al.Simulated long－term effects of different soil management regimes on the water balance in the Loess Plateau，China［J］.Field Crops Research，2007，100（2／3）：311－319.

［75］Jansson P E，Karlberg L.CoupModel－Coupled heat and Mass Transfer Model for Soil－Plant－Atmosphere System［R］.Stockholm：Royal Institute of Technology，Dept.of Land and Water Resources Engineering，2004.

［76］ Karlberg L，Alon B G，Jansson P E，et al.Modelling transpiration and growth in salinity－stressed tomato under different climatic conditions［J］.Ecological Modelling，2006，190（1／2）：15－40.

［77］ Zhang S L，Lovdahl L，Grip H，et al.Modelling the effects of mulching and fallow cropping on water balance in the Chinese Loess Plateau［J］.Soil ＆ Tillage Research，2007，93（2）：283－298.

［78］ Jansson P E，Moon D.A coupled model of water，heat and mass transfer using object orientation to improve flexibility and functionality［J］.Environmental Modelling and Software，2001，16（1）：37－46.

［79］ Jansson P E，Cienciala E，Grelle A，et al.Simulated evapotranspiration from the Norunda forest stand during the growing season of a dry year［J］.Agricultural and Forest Meteorology，1999，98／99（SI）：621－628.

［80］ Eckersten H，Blomb?ck K，K?tterer T，et al.Modelling C，N，water and heat dynamics in winter wheat under climate change in southern Sweden［J］.Agriculture，Ecosystems ＆ Environment，2001，86（3）：221－235.