劉 勇，胡 霞*，李宗超，孫貞婷，程亞倩，呂艷麗

(1 北京師范大學(xué)環(huán)境演變與自然災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室，北京 100875；2 北京師范大學(xué)減災(zāi)與應(yīng)急管理研究院，北京 100875)

基于CT的青海湖流域芨芨草草地土壤大孔隙特征分析①

劉 勇1,2，胡 霞1,2*，李宗超1,2，孫貞婷1,2，程亞倩1,2，呂艷麗1,2

(1 北京師范大學(xué)環(huán)境演變與自然災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室，北京 100875；2 北京師范大學(xué)減災(zāi)與應(yīng)急管理研究院，北京 100875)

為探究青海湖流域芨芨草草地的土壤大孔隙特征，以青海湖流域芨芨草草地為研究對(duì)象，在研究區(qū)芨芨草斑塊處和芨芨草斑塊間分別采集原狀土柱，每塊樣地選取3個(gè)重復(fù)樣品，通過對(duì)樣品進(jìn)行CT掃描并利用Fiji軟件解譯的方法分析土壤孔隙的結(jié)構(gòu)特征差異。結(jié)果表明：芨芨草斑塊土壤平均大孔隙數(shù)量、不同深度大孔隙數(shù)量、平均大孔隙度以及不同深度大孔隙度均大于芨芨草斑塊間土壤。芨芨草斑塊土壤大孔隙平均數(shù)量是芨芨草斑塊之間的2.8倍，平均大孔隙度是芨芨草斑塊間的12.1倍，芨芨草斑塊的土壤大孔隙度隨土壤深度增加呈先增加后減少的趨勢(shì)，斑塊間的大孔隙度呈遞減趨勢(shì)。

青海湖流域；芨芨草；CT；土壤大孔隙

青海湖是我國面積最大的內(nèi)陸咸水湖泊，湖泊面積 4 318 km2(2008 年)，流域總面積29 661 km2，其流域地形復(fù)雜，加上湖泊河流等水文因素作用，因此該地區(qū)土壤植被類型繁多[1]。

芨芨草屬禾本科芨芨草屬植物[2]，廣泛分布于青海湖盆地、共和盆地以及柴達(dá)木盆地等青海東部地區(qū)[3]，具有耐寒、耐旱、耐鹽、適應(yīng)性強(qiáng)等特點(diǎn)[4]，是重要的牧草、纖維植物和水土保持植物[5]。陳桂琛和彭敏[3]對(duì)青海地區(qū)芨芨草群落的特征及分布規(guī)律進(jìn)行研究后發(fā)現(xiàn)，該地區(qū)分布的芨芨草群落特征和生境條件具有草原植被特征和地帶性意義。水分是限制芨芨草草原生態(tài)系統(tǒng)發(fā)育的關(guān)鍵因子，小降水事件(<5 mm)是該生態(tài)系統(tǒng)水分補(bǔ)給的重要來源，而且能被該草原生態(tài)系統(tǒng)高效利用，避免干旱脅迫[6]。草地退化能夠?qū)е峦寥揽紫抖葴p少并間接導(dǎo)致0 ~ 5 cm土壤容積含水量極顯著升高，20 cm以下顯著降低[7]。一些研究發(fā)現(xiàn)芨芨草能顯著提高草原土壤速效養(yǎng)分含量以及土壤的陽離子代換量[8]，土壤理化性狀的改變會(huì)影響到生態(tài)系統(tǒng)功能和效益的發(fā)揮[9]。到目前為止，國內(nèi)對(duì)芨芨草的研究主要集中在芨芨草草原生態(tài)系統(tǒng)的空間分布格局、物種多樣性以及生態(tài)功能等方面[10–11]，對(duì)于芨芨草地土壤結(jié)構(gòu)的研究，尤其是土壤孔隙特征的研究相對(duì)較少。

