劉進輝，王雪芹，馬　洋1中國科學院新疆生態(tài)與地理研究所，烏魯木齊　8300112新疆策勒荒漠草地生態(tài)系統(tǒng)國家野外科學觀測研究站，策勒　8483003中國科學院大學，北京　100049

?

沙漠綠洲過渡帶檉柳灌叢沙堆-丘間地系統(tǒng)土壤養(yǎng)分空間異質性

劉進輝1，2，3，王雪芹1，*，馬洋1，2，3
1中國科學院新疆生態(tài)與地理研究所，烏魯木齊830011
2新疆策勒荒漠草地生態(tài)系統(tǒng)國家野外科學觀測研究站，策勒848300
3中國科學院大學，北京100049

摘要:以塔克拉瑪干沙漠南緣策勒沙漠綠洲過渡帶為研究區(qū)，選取植被蓋度依次為30%、15%—20%、10%和＜5%的4個典型樣地，對各樣地的檉柳灌叢沙堆-丘間地系統(tǒng)典型部位0—10 cm和10—20 cm土壤進行系統(tǒng)采樣，分析土壤枯落物、有機質、全效N、P、K和速效N、P、K含量，從風沙地貌的角度研究過渡帶沙漠化進程中土壤養(yǎng)分分布特征。結果表明:(1)隨著植被總蓋度的降低，灌叢沙堆表層0—10 cm土壤的枯落物、有機質、全N、全K、速效K含量呈降低趨勢，平均降幅為69.3%、37.0%、35.3%、8.3%和24.5%。全P含量沒有明顯的變化，速效N和速效P含量呈先減小后增大的趨勢;(2)從灌叢下→沙堆邊緣→丘間地→風影區(qū)，土壤枯落物、有機質、全N、全P、速效N、速效P、速效K含量表現(xiàn)出先減小后增大的趨勢，最大值位于灌叢下部位，最小值位于丘間地，但在植被蓋度＜5%的樣地，沙堆-丘間地系統(tǒng)各部位之間養(yǎng)分含量無顯著性差異;(3)與0—10 cm土層相比，10—20 cm土層的枯落物、有機質、全N、全P和速效P含量顯著降低，平均降幅為40.0%、27.0%、25.0%、4.5%和2.0%，全K、速效N和速效K含量明顯增加，平均增幅為4.8%、103.3%和12.1%。隨著植被蓋度的降低，10—20 cm土層的養(yǎng)分含量具有與0—10 cm土層相對一致的變化趨勢;(4)灌叢的生物反饋作用使得灌叢下部位具有明顯的養(yǎng)分富集效應，但隨著植被總蓋度的降低和風沙活動不斷加劇，非生物因素逐漸占據(jù)主導，灌叢的養(yǎng)分富集效應趨于減弱和消失。

關鍵詞:沙漠綠洲過渡帶;檉柳灌叢沙堆;枯落物;土壤養(yǎng)分

劉進輝，王雪芹，馬洋.沙漠綠洲過渡帶檉柳灌叢沙堆-丘間地系統(tǒng)土壤養(yǎng)分空間異質性.生態(tài)學報，2016，36(4):979-990.

Liu J H，Wang X Q，Ma Y.Spatial variation of soil nutrients of Tamarix ramosissima nebkhas and interdune areas in a desert-oasis ecotone.Acta Ecologica Sinica，2016，36(4):979-990.

灌叢沙堆是風沙流遇到灌叢阻攔，沙物質在灌叢及其周圍堆積而成的生物地貌類型［1］，廣泛分布于世界干旱、半干旱及半濕潤沙地和部分沙質海岸帶［2-5］。它的形成演化是風成沙與灌叢相互作用的結果［6］，如在塔里木盆地，檉柳灌叢沙堆受植物和風沙活動的共同作用［7］，形成沙物質與枯落物互層的交錯層理［8］。因此凡是能提供沙包形成基本條件的地區(qū)都有它的分布［9］，也是沙地和沙漠的重要組成部分［10］。

Klemmedson和Barth［11］在對荒漠生態(tài)系統(tǒng)中植物生物量與養(yǎng)分分配關系研究時，首次提出了“肥島”的概念。陳廣生等［12］將其定義為干旱、半干旱區(qū)灌木冠幅下限制性土壤資源的顯著聚集現(xiàn)象，其中的資源是指土壤水分、養(yǎng)分、微生物、動物及其由灌木或喬木等帶來的非生物環(huán)境等的總和。肥島現(xiàn)象在很多區(qū)域及其不同種類的灌叢下被發(fā)現(xiàn)［13-15］，普遍的研究都認為灌叢沙堆相比周圍空地具有較強的資源富集性［13，16-17］。在干旱區(qū)，水分及/或氮是最常見的限制植物生長的非生物因素［18］，養(yǎng)分資源的空間分布對干旱環(huán)境中植被格局的確定具有重要的作用［19］，養(yǎng)分有效性的提高對群落動態(tài)將產(chǎn)生深遠的影響［20］。Gibbens［21］對奇瓦瓦沙漠植物根系進行了研究，證實了該沙漠植物根系對上層土壤水分與養(yǎng)分存在激烈的競爭。雖然在干旱區(qū)表層土壤的含水量極低，但當水分較充分時，土壤養(yǎng)分會快速地變成限制沙漠生物過程的主要因素［22］。有研究發(fā)現(xiàn)，沙漠化逆轉與較高的土壤養(yǎng)分水平相關［23］。肥島是沙漠生態(tài)系統(tǒng)結構與功能的一個重要組成部分，因此有關肥島土壤養(yǎng)分特點的研究對于了解區(qū)域環(huán)境變化具有非常重要的意義［24］。尹傳華［25］研究發(fā)現(xiàn)，當生態(tài)系統(tǒng)退化為沙漠或半沙漠時，肥島效應趨于削弱。Okin等［26］的研究則認為，人為干擾導致干旱的灌木地土壤結皮和植被破壞，進而引發(fā)肥島的解體，而肥島的解體會對相鄰地區(qū)形成間接干擾。

