郭艷菊,馬曉靜,許愛(ài)云,史佳梅,劉萬(wàn)龍,許冬梅

寧夏大學(xué)農(nóng)學(xué)院, 銀川 750021

草地生態(tài)系統(tǒng)是陸地生態(tài)系統(tǒng)的重要組分,其能量流動(dòng)、物質(zhì)循環(huán)和自我修復(fù)對(duì)維系區(qū)域生態(tài)平衡、發(fā)展地區(qū)經(jīng)濟(jì)具有重要價(jià)值。然而,在氣候變化和人類活動(dòng)的雙重作用下,干旱、半干旱地區(qū)形成了以風(fēng)蝕為主要特征的土地沙化,其發(fā)生面積、危害程度已遠(yuǎn)遠(yuǎn)超出其他類型的土地退化[1—2]。土地沙化影響土壤養(yǎng)分的累積與分解、團(tuán)聚體的形成及植被生產(chǎn)力、植物群落結(jié)構(gòu),生態(tài)系統(tǒng)自我修復(fù)能力減弱[3—5]。干旱、半干旱地區(qū)降雨量少、蒸散量高、土壤保水能力差,加之難以得到地下水的補(bǔ)充,導(dǎo)致土壤水資源缺乏[6]。而土壤水分影響土壤侵蝕、碳氮轉(zhuǎn)化及植物的生長(zhǎng)、生理生態(tài)特性和群落結(jié)構(gòu),對(duì)沙地土壤的演化具有重要作用[7]。因此,探究沙化草地土壤水分的分異規(guī)律,對(duì)沙化草地生態(tài)系統(tǒng)的保護(hù)和修復(fù)具有重要的理論和現(xiàn)實(shí)意義。

土壤水分受土壤質(zhì)地、成土母質(zhì)、氣候條件、地形地貌、土地利用方式和植被覆蓋類型等因素的綜合作用,在垂直方向和水平方向表現(xiàn)出明顯的空間變異[8—11]。降水作為干旱、半干旱區(qū)草地生態(tài)系統(tǒng)唯一的水分來(lái)源,直接影響土壤水分含量及草地生產(chǎn)力[12]。在沙化草地生態(tài)系統(tǒng)中,地形通過(guò)起伏變化改變由降水格局產(chǎn)生的土壤水分空間分異規(guī)律,往往形成沙丘底部土壤水分含量最高、坡中次之、坡頂最低的空間分布格局[13]。土壤水分的空間變異影響植物的生長(zhǎng)發(fā)育及分布規(guī)律,而植被對(duì)土壤性狀也產(chǎn)生反饋?zhàn)饔?由植被控制的土壤水分分異是許多干旱、半干旱生態(tài)系統(tǒng)土壤水分變化的重要影響因素[6—7]。在沙地生態(tài)系統(tǒng)灌叢群落演替初期,土壤水分的空間異質(zhì)性決定灌木物種的分布格局,至群落相對(duì)穩(wěn)定的固定沙地階段,灌叢的分布格局一定程度上改變了土壤水分的空間分布特征,并對(duì)土壤水分的變化做出響應(yīng)[14]。此外,在干旱沙區(qū),灌木通過(guò)影響灌叢下土壤水分的空間分異,導(dǎo)致灌叢下草本植物的生物量呈環(huán)狀分布[15]。因此,沙化草地生態(tài)系統(tǒng)降水、徑流等水文過(guò)程決定土壤水分的分異,而由土壤水分變化驅(qū)動(dòng)的植被格局又反饋于土壤,引起土壤水分的再分配。

