魯為華，朱進忠，靳瑰麗

（1.石河子大學動物科技學院，新疆 石河子832003；2.新疆農(nóng)業(yè)大學草業(yè)與環(huán)境科學學院，新疆 烏魯木齊830052；3.新疆草地資源與生態(tài)重點實驗室，新疆 烏魯木齊830052）

幼苗的出苗格局往往和種子的形成以及傳播過程密切相關(guān)［1］，種子落粒過程要受到母株狀態(tài)的影響［2－4］。有些植物種子只能在短距離內(nèi)傳播甚至傳播距離為零。有學者認為短距離的種子傳播對于植物種群拓展而言是不成功的［5］。但實際上短距離擴散有利于在環(huán)境均一的條件下種群內(nèi)斑塊的形成，同時也有利于競爭種在異質(zhì)性環(huán)境中共存［6－9］。

種子傳播嚴格來講是屬于種群水平上發(fā)生的行為，單株傳播后形成的格局有可能由于母株個體數(shù)量太少而使得隨機效應被無形中放大。因此如何用同一種分析方法將多個個體納入到分析樣本中而實現(xiàn)對種子擴散過程的確切描述是一個現(xiàn)實問題。Thorsten Wiegand研究組利用“多樣方圖像整合技術(shù)”（combining the data from individual mapped replicate plots into mean，以下簡稱 CIM）成功實現(xiàn)了多個小尺度樣方的點格局分析，之后并以這種方法進行了大量的研究工作［10］。

近年來，眾多學者對絹蒿（Seriphidium）荒漠進行了大量的研究工作，包括絹蒿荒漠的植被動態(tài)變化［11］、絹蒿種子結(jié)構(gòu)特征及萌發(fā)特性［12，13］、退化過程中土壤和植被的變化規(guī)律［14，15］以及絹蒿生物學和生理學特性［16，17］、幼苗的群體變化規(guī)律［18］、成株的分布特征［19］等。但對絹蒿種子傳播后形成的幼苗空間分布格局未進行充分研究。因此，對絹蒿種群的種子傳播后發(fā)生的幼苗格局進行小尺度分析，揭示種子短距離擴散后形成的幼苗在不同尺度上的發(fā)生格局，期望從小尺度的空間分布格局中認識到絹蒿種群的某些生態(tài)學過程，為該類草地的合理利用和生態(tài)恢復提供理論依據(jù)。

1 材料與方法

1.1 試驗方法

首先在研究區(qū)由北向南設(shè)置3條樣線，統(tǒng)計包括絹蒿在內(nèi)的植物種類、密度、生物量、蓋度等指標，利用聚類分析方法將整個地段劃分為重度退化（SD，serious degradation）和極度退化（HD，heavy degradation）2個梯度［19］，重度退化區(qū)主要以叉毛蓬（Petrosimonia sibirica）為主，植被密度較大。極度退化區(qū)域主要以叉毛蓬為主，伴生有駱駝蓬（Peganum）等植物，植被稀疏，并有大量的裸地。兩區(qū)域內(nèi)絹蒿個體分布稀疏，可完全消除鄰接干擾效應。在2008年整個植物群落的生長旺季，對重度和極度退化2個樣地中大小級別相似的絹蒿成年植株和其周圍發(fā)生的幼苗坐標點進行詳細記錄，每個樣地取10個單株，2樣地共計20個單株。記錄時以單個母株為中心點，布置1m×1m的小樣方，盡量將所有幼苗控制在樣方范圍內(nèi)。以樣方正南正北方向，1m×1m樣方的左下角為坐標位置的參考點，按地上絹蒿幼苗的個體記錄，只要有絹蒿幼苗個體出現(xiàn)，就把該點坐標位置賦予絹蒿，并對每一絹蒿幼苗坐標位置進行標簽標注。

在對絹蒿幼苗定位的基礎(chǔ)上，按月統(tǒng)計樣方內(nèi)幼苗的發(fā)生數(shù)量以及死亡狀況，統(tǒng)計其數(shù)量的動態(tài)變化。

1.2 數(shù)據(jù)處理

采用生存分析的Kaplan-Meier來統(tǒng)計幼苗存活率大小，并利用Log-Rank秩相關(guān)以及Breslow檢驗來比較不同退化梯度下的幼苗存活能力差異［20］；利用Programita 2006軟件的O-ring函數(shù)進行點格局分析，同時結(jié)合軟件包內(nèi)的CIM技術(shù)將若干單株幼苗的點格局分析結(jié)果通過加權(quán)平均，將不同單株幼苗的分布狀況反映在相同尺度上［21］。該方法中，雙變量O-ring 統(tǒng)計的O12（r）可以用公式（1）計算，該公式通用于單變量和雙變量O-ring統(tǒng)計過程。

