黃土丘陵區(qū)小流域典型造林整地工程土壤水分特征曲線模擬

于 洋,衛(wèi) 偉,陳利頂,馮天驕,4,楊 磊,陳 蝶,4

1 中國水利水電科學(xué)研究院流域水循環(huán)模擬與調(diào)控國家重點(diǎn)實(shí)驗(yàn)室,北京 100048 2 中國科學(xué)院生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國家重點(diǎn)實(shí)驗(yàn)室,北京 100085 3 水利部水土保持生態(tài)工程技術(shù)研究中心,北京 100048 4 中國科學(xué)院大學(xué),北京 100049

水分是影響干旱半干旱生態(tài)系統(tǒng)植被恢復(fù)和生態(tài)建設(shè)的主要因素。在黃土高原半干旱小流域,為了實(shí)現(xiàn)有限水資源的高效利用,工程措施廣泛應(yīng)用于植被建設(shè),工程措施的開展也被認(rèn)為是流域減流減沙的主要手段[1- 3]。工程措施一方面通過改變徑流路徑,實(shí)現(xiàn)地表徑流就地?cái)r蓄,另一方面,通過改變微地形減少土壤蒸發(fā)面積,有效保持土壤水分[4]。由于土壤中的水分并不能被植物全部利用,植被生長能夠直接利用的部分為有效水,當(dāng)土壤基質(zhì)勢高于植物體內(nèi)基質(zhì)勢時(shí),水分從土壤進(jìn)入植物體,反之水分保留于土壤中,植物無法利用。水分特征曲線是用于描述不同基質(zhì)勢下土壤能夠保留水分含量的曲線。水土保持工程措施的主要功能就是充分利用有限的水分以維系植被生長,那么在利用工程措施進(jìn)行植被恢復(fù)的背景下,工程措施是否能提高植被對(duì)有效水分的利用？采用配對(duì)樣地實(shí)驗(yàn)的辦法,通過分析土壤水分特征曲線,對(duì)工程措施下植被有效水的體積質(zhì)量與相同植被類型的自然坡面進(jìn)行對(duì)比,是評(píng)價(jià)工程措施改良立地環(huán)境的有效途徑。

水分特征曲線能夠描述土壤水分能量和數(shù)量之間的關(guān)系,其可以用于研究土壤水分有效性和供水體積質(zhì)量的大小,對(duì)研究土壤水的滯留與運(yùn)移有重要作用。當(dāng)前,描述土壤水分特征曲線的經(jīng)驗(yàn)公式較多,比如Brooks-Corey模型(1964)、Gardner模型(1970)、Campbell模型(1974)以及van Genuchten模型(1980)等[5-6];人們基于不同模型在諸多生態(tài)系統(tǒng)開展了土壤水動(dòng)力學(xué)的相關(guān)研究,并通過對(duì)模型的參數(shù)進(jìn)行求解來分析參數(shù)的空間變異性[7-12]。受土壤自身理化性質(zhì)的變化以及外界因素的影響,土壤有效水分也會(huì)發(fā)生改變。在黃土高原,人們對(duì)不同土壤質(zhì)地以及不同土地利用方式下的土壤水力學(xué)特性進(jìn)行了報(bào)道,揭示了影響水分特性變異的主要因素[13-14]。工程措施改變了坡面原始下墊面結(jié)構(gòu),在改變徑流路徑的同時(shí)影響顆粒物的搬運(yùn)和堆積,但圍繞工程措施與植被長時(shí)間耦合條件下的土壤水分特征曲線模擬以及工程措施如何提高土壤有效水分,相關(guān)研究較少。基于此,本研究通過實(shí)測水分特征曲線,并與經(jīng)驗(yàn)?zāi)Ｐ拖嘟Y(jié)合進(jìn)行水分特征曲線模擬,以探討不同水土保持工程措施對(duì)土壤有效水分的影響,并為評(píng)價(jià)工程措施在植被恢復(fù)過程中的保持水土功能提供依據(jù)。

1 研究區(qū)概況