土壤大孔隙是土壤結(jié)構(gòu)特征的重要組成部分。大孔隙是在土壤漲縮、可溶性物質(zhì)溶解、凍融循環(huán)交替、耕種等物理過程及蚯蚓動(dòng)物活動(dòng)、植物根系生長(zhǎng)等生物過程的作用下，土體內(nèi)形成使水分和溶質(zhì)優(yōu)先遷移的物理孔隙[12]。研究表明，很多因素會(huì)對(duì)土壤大孔隙特征產(chǎn)生影響，凍融循環(huán)能夠顯著降低40 ~ 80 cm土層范圍內(nèi)的平均大孔隙面積以及Feret直徑[13]，小葉錦雞兒灌叢化對(duì)土壤的孔隙結(jié)構(gòu)有較大影響[14]，坡位能夠?qū)е峦寥来罂紫督Y(jié)構(gòu)特征產(chǎn)生差異[15]。土壤大孔隙直接影響土壤中水分的運(yùn)移方式和途徑[16]，形成孔隙流或優(yōu)勢(shì)流，加快淺層土壤中水分滲入土壤深處速率，提高土壤通透性和降水入滲率，改變土壤徑流形成過程和比例[17]。為評(píng)價(jià)大孔隙對(duì)水流及溶質(zhì)運(yùn)移特性的影響，需確定土壤中大孔隙的大小、數(shù)目、分布和連通性[18]。目前對(duì)大孔隙研究的方法主要有染色法、切片法、CT掃描法[19]，其中 CT掃描技術(shù)操作簡(jiǎn)單、迅速，可以保持土壤樣品原狀結(jié)構(gòu)，通過三維技術(shù)對(duì)土壤孔隙網(wǎng)絡(luò)結(jié)構(gòu)定量分析，是無損檢測(cè)土壤結(jié)構(gòu)的有效手段。

本文選擇青海湖流域的芨芨草草地為研究樣地，通過CT掃描技術(shù)確定芨芨草草地土壤的孔隙結(jié)構(gòu)、孔隙數(shù)量以及垂直分布規(guī)律，對(duì)比分析芨芨草斑塊和芨芨草斑塊間土壤孔隙特征，從而揭示芨芨草對(duì)土壤大孔隙特征的影響，為合理解釋土壤大孔隙對(duì)水分運(yùn)移以及優(yōu)勢(shì)物種芨芨草對(duì)水分利用的機(jī)理提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)域概況

研究區(qū)域選在青海湖流域沙柳河下游地區(qū)(圖1)，該地區(qū)氣候?qū)俑咴箨懶詺夂?，常年低溫、干旱少雨，太陽輻射?qiáng)烈，氣溫日較差大，年平均降水量為291 ~ 571 mm，年蒸發(fā)量為1 200 ~ 2 000 mm。沙柳河下游地區(qū)土壤為栗鈣土，主要植被類型為芨芨草、河谷灌木具鱗水柏枝以及農(nóng)作物(燕麥、青稞等)[20]。

圖1 青海湖流域樣地Fig. 1 Location of sample plots of Qinghai Lake valley

1.2 樣品采集

選擇芨芨草草原作為研究樣地，分別在芨芨草斑塊和芨芨草斑塊間用PVC管采集原狀土樣，并設(shè)置3個(gè)重復(fù)。PVC管直徑為10 cm，壁厚為4 mm，長(zhǎng)度為50 cm，將PVC管一端套住鐵帽，將其緩慢壓入采樣點(diǎn)土壤中，然后挖開旁邊土壤，將PVC管連同管內(nèi)土壤整個(gè)取出，空隙部分用海綿填充包裹，避免土壤散開，管內(nèi)土壤則為原狀土樣。通過室內(nèi)分析測(cè)定得到土壤基本理化性質(zhì)如表 1 所示。

表1 芨芨草草地土壤基本理化性質(zhì)Table 1 Soil properties within Achnatherum splendens patch

