馬劍平，王 飛,2*，郭樹江,3，紀(jì)永福,3，張瑩花，張?jiān)Ｄ?2，張衛(wèi)星，宋達(dá)成

(1.甘肅省治沙研究所，甘肅 蘭州 730070；2.甘肅省荒漠化與風(fēng)沙災(zāi)害防治國(guó)家重點(diǎn)實(shí)驗(yàn)室(培育基地)，甘肅 武威 733000；3.甘肅民勤荒漠草地生態(tài)系統(tǒng)國(guó)家野外科學(xué)觀測(cè)研究站，甘肅 民勤 733300)

土壤是植物生長(zhǎng)發(fā)育的基礎(chǔ)物質(zhì)，土壤粒徑分布是重要的土壤物理參數(shù)，能夠反映土壤風(fēng)化成土、退化過程及侵蝕程度[1]，與土壤水分[2]、結(jié)構(gòu)[3-4]、肥力狀況[5-6]關(guān)系密切，進(jìn)一步影響著土壤水分、養(yǎng)分截留與轉(zhuǎn)運(yùn)、植被生產(chǎn)力和生態(tài)恢復(fù)過程[7-8]。土壤是由不同大小、形狀的固體顆粒物構(gòu)成的多孔介質(zhì)，具有分形特征[9]。B.B.Mandelbort et al[10]提出分形理論，隨著分形理論的發(fā)展，該理論在土壤學(xué)研究中得到廣泛應(yīng)用。相關(guān)研究表明[11-13]，分形維數(shù)能夠反映土壤物理性質(zhì)，海拔、植物群落、土壤pH、水分等因素與土壤粒徑分布有關(guān)。

灌叢沙堆的發(fā)育過程是植被與荒漠環(huán)境相互作用的過程[14]，是我國(guó)干旱荒漠區(qū)特有地貌類型，在民勤綠洲-荒漠過渡帶分布廣泛。白刺(Nitrariatangutorum)是超旱生植物，耐風(fēng)沙、抗鹽堿，具有超強(qiáng)的適應(yīng)能力，由于其易分枝、耐沙埋、枝條沙埋后易著生不定根等諸多特性，因此常常形成大小不等的白刺灌叢沙堆。石羊河下游民勤綠洲外圍形成了我國(guó)干旱區(qū)獨(dú)具特色的以白刺灌叢沙堆與丘間地相間分布的景觀格局，對(duì)維護(hù)該區(qū)的生態(tài)平衡發(fā)揮著重要的作用。隨著白刺灌叢沙堆的演替發(fā)展，土壤結(jié)構(gòu)會(huì)有較大差異，但關(guān)于不同演替階段白刺灌叢沙堆土壤粒徑分布及其分形維數(shù)特征與土壤理化性質(zhì)的相互關(guān)系的研究較少。鑒于此，以民勤白刺灌叢沙堆為對(duì)象，利用野外調(diào)查與室內(nèi)試驗(yàn)相結(jié)合的方法，開展白刺灌叢沙堆演替過程中(發(fā)育階段、穩(wěn)定階段、衰退階段、嚴(yán)重衰退階段)土壤粒度組成、土壤分形維數(shù)、土壤理化性質(zhì)的研究，探析白刺灌叢沙堆演替過程中土壤粒徑分布及分形特征與土壤理化性質(zhì)的相互關(guān)系，以期為干旱荒漠區(qū)植被恢復(fù)可持續(xù)性提供理論和數(shù)據(jù)參考。

1 材料與方法

1.1 研究區(qū)概況