式中，n1為格局1（雙變量統(tǒng)計中的對象1）的點的數(shù)目；Rw1，i（r）表示以格局1中第i個點為圓心、半徑為r、寬為w的圓環(huán)；Point2［X］計算了區(qū)域X中格局2（雙變量統(tǒng)計中的對象2）的點的數(shù)目；Area［X］是區(qū)域X的面積。

為了將多個重復且進行相同處理的樣本等權(quán)納入到計算公式，根據(jù)Thorsten Wiegand的意見，在計算多個重復小樣方時，可以將公式（1）拓展為公式（2）。

公式（2）中，ij為第j個重復樣方的第i個點，n1j為格局1和重復j的點數(shù)目，N＝∑jn1j為所有重復樣方中格局1的點總數(shù)目，而公式（2）又可進一步簡化為公式（3）。

為檢驗O12（r）偏離零值的顯著程度，可以用Monte Carlo隨機模擬方法擬合它們95%的置信區(qū)間。用r作為橫坐標，上下包跡線作為縱坐標繪圖，置信區(qū)間一目了然。用種群實際分布數(shù)據(jù)（點圖）計算得到不同尺度下的O12（r）值若在包跡線以內(nèi)，則符合隨機分布；若在包跡線以上，則為聚集分布，如在包跡線以下則為均勻分布。

2 結(jié)果與分析

2.1 兩種退化梯度下單株產(chǎn)生幼苗數(shù)量動態(tài)變化及生存率分析

從2008年4月份開始，按月對所選定的植株以母株為中心定點觀測幼苗數(shù)量變化狀況，統(tǒng)計幼苗死亡數(shù)量和新出現(xiàn)的幼苗數(shù)量，并采用Log-Rank和Breslow檢驗來分析不同退化梯度下幼苗的生存力。

2種退化梯度下的幼苗數(shù)量存在差異（圖1），但其數(shù)量隨時間變化動態(tài)趨勢相似，在4月份有1個萌發(fā)高峰，完成整個萌發(fā)量的90%以上，5月份在數(shù)量上變化不大，但略微增長，之后由于大規(guī)模死亡而數(shù)量迅速下降。因此，不論在哪種條件下都呈現(xiàn)出“爆發(fā)性出苗”和“驟然性死亡”的特點。

為進一步比較2種退化梯度下幼苗的生存率大小，采用醫(yī)學統(tǒng)計中常用的生存率分析來比較樣本之間的生存率差異，2種退化梯度下的Kaplan-Meier生存率曲線自7月份開始出現(xiàn)差別（圖2），Log-Rank檢驗結(jié)果中，卡方值＝0.275 5，P＝0.599 6＞0.05，因此2種退化梯度下生存曲線不存在差別。同時，因為重度退化樣本觀察值n＝47，小于期望值E＝45.289 98，因此可以認為重度退化條件下幼苗死亡率稍高。進一步以觀察個體數(shù)量為權(quán)重進行Breslow檢驗結(jié)果，卡方值＝0.0 0 3 2，P＝0.954 7，說明在考慮觀察對象數(shù)量后其生存曲線也不存在顯著差異。說明在小尺度下由于單株對幼苗發(fā)生和存活提供了同質(zhì)的庇護條件，造成二者生存率不存在顯著差異。

圖1 兩種退化梯度下單株幼苗數(shù)量動態(tài)變化及幼苗死亡數(shù)量Fig.1 Quantity of survival seedling，mortality of seedling under two gradation degradation

圖2 兩種退化梯度下幼苗生存率曲線及Breslow檢驗結(jié)果Fig.2 Kaplan-Meier survival rate curve，Log-Rank＆ Breslow test results in SD and HD

2.2 幼苗空間分布格局分析

將幼苗的點坐標數(shù)據(jù)轉(zhuǎn)換為以母株為圓心，以最大傳播距離為半徑的點數(shù)據(jù)，以探討幼苗在母株周圍的分布狀態(tài)，并繪制雷達點分布圖，由于個體數(shù)量限制，這里在2個退化梯度中分別只選取3個較典型植株進行繪圖比較，具體結(jié)果見圖3。

2個退化梯度內(nèi)的3個單株5月所產(chǎn)生的幼苗數(shù)量是比較可觀的（圖3），但隨著時間進程（5－9月）其幼苗數(shù)量急劇下降，且存活的幼苗一般都距離母株很接近，重度退化梯度下，5月份幼苗與母株最遠距離可以達到近60cm左右，7月份則降至45cm左右，到9月份則縮短至25cm左右。極度退化條件下幼苗的分布也有類似的趨勢。