1.3 CT掃描與圖像處理

原狀土樣采用Light Speed VCT(64排多層螺旋CT掃描儀)掃描，該CT應(yīng)用在醫(yī)學(xué)領(lǐng)域，為保證掃描土壤樣品圖像的質(zhì)量和精度，需重新設(shè)定掃描參數(shù)。經(jīng)調(diào)整，峰值電壓為130 kV，檢測(cè)電流為200 mA，螺距為0.948︰1。掃描圖像矩陣為512 × 512，掃描層厚為0.625 mm/層，體素大小0.33 mm × 0.33 mm × 0.625 mm。

掃描圖像為tif格式，導(dǎo)入Fiji軟件，為了將土壤部分與PVC管壁分開，因此需要利用圓形工具對(duì)圖片進(jìn)行切割。經(jīng)多次對(duì)比分析確定合適閾值，對(duì)圖像進(jìn)行二值化處理，黑色代表孔隙，白色部分代表土壤或者其他介質(zhì)(圖 2)。根據(jù)二值化圖像提取大孔隙輪廓，然后對(duì)所有圖片進(jìn)行堆棧得到孔隙三維圖像，對(duì)二值化圖像進(jìn)行分析可以計(jì)算出孔隙面積、周長(zhǎng)等特征參數(shù)。

圖2 圖像處理過程Fig. 2 Image processing course

2 結(jié)果與分析

2.1 土壤孔隙三維形態(tài)特征

圖 3是青海湖流域芨芨草草地土壤大孔隙的三維圖。從圖中可以看出，芨芨草斑塊土壤的大孔隙明顯高于芨芨草斑塊之間土壤。而且，芨芨草斑塊土壤的大孔隙在整個(gè)土層都有分布，斑塊之間土壤的大孔隙主要分布在土層上部。

圖3 芨芨草斑塊和斑塊間土壤孔隙三維可視化圖Fig. 3 3D visualization of soil macropores within (up) and outside (down) Achnatherum splendens patch

芨芨草斑塊土壤的大孔隙在形態(tài)特征上都較連續(xù)；而芨芨草斑塊之間的土壤大孔隙在形態(tài)上較分散、不連續(xù)。一般情況下，蚯蚓等生物形成大孔隙且相對(duì)連續(xù)、呈管狀；而由根系形成的孔隙也相對(duì)比較連續(xù)、圓形，且隨深度增加大孔隙呈減小趨勢(shì)[21]，觀察芨芨草斑塊下土壤剖面和周邊環(huán)境，均未發(fā)現(xiàn)大量土壤動(dòng)物活動(dòng)，從土壤剖面觀測(cè)到芨芨草的主根主要分布在15 ~ 25 cm土層，三維圖中大孔隙分布呈向下擴(kuò)散狀且連續(xù)性較好，因此可以判斷該地土壤孔隙的形成主要是根系作用。從芨芨草斑塊間土壤孔隙三維圖可以看出，大孔隙在表層0 ~ 10 cm土壤分布較為密集，在20 ~ 30 cm深度處也有部分孔隙存在。

2.2 大孔隙數(shù)量

Warner等[22]在X射線攝影技術(shù)研究大孔隙特征的應(yīng)用中，將直徑比團(tuán)聚體間縫隙大的且孔徑大于1 mm的孔隙稱作大孔隙；Luxmoore等[23]以毛管模型為依據(jù)，推算出大孔隙當(dāng)量孔徑大于1 mm。因此本文選擇孔徑大于1 mm的孔隙作為土壤大孔隙，利用Fiji軟件，在二值化圖像中提取等效直徑大于1 mm的孔隙，可以計(jì)算出大孔隙數(shù)量。