試驗(yàn)研究區(qū)設(shè)在民勤綠洲-荒漠過渡帶，位于巴丹吉林沙漠東南緣(38°34′-39°38′N，102°53′-102°58′E)，屬于典型的干旱荒漠地區(qū)，海拔高度1 376～1 383 m。該區(qū)屬于溫帶大陸性干旱氣候，降水量少，多年平均降水量115.2 mm，集中于7-9月；多年平均蒸發(fā)量2 419.6 mm，是降水量的21倍；多年平均氣溫為7.7 ℃；熱量資源充足，日照時(shí)間長(zhǎng)，年平均日照時(shí)數(shù)達(dá)2 832.1 h；風(fēng)大沙多，年平均風(fēng)速2.5 m·s-1；土壤類型以風(fēng)沙土為主；現(xiàn)有的植被主要包括天然和人工2種類型，主要的植物有梭梭(Haloxylonammodendron)、白刺、沙拐棗(Calligonummongolicum)、鹽生草(Halogetonglomeratus)、豬毛菜(Salsolacollina)、沙米(Agriophyllumsquarrosum)等。

1.2 試驗(yàn)設(shè)計(jì)與方法

選擇白刺灌叢沙堆集中分布地帶，根據(jù)野外實(shí)地調(diào)查和已有的劃分標(biāo)準(zhǔn)[15]，結(jié)合白刺灌叢沙堆形態(tài)、土壤及植被狀況，選擇4個(gè)不同演替階段白刺灌叢沙堆為研究對(duì)象，主要包括發(fā)育階段、穩(wěn)定階段、衰退階段、嚴(yán)重衰退階段，每個(gè)演替階段沙堆的大小、植被生長(zhǎng)狀況、風(fēng)積狀況相對(duì)一致，以減小試驗(yàn)誤差，每個(gè)演替階段3次重復(fù)。對(duì)樣地內(nèi)白刺灌叢沙堆高度、長(zhǎng)度以及植被、土壤狀況進(jìn)行調(diào)查(表1)。

表1 白刺灌叢沙堆形態(tài)、植被及土壤特征

1.2.1 土壤樣品采集 于2019年夏季在不同演化階段白刺灌叢沙堆頂部的中心位置采集土樣，將土層共劃分為4個(gè)層次：0～10、10～20、20～40、40～60 cm，分別采集2份土壤樣品，一份用于土壤含水量測(cè)定；另一份裝入自封袋，帶回實(shí)驗(yàn)室，測(cè)定土壤粒度及養(yǎng)分含量。

1.2.2 指標(biāo)測(cè)定 采集的土壤在室內(nèi)自然風(fēng)干后分成2部分：一部分用于土壤粒度測(cè)定；另一部分用于土壤理化性質(zhì)的測(cè)定。土壤粒度用馬爾文激光粒度儀Malvern Mastersizer 2000測(cè)定，各粒級(jí)含量以體積百分含量表示。試驗(yàn)結(jié)果以美國(guó)制土壤粒徑分級(jí)標(biāo)準(zhǔn)輸出：黏粒(0～2 μm)、粉粒(2～50 μm)、極細(xì)砂(50～100 μm)、細(xì)砂(100～250 μm)、中砂(250～500 μm)、粗砂(500～1 000 μm)和極粗砂(1 000～2 000 μm)。分形維數(shù)采用土壤粒度體積含量數(shù)據(jù)進(jìn)行計(jì)算，即獲取土壤粒度體積分形維數(shù)D[16]：

(1)

式中：r是粒徑；Ri表征粒徑劃分中第i級(jí)粒徑；V(r

共測(cè)定了6個(gè)土壤因子。具體測(cè)定方法如下：采用烘干法測(cè)定土壤含水量(%)；采用重鉻酸鉀氧化-外加熱法測(cè)定[18]土壤有機(jī)質(zhì)含量(%)；采用半微量凱氏法測(cè)定[18]土壤全N含量(%)；采用碳酸氫鈉提取-鉬銻抗比色法測(cè)定[18]土壤速效P含量(mg·100 g-1)；采用電導(dǎo)儀法測(cè)定[19]土壤電導(dǎo)率μS；采用電極法測(cè)定[19]土壤pH。

1.3 數(shù)據(jù)分析

用SPSS19.0軟件對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)與分析。用單因素(One-way ANOVA)和Duncan法進(jìn)行方差分析與多重比較，回歸分析進(jìn)行線性擬合，用Pearson法對(duì)各指標(biāo)間的相關(guān)性進(jìn)行分析。