塔克拉瑪干沙漠南緣氣候極端干旱，風沙活動頻繁，不同植被條件下灌叢沙堆具有不同的蝕積特點。雖然對綠洲-沙漠過渡帶灌木引起的資源島現(xiàn)象已有報道［25，27］，但基于沙漠化進程中不同植被覆蓋條件下灌叢沙堆-丘間地系統(tǒng)土壤養(yǎng)分分異規(guī)律尚無系統(tǒng)的觀測和研究。鑒于此，本項研究以策勒綠洲沙漠過渡帶為研究區(qū)，在對不同植被覆蓋條件下檉柳灌叢沙堆-丘間地系統(tǒng)蝕積特征有所認識的前提下，重點研究土壤養(yǎng)分的空間異質性變化，揭示過渡帶沙漠化進程中檉柳灌叢沙堆-丘間地系統(tǒng)土壤養(yǎng)分分異規(guī)律，從風沙地貌的角度探討灌叢沙堆蝕積過程對土壤養(yǎng)分分布的影響。研究結果有助于深入理解綠洲-沙漠過渡帶土地退化成因和過程，可為綠洲外圍灌叢沙堆的科學保育和生態(tài)安全提供理論依據(jù)。

1　區(qū)域概況

研究區(qū)位于塔克拉瑪干沙漠南緣策勒綠洲西部過渡帶，介于80°03'—82°10'E，35°17'—39°30'N之間。該區(qū)域氣候極端干旱，年均降水量35 mm，蒸發(fā)量2595.3 mm，干燥度20.8，屬暖溫帶大陸性荒漠氣候。多年平均氣溫11.9℃，1月平均氣溫－5.8℃，7月平均氣溫25.1℃，極端最高氣溫41.9℃，極端最低氣溫－23.9℃。這里風沙活動頻繁，以西風為主，其頻率占62.43%—76.25%，年平均風速1.9 m/s，最大可達12.1 m/s。春季多大風，年平均3—9次，8級以上大風年平均40 d。土壤以風沙土和棕漠土為主，主要由極細沙和粉沙組成，土壤質地輕，含水量低，易于形成風沙流。天然植被類型單調，群落結構簡單，疏葉駱駝刺(Alhagi sparsifolia)、多枝檉柳(Tamarix ramosissima)和花花柴(Karelinia caspic)為主要建群種，部分可形成大小不等的植物沙包［28］。該區(qū)域沙漠化問題在綠洲西部有集中體現(xiàn)，受不同程度人類活動的影響，形成植物組合與蓋度呈現(xiàn)明顯差異的4種植被類型(表1)，分別與不同程度的沙漠化土地相對應，其間均有檉柳灌叢沙堆分布。楊帆等人［29］通過詳細的灌叢沙堆蝕積分布研究發(fā)現(xiàn)，不同植被覆蓋條件下檉柳灌叢沙堆-丘間地系統(tǒng)地表蝕積過程也出現(xiàn)明顯差異:在植被覆蓋較好的非沙漠化區(qū)域，地表以微弱積沙過程為主，主要是由于良好的植被攔截了相當一部分遠程輸送的沙塵;在植被條件居中的輕度和中度沙漠化區(qū)域，則出現(xiàn)地表物質的原地再分配，沙堆迎風坡和沙堆兩側局部發(fā)生風蝕，其它部位以積沙為主;在植被遭受嚴重破壞的重度沙漠化區(qū)域，除沙堆中上部外的大部分地表均以嚴重風蝕為顯著特點。不同程度沙漠化地表植被狀況和蝕積特征描述見表1。

表1　不同植被類型中檉柳灌叢沙堆地表蝕積特征描述Table1　The description of Tamarix ramosissima nebkha surface erosion feature in different vegetation types

2　研究方法

2.1實驗設計

基于大范圍野外調查，在策勒綠洲西部過渡帶確定4種不同植被類型的典型樣地，各樣地分別選擇一個大小適中、相對獨立的檉柳灌叢沙堆。其中，沙堆A位于植被蓋度為30%左右的疏葉駱駝刺群落中，沙堆B位于植被蓋度為15%—20%左右的疏葉駱駝刺-多枝檉柳群落中，沙堆C位于植被蓋度為10%的多枝檉柳-疏葉駱駝刺-花花柴群落中，沙堆D位于植被蓋度＜5%的多枝檉柳群落中。根據(jù)主風向、沙堆部位及蝕積分布狀況，將各灌叢沙堆-丘間地劃分為4個部位，依次為沙堆灌叢下(Under shrub)、沙堆邊緣區(qū)(Margin of nebkha)、沙堆風影區(qū)(Leeward shadow)和丘間地(Interdune)。其中沙堆灌叢下包括1、2、3、4、5號共5個采樣點，沙堆邊緣區(qū)包括6、8、9號共3個采樣點，沙堆風影區(qū)包括7、11、14號共3個采樣點，丘間地包括10、12、13、15、16、17、18號共7個采樣點(圖1)。

