唐 菡,諶 蕓,劉梟宏,何丙輝,李 勇, 強嬌嬌,李 鐵
西南大學資源環境學院,巖溶環境重慶市重點實驗室, 重慶 400715
喀斯特坡地不同于其他地區,具有獨特的“二元三維”[1]地質構造。地表水土流失嚴重,加之成土速率低、土層薄、土體與母巖之間存在明顯軟硬界面的特點[1],一旦當地植被受到人為干擾破壞,水土流失加劇,原本處于平衡狀態的生境將向石漠化方向演變[2]。工程措施(如依自然地形建蓄水池、用擋墻維護落水洞等[3])與植物措施(如關鍵點位布置植物、植物籬等)結合的方式[4]最能有效緩解水土流失加劇的困境。地埂籬為植物籬中的一種,最初用于集水,20世紀初用于水土保持,主要通過根系及根對土體的加強達到固土作用[5]。草本植物須根量大且根系集中分布于淺土層,其對薄土的固持效果較灌、喬木顯著,這是由于深粗根(d>2.0 mm)在土體中起抗滑樁和扶壁的作用,而淺細根(d≤2.0 mm)主要為三維加筋、固結土壤的作用[6- 7]。地埂籬根系主要通過三方面加強土體:1)通過穿插于土體中增大根與土粒、土粒與土粒間的摩擦力;2)根系分泌物促進土粒膠結作用,增加土體粘聚力和抗蝕性能;3)根系自身的抗拉、抗剪性能也能在一定程度上增強土體的抗剪、抗沖力[8- 10]。
拉巴豆(DolichoslablabL.)為豆科蝶形花亞科菜豆族扁豆屬,是熱帶、亞熱帶溫暖地區優良的高營養高產牧草和蔬菜[11],有較高的經濟價值,較強的環境適應性和良好的水土保持功能。我國西南喀斯特地區地勢偏遠、環境惡劣、貧困率高,將拉巴豆應用于喀斯特坡地,既可帶動經濟發展又可保持水土。目前,西南地區雖已推廣種植,但有關其生態效應的研究僅涉及根系生物特性[12],根系的力學特性研究方面存在空白。已有研究證明不同植物根系具有不同的力學特性,即使同一種植物,在不同環境或于不同徑級也存在較大差異[13]。喀斯特植物根系主要受限于土壤厚度和水分,根系形態特征隨著限制條件改變的同時,其化學成分和力學特性也隨之改變[13]。有關根系力學特征的研究表明根系抗拉特性(抗拉強度、延伸率、彈性模量等)與化學成分、微觀結構皆存在一定的相關性[13- 18],但少見以巖溶地貌為背景的根系力學特性方面的研究[6]。本文對喀斯特坡地拉巴豆地埂籬根及根-土復合體在不同坡位的形態分布、力學和化學成分特征進行了綜合分析,從不同坡位、不同徑級根系之間各指標的相關性和差異性著手,探討根及根-土復合體的力學特性及影響因素,旨在探索拉巴豆地埂籬根系的水土保持潛能及在喀斯特地區推廣價值,為喀斯特坡地科學實施水土保持措施提供理論支撐。
研究區位于重慶市酉陽土家族苗族自治縣泔溪鎮龍潭槽谷,地處東經108°58′,北緯28°58′(具體位置見圖1)。槽谷由背斜發育,谷底及順層坡一側為碳酸鹽巖,黃壤,坡度30°左右。該區年均溫14.6℃,年均降雨量1200 mm,年均日照時長1131 h,無霜期261 d,雨熱同期,降雨主要集中于5—9月,屬典型的中亞熱帶季風濕潤氣候。采樣區處于順層坡人工林坡面內,整個坡面植被覆蓋率較低,主要為花椒經濟林,另外存在少量抗旱、耐瘠薄野生喬灌草,如馬尾松(Pinusmassoniana)、紅背山麻桿(Alchorneatrewioides)、沿階草(Ophiopogonbodinieri)等。