進一步利用Programita 2006軟件的O-ring函數(shù)和CIM技術(shù)對幼苗在不同尺度上的分布狀態(tài)進行分析（圖4）。重度退化條件下，5月份幼苗密度較大，因此在小尺度上聚集強度較大，其密度函數(shù)最大值出現(xiàn)在半徑為6 cm的圓環(huán)上，其值為0.022 2，之后呈直線下降，但在半徑小于35cm之前的一系列圓環(huán)上，仍然呈現(xiàn)聚集分布，在到達40cm尺度之后轉(zhuǎn)入隨機或者均勻分布。到7月份，幼苗數(shù)量由于環(huán)境篩選而數(shù)量減少，因此聚集強度與5月相比顯著下降，最大聚集強度出現(xiàn)在半徑為6cm的圓環(huán)上，其函數(shù)最大值為0.016 67，并且偏離聚集分布的尺度提前至33cm；而到9月份，幼苗數(shù)量進一步減少，其格局分布雖然是在小尺度上呈聚集分布，但是曲線出現(xiàn)明顯的波谷，并且聚集強度總體上進一步減小，其密度函數(shù)最大值出現(xiàn)在半徑為7cm的圓環(huán)上，其值為0.013 02，之后在22cm尺度上轉(zhuǎn)入隨機分布。

極度退化條件下，由于幼苗數(shù)量隨著退化程度增加其數(shù)量較少，因此，5月份聚集強度與重度退化比較顯著降低，最大聚集強度出現(xiàn)在半徑為4cm的圓環(huán)上，其值為0.015 3，在到達28cm時，轉(zhuǎn)入隨機分布。到7月份，小尺度內(nèi)聚集強度增加，其O（r）最大值仍然在4cm的圓環(huán)上，其值為0.019 1，在25cm時則轉(zhuǎn)入隨機分布，37 cm時則轉(zhuǎn)入均勻分布。9月份，曲線整體趨勢未發(fā)生明顯變化，但在4cm圓環(huán)上其聚集強度與7月相比略微減小，其O（r）值為0.019 0，在19cm尺度上則轉(zhuǎn)入隨機分布。從O（r）值變化曲線可以看出，絹蒿幼苗種子傳播距離分布在5～30cm內(nèi)，其種子集中落在5～10cm的尺度內(nèi)。當然，不同形體大小的植株種子落粒距離可能會有所不同，但本試驗中所選取的植株種子基本在此范圍內(nèi)。隨時間推移，母株周圍的種子有可能因為其他因素被傳播到更遠的距離，其潛在生態(tài)位空間有可能會加大。

圖3 重度（A）和極度（B）退化3個單株產(chǎn)生幼苗5、7、9月的整合點分布圖（坐標單位：cm）Fig.3 CIM map of three maternal Seriphidiumseedling in three mounth（Unit：cm）

圖4 重度和極度退化條件下3個時間點3個單株幼苗經(jīng)CIM整合后的點格局分布圖Fig.4 Combined spatial patterns of emergence seedling on three maternal piant by CIM in temporal sequence on SD and HD

在2種退化條件下，各時間點上幼苗在空間上的分布格局變化趨勢相似，但是經(jīng)過夏季干旱及高溫等環(huán)境因子的強烈篩選之后，距離母株較遠的幼苗大部分已經(jīng)死亡，因此在9月份的分布格局與初始格局已經(jīng)有了明顯不同，其有效傳播距離尺度縮小，在重度退化區(qū)域已經(jīng)縮小到25cm左右，極度退化區(qū)域縮小至20cm左右。從時空尺度來看，絹蒿幼苗對母株具有很強的依賴性，其格局變化總是朝有利于自身生存的方向發(fā)展。

3 結(jié)論與討論

絹蒿幼苗在較早時間萌發(fā)是一種對中亞荒漠干旱氣候的適應性對策，一般在4－5月份土壤中的種子就基本全部萌發(fā)，盡早萌發(fā)無形中增加了幼苗適合度，因為這種對策使得幼苗個體盡早獲得有限的資源并因此而增加其生存的幾率并對之后的生長、繁殖等生活史的一系列過程產(chǎn)生積極影響［22，23］。從其生存率來看，大量的幼苗集中萌發(fā)為環(huán)境篩選提供了足夠數(shù)量的個體，使得最終還有一部分幼苗存活下來，這將對種群的更新和繁衍提供足夠的物質(zhì)基礎(chǔ)。