2 結(jié)果與分析

2.1 白刺灌叢沙堆演替過程中土壤粒度特征

由表2可知，民勤干旱荒漠區(qū)發(fā)育階段白刺灌叢沙堆土壤以細(xì)砂含量最高，為54.79%，其次為中砂，為19.30%；其他3種演替階段白刺灌叢沙堆土壤以細(xì)砂和極細(xì)砂為主，分別為46.71%～64.06%和17.70%～36.29%。經(jīng)方差分析，除極粗砂外，不同演替階段白刺灌叢沙堆土壤各粒級(jí)含量差異極顯著(P<0.01)。黏粒、粉粒、極細(xì)砂含量在0～60 cm土層衰退階段>嚴(yán)重衰退階段>穩(wěn)定階段>發(fā)育階段，發(fā)育階段黏粒、極細(xì)砂含量顯著(P<0.05)低于其他演替階段；細(xì)砂含量在0～60 cm土層穩(wěn)定階段>嚴(yán)重衰退階段>發(fā)育階段>衰退階段，差異極顯著(P<0.01)，其中穩(wěn)定階段細(xì)砂含量顯著(P<0.05)高于發(fā)育與衰退階段；中砂含量在0～60 cm土層發(fā)育階段>穩(wěn)定階段>衰退階段>嚴(yán)重衰退階段，發(fā)育階段中砂含量顯著(P<0.05)高于其他演替階段；粗砂含量在0～60 cm土層發(fā)育階段>衰退階段>嚴(yán)重衰退階段>穩(wěn)定階段，差異顯著(P<0.05)，其中發(fā)育階段粗砂含量顯著(P<0.05)高于嚴(yán)重衰退階段與穩(wěn)定階段；極粗砂含量在0～60 cm土層不同演替階段差異不顯著(P>0.05)。綜上所述，隨著白刺灌叢沙堆不斷演替發(fā)展，黏粒、粉粒、極細(xì)砂含量逐漸增多，到嚴(yán)重衰退階段降低；細(xì)砂含量從發(fā)育階段到穩(wěn)定階段呈升高趨勢(shì)，穩(wěn)定階段達(dá)到最大，隨后呈下降趨勢(shì)；中砂、粗砂、極粗砂含量整體呈降低趨勢(shì)，發(fā)育階段含量最高。

表2 白刺灌叢沙堆土壤粒度組成

發(fā)育階段、穩(wěn)定階段白刺灌叢沙堆土壤各粒級(jí)在不同土層之間差異不顯著(P>0.05)；衰退階段0～10 cm土層黏粒、極細(xì)砂含量顯著(P<0.05)高于40～60 cm土層，而40～60 cm土層中砂、粗砂含量顯著(P<0.05)高于0～10 cm土層；嚴(yán)重衰退階段0～10 cm土層黏粒、粉粒含量顯著(P<0.05)高于其他土層，而40～60 cm土層細(xì)砂含量顯著(P<0.05)高于0～40 cm土層。

4個(gè)演替階段土壤分形維數(shù)變化范圍為1.665～2.242，隨著白刺灌叢沙堆的演替，分形維數(shù)呈先升高后降低趨勢(shì)，衰退階段達(dá)到最大，發(fā)育階段白刺灌叢沙堆0～60 cm土壤分形維數(shù)顯著(P<0.05)低于其他3個(gè)演替階段，其他3個(gè)演替階段差異不顯著(P>0.05)。同時(shí)，4個(gè)演替階段中表層0～10 cm分形維數(shù)高于其他土層。

2.2 土壤粒度分形維數(shù)與粒徑分布的關(guān)系

