貝殼堤島3種植被類型的土壤顆粒分形及水分生態特征

夏江寶，張淑勇，王榮榮，趙艷云，孫景寬，劉京濤，劉 慶

(1. 濱州學院山東省黃河三角洲生態環境重點實驗室， 濱州 256603；2. 山東省黃河三角洲野生植物資源開發利用工程技術研究中心， 濱州 256603； 3. 山東農業大學林學院， 泰安 271018)

貝殼堤島3種植被類型的土壤顆粒分形及水分生態特征

夏江寶1,2，張淑勇1,3，*，王榮榮1,3，趙艷云1，孫景寬1，劉京濤1，劉 慶1

(1. 濱州學院山東省黃河三角洲生態環境重點實驗室， 濱州 256603；2. 山東省黃河三角洲野生植物資源開發利用工程技術研究中心， 濱州 256603； 3. 山東農業大學林學院， 泰安 271018)

運用土壤粒徑質量分布原理與分形學理論，以黃河三角洲貝殼堤島的3種植被類型為研究對象，以裸地為對照，測定分析土壤顆粒分形維數、粒徑組成和水分物理參數，探討不同植被類型的土壤顆粒分形特征及其影響因素。結果表明：貝殼砂土壤中粗砂粒含量最高；其次是細砂粒，而石礫和粉粘粒含量較低。灌木林地和草地具有降低石礫、粗砂粒，增加細砂粒和粉粘粒含量的作用。不同植被類型土壤顆粒分形維數均值在1.5845—1.9157之間，大小依次為酸棗林、杠柳林和草地，表層高于20—40cm土層。酸棗林、杠柳林及草地0—40cm土壤容重均值分別比裸地低23.87%，14.51%和10.47%；總孔隙度均值分別比裸地增加16.96%，16.71%和1.31%。植被恢復措施對貝殼砂表層的孔隙結構、疏松程度改善較好，草地及灌木林地的蓄水性能均表現為0—20cm高于20—40cm。土壤顆粒分形維數與粉粘粒含量、毛管孔隙度、總孔隙度、飽和蓄水量、吸持蓄水量等呈極顯著正相關，與細砂粒含量呈顯著正相關，與粗砂粒含量和容重呈極顯著負相關，與石礫含量、非毛管孔隙度、滯留蓄水量的相關性不大。從土壤分形維數及其水分生態特征來看，貝殼堤島3種植被類型的改良土壤物理性質及蓄水保土功能表現為灌木林好于草地，其中酸棗林好于杠柳林，0—20cm好于20—40cm土層。

分形維數；土壤粒徑分布；土壤水分；植被類型；貝殼堤島

土壤顆粒分形特征不但表征土壤粒徑的大小組成及孔隙分布[1- 2]，還能反映土壤水力學特征[3]、土壤質地均勻程度及通氣透水性[4- 5]及土壤肥力[6]等特性。利用粒徑的重量分布描述土壤顆粒組成的分形維數具有求解精確、簡便的特點[7- 9]，因此，這一方法得到廣泛應用。目前對土壤顆粒分形特征的研究已有傳統的分形維數與不同土壤質地關系的分析[8- 9]集中到某一具體立地類型、同種土壤質地下的不同土地利用方式或植被恢復措施效益分析及土壤質量評價等方面，如用土壤顆粒分形特征反映科爾沁沙地農田沙漠化演變過程[10]、庫布齊沙漠沙柳沙障構建方式優劣[11]、沂蒙山區的植被恢復效果[12]、土石山區林地土壤質量評價[5,13]、黃河三角洲灘地[6]及黃土丘陵溝壑區[14]土地利用類型比較、退耕還湖安慶沿江濕地土壤演變狀況[15]等，同時也從單一的土壤顆粒分形向不同土壤類型的多重分形維數轉變[14,16- 19]。

貝殼堤是淤泥質或粉砂質海岸所特有的一種灘脊類型，黃河三角洲的濱州市無棣、沾化縣境內分布著兩道貝殼堤，與天津、河北的貝殼堤相連，組成規模宏大、國內獨有的貝殼灘脊海岸，與美國圣路易斯安娜州和南美蘇里南的貝殼堤并稱為世界三大古貝殼堤，在世界第四紀地質和海岸地貌研究中占有極其重要的位置[20]。貝殼堤島這一獨特的生態系統，引起了眾多學者的關注，目前對黃河三角洲貝殼堤島的研究主要集中在貝殼堤島脆弱生態系統特征及其保護管理對策[20- 21]、植被及微生物分布特征[22]、典型灌草生理生態特征[23- 24]及貝殼沙中微量元素含量和形態特征[25]等方面，而對該區域主要灌草植被恢復措施下的貝殼砂土壤顆粒分形特征及其影響因素的分析尚未見報道。鑒于此，本文以黃河三角洲貝殼堤島的杠柳(PeriplocasepiumBunge)林、酸棗(Ziziphusjujubavar.spinosaHu)林及砂引草(MesserschmidiasibiricaLinn.)草地3種植被類型為研究對象，并以裸地為對照，運用土壤單重分形學原理與方法，測定分析貝殼堤島不同植被類型的土壤顆粒分形維數、顆粒組成、容重、孔隙度及蓄水性能等指標，探討不同植被類型對土壤物理結構的改良作用，闡明貝殼砂植被恢復措施下的土壤顆粒分形特征及其影響因素，以其為貝殼砂生境下的土壤顆粒分形學機制奠定基礎，為黃河三角洲貝殼堤島灌草種類選擇及模式構建提供理論依據和技術參考。