寧夏東部風(fēng)沙區(qū)地處我國(guó)北方農(nóng)牧交錯(cuò)的生態(tài)脆弱帶,為黃土高原向鄂爾多斯臺(tái)地的過(guò)渡地帶,也是半干旱向干旱氣候的過(guò)渡帶,在全球氣候變化和人為干擾因素共同作用下,植被生態(tài)系統(tǒng)穩(wěn)定性較低和土壤不穩(wěn)定,生態(tài)演替往往呈現(xiàn)出雙向性,以往對(duì)寧夏東部風(fēng)沙區(qū)沙化草地的研究主要集中在沙化草地類型的劃分[16]、植物光合生理生態(tài)[17]、優(yōu)勢(shì)物種替代[18]及土壤養(yǎng)分[19—21]、微生物數(shù)量[22]和酶活性[23]等對(duì)草地沙化的響應(yīng),針對(duì)沙化草地小尺度范圍內(nèi)植被和土壤的空間異質(zhì)性也有一些研究[14,24]。然而,針對(duì)動(dòng)態(tài)變化的沙化草地生態(tài)系統(tǒng)其植被和土壤水分在沙化演替過(guò)程中的空間變異規(guī)律及其互饋關(guān)系如何,目前還鮮有報(bào)道。本文基于草地生態(tài)系統(tǒng)在寧夏東部風(fēng)沙區(qū)的生態(tài)地位,采用經(jīng)典統(tǒng)計(jì)學(xué)和地統(tǒng)計(jì)學(xué)方法,研究自潛在沙化草地至重度沙化草地土壤水分、植被蓋度和物種豐富度指數(shù)的變化,探討寧夏東部風(fēng)沙區(qū)草地沙化演替過(guò)程中土壤水分和植被的空間分異規(guī)律及其相互關(guān)系,以期為草地沙化演變過(guò)程中生物、非生物因子的定量化測(cè)度提供依據(jù),對(duì)明晰草地沙化敏感區(qū)域及干旱風(fēng)沙區(qū)生態(tài)環(huán)境保護(hù)和建設(shè)具有指導(dǎo)意義。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于寧夏東部鹽池縣哈巴湖自然保護(hù)區(qū)(E107°5′45″—107°6′7″,N37°43′34″—37°44′30″),毛烏素沙地西南緣(圖1),地貌為緩坡丘陵,平均海拔1501 m,年均降水量302.38 mm,70%以上集中在6—9月份,年均潛在蒸發(fā)量2131.7 mm。年平均日照時(shí)數(shù)2876 h,1月平均氣溫-6.76 ℃,7月平均氣溫23.91 ℃,年平均氣溫9.93 ℃,年平均無(wú)霜期162 d。土壤類型主要為灰鈣土、風(fēng)沙土,土壤質(zhì)地為沙壤和粉沙壤。研究區(qū)地處中我國(guó)溫帶草原的過(guò)渡地帶,屬歐亞草原區(qū),植被組成以黑沙蒿(Artemisiaordosica)、短花針茅(Stipabreviflora)、豬毛蒿(Artemisiascoparia)、阿爾泰狗娃花(Heteropappusaltaicus)、賴草(Leymussecalinus)、糙隱子草(Cleistogenessquarrosa)等為主。

圖1 研究區(qū)地理位置

1.2 野外調(diào)查與樣品采集方法

于2018年7月選取哈巴湖自然保護(hù)區(qū)典型區(qū)域?yàn)檠芯繉?duì)象,依據(jù)《天然草地退化、沙化、鹽漬化的分級(jí)指標(biāo)》(GB19733—2003),將研究區(qū)劃分為潛在沙化草地、輕度沙化草地、中度沙化草地和重度沙化草地4種沙化草地類型[25]。采用樣線法,自潛在沙化草地(北)至重度沙化草地(南)設(shè)置3條平行樣線,樣線長(zhǎng)為3.6 km,相鄰兩條樣線之間的間隔距離為500 m,在每條樣線上每隔200 m設(shè)置樣方(草本植物1 m×1 m,灌木10 m×10 m),共計(jì)57個(gè),調(diào)查植物物種組成、密度和蓋度(圖2)。同時(shí),在每個(gè)樣方,去除地表覆蓋物后,挖掘土壤剖面,采用多點(diǎn)混合法按0—10 cm、10—20 cm、20—30 cm、30—40 cm、40—50 cm、50—60 cm、60—70 cm、70—80 cm、80—90 cm和90—100 cm分層采集土壤樣品,將同一土層樣品混合后帶回實(shí)驗(yàn)室,采用烘干法測(cè)定土壤水分(SM)。土壤樣品采集保證在當(dāng)?shù)赜行Ы涤?>5 mm降水)后的 20 d左右集中進(jìn)行,以消除降水等因素對(duì)土壤水分的影響。

