李 潮，謝應忠，許冬梅，李 俠

(寧夏大學農(nóng)學院，寧夏 銀川 750021)

種間關系指不同物種種群之間相互作用所形成的關系，是植物種群間相互聯(lián)系、相互影響的反映，為植物群落結(jié)構(gòu)和數(shù)量的重要特征之一[1]。在群落生態(tài)學研究中，通過對種間關系進行分析，可以確定趨向于生活在一起的物種，進而確定群落邊界。因此，測定植物群落的種間關系，能夠揭示植物種群間在空間上共同出現(xiàn)的程度，從而辨析植物是獨立分布，還是相互吸引或相互排斥[2]。近年來，有關寧夏東部風沙區(qū)荒漠草原植物群落數(shù)量分析的研究主要集中在生態(tài)位測度[3-4]和群落排序[5]等方面，但對植物群落種間關系方面的研究未有深入報道。本研究通過實地調(diào)查，定量研究寧夏東部風沙區(qū)荒漠草原植物群落的種間關系，旨在闡明群落的組成和動態(tài)，完善該地區(qū)群落生態(tài)學的研究，可進一步為區(qū)域生態(tài)環(huán)境建設、生態(tài)安全維護提供理論指導，對該地區(qū)草地生態(tài)系統(tǒng)的可持續(xù)利用具有實用意義。

1 材料與方法

1.1研究區(qū)自然概況 研究區(qū)位于寧夏回族自治區(qū)鹽池縣(106°30′-107°47′ E，37°04′-38°10′ N)，海拔1 295～1 951 m。地勢南高北低，北臨毛烏素沙地，南靠黃土高原，北側(cè)為鄂爾多斯緩坡丘陵，南側(cè)為黃土丘陵。

鹽池縣屬于典型的中溫帶大陸性氣候，光能豐富、熱量偏少，年平均氣溫7.5 ℃，年均無霜期165 d。年平均風速2.8 m·s-1，冬春風沙天氣較多。多年平均降水量280 mm，自西北向東南遞增，其中70%降水集中在6-9月，且降水年際變率大。灰鈣土為主要的土壤類型，其次是風沙土和黑壚土，土壤質(zhì)地以沙土、沙壤和粉沙壤為主[5]。植被以荒漠草原為主，風蝕作用強烈。

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

1.2.1植物群落調(diào)查 2012年7月下旬在寧夏鹽池縣進行野外調(diào)查。采用隨機取樣法進行取樣，結(jié)合樣方法與樣線法對植物群落進行調(diào)查，樣線沿流動沙丘外圍向里布設，在樣線上根據(jù)群落變化布設樣方，樣方面積為1 m×1 m。共布設樣方119個，測定植被蓋度、密度、多度、頻度，同時記錄樣方所在的地理位置、微地形等。

重要值(IV)計算公式為:IV=(相對蓋度+相對密度+相對高度+相對頻度)×100/4[3]。

1.2.2分析方法 依據(jù)調(diào)查結(jié)果，剔除頻度<10%的物種后[6]，以主要物種的重要值作為數(shù)量指標，用于種間關聯(lián)的相關分析。

Jaccard關聯(lián)指數(shù)：按照物種在樣方內(nèi)存在與否，獲得種-樣方的“0、1”二元矩陣，“0”代表不存在，“1”代表存在，用于計算種間關聯(lián)。依據(jù)2×2列聯(lián)表的a、b、c值計算Jaccard關聯(lián)指數(shù)(JI)，JI=a/(a+b+c)。式中，a為兩個種都存在的樣方數(shù)，b為B種存在而A種不存在的樣方數(shù)，c為A種存在而B種不存在的樣方數(shù)。

種間關聯(lián)的2×2列聯(lián)表及χ2檢驗：在2×2列聯(lián)表的基礎上，計算χ2值，可以判斷兩個物種關聯(lián)與否。

數(shù)據(jù)處理用Microsoft Ecxel 2007及SPSS 17.0數(shù)據(jù)處理軟件完成。

2 結(jié)果與分析