圖1 研究區位置圖Fig.1 Study area location map黃色區域表示順層坡研究區
依據自然條件差異,坡面可劃分為上、中、下三個坡位,上坡巖石裸露率及坡度明顯大于中、下坡,而下坡土層最厚,相鄰坡位間海拔相差約為50 m。2018年4月于坡面的上、中、下坡位地埂位置撒播拉巴豆,播種密度0.03 kg/m2(出苗率80%左右),每個坡位籬帶長約200 m、寬約20 cm,常規管護。
全根根系的采集:2018年7月中旬取植株根系,采用挖掘法,每個坡位取5株標準株(標準株的選擇參見[19]),為保證根系完整性,連根帶土一并帶回實驗室。
根-土復合體的采集:每個坡位選取3個對照點(未播拉巴豆地段)和3個地埂籬地點(選點位于拉巴豆標準株附近)。對照裸地取樣時將土壤表面的植株、枯落物等雜質清理干凈,避開雜草根系用鋁盒、容重環刀(底面積20 cm2,高5 cm)、抗剪環刀(底面積30 cm2,高2 cm)取土。地埂籬地取樣時先將拉巴豆地上部分剪去,扒掉2 cm土層(0—2 cm土層無根系生長),每個抗剪環刀均以植株為中心取樣;采集抗沖環刀(長×寬×高:20 cm×10 cm×10 cm)時選取植株密度相近的區域。對照裸地共計36個抗剪土樣,地埂籬地共計36個抗剪土樣和9個抗沖樣。
1.3.1土壤基礎物理性質
土壤容重、孔隙度(總孔隙度、毛管孔隙度和非毛管孔隙度)、飽和含水量、田間持水量和自然含水率的測定采用環刀法[20]。
1.3.2根系形態指標
所取根土樣品置于0.05 mm的網篩中用細弱的水沖洗,分離根土晾干表面水分,整個過程輕柔避免損傷根系。根-土復合體中的根在抗剪、抗沖試驗后再處理,所有根系試樣盡快測定以得到最精確數據。采用EPSON(PERFECTION C700)掃描儀進行根系灰度掃描,WinRHIZO(Pro.2009)根系分析系統進行分析,具體操作步驟和計算過程參考[10]。
1.3.3根系化學成分
將全根根系樣按0.0 1.3.4單根抗拉特性指標 從全根根系中剪取順直、直徑較均一且完好無損的根備用。拉伸試驗夾具鉗口處應力集中,為避免根皮受夾具應力破壞產生脫皮現象或皮斷裂讓試驗機誤記錄為根斷裂,在根兩頭包上醫用橡皮膠,可分散夾具周邊應力增加摩擦力提高試驗成功率(拉伸斷裂點位于根樣中部視為試驗成功),包上膠帶后根系樣本標距為50 mm[24]。采用測力范圍0—5 kN,加載速率范圍0—60 MPa/s,最小分度值1×10-5N的微機控制電子萬能試驗機,夾具型號DSA502A。樣本試驗前用精度為0.01 mm的數顯游標卡尺(型號CD- 6 ASX)測量根系上中下3個點位的直徑取均值即為平均直徑,拉伸速率設置為7 mm/min(速率設置參考[25]),每次試驗微機自動記錄所有數據。 1.3.5根-土復合體力學指標 土壤抗剪強度用南京土壤儀器廠生產的ZJ型應變控制式直剪儀測定并根據庫倫定律計算內摩擦角φ和粘聚力c。沖刷試驗的沖刷槽坡度定為30°,試驗過程及測定方法見[10]。 運用Microsoft Office Excel 2013和 SPSS 22.0軟件進行數據統計分析。進行了差異性檢驗(One-way ANOVA)、協方差分析(Multivariable)和相關分析(Spearman)。 2.1.1根系形態特征 對各坡位全根根系進行形態指標統計以及差異性分析,得出(表1):整體上d≤1.0 mm根系占優勢,全根根系平均重量僅1 mg左右。各全根指標不同坡位之間均無顯著性差異,同一徑級指標在不同坡位,僅下坡根尖數(Root tips,RT)顯著大于上、中坡,同一坡位0.0 全根指標除RT外均于中坡出現最大值。根長(Root length,RL)和根表面積(Root surface area,RSA),下坡大于上坡;RT呈現下坡>中坡>上坡的特征。徑級指標顯示:0.0 結合采樣情況和隨機選取的根系灰度掃描圖2和表1、表2可知:拉巴豆側根發達(尤其是上坡)、大部分向土層深厚的區域延展,主根較短且深度小于15 cm,細根量大;中、下坡細根量較上坡大且RL、RSA、RT和根的分支均較上坡增加;上坡根系無根瘤菌而下坡根瘤菌明顯。上坡土壤物理性質較差,容重大且與中、下坡有顯著差異,總孔隙度、飽和含水量和自然含水率顯著小于中、下坡。由此可見土壤物理性質對根系生長產生了一定的影響,透氣性差且含水量低的上坡土壤的根系各項指標值均較小。 表1 不同坡位拉巴豆根系全根與徑級指標 表中數據為平均值±標準誤差(n=5),不同小寫字母表示同一指標不同徑級間存在顯著差異(P<0.05);不同大寫字母表示同一指標不同坡位間存在顯著差異(P<0.05) 表2 順層坡不同坡位地埂籬土壤的基本物理性質 表中數據為平均值±標準誤差(n=9),不同小寫字母表示同一指標不同坡位間存在顯著差異(P<0.05) 圖2 拉巴豆根系灰度掃描圖Fig.2 Grayscale scan of Dolichos lablab L.A1、A2表示上坡隨機根系; B1、B2表示中坡根系; C1、C2表示下坡根系 2.1.2根系化學成分特征 對根系纖維素(Cellulose,C)、半纖維素(Hemicellulose,H)、木質素(Lignin,L)含量進行測量,結果表明:纖維素含量均值中坡>上坡>下坡,木質素含量排序與纖維素相反,半纖維素含量于中坡值最大,上、下坡均值相近,但各化學成分于各坡位之間均無顯著性差異。 將直徑與這三個指標進行回歸分析,回歸方程見圖3。纖維素含量與直徑的擬合曲線為冪函數關系,半纖維素、木質素含量與直徑均為二次多項式函數關系。相關性分析得出,僅0.0 圖3 不同徑級根系化學成分含量與平均直徑的關系Fig.3 Relationships between root chemical composition content and different diameters 試驗成功的拉巴豆根系直徑范圍為0.19—0.93 mm,平均值為0.46 mm,試驗成功率65.52%。極限抗拉力范圍為1.06—18.56 N,均值6.77 N。抗拉強度范圍6.97—57.50 MPa,均值29.57 MPa。抗拉特性與直徑的擬合方程均為冪函數形式,其指數表征極限抗拉力隨直徑的增長率、衰減率,系數可被視為比例因子[26]。擬合方程中不同坡位的指數和系數均有一定的差異,極限抗拉力與直徑的擬合方程顯示各坡位增長率相似,但下坡比例因子較上、中坡低;抗拉強度與直徑擬合方程顯示中坡較上、下坡衰減率高0.80、1.67倍,比例因子相似。雖同為一種植物,但其根系力學特性由于環境影響造成根系化學成分含量或粗細根比例不一致,從而導致不同坡位根系擬合方程常量的差異。 2.2.1單根極限抗拉力及影響因素 擬合結果顯示(圖4),極限抗拉力與直徑正相關,同一直徑的極限抗拉力中坡>上坡>下坡。將直徑和對應根系極限抗拉力進行自然對數變換,以直徑為協變量進行協方差分析,結果表明直徑對極限抗拉力無顯著影響。極限抗拉力上坡與中、下坡均有顯著性差異。 回歸分析表明(僅對存在顯著相關關系的變量進行說明,表3): 0.5 2.2.2單根抗拉強度及影響因素 抗拉強度與直徑負相關(圖4),協方差分析表明直徑對抗拉強度有顯著性影響(P<0.05)。 0.0 乳化劑是微乳劑配方的關鍵性因素.一般的非離子型表面活性劑HLB值對溫度的變化很敏感,離子型表面活性劑的親水親油平衡對溫度變化不敏感,但其在低溫下的溶解度顯著降低[5].實驗中乳化劑用非離子型和離子型表面活性劑復配,其HLB值在13~15之間. 回歸分析表明(表3):0.0 圖4 抗拉特性與直徑在不同坡位的關系Fig.4 Relationships between maximum tension, tensile strength and mean diameter of Dolichos lablab L. at different slope portions 徑級Diameter根系化學成分Root chemical compositions極限抗拉力Maximum tension F/N抗拉強度Tensile strength T/MPa回歸函數R2Pearson回歸函數R2Pearson0.0 *表示在0.05的水平上顯著相關,**表示在0.01的水平上極顯著相關 2.3.1復合體中根的形態和化學成分特征 表4表明:根系形態指標中僅抗剪土樣根表面積密度(Root surface area density,RSAD)上坡顯著大于中、下坡,根重密度(Root weight density,RWD)上坡顯著大于中坡。抗沖土樣坡位間的根系形態指標均無顯著性差異。根-土復合體抗剪土樣的根系形態指標均為上坡>下坡>中坡而抗沖土樣形態指標于中、下坡較高。具體表現為抗沖土樣根系RSAD于下坡值最大,根體積密度(Root volume density,RVD)、RWD于中坡值最大;抗剪土樣根系指標于各坡位呈現與抗沖土樣根系指標和全根根系指標相反的規律。由此可見,不同坡位土層厚度的差異,影響著根系生長分布。較厚土層根系舒展從而更加分散、扎根較深,較薄土層(如上坡),根系相對集中于表土層,所以上坡抗剪土樣根量較中、下坡大,各指標值亦較高。 抗剪土樣僅上坡木纖比顯著大于下坡,其他指標均無顯著性差異。抗沖土樣各指標均無顯著性差異。總體上抗剪和抗沖土樣中、下坡根系纖維素含量均較上坡高;而抗剪與抗沖土樣木質素含量排序分別為中坡>上坡>下坡、下坡>中坡>上坡;抗剪土樣木纖比和木質素含量排序一致。結合全根根系的形態指標再次表明,纖維素含量與根系徑級有一定的相關性,中、下坡細根量大,其纖維素含量也表現出更高的水平。 表4 不同坡位抗剪、抗沖土樣根系形態與化學成分指標 表中數據為平均值±標準誤差(n=3),不同小寫字母表示同一指標不同坡位間存在顯著差異(P<0.05) 2.3.2復合體的抗剪強度指標及影響因素 表5中地埂籬地的內摩擦角、粘聚力(下坡除外)均大于對照裸地,僅上坡粘聚力與對照裸地存在顯著性差異。中坡粘聚力顯著小于上、下坡且減幅高達1倍以上。數值上內摩擦角排序:中坡>下坡>上坡,坡位間波動幅度較小,對照裸地與地埂籬地均值分別為85.42°、86.08°;粘聚力排序:上坡>下坡>中坡,對照裸地與地埂籬地均值分別為27.60 kPa、36.92 kPa,地埂籬地較對照粘聚力的增幅最高達45.67%。 相關性分析得出(表6):上坡復合體內摩擦角與RVD、RWD呈顯著和極顯著正相關,上、中坡粘聚力與RVD呈顯著負相關。RLD和RSAD,RVD和RWD總是與抗剪指標有一致或相似的相關關系,如各坡位的內摩擦角與RLD、RSAD都呈負相關,與RVD、RWD均呈正相關。 表5 不同坡位根-土復合體抗剪強度指標和抗沖指標 表中數據為平均值±標準誤差(n=3),不同小寫字母表示同一指標不同坡位間存在顯著差異(P<0.05);不同大寫字母表示同一指標的對照與地埂籬地間存在顯著差異(P<0.05) 圖5 不同坡位土壤抗沖指數動態變化 Fig.5 Dynamic changes of Anti-scourability at different slope portions 2.3.3復合體的抗沖指數及影響因素 抗沖指數均值排序(表5):下坡>中坡>上坡,坡位之間的數值變幅在0.34—5.31之間。差異性檢驗表明,在抗沖指數動態變化過程中t=4 min時各坡位之間有顯著性差異。沖刷過程中各坡位復合體抗沖指數均呈波狀變化(圖5),沖刷初期抗沖指數值最低,上、下坡抗沖指數峰值較中坡延遲4 min,各坡位均呈先增后減再平緩波動上升趨勢。沖刷初期由于表層土壤基本無根系且多浮土所以流失土量大;隨著時間增加抗沖指數先增后減,這個過程主要是根系起到的緩沖作用;隨后泥沙量漸漸變少,抗沖指數基本處于緩升趨勢。整體來看,中、下坡復合體抗沖性能較優,上坡抗沖指數小其波動幅度也較小。 抗沖指數受RSAD和纖維素影響(表6),中、下坡抗沖指數與RSAD呈顯著正相關、與纖維素呈顯著或極顯著負相關關系。纖維素和木質素含量與抗沖指標有相反的相關關系,如抗沖指數與纖維素含量呈負相關,卻與木質素含量呈正相關。 表6 不同坡位根-土復合體抗剪強度指標和抗沖指標與根系形態指標、化學成分相關關系 Table 6 Correlation analysis between shear strength parameters, anti-scourability, and root morphology paratemns, chemical compositions at different slope portions 指標Parameters根長密度Root length density根表面積密度Root surface area density根體積密度Root volume density根重密度Root weight density纖維素Cellulose木質素Lignin木纖比Wood fiber ratio半纖維素Hemicellulose內摩擦角上坡-0.383-0.6980.997*0.998**0.4400.471-0.096-0.781Internal friction angle中坡-0.993-0.8410.7350.897-0.9870.158-0.731-0.753下坡-0.620-0.5160.7800.835-0.0210.030-0.2130.988粘聚力上坡0.6240.869-0.980*-0.960-0.641-0.2480.3290.907Cohesion中坡0.6190.977-0.971*-0.946-0.947-0.3170.610-0.850下坡0.6020.920-0.937-0.908-0.834-0.8610.940-0.660抗沖指數上坡-0.5900.7410.8070.739-0.906*0.7620.6630.140Anti-scourability中坡-0.2000.985*0.9640.971-0.990**0.9550.999*-0.987下坡-0.3900.833*0.8100.775-0.916*0.8850.946-0.116 *表示在0.05的水平上顯著相關,**表示在0.01的水平上極顯著相關 喀斯特地區植物根系為適應土壤和水分異質性強(高位、陡坡土薄缺水)的環境,形成了自身響應機制[7,27]。研究區中拉巴豆根系的形態與分布在不同坡位產生了空間變異。全根和徑級指標均表現為中、下坡大于上坡;0.0 喀斯特坡地土儲量的空間變異可通過根系在土體中的分布深度與集中程度反演獲知。上坡地勢高、水土流失動力潛能大且侵蝕力較強,所以土層最薄,中坡接收上坡流失物質后土層增厚,形成與下坡相似甚至更優的土壤條件。本研究得出抗剪土樣根系總體指標排序均為上坡>下坡>中坡,看似有悖于全根根系形態指標規律,實則反映不同坡位根系分布隨土儲量的變化:上坡土層薄,根系相對集中;中坡土層厚,根系相對舒展分散;下坡與中坡土層厚度相似,根系亦相對分散。 結合土壤物理性質發現,根系與土壤的交互作用形成了一個正反饋機制:根量較大、細根量多、形態指標較優的中、下坡,土壤容重、總孔隙度、飽和含水量和自然含水率均較上坡優,土壤物理性質優良又促進根系的分化發育。由此可見,植物根系形態特征與地勢、土壤厚度、土壤水等有一定的關聯[13,29]。 草本根系抗拉試驗成功率一般在60%—75%[17,25]范圍內,拉巴豆根系試驗成功率為65.52%,相較灌喬木根系成功率高[18,30- 31]。拉巴豆平均抗拉強度29.57 MPa,約為狗牙根(d≤1.5 mm)的1.16倍、馬唐(d≤1.5 mm)的2.03倍[17,32],可見拉巴豆根系有較好的生物力學特性。不同坡位根系抗拉特性指標與直徑的擬合方程雖均為冪函數相關,但其常數存在一定差異,這種現象可歸結為環境影響了粗、細根的比例和根系化學成分含量,同時試驗誤差產生相較真實值一定程度的偏移[30];不同徑級根系抗拉特性指標與化學成分含量的回歸方程常數的差異也存在實驗誤差影響,具體影響因素和變化機制還有待驗證。 本研究中拉巴豆根系極限抗拉力隨直徑增大而增大,抗拉強度隨直徑減小而增強,這與大部分草、灌根系力學特性一致[8,25]。極限抗拉力與木質素含量和木纖比均正相關,抗拉強度與纖維素含量正相關、與木質素含量和木纖比均負相關。大量研究得出同樣結論[8,17,32],而某些研究卻指出抗拉強度與纖維素含量存在負相關關系[15],以及抗拉特性與化學成分無相關關系的結論[33]。關于抗拉特性與化學成分關系不一的本質原因目前存在這樣的解釋:根系維管束中木質部與韌皮部的化學成分含量存在差異,且在增粗生長的動態變化過程中木質部與韌皮部的分化速率不同從而導致不同物種或同一物種不同徑級抗拉特性與化學成分關系的變化[15]。還有研究指出環境因素間接影響木質部中成分含量[13],生長于山脊、坡地等分散地形的根系相較平地、溝谷地等匯聚地形的強度更大,其原理在于較大應力環境下根系木質部的木纖比降低、纖維素局部積累,從而使根系抗拉強度變大[13,33]。 根-土復合體相較對照裸地內摩擦角和粘聚力均增大,內摩擦角排序:中坡>下坡>上坡,粘聚力:上坡>下坡>中坡。大量研究指出草本根系能有效提高土體抗剪強度[9,34]。剪切過程,根系與土顆粒表面摩擦力及顆粒間嵌入作用產生的咬合力影響內摩擦角大小,分泌物與土壤顆粒間膠結作用影響粘聚力大小[9],所以內摩擦角、粘聚力一般隨根長密度、根表面積密度的增大而增大[8- 9]。本研究顯示上坡復合體較另兩坡位粘聚力更大,這與根系形態指標變化一致:上坡復合體根系形態指標均較大,根體積密度、根重密度尤其突出(均較中、下坡大1.5倍以上,且有顯著性差異)。但相關分析得出上坡根系形態指標與粘聚力無顯著相關,甚至負相關(根體積密度),這表明粘聚力并非受根系影響。本研究中拉巴豆生長時間較短,單位土體根量少,對土體的抗剪性能影響未達顯著水平,或土壤類型差異導致粘聚力受其他因素影響更突出。已有研究證明含水率也顯著影響粘聚力,隨含水率增加,粘聚力呈先增后減的變化趨勢(達峰值時含水率約為20%[35])。由中坡根量大但含水率太高(接近峰值含水率的兩倍)最終粘聚力表現最差,猜想本研究的粘聚力受含水率影響更大。 抗沖指數與形態指標和化學成分相關分析表明,其僅受根表面積、纖維素影響。相關研究表明不同背景環境或土壤類型影響抗沖指數的因素存在差異,土壤物理性質(容重、總孔隙度、細砂粒含量等)[36]、有機質[9]、毛根量[34]、毛根形態指標(根長、根表面積等)[37]均被證明與抗沖指數相關。本研究抗沖指數呈波動變化并最終平穩上升,這與楊玉梅等[37]在四川雅安黃壤研究區的結論類似,諶蕓等[10]于紫色土區得出抗沖指數動態變化可用對數函數或二次函數較好擬合,不同土壤類型為產生差異的主要原因。 1)拉巴豆根系全根指標中根長、根表面積、根體積、根重均表現為中坡>下坡>上坡; 0.0 2)抗拉強度與直徑呈冪函數負相關;抗拉強度與纖維素呈對數函數正相關、與木質素和木纖比均呈對數函數負相關;極限抗拉力與木質素呈對數函數正相關、與木纖比呈線性函數正相關。總體上抗拉強度與較細徑級、極限抗拉力與較粗徑級的化學成分相關度更高。 3)復合體根系形態指標均于上坡的抗剪土樣出現最大值,而抗沖土樣于中、下坡值較大。僅上坡抗剪土樣的根體積密度和根重密度顯著大于中、下坡。化學成分指標中抗剪、抗沖土樣纖維素含量均于中、下坡值較高。僅上坡抗剪土樣的木纖比顯著大于中、下坡。 4)復合體內摩擦角、粘聚力均較對照裸地增大,但僅上坡位的粘聚力顯著大于對照。內摩擦角隨根體積密度、根重密度增大而增大,粘聚力隨根體積密度增大而減小。抗沖指數隨根表面積密度增大而增大、隨纖維素含量增大而減小。1.4 數據處理
2 結果與分析
2.1 根系形態與化學成分特征




2.2 單根極限抗拉力與抗拉強度


2.3 根-土復合體抗剪/沖性能特征




3 討論
3.1 拉巴豆根系形態分布及影響因素
3.2 拉巴豆根系抗拉特性及影響因素
3.3 拉巴豆根-土復合體力學特性及影響因素
4 結論