圖2 樣方分布示意圖

1.3 數(shù)據(jù)計(jì)算與分析

物種豐富度指數(shù)(R)采用Margalef指數(shù),計(jì)算公式如下:

R=(S-1)/lnN

式中,N為樣方中全部物種的個(gè)體總數(shù)；S為樣方中物種數(shù)。

采用Excel 2010進(jìn)行數(shù)據(jù)的基礎(chǔ)處理；采用SPSS 22.0進(jìn)行經(jīng)典統(tǒng)計(jì)；采用R 4.0.2進(jìn)行Pearson相關(guān)分析；利用GS+ 9.0進(jìn)行土壤水分、植被蓋度及物種豐富度指數(shù)的空間異質(zhì)性分析。半變異函數(shù)計(jì)算公式如下:

式中,r(h)為半變異函數(shù),h為兩樣點(diǎn)間空間間隔距離,N(h)為間隔距離為h時(shí)的樣點(diǎn)對(duì)的總數(shù),Z(xi)代表樣點(diǎn)xi處的實(shí)測(cè)值,Z(xi+h)代表與xi間距為h時(shí)樣點(diǎn)的實(shí)測(cè)值,半變異函數(shù)有基臺(tái)值(C0+C)、塊金值(C0)、塊金系數(shù)(C0/(C0+C))和變程(A0)4個(gè)重要參數(shù)。塊金值(C0)表示由于測(cè)量誤差及小于最小取樣尺度引起的隨機(jī)變異,其值越大表明較小尺度上的某些過(guò)程越不能忽視；基臺(tái)值(C0+C)反映系統(tǒng)內(nèi)總的變異,包括結(jié)構(gòu)性變異和隨機(jī)性變異兩部分。塊金系數(shù)(C0/(C0+C))表示隨機(jī)部分引起的變異占總變異的比例,C0/(C0+C) 小于25%表明在空間上呈強(qiáng)烈的相關(guān)性,主要受結(jié)構(gòu)性因素的影響；介于25%和75%之間,表明在空間上呈中等相關(guān)性；大于75%時(shí),在空間上呈弱相關(guān)性,主要受隨機(jī)性因素的影響。變程A0也稱空間相關(guān)范圍,能夠準(zhǔn)確反映區(qū)域化變量自相關(guān)范圍的大小[26]。

利用普通克里格(Ordinary Kriging)法進(jìn)行空間插值估計(jì),并進(jìn)行交叉實(shí)證檢驗(yàn)。依據(jù)克里格插值結(jié)果,繪制土壤水分及植被蓋度、物種豐富度的空間分布圖。

2 結(jié)果

2.1 沙化草地土壤水分、植被蓋度及物種豐富度指數(shù)的描述性統(tǒng)計(jì)分析

自潛在沙化草地至重度沙化草地,0—100 cm土層土壤水分含量分布范圍為0.82%—28.22%,其中以70—100 cm深層土壤變異較大,極大值和極小值的差值范圍為22.82%—25.58%；隨剖面深度的增加,土壤水分含量總體呈先升高后降低的趨勢(shì),以70—80 cm和80—90 cm土層最高。植被蓋度和物種豐富度指數(shù)的變化范圍分別為41.00%—93.00%和0.82—2.80(圖3,表1)。描述性統(tǒng)計(jì)分析結(jié)果表明,0—100 cm各土層土壤水分的變異系數(shù)介于0.36—0.48,植被蓋度和物種豐富度指數(shù)的變異系數(shù)分別為0.20和0.26,均屬于中等變異。各土層土壤水分、植被蓋度及物種豐富度指數(shù)的偏度均為正值,表明數(shù)據(jù)在右側(cè)分散。土壤水分的峰度介于-0.67—3.30,表明數(shù)據(jù)分布較為陡峭,存在極值；植被蓋度及物種豐富度指數(shù)的峰度均為負(fù)值,數(shù)據(jù)分布的峰都低于標(biāo)準(zhǔn)正態(tài)。采用Kolmogorov-Smirnov(K-S) 法檢驗(yàn),將K-S值小于0.05的數(shù)據(jù)進(jìn)行了自然對(duì)數(shù)轉(zhuǎn)換,轉(zhuǎn)換后的數(shù)據(jù)均符合正態(tài)分布 (P>0.05)。