2.1寧夏東部風沙區(qū)荒漠草原植物群落特征 在所有調(diào)查樣方中，共鑒定出74種植物，分屬21科57屬，剔除頻度<10%的物種后，剩余21種植物(表1)。其中，多年生草本占52.4%，一年生草本占33.3%，半灌木占14.3%。種數(shù)最多的3個科是禾本科、豆科和藜科。有一部分物種具有較為廣泛的地理分布范圍，如賴草在70個樣地中出現(xiàn)，絲葉山苦荬出現(xiàn)在68個樣地，黑沙蒿出現(xiàn)在62個樣地。

2.2寧夏東部風沙區(qū)荒漠草原植物群落的種間關系 寧夏東部風沙區(qū)荒漠草原植物群落21個主要物種210個物種對數(shù)間的χ2檢驗、Pearson相關系數(shù)和Spearman秩相關系數(shù)的檢驗結(jié)果顯示，χ2檢驗共有21個物種對數(shù)為顯著正關聯(lián)，25個物種對數(shù)為顯著負關聯(lián)；Pearson相關系數(shù)有12個物種對數(shù)為顯著正相關，6個物種對數(shù)為顯著負相關；Spearman秩相關系數(shù)有21個物種對數(shù)為顯著正相關，10個物種對數(shù)為顯著負相關(表2)。

2.2.1χ2檢驗與Jaccard關聯(lián)指數(shù) 校正后的χ2檢驗值排出的半矩陣圖表達的是種間聯(lián)結(jié)的程度，這有助于對群落中植物種間關系的認識。圖1所示，正關聯(lián)的物種對數(shù)為103個，負關聯(lián)的物種對數(shù)為106個，無關聯(lián)的物種對數(shù)為1個。正、負關聯(lián)物種對數(shù)分別占總物種對數(shù)的49.05%和50.48%。物種間呈正關聯(lián)可能是因一個物種依賴于另一個物種而存在，或兩者受生物因子和非生物因子影響而共存。負關聯(lián)則可能是因競爭、空間排擠或他感作用，或不同的環(huán)境要求[7]。負關聯(lián)物種對數(shù)在植物群落中的出現(xiàn)率較大，說明該區(qū)域的植物群落不穩(wěn)定，當受到外界因素干擾時，容易發(fā)生波動，群落仍處于發(fā)展演替之中。

χ2檢驗物種對數(shù)間呈極顯著、顯著關聯(lián)的物種對數(shù)有46個，占總物種對數(shù)的21.9%。其中物種對數(shù)為顯著正關聯(lián)共有21個，占總物種對數(shù)的10.0%，物種對數(shù)為顯著負關聯(lián)共有25個，占總物種對數(shù)的11.9%，呈顯著正、負關聯(lián)的物種對數(shù)相差較小。群落系統(tǒng)發(fā)育或群落演替過程愈成熟，其植物種類組成也愈穩(wěn)定，種間關系也愈趨向正關聯(lián)，以保證物種間的共存[8]。物種對數(shù)間呈極顯著、顯著關聯(lián)的物種對數(shù)不多，說明該地區(qū)植物種群之間的關聯(lián)性較弱。

Jaccard關聯(lián)指數(shù)的半矩陣圖表達的是物種共同出現(xiàn)的百分率。由圖2看出，JI為0的物種對數(shù)為22個，占總物種對數(shù)的10.5%；JI值在0～0.2的物種對數(shù)有150個，占總物種對數(shù)的71.4%，說明植物群落的主要物種間關聯(lián)程度很小，存在相當程度的獨立分布格局。同時，JI值大于0.4的物種對數(shù)為2個，僅占總物種對數(shù)的0.95%，這也進一步說明該地區(qū)植物群落的種間關聯(lián)性相對較弱。