由表3可知，除細(xì)砂、粗砂、極粗砂以外，土壤分形維數(shù)與其他粒級(jí)顆粒百分含量間呈一定的線性關(guān)系。分形維數(shù)與黏粒、粉粒、極細(xì)砂、中砂的擬合程度均較好，決定系數(shù)介于0.573～0.783，且分形維數(shù)與黏粒、粉粒、極細(xì)砂呈極顯著正相關(guān)(P<0.01)，表明土壤中細(xì)顆粒物質(zhì)越多，分形維數(shù)值越大；與中砂呈極顯著負(fù)相關(guān)(P<0.01)，即土壤中砂粒含量越高，分形維數(shù)值越小。

表3 分形維數(shù)與各參數(shù)間線性擬合關(guān)系

2.3 不同演替階段白刺灌叢沙堆土壤理化性質(zhì)

由表4可知，不同演替階段白刺灌叢沙堆0～60 cm土壤含水量隨土層增加整體呈增加趨勢(shì)。經(jīng)方差分析，衰退階段白刺灌叢沙堆土壤含水量顯著(P<0.05)高于發(fā)育、嚴(yán)重衰退階段，與穩(wěn)定階段差異不顯著(P>0.05)。不同演替階段白刺灌叢沙堆土壤有機(jī)質(zhì)、全氮、速效磷含量均表現(xiàn)出隨土層增加而降低的趨勢(shì)，其表層(0～10 cm)含量最大，表聚性明顯。經(jīng)方差分析，4個(gè)演替階段白刺灌叢沙堆0～60 cm土壤有機(jī)質(zhì)、全氮、速效磷含量差異顯著(P<0.05)，變化規(guī)律為衰退階段>嚴(yán)重衰退階段>穩(wěn)定階段>發(fā)育階段，其中衰退階段顯著(P<0.05)高于穩(wěn)定、發(fā)育階段，與嚴(yán)重衰退階段差異不顯著(P>0.05)，而穩(wěn)定與發(fā)育階段差異不顯著(P>0.05)。4個(gè)演替階段白刺灌叢沙堆0～60 cm土壤電導(dǎo)率差異不顯著(P>0.05)，變化規(guī)律為衰退階段>嚴(yán)重衰退階段>穩(wěn)定階段>發(fā)育階段，其中衰退階段顯著(P<0.05)高于發(fā)育階段。4個(gè)演替階段白刺灌叢沙堆0～60 cm土壤pH差異顯著(P<0.05)，變化規(guī)律與其他指標(biāo)相反，發(fā)育階段顯著(P<0.05)高于衰退階段。4個(gè)演替階段白刺灌叢沙堆不同土層有機(jī)質(zhì)、速效磷、pH差異不顯著(P>0.05)；發(fā)育、穩(wěn)定、衰退階段白刺灌叢沙堆不同土層全氮含量、電導(dǎo)率差異不顯著(P>0.05)，而嚴(yán)重衰退階段表層(0～10 cm)顯著(P<0.05)高于其他土層。

表4 白刺灌叢沙堆土壤理化性質(zhì)

2.4 土壤理化性質(zhì)、土壤粒度組成及分形維數(shù)間的相關(guān)關(guān)系

由表5可知，黏粒、粉粒及極細(xì)砂含量與有機(jī)質(zhì)、全N、速效P含量呈極顯著正相關(guān)(P<0.01)，且與含水量、電導(dǎo)率、pH值的相關(guān)性不顯著(P>0.05)；細(xì)砂含量與有機(jī)質(zhì)、全N含量呈極顯著負(fù)相關(guān)(P<0.01)，與pH值呈極顯著正相關(guān)(P<0.01)；中砂含量與有機(jī)質(zhì)、全N含量呈極顯著負(fù)相關(guān)(P<0.01)，與速效P含量呈顯著負(fù)相關(guān)(P<0.05)；粗砂含量與全N含量呈顯著負(fù)相關(guān)(P<0.05)；極粗砂含量與各指標(biāo)之間相關(guān)性不顯著(P>0.05)。分形維數(shù)與有機(jī)質(zhì)、全N含量呈極顯著正相關(guān)(P<0.01)，與速效P含量呈顯著正相關(guān)(P<0.05)。綜合分析表明，白刺灌叢沙堆土壤粒度的形成受到土壤有機(jī)質(zhì)、全N、速效P及pH的影響較大(P<0.05)。

