湯木子

(四川大學(xué) 建筑與環(huán)境學(xué)院, 四川 成都 610065)

陸地及水域等構(gòu)成了紛繁復(fù)雜的生態(tài)系統(tǒng)，其中又包含數(shù)量眾多的生態(tài)子系統(tǒng)。濕地是地球表層中由水陸生態(tài)系統(tǒng)相互作用形成的過(guò)渡性自然綜合體，具有獨(dú)特的植物群落結(jié)構(gòu)特征和較高的生物多樣性，是生態(tài)系統(tǒng)中至關(guān)重要的組成部分，為動(dòng)植物提供了豐富的資源和環(huán)境[1-3]。從生態(tài)環(huán)境的角度來(lái)講，由于濕地的構(gòu)成要素眾多且存在顯著的交互作用，因此其復(fù)雜性和多樣性尤為突出，對(duì)于水陸生態(tài)而言更是如此，作為人類生存發(fā)展的載體，動(dòng)植物的生長(zhǎng)發(fā)育均受其直接制約，而濕地植物群落的結(jié)構(gòu)與特征是濕地土壤與微生物因子的綜合效應(yīng)[4-5]。立足于濕地生態(tài)環(huán)境，其中關(guān)鍵的構(gòu)成部分離不開土壤和微生物，絕大多數(shù)植被的生長(zhǎng)發(fā)育離不開土壤，在物質(zhì)和能量交換過(guò)程中離不開微生物，從而形成了土壤—植物局地生態(tài)系統(tǒng)[6]。作為土壤生態(tài)的重要構(gòu)成部分，微生物對(duì)于凋落物的分解發(fā)揮著無(wú)可替代的效果，同時(shí)在酶的參與下進(jìn)一步提升了分解效果，二者是息息相關(guān)的，此外，環(huán)境因子也發(fā)揮著明顯的制約效應(yīng)[7-9]。土壤微生物在凋落物等腐殖質(zhì)的降解轉(zhuǎn)化過(guò)程中尤為關(guān)鍵，直接制約著土壤肥力及活性，因此，在環(huán)境因子中微生物常常作為關(guān)鍵指標(biāo)對(duì)生態(tài)系統(tǒng)起著決定性作用，土壤與微生物及植被之間的關(guān)系具有重要的現(xiàn)實(shí)意義[10-11]。

若爾蓋濕地屬于典型的高原沼澤濕地，且分布面積較廣，由于地處青藏東緣，地理位置的特殊性也造就了其特殊地位：不僅是黃河源頭區(qū)域，而且還是長(zhǎng)江源頭區(qū)域，成為典型的“高原之腎”；由于其蓄水效應(yīng)等顯著，也具有“黃河蓄水池”的美譽(yù)，在水質(zhì)調(diào)節(jié)方面具有重要作用。由于面積分布廣、地理位置特殊，生態(tài)分布的多樣性顯著，在局地生態(tài)碳循環(huán)方面效果顯著，同時(shí)具有重要的經(jīng)濟(jì)價(jià)值[12-13]。對(duì)于該區(qū)域而言，由于獨(dú)特的生態(tài)環(huán)境，植被群落結(jié)構(gòu)復(fù)雜，物種較為豐富[14-15]，對(duì)于該區(qū)域碳平衡具有關(guān)鍵的調(diào)節(jié)作用；對(duì)于濕地的植物群落結(jié)構(gòu)而言，通過(guò)連續(xù)的碳輸入影響著土壤碳的固存，而土壤微生物在這過(guò)程中發(fā)揮著重要作用，主要是通過(guò)呼吸作用來(lái)影響土壤碳的固存。