研究區(qū)地處甘肅定西龍灘流域(104°27′—104°31′ E,35°43′—35°46′ N),流域面積16.1 km2,海拔1840—2260 m,多年平均降水量386.3 mm,屬于典型半干旱黃土丘陵區(qū),降水量季節(jié)分配不均,春季降水較為稀少,主要集中在7—9月份。流域年均氣溫為6.8 ℃,平均無霜期152 d,潛在蒸發(fā)量為1649.0 mm,年均相對(duì)濕度72%。流域土壤類型為黃綿土,土質(zhì)均一。流域內(nèi)主要土地利用類型包括撂荒地、灌木林地、喬木林地、農(nóng)地、人工草地。同時(shí)分布著零星的天然草地。流域內(nèi)的人工植被以檸條(Caraganakorshinskii)、側(cè)柏(Platycladusorientalis)、山杏(Prunusarmeniaca)、油松(Pinustabulaeformis)、紫花苜蓿(Medicagosativa)為主,天然植被則以多年生草本為主,主要種類包括賴草(Leymussecalinus)、長芒草(Stipabungeana)、阿爾泰狗娃花(Heteropappusaltaicus)等。自20世紀(jì)50年代起,流域內(nèi)開展了大面積植被恢復(fù)工程,取得了較好的生態(tài)效益。水平溝、水平階、魚鱗坑、反坡臺(tái)等水土保持工程措施廣泛應(yīng)用于植被建設(shè),其中魚鱗坑、反坡臺(tái)兩類工程措施的應(yīng)用已有四十逾年的歷史,采用水平溝、水平階兩類工程措施進(jìn)行植被恢復(fù)的時(shí)間也有二十年以上,并采用人工植被檸條、側(cè)柏、油松、山杏與不同工程措施相結(jié)合進(jìn)行植被恢復(fù),形成了人工植被與工程措施相結(jié)合的建設(shè)體系。

2 研究方法

2.1 試驗(yàn)地基本特征

于2015年生長季(5月—10月),以流域內(nèi)開展植被恢復(fù)時(shí)間最長的油松魚鱗坑和油松反坡臺(tái)為研究對(duì)象(1972年始建),同時(shí)以油松自然坡面樣地為對(duì)照,采用微創(chuàng)式零干擾小區(qū)搭建技術(shù)布設(shè)油松魚鱗坑和油松反坡臺(tái)監(jiān)測小區(qū)(圖1),搭建完畢后進(jìn)行土壤樣品的采集,并采用每木檢尺對(duì)植被進(jìn)行調(diào)查,記錄不同植被的樹高、胸徑和冠幅,試驗(yàn)地基本特征見表1。

表1 試驗(yàn)地基本特征

圖1 魚鱗坑、反坡臺(tái)及對(duì)照樣地示意圖Fig.1 Diagram of fish scale pits, adversely graded tableland and control site

2.2 土壤水分特征曲線測定

野外樣品采樣通過剖面法進(jìn)行,分別在油松魚鱗坑和油松反坡臺(tái)坡面小區(qū)挖2 m長、1.5 m寬、1 m深的剖面,按0—20、20—40、40—60和60—80 cm 分層取樣,每層重復(fù)3次。采用沙箱法和壓力膜法測定土壤水分特征曲線,其中沙箱測定范圍為0—80 cm水柱,將裝有土樣的環(huán)刀放置砂箱中,測完基質(zhì)勢為-80 cm時(shí)將環(huán)刀取出。然后將環(huán)刀土置于壓力板上,進(jìn)行24 h飽和后放入壓力鍋中由低壓到高壓依次進(jìn)行脫水測試,每次加壓后至水不再流出時(shí)即達(dá)到平衡。平衡后取出土樣,迅速稱重,再將土樣重新裝入壓力室內(nèi)。重復(fù)上述過程,直至最后一個(gè)壓力,取出土樣后用烘干法測其土壤含水量,通過計(jì)算得到土壤體積含水量與土壤基質(zhì)勢的關(guān)系[15]。

