河西走廊綠洲–荒漠過渡帶濕生與鹽生生態型蘆葦種群空間分布格局分析

芝祥紅，獨肖艷，李朝周,3，焦 健，劉 鑫

（1.甘肅農業大學生命科學技術學院，甘肅蘭州730070；2.甘肅農業大學林學院，甘肅蘭州730070；3.甘肅省作物遺傳改良和種質創新重點實驗室，甘肅蘭州730070）

河西走廊地處巴丹吉林沙漠南緣，綠洲–荒漠過渡地帶是維持綠洲穩定和防止綠洲沙漠化的前沿陣地，過渡帶植被生長與健康狀況決定著綠洲的衰退與擴張，對綠洲的穩定性起決定性作用。近幾十年來在人類作用和區域氣候變化雙重驅動下，該區域地下水位逐年下降，土壤旱化、鹽漬化，植被生長受到限制[1]，急需進一步加強區域植被保護工作，維護生態安全。

蘆葦(Phragmites australis)是根莖型克隆禾草，作為河西走廊地區植被的主要建群種之一廣泛分布于綠洲農田外圍、鹽堿地及沙丘上，約占耕地總面積的50%[2]，蘆葦作為典型的克隆植物，萌生能力強，眾多根莖串聯交織，增加了下墊面的粗糙度，能降低風速，減輕風沙危害，有效遏制河西走廊沙漠化的擴張，對綠洲的穩定作用顯著。以往國內對蘆葦的研究多集中于黃河三角洲地區和松嫩平原蘆葦種群，研究內容主要涉及遺傳多樣性[3-4]、生理生化特征[5-6]、生態學特征[7-9]及種群空間分布格局[10-11]等，河西走廊地區蘆葦的研究多集中于遺傳基因[12-13]和生理生化特性[5,14]，研究表明，該區域蘆葦源于濕地生境，長期生長于鹽堿化或沙丘土壤中，其形態特征與各自生長環境相適應，產生了的極大變異，有鹽生蘆葦、濕生蘆葦和沙生蘆葦不同生態型的分化[14]，目前對不同生態型蘆葦種群空間格局的研究鮮見報道。

研究種群空間分布格局的方法很多，不同的空間分布格局研究方法是從不同的角度來度量種群的聚集特征。擴散系數(C)、Moristia指數(Iδ)均不受樣方大小的影響；平均擁擠度(m*)受到種群密度的影響，可以表征種群在空間分布即樣方之間的差異，還表征單位樣方內種群的數量和密集程度；聚塊性指數(C*)只考慮空間格局本身的性質，不受種群密度的影響，能作為衡量集群程度的指標，表征種群對空間的侵占、擴展程度，Green指數(IG)可以比較種群的聚集程度，體現種群聚集特征差異[15]。在聚集強度分析中，需要用不同參數對格局強度的計算結果進行對比分析，以檢驗結果的可靠性。

本研究在甘肅省臨澤縣綠洲–荒漠過渡帶設置樣地，選取有蘆葦種群天然分布的鹽堿、濕地設置樣地，采用具有相同本質的或較大相似性的擴散系數、負二項參數、Moristia指數、平均擁擠度、聚塊性指標、Green指數6種指標及Iwao回歸分析法，綜合分析在濕地和鹽堿生境生長的不同尺度及不同高度級蘆葦種群分布格局及聚集強度，為研究種群個體行為、種群擴散型及逆境適應提供一定信息，為綠洲生態系統的穩定性及植被恢復提供理論依據。

1 研究區概況

以位于巴丹吉林沙漠南緣的甘肅省臨澤縣為定位研究區。臨澤縣處于99°51′?100°30′E，38°57′?39°42′N，河西走廊中部，是典型的溫帶大陸性荒漠氣候，年均氣溫7.7℃，多年平均降水量118.4 mm，年蒸發量1 830.4 mm，晝夜溫差大，光照和水分配置極不均衡，土壤次生鹽漬化嚴重[16]。

2 研究方法

2.1 樣地調查與土壤理化性質測定

本研究分別選取有蘆葦天然分布的濕地和鹽堿生境設置樣地，其土壤理化性質如表1所列。

1)濕地生境(H1)：調查樣地設于臨澤縣平川鎮三一村黑河岸邊，該生境生長濕生蘆葦，主要伴生種為苜蓿(Medicago sativa)、野艾蒿(Artemisia lavandulaefolia)等。