隨著工業(yè)進(jìn)程的不斷深入，人類活動(dòng)干預(yù)也日漸突出，對(duì)該區(qū)域的生態(tài)環(huán)境產(chǎn)生了較大的影響：①造成了濕地分布面積呈大幅的下降；②濕地水質(zhì)質(zhì)量明顯下降，原有生態(tài)環(huán)境和棲息地被大幅改變，生物多樣性逐年降低，植被生產(chǎn)力也明顯降低，導(dǎo)致各生態(tài)系統(tǒng)服務(wù)功能衰減[16]。由于若爾蓋濕地區(qū)域地處兩大流域的源頭，其重要地位尤為顯著，植被多樣性穩(wěn)定直接關(guān)乎整個(gè)濕地的多樣性?；诖?，本文2016—2019年研究了若爾蓋高原濕地不同植物群落結(jié)構(gòu)特征與土壤微生物群落多樣性的長(zhǎng)期變化特征及其二者之間的演化關(guān)系，以期為該區(qū)的植被恢復(fù)和生態(tài)環(huán)境保護(hù)提供有益的借鑒。

1 研究區(qū)概況

若爾蓋縣位于四川省西北部，地處青藏高原東北部的若爾蓋高原上，地理位置介于東經(jīng)102°08′—103°39′，北緯32°56′—34°19′之間。在行政上與瑪曲縣、碌曲縣、若爾蓋縣、阿壩縣和紅原縣相鄰。試驗(yàn)地點(diǎn)位于四川阿壩若爾蓋縣濕地國(guó)家級(jí)自然保護(hù)區(qū)管理局內(nèi)。本研究選取若爾蓋高原濕地不同植物群落，從2016年開始進(jìn)行連續(xù)4 a的定位觀測(cè)。該區(qū)域面積較廣，具有典型的高原沼澤濕地特征，其緊鄰青藏高原的獨(dú)特地理分布，屬于3 000 m以上的高海拔濕地，擁有獨(dú)特的濕地生態(tài)特征，在中國(guó)濕地分布中顯得尤為獨(dú)特。受制于所在地理分布，大陸性高寒季風(fēng)氣候特點(diǎn)突出，導(dǎo)致晝夜溫差較大，結(jié)合近年來(lái)該區(qū)域氣象資料得知，最低溫度達(dá)到了-33 ℃，1月溫度相對(duì)較低，約為-10 ℃，在冬季酷寒天氣影響下形成了明顯的凍土層；即使是進(jìn)入溫度較高的夏季，其均溫也僅為10 ℃左右，最高氣溫約為25 ℃。該區(qū)域降雨量并不高，僅為600 mm左右，全年日照時(shí)長(zhǎng)約為2 400 h，植被分布方面優(yōu)勢(shì)物種主要有木里苔草(Carexmuliensis)、烏拉苔草(Carexmeyeriana)，藏嵩草(Kobresiatibetica)、花葶驢蹄草(Calthascaposa)等。

2 材料與方法

2.1 試驗(yàn)方法

本試驗(yàn)選取的濕地植被包括木里苔草、烏拉苔草、藏嵩草、花葶驢蹄草4種植物群落。為了最大程度地降低試驗(yàn)誤差，對(duì)各個(gè)群落設(shè)置3個(gè)樣地，要求長(zhǎng)、寬均為20 m，在各個(gè)樣地設(shè)置樣方15個(gè)，長(zhǎng)、寬均為1 m。正式試驗(yàn)時(shí)間為2016—2019年，連續(xù)4 a定位觀測(cè)。土壤的采樣于每年的8月進(jìn)行，通過(guò)5點(diǎn)混合法，將樣品進(jìn)行均勻混合后一部分在-80 ℃下保存，剩余部分去除雜質(zhì)后進(jìn)行風(fēng)干處理，而后進(jìn)行0.25 mm篩過(guò)篩處理，選取其中50 g，磨細(xì)后過(guò)0.05 mm篩，然后將其用聚乙烯袋保存，以備進(jìn)行指標(biāo)測(cè)定，為了降低試驗(yàn)誤差，進(jìn)行3次重復(fù)測(cè)定。

α多樣性計(jì)算如下[17]：

Margalef豐富度(S)：S=(N-1)/lnN

(1)

Wiener多樣性(H)：H=-∑(PilnPi)

(2)

Simpson優(yōu)勢(shì)度(D)：D=1-∑(Pi)2

(3)

Pielou均勻度(JP)： JP=H/lnS

