基于CT掃描研究青海湖流域高寒草甸不同坡位土壤大孔隙結構特征①

胡 霞,李宗超，劉 勇，孫貞婷，呂艷麗

(1 北京師范大學環境演變與自然災害教育部重點實驗室，北京 100875；2 北京師范大學減災與應急管理研究院，北京 100875)

?

基于CT掃描研究青海湖流域高寒草甸不同坡位土壤大孔隙結構特征①

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

(1 北京師范大學環境演變與自然災害教育部重點實驗室，北京 100875；2 北京師范大學減災與應急管理研究院，北京 100875)

摘 要：以青海湖流域的高寒草甸土壤為研究對象，對該高寒草甸的坡上、坡中和坡下土壤分別取原狀土柱進行CT掃描，利用Fiji 軟件分析土壤大孔隙結構特征的差異。結果表明：坡位對土壤的大孔隙結構有較大影響，坡上土壤的大孔隙數量、大孔隙度和大孔隙等效直徑均大于坡中土壤，坡中土壤的大孔隙參數大于坡下土壤，坡上土壤的平均大孔隙度是坡中和坡下土壤的15.5和46.5倍。坡上土壤的大孔隙主要分布在150～400 mm 土層深度，而坡中土壤的大孔隙主要分布在0～150 mm深度，坡下土壤的大孔隙主要在0～200 mm深度分布。坡上和坡中土壤的大孔隙形成主要是在土壤團聚體的作用下形成，植物根系在坡下土壤的大孔隙形成中占主導作用。

關鍵詞：高寒草甸；坡位；大孔隙；青海湖；根系

青海湖流域位于青藏高原東北部，是我國西北干旱區、西南高寒區和東部季風區的交匯區，青海湖流域受氣候變化和人類活動的影響，青海湖流域生態惡化，出現草場退化、土地沙漠化面積擴大、水土流失嚴重、生物多樣性減少等一系列生態問題[1]。青海湖流域的土地退化面積超過6.9×105hm2，以每年3% 的速度退化[1]。土地退化會影響土壤結構特征的變化[2-3]，其中，孔隙結構是描述土壤結構特征的重要參數[4]。大孔隙是在土壤漲縮、可溶性物質溶解、凍融循環交替、耕種等物理過程及蚯蚓等動物活動、植物根系生長等生物過程的作用下，土體內形成的使水分和溶質優先遷移的物理孔隙[5-6]。大孔隙的存在可形成大孔隙流或優勢流，使淺層土壤中的水分快速滲入土壤深處或地下水中，加快了地下水響應速度，增加土壤通氣性，提高降雨的入滲率，改變了坡地徑流的形成過程和不同徑流成分的比例[7-10]。

坡面水土流失是導致土壤質量退化及坡面生產力下降的重要原因[11-12]，并造成不同坡位土壤性質變異[13]。研究坡面土壤特征，探討坡面土壤質量退化原因，對其進行植被恢復，將有助于控制土壤退化，提高土地生產力，減緩水土流失，加速退化生態系統的恢復與重建，改善生態環境。目前關于坡面土壤方面的研究主要集中在坡地養分流失規律和形態[14-15]、退耕地土壤性質變化[16]、植被恢復[17]、土壤水文效應[18]等。但是有關不同坡位對土壤結構的影響研究相對較少。本研究選擇青海湖流域高寒草甸的坡上、坡中和坡下土壤，通過CT掃描分析高寒地區不同坡位土壤的大孔隙三維結構特征、孔隙數量、孔隙度特征，研究了青海湖流域高寒草甸不同坡位土壤的大孔隙類型和特征。

1 材料與方法

1.1 研究區概況

青海湖流域位于青藏高原東北部，流域面積29 661 km2，海拔 3 194～5 174 m(圖1)。由于處在中國東部季風區、西北部干旱區和西南部青藏高原高寒區的交匯地帶，再加上自身的湖泊效應，青海湖流域形成了明顯的高寒、少雨、多風、太陽輻射強烈、氣溫日較差大的地區氣候特點。流域內的自然植被包括高寒灌叢、河谷灌叢、沙生灌叢、干旱草甸、高寒草原、溫性草原等。

1.2 樣品采集