圖1　土壤采樣點位置Fig.1The location of soil sampling灌叢沙堆:是指叢下積土的范圍;灌叢下:是指灌叢基部分布(或影響)范圍

2.2土樣采集

各樣點分0—10 cm和10—20 cm兩層取樣。將直徑50 mm、高100 mm的環(huán)刀垂直壓入土中，使環(huán)刀與土壤表層相平齊，用環(huán)刀托取出環(huán)刀后將土樣裝入自封袋，每個樣點土樣為3次環(huán)刀取樣混合而成。10—20 cm土樣采集時，用鐵鍬鏟去表層10 cm厚土壤，按相同方法采集10—20 cm的土樣。所有土樣帶回實驗室進行測試和分析。

2.3指標測定

將各土樣過0.5 mm的篩子，稱取枯落物的重量，根據(jù)環(huán)刀尺寸計算沙堆各部位枯落物百分含量。過篩后的土樣全部送實驗室進行有機質、全N、全P、全K、速效N、速效P和速效K共7個指標的測定。有機質采用重鉻酸鉀法測定，全氮采用凱氏法利用FOSS全自動定氮儀測定，全磷用DU800紫外分光光度計測定、全鉀用火焰光度計測定，速效氮用堿解擴散法測定，速效磷用鉬銻抗比色法測定，速效鉀用中性醋酸銨浸提法測定。

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

用SPSS 21.0軟件對各項指標進行描述統(tǒng)計分析和數(shù)據(jù)正態(tài)性檢驗(one sample Kolomogorov-Semimov，即K-S檢驗)。對同一灌叢沙堆不同部位之間以及不同沙堆相同部位之間的枯落物含量及土壤養(yǎng)分各指標數(shù)據(jù)進行單因素方差分析(One-way ANOVA)和多重比較(LSD)，并用Excel軟件進行分析和作圖。

3　研究結果

3.1枯落物含量變化規(guī)律

隨著植被蓋度的下降，0—10 cm土層枯落物含量總體呈降低的趨勢，而沙堆不同部位降低幅度又有不同(圖2)。對于沙堆灌叢下部位，沙堆A枯落物含量為2.21%，到沙堆D降低至0.94%，降幅達到57%;對于沙堆邊緣區(qū)，沙堆A枯落物含量為1.16%，至沙堆D則降低到0.38%，降幅達到67%;對于丘間地，枯落物含量在沙堆A為0.53%，到沙堆D則降至0.13%，最大降幅為75%;對于沙堆風影區(qū)，則從沙堆A的0.5%降至沙堆D 的0.11%，最大降幅為78%。就單個沙堆不同部位而言，沙堆枯落物含量以灌叢下最大，其次為沙堆邊緣區(qū)，丘間地和風影區(qū)較差，其中沙堆灌叢下顯著高于其它3個部位。將4個沙堆丘間地枯落物與灌叢下相比，沙堆A、B、C和D的降幅依次為為77%、81%、91%和81%。與0—10 cm土層相比，10—20 cm土層枯落物含量普遍降低，降幅為40%左右，其中灌叢下變化最為明顯。無論是不同沙堆間還是同一沙堆不同部位間，10—20 cm土層與0—10 cm土層相比，均具有相對一致的變化趨勢。

圖2　 土壤有機質含量變化規(guī)律Fig.2The changes of soil litter content(mean±S.E.) Different lowercase letters indicate significant statistical differences(P＜0.05);U:灌叢下Under shrub;M:沙堆邊緣Margin of nebkha;I:丘間地Interdune;L:風影區(qū)Leeward shadow;A、B、C和D依次代表植被蓋度為30%、15%—20%、10%和＜5%的典型樣地內所選的檉柳灌叢沙堆

3.2土壤有機質含量變化規(guī)律

隨著植被蓋度的下降，0—10 cm土層有機質含量總體呈逐漸降低的趨勢，而沙堆不同部位降低幅度又有不同(圖3)。對于沙堆灌叢下部位，沙堆A有機質含量為4.67 g/kg，到沙堆D降低至2.26 g/kg，降幅達到52%;對于沙堆邊緣區(qū)，沙堆A有機質含量為3.62 g/kg，至沙堆D則降低到2.31 g/kg，降幅達到36%;對于丘間地，有機質含量在沙堆A為2.50 g/kg，到沙堆D則降至1.93 g/kg，最大降幅為23%;對于沙堆風影區(qū)，有機質含量先降低后升高，最小值位于沙堆C，含量為1.93 g/kg。就單個沙堆不同部位而言，沙堆有機質含量以灌叢下最大，其次為沙堆邊緣區(qū)和風影區(qū)，丘間地最小，其中沙堆灌叢下顯著高于其它3個部位。與0—10 cm土層相比，沙堆A、B、C各部位10—20 cm土層有機質含量普遍降低，其中以灌叢下和沙堆邊緣兩部位比較明顯，平均降幅分別為23.3%和30.75%。沙堆D各部位10—20 cm土層較0—10 cm有所增加，平均增幅為9.66%。無論是不同沙堆間還是同一沙堆不同部位間，10—20 cm土層與0—10 cm土層相比，差異顯著性較枯落物含量均有減弱。