1 研究地概況

圖1 貝殼堤島不同植被類型的土壤采樣點示意圖Fig.1 Schematic map of sampling sites of soil under different vegetation types in Shell Island

本研究選擇在山東濱州市無棣縣中東部濱海低地的汪子島(38°02′51″—38°21′06″N，117°46′58″—118°05′43″E)，總面積約435.4km2，位于濱州貝殼堤島與濕地國家級自然保護區的緩沖區內(圖1)。該保護區屬于暖溫帶東亞季風大陸性半濕潤氣候區，多年平均降水量為550mm，主要集中在6—9月份，占全年降水量的71%以上；多年平均蒸發量為2430.6 mm，蒸降比為4.4。多年平均氣溫為12.36℃，極端最低氣溫為-25.3℃，極端最高氣溫為37.5℃；多年平均日照2849 h/a，平均無霜期205d。貝殼堤島地勢平坦，海拔一般在5m以下，潛水水位淺，礦化度高。貝殼砂平均厚度達1.0—2.5m，局部達3—4m，土壤類別主要是貝殼砂土類和濱海鹽土類，向海側和向陸側以濱海鹽漬土為主，成土母質由風積物和鈣質貝殼土壤化組成。植被類型以草本、灌木為主，現有樹種較少，受自然因素和平島挖砂等人類活動的影響，貝殼堤島植被主要種如杠柳、酸棗、檉柳(TamarixchinensisLour)等形成的天然次生灌木林在逐步減少，草本以砂引草、獐毛(Aeluropussinessis(Debeaux)Tzvel)、二色補血草(Limoniumbicolor(Bunge) O.Ktunze)等為主。

2 研究材料與方法

2.1 研究樣地

在貝殼堤島灌草植被集中分布的灘脊地帶，選擇生境條件一致的酸棗、杠柳天然次生灌木林，及以砂引草為主的草地為試驗樣地，并以相同地段的貝殼砂裸地作為對照。酸棗林平均樹高1.65m，平均基莖1.38cm，林分郁閉度0.85，覆蓋度90%，樹齡平均為8a生，林下草本以青蒿(ArtemisiacarvifoliaBuch.)、地膚 (Kochiascoparia(Linn.)Schrad) 、沙打旺(AstragalusadsurgensPall.)為主。杠柳林平均樹高1.36m，平均基莖1.21cm，林分郁閉度0.75，覆蓋度85%，樹齡平均為6a生，林下草本以狗尾草(Setariaviridis(Linn.)Beauv.)、鵝絨藤(Cynanchumchinense(Thunb.)Mak.)為主。砂引草為多年生草本，主要借根狀莖的延伸進行無性繁殖，水分適宜時，也能用正常的種子繁殖，株高平均為0.38m，覆蓋度65%，伴生種有青蒿和地膚。

2.2 土壤樣品的采集與測定

在每種植被類型內設置3個面積為10m×10m的樣地，在每個樣地內按S形選取5個測點，取0—20cm以及20—40cm土層的土壤樣品，把同一樣地5個土壤樣品分層混勻后，進行風干處理。烘干法測定土壤含水量，環刀浸水法測定土壤容重和孔隙度等參數[26]，并由公式計算一定土層深度內的吸持蓄水量、滯留蓄水量和飽和蓄水量[27]，本研究按0.2 m深度計算。樣品風干處理后，采用機械篩分法與比重計法測定土壤粒徑質量分布，粒徑分級標準依據中國制和參照文獻[12]。基于Turcotte[28]和Sperry[29]等土壤顆粒分形維數計算公式，利用Tyler等[3]和楊培嶺等[7]推導出的新模型公式，土壤顆粒分形維數計算公式可表示如下：

(1)

式中,D為土壤顆粒分形維數；di為兩相鄰粒級di與di+1間土粒平均直徑(mm)；dmax為最大粒級土粒平均直徑(mm)；Wi為直徑小于di的累積質量(g)；W0為土壤樣品總質量(g)。具體應用時，首先求出土壤樣品不同粒徑di的lg(di/dmax)和lg(Wi/W0)值，并將兩者進行線性擬合分析求得斜率K，則土壤分形維數為D=3-K。具體的土壤粒徑計算原理及推導過程詳見文獻[8,11]。采用SPSS13.0統計軟件中的one-way ANOVA、LSD(α=0.05， 0.01和0.001)和Peareson相關分析方法分別進行方差分析、多重比較和相關性分析。

3 結果與分析

3.1 土壤顆粒組成與分形維數

3.1.1 不同植被類型的土壤顆粒組成