選擇青海湖流域的高寒草甸作為研究樣地，樣地位于沙柳河中游，地理位置100.035°E，37.575°N(圖1)。分別在高寒草甸的坡上、坡中和坡下采集原狀土柱，采樣設置3個重復。坡上采樣點的經緯度分別是100.011°E，37.596°N，其主要植被是小嵩草、金露梅、珠芽蓼、燈芯草、黑褐苔草、冷地早熟禾等；坡中采樣點的經緯度分別是100.0074°E，37.595°N，其主要植被是金露梅、高山嵩草、青藏苔草、釘柱委陵菜、珠芽蓼、重齒風毛菊、美麗風毛菊等；坡下采樣點的經緯度分別是100.0051°E， 37.596°N，其主要植被是高山嵩草、青藏苔草、燈芯草、腎形子黃耆、美麗風毛菊、珠芽蓼、扁囊苔草等。土樣采集用內徑100 mm、管壁厚度4 mm、長度500 mm的圓柱狀PVC管，將PVC管的一端打磨成刀口，沿著事先挖好的剖面輕輕敲入土壤中，采樣深度0～500 mm，將其壓入采集地的土壤中，然后整個取出，用海綿包裹。所得PVC管內的原狀土稱作土芯。表1是坡上、坡中和坡下的海拔、植被類型和土壤的基本理化性質。

圖1　研究區域Fig.1 Location of the experimental sites

表1　不同坡位土壤的基本理化性質Table 1 Soil properties of different slope sites

1.3 CT掃描

采用GE公司生產的LightSpreed VCT(64排螺旋CT掃描儀)掃描原狀土柱，掃描電壓為140 kV，電流為200 mA，掃描厚度為0.625 mm，分辨率為0.3 mm，圖像矩陣為512 × 512。掃描前對原狀土柱進行定位，確保土芯在螺旋射線管的中心。

1.4 數據處理

首先通過人工制作大孔隙的方法獲得閾值。選取一根內徑為10 mm 的PVC 管垂直放在內徑100 mm 的PVC 管中央，其周圍裝滿回填土，做成回填土柱。用CT 掃描儀掃描土柱，掃描完的圖片在Fiji 軟件中處理人工制作的大孔隙。設置一個閾值，計算人工大孔隙的大小，與實際的大孔隙大小相比較，如果計算的數值與實際的大小相差較大，再重新設置一個閾值計算人工孔隙大小，直到它們的差值不超過1%，得到圖像分割的閾值(圖2)。

圖2　人工大孔隙的灰度直方圖Fig.2 Gray histogram of artificial macroporosity

將提取的jpg 格式圖像導入Fiji 軟件，用圓形工具提取感興趣區域(直徑為49.11 mm)，避免切入縫隙。用中值濾波法消除圖像噪音，并根據所得到的分割閾值利用全局閾值法處理圖像，得到二值圖像。利用Fiji 軟件對二維圖像的二值圖分析得到孔隙面積、周長以及孔隙圖，利用Fiji 軟件編寫腳本程序，對所有圖像進行批量處理。對圖像去背景處理后，通過3D viewer 查看器通過體積渲染得到三維大孔隙圖，從而更好地觀察孔隙的宏觀結構。

2 結果

2.1 青海湖流域不同坡位土壤大孔隙的三維形態特征

圖3是青海湖流域不同坡位(坡上、坡中和坡下)土壤大孔隙的三維圖。從圖3中可以看出，不同坡位土壤的大孔隙有很大差異。隨著坡位的下降，土壤的大孔隙數量呈現下降趨勢，即坡上土壤的大孔隙明顯高于坡中土壤，坡中土壤的大孔隙明顯高于坡下土壤。而且，坡上土壤的大孔隙主要位于土層中部，坡中和坡下土壤的大孔隙主要分布在土層上部。

Katuwal 等[19]研究表明，生物孔隙(蚯蚓洞穴和根系管道)一般呈管狀，有方向性；Luo等[20]發現由蚯蚓引起的大孔隙一般具有連續性、相對比較大、呈管狀；由根系引起的大孔隙一般具有連續性、圓形，隨著深度的增加孔隙的大小呈下降趨勢；而由團聚體形成(凍融交替或干濕交替)的大孔隙一般很小、分散分布、連續性較差。本研究中的高寒草甸的坡上、坡中土壤的大孔隙在形態特征上都較分散、不連續；而坡下土壤的大孔隙在形態上較細，相對較連通。因此，坡上和坡中土壤的大孔隙主要是在土壤團聚體作用下形成；而坡下土壤的大孔隙主要是根系作用下形成。

2.2 不同坡位土壤大孔隙的數量、孔隙度及孔隙直徑特征

表2是不同坡位土壤大孔隙數量、平均直徑、大孔隙度以及大孔隙度在不同土層的分布。表2數據顯示，青海湖流域高寒草甸坡上土壤的平均大孔隙數量是14，坡中和坡下的大孔隙數量分別是3和2，而且，坡上、坡中和坡下土壤的平均大孔隙數量之間都沒有顯著差異。坡上土壤的大孔隙數量是坡中和坡下土壤的大孔隙數量的4.6和7倍。