圖3　 土壤有機質含量變化規(guī)律Fig.3The changes of soil organic matter content(mean±S.E.)Different lowercase letters indicate significant statistical differences(P＜0.05);U:灌叢下Under shrub;M:沙堆邊緣Margin of nebkha;I:丘間地Interdune;L:風影區(qū)Leeward shadow;A、B、C和D依次代表植被蓋度為30%、15%—20%、10%和＜5%的典型樣地內所選的檉柳灌叢沙堆

3.3土壤全N、全P、全K含量變化規(guī)律

隨著植被蓋度的下降，0—10 cm土層土壤全N和全K含量在灌叢下、沙堆邊緣、丘間地3部位呈逐漸降低的趨勢，而不同部位降低幅度又有不同(圖4，圖5，圖6)。對于沙堆灌叢下部位，沙堆A全N和全K含量分別為0.26 g/kg和16.27 g/kg，到沙堆D分別降低至0.13 g/kg和14.73 g/kg，降幅分別為50%和9.4%;對于沙堆邊緣區(qū)，沙堆A全N和全K含量分別為0.20 g/kg和15.94 g/kg，至沙堆D則降低到0.14 g/kg和14.53 g/kg，降幅為30%和8.8%;對于丘間地，全N和全K含量在沙堆A分別為0.15 g/kg和15.39 g/kg，到沙堆D則降至0.11 g/kg和14.33 g/kg，降幅分別為26%和6.8%。對于沙堆風影區(qū)，全N含量先降低后升高，最小值位于沙堆B，含量為0.12 g/kg，沙堆D大于沙堆B和沙堆C。全K含量先升高后降低，最大值位于沙堆B，含量為16.22 g/kg。從沙堆A→沙堆B→沙堆C→沙堆D，0—10 cm土層土壤全P含量變化趨勢不明顯。就單個沙堆不同部位而言，沙堆A在灌叢下全N和全P含量與其它三部位之間具有顯著差異，沙堆B、C、D各部位之間差異不顯著，4個沙堆各部位全N和全K含量最小值均位于丘間地。與0—10 cm土層相比，10—20 cm土層全N和全P含量有所降低，平均降幅分別為25%和4.5%，土層全K含量略有增加，平均增幅為4.8%。無論是不同沙堆間還是同一沙堆不同部位間，10—20 cm土層與0—10 cm土層相比，土壤全N、全P、全K含量具有相對一致的變化趨勢。

圖4　 土壤全氮含量變化規(guī)律Fig.4The changes of soil total nitrogen content(mean±S.E.)Different lowercase letters indicate significant statistical differences(P＜0.05);U:灌叢下Under shrub;M:沙堆邊緣Margin of nebkha;I:丘間地Interdune;L:風影區(qū)Leeward shadow;A、B、C和D依次代表植被蓋度為30%、15%—20%、10%和＜5%的典型樣地內所選的檉柳灌叢沙堆

3.4土壤速效N、速效P、速效K含量變化規(guī)律

圖5　土壤全磷含量變化規(guī)律Fig.5The changes of soil total phosphorus content(mean±S.E.)Different lowercase letters indicate significant statistical differences(P＜0.05);U:灌叢下Under shrub;M:沙堆邊緣Margin of nebkha;I:丘間地Interdune;L:風影區(qū)Leeward shadow;A、B、C和D依次代表植被蓋度為30%、15%—20%、10%和＜5%的典型樣地內所選的檉柳灌叢沙堆

圖6　土壤全鉀含量變化規(guī)律Fig.6The changes of soil total potassium content(mean±S.E.)Different lowercase letters indicate significant statistical differences(P＜0.05);U:灌叢下Under shrub;M:沙堆邊緣Margin of nebkha;I:丘間地Interdune;L:風影區(qū)Leeward shadow;A、B、C和D依次代表植被蓋度為30%、15%—20%、10%和＜5%的典型樣地內所選的檉柳灌叢沙堆

圖7　土壤速效氮含量變化規(guī)律Fig.7The changes of soil available nitrogen content(mean±S.E.)Different lowercase letters indicate significant statistical differences(P＜0.05);U:灌叢下Under shrub;M:沙堆邊緣Margin of nebkha;I:丘間地Interdune;L:風影區(qū)Leeward shadow;A、B、C和D依次代表植被蓋度為30%、15%—20%、10%和＜5%的典型樣地內所選的檉柳灌叢沙堆

圖8　土壤速效磷含量變化規(guī)律Fig.8The changes of soil available phosphorus content(mean±S.E.)Different lowercase letters indicate significant statistical differences(P＜0.05);U:灌叢下Under shrub;M:沙堆邊緣Margin of nebkha;I:丘間地Interdune;L:風影區(qū)Leeward shadow;A、B、C和D依次代表植被蓋度為30%、15%—20%、10%和＜5%的典型樣地內所選的檉柳灌叢沙堆