由表1可知，貝殼砂土壤中，粗砂粒含量最高，在50.15%—72.30%之間，平均為61.31%；其次為細砂粒，含量在8.65%—32.65%之間，平均為19.97%；而石礫和粉粘粒含量相對較低，石礫含量在12.65%—20.13%之間，平均為17.33%；粉粘粒含量在0.25%—4.92%之間，平均僅為1.39%。表明貝殼砂生境內粗砂粒、細砂粒含量顯著高于石礫、粉粘粒，具有壤質砂土特點，屬于多礫質粗砂土的范疇[8]。不同植被類型土壤粒徑的質量分布差異極顯著(F=380.449，Sig.=0.000，Plt;0.001)，3種植被類型0—40cm貝殼砂土壤剖面中石礫含量差異極顯著(F=12.234，Sig.=0.000，Plt;0.001)，均值大小依次為草地lt;杠柳林lt;酸棗林lt;裸地，分別比裸地低20.43%、16.87%和3.98%。粗砂粒含量差異極顯著(F=11.205，Sig.=0.000，Plt;0.001)，均值大小依次為酸棗林lt;杠柳林lt;草地lt;裸地，分別比裸地低15.85%、11.78%和2.76%。細砂粒含量差異極顯著(F=8.740，Sig.=0.000，Plt;0.001)，均值大小依次為杠柳林gt;酸棗林gt;草地gt;裸地，分別比裸地高71.37%、62.10%和38.33%。粉粘粒含量差異極顯著(F=14.727，Sig.=0.000，Plt;0.001)，均值大小依次為酸棗林gt;杠柳林gt;草地gt;裸地，分別是裸地的8.26、4.08和2.18倍。表明灌木林地和草地具有減小石礫和粗砂粒含量，增加細砂粒和粉粘粒含量的作用，即植被恢復措施，具有使貝殼砂由粗粒徑向細粒徑轉變的效能。在垂直結構上，不同植被類型不同顆粒分布表現出一定的差異，占主要成分的粗砂粒均表層低于20—40cm，細砂粒除草地差異不顯著外(Pgt;0.05)，其他均表現為表層高于20—40cm。占百分比最少的粉粘粒含量均表現為表層高于20—40cm，石礫含量除杠柳林外，其他均表現為表層高于20—40cm。

表1 各植被類型土壤中不同粒徑范圍土壤顆粒質量與總質量百分比

3.1.2 土壤顆粒分形維數與土壤粒級分布的關系

3種植被類型及裸地0—40cm土壤顆粒分形維數差異極顯著(F=24.70，Sig.=0.000，Plt;0.001)，均值大小表現為酸棗林gt;杠柳林gt;草地gt;裸地(圖2)，與裸地相比，分別增加36.52%，23.67%，12.92%。在垂直結構上，土壤顆粒分形維數均表現為表層高于20—40cm土層，差異均極顯著(Plt;0.001)。不同土壤層次顆粒分形維數大小均表現為酸棗林gt;杠柳林gt;草地gt;裸地，0—20cm分形維數分別比裸地高48.39%、32.75%和15.69%；20—40cm分形維數分別比裸地高23.96%、14.07%和9.98%。表明不同的植被恢復措施對貝殼砂土壤分形維數影響較大，并且隨著貝殼砂深度的不同，其分形維數也表現出較大差異。

土壤顆粒分形維數對各個粒級土粒含量的反映程度不同，為確定分形維數與各粒級含量的關系，對分形維數與石礫、粗砂粒、細砂粒和粉粘粒的含量進行相關性分析(圖3)。土壤顆粒分形維數與粉粘粒含量呈極顯著正相關(相關系數r=0.940，Plt;0.01)，與細砂粒呈顯著正相關(r=0.771，Plt;0.05)，與粗砂粒含量呈極顯著負相關(r=-0.947，Plt;0.01)，與石礫含量的相關性不顯著(r=0.417，Pgt;0.05)。可見，貝殼砂生境下，土壤顆粒分形維數隨粉粘粒以及細砂粒含量的增加而增加，隨粗砂粒含量的增加而減少，其中，對土壤顆粒分形維數影響程度較大的是粗砂粒和粉粘粒含量，其次是細砂粒含量，石礫含量較小。

圖2 不同植被類型下的土壤顆粒分形維數Fig.2 Soil particle fractal dimension in different vegetation types

圖3 土壤顆粒分形維數與粒徑組成的相關性Fig.3 Relationships between soil fractal dimension and soil particle content

3.2 土壤顆粒分形維數與土壤容重和孔隙度

各植被類型下的土壤容重均低于裸地(表2)，差異顯著(F=17.212，sig.=0.001，Plt;0.05)，其中0—40cm土壤容重均值酸棗林、杠柳林、草地，分別比裸地低23.87%，14.51%和10.47%。總孔隙度差異顯著(F=9.607，sig.=0.005，Plt;0.05)，均值大小為酸棗林gt;杠柳林gt;草地，分別比裸地高16.96%，16.71%和1.31%，表明酸棗林的通氣透水性能較好，其次為杠柳林，草地較差。酸棗林和杠柳林的毛管孔隙度均值較高，差異不顯著(Pgt;0.05)，草地的毛管孔隙度與裸地接近，表明酸棗林和杠柳林土壤中有效水的貯存容量較大，利于植被根系對水分的有效利用，而草地在維持自身生長發育所貯存水分的潛能相對較低。非毛管孔隙度均值最大的為杠柳林，比裸地增加74.17%，涵養水源潛能相對較好；而酸棗林和草地差異不顯著(Pgt;0.05)，僅比裸地增加22.50%，12.71%。在垂直結構上，土壤容重表現為表層低于20—40cm土層，總孔隙度和毛管孔隙度則與之相反。

表2 不同植被類型的土壤容重和孔隙度特征

圖4 土壤顆粒分形維數與土壤容重和孔隙度的相關性Fig.4 Relationships between soil fractal dimension and soil density and porosity

圖5 土壤顆粒分形維數與土壤蓄水量的相關性Fig.5 Relationships between soil fractal dimension and soil water-storage capacity

