蘭道云,畢華興,2,3,4,5,趙丹陽,王 寧,云慧雅,王珊珊,崔艷紅
(1.北京林業(yè)大學(xué)水土保持學(xué)院,北京 100083;2.山西吉縣森林生態(tài)系統(tǒng)國家野外科學(xué)觀測研究站,北京 100083;3.水土保持國家林業(yè)局重點(diǎn)實(shí)驗(yàn)室(北京林業(yè)大學(xué)),北京 100083;4.北京市水土保持工程技術(shù)研究中心(北京林業(yè)大學(xué)),北京 100083;5.林業(yè)生態(tài)工程教育部工程研究中心(北京林業(yè)大學(xué)),北京 100083)
晉西黃土區(qū)嚴(yán)重的土壤侵蝕導(dǎo)致大量的土壤和養(yǎng)分流失,使該地區(qū)的自然土壤遭到破壞。植被恢復(fù)是改善土壤環(huán)境最有效的方法。國內(nèi)外研究表明,植被的生長發(fā)育與土壤性狀的演變之間是互饋?zhàn)饔?,土壤為植物生長提供水分和養(yǎng)分,其含量不僅影響植被個體發(fā)育,也進(jìn)一步?jīng)Q定植物群落的類型和分布,同時植被也是土壤有機(jī)質(zhì)最主要的來源,對土壤理化性質(zhì)有著深刻的影響。土壤是生態(tài)系統(tǒng)諸多生態(tài)過程的載體,植物的生長對土壤結(jié)構(gòu)與養(yǎng)分狀況起著關(guān)鍵作用。
油松()是黃土區(qū)主要造林樹種之一,具有保育土壤、涵養(yǎng)水源、固碳釋氧等功能。保育土壤功能是指森林中活地被物和凋落物截留降水、降低水滴對表土的沖擊和地表徑流的侵蝕作用;同時林木根系固持土壤,從而減少土壤肥力損失,以及改善土壤結(jié)構(gòu)的功能,是森林生態(tài)系統(tǒng)中重要的服務(wù)功能之一。油松林地的固土保肥能力反映了油松林地抵抗水土流失、改善土壤肥力的能力。近年來諸多學(xué)者對林地的保育土壤功能進(jìn)行了研究,田寧寧等使用固土量公式和保肥量公式對晉西黃土區(qū)退耕還林地的保育土壤功能進(jìn)行評價;張向峰等使用模糊綜合評價法對縉云山水源涵養(yǎng)林的保育土壤功能進(jìn)行了評價。目前對林地保育土壤功能客觀評價方法的研究還較少,熵權(quán)法是一種客觀計算各指標(biāo)權(quán)重的數(shù)學(xué)方法,能夠剔除主觀因素的影響,客觀地判斷評價指標(biāo)的相對重要性。在面對土壤保育功能指標(biāo)本身涉及到的數(shù)據(jù)多且復(fù)雜,熵權(quán)法解決了因各層次內(nèi)涵不同、不能相加的弊端。近年來,熵權(quán)法已成功運(yùn)用于土壤抗蝕性評價和水源涵養(yǎng)評價,說明熵權(quán)法在保育土壤方面也同樣具有可行性與科學(xué)性。
關(guān)于不同林分密度對林分土壤及其影響因素早有報道,主要集中在不同地區(qū)、不同植被類型對其土壤養(yǎng)分、物理性質(zhì)或者土壤抗蝕性單獨(dú)進(jìn)行研究,缺少有關(guān)林分密度對保育土壤功能的影響及綜合評價方面的研究。本文以晉西黃土區(qū)4種不同密度的油松人工林為研究對象,對研究區(qū)的土壤理化性質(zhì)、坡面產(chǎn)流產(chǎn)沙進(jìn)行分析,運(yùn)用熵權(quán)法客觀評價油松人工林的保育土壤功能,擬從林分密度對晉西黃土區(qū)油松人工林保育土壤功能角度出發(fā),選擇林分經(jīng)營密度。
研究區(qū)位于山西省臨汾市吉縣蔡家川流域(36°14′27′′—36°18′23′′N,110°39′45′′—110°47′45′′E),流域面積40 km,海拔900~1 590 m,屬黃土高原東南部暖溫帶半濕潤區(qū),年平均氣溫10 ℃,多年平均降水量574.8 mm,年平均蒸發(fā)量1 726.3 mm。屬于典型的黃土殘塬溝壑區(qū),土壤類型為微堿性褐土,黃土母質(zhì)。該地區(qū)主要喬木樹種有油松(Carr.)、刺槐(L.)、側(cè)柏((L.)Franco)、山楊(Dode)、遼東櫟(Mary)等;主要灌木樹種有沙棘(Linn)、虎榛子(Decne.)、黃刺玫(Lindl)、胡枝子(Turcz.)等;主要的草本植物有艾蒿(H.Lév.&Vaniot)、草地早熟禾(L.)、黃花蒿(Linn.)等。
1.2.1 樣地選擇 研究區(qū)于1991年開始大面積栽植油松?;诖罅繉?shí)際調(diào)查,在研究區(qū)選擇林齡為29年、成土母質(zhì)基本相同、坡度、坡位、坡向等立地條件基本一致的油松林為研究對象,設(shè)置低密度(925株/hm)、中密度(1 325,1 750株/hm)和高密度(2 250株/hm)3個密度梯度,每個密度梯度3個重復(fù)。布設(shè)20 m×20 m標(biāo)準(zhǔn)樣地,并對樣地內(nèi)的喬木每木檢尺,記錄樣地內(nèi)所有樹木的胸徑、樹高等生長指標(biāo);根據(jù)“五點(diǎn)法”布設(shè)5 m×5 m灌木樣方5個,1 m×1 m草本樣方5個,記錄林下植被的生長情況。使用地質(zhì)羅盤測定坡向和坡度;使用GPS數(shù)據(jù)采集儀定位測定海拔高度。不同密度油松人工林地基本特征見表1。
表1 不同密度油松人工林樣地基本特征
1.2.2 土壤取樣與測定方法 在不同密度油松林地(各3塊樣地)中心點(diǎn)處選擇1個土壤剖面采樣點(diǎn),沿著土壤剖面0—10,10—20,20—30,30—40,40—50,50—60 cm分層取樣,每層3個重復(fù),帶回實(shí)驗(yàn)室,待其風(fēng)干后研磨過篩裝袋,用于測定土壤化學(xué)指標(biāo)。土壤容重、孔隙度等采用環(huán)刀法測定;土壤化學(xué)性質(zhì)的測定指標(biāo)主要有土壤有機(jī)質(zhì)、全氮、全磷、銨態(tài)氮、硝態(tài)氮和速效磷等。土壤有機(jī)質(zhì)采用重鉻酸鉀—外加熱法測定,土壤全氮采用半微量凱式法測定,土壤全磷采用硫酸—混合加速劑消煮—蒸餾法測定,土壤銨態(tài)氮、硝態(tài)氮采用堿解擴(kuò)散法和氫氧化鈉熔融—鉬銻抗比色法測定,土壤速效磷采用NHF—HCL浸提—鉬銻抗比色法測定。
在油松林地坡度相近的坡面布設(shè)8個標(biāo)準(zhǔn)徑流小區(qū)(每種密度各2個),對4種密度油松林地進(jìn)行天然降雨條件下產(chǎn)流產(chǎn)沙的觀測,并記錄場降雨條件下的降雨量(林內(nèi)降雨)、徑流量和泥沙含量。
1.2.3 數(shù)據(jù)處理與分析 試驗(yàn)取得的數(shù)據(jù)采用Excel 2019進(jìn)行處理,各指標(biāo)間數(shù)據(jù)分析采用SPSS 26.0軟件分析,其中以林分密度為自變量,以不同層次容重、總孔隙度、毛管孔隙度、非毛管孔隙度、有機(jī)質(zhì)、全氮、全磷、銨態(tài)氮、硝態(tài)氮和速效磷作為因變量進(jìn)行單因素方差分析(One-way Anova),采用Duncan多重比較進(jìn)行差異性分析。
1.2.4 保育土壤功能評價方法 依據(jù)保育土壤功能的定義,結(jié)合前人的研究成果,按照科學(xué)性、可操作性的原則,選取能夠反映油松人工林保育土壤功能(土壤結(jié)構(gòu)改良、土壤肥力改善、截流減沙)的12指標(biāo),包括容重、總孔隙度、毛管孔隙度、非毛管孔隙度、有機(jī)質(zhì)、全氮、全磷、銨態(tài)氮、硝態(tài)氮、速效磷、徑流量和泥沙量,建立油松人工林保育土壤功能評價體系。采用熵權(quán)法(the entropy weight method, EWM)評價林分保育土壤能力,該方法屬于客觀賦權(quán)法,通過調(diào)查數(shù)據(jù)結(jié)果提供的信息大小來確定權(quán)重。一般對于某項(xiàng)指標(biāo)可以用熵值來判斷其離散程度,其信息熵越小,離散程度越大,說明該指標(biāo)的權(quán)重越大。利用信息熵計算出各指標(biāo)的權(quán)重,為多指標(biāo)綜合評價提供依據(jù)。
(1)數(shù)據(jù)標(biāo)準(zhǔn)化處理。由于選取的評價指標(biāo)單位不同,需要進(jìn)行無量綱化處理,具體分為2種情況:
式中:為第個評價對象對應(yīng)第個指標(biāo);為標(biāo)準(zhǔn)化后對應(yīng)的指標(biāo)值。
(4)信息熵冗余度計算公式為:=1-。
式中:為指標(biāo)個數(shù)。
土壤容重和孔隙度是反映土壤物理結(jié)構(gòu)的重要參數(shù),容重是土壤緊實(shí)度指標(biāo),反映土壤透水通氣性的大小。由表2可知,在垂直剖面上,隨著土層深度增加,不同密度油松林地土壤容重均呈增大趨勢,深層可達(dá)1.30 g/cm。在0—60 cm土層,油松林地土壤平均容重為1.21 g/cm,變化范圍為0.94~1.35 g/cm。不同密度油松林地土壤平均容重從小到大依次為1 750株/hm<1 325株/hm<2 250株/hm<925株/hm。差異性檢驗(yàn)結(jié)果表明,中密度(1 325,1 750株/hm)油松林地土壤容重與低密度(925株/hm)和高密度(2 250株/hm)油松林地間均存在顯著差異(<0.05)。隨著密度的增大,土壤容重的總體趨勢表現(xiàn)為先減小后增大。
表2 不同密度油松人工林地土壤容重
土壤孔隙可以為植物提供水和空氣,包括毛管孔隙和非毛管孔隙2種,體現(xiàn)了土壤貯蓄和調(diào)節(jié)水分能力的大小。在垂直剖面上,不同密度油松林地土壤總孔隙度、毛管孔隙度和非毛管孔隙孔隙度均隨土層深度的增加而減小(表3)。0—60 cm土層土壤平均總孔隙度為52.04%,變化范圍為49.41%~55.69%,不同密度油松林地土壤平均總孔隙度從小到大依次為925株/hm<1 750株/hm<1 325株/hm<2 250株/hm。差異檢驗(yàn)結(jié)果表明,密度為1 750株/hm油松林地土壤平均總孔隙度與其他密度油松林地之間差異顯著(<0.05)。0—60 cm土層土壤平均毛管孔隙度為48.13%,變化范圍為44.78%~51.83%,不同密度油松林地土壤平均毛管孔隙度從小到大依次為925/hm<1 750株/hm<1 325株/hm<2 250株/hm,差異檢驗(yàn)結(jié)果表明,密度為1 750株/hm的油松林地土壤平均毛管孔隙度和925株/hm油松林地之間無顯著差異,而與1 325,2 250株/hm油松林地之間差異顯著(<0.05)。在0—60 cm土層土壤平均非毛管孔隙度為3.91%,變化范圍為3.01%~4.64%,不同密度油松林地平均非毛管孔隙度從小到大依次為925株/hm<1 325株/hm<2 250株/hm<1 750株/hm,差異檢驗(yàn)結(jié)果表明,密度為1 750株/hm油松林地土壤平均非毛管孔隙度與其他林分密度之間差異顯著(<0.05)。
表3 不同密度油松人工林地土壤孔隙狀況
表4為不同密度油松林地全量養(yǎng)分狀況。土壤全量養(yǎng)分主要包括全氮和全磷,其含量主要受到森林生態(tài)系統(tǒng)植被養(yǎng)分循環(huán)過程和土壤成土母質(zhì)的影響。土壤有機(jī)質(zhì)是土壤養(yǎng)分的主要來源,影響土壤的潛在生產(chǎn)力,是土壤肥力狀況的重要指標(biāo),一般指存在于土壤中的所有含碳有機(jī)物。0—60 cm土層不同密度油松林地土壤平均有機(jī)質(zhì)含量為7.29 g/kg,變化范圍為3.04~15.94 g/kg,且隨著土層深度的增加而減少。不同密度油松林地土壤平均有機(jī)質(zhì)含量從大到小依次是1 325株/hm>1 750株/hm>925株/hm>2 250株/hm。差異性檢驗(yàn)結(jié)果表明,密度為1 750株/hm的油松林地土壤平均有機(jī)質(zhì)含量分別與925,2 250株/hm油松林地差異顯著(<0.05)。
表4 不同密度油松人工林地土壤全量養(yǎng)分狀況
0—60 cm土層土壤平均全氮含量為0.32 g/kg,變化范圍為0.16~0.55 g/kg。不同密度油松林土壤平均全氮含量從大到小依次為1 750株/hm>1 325株/hm>925株/hm>2 250株/hm。差異檢驗(yàn)結(jié)果表明,林分密度為1 750株/hm油松林地土壤平均全氮含量與其他林分密度之間差異顯著(<0.05)。0—60 cm土層土壤平均全磷含量為0.51 g/kg,變化范圍0.48~0.54 g/kg,不同密度油松林土壤平均全磷含量從大到小依次為1 750株/hm>2 250株/hm>1 325株/hm>925株/hm。差異檢驗(yàn)結(jié)果表明,密度為1 750株/hm油松林地土壤全磷含量與其他林分密度之間差異顯著(<0.05)。土壤全量養(yǎng)分含量最大的為密度1 750株/hm油松林地,且隨著密度增大呈先增大后減少的趨勢。
土壤速效類養(yǎng)分(銨態(tài)氮、硝態(tài)氮、速效磷)含量高低直接關(guān)系土壤肥力的高低。0—60 cm土層土壤銨態(tài)氮平均含量為10.63 mg/kg,變化范圍為7.76~14.71 mg/kg(表5)。不同密度油松林地土壤銨態(tài)氮平均含量從大到小依次為2 250株/hm>1 325株/hm>1 750株/hm>925株/hm。差異性檢驗(yàn)結(jié)果表明,4種不同密度油松林地之間土壤銨態(tài)氮平均含量差異顯著(<0.05)。0—60 cm土層土壤硝態(tài)氮平均含量為13.33 mg/kg,變化范圍為7.76~18.03 mg/kg,不同密度油松林地土壤硝態(tài)氮平均含量從大到小依次為1 750株/hm>1 325株/hm>925株/hm>2 250株/hm,差異檢驗(yàn)結(jié)果表明,密度為1 750株/hm油松林地土壤硝態(tài)氮平均含量與其他林分密度之間差異顯著(<0.05)。0—60 cm土層土壤速效磷平均含量為2.89 mg/kg,變化范圍為2.19~3.42 mg/kg,各密度林分速效磷平均含量從大到小依次為1 750株/hm>925株/hm>1 325株/hm>2 250株/hm,密度為1 750株/hm的油松林地分別與1 325,2 250株/hm油松林地之間土壤平均速效磷差異顯著(<0.05)。密度為1 750株/hm油松林地土壤速效類養(yǎng)分含量最大。
表5 不同密度油松人工林地速效養(yǎng)分狀況
在林地上由于活地被物對降雨截流和枯落物的攔截,降雨并不能直接到達(dá)林地表面形成徑流,只有滿足活地被物截流量和枯落物的吸水量對降雨的再分配后,降雨才能到達(dá)土壤表面,對土壤表面產(chǎn)生沖擊,當(dāng)?shù)竭_(dá)地面的凈雨強(qiáng)度大于土壤的滲透強(qiáng)度,地表產(chǎn)生徑流,對土壤產(chǎn)生侵蝕。從表6可以看出,油松林地徑流量總體表現(xiàn)為隨著林分密度增大呈先減小后增大趨勢。4種不同密度油松人工林徑流量為7.90~9.01 L/100 m,不同密度油松林徑流量從小到大依次為1 750株/hm<2 250株/hm<1 325株/hm<925株/hm,其中1 750株/hm油松林地的地表徑流量最小。泥沙量表現(xiàn)為隨著林分密度增大呈先減小后增大趨勢,4種不同密度油松林泥沙量為38.49~121.88 g/100 m,不同密度油松林泥沙量從小到大依次為1 750株/hm<2 250株/hm<1 325株/hm<925株/hm,1 750株/hm密度的油松林地泥沙量最小。
表6 不同密度油松人工林降雨徑流泥沙含量
森林生態(tài)系統(tǒng)的保育土壤功能主要體現(xiàn)在活地被物對土壤結(jié)構(gòu)改良、土壤肥力改善及截流減沙的能力,表征林分保育土壤能力的評價指標(biāo)較多,指標(biāo)間相互依賴或排斥,且林分保育土壤能力受到各指標(biāo)的影響大小不同,因此綜合采用多個指標(biāo)評價林分保育土壤能力時需要對指標(biāo)權(quán)重進(jìn)行賦值,本研究運(yùn)用熵權(quán)法對不同密度油松林地指標(biāo)權(quán)重賦值(表7)。
由表7可知,各項(xiàng)指標(biāo)的權(quán)重為0.059 9~0.110 1。通過熵權(quán)法計算得出的結(jié)果,結(jié)合權(quán)重與各指標(biāo)標(biāo)準(zhǔn)化后數(shù)值,計算各密度林分保育土壤能力的綜合評價結(jié)果為1 750株/hm(0.77)>1 325株/hm(0.52)>2 250株/hm(0.45)>925株/hm(0.25)。隨著林分密度增大,油松林保育土壤的功能先增強(qiáng)后減弱,密度為1 750株/hm的林地土壤結(jié)構(gòu)、土壤肥力和林地截流減沙能力的效果較好。
表7 不同密度油松人工林各指標(biāo)權(quán)重
通過對晉西黃土區(qū)4種不同密度油松林地土壤物理特性進(jìn)行研究發(fā)現(xiàn),隨著土層深度增加,土壤容重逐漸增大;隨著林分密度增大,土壤容重先減小后增大,1 750株/hm密度下油松林地土壤容重最小。陳莉莉等對秦嶺不同密度油松土壤特性的研究表明,林分密度對土壤容重具有顯著影響,且隨林分密度增大土壤容重先減小后增大,與本研究結(jié)果相同。一般情況下,土壤中大小孔隙同時存在,總孔隙度在50%左右,其中非毛管孔隙占1/5~2/5為最好,此種情況下土壤孔隙通透性比較協(xié)調(diào)。本研究中4種林分密度林地土壤總孔隙度具有顯著的差異性,且都在50%左右,但非毛管孔隙度偏小(3.41%~4.40%),原因在于土壤孔隙狀況除了受林分密度影響外,還可能與林下植被、土壤動物有關(guān)。過小或過大的林分密度均使土壤物理性質(zhì)降低,這可能是由于土壤物理性質(zhì)受土壤質(zhì)地、結(jié)構(gòu)和有機(jī)質(zhì)等眾多因素的影響,相較于低密度林分而言,中密度的林分植被根系對土壤結(jié)構(gòu)具有更好的改善作用,林下水熱條件適中,枯落物層豐厚,有利于水分和養(yǎng)分循環(huán)。但植被對土壤物理性質(zhì)的改善作用有限,林分密度過大時,林木間競爭增強(qiáng),反而不利于土壤發(fā)育。
不同密度油松林地間土壤養(yǎng)分含量均存在差異,且隨著林分密度增大,土壤有機(jī)質(zhì)和土壤氮、磷元素均表現(xiàn)為先增多后減少的趨勢。其原因可能是林分密度通過改變光照和土壤微環(huán)境影響土壤有機(jī)質(zhì)的輸入輸出,同時有機(jī)質(zhì)的分解產(chǎn)生氮元素和磷元素,以供植物生長。隨著林分密度增大,枯落物增多,凋落物的腐殖化過程能分解更多的有機(jī)質(zhì)、氮元素和磷元素,土壤養(yǎng)分含量也隨之增加。但林分密度過大時,針葉凋落物分解較慢,土壤養(yǎng)分的消耗速率大于補(bǔ)償速率,導(dǎo)致土壤養(yǎng)分含量減小,這與任麗娜等的研究結(jié)果基本一致。土壤養(yǎng)分除全磷外,均存在表聚現(xiàn)象,即富集在0—30 cm土層,全磷含量隨土層深度增加無明顯變化,與王巖松等的研究結(jié)果一致。這可能是由于磷元素主要來自于巖石的風(fēng)化淋溶,受到土壤母質(zhì)的影響,在土壤中分布比較均勻;同時,林分密度還對坡面徑流泥沙產(chǎn)生影響,王曉康等的研究表明,徑流泥沙量隨林分密度的增大呈先減小后增大趨勢,與本研究結(jié)果一致。林分密度增大,郁閉度也隨之增大,因?yàn)榱謨?nèi)生境狀況改變,植物對光照、水分、養(yǎng)分條件的利用方式發(fā)生變化,從而導(dǎo)致林下物種組成和結(jié)構(gòu)差異;林分密度過大時,灌草覆蓋度不一定大,強(qiáng)降雨條件下對土壤侵蝕的防治作用較差。
保育土壤功能是森林生態(tài)系統(tǒng)重要的生態(tài)服務(wù)功能之一,盧妮妮等對杉木及其混交林保育土壤功能研究時,以土壤養(yǎng)分的總含量作為林分保育土壤的指標(biāo),對土壤養(yǎng)分元素運(yùn)用主成分分析法進(jìn)行評價;賀山峰等對科爾沁沙地小葉錦雞兒群落保育土壤功能的研究中,通過對比栽植前后理化性質(zhì)對保育土壤功能進(jìn)行判斷,但并未量化指標(biāo)。與前2位學(xué)者的研究不同,本研究在土壤物理結(jié)構(gòu)和土壤養(yǎng)分的基礎(chǔ)上加入徑流泥沙量這一新的監(jiān)測指標(biāo),采用熵權(quán)法對土壤保育功能進(jìn)行量化評價,為探究保育土壤功能的指標(biāo)選取及量化評價提供了嘗試。本研究結(jié)果顯示,隨著林分密度的增大,保育土壤功能先增強(qiáng)后減弱,與已有研究結(jié)果一致。評價結(jié)果表明,中等種植密度(1 750株/hm)比較適合作為研究的經(jīng)營密度。也有研究分別在土壤物理性質(zhì)、養(yǎng)分循環(huán)、徑流泥沙方面對不同密度的油松進(jìn)行探索表明,中等密度(1 500,1 675株/hm)有利于提高土壤物理結(jié)構(gòu)和養(yǎng)分積累,與本研究結(jié)果一致。賀康寧等的研究也表明,1 900株/hm的油松林具有較好的攔蓄徑流、固持土壤能力。綜上所述,密度為1 750株/hm的油松人工林作為經(jīng)營密度能夠改善土壤結(jié)構(gòu),提高土壤肥力,減少徑流泥沙,并且能為晉西黃土區(qū)的水土保持林經(jīng)營管理提供一定參考。
(1)不同密度油松林地之間土壤容重、土壤孔隙度、有機(jī)質(zhì)、土壤全量養(yǎng)分和速效養(yǎng)分含量均存在顯著差異(<0.05),且林分密度變化不同程度地改變土壤理化性質(zhì)和徑流泥沙等12項(xiàng)指標(biāo)。隨著林分密度增大,土壤結(jié)構(gòu)改良狀況和土壤肥力整體呈先上升后下降的趨勢,徑流泥沙量呈先減小后增大趨勢。
(2)通過熵權(quán)法分析評價認(rèn)為,密度為1 750株/hm油松林地得分最高。中等種植密度是研究區(qū)比較合適的營林密度,有利于改善土壤結(jié)構(gòu),增強(qiáng)土壤肥力,提高林地截流減沙能力,使林地具有更好的保育土壤功能。針對研究區(qū)現(xiàn)狀,建議在油松人工林管護(hù)過程中,對林分密度進(jìn)行調(diào)控,以達(dá)到加強(qiáng)林分改良土壤、截流減沙的目的。
(3)保育土壤功能是多個因素的綜合反映,其變化是一個長期的過程,受林分密度、林分生長、外界干擾等多因素影響。今后可進(jìn)一步對林分保育土壤功能進(jìn)行探究,進(jìn)行長時間序列的定位觀測,綜合考慮多方面因素,找出影響林分保育土壤功能的關(guān)鍵因子。