2.2.2Pearson 相關分析與Spearman秩相關分析χ2檢驗只涉及種的存在與否，不能反映物種對數(shù)間的數(shù)量變化關系，因此，采用Pearson 相關系數(shù)和Spearman 秩相關系數(shù)進一步定量分析植物群落的種間關系。通過Pearson相關分析，呈顯著相關的物種對數(shù)為18個，其中呈極顯著正相關的6個，呈顯著正相關的6個，呈極顯著負相關的兩個，呈顯著負相關4個(圖3、表2)。Spearman秩相關分析表明，呈顯著相關的物種對數(shù)為31個，其中呈極顯著正相關的18個，呈顯著正相關的3個，呈極顯著負相關的5個，呈顯著負相關5個(圖4、表2)。Spearman秩相關分析和χ2檢驗得出的呈顯著相關的物種對數(shù)均多于Pearson相關分析的結(jié)果(表2)，原因在于，基于Pearson相關系數(shù)的統(tǒng)計分析對于數(shù)據(jù)的分布類型十分敏感，只有在數(shù)據(jù)呈近似正態(tài)分布的時候才能被采用，否則就可能導致錯誤的結(jié)論。而自然界大多數(shù)物種多呈簇生分布，所以Pearson相關分析的結(jié)果帶有一定的局限性[9]。

比較圖1和圖4，在Spearman秩相關分析和χ2檢驗中，呈正相關的物種對數(shù)的數(shù)量相同，但χ2檢驗中呈顯著正相關的物種對數(shù)在Spearman秩相關分析中均為極顯著正相關。而對于呈顯著負相關的物種對數(shù)的數(shù)量，χ2檢驗明顯多于Spearman秩相關分析，這是由于χ2檢驗是定性分析，只涉及物種在樣地內(nèi)存在與否，不涉及物種的數(shù)量特征，不能非常準確地描述物種對數(shù)間的關聯(lián)性強弱。因此，本研究只對Spearman秩相關分析結(jié)果進行說明。

表2 寧夏東部風沙區(qū)荒漠草原21種主要植物種間關系的3種檢驗結(jié)果Table 2 Results of interspecific relationship produced by χ2-test, Pearson’s correlation coefficient, Spearman’s rank correlation coefficient, respectively

圖1 寧夏東部風沙區(qū)荒漠草原植物群落種間關聯(lián)性的χ2檢驗半矩陣圖Fig.1 The semi-matrix of χ2-test among 21 major species

Spearman秩相關分析呈極顯著正相關的物種對數(shù)為18個，如紫翅豬毛菜、地錦、狗尾草、尖頭葉藜兩兩物種之間都呈極顯著正相關。這是因為這4種植物都為一年生草本，它們之間相互依賴，資源利用方式相似，對生境的適應和反應相同。而物種對數(shù)間呈極顯著負相關，表明兩物種由于空間排擠、競爭，以及不同的資源利用方式，對生境要求不同，出現(xiàn)生態(tài)位的分離[10]。

3 討論

3.1種間關聯(lián)測定 種間關聯(lián)測定有多種方法，主要有χ2檢驗、共同出現(xiàn)百分率法、方差比率法、關聯(lián)系數(shù)、分布相似性百分率法、相關系數(shù)、點相關系數(shù)等[11]。本研究應用的是χ2檢驗和Jaccard關聯(lián)指數(shù)兩種方法。其中，在χ2檢驗中呈極顯著和顯著正關聯(lián)的物種對數(shù)，其JI值均較大。如絲葉山苦荬-狗尾草的JI值是0.41，χ2檢驗中也呈極顯著關聯(lián)，這說明Jaccard聯(lián)結(jié)指數(shù)和χ2檢驗在反映物種間的關聯(lián)程度方面較為一致。但是，JI值的運用也有局限性，如賴草-絲葉山苦荬的χ2檢驗值很小，物種對數(shù)間并沒有表現(xiàn)出顯著關聯(lián)性，而其JI值是0.41。在測定正、負關聯(lián)的過程中，賴草-絲葉山苦荬是所有物種對數(shù)里唯一一個無關聯(lián)的物種對，其ad=bc，說明賴草、絲葉山苦荬兩個物種隨機出現(xiàn)，無關聯(lián)。同時從χ2檢驗的數(shù)學角度來看，這兩個種出現(xiàn)的概率恰與概率統(tǒng)計的期望值相等，說明賴草-絲葉山苦荬間完全獨立。以上表明χ2檢驗方法和Jaccard關聯(lián)指數(shù)應該同時使用才能更準確、更合理地描述和測定種間關聯(lián)程度。