由圖4可知，土壤顆粒分形維數與土壤容重呈極顯著負相關(r=-0.928，Plt;0.01)，與毛管孔隙度、總孔隙度均呈極顯著正相關(Plt;0.01)，相關系數分別為0.883，0.857，與非毛管孔隙度的相關性不顯著(r=0.350，Pgt;0.05)。可見，土壤顆粒分形維數對土壤容重和孔隙度狀況的反映程度不一樣，其中反映程度最大的是土壤容重，其次是毛管孔隙度，總孔隙度次之。表明貝殼砂生境下土壤容重越小，總孔隙度和毛管孔隙度越大，土壤顆粒分形維數越大，但土壤顆粒分形維數難以反映非毛管孔隙度的大小。

3.3 土壤顆粒分形維數與土壤蓄水性能

不同植被類型下土壤飽和蓄水量(F=9.628，sig.=0.006，Plt;0.05)、吸持蓄水量(F=7.606，sig.=0.010，Plt;0.05)差異均顯著，均值大小均表現為酸棗林gt;杠柳林gt;草地(表3)，其中飽和蓄水量分別比裸地高16.94%、16.70%和1.30%，吸持蓄水量分別比裸地高16.63%，13.49%和0.66%。滯留蓄水量差異極顯著(F=25.011，sig.=0.000，Plt;0.001)，均值大小表現為杠柳林gt;酸棗林gt;草地，分別是裸地的1.74、1.23和1.13倍，表明酸棗林供給植物有效水利用較好，杠柳林涵養水源潛能較好，草地貯存水分的能力較差。在垂直結構上，飽和蓄水量及吸持蓄水量均表現為土壤表層大于20—40cm。

表3 不同植被類型的土壤蓄水性能

土壤顆粒分形維數與飽和蓄水量(r=0.855，Plt;0.01)、吸持蓄水量(r=0.881，Plt;0.01)均呈極顯著正相關，與滯留蓄水量的相關性不明顯(r=0.340，Pgt;0.05)(圖5)。表明土壤顆粒分形維數對土壤蓄水指標的反映程度不一樣，其中反映程度最大的是飽和蓄水量，其次是吸持蓄水量。即飽和蓄水量和吸持蓄水量越大，土壤顆粒分形維數越大，但土壤顆粒分形維數對滯留蓄水量的反映程度不高。

4 討論與結論

4.1 貝殼砂土壤顆粒分形維數與土壤顆粒組成的關系

土壤顆粒組成與分形維數，除了與成土母質、土壤質地、物理化學風化過程有關外[8- 9,30]，還與土地利用方式[6,14]、不同植被類型[5- 6,12]有一定關系。在相同地段、相同生境內，植被可通過地表覆蓋、枯枝落葉層的攔截降雨、減弱風蝕、水蝕等過程，保存了地表層的細沙粒及粉粘粒含量；同時枯落物形成的腐殖質層及植物根系本身穿插、殘體及分泌物均可有效改變土壤物理結構和化學性能，從而影響土壤的顆粒組成及分形維數[5- 6,12,16]。黃河三角洲貝殼堤島灘脊地帶的粗砂粒含量最高，其次為細砂粒，石礫和粉粘粒含量較低。研究表明，屬于潮土類型的黃河三角洲灘地有林地粉粒含量最高，但棉花地和荒草地砂粒含量最高[6]；屬于棕壤類型的冀北山地森林土壤中粗粉粒和沙粒含量顯著高于細黏粒、細粉粒，混交林的細粒含量均高于純林[13]。可見，貝殼砂生境內土壤主要來源于風化的貝殼，顆粒組成相對較粗，但灌木林地及草地具有降低貝殼砂石礫、粗砂礫，增加細砂粒和粉粘粒的作用，并且0—20cm土層表現明顯，即植被措施具有較好的增加貝殼砂細顆粒含量的作用，并且灌木林地好于草地，這與其植被覆蓋度高、生物量大，根系發達有一定關系。

土壤粒徑大小及含量對土壤顆粒間的組合、孔隙大小、數量及幾何形態都起著決定作用[3,9,19]。我國不同質地類型內土壤分形維數按砂土類-壤土類-粘壤土類-粘土類四大類質地順序，依次增大，測量范圍值在1.834—2.904之間，其中壤質砂土最低，范圍值在1.834—2.641之間[8]。黃河三角洲灘地荒草地及有林地土壤顆粒分形維數在2.4657—2.6798之間[6]，山地森林土壤顆粒分形維數在2.0570—2.3739之間[13]，而貝殼砂土壤顆粒分形維數為1.3632—2.1416，明顯低于壤土、粘壤土及粘土類，表層高于20—40cm土層，3種植被類型的分形維數大小依次為酸棗林、杠柳林、草地，均值為1.7452，接近壤質砂土的測量最低值1.834[8]，可見與山地森林或黃河灘地土壤相比，貝殼堤島的土壤粗顆粒較多，分形維數偏低。但貝殼堤島的灌木林及荒草地仍提高了相同生境下的分形維數，由裸地的1.4032增至酸棗林下的1.9157，增幅為36.52%，表明灌木林下貝殼砂變細，易形成良好的土壤結構，分形維數有增大趨勢，0—20cm土層這種改善作用較為明顯，酸棗林分形維數表層最高，這與其覆蓋度最大，分解層枯枝落葉豐厚易增加表層土壤養分含量和粉粘粒含量有一定關系。同時樹木根系的生長影響土壤的物理化學以及生物學性質，加快了土壤風化的速度及腐殖質的形成，有利于細砂粒物質的固定[5,12,16]。而草地覆蓋度較低，枯落物儲量低，細小顆粒易被風蝕掉，草本植物根系較淺且不發達，土壤易粗粒化，分形維數較低。在無植被覆蓋的地方，因風蝕而引起細顆粒和營養物質被吹蝕，土壤有粗粒化變重趨勢，粗砂粒含量最多，保水性能差，植被生長困難，易惡變為嚴重的退化沙化質地[10- 11]。