隨著植被蓋度的減小，0—10 cm土層土壤速效N含量在灌叢下和沙堆邊緣兩部位呈先減小后增大的趨勢，丘間地和風影區(qū)兩部位呈先減小后增大再減小的趨勢(圖7，圖8，圖9)。土壤速效P含量總體呈先減小后增大的規(guī)律，沙堆A 4個部位速效P含量均大于其它3個沙堆。土壤速效K含量總體呈逐漸降低的趨勢，而沙堆不同部位降低幅度又有不同。對于沙堆灌叢下部位，沙堆A土壤速效K含量為473.9 mg/kg，到沙堆D降低至311.31 mg/kg，降幅達到34.3%;對于沙堆邊緣區(qū)，沙堆含量為402.5 mg/kg，至沙堆D則降低到272.84 mg/kg，降幅達到32.21%;對于丘間地，速效K含量在沙堆A為287.49 mg/kg，到沙堆D則降至235.72 mg/ kg，降幅為18%;對于沙堆風影區(qū)，則從沙堆A的316 mg/kg降至沙堆D的273.9 mg/kg，降幅為13.32%。就單個沙堆而言，沙堆C灌叢下土壤速效N含量與丘間地、風影區(qū)兩部位之間具有顯著差異，其余沙堆各部位之間沒有顯著差異。土壤速效P含量4個沙堆4個部位之間均沒有顯著性差異。4個沙堆從灌叢下→沙堆邊緣→丘間地→風影區(qū)速效K含量大體呈先減小后增大的趨勢，最小值位于丘間地。

圖9　 土壤速效鉀含量變化規(guī)律Fig.9The changes of soil available potassium content(mean±S.E.) Different lowercase letters indicate significant statistical differences(P＜0.05);U:灌叢下Under shrub;M:沙堆邊緣Margin of nebkha;I:丘間地Interdune;L:風影區(qū)Leeward shadow;A、B、C和D依次代表植被蓋度為30%、15%—20%、10%和＜5%的典型樣地內所選的檉柳灌叢沙堆

與0—10 cm土層相比，10—20 cm土層土壤速效N含量明顯增加，而沙堆不同部位增加幅度又有不同。對于沙堆灌叢下部位，0—10 cm土層速效N平均含量為36.75 mg/kg，到10—20 cm土層增加至47.5 mg/kg，增幅達到29.25%;對于沙堆邊緣區(qū)，0—10 cm土層含量為35 mg/kg，至10—20 cm土層增加到75 mg/kg，增幅達到114%;對于丘間地，0—10 cm土層含量為45 mg/kg，到10—20 cm土層增加至97.25 mg/kg，最大增幅為116%;對于沙堆風影區(qū)，則從0—10 cm土層的46.5 mg/kg增至10—20 cm土層的118.25 mg/kg，最大增幅為154%。10—20 cm土層速效N含量，在4個沙堆均表現(xiàn)出從灌叢下→沙堆邊緣→丘間地→風影區(qū)土壤逐漸增大的趨勢。土壤速效K含量也有所增加，平均增幅為12.06%，其中以丘間地增加比較明顯，增幅為22%。土壤速效P含量略有降低，平均降幅為2%，其中沙堆A降低比較明顯，降幅為11.4%。不同沙堆間土壤速效P含量與0—10 cm土層相比，具有相對一致的變化趨勢。

4　討論

本項研究發(fā)現(xiàn)，隨著植被蓋度的降低，沙堆各部位枯落物及有機質含量呈降低趨勢，這主要是由于植被覆蓋條件和風沙活動相互作用的結果。楊東亮等人［30］對研究區(qū)典型下墊面風沙活動的詳細研究發(fā)現(xiàn)，駱駝刺群落分布區(qū)風沙活動強度明顯比其它區(qū)域弱，一方面植被條件較好會產(chǎn)生較多的枯落物覆蓋地表，另一方面風沙活動弱，致使枯落物和有機質能夠在灌叢下大量的累積。隨著植被條件變差，對風力的阻滯作用減弱，風沙活動強度增加，地表由積沙狀態(tài)變?yōu)轱L蝕狀態(tài)，勢必造成枯落物和有機質累積量的減少，肥島效應趨于消失。

有學者［31-36］對檉柳灌叢沙堆表面壓力及流場分布特征研究發(fā)現(xiàn)，沙堆表面風蝕主要發(fā)生在沙堆迎風坡及兩側區(qū)域，而沙堆背風區(qū)域和灌叢保護的沙堆頂部是積沙發(fā)生的主要區(qū)域。本項研究發(fā)現(xiàn)就單個灌叢沙堆而言，枯落物和有機質含量以灌叢下最大，其次為沙堆邊緣，再次為風影區(qū)和丘間地，這與沙堆各部位的蝕積狀況相對應。土壤有機質、全N、全K、速效K含量與枯落物含量表現(xiàn)出相同的變化趨勢，說明枯落物是該研究區(qū)表層土壤養(yǎng)分的主要來源，枯落物的累積對養(yǎng)分的富集有很大的影響。土壤有機質、全N、全K、速效P、速效K含量最大值都位于灌叢下部位，這除了枯落物分解的貢獻外，還與灌叢對風速的降低和降塵的捕獲有關。灌叢沙堆4個部位土壤有機質、全N、全P、全K、速效K含量的最小值都位于丘間地，主要是由于丘間地風蝕作用較強烈，致使土壤細粒物質和養(yǎng)分流失。

隨著植被蓋度的降低風蝕作用普遍增強，在重度沙漠化區(qū)域表現(xiàn)最為強烈，沙物質處于頻繁的交換狀態(tài)［15］，所以本研究發(fā)現(xiàn)沙堆D土壤各養(yǎng)分指標4個部位之間沒有顯著性差異。Li等［37］研究發(fā)現(xiàn)直徑小于50 um的土壤微粒含有較高的有機碳和全氮，本項研究中灌叢沙堆風影區(qū)部位沙物質主要以沉積為主，細粒物質最容易在此部位沉積，灌叢沙堆土壤有機質、全N、全P、速效N、速效P、速效K含量從灌叢下→沙堆邊緣→丘間地→風影區(qū)表現(xiàn)出先減小后增大的趨勢。由于沙堆D所在區(qū)域風沙活動頻繁，導致灌叢沙堆沙物質也處于頻繁的交換運移狀態(tài)，灌叢沙堆10—20 cm土層有機質、全N、速效N、速效P、速效K含量相比0—10 cm土層各部位之間差異性明顯減弱。