2.3 水分特征曲線模型選取

Brook-Corey模型、Gardner模型和van Genuchten模型是國內(nèi)外學(xué)者較常使用且擬合較好的土壤水分特征曲線模型,3種模型均可通過多元非線性回歸的方式求解參數(shù)值。模型的表現(xiàn)形式如下：

Gardner:

h=aθ-b

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

土壤飽和導(dǎo)水率采用RETC軟件通過實(shí)測機(jī)械組成及土壤容重計(jì)算而得。采用SAS9.2對(duì)不同工程措施下土壤水分特征曲線模型進(jìn)行非參數(shù)求解,采用混合線性模型分析不同工程措施土壤機(jī)械組成與飽和導(dǎo)水率。相關(guān)數(shù)據(jù)采用平均值±標(biāo)準(zhǔn)差(Mean±SD)的方式來表示,同一深度不同工程措施之間顯著性采用最小顯著差數(shù)法(LSD)進(jìn)行檢驗(yàn)(顯著性水平P< 0.05),繪圖采用Sigmaplot12.5完成。

3 結(jié)果與分析

3.1 土壤機(jī)械組成

油松魚鱗坑、油松反坡臺(tái)及對(duì)照樣地土壤機(jī)械組成變化見表2,由表可知,不同工程措施土壤機(jī)械組成之間差異顯著(P<0.05),且不同深度表現(xiàn)不同。具體來說,在0—20 cm,油松反坡臺(tái)土壤粘粒(11.58%)、粉粒含量(73.48%)顯著高于油松魚鱗坑(9.73%、68.45%)以及坡面對(duì)照(9.71%、67.99%),但油松魚鱗坑與對(duì)照樣地之間粘粒含量之間差異不顯著(P>0.05),同時(shí),兩者粉粒含量顯著高于坡面對(duì)照(P<0.05)。20—40 cm,魚鱗坑(71.61%)和反坡臺(tái)(71.73%)粘粒含量顯著高于對(duì)照(68.60%),同時(shí),粉粒含量兩者均顯著高于對(duì)照(P<0.05),但兩者之間無顯著差異(P>0.05)。40—60 cm,魚鱗坑(72.52%)和反坡臺(tái)(70.54%)之間粉粒含量無顯著差異(P>0.05),但兩者顯著高于對(duì)照(67.29%)。粘粒含量反坡臺(tái)(10.26%)顯著高于魚鱗坑(9.95%),同時(shí)兩者均顯著高于對(duì)照(9.76%)。

表2 不同工程措施土壤粒徑分布

不同小寫字母表示統(tǒng)一工程措施不同深度之間差異顯著(P<0.05),不同大寫字母表示同一深度不同工程措施之間差異顯著(P<0.05),數(shù)值表現(xiàn)為平均值±標(biāo)準(zhǔn)差

3.2 飽和導(dǎo)水率

油松魚鱗坑、油松反坡臺(tái)和坡面對(duì)照之間土壤飽和導(dǎo)水率差異顯著(P<0.05),不同深度表現(xiàn)不同(圖2)。具體來說,0—20 cm魚鱗坑和反坡臺(tái)飽和導(dǎo)水率無顯著差異(P>0.05),但兩者均顯著高于坡面對(duì)照(P<0.05),這也與 20—40 cm相同。油松魚鱗坑在 40—80 cm 與對(duì)照之間飽和導(dǎo)水率不顯著(P>0.05),60—80 cm 三者之間飽和導(dǎo)水率無顯著差異(P>0.05)。

圖2 不同工程措施土壤飽和導(dǎo)水率 Fig.2 Soil saturated hydraulic conductivity of different engineering measures 不同小寫字母表示同一深度不同工程措施之間差異顯著,P<0.05

3.3 土壤水分特征曲線模擬