科爾沁沙地農田沙漠化演變過程中[10]，黃河三角洲灘地不同土地利用方式下[6]，山地典型森林植被下[5,13]的土壤顆粒分形維數與土壤中砂粒、粗粉粒含量呈顯著負相關，而與粘粒、粉粒含量呈顯著的正相關，這與本研究結果基本一致。貝殼砂土壤顆粒分形維數與粉粘粒含量呈極顯著正相關，與細砂粒呈顯著正相關，與粗砂粒含量呈顯著負相關，但與大粒徑石礫含量相關性不顯著，與相關研究表明土壤質地由粗到細，分形維數由小到大的結論一致[3- 4,8- 9]，可見土壤顆粒分形維數在描述貝殼砂這一特殊土壤質地是可行的。土壤顆粒分形維數對各個粒級土粒含量的反映程度有一定差異，總體表現為土壤顆粒組成或團粒組成的分形維數是隨著土壤質地變細而增大，隨粗粒徑砂粒含量的增加而變小，這一規律隨土壤質地、土地利用方式或植被類型的不同略有差異，如冀北山地森林與黏粒含量相關性不顯著[13]，重慶四面山林地土壤分形維數與粉粒含量的相關性不顯著[5]，這可能與成土母質、土壤粒徑分布范圍、含量及質地均勻程度有關。

4.2 貝殼砂土壤顆粒分形維數與土壤物理參數的關系

結構良好的土壤容重為1.25—1.35 g/cm3，水氣關系協調的土壤總孔隙度為40%—50%。貝殼堤島灌木林地及草地土壤容重在1.04—1.44 g/cm3，均低于裸地，總孔隙度為44.53%—56.66%，表明植被覆蓋對貝殼砂土壤物理結構改善較好，灌木林地0—20cm總孔隙度達54.20%以上，土壤容重表層低于20—40cm土層，這與表層枯落物回歸土壤及腐殖質層的形成有一定關系，較高的覆蓋度，在一定程度上減弱了風蝕及降雨對土壤的沖刷淋蝕，且較多的殘次根系使毛管孔隙度增大，在一定程度上改善了土壤通氣狀況[5,12]。裸地總孔隙度偏低，容重偏大，可能與泥質海岸上潮時攜帶的淤泥堆積有一定關系，而有植被覆蓋的地方，除了植被本身的改良土壤效應外，向海側的酸棗、杠柳及檉柳等沖浪林帶能有效抵擋潮灘淤泥。

山地森林土壤分形維數與土壤容重呈正相關，與總孔隙度呈負相關關系[13]。一般而言，土壤分形維數越高，土壤質地越粘重，通透性越差；分形維數越低，土壤結構越松散[9]。但來源于風化未全的貝殼砂，土壤顆粒較粗，孔隙度較大，不利于貯存水分，因此，植被恢復措施下的貝殼砂生境則呈現與之相反的變化規律，土壤顆粒分形維數與土壤容重呈極顯著負相關，與毛管孔隙度、總孔隙度呈極顯著正相關，與非毛管孔隙度相關性不顯著。即貝殼砂生境下，土壤孔隙度大、通氣透水性能好，則分形維數也越高；而土壤變得密實，容重增大，則分形維數也越低；與山地森林表現規律不一致[13]，這主要與貝殼砂土壤中粗粒徑含量相對較高、分形維數相對較低有一定關系，這也表明隨著土壤顆粒組成、粒徑大小及含量的不同，土壤顆粒分形維數與土壤的疏松程度、通氣性能表現出一定的閾值效應。

4.3 貝殼砂土壤顆粒分形維數與土壤蓄水性能的關系

吸持蓄水量為毛管持水，主要用來貯存植物生理用水，滯留蓄水量為非毛管持水，多反映植被的涵養水源功能，飽和蓄水量為吸持貯存與滯留貯存的總和，可反映植被減少地表徑流和防止土壤侵蝕的功能[26,31]。黃河三角洲灘地有林地0—40cm飽和蓄水量在82.69—102.85mm，吸持蓄水量在79.07—93.34mm，滯留蓄水量在2.02—9.51mm[31]；貝殼堤島灌木林及草地0—40cm飽和蓄水量在89.06—113.32mm，吸持蓄水量在82.80—104.36mm，滯留蓄水量在4.56—8.96mm。可見，從保持水土、植物有效水利用及水源涵養的角度來看，貝殼堤島灘脊地帶好于黃河灘地，貝殼砂具有一定涵蓄水分的潛力，并且灌木林的蓄水性能好于草地。植被恢復措施對貝殼砂表層的孔隙結構、疏松程度改善較好，通氣、透水性能比較協調，因此灌木林地及草地的蓄水性能均表現為0—20cm高于20—40cm，這與黃河灘地的白蠟林、刺槐林結果類似[31]。