表1 寧夏東部風(fēng)沙區(qū)沙化草地不同土層土壤水分、植被蓋度及物種豐富度指數(shù)的描述性統(tǒng)計(jì)

圖3 寧夏東部風(fēng)沙區(qū)沙化草地0—100 cm 土層土壤水分、植被蓋度及物種豐富度指數(shù)的分布

2.2 沙化草地土壤水分、植被蓋度及物種豐富度指數(shù)的空間異質(zhì)性

半變異函數(shù)模型最優(yōu)擬合表明,0—10 cm土層土壤水分和物種豐富度指數(shù)的最優(yōu)模型分別為球狀模型和指數(shù)模型,10—100 cm各土層土壤水分和植被蓋度的最優(yōu)模型均為高斯模型(圖4),決定系數(shù)為0.593—0.941,且殘差較小,曲線擬合效果較好。各土層土壤水分和物種豐富度指數(shù)的塊金系數(shù)為26.31%—41.86%,介于25%—75%之間,表現(xiàn)為中等的空間自相關(guān)性。植被蓋度的塊金系數(shù)為2.50%,小于25%,表現(xiàn)出強(qiáng)烈的空間自相關(guān)性。各土層土壤水分、植被蓋度及物種豐富度指數(shù)的自相關(guān)范圍較大,變程為1453—8110 m (表2) 。0—10 cm表層土壤水分的塊金值、變程及塊金系數(shù)均較高,隨土層的加深,土壤水分的塊金值呈波動(dòng)性下降趨勢(shì),基臺(tái)值則無(wú)明顯變化規(guī)律。

表2 寧夏東部風(fēng)沙區(qū)沙化草地土壤水分、植被蓋度及物種豐富度指數(shù)的半變異函數(shù)理論模型及相關(guān)參數(shù)

圖4 寧夏東部風(fēng)沙區(qū)沙化草地0—100 cm土層土壤水分、植被蓋度及物種豐富度指數(shù)的半變異函數(shù)r(h)-h散點(diǎn)圖

2.3 沙化草地土壤水分、植被蓋度及物種豐富度指數(shù)的空間分布特征

采用普通 Kriging 法對(duì)草地土壤水分、植被蓋度及物種豐富度指數(shù)進(jìn)行插值繪圖(圖5)?？梢钥闯?自潛在沙化草地至重度沙化草地,0—100 cm各土層土壤水分含量和植被蓋度的空間插值圖呈條帶狀和斑塊狀的梯度變化,物種豐富度指數(shù)的空間插值圖則呈較為明顯的斑塊分布；自潛在沙化草地至重度沙化草地,土壤水分含量、植被蓋度及物種豐富度指數(shù)總體表現(xiàn)為逐漸降低的趨勢(shì),且各土層土壤水分西部比東部變化較為平緩。隨土層加深,土壤水分的高值區(qū)和低值區(qū)呈先升高后降低,高值區(qū)由東向西北方和低值區(qū)由南向北變化。

從土壤水分含量的分異看,0—10 cm土層土壤水分克里格插值圖以斑塊狀居多,最高值出現(xiàn)在東側(cè),等值線彎曲,有明顯的“高地”和“低谷”,空間變異程度較高；其余土層土壤水分等值線較平直,過(guò)渡較為平滑,20—60 cm土層土壤水分總體表現(xiàn)為東北高西南低。

2.4 沙化草地各土層土壤水分、植被蓋度及物種豐富度指數(shù)的相互關(guān)系

相關(guān)分析表明(圖6),寧夏東部風(fēng)沙區(qū)沙化草地植被蓋度與0—40 cm各土層土壤水分含量呈顯著正相關(guān)(P<0.05),與40—100 cm各土層土壤水分含量呈極顯著正相關(guān)(P<0.01),相關(guān)系數(shù)為0.26—0.55；與物種豐富度指數(shù)相關(guān)不顯著(P>0.05)。物種豐富度指數(shù)與0—50 cm各土層土壤水分含量呈顯著正相關(guān)(P<0.05)。除0—10 cm土層與60—70 cm、70—80 cm土層,10—20 cm土層與70—80 cm土層之間土壤水分含量相關(guān)不顯著外(P>0.05),其余各土層之間土水分含量均呈顯著或極顯著正相關(guān)(P<0.05或P<0.01),且相鄰?fù)翆又g土壤水分含量相關(guān)系數(shù)較高。