表2數據還顯示，坡上土壤的平均大孔隙度是4.65%，坡中和坡下分別是0.30% 和0.10%，坡上土壤的大孔隙度是坡中和坡下土壤的大孔隙度的15.5 和46.5倍。而且，坡上、坡中和坡下土壤的平均大孔隙度之間都沒有顯著差異。坡上土壤在0～100、100～200、200～300、300～400、400～500 mm土層深度，土壤的大孔隙度分別是0.54%、3.56%、17.44%、6.44% 和0.24%；坡中土壤的大孔隙度分別是1.02%、0.21%、0.01%、0.01% 和0.01%；坡下土壤的大孔隙度分別是0.25%、0.17%、0.01%、0.01% 和無。而且，在同一土層深度，坡上、坡中和坡下土壤的大孔隙度沒有明顯差異。

圖4是不同坡位土壤的大孔隙度隨著深度的變化?？梢姡律贤寥赖拇罂紫抖让黠@大于坡中土壤的大孔隙度，坡中土壤的大孔隙度大于坡下土壤的大孔隙度。坡上土壤的大孔隙度在0～150 mm深度呈現緩慢下降的趨勢，在150～400 mm深度，土壤的大孔隙度呈現急劇上升的趨勢；對于坡中土壤，在土層深度0～150 mm，大孔隙度呈現緩慢下降的趨勢，在150～400 mm深度，土壤的大孔隙度保持穩定；坡下土壤的大孔隙度隨著土壤深度的變化趨勢與坡中土壤的變化趨勢類似，在0～200 mm深度，土壤的大孔隙度呈現緩慢下降的趨勢。

從表2中可以看出，坡上土壤的平均大孔隙直徑是1.63 mm，坡中和坡下土壤的大孔隙直徑是0.76 mm和0.40 mm，且坡上、坡中和坡下土壤的大孔隙平均直徑沒有顯著差異。圖5是坡上、坡中和坡下土壤不同直徑大孔隙所占比例。從圖5中可以看出，3個大孔隙直徑(0＜D＜1 mm，1＜D＜2 mm和2＜D＜3 mm)在所有大孔隙直徑中占所有大孔隙的85% 以上，本研究結果與前人的研究結果一致，土壤中大孔隙直徑主要是2 mm和3 mm直徑，這兩個直徑大孔隙占所有大孔隙的50% 以上[21]。

坡上土壤的大孔隙度在0～150 mm深度較小，在150～400 mm孔隙度急劇增加，說明坡上土壤在表面主要以坡面徑流為主，水分入滲后主要以大孔隙流的形式向下運移。坡中和坡下部分的大孔隙度主要集中于土壤表層(0～200 mm)，說明坡中和坡下表層土壤入滲能力較強，水分主要以大孔隙流的形式向下運移，土壤入滲后可形成壤中流。而較深層(200～450 mm)土壤孔隙度較低，說明水分在較深層土壤運移主要以基流的形式進行，流速較慢。

圖3　不同坡位土壤大孔隙的三維圖Fig.3 3D visualization of soil macropore networks in soil columns in different slope sites

表2　不同坡位土壤的大孔隙特征Table 2 Numbers,diameters and porosities of CT-measured macropores at different soil depths and slope sites

圖4　不同坡位土壤大孔隙度隨著深度的變化Fig.4 Macroporosity distribution along soil column depth in different slope sites

圖5　不同坡位土壤大孔隙的直徑分布Fig.5 Diameter distribution of macropores in different slope sites

3 結論

青海湖流域高寒草甸坡位對土壤大孔隙的影響較大，坡上土壤的大孔隙度明顯大于坡中和坡下土壤的大孔隙度，坡中土壤的大孔隙度大于坡下土壤的大孔隙度。坡上土壤的大孔隙主要位于150～400 mm深度土層，坡中土壤的大孔隙主要位于0～150 mm深度土層，坡下土壤的大孔隙在深度0～200 mm深土層度分布較多。坡上土壤的平均大孔隙直徑明顯大于坡中和坡下土壤。青海湖流域高寒草甸坡上和坡中土壤的大孔隙主要是在土壤團聚體的作用下形成，植物的根系在坡下土壤的大孔隙形成中占主要作用。

參考文獻：

[1]肖青寧.青海湖環湖地區生態環境問題及建議[J].青海草業,2002,11(4):20-23

[2]Sharma P,Aggarwal K.Soil structure under different land uses[J].Catena,1984,11:197-200