分別g使用三種土壤水分特征曲線模型對(duì)油松不同工程措施土壤水分特征曲線進(jìn)行參數(shù)擬合,得到每個(gè)模型參數(shù)的模擬值以及判別系數(shù)(表3),根據(jù)不同模型的擬合結(jié)果,可知不同工程措施不同深度,各模型擬合結(jié)果并不一致。由表3可知,對(duì)于油松魚鱗坑來說,van Genuchten模型在0—20、20—40 cm和40—60 cm深度的模擬精度較好,判別系數(shù)分別為0.950、0.927和0.898,Gardner模型對(duì)60—80 cm深度的模擬精度高于其他兩組模型,對(duì)油松反坡臺(tái)來說,Gardner模型對(duì)40—60 cm的模擬精度高于其他兩組模型,除40—60 cm外,van Genuchten模型的模擬精度均高于0.900。由此可知,3種模型均能夠很好的擬合不同工程措施土壤水分特征曲線。此外,通過對(duì)不同工程措施實(shí)測值與模型模擬的擬合結(jié)果進(jìn)行比較(圖3、圖4和圖5),具體反映在曲線上可表現(xiàn)為低吸力條件下曲線比較陡直,在中高吸力段區(qū)間內(nèi)的曲線漸趨平緩。與坡面對(duì)照相比,魚鱗坑和反坡臺(tái)表現(xiàn)出較好的保水功能。

表3 土壤水分特征曲線參數(shù)擬合

圖3 油松反坡臺(tái)不同深度土壤水分特征曲線模擬Fig.3 Soil water retention curve of different depth of adversely graded tableland- Pinus tabulaeformis

圖4 油松魚鱗坑不同深度土壤水分特征曲線模擬Fig.4 Soil water retention curve of different depth of fish-scale pits- Pinus tabulaeformis

圖5 對(duì)照樣地不同深度土壤水分特征曲線模擬Fig.5 Soil water retention curve of different depth of control plot

3.4 土壤有效水

采用田間持水量、有效水含量等指標(biāo)來表征不同工程措施措施的土壤持水能力,反坡臺(tái)、魚鱗坑與自然坡面水分有效性常數(shù)見表4。由表可知,魚鱗坑、反坡臺(tái)的土壤水分常數(shù)較低,0—20 cm,魚鱗坑(0.136 cm3/cm3)、反坡臺(tái)(0.135 cm3/cm3)土壤有效水含量均高于對(duì)照(0.106 cm3/cm3)。結(jié)合不同深度(0—80 cm)兩類工程措施土壤有效水含量,較之對(duì)照(0.118 cm3/cm3),魚鱗坑(0.138 cm3/cm3)和反坡臺(tái)(0.129 cm3/cm3)平均有效水分分別提高15%和9%。

表4 不同工程措施土壤有效水含量

4 討論

人們通常使用模型根據(jù)土壤機(jī)械組成、土壤容重等性質(zhì)分層填充的重塑土進(jìn)行測定和推導(dǎo),反映去除影響土壤性質(zhì)因素后的土壤水力學(xué)特征[16]。土壤質(zhì)地、容重以及土壤結(jié)構(gòu)是影響水分特征曲線的主要因素,其中土壤質(zhì)地對(duì)水分特征曲線的影響最為明顯,魚鱗坑、反坡臺(tái)等工程措施重塑下墊面微地形,改變了徑流路徑,影響了土壤顆粒的分離、搬運(yùn)和沉積過程,工程措施開展后,在與植被長期作用下對(duì)土壤表面積及土壤孔隙產(chǎn)生影響,改變了土壤質(zhì)地,從而影響水分特征曲線[17-18]。由于土壤含水量主要是指土壤中水分含量的總體值,而真正能夠被植物所吸收利用的水量以及其利用水分的難以程度取決于有效水分含量,植物對(duì)土壤水分利用的有效程度取決于土水勢的高低,而不能單從含水量角度進(jìn)行分析。所以在評(píng)價(jià)工程措施涵養(yǎng)水源的生態(tài)效益時(shí),不僅僅需要關(guān)注土壤水分的變化,更需關(guān)注工程措施開展后土壤有效水分的變化。