表5 土壤理化性質(zhì)與土壤粒度組成、分形維數(shù)間的相關(guān)關(guān)系

3 結(jié)論與討論

3.1 結(jié)論

民勤白刺灌叢發(fā)育階段沙堆土壤以細(xì)砂(54.79%)含量最高，其次為中砂(19.30%)；其他3種演替階段白刺灌叢沙堆土壤以細(xì)砂和極細(xì)砂為主，分別為46.71%～64.06%和17.70%～36.29%，隨著白刺灌叢沙堆的演替，土壤細(xì)顆粒呈先增加后降低的趨勢(shì)，衰退階段土壤黏粒、粉粒、極細(xì)砂含量最高。分形維數(shù)與黏粒、粉粒、極細(xì)砂呈極顯著正相關(guān)(P<0.01)，與中砂呈極顯著負(fù)相關(guān)(P<0.01)，表明黏粉粒、極細(xì)砂含量越多，土壤分形維數(shù)越高。土壤分形維數(shù)與土壤養(yǎng)分之間存在顯著相關(guān)性，因此可以用分形維數(shù)來反映土壤質(zhì)地及肥力狀況。白刺灌叢沙堆充分發(fā)揮了防風(fēng)固沙作用以及對(duì)土壤的保育功能，有利于改善土壤結(jié)構(gòu)和肥力，為荒漠生態(tài)系統(tǒng)的恢復(fù)和重建提供物質(zhì)保障。

3.2 討論

3.2.1 白刺灌叢沙堆演替過程中土壤粒度組成與分形維數(shù) 土壤顆粒組成受母質(zhì)特征及環(huán)境變化的影響[20]，在一定程度上決定土壤的基本性狀。民勤荒漠-綠洲過渡帶生態(tài)環(huán)境脆弱，近年來，由于土壤水分收支不平衡導(dǎo)致白刺沙堆退化[21]。本研究中，白刺灌叢沙堆在演替過程中，土壤粒度組成仍以細(xì)砂為主，除極粗砂外，不同演替階段白刺灌叢沙堆土壤各粒級(jí)含量差異極顯著(P<0.01)。其中從發(fā)育階段到穩(wěn)定階段粒度組成變化最大，除粉粒與極粗砂外，其他土壤各粒級(jí)含量存在顯著差異(P<0.05)，表現(xiàn)為：白刺灌叢沙堆穩(wěn)定后，黏粒、極細(xì)砂、細(xì)砂含量增加，中砂與粗砂含量降低。到衰退階段，黏粒、粉粒、極細(xì)砂含量增加，細(xì)砂與中砂含量降低，但粗砂與極粗砂含量升高，與穩(wěn)定階段相比，黏粒、粉粒增加顯著，細(xì)砂顯著降低(P<0.05)。嚴(yán)重衰退階段，白刺灌叢沙堆活化，黏粒、粉粒、極細(xì)砂含量降低，砂粒含量增加。本研究中，白刺灌叢沙堆處衰退階段時(shí)，黏粒、粉粒含量相對(duì)較高，砂粒含量降低，主要是因?yàn)樵撾A段處于白刺灌叢沙堆衰退初期，與穩(wěn)定階段相比，白刺存在一定的枯枝率，但是覆蓋度仍較高，以及沙堆表層結(jié)皮較厚、堅(jiān)固，有效地減弱了風(fēng)力對(duì)灌叢沙堆土壤黏粒、粉粒的吹蝕。從整個(gè)土層垂直剖面來看，發(fā)育、穩(wěn)定階段白刺灌叢沙堆土壤黏粒、粉粒含量表現(xiàn)為“高-低-高”，衰退階段與嚴(yán)重衰退階段表現(xiàn)為“高-低”。