(4)

β多樣性計(jì)算如下[15]：

Cody指數(shù)=(a+b)/2

(5)

Sorenson指數(shù)=2j/(a+b)

(6)

式中：N表示樣方中物種數(shù)目；Pi表示重要值，計(jì)算方法為：(相對(duì)蓋度+相對(duì)多度+相對(duì)高度)/3；a表示樣方1的物種數(shù)目；b表示樣方2的物種數(shù)目；j表示樣方1和2的共有物種數(shù)目。

2.2 土壤微生物群落功能多樣性

土樣養(yǎng)分、微生物量和微生物數(shù)目參考《土壤農(nóng)化分析》[18]。

Biolog-Eco板中的碳源共有31種，為了對(duì)微生物的利用水平進(jìn)行較準(zhǔn)確的探究，本研究借助于平均顏色變化率(AWCD)這一研究方式。AWCD值越大，說(shuō)明土壤微生物對(duì)碳源的利用水平越高，且微生物具有更大的豐富度，該值的計(jì)算公式如下[18]：

AWCD=∑〔(Ci-R)/31〕

(7)

式中：第i孔的吸光值用Ci代表,而對(duì)照孔則用R表示;若(Ci-R)<0,則計(jì)算過(guò)程中將之記為0,即(Ci-R)≥0。

土壤微生物作為土壤肥力保持的重要參與者，其群落多樣性的測(cè)量常通過(guò)以下指數(shù)測(cè)定實(shí)現(xiàn)：Shannon多樣性指數(shù)(H)，McIntosh(Pi)和Simpson(Ds)指數(shù)[18]。

AWCD=∑(Ci-Ri)/n

(8)

H=-∑Pi(lnPi)

(9)

Pi=(Ci-Ri)/∑(Ci-Ri)

(10)

S=被利用碳源的總數(shù)

(當(dāng)吸光值≥0.25的情況下微孔總和)

D=1-∑Pi

(11)

式中：Ci為第i孔的吸光值；R為對(duì)照孔的吸光值。對(duì)(Ci-R)<0的孔，計(jì)算中記為0，即(Ci-R)≥0。Pi為第i孔的相對(duì)吸光值與所有整個(gè)微平板的相對(duì)吸光值總和的比值。

2.3 統(tǒng)計(jì)分析

采用SPSS 21和Excel 2010進(jìn)行數(shù)據(jù)整理和分析，其中所有的值均為平均值±標(biāo)準(zhǔn)偏差表示，采用單因素方差分析(LSD方法)比較不同年份各個(gè)指標(biāo)之間的差異顯著性(p<0.05)，用Pearson相關(guān)系數(shù)法檢驗(yàn)各個(gè)指標(biāo)之間的相關(guān)性。采用CANACO 5.0進(jìn)行排序分析。

3 結(jié)果與分析

3.1 若爾蓋高原濕地植物群落多樣性特征

由表1可知，對(duì)于α多樣性，2016—2019年若爾蓋高原濕地植被Shannon-Wiener多樣性指數(shù)、Margalef豐富度指數(shù)、Pielou均勻度指數(shù)變化趨勢(shì)一致，均表現(xiàn)為：烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草，隨著年份的增加而逐漸增加趨勢(shì)；Simpson優(yōu)勢(shì)度指數(shù)表現(xiàn)為：烏拉苔草<木里苔草<藏嵩草<花葶驢蹄草，呈現(xiàn)出隨著年份的增加而逐漸增加趨勢(shì)，不同年份而Shannon-Wiener多樣性指數(shù)差異不顯著(p>0.05)。對(duì)于β多樣性，2016—2019年若爾蓋高原濕地植被Cody指數(shù)呈一致的變化規(guī)律，均表現(xiàn)為：烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草，呈現(xiàn)出隨著年份的增加而逐漸增加趨勢(shì)；而Sorenson指數(shù)呈出相反的變化趨勢(shì)，隨著年份的增加而逐漸增加趨勢(shì)，不同年份Sorenson指數(shù)差異不顯著(p>0.05)。