土壤孔隙大小分布決定著土壤持水性能，土壤顆粒分形維數在一定程度上能夠反映土壤的蓄水特性[5,13]。冀北山地森林土壤顆粒分形維數與最大持水量呈負相關，與田間持水量相關性不顯著[13]；但重慶四面山林地下的土壤顆粒分形維數與飽和含水量、毛管持水量和田間持水量均呈正相關，與非毛管孔隙度呈負相關[5]，這與本研究結果類似，貝殼砂生境下土壤顆粒分形維數與飽和蓄水量、吸持蓄水量呈極顯著正相關，而與滯留蓄水量相關性不大，表明土壤顆粒分形維數隨著成土母質、植被類型的不同，對土壤蓄水性能的反映表現出一定的差異。分形維數高的土壤中細沙粒或粘粒含量相對較高，而粘粒含量高的土壤利于土壤團聚體的形成，可改善土壤通氣、透水性能，增強土壤毛管孔隙度；同時單位土粒表面積越大，土壤對水分子的吸附力越大，土壤的持水性能易增強[5]。從土壤顆粒分形及其水分生態特征來看，貝殼堤島3種植被類型的改良土壤物理性質及蓄水保土功能表現為灌木林好于草地，其中酸棗林好于杠柳林。

[1] Tyler S W, Wheatcraft S W. Fractal scaling of soil particle size distributions: analysis and limitations. Soil Science Society of American Journal, 1992, 56(2): 362- 369.

[2] Bittelli M, Campbell G S, Flury M. Characterization of particle-size distribution in soils with a fragmentation model. Soil Science Society of America Journal, 1999, 63(4): 782- 788.

[3] Tyler S W, Wheatcraft S W. Application of fractal mathematics to soil water retention estimation. Soil Science Society of American Journal, 1989, 53(4): 987- 996.

[4] Wu C Z, Hong W. Fractal features of soil aggregate structure under different management model. Acta Pedologica Sinica, 1999, 36(2): 162- 167.

[5] Wang X, Zhang H J, Cheng J H, Zhang K, Zhang J W, Song N, Sun L. Fractal characteristics and related affecting factors of particle size distribution of different forest soil in Simian Mountains Chongqing. Journal of Soil and Water Conservation, 2011, 25(3): 154- 159.

[6] Lü S Q, Gao P, Geng G P, Zhang J, Xia J B. Characteristics of soil particles and their correlation with soil organic matter in low lands of the Yellow River Delta. Journal of Soil and Water Conservation, 2011, 25(6): 134- 138.

[7] Yang P L, Luo Y P, Shi Y C. Soil fractal character token by particle mass distribution. Chinese Science Bulletin, 1993, 38(20): 1896- 1899.

[8] Li D C, Zhang T L. Fractal features of particle size distribution of soils in China. Soil and Environmental Sciences, 2000, 9(4): 263- 265.

[9] Huang G H, Zhan W H. Fractal property of soil particle size distribution and its application. Acta Pedologica Sinica, 2002, 39(4): 490- 497.

[10] Su Y Z, Zhao H L. Fractal features of soil particle size distribution in the desertification process of the farmland in Horqin Sandy Land. Acta Ecologica Sinica, 2004, 24(1): 71- 74.

[11] Li H L, Wan L L, Dong Z, Liu Z, Wang L Y. Effects of sand barriers ofSalixPsammophilaon soil particle size and fractal dimension. Chinese Journal of Soil Science, 2012, 43(3): 540- 545.

[12] Liu X, Zhang G C, Heathman G C, Wang Y Q, Huang C H. Fractal features of soil particle-size distribution as affected by plant communities in the forested region of Mountain Yimeng, China. Geoderma, 2009, 154(1/2), 123- 130.

[13] Liu Y, Chen B, Yang X B, Zhao X M, Wang Y M. Fractal characteristics of soil particles of typical forest in north mountain of Hebei Province. Journal of Soil and Water Conservation, 2012, 26(3): 159- 163, 168.

[14] Wang D, Fu B J, Chen L D, Zhao W W, Wang Y F. Fractal analysis on soil particle size distributions under different land-use types: a case study in the loess hilly areas of the Loess Plateau, China. Acta Ecologica Sinica, 2007, 27(7): 3081- 3089.

[15] Zhang P J, Zhao Y Q. Fractal features of soil particles of riverain wetlands recovered from farmlands in Anqing, Anhui Province. Journal of Ecology and Rural Environment, 2012, 28 (2): 128- 132.

[16] Yu D M, Hu X L, Zhang G C, Liu X, Yao X Y, Hu X L. Multifractal analysis on soil particle size distribution for different vegetation types in Jiangzihe small watershed. Science of Soil and Water Conservation, 2011, 9(5): 79- 85.

[17] Martín M A, Rey J M, Taguas F J. Anentropy based heterogeneity index form as size distributions in earth science. Ecological Modelling, 2005, 182(3/4): 221- 228.

[18] Martín M A, Montero E. Laser diffraction and multifractal analysis for the characterization of dry soil volume-size distributions. Soil and Tillage Research, 2002, 64(1/2): 113- 123.

[19] Filgueira R R, Fournier L L, Cerisola C I, Gelati P, Garcia M G. Particle-size distribution in soils: a critical study of the fractal model validation. Geoderma, 2006, 134(3/4): 327- 334.

[20] Tian J Y, Xie W J, Sun J K. Current status of vulnerable ecosystem of shell islands and protection measures in Yellow River Delta. Environmental science and management, 2009, 34(8): 138- 143.

[21] Liu Z J, Zhang X L, Li P, Li P Y, Xu Y Q. Protectionstatus and management countermeasures of Binzhou shell barrier islands and wetlands system. Ocean Development and Management, 2010, 27(1): 65- 68.

[22] Zhao Y Y, Hu X M, Liu J T. Distribution characteristics of soil microorganisms and its relationship with vegetation in shell bars. Bulletin of Soil and Water Conservation, 2012, 32(2): 267- 270.