在生物過程與非生物過程的交互作用下，干旱區(qū)土壤資源具有較強的空間異質性，肥島現(xiàn)象正是這種異質性的充分體現(xiàn)［38］。生物作用大于非生物作用是荒漠生態(tài)系統(tǒng)中土壤資源在灌叢下富集的重要原因［39］。本研究表明，在植被覆蓋較好的條件下，灌叢植物對降塵的捕獲、對風速的明顯降低，致使枯落物和有機質累積。隨著植被的大量破壞，沙漠化程度加劇，非生物作用大于生物作用時肥島現(xiàn)象則趨于消失。本項研究沙堆D養(yǎng)分含量的變化規(guī)律都證實了這種說法。

荒漠化土壤養(yǎng)分的變化是荒漠化過程中能流、物流中的重要一環(huán)［11］。土壤養(yǎng)分的供給與植物生長之間的平衡關系決定了生態(tài)環(huán)境演化的過程，荒漠化過程是這種平衡關系被破壞的表現(xiàn)。在脆弱的生態(tài)環(huán)境中，土壤養(yǎng)分與植被之間存在著密切的聯(lián)系，土壤養(yǎng)分的變化直接影響荒漠化的程度與發(fā)展［40］。有學者研究指出長期的放牧會引起水分、氮及其他土壤資源的異質性，這種異質性會促進荒漠灌木的入侵，導致土壤資源在灌木下的長期累積，灌叢間裸地的土壤則因侵蝕而流失，這種物理與生物的雙重作用導致灌叢島的形成，同時會導致土地退化［16，41］。李新榮［42］研究發(fā)現(xiàn)灌木的定居會引起土壤資源的空間異質性變化，在流動沙地的恢復治理過程中可以利用這一原理打破流動沙地質地均一、養(yǎng)分貧瘠的特點，將有限的資源集中在灌木叢下，促進灌木的生長從而實現(xiàn)流沙地的固定。這也是在荒漠化治理過程中選擇旱生灌木種作為先鋒植物的原因。本項研究發(fā)現(xiàn)在策勒綠州沙漠過渡帶荒漠化過程中，隨著植被蓋度的降低，風沙活動不斷加劇，致使養(yǎng)分呈現(xiàn)有的這種分布規(guī)律，但是這種養(yǎng)分的分布對該研究區(qū)的植被將產(chǎn)生怎樣的影響有待進一步的研究。

5　結論

(1)隨著植被蓋度的降低，灌叢沙堆表層0—10 cm各部位土壤枯落物和有機質呈降低趨勢，平均降幅為69.3%和37.0%。就單個灌叢沙堆而言，枯落物和有機質含量以灌叢下最大，其次為沙堆邊緣和風影區(qū)，丘間地最小。10—20 cm土層枯落物和有機質含量相對于0—10 cm土層明顯降低，平均降幅為40.0%和27.0%，沙堆各部位之間的差異顯著性均有減弱，但兩個土層的變化趨勢保持一致。

(2)隨著植被蓋度的降低，灌叢沙堆表層0—10 cm土壤全N含量在灌叢下、沙堆邊緣、丘間地三部位呈逐漸降低的趨勢，平均降幅為50.0%、30.0%和26.0%。在沙堆風影區(qū)全N含量呈先降低后升高的趨勢。灌叢沙堆土壤全K含量呈降低趨勢，平均降幅為8.3%。土壤全P含量沒有明顯的變化趨勢。就單個灌叢沙堆而言，全N、全K含量以灌叢下最大，其次為沙堆邊緣和風影區(qū)，丘間地最小。與0—10 cm土層相比，10—20 cm土層全N和全P含量有所降低，平均降幅為25.0%和4.5%，土層全K含量略有增加，平均增幅為4.8%。但兩個土層的變化趨勢保持一致。

(3)隨著植被蓋度的降低，灌叢沙堆表層0—10 cm土壤速效N和速效P含量各部位大體呈先減小后增大的趨勢，各部位之間沒有顯著的差異性，土壤速效K呈降低趨勢，灌叢下、沙堆邊緣、丘間地和風影區(qū)各部位的降幅為34.3%、32.2%、18.0%、13.3%。與0—10 cm土層相比，10—20 cm土層土壤速效N和速效K含量明顯增加，平均增幅為103.3%和12.1%，風影區(qū)速效N含量明顯高于其他3個部位。土壤速效P含量略有降低，平均降幅為2%。

(4)在荒漠綠洲過渡帶，輕度沙漠化和中度沙漠化區(qū)域，灌叢的生物反饋作用使得灌叢下部位具有明顯的養(yǎng)分富集效應，但隨著植被總蓋度的降低，風沙活動不斷加劇，導致養(yǎng)分富集效應逐漸減弱，肥島效應消失。

參考文獻(References):

［1］Tengberg A.Nebkha dunes as indicators of wind erosion and land degradation in the Sahel zone of Burkina Faso.Journal of Arid Environments，1995，30(3):265-282.