表1 2016-2019年若爾蓋高原濕地植物群落多樣性

3.2 若爾蓋高原濕地不同植物群落土壤養(yǎng)分

由表2可知，2016—2019年若爾蓋高原濕地植被土壤有機(jī)碳、全氮、全磷和全鉀含量呈一致的變化規(guī)律，均表現(xiàn)為：烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草，隨著年份的增加而逐漸增加趨勢(shì)。2018年和2019年土壤有機(jī)碳和全氮含量差異不顯著(p>0.05)，二者顯著高于2016年和2017年(p<0.05)。土壤全磷在不同年份差異均不顯著(p>0.05)。其中土壤有機(jī)碳變化范圍在8.15～17.01 g/kg之間，全氮變化范圍在1.25～1.83 g/kg之間，全磷變化范圍在0.69～0.89 g/kg之間，全鉀變化范圍在31.09～41.23 g/kg之間。

表2 2016-2019年若爾蓋高原濕地不同植物群落土壤養(yǎng)分

3.3 若爾蓋高原濕地不同植物群落土壤微生物量

由表3可知，若爾蓋高原濕地不同植物群落土壤微生物量存在一定差異。2016—2019年土壤微生物量碳變化范圍在231.9～346.2 mg/kg之間，整體表現(xiàn)為：烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草，隨著年份的增加呈增加趨勢(shì)，不同年份土壤微生物量碳含量差異均顯著(p<0.05)，而不同年份土壤微生物量磷含量差異不顯著(p>0.05)。土壤微生物量氮變化范圍在65.3～101.5 mg/kg之間，整體表現(xiàn)為:烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草。微生物量磷變化范圍在15.2～26.2mg/kg之間，整體表現(xiàn)為:烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草。不同年份之間若爾蓋高原濕地不同植物群落土壤微生物量差異不顯著(p>0.05)。

表3 2016-2019年若爾蓋高原濕地不同植物群落土壤微生物量

3.4 若爾蓋高原濕地不同植物群落土壤微生物群落結(jié)構(gòu)

由表4可知，若爾蓋高原濕地不同植物群落土壤細(xì)菌數(shù)量占絕對(duì)優(yōu)勢(shì)，占微生物總數(shù)的70%左右，其次是放線菌，真菌最少。2016—2019年若爾蓋高原濕地植物群落土壤細(xì)菌、真菌和放線菌數(shù)量呈一致的變化規(guī)律，均表現(xiàn)為:烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草，隨著年份的增加而逐漸增加趨勢(shì)；并且不同年份土壤細(xì)菌數(shù)量、真菌數(shù)量和放線菌數(shù)量差異均顯著(p<0.05)。

表4 2016-2019年若爾蓋高原濕地不同植物群落土壤微生物數(shù)量

由表5可知，若爾蓋高原濕地不同植物群落土壤微生物群落功能多樣性指數(shù)存在一定差異，2016—2019年土壤微生物物種豐富度指數(shù)(H)變化范圍在1.78～4.05之間，整體表現(xiàn)為：烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草，隨年份的增加呈增加趨勢(shì)；土壤微生物均勻度指數(shù)(E)變化范圍在0.31～0.81之間，整體表現(xiàn)為：烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草，隨年份的增加呈增加趨勢(shì)，并且不同年份土壤微生物均勻度指數(shù)(E)差異不顯著(p>0.05)；優(yōu)勢(shì)度指數(shù)(D)變化范圍變化范圍在0.42～0.62之間，整體表現(xiàn)為：烏拉苔草<木里苔草<藏嵩草<花葶驢蹄草，隨年份的增加呈增加趨勢(shì)，并且不同年份土壤微生物優(yōu)勢(shì)度指數(shù)(D)差異不顯著(p>0.05)；碳源利用豐富度指數(shù)(S)變化范圍在5.77～16.43之間，整體表現(xiàn)為：烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草，隨年份的增加呈增加趨勢(shì)。

表5 2016-2019年若爾蓋高原濕地不同植物群落土壤微生物群落多樣性

3.5 植物群落多樣性與土壤微生物之間的相關(guān)分析