圖2 寧夏東部風沙區(qū)荒漠草原植物群落種間關聯(lián)性的Jaccard關聯(lián)指數(shù)半矩陣圖Fig.2 The semi-matrix of Jaccard association index among 21 major plant species

圖3 寧夏東部風沙區(qū)荒漠草原植物群落種間關聯(lián)性的Pearson相關系數(shù)半矩陣圖Fig.3 The semi-matrix of Pearson’s correlation coefficient among 21 major species

圖4 寧夏東部風沙區(qū)荒漠草原植物群落種間關聯(lián)性的Spearman秩相關系數(shù)半矩陣圖Fig.4 The semi-matrix of Spearman’s rank correlation coefficient among 21 species

單獨運用JI值的另外一個局限是在χ2檢驗中呈顯著負關聯(lián)的物種對，JI值沒有表現(xiàn)出相應數(shù)值較大的特征，這與χ2檢驗中呈顯著正關聯(lián)物種對的JI值較大的情況不同。這種情況可能是由于Jaccard關聯(lián)指數(shù)公式的分子為a，a表示含有兩個種A和種B的樣方數(shù)，當物種對數(shù)間呈正關聯(lián)時a較大，而當物種對數(shù)間為負關聯(lián)時，a較小，b、c較大，而分母(a+b+c)的值不管是正關聯(lián)還是負關聯(lián)，始終變化不大。因此，當物種對數(shù)間為正關聯(lián)時，a較大，所以JI值也較大，可以較好地反映正關聯(lián)的強弱。而當物種對數(shù)間為負關聯(lián)時，a較小，其JI值也較小，所以，Jaccard聯(lián)接指數(shù)不能較好地反映物種對數(shù)間負關聯(lián)性的顯著情況。同時，χ2檢驗也有其局限性，因為χ2檢驗使用的是存在與否的二元數(shù)據(jù)，不能區(qū)分關聯(lián)程度的大小?？傮w來看，Jaccard聯(lián)接指數(shù)和χ2檢驗都有助于客觀測定兩物種種間關聯(lián)性，一定程度上可以說明物種對數(shù)間的相互關系，并表明其關聯(lián)程度。因此，在實際運用中，要結(jié)合兩種方法使用才能更好地反映種間關系[12]。

3.2種間相關測定 對于χ2檢驗定性分析的局限性，Spearman秩相關分析方法可以有效規(guī)避這一部分。如短花針茅-紫翅豬毛菜、短花針茅-尖頭葉藜、短花針茅-霧冰藜在χ2檢驗中均為顯著負相關關系，而在Spearman秩相關分析結(jié)果里并無顯著相關關系。這可能是由于短花針茅在χ2檢驗中使用的二元數(shù)據(jù)，與Spearman秩相關分析中使用的重要值數(shù)據(jù)有較大差異，因為不同樣地中的重要值數(shù)據(jù)不相等且數(shù)值差距也較大，所以分析過程中所使用數(shù)據(jù)的差異導致了兩種方法出現(xiàn)不同的分析結(jié)果。重要值的數(shù)量特征更細化、更確切地描述了群落的特征，所以使用Spearman秩相關分析方法可以更合理地對種間相關性進行測定。種間相關可以彌補χ2檢驗利用二元數(shù)據(jù)進行分析的不足，尤其是避免了伴生種、偶見物種對種間關系研究的影響[9]。