絹蒿種子落粒對母株具有強烈的依賴性，通過多尺度變換的點格局分析發(fā)現(xiàn)，在不同尺度上，幼苗的分布格局是出現(xiàn)變化的，在小尺度上，總是呈現(xiàn)出強烈的聚集分布，而隨著尺度增加，其聚集強度逐漸減小，隨著時間變化，聚集強度也在發(fā)生變化，也就是說，種子向幼苗以及成株轉(zhuǎn)變過程中在時間尺度上面臨著強大的環(huán)境選擇壓力，其分布格局總是傾向于使幼苗建成成功率增加的方向轉(zhuǎn)換。

［1］Walker L R，del Moral R，Primary Succession and Ecosystem Rehabilitation［M］.Cambridge，UK：Cambridge University Press，2003：442.

［2］Harper J L.The role of predation in vegetational diversity［J］.Brookhaven Symposia in Biology，1969，22：48-62.

［3］Crawley.Plant-herbivore dynamics［A］.In：Crawley M.Plant Ecology［M］.Oxford：Blackwell Science，1997：401-474.

［4］Kitajima K，F(xiàn)enner M.The Ecology of Regeneration in Plant Communities（2nd edition）［C］.Wallingford：CAB International Publishing，2000：331-360.

［5］Cousens R D，Wiegand T.Small-scale spatial structure within patterns of seed dispersal［J］.Oecologia，2008，158（3）：437-448.

［6］Snyder R E，Chesson P.Local dispersal can facilitate coexistence in the presence of permanent spatial heterogeneity［J］.Ecology Letters，2003，6：301-309.

［7］Murrell D J，Law R.Heteromyopia and the spatial coexistence of similar competitors［J］.Ecology Letters，2003，6：48-59.

［8］Green D G.Connectivity and complexity in ecological systems［J］.Pacific Conservation Biology，1994，1：194-200.

［9］劉鐘齡，王煒，梁存柱，等.內(nèi)蒙古草原植被在持續(xù)牧壓下退化演替的模式與診斷［J］.草地學報，1998，6（4）：244-251.

［10］De Luis Martín，Ravents J，Wiegand T.Temporal and spatial differentiation in seedling emergence may promote species coexistence in Mediterranean fire-prone ecosystems［J］.Ecography，2008，31：620-629.

［11］劉長娥，許鵬，安沙舟.土質(zhì)溫性荒漠草地植物生育節(jié)律分析［J］.新疆農(nóng)業(yè)科學，2006，43（1）：11-13.

［12］張連義.伊犁絹蒿種子天然種衣吸水保水性能研究［J］.植物學報，1998，8：775-777.

［13］孫會忠，宋月芹，賀學禮.絹蒿屬8種植物果實微形態(tài)特征及其分類學意義［J］.武漢植物學研究，2007，25（5）：440-444.

［14］顧愛星，范燕敏，武紅旗，等.天山北坡退化草地土壤環(huán)境與微生物數(shù)量的關(guān)系［J］.草業(yè)學報，2010，19（2）：116-123.

［15］劉洪來，朱進忠，靳瑰麗.伊犁絹蒿荒漠不同退化階段草地植物多樣性分析［J］.新疆農(nóng)業(yè)科學，2007，44（5）：632-636.

［16］孫宗玖，安沙舟，許鵬.伊犁絹蒿構(gòu)件動態(tài)變化研究［J］.草地學報，2007，15（5）：454-459.

［17］孫宗玖，安沙舟，許鵬.伊犁絹蒿體內(nèi)可塑性營養(yǎng)物質(zhì)動態(tài)變化及分配特征的研究［J］.草業(yè)學報，2008，17（2）：151-156.

［18］魯為華，朱進忠，王東江，等.天山北坡圍欄封育條件下伊犁絹蒿幼苗分布格局及數(shù)量動態(tài)變化規(guī)律研究［J］.草業(yè)學報，2009，18（4）：17-26.

［19］魯為華，朱進忠，王東江，等.天山北坡兩種退化梯度下伊犁絹蒿種群空間分布的點格局分析［J］.草業(yè)學報，2009，18（5）：142-149.

［20］Kaplan E L，Meier P.Nonparametric estimation from incomplete observations［J］.Journal of the American Statistical Association，1958，53：457-481.

［21］Wiegand T，Moloney K.Rings，circles and null models for point pattern analysis in ecology［J］.Oikos，2004，104：209-229.

［22］Pausas J G.Plant functional traits in relation to fire in crown-fire ecosystems［J］.Ecology，2004，85：1085-1100.

［23］Bolker B M，Pacala S W.Spatial moment equations for plant competition：understanding spatial strategies and the advantages of short dispersal［J］.American Naturalist，1999，153：575-602.