表 2是芨芨草斑塊和斑塊之間的土壤平均大孔隙數(shù)量及不同深度土壤大孔隙平均數(shù)量。從表2可以看出，芨芨草斑塊的土壤大孔隙數(shù)量明顯大于芨芨草斑塊之間土壤，芨芨草斑塊下的大孔隙平均數(shù)量是芨芨草斑塊之間的2.8倍，其數(shù)值分別是31和11個(gè)。芨芨草斑塊下土壤大孔隙主要分布在 0 ~ 20 cm土層，尤其在0 ~ 10 cm深度土層大孔隙數(shù)量最多，占總大孔隙數(shù)量的41.17%，且隨土壤深度的增加大孔隙數(shù)量逐漸減少；芨芨草斑塊間土壤大孔隙數(shù)量在土壤表層最多，占總大孔隙數(shù)量的47.30%，在20 ~ 30 cm土層大孔隙數(shù)量要多于10 ~ 20 cm和30 ~ 40 cm土層，因此，芨芨草斑塊間的土壤大孔隙數(shù)量在土層中呈現(xiàn)先增加后減少的現(xiàn)象。

表2 芨芨草斑塊和芨芨草斑塊間土壤的大孔隙數(shù)量Table 2 Numbers of soil macropores within and outisde Achnatherum splendens patch

2.3 大孔隙度特征

芨芨草斑塊各個(gè)土層的土壤大孔隙度均大于斑塊間土壤大孔隙度，且平均大孔隙度是芨芨草斑塊間的12.1倍，其中芨芨草斑塊在0 ~ 10、10 ~ 20、20 ~ 30、30 ~ 40 cm土層的大孔隙度分別是斑塊間的28.1、6.67、1.1、2.1倍。 在土壤0 ~ 10 cm深度，芨芨草斑塊和芨芨草斑塊之間的大孔隙度存在顯著性差異，且芨芨草斑塊的土壤大孔隙度明顯大于芨芨草斑塊間；在土壤10 ~ 20、20 ~ 30和30 ~ 40 cm深度，芨芨草斑塊和芨芨草斑塊之間的大孔隙度沒有顯著性差異，但是前者的大孔隙度明顯大于后者。

圖4是土壤大孔隙度隨土壤深度變化圖，從圖中可以看出，芨芨草斑塊3個(gè)重復(fù)的土壤大孔隙度隨土壤深度變化的趨勢(shì)相同，隨著土壤深度的增加大孔隙度先增加后減少，重復(fù)2和重復(fù)3大孔隙度的峰值出現(xiàn)在3 cm處，重復(fù)1的峰值出現(xiàn)在8.8 cm處，重復(fù)2和重復(fù)3的大孔隙度下降，在深度10 cm處以下趨于平緩，重復(fù)1的大孔隙度在土深20 cm處趨于平緩。

表3 芨芨草斑塊和芨芨草斑塊間土壤大孔隙度(%)Table 3 Mean macroporosity and macroporosities in different depths within and outside Achnatherum splendens patch

圖4 芨芨草斑塊和芨芨草斑塊間土壤的土壤大孔隙度隨土壤深度變化特征Fig. 4 Macroporosity distributions in soil profiles within (left) and outside (right) Achnatherum splendens patches

芨芨草斑塊間的 3個(gè)重復(fù)的土壤大孔隙度隨土壤深度的變化一致，隨土壤深度增加大孔隙度呈減少趨勢(shì)，3個(gè)重復(fù)在0 ~ 2.5 cm的范圍內(nèi)大孔隙度隨深度的增加迅速遞減，重復(fù)1在16 ~ 32 cm范圍內(nèi)的大孔隙度增加并出現(xiàn)波動(dòng)。對(duì)比芨芨草斑塊和斑塊之間的大孔隙度隨土壤深度變化圖，能夠明顯看出芨芨草斑塊下的大孔隙度較大且隨土壤深度變化呈先增加后減少的趨勢(shì)；而芨芨草斑塊間的土壤大孔隙度較小且隨土壤深度變化呈先遞減后平穩(wěn)的趨勢(shì)。從芨芨草斑塊與芨芨草斑塊間土壤大孔隙特征的差別來看，主要原因可能是植被的根系作用。