3.3種間關系與群落演替 對于草地植物群落的種間關系與群落演替之間的關系，很多學者已對此做了研究。劉玨宏等[6]發(fā)現(xiàn)內(nèi)蒙古錫林郭勒草原研究區(qū)處于演替中期，未形成總體關聯(lián)性，呈正相關狀態(tài)。種間關聯(lián)較為松散，總體處于相對共存的狀態(tài)。房飛等[13]對高寒草原植物群落進行研究，發(fā)現(xiàn)原始生態(tài)種之間多呈正關聯(lián)，處于退化階段的物種則表現(xiàn)出無關聯(lián)或負關聯(lián)，一定程度上反映了群落演替的逆向過程。在本研究中，物種對數(shù)間呈顯著、極顯著相關的物種對數(shù)遠低于呈不顯著相關的對數(shù)，說明該區(qū)域植物群落的種間關系松散。物種對數(shù)間呈顯著正相關或顯著負相關，說明群落在長期的演替過程中，各種群之間既相互依存又相互排斥。一部分種群占據(jù)適宜的生態(tài)位，和諧相處，而另一部分種群獨立性較強或進行競爭[8]，所以該區(qū)域植物群落尚處于不斷的發(fā)展演替之中。通過對群落種間關系進行研究，可以找出群落內(nèi)和群落間的聯(lián)系，有助于認識植物群落的分布規(guī)律和動態(tài)演替過程，并為區(qū)域生態(tài)環(huán)境建設和合理利用草地資源提供依據(jù)。

4 結(jié)論

寧夏東部風沙區(qū)荒漠草原植物群落的種間關聯(lián)性較弱，旱生叢狀禾草、雜類草和一年生草較多，并有大量沙生植物，群落穩(wěn)定性低。物種對數(shù)間呈負關聯(lián)物種對數(shù)多于呈正關聯(lián)物種對數(shù)，呈極顯著和顯著相關的物種對數(shù)所占比例較小，僅為14.76%(Spearman秩相關分析結(jié)果)，物種對間競爭水平低，群落結(jié)構(gòu)簡單化，群落種間關系松散，物種趨向于獨立分布，仍處于群落演替之中。

[1] 王伯蓀.植物種群學[M].廣州：中山大學出版社，1989：115-120.

[2] 吳征鎰.中國植被[M].北京：科學出版社，1980.

[3] 張晶晶，許冬梅.寧夏荒漠草原不同封育年限優(yōu)勢種群的生態(tài)位特征[J].草地學報，2013，21(1):73-78.

[4] 程中秋，張克斌，劉建，等.寧夏鹽池荒漠草原區(qū)天然草地植物生態(tài)位研究[J].水土保持研究，2011，18(3):37-40.

[5] 方楷，宋乃平，安慧，等.寧夏鹽池荒漠草原植被的數(shù)量分類和排序[J].生態(tài)學雜志，2011，30(12):2719-2725.

[6] 劉玨宏，高慧，張麗紅，等.內(nèi)蒙古錫林郭勒草原大針茅——克氏針茅群落的種間關聯(lián)特征分析[J].植物生態(tài)學報，2010，34(9):1016-1024.

[7] 李博.生態(tài)學[M].北京：高等教育出版社，2000：122-124.

[8] 杜道林，劉玉成，李睿.縉云山亞熱帶烤樹林優(yōu)勢種群間聯(lián)結(jié)性研究[J].植物生態(tài)學報，1995，19(2):149-157.

[9] 張峰，上官鐵梁.山西翅果油樹群落種間關系的數(shù)量分析[J].植物生態(tài)學報，2000，24(3):351-355.

[10] 張金屯，焦蓉.關帝山神尾溝森林群落木本植物種間聯(lián)結(jié)性與相關性研究[J].植物研究，2003，23(4):458-463.

[11] 彭少麟，周厚誠，郭少聰，等.鼎湖山地帶性植被種間關系研究[J].植物學報，1999，41:1239-1244.

[12] 李軍玲，張金屯，袁建英.關帝山亞高山灌叢群落和草甸群落優(yōu)勢種的種間關系[J].草地學報，2004，12(2):113-119.

[13] 房飛，胡玉昆，張偉，等.高寒草原植物群落種間關系的數(shù)量分析[J].生態(tài)學報，2012，32(6):1898-1907.