[3]李宗超,胡霞.小葉錦雞兒灌叢化對退化沙質草地土壤孔隙特征的影響[J].土壤學報,2015,52(1):242-248

[4]Brewer.Fabric and Mineral Analysis of Soils[M].Company Huntington,NY:Krieger Publishing,1976

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

[6]馮杰,郝振純.水及溶質在有大孔隙的土壤中運移機制研究進展[J].河海大學學報,2002,30(2):63-69

[7]Germann P,Beven K.Water flow in soil macropores.I.An experimental approach[J].Journal of Soil Science,1981,32:1-13

[8]Moidrup P,Olesen T,Sconning P,et al.Predicting the gas diffusion coefficient in undisturbed soil from soil water characteristics[J].Soil Science Society of America Journal,2000,64(8):94-100

[9]Coppola A,Kutílek M,Frind E O.Transport in preferential flow domains of the soil porous system:Measurement,interpretation,modelling,and upscaling[J].Journal of Contaminant Hydrology,2009,104:1-3

[10]盛豐,張利勇,王康.土壤大孔隙發育特征對水和溶質輸移的影響[J].土壤,2015,47(5):1 007-1 013

[11]史衍璽,唐克麗.人為加速侵蝕下土壤質量的生物學特性變化[J].水土保持學報,1998,4(1):28-31

[12]張少良,張興義,劉曉冰,等.典型黑土區不同坡位剖面土壤速效鉀空間分布規律研究.土壤,2014,46(2):218-224

[13]Miller P M,Singer M J,Nielsen D R.Spatial variability of wheat yield and soil properties on complex hills[J].Soil Science Society of America Journal,1988,52:1 133-1 141

[14]Flanagan D C,Foster G R.Storm pattern effect on nitrogen and phosphorus losses in surface runoff[J].Transactions of the ASAE,1989,32:535-544

[15]Sharpley A N,Smith S J,Jones O R,et al.The transport of bio2available phosphorus inagricultural runoff[J].Journal Environment Quality,1992,21:30-351

[16]彭文英,張科利,陳瑤,等.黃土坡耕地退耕還林后土壤性質變化研究[J].自然資源學報,2005,20(2):272-278

[17]萬雪琴,胡庭興,張健,等.坡耕地退耕還林后的植被恢復[J].林業科學,2005,41(2):191-194

[18]王國梁,劉國彬,常欣,等.黃土丘陵區小流域植被建設的土壤水文效應[J].自然資源學報,2002,3(17):339-345

[19]Katuwal S,Norgaard T,Moldrup P,et al.Linking air and water transport in intact soils to macropore characteristics inferred from X-ray computed tomography[J].Geoderma,2015,237-238:9-20

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

[21]Iversen B V,Lamande M,Torp S B,et al.Macropores and macropore transport:Relating basic soil properties to macropore density and soil hydraulic properties[J].Soil Sci,2012,177(9):535-542

Soil Macropores of Alpine Kobresia Meadow in Different Slope Positions in the Qinghai Lake Watershed

HU Xia1,2,LI Zongchao1,2,LIU Yong1,2,SUN Zhenting1,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)

Abstract:The objective of this study is to quantify macropore structures of Alpine Kobresia meadow soils in different slope positions in the Qinghai Lake Watershed.A total of nine soil cores(0-50 cm deep)were taken respectively at the upper-slope,mid-slope and lower-slope sites with 3 replicates in each site,and the cores were scanned with a GE HISPEED FX/I medical scanner.Numbers of macropores,macroporosities and macropore equivalent diameters were interpreted by using X-ray computed tomography.The results indicated that slope position could significantly influence soil macropores.Soils on upper-slope had greater macroporosity,deeper and longer macropores than other two sites.Macroporosity was 15.5-46.5 times greater in soils on upper-slope than other two sites,and macropores were distributed mainly in the 150-400 mm layer,but mainly in the 0-150 mm and 0-200 mm soil layer in the mid-slope and the lower-slope,respectively.The form of macropores in lower slope can be attributed to root growth and decay,while can be attributed to aggregate function in the upper and middle slopes.

Key words:Alpine Kobresia meadow; Slope position; Macropore; Qinghai Lake; Root

作者簡介：胡霞(1978—)，女，江蘇淮安人，博士，副教授，主要從事土壤物理及其生態功能研究。E-mail:huxia@bnu.edu.cn

基金項目：①國家自然科學基金項目(41471018)和國家自然科學基金重點項目(41130640)資助。

DOI:10.13758/j.cnki.tr.2016.01.027

中圖分類號：S157.1