由表6可知，Margalef指數(shù)和Shannon-Wiener指數(shù)與細(xì)菌數(shù)量、微生物量碳、微生物量氮、物種豐富度指數(shù)(H)、碳源利用豐富度指數(shù)(S)、優(yōu)勢(shì)度指數(shù)(D)和均勻度指數(shù)(E)呈顯著正相關(guān)(p<0.05)；Simpson指數(shù)與細(xì)菌數(shù)量、微生物量碳、微生物量氮、物種豐富度指數(shù)(H)、碳源利用豐富度指數(shù)(S)和均勻度指數(shù)(E)呈顯著負(fù)相關(guān)(p<0.05)；Pielou指數(shù)與細(xì)菌數(shù)量、微生物量碳、微生物量氮、物種豐富度指數(shù)(H)、碳源利用豐富度指數(shù)(S)呈顯著正相關(guān)(p<0.05)；Cody指數(shù)與細(xì)菌數(shù)量、微生物量碳、微生物量氮、物種豐富度指數(shù)(H)、碳源利用豐富度指數(shù)(S)、優(yōu)勢(shì)度指數(shù)(D)和均勻度指數(shù)(E)呈顯著正相關(guān)(p<0.05)；Sorenson指數(shù)細(xì)菌數(shù)量、微生物量碳、微生物量氮、物種豐富度指數(shù)(H)、碳源利用豐富度指數(shù)(S)、優(yōu)勢(shì)度指數(shù)(D)呈顯著負(fù)相關(guān)(p<0.05)。

表6 若爾蓋高原濕地植物群落多樣性與土壤微生物之間的相關(guān)分析

此外，采用冗余分析RDA研究土壤養(yǎng)分、土壤微生物數(shù)量和微生物群落多樣性的排序關(guān)系，將植物群落多樣性作為響應(yīng)變量，土壤養(yǎng)分、土壤微生物數(shù)量和微生物群落多樣性為解釋變量，通過(guò)排序分析更直觀地展現(xiàn)不同的環(huán)境因子所產(chǎn)生的不同影響(圖1)。其中箭頭象限代表著相關(guān)性，箭頭長(zhǎng)度代表相關(guān)程度，其夾角越大則說(shuō)明其相關(guān)性越弱；通過(guò)對(duì)比發(fā)現(xiàn)，前兩個(gè)排序軸的累積解釋率達(dá)到了55.45%。因此，土壤養(yǎng)分、土壤微生物數(shù)量和微生物群落多樣性很好地解釋了植物群落多樣性的變化。此外，通過(guò)研究發(fā)現(xiàn)，土壤微生物多樣性越豐富，植物多樣性越豐富，二者之間存在明顯的正相關(guān)關(guān)系。

注：圖中實(shí)心箭頭和開口箭頭分別代表響應(yīng)變量(植物多樣性)和解釋變量(土壤養(yǎng)分、土壤微生物數(shù)量和生物群落多樣性)。

4 討論與結(jié)論

4.1 討 論