[23] Xia J B, Tian J Y, Zhang G C, Li T. Photosynthetic and physiological characteristics of three shrubs species in shell islands of Yellow River Delta. Acta Botanica Boreali-Occidentalia Sinica, 2009, 29(7): 1452- 1459.

[24] Li T, Liu Q, Tian J Y, Sun J K. Reponses on salt stress on characteristics of growth and protective enzymes ofLimoniumBicolorin shell islands of the Yellow River Delta. Bulletin of Soil and Water Conservation, 2010, 30(1): 85- 88.

[25] Liu Q, Sun J K, Tian J Y, Lv Z J, Xu S. The characteristics of microelements content and forms of seashell islands in Yellow River Delta. Journal of Soil and Water Conservation, 2009, 23(4): 204- 207, 212.

[26] Institute of Soil Science, Chinese Academy of Sciences. ed. Soil Physical and Chemical Analysis. Shanghai: Shanghai Science and Technology Press, 1978: 511- 512, 522- 524.

[27] Xia J B，Xie W J，Lu Z H, Jia Q, Dong L J. Effects of different irrigation modes with reclaimed water on soil eco-hydrological characteristics of reed land. Acta Ecologica Sinica, 2010, 30(15): 4137- 4143.

[28] Turcotte D L. Fractal fragmentation. Journal of Geography Research, 1986, 91(12): 1921- 1926.

[29] Sperry J S, Hacke U G. Desert shrub water relations with respect to soil characteristics and plant functional type. Functional Ecology, 2002, 16(3): 367- 378.

[30] Montero E. Renyi dimensions analysis of soil particle-size distributions. Ecological Modelling, 2005, 182(3/4): 305- 315.

[31] Xu J W, Li C R, Xia J B, Liu L J, Wang Y H. Soil hydrological properties of different forest types in lowlands of the Yellow River Delta. Journal of Soil and Water Conservation, 2009, 23(1): 173- 176.

參考文獻：

[4] 吳承禎, 洪偉. 不同經營模式土壤團粒結構的分形特征研究. 土壤學報, 1999, 36(2): 162- 167.

[5] 王賢, 張洪江, 程金花, 張焜, 張靜雯, 宋楠, 孫龍.重慶四面山幾種林地土壤顆粒分形特征及其影響因素. 水土保持學報, 2011, 25(3): 154- 159.

[6] 呂圣橋, 高鵬, 耿廣坡, 張杰, 夏江寶. 黃河三角洲灘地土壤顆粒分形特征及其與土壤有機質的關系. 水土保持學報, 2011, 25(6): 134- 138.

[7] 楊培嶺, 羅元培, 石元春. 用粒徑的重量分布表征的土壤分形特征. 科學通報, 1993, 38(20): 1896- 1899.

[8] 李德成, 張桃林. 中國土壤顆粒組成的分形特征研究. 土壤與環境, 2000, 9(4): 263- 265.

[9] 黃冠華, 詹衛華. 土壤顆粒的分形特征及其應用. 土壤學報, 2002, 39(4): 490- 497.

[10] 蘇永中, 趙哈林. 科爾沁沙地農田沙漠化演變中土壤顆粒分形特征. 生態學報, 2004, 24(1): 71- 74.

[11] 李紅麗, 萬玲玲, 董智, 劉振, 王麗英. 沙柳沙障對沙丘土壤顆粒粒徑及分形維數的影響. 土壤通報, 2012, 43(3): 540- 545.

[13] 劉陽, 陳波, 楊新兵, 趙心苗, 王永明. 冀北山地典型森林土壤顆粒分形特征. 水土保持學報, 2012, 26(3): 159- 163, 168.

[14] 王德, 傅伯杰, 陳利頂, 趙文武, 汪亞峰. 不同土地利用類型下土壤粒徑分形分析——以黃土丘陵溝壑區為例. 生態學報, 2007, 27(7): 3081- 3089.

[15] 張平究, 趙永強. 退耕還湖后安慶沿江濕地土壤顆粒分形特征. 生態與農村環境學報, 2012, 28 (2): 128- 132.

[16] 于東明，胡小蘭，張光燦，劉霞，姚孝友，胡續禮. 江子河小流域不同植被類型土壤粒徑的多重分形特征. 中國水土保持科學, 2011, 9(5): 79- 85.

[20] 田家怡, 謝文軍, 孫景寬. 黃河三角洲貝殼堤島脆弱生態系統破壞現狀及保護對策. 環境科學與管理, 2009, 34(8): 138- 143.

[21] 劉志杰, 張曉龍, 李萍, 李培英, 徐元芹. 濱州貝殼堤島與濕地系統保護現狀及其管理對策. 海洋開發與管理, 2010, 27(1): 65- 68.

[22] 趙艷云, 胡相明, 劉京濤. 貝殼堤地區微生物分布特征及其與植被分布的關系. 水土保持通報, 2012, 32(2): 267- 270.

[23] 夏江寶, 田家怡, 張光燦, 李田. 黃河三角洲貝殼堤島3種灌木光合生理特征研究. 西北植物學報, 2009, 29(7): 1452- 1459.

[24] 李田, 劉慶, 田家怡, 孫景寬. 黃河三角洲貝殼堤島二色補血草生長和保護酶特性對鹽脅迫的響應. 水土保持通報, 2010, 30(1): 85- 88.

[25] 劉慶, 孫景寬, 田家怡, 呂振江, 許帥. 黃河三角洲貝殼堤島貝殼沙中微量元素含量及形態特征. 水土保持學報, 2009, 23(4): 204- 207, 212.