3 結(jié)論

從大孔隙的三維形態(tài)上來看，芨芨草斑塊土壤的大孔隙在形態(tài)特征上較連續(xù)，而芨芨草斑塊之間的土壤大孔隙在形態(tài)上較分散、不連續(xù)。芨芨草斑塊土壤大孔隙平均數(shù)量和大孔隙度明顯高于芨芨草斑塊間土壤，且芨芨草斑塊土壤的大孔隙的連通度高于斑塊間土壤。芨芨草斑塊土壤大孔隙主要分布在0 ~ 20 cm土層，尤其在0 ~ 10 cm深度土層大孔隙數(shù)量最多，芨芨草斑塊間土壤大孔隙數(shù)量在土壤表層最多。芨芨草斑塊和芨芨草斑塊之間的10 ~ 40 cm深度土壤沒有顯著性差異，但前者大孔隙度明顯大于后者。芨芨草斑塊土壤大孔隙度隨著土壤深度的增加先增加后減少，芨芨草斑塊間土壤大孔隙度隨土壤深度增加呈減少趨勢(shì)。

[1] 劉吉峰, 李世杰, 秦寧生, 等. 青海湖流域土壤可蝕性K 值研究[J]. 干旱區(qū)地理, 2006, 29(3): 321–326

[2] 劉尚武, 中國科學(xué)院, 西北高原生物研究所, 等. 青海植物志. 第二卷[M]. 西寧: 青海人民出版社, 1999

[3] 陳桂琛, 彭敏. 青海省芨芨草草原的群落特征及其分布規(guī)律[J]. 西北植物學(xué)報(bào), 1993, 13(2): 154–162

[4] 張雨軍. 對(duì)芨芨草草場(chǎng)資源在鄂爾多斯市開發(fā)利用的建議[J]. 內(nèi)蒙古草業(yè), 2002, 14(2): 28–29

[5] 衛(wèi)東, 王彥榮. 芨芨草種子發(fā)芽檢驗(yàn)方法的研究[J]. 草業(yè)科學(xué), 1998, 15(4): 29–32

[6] Ehleringer J, Dawson T. Water uptake by plants: Perspectives from stable isotope composition[J]. Plant, Cell & Environment, 1992, 15(9): 1073–1082

[7] 張法偉, 郭竹筠, 李以康, 等. 青海湖芨芨草干草原植被退化對(duì)土壤溫濕特征的影響[J]. 干旱區(qū)研究, 2013, 30(2): 219–225

[8] 王庫, 于天富. 水土保持植物——芨芨草對(duì)土壤養(yǎng)分的影響[J]. 土壤, 2002, 34(3): 170–172

[9] 郭忠升, 邵明安. 半干旱區(qū)人工林草地土壤旱化與土壤水分植被承載力[J]. 生態(tài)學(xué)報(bào), 2003, 23(8): 1640–1647

[10] 周國英, 陳桂琛, 魏國良, 等. 青海湖地區(qū)芨芨草群落主要種群分布格局研究[J]. 西北植物學(xué)報(bào), 2006, 26(3): 579–584

[11] 周國英, 陳桂琛, 趙以蓮, 等. 青海湖地區(qū)芨芨草群落特征及其物種多樣性研究[J]. 西北植物學(xué)報(bào), 2003, 23(11): 1956–1962

[12] 馮杰, 郝振純. CT 掃描確定土壤大孔隙分布[J]. 水科學(xué)進(jìn)展, 2002, 13(5): 611–617

[13] 王恩姮, 盧倩倩, 陳祥偉. 模擬凍融循環(huán)對(duì)黑土剖面大孔隙特征的影響[J]. 土壤學(xué)報(bào), 2014, 51(3): 490–496