對(duì)于濕地不同的植物群落，不僅受到土壤養(yǎng)分等狀況的直接影響，還受到土壤微生物等調(diào)控因子的影響，主要是由于土壤有機(jī)質(zhì)的合成與分解離不開微生物的代謝，土壤微生物在土壤養(yǎng)分循環(huán)和能量交換過(guò)程中發(fā)揮著重要作用，再加上土壤微生物自身處于不斷的生長(zhǎng)和繁殖動(dòng)態(tài)中，因此，土壤微生物對(duì)于植被群落結(jié)構(gòu)的影響是處于不斷變化中[19-23]。此外，環(huán)境因子(水、熱條件等)也影響著植物群落多樣性和結(jié)構(gòu)分布特征。在土壤微生物的作用下，土壤理化特性發(fā)生改變，相應(yīng)地影響著土壤微生物新陳代謝活動(dòng)，影響其土壤有機(jī)質(zhì)降解效率[24-27]。通過(guò)連續(xù)4 a的觀測(cè)分析得知，對(duì)于α多樣性，Shannon多樣性、Margalef豐度、Pielou均勻度指數(shù)，雖然不同植物群落表現(xiàn)的有所差異，但是整體而言，變化趨勢(shì)相一致，均表現(xiàn)為：在年份不斷增加的情況下呈現(xiàn)上升趨勢(shì)，最高的是烏拉苔草，其次是木里苔草，而花葶驢蹄草最低；對(duì)于Simpson指數(shù)，呈相反的變化規(guī)律。對(duì)于β多樣性，在年份不斷增加的情況下，Cody指數(shù)呈現(xiàn)上升趨勢(shì)，且最高的是烏拉苔草，其次是木里苔草，而花葶驢蹄草最低；而Sorenson指數(shù)呈相反的變化趨勢(shì)[28-29]。綜合來(lái)看，若爾蓋濕地植被多樣性隨著年份的增加呈現(xiàn)整體的上升趨勢(shì)，尤其是烏拉苔草群落，具有較穩(wěn)定的生態(tài)群落分布；從Sorenson指數(shù)的角度來(lái)說(shuō)，該指數(shù)的下降說(shuō)明物種不斷增多，這也是自然環(huán)境及人為影響綜合作用的結(jié)果。

若爾蓋高原濕地不同植物群落土壤細(xì)菌數(shù)量占絕對(duì)優(yōu)勢(shì)，占微生物總數(shù)的70%左右，其次是放線菌，真菌最少；其中細(xì)菌數(shù)目、真菌數(shù)目和放線菌數(shù)目均表現(xiàn)為：烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草，隨著年份的增加而逐漸增加趨勢(shì)。對(duì)于烏拉苔草群落，大量的凋落物形成了大量的腐殖質(zhì)來(lái)源，在微生物作用下，其被分解為有機(jī)質(zhì)，為土壤養(yǎng)分的積累奠定了基礎(chǔ)，同時(shí)，凋落物能夠有效促進(jìn)土壤水分的保持，利于涵養(yǎng)降雨等水分，從而在促進(jìn)自身生長(zhǎng)的同時(shí)提升微生物新陳代謝的條件[30-31]，最終微生物群落多樣性呈一致的變化趨勢(shì)。相關(guān)性分析表明，土壤微生物群落多樣性與植被多樣性具有顯著的相關(guān)性(p<0.05)由此說(shuō)明微環(huán)境因子中微生物群落多樣性顯著影響濕地植被多樣性。

4.2 結(jié) 論

若爾蓋高原濕地植物群落物種Shannon-Wiener多樣性指數(shù)、Margalef豐富度指數(shù)、Pielou均勻度指數(shù)呈一致的變化規(guī)律，均表現(xiàn)為：烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草，隨著年份的增加而逐漸增加趨勢(shì)；植物群落土壤養(yǎng)分含量(土壤有機(jī)碳、全氮、全磷和全鉀)、微生物量碳、氮、磷和土壤微生物群落多樣性指數(shù)基本表現(xiàn)出相同的規(guī)律(烏拉苔草>木里苔草>藏嵩草>花葶驢蹄草)，隨著年份的增加而逐漸增加趨勢(shì)，局部有所波動(dòng)。若爾蓋高原濕地植被群落物種多樣與土壤微生物群落多樣性的變化相一致，二者具有正相關(guān)關(guān)系，若爾蓋高原濕地地上植被與地下微生物群落多樣性之間強(qiáng)耦合關(guān)系，并且這種關(guān)系并不隨時(shí)間而改變，二者具有對(duì)環(huán)境變化響應(yīng)的統(tǒng)一性。土壤微生物環(huán)境因子改變著土壤養(yǎng)分和微生物特性，另一方面土壤環(huán)境為微生物活動(dòng)創(chuàng)造更為適宜的新陳代謝環(huán)境，二者具有明顯的相互促進(jìn)性，在濕地生態(tài)保護(hù)過(guò)程中要充分注重微生物和養(yǎng)分的協(xié)調(diào)，合理控制氮、磷等含量，發(fā)揮微生物效應(yīng)。