圖6 油松魚鱗坑、油松反坡臺(tái)與對(duì)照樣地土壤機(jī)械組成Fig.6 Soil texture distribution of fish-scale pits- Pinus tabulaeformis, adversely graded tableland- Pinus tabulaeformis and cotrol site

圖7 油松魚鱗坑與油松反坡臺(tái)工程措施示意圖 Fig.7 Schematic diagram of the of fish scale pits- Pinus tabulaeformis and adversely graded tableland- Pinus tabulaeformis

對(duì)魚鱗坑、反坡臺(tái)土壤水分特征曲線測定結(jié)果表明,水分特征曲線在低吸力段比較陡直,而在中吸力段和高吸力段區(qū)間內(nèi)逐漸平緩。這主要是由于在低吸力階段的土壤水勢范圍,土壤釋放水分過程取決于土壤結(jié)構(gòu)較粗的孔隙分布。因水吸力較小,脫水過程主要靠毛管力起作用。在高吸力段區(qū)間時(shí),水吸力較大,主要是土壤顆粒表面吸附力起作用。通過對(duì)不同工程措施土壤機(jī)械組成的分析可知(表2),研究區(qū)內(nèi)機(jī)械組成主要分布在粉砂粒徑中(2—2000 μm),較之粉粒,砂粒所占比重也較小,粘粒含量最低,僅在10%上下,所以坡面分布的主要為粉砂質(zhì)壤土,而非黃土母質(zhì)發(fā)育較為普遍的砂質(zhì)壤土。較之對(duì)照,工程措施改變了土壤質(zhì)地,魚鱗坑與反坡臺(tái)長期作用提高了土壤粉粒含量(圖6)。黃土土質(zhì)通常顆粒較粗且發(fā)育緩慢,很難形成好的土壤結(jié)構(gòu)和顆粒組成,導(dǎo)致其持水能力和保水能力較差,工程措施能夠加速土壤顆粒細(xì)化,改善土壤質(zhì)地,從而改善土壤持水能力,在魚鱗坑、反坡臺(tái)等工程措施長期作用下,土壤有效水分得到了提高。結(jié)合對(duì)土壤飽和導(dǎo)水率的計(jì)算結(jié)果,0—20 cm油松魚鱗坑和油松反坡臺(tái)飽和導(dǎo)水率與自然坡面之間差異顯著(P<0.05),但兩者間無顯著差異(P>0.05),由此可見,工程措施對(duì)土壤環(huán)境的改良主要體現(xiàn)在對(duì)表層的影響。通過計(jì)算水分有效性常數(shù),結(jié)果顯示工程措施能夠提高土壤有效水分,魚鱗坑和反坡臺(tái)有效水分較之對(duì)照分別提高15%和9%,對(duì)表層土壤水分有效性的改善更為顯著。魚鱗坑、反坡臺(tái)均能夠減緩坡度,改變上方來水路徑,削減徑流沖刷力并促使降水就地入滲,雖然兩者均能提高土壤有效水分,但因兩種工程措施自身結(jié)構(gòu)差異,致使研究區(qū)內(nèi)魚鱗坑和反坡臺(tái)對(duì)有效水分提高程度不同(圖7)。與此同時(shí),對(duì)魚鱗坑和反坡臺(tái)土壤有效水分測定結(jié)果顯示,兩種工程措施的土壤有效水分較低,主要由于研究區(qū)土壤孔隙較大,粒間孔隙粗,毛管力微弱,施加較小的吸力,大孔隙中的水分就被排出;而保持在小孔隙中的水分只有在較大吸力下才能被排出,這也是研究區(qū)土壤持水量較低的內(nèi)在原因。