2)鹽堿生境(H2)：位于臨澤縣鴨暖鎮馬營村鹽堿荒漠區，該生境生長鹽生蘆葦，主要伴生種為冰草(Agropyron cristatum)等。

2019年10月，在濕地和鹽堿生境，選取無顯著人為干擾、地形變化不大的地段設置20 m ×20 m的蘆葦種群永久性調查樣地，采用相鄰格子法將樣地分割成100個2 m ×2 m的小樣方，計數小樣方內蘆葦分株數目及其株高。

用土鉆(直徑2 cm)在各生境樣地進行土壤取樣，深度為1 m，5次重復。采集土壤樣品混勻后過2 mm篩備用，土壤含水量用烘干法測定；土壤含鹽量用質量法測定[17]，土水比1?5；pH用電位法測定；有機質含量采用重鉻酸鉀容量法測定；堿解氮含量采用堿解擴散法測定；速效磷含量采用0.5 mol·L?1NaHCO3浸提–鉬銻抗吸光光度法測定[18-19]。

表1 濕生和鹽堿生境土壤理化性質Table1 Soil physical and chemical characteristicsin wetland and salinehabitats

2.2 高度級劃分

根據蘆葦生長季末的調查結果，濕生與鹽生蘆葦株高差異明顯，本研究初步嘗試做了濕生與鹽生蘆葦高度級的劃分。濕生蘆葦高度級劃分為：Ⅰ，H< 140 cm；Ⅱ，140 cm ≤ H<210 cm；Ⅲ，210 cm ≤H< 280 cm；Ⅳ，280 cm ≤ H<350 cm；Ⅴ，350 cm ≤H<420 cm；Ⅵ，H ≥420 cm；鹽生蘆葦高度級劃分：Ⅰ，H<10 cm；Ⅱ，10 cm ≤H<15 cm；Ⅲ，15 cm ≤H<20 cm；Ⅳ，20 cm ≤H<25 cm；Ⅴ，25 cm ≤H<30 cm；Ⅵ，H ≥30 cm。

2.3 種群的分布格局測定

以2 m ×2 m的樣方為基本單位，按相鄰位置組合為4、8、16、20、40和100 m26個尺度，分別測定種群的分布格局。采用擴散系數(C)、負二項參數(K)、Moristia指數(Iδ)、平均擁擠度(m*)、聚塊性指數(C*)、Green指數(IG)6種聚集指標和Iwao模型來判定蘆葦種群的格局分布類型，分析不同取樣尺度及不同高度級蘆葦種群的分布格局及聚集強度。各指標常用計算方法如下[20-22]。

1)擴散系數(C)：

式中：v為樣本方差，m為樣本均值，當C< 1時為均勻分布；C=1時為隨機分布；C> 1時為聚集分布。

2)負二項參數(K)：

式中：K值用于度量聚集程度，K值越小，聚集程度越高，當K>8時逼近泊松分布。

3)Moristia指數(Iδ)：

式中：n為樣方數，N為個體總數，ni是第i個樣方的個體數，當Iδ<1時為均勻分布，Iδ= 1時為隨機分布，Iδ>1時為聚集分布。

4)平均擁擠度(m*)：

式中：m*反映了樣方內生物個體的擁擠程度，數值越大聚集強度越大。

5)聚塊性指數(C*)：

當C*<1時為均勻分布；C*= 1時為隨機分布；C*>1時為聚集分布。

6)Green指數(IG)：

當IG< 0時為均勻分布；IG= 0時為隨機分布；IG> 0時為聚集分布。

7)Iwao回歸分析法：

當α=0時，分布的基本成分為單個個體；當α>0時，分布的基本成分為個體群；當α< 0時，個體間相互排斥。α和β的聯合意義為：當α=0，β= 1時，為隨機分布；當α> 0，β= 1或α= 0，β> 1或α> 0，β>1時，呈聚集型分布；α<0，β> 1時，種群為均勻分布。在實際應用中α和β不會剛好為0和1，本研究采用下式對隨機偏離程度進行F檢驗[23-24]。

若F≥Fα(2，n?2)，則種群為聚集分布；反之則為隨機分布。

2.4 數據分析

用Excel 2007繪制圖表，SPSS 19.0軟件對數據進行分析處理。

3 結果與分析

3.1 濕生與鹽生蘆葦高度級結構