［2］穆桂金.塔克拉瑪干沙漠灌草丘的發(fā)育特征及環(huán)境意義.干旱區(qū)研究，1994，11(1):34-40.

［3］Nickling W G，Wolfe S A.The morphology and origin of Nabkhas，region of Mopti，Mali，West Africa.Journal of Arid Environments，1994，28 (1):13-30.

［4］Mares M A.Encyclopedia of Deserts.Norman:University of Oklahoma Press，1999:189-189.

［5］Tengberg A，Chen D L.A comparative analysis of nebkhas in central Tunisia and northern Burkina Faso.Geomorphology，1998，22(2):181-192.

［6］Ardon K，Tsoar H，Blumberg D G.Dynamics of nebkhas superimposed on a parabolic dune and their effect on the dune dynamics.Journal of Arid Environments，2009，73(11):1014-1022.

［7］Li Z Z，Wu S L，Chen S J，Chen X L，Jin J H，Liu Q.Bio-geomorphologic features and growth process of Tamarix nabkhas in Hotan River Basin，Xinjiang.Journal of Geographical Sciences，2010，20(2):205-218.

［8］Qong M，Takamura H，Hudaberdi M.Formation and internal structure of Tamarix cones in the Taklimakan Desert.Journal of Arid Environments，2002，50(1):81-97.

［9］武勝利，李志忠，肖晨曦，孫秋梅，劉麗梅.灌叢沙堆的研究進展與意義.中國沙漠，2006，26(5):734-738.

［10］朱震達，陳廣庭.中國土地沙質荒漠化.北京:科學出版社，1994，24-24.

［11］Klemmedson J O，Barth R C.Distribution and balance of biomass and nutrients in desert shrub ecosystems.US/IBP Desert Biome Research Memo.Logan:Utah State University Press，1975:18-18.

［12］陳廣生，曾德慧，陳伏生，范志平，耿海東.干旱和半干旱地區(qū)灌木下土壤“肥島”研究進展.應用生態(tài)學報，2003，14(12):2295-2300.

［13］Charley J L，West N E.Plant-induced soil chemical patterns in some shrub-dominated semi-desert ecosystems of Utah.Journal of Ecology，1975，63(3):945-963.

［14］Schlesinger W H，Raikes J A，Hartley A E，Cross A F.On the spatial pattern of soil nutrients in desert ecosystems.Ecology，1996，77(2): 364-374.

［15］Li P X，Wang N，He W M，Krüsi B O，Gao S Q，Zhang S M，Yu F H，Dong M.Fertile islands under Artemisia ordosica in inland dunes of northern China:Effects of habitats and plant developmental stages.Journal of Arid Environments，2008，72(6):953-963.

［16］Schlesinger W H，Pilmanis A M.Plant-soil interactions in deserts.Biogeochemistry，1998，42(1/2):169-187.

［17］Zaady E.Seasonal change and nitrogen cycling in a patchy Negev Desert:a review.Arid Land Research and Management，2005，19(2):111-124.

［18］Smith W K，Monson R K，Anderson J E.Physiological Ecology of North American Desert Plants:with 86 figures and 21 tables.New York，NY: Springer，1997.

［19］Aguiar M R，Sala O E.Patch structure，dynamics and implications for the functioning of arid ecosystems.Trends in Ecology and Evolution，1999，14(7):273-277.

［20］Dhillion S S.Environmental heterogeneity，animal disturbances，microsite characteristics，and seedling establishment in a Quercus havardiicommunity.Restoration Ecology，1999，7(4):399-406.

［21］Gibbens R P，Lenz J M.Root systems of some Chihuahuan Desert plants.Journal of Arid Environments，2001，49(2):221-263.

［22］West N E，Skujins J.Nitrogen in Desert Ecosystems.Stroudsburg，PA:Dowden，Hutchinson and Ross Inc.，1978.

［23］Allington G R H，Valone T J.Reversal of desertification:the role of physical and chemical soil properties.Journal of Arid Environments，2010，74 (8):973-977.

［24］Titus J H，Nowak R S，Smith S D.Soil resource heterogeneity in the Mojave Desert.Journal of Arid Environments，2002，52(3):269-292.

［25］尹傳華.塔克拉瑪干沙漠北部檉柳灌叢肥島、鹽島的特征及其與荒漠化過程的關系［D］.北京:中國農業(yè)大學，2007.

［26］Okin G S，Murray B，Schlesinger W H.Degradation of sandy arid shrubland environments:observations，process modelling，and management implications.Journal of Arid Environments，2001，47(2):123-144.

［27］Li J，Zhao C，Zhu H，Li Y，Wang F.Effect of plant species on shrub fertile island at an oasis-desert ecotone in the South Junggar Basin，China.Journal of Arid Environments，2007，71(4):350-361.

［28］楊帆，王雪芹，楊東亮，韓章勇，鄭新倩.風沙侵襲影響下綠洲邊緣植物群落點格局分析.干旱區(qū)研究，2012，29(5):763-769.

［29］楊帆，王雪芹，楊東亮，韓章勇.不同沙源供給條件下檉柳灌叢與沙堆形態(tài)的互饋關系—以策勒綠洲沙漠過渡帶為例.生態(tài)學報，2012，32(9):2707-2719.

［30］楊東亮.策勒綠洲-沙漠過渡帶典型下墊面風沙活動研究［D］.北京:中國科學院研究生院，2011.