土壤粒度分形維數(shù)能夠反映土壤顆粒分布及質(zhì)地的均勻程度[22-23]，一般來說土壤分形維數(shù)越大，土壤抗侵蝕能力越強(qiáng)[24]。本研究中，4個(gè)演替階段土壤分形維數(shù)變化范圍為1.665～2.242，發(fā)育階段白刺灌叢沙堆0～60 cm土壤分形維數(shù)顯著(P<0.05)低于其他3個(gè)演替階段，其他3個(gè)演替階段差異不顯著(P>0.05)，說明白刺灌叢沙堆對(duì)土壤分形維數(shù)的影響較大，發(fā)育階段白刺灌叢沙堆植被覆蓋度相對(duì)較低，風(fēng)蝕嚴(yán)重，土壤黏粉粒被風(fēng)蝕。結(jié)果表明，分形維數(shù)與黏粒、粉粒、極細(xì)砂呈極顯著正相關(guān)(P<0.01)，與中砂呈極顯著負(fù)相關(guān)(P<0.01)，表明土壤粒度分形維數(shù)較好地反映了白刺灌叢沙堆的顆粒粗細(xì)狀況，且分形維數(shù)隨著黏粒、粉粒與極細(xì)砂含量的增加而增大，隨著中砂含量的增加而減小。因此，土壤分形維數(shù)在一定程度上能夠反映土壤粒徑的分布狀況和土壤質(zhì)地的均勻程度，并作為判斷土壤質(zhì)地差異的重要指標(biāo)。

3.2.2 分形維數(shù)與土壤顆粒組成與理化性質(zhì)的關(guān)系 土壤顆粒組成影響土壤養(yǎng)分的供應(yīng)與植物的生長(zhǎng)發(fā)育。土壤養(yǎng)分含量主要與土壤細(xì)顆粒含量有關(guān)，即土壤細(xì)顆粒物質(zhì)的含量更能夠反映土壤質(zhì)量水平[25]。本研究發(fā)現(xiàn)，黏粒、粉粒及極細(xì)砂含量與有機(jī)質(zhì)、全N、速效P含量呈極顯著正相關(guān)(P<0.01)；細(xì)砂含量與有機(jī)質(zhì)、全N含量呈極顯著負(fù)相關(guān)(P<0.01)，與pH值呈極顯著正相關(guān)(P<0.01)；中砂含量與有機(jī)質(zhì)、全N含量呈極顯著負(fù)相關(guān)(P<0.01)，與速效P含量呈顯著負(fù)相關(guān)(P<0.05)；粗砂含量與全N含量呈顯著負(fù)相關(guān)(P<0.05)；極粗砂含量與各指標(biāo)之間相關(guān)性不顯著(P>0.05)。這表明，土壤黏粒、粉粒及極細(xì)砂含量與土壤養(yǎng)分含量密切相關(guān)，這與前人的研究結(jié)果相一致[25]。本研究中，分形維數(shù)與有機(jī)質(zhì)、全N含量呈極顯著正相關(guān)(P<0.01)，與速效P含量呈顯著正相關(guān)(P<0.05)，同時(shí)分形維數(shù)與黏粒、粉粒、極細(xì)砂呈極顯著正相關(guān)(P<0.01)，表明土壤顆粒組成特征、分形維數(shù)和養(yǎng)分含量三者之間密切相關(guān)，土壤中細(xì)顆粒含量增多有利于土壤養(yǎng)分的積累，在一定程度上土壤分形維數(shù)可以作為土壤質(zhì)量演變的評(píng)價(jià)指標(biāo)。土壤分形維數(shù)與含水量、pH間均無顯著相關(guān)性(P>0.05)，這與羅雅曦等[26]的研究結(jié)果相一致，說明土壤顆粒大小不是土壤含水量與pH值變化的主要影響因素。土壤分形維數(shù)與電導(dǎo)率無顯著相關(guān)性(P>0.05)，這與蔣嘉瑜[27]的研究結(jié)果一致，可能與研究區(qū)惡劣的自然環(huán)境狀況有關(guān)，使鹽堿程度受土壤顆粒組成影響比較小。