[14] 李宗超, 胡霞. 小葉錦雞兒灌叢化對(duì)退化沙質(zhì)草地土壤孔隙特征的影響[J]. 土壤學(xué)報(bào), 2015, 42(1): 242–248

[15] 胡霞, 李宗超, 劉勇, 等. 基于 CT掃描研究青海湖流域高寒草甸不同坡位土壤大孔隙結(jié)構(gòu)特征[J]. 土壤, 2016, 48(1): 180–185

[16] Dunkerley D. Assessing the influence of shrubs and their interspaces on enhancing infiltration in an arid australian shrubland[J]. The Rangeland Journal, 2000, 22(1): 58–71

[17] Germann P, Beven K. Water flow in soil macropores i. An experimental approach[J]. Journal of Soil Science, 1981, 32(1): 1–13

[18] 馮杰, 郝振純. CT 在土壤大孔隙研究中的應(yīng)用評(píng)述[J].灌溉排水, 2000, 19(3): 71–76

[19] 劉偉, 區(qū)自清, 應(yīng)佩峰. 土壤大孔隙及其研究方法[J].應(yīng)用生態(tài)學(xué)報(bào), 2001, 12(3): 465–468

[20] 吳華武, 李小雁, 蔣志云, 等. 基于dD和d~(18)O的青海湖流域芨芨草水分利用來源變化研究[J]. 生態(tài)學(xué)報(bào), 2015, 35(24): 8174–8183

[21] Luo L, Lin H, Li S. Quantification of 3-d soil macropore networks in different soil types and land uses using computed tomography[J]. Journal of Hydrology, 2010, 393(1): 53–64

[22] Warner G S, Nieber J L, Moore I D, et al. Characterizing macropores in soil by computed tomography[J]. Soil Science Society of America Journal, 1989, 53(3): 653–660

[23] Luxmoore R J, Jardine P M, Wilson G V, et al. Physical and chemical controls of preferred path flow through a forested hillslope[J]. Geoderma, 1990, 46(1): 139–154

Soil Macropore Characteristics of Achnatherum splendens Steppe in the Qinghai Lake Watershed

LIU Yong1,2, HU Xia1,2*, LI Zongchao1,2, SUN Zhenting1,2, CHENG Yaqian1,2, LV Yanli1,2
(1 Key Laboratory of Environmental Change and Natural Disaster of Ministry of Education, Beijing Normal University, Beijing 100875, China; 2 Academy of Disaster Reduction and Emergency Management, Beijing Normal University, Beijing 100875, China)

Soil macropore characteristics of Achnatherum splendens Steppe in the Qinghai Lake Watershed of northeastern Qinghai-Tibet Plateau by using computed tomography (CT) was quantified. A plot within a Achnatherum splendens patch and a plot outside a Achnatherum splendens patch were selected, and 3 soil cores (0 – 50 cm deep) were taken as replicates from each plot and scanned with a GE HISPEED FX/I medical scanner. The results indicated that the number of macropores and the macroporosity of soil within the Achnatherum splendens patch were larger than those of soil outside the patch, the mean numbers of macropore and macroporosity of soil within the Achnatherum splendens patch were 2.8 times and 12.1 times greater than those in soil outside the patch. With the increase of soil depth, the macroporosity increased first and then decreased within the Achnatherum splendens patch, while it decreased continuously in soil outside the patch.

Achnatherum splendens; CT; Soil macropores

S157.1

A

10.13758/j.cnki.tr.2017.01.027

國家自然科學(xué)基金項(xiàng)目(41471018)和國家自然科學(xué)基金重點(diǎn)項(xiàng)目(41130640)資助。

* 通訊作者(huxia@bnu.edu.cn)

劉勇(1990—)，男，山東棗莊人，碩士研究生，主要從事土壤物理及其生態(tài)功能研究。E-mail: esc1909@163.com