[26] 中國科學院南京土壤研究所編. 土壤理化分析. 上海: 上海科學技術出版社, 1978: 511- 512, 522- 524.

[27] 夏江寶, 謝文軍, 陸兆華, 賈瓊, 董立杰. 再生水澆灌方式對蘆葦地土壤水文生態特性的影響. 生態學報, 2010, 30(15): 4137- 4143.

[31] 許景偉, 李傳榮, 夏江寶, 劉立杰, 王月海. 黃河三角洲灘地不同林分類型的土壤水文特性. 水土保持學報, 2009, 23(1): 173- 176.

WaterecologyandfractalcharacteristicsofsoilparticlesizedistributionofthreetypicalvegetationsinShellIsland

XIA Jiangbao1,2, ZHANG Shuyong1,3,*, WANG Rongrong1,3, ZHAO Yanyun1, SUN Jingkuan1，LIU Jingtao1, LIU Qing1

1BinzhouUniversity,ShandongProvincialKeyLaboratoryofEco-EnvironmentalScienceforYellowRiverDelta,Binzhou256603,China2ShandongProvincialEngineeringandTechnologyResearchCenterforWildPlantResourcesDevelopmentandApplicationofYellowRiverDelta,Binzhou256603,China3CollegeofForestry,ShandongAgriculturalUniversity,Taian271018,China

Little information is available about soil particle-size distribution in shell ridge ecosystems. Based on the principle of quality distribution of soil particle size and fractal theory,PeriplocasepiumBunge,Ziziphusjujubavar.spinosaHu andMesserschmidiasibiricaLinn. growing on Shell Island, located on the Yellow River Delta, were analyzed in contrast with bare soil to determine the fractal dimensions of soil particles, soil particle-size distribution and soil water physical characteristics, in order to explore the effects of different vegetation types on soil pore structure and hydrological physical characteristics, and the effect of fractal mechanisms of different vegetation types on soil particle composition in shell sand habitats. The results showed that: 1) In shell sand soil, the distribution of soil particles among three typical vegetations was as follows: coarse sand, fine sand, gravel and silt-clay. And the mean value of coarse, fine sand, gravel and silt-clay was 61.31%, 19.97%, 17.33%, 1.39%, respectively, indicating that grassland and shrub could decrease the particle content of gravel and coarse sand, and improve the content of soil fine sand and silt-clay, which performed well in 0—20cm soil layer. With the vegetation restoration measures, the mean value of fractal dimension of shell sand varied from 1.5845 to 1.9157. The sequence of fractal dimension was as follows:Z.jujubavar.spinosaHu,P.sepiumBunge and Grassland, and the surface soil layer was larger than that in 20—40cm soil layer. Fractal dimension had a very significant positive correlation with silt-clay (r=0.940，Plt;0.01), a significant positive correlation with fine sand (r=0.771，Plt;0.05), and a very significant negative correlation with coarse sand (r=-0.947，Plt;0.01), but the correlation with gravel was not obvious. 2) The mean value of soil density in 0—40cm ofZ.jujubavar.spinosaHu,P.sepiumBunge and Grassland were lower by 23.87%，14.51% and 10.47%, respectively, compared with bare land (1.48 g/cm3); the mean value of total porosity of soil increased by 16.96%,16.71%, 1.31%, respectively, compared with bare land (45.39%), and the degree of porosity in the surface soil layer was higher than that in 20—40cm soil layer. Fractal dimension had a very significant negative correlation with soil density, and a very significant positive correlation with capillary porosity and total porosity, but there was no significant correlation with non-capillary porosity. In shell sand habitats, the ventilation permeability was better with soil porosity decreasing, fractal dimension as well; but the higher soil-compacting degree, the bigger soil density, fractal dimension was lower. 3) Vegetation restoration measures could well improve pore structure of shell sand surface and the degree of porosity, as well as ventilation, pervious performance. Soil water storage capacity of grasslands and shrub land all performed better in 0—20cm than that in 20—40cm. Under shell sand habitats, fractal dimension had a very significant positive correlation with total water-storage capacity and capillary water-storage capacity, but the correlation with non-capillary water-storage capacity was not obvious. Judging from the fractal characteristics and soil hydrological physical characteristics, we concluded that the sequence of improving the capacity of soil physical properties and soil and water conservation of three typical vegetations in Shell Island was that shrub,was better than grassland, among which,Z.jujubavar.spinosaHu was better thanP.sepiumBunge, and soil layer in 0—20cm than that in 20—40cm.

fractal dimension; soil particle-size distribution; soil water; vegetation type; Shell Island

國家自然科學基金資助項目(31100468); 山東省高校科研發展計劃資助項目(J13LC03); 國家科技支撐資助項目(2009BADB2B05); 山東省黃河三角洲生態環境重點實驗室開放基金資助項目(2012KFJJ04)

2012- 07- 14;

2012- 11- 19

*通訊作者Corresponding author.E-mail: zhsyong@126.com

10.5846/stxb201207140994

夏江寶,張淑勇,王榮榮,趙艷云,孫景寬,劉京濤,劉慶.貝殼堤島3種植被類型的土壤顆粒分形及水分生態特征.生態學報,2013,33(21):7013- 7022.

Xia J B, Zhang S Y, Wang R R, Zhao Y Y, Sun J K, Liu J T, Liu Q.Water ecology and fractal characteristics of soil particle size distribution of three typical vegetations in Shell Island.Acta Ecologica Sinica,2013,33(21):7013- 7022.