各高度級蘆葦存活個體數分析表明，濕地生境蘆葦種群數量隨高度級的增加呈現先上升后降低的趨勢，Ⅱ～Ⅳ高度級個體數占總個體數的88.63%，Ⅲ級存活個體數最多，達到4436株，占總個體數的39.92%；Ⅰ級和Ⅵ級個體數均比較少，分別為170和174株，占總個體數的1.53%和1.57%(圖1)。鹽堿生境蘆葦種數量隨高度級的增加同樣呈現先上升后降低的趨勢，Ⅲ級存活個體數最多，達到5103株，占總個體數的32.85%。Ⅰ級和Ⅵ級存活個體數均比較少，分別為762和744株，占總個體數的4.91%和4.79%(圖1)。兩生態型蘆葦種群高度級結構總體趨勢變化相同，且最低高度級與最高高度級占比均較小，兩者之和不足10%，存活個體數最多的高度級均為Ⅲ級。

圖1濕生與鹽生蘆葦種群高度級結構Figure 1 Height class structure of wet Phragmites australis and saline Phragmites australis populations

3.2 不同取樣尺度濕生與鹽生蘆葦種群空間分布格局聚集指數

采用6種聚集指數對濕生與鹽生蘆葦種群進行6個取樣尺度的空間分布格局分析，結果表明(表2)，濕生與鹽生蘆葦的C和C*值均大于1，IG值均大于0，表明兩種態型蘆葦在各取樣尺度均為聚集分布。Moristia指數(Iδ)及負二項參數(K)可以度量群落聚集強度。Iδ值越大聚集程度越高，但K值越大，聚集程度就越弱。以此為依據，濕生與鹽生蘆葦Iδ值隨取樣面積的增加呈減小趨勢，鹽生蘆葦K值隨取樣面積的增加呈增加趨勢，濕生蘆葦K值基本隨取樣面積的增加呈增加趨勢，但取樣面積超過20 m2K值略有不顯著降低，說明兩生態型蘆葦各自在不同取樣尺度間聚集程度相當。各取樣尺度鹽生蘆葦的Iδ值普遍小于濕生蘆葦，說明濕生蘆葦聚集程度大于鹽生蘆葦。鹽生蘆葦K值遠大于濕生蘆葦，且均大于8，說明各取樣尺度鹽生蘆葦聚集程度小于濕生蘆葦，且趨近泊松分布，有一定的隨機趨勢。

表2 不同取樣尺度濕生與鹽生蘆葦種群分布格局聚集指數Table2 Population distribution pattern aggregation index of wet Phragmites australis and saline Phragmites australis at different sampling scales

3.3 不同取樣尺度濕生與鹽生蘆葦種群空間分布格局的Iwao回歸模型分析

本研究以2 m ×2 m為基礎樣方，所以對蘆葦種群用Iwao方程進行回歸時，不對取樣面積為4 m2的種群建立回歸方程。濕生與鹽生蘆葦基于各取樣尺度得到的回歸方程(表3)，其相關系數R均大于0.8，說明回歸模型適合解釋該種群分布類型。各回歸方程中α值均大于1，表明分布的基本成分為個體群。取樣尺度為8和16 m2的濕生蘆葦和各取樣尺度的鹽生蘆葦回歸方程的β值均近似為1。應用Iwao模型隨機偏離程度進行檢驗，結果表明各取樣尺度濕生與鹽生蘆葦F值均大于F0.05(2，n?2)，即Iwao指數測定結果進一步說明濕生與鹽生蘆葦種群空間分布格局為聚集分布型。可見，用聚集度指標測定蘆葦種群空間分布格局的結果具有良好的一致性。

3.4 不同高度級濕生與鹽生蘆葦種群聚集指標分析

對濕生與鹽生蘆葦各高度級種群的空間格局進行綜合性分析(表4)，結果表明，濕生與鹽生蘆葦的C和C*值均大于1，I、Ca和IG值均大于0，表明兩種生態型蘆葦各高度級種群均為聚集分布。濕生與鹽生蘆葦Iδ值隨株高的增加呈先減小后增大的趨勢，濕生蘆葦Iδ值均大于鹽生蘆葦，且Iδ值在Ⅴ級和Ⅵ級達到高峰，分別為2.066和2.102，說明濕生蘆葦各高度級聚集程度大于鹽生蘆葦，尤其以Ⅴ和Ⅵ級聚集強度最大。濕生蘆葦K值隨株高的增加呈先增加后減小趨勢，鹽生蘆葦呈波動變化。鹽生蘆葦K值遠大于濕生蘆葦，且除Ⅰ和Ⅴ級外均大于8，由此可知，各高度級鹽生蘆葦聚集程度小于濕生蘆葦，高度級Ⅱ、Ⅲ、Ⅳ和Ⅵ級趨近泊松分布，有一定的隨機趨勢。