［31］武勝利，李志忠，惠軍，肖晨曦，孫秋梅.灌叢沙堆表面壓力分布特征的實驗研究.干旱區(qū)地理，2006，29(6):790-796.

［32］李志忠，武勝利，肖晨曦，劉麗梅，孫秋梅，張中華.新疆和田河流域灌叢沙堆風洞流場的實驗研究(Ⅱ).中國沙漠，2007，27(1): 15-19.

［33］李志忠，武勝利，肖晨曦，劉麗梅，孫秋梅，張中華.新疆和田河流域灌叢沙堆風洞流場的實驗研究(Ⅰ).中國沙漠，2007，27(1): 9-14.

［34］Judd M J，Raupach M R，F(xiàn)innigan J J.A wind tunnel study of turbulent flow around single and multiple windbreaks，part I:velocity fields.Boundary-Layer Meteorology，1996，80(1/2):127-165.

［35］凌裕泉，屈建軍，金炯.稀疏天然植被對輸沙量的影響.中國沙漠，2003，23(1):12-17.

［36］李志忠，武勝利，王曉峰，何牡丹，葛琳，穆赫塔爾，徐國強.新疆和田河流域檉柳沙堆的生物地貌發(fā)育過程.地理學報，2007，62(5): 462-470.

［37］Li J R，Okin G S，Epstein H E.Effects of enhanced wind erosion on surface soil texture and characteristics of windblown sediments.Journal of Geophysical Research:Biogeosciences，2009，114(G2):1-9.

［38］Wezel A，Rajot J L，Herbrig C.Influence of shrubs on soil characteristics and their function in Sahelian agro-ecosystems in semi-arid Niger.Journal of Arid Environments，2000，44(4):383-398.

［39］Garner W，Steinberger Y.A proposed mechanism for the formation of fertile islands in the desert ecosystem.Journal of Arid Environments，1989，16(3):257-262.

［40］呂貽忠，馬興旺.荒漠化土壤養(yǎng)分變化的影響因素研究進展.生態(tài)環(huán)境，2003，12(4):473-477.

［41］Schlesinger W H，Reynolds J F，Cunningham G L，Huenneke L F，Jarrell W M，Virginia R A，Whitford W G.Biological feedbacks in global desertification.Science，1990，247(4946):1043-1048.

［42］李新榮.干旱沙區(qū)土壤空間異質性變化對植被恢復的影響.中國科學(D輯:地球科學)，2005，35(4):361-370.

Spatial variation of soil nutrients of Tamarix ramosissima nebkhas and interdune areas in a desert-oasis ecotone

LIU Jinhui1，2，3，WANG Xueqin1，*，MA Yang1，2，3

1 Xinjiang Institute of Ecology and Geography，Chinese Academy of Sciences，ürümqi 830011，China

2 Cele National Station of Observation＆Research for Desert Grassland Ecosystem，Cele 848300，China

3 University of Chinese Academy of Sciences，Beijing 100049，China

Abstract:In this study four kinds of typical plots，with vegetation cover of 30%，15%—20%，10%，and＜5%，respectively，were selected to investigate the spatial variation of soil nutrients in Tamarix ramosissima nebkhas and interdune areas in a Qira oasis-desert ecotone at the southern rim of the Taklimakan desert.Soil samples from two layers (0—10 cm and 10—20 cm)of the nebkha and interdune in the four plots were collected for analysis of the litter，soil organic matter，total N，total P，total K，available N，available P，and available K content.The following results were obtained:First，with decreasing vegetation cover from Plot 1 to Plot 4，litter content，soil organic matter content，total Nbook=980,ebook=92content，total K content，and available K content in the top(0—10 cm)layer of the nebkhas decreased by 69.3%，37.0%，35.3%，8.3%，and 24.5%，respectively.The available N and available P content decreased at first，and then increased with decrease in vegetation cover，while total P content showed no obvious change.Second，the litter content，soil organic matter content，total N content，total P content，available N content，available P content，and available K content in the top layer first decreased，and then increased from under shrub to the margin of the nebkha，to the interdune，to the leeward shadow.The maximum values were mainly found under shrub，while the minimum values were mainly found in the interdune regions.For each soil nutrient index，there was no significant difference between the four positions in Plot 4 (vegetation cover is＜5%).Third，compared with the top soil layer，there was a significant reduction in litter content，organic matter content，total N content，total P content，total K content，available P content，and available K content in the bottom(10—20 cm)soil layer.The above corresponding mean values were lower by 40.0%，27.0%，25.0%，4.5%，and 2.0%，respectively，while total K content，available N content，and available K content were higher by 4.8%，103.3%，and 12.1%，respectively.With decreasing vegetation coverage，the changes in soil nutrient content in the two soil layers showed a consistent trend.Fourth，in non-desertified and lightly desertified areas of the oasis-desert ecotone，the biological feedback of shrubs led to obvious nutrient enrichment effects around them.However，as the total vegetation cover decreased and the intensity of desertification increased，non-biological factors began to play leading roles，and the nutrient enrichment effects of shrubs tended to disappear gradually.

Key Words:desert-oasis ecotone;Tamarix ramosissima nebkhas;litter content;soil nutrient

*通訊作者

Corresponding author.E-mail:xqwang@ms.xjb.ac.cn

收稿日期:2014-06-03;網(wǎng)絡出版日期:2015-07-09

基金項目:國家自然科學基金項目(41371042);國家科技支撐計劃(2014BAC14B02)

DOI:10.5846/stxb201406031141