圖6 寧夏東部風(fēng)沙區(qū)沙化草地0—100 cm各土層土壤水分、植被蓋度及物種豐富度指數(shù)的相關(guān)分析

3 討論

土地沙化對(duì)區(qū)域生態(tài)、水文和農(nóng)業(yè)等均造成嚴(yán)重危害,一直是生態(tài)環(huán)境治理的重點(diǎn)問(wèn)題之一[27]。水分是干旱半干旱地區(qū)草地植物生長(zhǎng)的重要制約因素[28]。寧夏東部風(fēng)沙區(qū)沙化草地0—100 cm各土層土壤水分含量總體處于較低水平。這是由于當(dāng)?shù)亟邓偾易兟蚀?加之土質(zhì)疏松,強(qiáng)烈的風(fēng)蝕加速土壤蒸發(fā),土壤保水能力較差[10,29]。植被蓋度自潛在沙化草地至重度沙化草地逐漸降低。在草地沙化過(guò)程中,伴隨著土壤細(xì)顆粒的吹蝕,土壤結(jié)構(gòu)變得更為松散,土壤養(yǎng)分損失,保水能力減弱,限制了植物的生長(zhǎng)和發(fā)育,植被發(fā)生逆向演替,群落向著更為耐旱耐貧瘠的方向演替,植被蓋度、密度及生物量等降低[30—32]。隨草地沙化程度的加重,植被蓋度、物種豐富度指數(shù)及各土層土壤含水量總體呈下降趨勢(shì),與以往的研究結(jié)果一致[32—33]。

植被及其生存的立地環(huán)境在任何尺度上都是時(shí)空連續(xù)的變異體,具有一定的空間特征[34]。寧夏東部風(fēng)沙區(qū)沙化草地0—100 cm各土層土壤水分含量的空間變異性表現(xiàn)為中等變異,與以往研究得出的結(jié)果相似[6,10,30]。隨土層的加深,土壤水分含量呈先升高后降低的趨勢(shì),以60—90 cm各土層較高,且表層土壤水分含量的變異系數(shù)較大,可能是研究區(qū)地處毛烏素沙地南緣,強(qiáng)烈的蒸發(fā)、風(fēng)蝕及水分入滲作用導(dǎo)致表層土壤水分含量較低,并使得空間異質(zhì)性增加[35]。隨土層深度的增加,影響土壤水分運(yùn)移的作用力發(fā)生改變,土壤水分含量增加[36]。

空間異質(zhì)性是隨機(jī)性因素和結(jié)構(gòu)性因素共同作用的結(jié)果,塊金系數(shù)(C0/C0+C)可表明隨機(jī)因素引起的空間變異占總變異的比例。寧夏東部風(fēng)沙區(qū)沙化草地物種豐富度指數(shù)及各土層土壤水分含量的塊金系數(shù)為26.31%—41.86%,表現(xiàn)為中等的空間自相關(guān)性,表明由地形地貌、成土母質(zhì)及植被等自然要素引起的結(jié)構(gòu)性因素和人類干擾等要素引起的隨機(jī)因素在其空間總變異中均發(fā)揮作用[26]。表層土壤水分較其他土層的塊金系數(shù)大,強(qiáng)烈的風(fēng)蝕通過(guò)影響水分循環(huán)及植被覆蓋對(duì)降雨和表層地下水的截留作用導(dǎo)致表層水分含量的變化更趨于不穩(wěn)定[37]。植被蓋度的塊金系數(shù)為2.50%,小于25%,表現(xiàn)出強(qiáng)烈的空間自相關(guān)性。植物是土壤中有機(jī)質(zhì)及各種營(yíng)養(yǎng)元素的主要來(lái)源,同時(shí),土壤作為植物生長(zhǎng)最重要的環(huán)境因子,為植物的生長(zhǎng)提供機(jī)械支持和物質(zhì)基礎(chǔ)[34]。植被蓋度、物種豐富度指數(shù)及各土層土壤水分含量的變程為1453—8110 m,說(shuō)明其生態(tài)學(xué)過(guò)程在較大的尺度上發(fā)揮作用,連續(xù)性較好。