表3 不同取樣尺度濕生與鹽生蘆葦分布格局的Iwao方程模擬結果Table 3 Iwao equation simulation results of distribution patterns of wet Phragmites australis and saline Phragmites australis at different sampling scales

表4 不同高度級濕生與鹽生蘆葦種群分布格局聚集指數Table 4 Distribution pattern aggregation index of wet Phragmites australis and saline Phragmites australis populations at different height classes

4 討論

種群分布格局是物種與環境因素相互作用的結果，與物種的生境(包括土壤、地形、氣候、水分和養分等)有密切的聯系[20,25]。在氣候一致的區域內，微環境會影響小尺度的植被格局[26-27]。種群聚集分布體現了種群內部正向(相互有利)的生態關系，均勻分布體現了種群內部負向(相互排斥)的生態關系，隨機分布則表明種群內部無明確的生態關系[28]。有研究表明，自然種群在生態脆弱地帶多服從聚集分布，而隨機分布較為少見[29-30]。蘆葦具有發達的根狀莖，水平擴展能力極強，河西走廊地區蘆葦生長環境較為惡劣，該地區蘆葦通過營養繁殖進行種群更新，擴展種群的生態位空間，經過長時間的演化，形成聚集分布的空間格局，來相互庇護抵御不利環境。此外，種群分布格局與物種的生態生物學特性和種間的競爭排斥密切相關[31]。本研究在沙漠–綠洲過渡地帶，氣候條件一致的濕生和鹽生生境進行蘆葦種群的研究，結果表明，濕生與鹽生蘆葦在各取樣尺度上均為集群分布，其中鹽生蘆葦有一定的隨機分布趨勢。這與張瑩花等[32]和趙艷云等[11]對蘆葦種群的研究結果相吻合。張瑩花等[32]對石羊河中游河岸蘆葦空間分布格局的研究表明，蘆葦種群表現出集群分布特征，且距離河岸近的蘆葦聚集程度大于距離河岸遠的蘆葦。趙艷云等[11]對渤海灣貝殼堤濕地蘆葦種群與蒙古蒿(Artemisia mongolica)種群空間分布格局的研究表明，在小尺度上，蘆葦種群和蒙古蒿種群聚集分布尤為明顯。

蘆葦種群空間分布格局對土壤水分的響應要強于對土壤鹽分的響應。蘆葦起源于濕地，朱學藝等[33]認為，河西走廊區域的蘆葦演化規律為濕地水位下降，濕生蘆葦登陸形成輕度鹽化草甸蘆葦，水位繼續下降，形成高度鹽漬化的次生裸地，于是鹽堿蘆葦出現。貢璐等[34]對塔里木盆地南緣旱生蘆葦的研究表明，水分對蘆葦的影響大于鹽分，水分是蘆葦生態特征變化的關鍵驅動因子。夏貴菊等[35]對銀川平原蘆葦的研究表明，土壤含水量是影響銀川平原蘆葦生態特征的最關鍵因子，其次是pH和全鹽含量。焦德志等[36]對扎龍濕地蘆葦的研究顯示，影響蘆葦分株生態可塑性變異的主要驅動因子是土壤含水量，土壤鹽度的增加對蘆葦的密度、高度、莖粗、生物量以及根莖養分儲備都會產生負效益。張爽等[37]用不同鹽度水灌溉對蘆葦生長的研究表明，隨著鹽濃度增加，蘆葦株高、莖粗、鮮重、干重、纖維素含量均表現為降低趨勢。李嘗君等[38]對艾比湖濕地自然保護區克隆植物群落的研究說明花花柴(Karelinia caspiacaspica)和羅布麻(Apocynum venetum)更傾向于通過高度的集群分布來適應水鹽脅迫環境。王梅等[39]對鹽穗木(Halostachys caspica)種群空間格局的研究表明，隨著地下水埋深的增加，種群聚集強度增加。結合以上幾位學者的研究，本研究認為，濕生生境水分不再是蘆葦生長的限制因子，種群聚集度比較高；鹽生生境由于水鹽脅迫，水分成為蘆葦種群發展的限制因子，蘆葦為了滿足自身對水分、養分的需求，有隨機分布的傾向，進而獲得高大的空間及資源來維持其生存和發展，土壤水分的影響較土壤鹽分更利于種群聚集強度。