利用水分特征曲線既能進(jìn)行土壤水吸力與土壤含水率之間的換算,也能夠間接反映土壤孔隙大小的分布,有利于在應(yīng)用物理、數(shù)學(xué)方法對(duì)水分運(yùn)動(dòng)進(jìn)行定量分析時(shí)提供有關(guān)參數(shù)。選擇Brooks-Corey、Gardner及van Genuchten模型對(duì)土壤水分特征曲線的實(shí)測值和模擬值進(jìn)行擬合的結(jié)果顯示,不同工程措施的水分特征曲線在不同的深度形狀和趨勢基本一致,van Genuchten對(duì)0—20 cm 擬和結(jié)果的決定系數(shù)均在0.95以上,高于其他兩種模型。Brooks-Corey模型形式較Gardner模型復(fù)雜,參數(shù)較多,需考慮飽和含水率、殘余含水率和進(jìn)氣吸力對(duì)水分特征曲線的影響,Gardner模型參數(shù)較少,形式簡單,運(yùn)算和求解較為方便。三種模型也被廣泛的應(yīng)用于土壤水分特征曲線的擬合研究中[19-22]。van Genuchten模型是最常用的描述土壤水分特征曲線的模型,該模型適用的土壤質(zhì)地范圍和土壤含水量范圍均較為廣泛,是推求土壤水分運(yùn)動(dòng)參數(shù)的常用模型。總體來看,van Genuchten模型擬合精度變化較小,但參數(shù)較多。相關(guān)學(xué)者基于不同水吸力下的土壤含水量和RETC主要涉及的4種模型,分別確定了van Genuchten為最優(yōu)的土壤水分特征曲線擬合模型[23- 24]。從運(yùn)算角度看,三種模型均可適用于長期開展水土保持工程措施后土壤水分特征曲線的模擬。采用非線性擬合求解參數(shù)的結(jié)果表明,van Genuchten模型的決定系數(shù)較高。

水分特征曲線是研究水分入滲、蒸發(fā)等生態(tài)水文過程的關(guān)鍵,對(duì)于探究土壤水分的貯存、保持、供應(yīng)以及土壤-植物-大氣連續(xù)體中的水分動(dòng)態(tài)具有重要意義。土壤是工程措施與植被影響的直接環(huán)境,兩者共同影響土壤屬性[25],水分是開展植被恢復(fù)的主要因子,所以通過實(shí)測土壤水分特征曲線并定量分析土壤有效水分的變化,能夠?yàn)樵u(píng)價(jià)不同工程措施保持水土的生態(tài)服務(wù)功能提供依據(jù)。當(dāng)前,工程措施被廣泛應(yīng)用于植被建設(shè)中,通過重塑下墊面微地形,并在與植被長期作用下改變了土壤質(zhì)地(圖6),土壤質(zhì)地是影響水分特征曲線的關(guān)鍵因素之一。較之對(duì)照,魚鱗坑和反坡臺(tái)顯著提高了土壤粘粒和粉粒含量,兩者土壤有效水分分別提高15%和9%,研究區(qū)內(nèi)工程措施的開展對(duì)生態(tài)恢復(fù)具有積極作用。

5 結(jié)論

1)黃土半干旱小流域,魚鱗坑和反坡臺(tái)兩種工程措施長期作用改變了土壤機(jī)械組成,反坡臺(tái)(10.63%)、魚鱗坑(9.78%)粘粒含量顯著高于對(duì)照樣地(9.66%),粉粒含量反坡臺(tái)(71.42%)與魚鱗坑(70.74%)顯著高于對(duì)照樣地(67.85%)。

2)工程措施有效提高土壤飽和導(dǎo)水率,與對(duì)照樣地相比,魚鱗坑和反坡臺(tái)表現(xiàn)出較好的保水功能。魚鱗坑、反坡臺(tái)均能夠提高土壤有效水分,魚鱗坑和反坡臺(tái)有效水分較之對(duì)照分別提高15%和9%,且兩種工程措施對(duì)表層土壤有效水分改善更加顯著。

3)采用Brooks-Corey、Gardner及van Genuchten模型對(duì)不同深度土壤水分特征曲線實(shí)測值和模擬值進(jìn)行擬合,兩種工程措施的水分特征曲線不同深度形狀和趨勢基本一致,非線性擬合求解參數(shù)的結(jié)果表明van Genuchten模型決定系數(shù)較高。