寧夏東部風(fēng)沙區(qū)沙化草地各土層土壤水分和植被蓋度的空間插值圖呈條帶狀和斑塊狀的梯度變化。任何生態(tài)系統(tǒng)的發(fā)展都具有雙向性,在自然因素和人為因素的共同作用下,可能發(fā)生正向或逆向演替,研究區(qū)由于生境的脆弱、敏感性,在內(nèi)外因素的擾動(dòng)下更趨不穩(wěn)定,因此植被蓋度和不同土層土壤水分含量自潛在沙化草地至重度沙化草地逐漸降低的條帶狀梯度變化過(guò)程中出現(xiàn)了局部升高的斑塊格局。物種豐富度指數(shù)則呈現(xiàn)較為明顯的斑塊格局,高值區(qū)出現(xiàn)在潛在沙化草地至輕度沙化草地、輕度沙化草地至中度沙化草地之間的過(guò)渡區(qū)域,體現(xiàn)了不同生態(tài)類型過(guò)渡和轉(zhuǎn)換過(guò)程中的邊緣效應(yīng)；同時(shí),也表明草地沙化演替過(guò)程中,植被物種發(fā)生替代,一些抗風(fēng)蝕、耐沙埋的灌木、半灌木逐步占據(jù)優(yōu)勢(shì)地位[38—40]。土壤水分的空間分布受地形、土壤結(jié)構(gòu)、質(zhì)地及局部微氣候等環(huán)境因子的直接影響,同時(shí),由于植物與土壤之間的互饋?zhàn)饔?植物個(gè)體發(fā)育、種群及群落結(jié)構(gòu)等生物因子也對(duì)土壤水分產(chǎn)生影響[34,41]。伴隨著草地的沙化,由風(fēng)蝕作用產(chǎn)生的土壤粗?；瘜?dǎo)致土壤質(zhì)地變差、結(jié)構(gòu)更為松散,加之流沙堆積,局地微地形發(fā)生變化,引起土壤水分的再分配[29]。其次,自潛在沙化草地至重度沙化草地,群落發(fā)生逆向演替,結(jié)構(gòu)趨于簡(jiǎn)單,優(yōu)勢(shì)物種更替,群落蓋度、密度及生物量等的變化導(dǎo)致土壤養(yǎng)分輸入、微生物活動(dòng)、水分蒸發(fā)蒸騰作用及土壤水分的運(yùn)移過(guò)程發(fā)生改變[33]。此外,沙化草地生態(tài)系統(tǒng)中灌木、半灌木的出現(xiàn),致使灌叢周圍在風(fēng)蝕作用下形成土壤水分、養(yǎng)分等的富集,呈現(xiàn)明顯的“沃島效應(yīng)”,進(jìn)而引起土壤水分的分異[42—43]。

4 結(jié)論

寧夏東部風(fēng)沙區(qū)沙化草地土壤水分含量總體處于較低水平,在0—100 cm土層內(nèi)隨剖面深度的增加呈先上升后下降的趨勢(shì),以70—90 cm土層最高。各土層土壤水分含量、植被蓋度及物種豐富度指數(shù)均屬于中等變異,受地形地貌、土壤母質(zhì)等結(jié)構(gòu)性因素和人類擾動(dòng)等隨機(jī)性因素的共同作用,其空間分布呈現(xiàn)出條帶狀和斑塊狀的梯度變化格局,且連續(xù)性較好,自潛在沙化草地至重度沙化草地,總體表現(xiàn)為逐漸降低的趨勢(shì)?；诟鲄?shù)的空間分布格局,在草地沙化演替過(guò)程中不同沙化草地類型之間的過(guò)渡和轉(zhuǎn)換區(qū)是土壤水分和植被變化的敏感區(qū)域,應(yīng)作為草地沙化防治的關(guān)鍵區(qū)域加強(qiáng)管理。

致謝：試驗(yàn)過(guò)程中得到寧夏鹽池縣科技局及彭文棟、趙永升、李永康等的大力支持和幫助。