研究表明，濕生與鹽生蘆葦各自在不同取樣尺度間聚集程度相當。這與程煜等[40]對檫樹(Sassafras tsumu)種群在不同取樣尺度的分布格局研究結果相似。而冷艷芝等[41]對馬鹿河紅椿(Toona ciliata)種群分布格局的研究表明，隨取樣面積加大，紅椿種群分布的聚集性逐漸減弱。張瑩花等[32]對石羊河中游河岸蘆葦群落空間格局研究發現，距離河岸遠的蘆葦在小尺度上為集群分布，較大尺度上為隨機分布。可見，植物種群分布格局并不是一成不變的，而是會依據各自的生物學特性隨尺度變化而導致的環境異質性產生動態響應。

目前，關于種群空間分布的研究大多基于距離進行分析，對于蘆葦來說，株高是一個可以反映植株生長狀況的指標，因此本研究基于株高分析蘆葦的空間分布格局。結果表明，各高度級濕生蘆葦聚集程度大于鹽生蘆葦，尤其以Ⅴ和Ⅵ級聚集強度最大。濕生蘆葦處于相對適宜生境，故分株有足夠的營養物質來完成高生長，本次取樣時間處于蘆葦生長季后期，該時期蘆葦分株基本已發育完全，故累積于Ⅴ和Ⅵ級，產生其聚集強度最大的現象。鹽生蘆葦Ⅱ、Ⅲ、Ⅳ、Ⅵ級有一定的隨機趨勢，這與不同取樣尺度鹽生蘆葦研究的結論吻合，而Ⅰ、Ⅴ級沒有隨機趨勢。影響株高的因素非常復雜，本研究中發現土壤含水量、含鹽量、有機質含量等理化性質均與株高有顯著性相關關系，此外，其生物學特性和伴生種群的影響也不可忽視，有待于今后做更深層次的研究分析。

5 結論

濕生與鹽生蘆葦在各取樣尺度上均為集群分布，其中鹽生蘆葦有一定隨機分布的趨勢。濕生與鹽生蘆葦Moristia指數(Iδ)值隨取樣面積的增加呈減小趨勢，負二項參數(K)值隨取樣面積的增加呈增加趨勢，但差異并不大，說明兩生態型蘆葦各自在不同取樣尺度間聚集程度相當。植物種群分布格局并不是一成不變的，而是會依據各自的生物學特性隨尺度變化而導致的環境異質性產生動態響應。濕生蘆葦聚集強度大于鹽生蘆葦，說明土壤水分的較土壤鹽分更利于增加種群聚集強度。濕生與鹽生蘆葦高度級結構基本一致，隨高度級的增加呈現先上升后降低的趨勢，均為Ⅲ級個體最多，壯齡階段的蘆葦始終占絕對比重，說明齡級結構比較穩定。各高度級濕生蘆葦聚集強度大于鹽生蘆葦，尤其以Ⅴ和Ⅵ級濕生蘆葦聚集強度最大。鹽生蘆葦Ⅱ、Ⅲ、Ⅳ、Ⅵ級種群分布格局有一定的隨機趨勢，而Ⅰ、Ⅴ級沒有隨機趨勢。

河西走廊地區蘆葦生長環境較為惡劣，該地區蘆葦通過營養繁殖進行種群更新，擴展種群的生態位空間，經過長時間的演化，形成聚集分布的空間格局，來相互庇護抵御不利環境。鹽生蘆葦可能因為種內個體的競爭、生化他感作用，以及生境作用呈現一定程度的隨機分布趨勢，有待于今后深入研究加以驗證。濕生生境水分不再是蘆葦生長的限制因子，種群聚集度比較高；鹽生生境由于水鹽脅迫，土壤水分虧缺和高含鹽量成為鹽生蘆葦種群發展的限制因子，蘆葦為了滿足自身對水分、養分的需求，有隨機分布的傾向，進而獲得較大的空間及資源來維持其生存和發展。