李建慶,梅增霞,楊忠岐

1 濱州學院, 山東省黃河三角洲生態(tài)環(huán)境重點實驗室, 濱州　256603　　2 中國林業(yè)科學研究院森林生態(tài)環(huán)境與保護研究所,國家林業(yè)局森林保護學重點實驗室, 北京　100091

?

不同林分白蠟樹云斑白條天牛種群空間格局地統(tǒng)計學分析

李建慶1,梅增霞1,楊忠岐2,*

1 濱州學院, 山東省黃河三角洲生態(tài)環(huán)境重點實驗室, 濱州2566032 中國林業(yè)科學研究院森林生態(tài)環(huán)境與保護研究所,國家林業(yè)局森林保護學重點實驗室, 北京100091

云斑白條天牛是我國重要林業(yè)害蟲,在黃河三角洲地區(qū)對白蠟樹造成嚴重危害。利用地統(tǒng)計學方法,通過調(diào)查城市行道綠化林、庭院綠化林、公路林和片林4種不同林分類型的白蠟樹林地,以產(chǎn)卵刻槽、排糞孔和羽化孔為調(diào)查指標,對危害白蠟云斑白條天牛種群卵、幼蟲、成蟲(蛹)的空間格局進行了研究。結(jié)果表明：除片林樣點成蟲(蛹)的半變異函數(shù)模型為指數(shù)模型外,危害白蠟云斑白條天牛種群卵、幼蟲、成蟲(蛹)在不同林分類型下的半變異函數(shù)模型均為球形模型,其指明空間分布均為聚集分布,種群的聚集分布習性與林分類型相關(guān)性不強。但不同林分類型云斑白條天牛種群的空間相關(guān)性地統(tǒng)計學參數(shù)存在差異,片林調(diào)查樣點的卵、幼蟲和成蟲(蛹)的變程a、塊金值C0和基臺值C0+C均小于其他3種林分類型的統(tǒng)計值,片林樣點成蟲(蛹)的半變異函數(shù)模型不同于其他林分類型為指數(shù)模型。

云斑白條天牛；空間分布；地統(tǒng)計學；白蠟樹

白蠟樹是一優(yōu)良園林綠化樹種,由于其耐鹽堿能力強,在土壤鹽漬化較為嚴重的黃河三角洲地區(qū)廣泛種植。云斑白條天牛Batoceralineolata,又名云斑天牛,屬鞘翅目天牛科,是我國重要林業(yè)害蟲,在黃河三角洲地區(qū)對白蠟樹的危害尤其嚴重,在當?shù)貫I州和東營城區(qū)管理不善的路段,白蠟行道林受害株率可達88.24%[1],筆者調(diào)查路段甚至達100%,受害嚴重的白蠟樹蟲口數(shù)可達32個[2],云斑白條天牛的危害已經(jīng)限制了當?shù)匕紫灅涞臄U大種植。因此,研究危害白蠟樹云斑白條天牛種群空間格局、動態(tài)變化和增長趨勢,對指導云斑白條天牛的防治具有重要意義。

為做好白蠟樹云斑白條天牛(Batoceralineolata)的防治工作,需要進一步深入研究其空間布局等種群分布特點。利用地統(tǒng)計學的方法來分析云斑白條天牛的種群空間格局,可彌補傳統(tǒng)統(tǒng)計方法的不足,最大限度分析利用其調(diào)查數(shù)據(jù),獲取精確的種群信息,可有效的避免系統(tǒng)誤差,更好的指導生產(chǎn)防治[3-4]。李友常、柳林俊等利用地統(tǒng)計學分析了楊樹光肩星天牛Anoplophoraglabripennis和青楊天牛Saperdapapulnea的空間格局[5-6],但由于天牛類害蟲蛀干危害的特殊性,利用地統(tǒng)計學方法來研究天?？臻g格局的報道還不多,云斑白條天牛種群空間格局地統(tǒng)計學分析的研究報道尚未見到。本研究利用地統(tǒng)計學方法研究了不同林分白蠟樹云斑白條種群不同蟲態(tài)空間格局,以期為種群預(yù)測預(yù)報和制定合理防治策略提供理論指導。

1　材料與方法

1.1試驗地概況

試驗地設(shè)在黃河三角洲地區(qū)的山東省東營和濱州市,地理位置為36°41′—38°16′N,117°15′—119 °10′E,屬北溫帶濕潤氣候區(qū),一年四季分明,年平均氣溫在11.7—12.6℃之間。該區(qū)域地勢平坦,海拔較低且蒸發(fā)量較大,土壤次生鹽堿化嚴重,耐鹽堿的白蠟樹是當?shù)刂匾G化樹種,大量用于城區(qū)和道路園林綠化。

試驗地栽植樹種為絨毛白蠟,主要林分類型為城市行道綠化林、庭院綠化林、公路林和片林,依據(jù)林分類型選取試驗地4塊(表1)。選取的試驗地代表了不同胸徑大小的白蠟林地,片林由于栽植密度高,種間競爭大,生長慢,平均胸徑8.69 cm,城市行道綠化林、庭院綠化林和公路林為道路兩側(cè)單行種植,種間競爭小,生長快,平均胸徑在14.38—18.27 cm之間。選取的試驗地云斑白條天牛危害率較高,有蟲株率在22.93%—100%之間,但白蠟樹處于天牛入侵危害的前期和中期,外在長勢基本正常。

表1　白蠟調(diào)查樣地基本情況統(tǒng)計表

表中數(shù)據(jù)為平均數(shù)±標準誤

1.2調(diào)查方法

分別選取城市行道綠化林、庭院綠化林、公路林和片林各1塊,采取逐株取樣調(diào)查的方法,統(tǒng)計每株白蠟上的云斑白條天牛產(chǎn)卵刻槽、排糞孔和羽化孔的數(shù)量,分別代表云斑白條天牛的卵、幼蟲和成蟲(或蛹)的蟲口數(shù)。排糞孔調(diào)查時間為5月中旬,產(chǎn)卵刻槽和羽化孔的調(diào)查時間為7月上旬。

1.3空間格局分析方法[7]

半方差函數(shù)分析種群空間格局,實質(zhì)是分析表征種群數(shù)量的變量空間變異特征和程度。對觀察數(shù)據(jù)序列Z(xi),i=1,2,…,n,樣本半變異函數(shù)值γ(h)的公式如下：

式中,N(h)是被h分割的數(shù)據(jù)對(xi,xi+h)的對數(shù),Z(xi)和Z(xi+h)分別是點xi和xi+h處樣本的測量值,h是分隔兩樣點的距離。

半方差函數(shù)曲線圖的形狀反映了變量空間分布的結(jié)構(gòu)及相關(guān)類型,以及空間相關(guān)范圍的大小,一般常用模型為球形模型、指數(shù)模型、高斯模型等,這3種模型所對應(yīng)的空間分布均可為聚集分布,其理論模型及公式如下：

球狀模型

指數(shù)模型

高斯模型

不同理論模型的變異函數(shù)方程中,C0為塊金值其大小反映變量的隨機程度。C+C0為基臺值其大小可反映變量變化幅度的大小。a為變程表示以a為半徑的鄰域內(nèi)的任何其它Z(xi)與Z(xi+h)間存在的空間相關(guān)性。C/(C+C0)為空間結(jié)構(gòu)比率,反映了空間結(jié)構(gòu)系數(shù)在樣點變量差異中所占的比率,大小可判定系統(tǒng)內(nèi)變量的空間相關(guān)程度,介于25%—75%之間,有中等程度的空間相關(guān)性,大于75%則具有較強的空間相關(guān)性。

1.4數(shù)據(jù)統(tǒng)計與分析

本文數(shù)據(jù)應(yīng)用DPS軟件的地理統(tǒng)計模塊作地統(tǒng)計學分析。

2　結(jié)果與分析

2.1云斑白條天牛的危害與分布

云斑白條天牛在黃河三角洲地區(qū)主要危害白蠟樹,2年1代,跨3個年度,以幼蟲在樹干內(nèi)部蛀干危害,通過排糞孔排出蟲糞和木屑,第2年秋在樹干內(nèi)化蛹越冬,第3年春天羽化在樹干上咬一羽化孔,鉆出樹體,取食白蠟?zāi)壑l補充營養(yǎng),然后在樹干上咬一產(chǎn)卵刻槽產(chǎn)卵。因此,以產(chǎn)卵刻槽、排糞孔和羽化孔3個指標,分別調(diào)查云斑白條天牛不同蟲態(tài)卵、幼蟲、成蟲(蛹)的種群數(shù)量,其基本統(tǒng)計學特征見表2。由表1和表2的數(shù)據(jù)可見,片林的平均胸徑最小,受害率也最輕,其株均產(chǎn)卵刻槽、排糞孔和羽化孔的數(shù)量均明顯低于城市行道綠化林、庭院綠化林和公路林,這也表明白蠟樹林分類型與云斑白條天牛的危害密切相關(guān),在實際栽植環(huán)境中,城市行道綠化林、庭院綠化林和公路林的立地條件較為類似,而片林栽植密度較大,與其他林分類型差別較大,種群空間分布特征可能存在差異。

2.2卵種群數(shù)量的半變異函數(shù)與空間格局

對4個不同林分類型白蠟林地產(chǎn)卵刻槽的分布進行半變異函數(shù)理論模型擬合,擬合函數(shù)模型及參數(shù)見表3,擬合曲線見圖1。各林分白蠟林地,云斑白條天牛卵種群數(shù)量的半變異函數(shù)擬合模型均為球形模型,半變異函數(shù)的決定系數(shù)R2介于0.70—0.88,平均值為0.75,決定系數(shù)值均較大,說明擬合模型的擬合程度均較高。根據(jù)擬合模型可知,危害白蠟云斑白條天牛卵種群數(shù)量(產(chǎn)卵刻槽)的空間分布有明顯的空間結(jié)構(gòu),均為球形模型,其指明的空間分布型為聚集分布。

表2　危害白蠟云斑白條天牛統(tǒng)計指標調(diào)查數(shù)據(jù)的基本統(tǒng)計參數(shù)

圖1　不同林分白蠟林地云斑白條天牛產(chǎn)卵刻槽的半變差異函數(shù)Fig.1　Variogram of B. lineolata egg niche at different type ash forest

從擬合函數(shù)參數(shù)的塊金值C0看,最小塊金值1.55,表明引起變量的隨機程度較小,最大塊金值為62.58,表明引起變量的隨機程度較大。從變程a看,在8.63—75.08 m之間,其空間依賴范圍均在研究尺度之內(nèi),當間隔距離在變程之內(nèi)時,具有明顯的空間依賴性,當間隔距離超過變程時,半方差函數(shù)趨于穩(wěn)定,不存在空間相關(guān)性。從基臺值C0+C看,最小為5.85,變量的變化幅度較小,最大為85.44,表明變量的變化幅度較大。從空間結(jié)構(gòu)比率C/(C0+ C)看,最小值為19.81%,最大為73.53%,平均為37.59%,介于25%—75%之間,有中等程度的空間相關(guān)性,其中庭院綠化林小于25%,空間相關(guān)性不高。

表3　危害白蠟云斑白條天牛各調(diào)查指標的半變異函數(shù)模型及參數(shù)

*表示差異顯著,**表示差異極顯著

2.3幼蟲種群數(shù)量的半變異函數(shù)與空間格局

對4個不同林分類型白蠟林地排糞孔的分布進行半變異函數(shù)理論模型擬合,擬合函數(shù)模型及參數(shù)見表3,擬合曲線見圖2。各林分白蠟林地,云斑白條天牛幼蟲種群數(shù)量的半變異函數(shù)擬合模型均為球形模型,半變差函數(shù)的決定系數(shù)R2介于0.64—0.85,平均值為0.75,決定系數(shù)值均較大,說明擬合模型的擬合程度均較高。根據(jù)擬合模型可知,危害白蠟云斑白條天牛幼蟲種群數(shù)量(排糞孔)的空間分布有明顯的空間結(jié)構(gòu),均為球形模型,其指明的空間分布型均為聚集分布。

從擬合函數(shù)參數(shù)的塊金值C0看,最小塊金值2.92,表明引起變量的隨機程度較小,最大塊金值為16.40,表明引起變量的隨機程度較大。從變程a看,在15.49—67.02 m之間,其空間依賴范圍均在研究尺度之內(nèi),當間隔距離在變程之內(nèi)時,具有明顯的空間依賴性,當間隔距離超過變程時,半方差函數(shù)趨于穩(wěn)定,不存在空間相關(guān)性。從基臺值C0+C看,最小為3.64,變量的變化幅度較小,最大為24.32,表明變量的變化幅度較大。從空間結(jié)構(gòu)比率C/(C0+C)看,最小值為16.29%,最大為35.81%,平均為26.18%,介于25%—75%之間,有中等程度的空間相關(guān)性,庭院綠化林和片林小于25%,空間相關(guān)性不強。

圖2　不同林分白蠟林地云斑白條天牛排糞孔的半變差異函數(shù)Fig.2　Variogram of B. lineolata frass holes at different type ash forest

2.4成蟲(蛹)種群數(shù)量的半變異函數(shù)與空間格局

對4個不同林分類型白蠟林地羽化孔的分布進行半變異函數(shù)理論模型擬合,擬合函數(shù)模型及參數(shù)見表3,擬合曲線見圖3。各林分白蠟林地,云斑白條天牛成蟲(蛹)種群數(shù)量的擬合模型除片林樣點(為指數(shù)模型)外均為球形模型,半變差函數(shù)的決定系數(shù)R2介于0.74—0.80,平均值為0.77,決定系數(shù)值均較大,說明擬合模型的擬合程度均較高。根據(jù)擬合模型可知,白蠟上云斑白條天牛成蟲(蛹)種群數(shù)量(羽化孔)的空間分布有明顯的空間結(jié)構(gòu),球形模型和指數(shù)模型指明的空間分布型均為聚集分布。片林樣點羽化孔的擬合模型為指數(shù)模型,主要原因為羽化孔種群密度較低僅為0.43個/株,結(jié)合實際調(diào)查的分布情況,數(shù)量不多的羽化孔主要分布于片林周邊林緣的樹木上,處于片林內(nèi)部的樹木基本無分布,導致片林羽化孔的分布較其他林分類型要相對集中。

從擬合函數(shù)參數(shù)的塊金值C0看,最小塊金值0.66,表明引起變量的隨機程度較小,最大塊金值為2.66,表明引起變量的隨機程度較大。從變程a看,在6.95—31.24 m之間,其空間依賴范圍均在研究尺度之內(nèi),當間隔距離在變程之內(nèi)時,具有明顯的空間依賴性,當間隔距離超過變程時,半方差函數(shù)趨于穩(wěn)定,不存在空間相關(guān)性。從基臺值C0+C看,最小為0.88,變量的變化幅度較小,最大為4.37,表明變量的變化幅度較大。從空間結(jié)構(gòu)比率C/(C0+C)看,最小值為18.53%,最大為55.84%,平均為34.68%,介于25%—75%之間,有中等程度的空間相關(guān)性,其中城市行道綠化林小于25%,空間相關(guān)性不強。

圖3　不同林分白蠟林地云斑白條天牛羽化孔的半變差異函數(shù)Fig.3　Variogram of B. lineolata emergence holes at different type ash forest

3　結(jié)論與討論

通過對不同林分白蠟林地調(diào)查數(shù)據(jù)的地統(tǒng)計學分析,表明危害白蠟云斑白條天牛的各個蟲態(tài)在不同林分類型下的空間格局均呈聚集分布,卵、幼蟲、成蟲(蛹)在城市行道綠化林、庭院綠化林、公路林和片林環(huán)境中均表現(xiàn)為聚集分布,聚集分布的習性與其生存的林地環(huán)境關(guān)系不大。

本研究調(diào)查的4種林分類型中, 城市行道綠化林、庭院綠化林和公路林的立地條件較為類似,片林的栽植密度大,環(huán)境條件明顯不同,因此在一些統(tǒng)計指標上與其他3種林分類型表現(xiàn)出較大差異,片林調(diào)查樣點的卵、幼蟲和成蟲(蛹)的變程a、塊金值C0、基臺值C0+C均小于其他3種林分類型的統(tǒng)計值。同時,僅在片林環(huán)境下成蟲(蛹)的半變異函數(shù)模型為指數(shù)模型,除此之外,不同蟲態(tài)在不同林分類型下的半變異函數(shù)模型均為球形模型。這均與片林條件下云斑白條天牛的危害特點和分布情況有關(guān),片林受害率較輕,且分布不均,主要為害在林緣樹木,林地內(nèi)部樹木受害較少。

空間格局是昆蟲種群的重要特征之一,它是由昆蟲種群的生物學特性與特定生境條件協(xié)同進化的結(jié)果。研究昆蟲的空間格局有助于了解昆蟲的生態(tài)特性,對提高抽樣技術(shù)和揭示種群消長規(guī)律、預(yù)測預(yù)報及制訂害蟲防治策略都具有重要意義[8-10]。由于云斑白條天牛一生只有一個排糞孔,以排糞孔來統(tǒng)計其種群數(shù)量是科學的,據(jù)此,通過傳統(tǒng)生物統(tǒng)計學技術(shù)研究表明云斑白條天牛幼蟲在楊樹、桑樹和核桃上的空間分布型為聚集分布[11-14],但傳統(tǒng)生物統(tǒng)計方法忽略了研究對象的空間位置,且對取樣單位的大小依賴較大,而地統(tǒng)計學能夠最大限度地利用野外調(diào)查所提供的各種信息,不僅可對區(qū)域化變量的空間相關(guān)性進行定量化描述,還可對空間分布進行模擬,能獲得更多的統(tǒng)計結(jié)果[3,15-18]。本研究通過地統(tǒng)計學的方法,除獲得了不同林分類型白蠟林地云斑白條天牛卵、幼蟲、成蟲(蛹)的空間分布型外,還得出各個蟲態(tài)的空間擬合模型、空間依賴性的變程等信息,對指導云斑白條天牛蟲情測報和防治更具指導價值。

[1]楊忠岐, 李建慶, 梅增霞, 王小藝, 景志高, 孟向東. 釋放花絨寄甲防治危害白蠟的云斑天牛. 林業(yè)科學, 2011, 47(12): 78-84.

[2]梅增霞, 李建慶. 黃河三角洲地區(qū)云斑白條天牛的發(fā)生危害與防治. 濱州學院學報, 2014, 30(3): 58-62.

[3]潘杰, 王濤, 宗世祥, 溫俊寶, 駱有慶. 紅脂大小蠹種群空間格局地統(tǒng)計學分析及抽樣技術(shù). 生態(tài)學報, 2011, 31(1): 195-202.

[4]Yünsel T Y, Ersoy A, ?etin M. Geostatistical analysis of spatial distribution of salt bed thickness. Journal of Mining Science, 2002, 38(6): 565-572.

[5]李友常, 夏乃斌, 屠泉洪, 駱有慶, 溫俊寶. 楊樹光肩星天牛種群空間格局的地統(tǒng)計學研究. 生態(tài)學報, 1997, 17(4): 393-401.

[6]柳林俊. 雁北地區(qū)青楊天牛空間分布型的地統(tǒng)計學研究. 陜西林業(yè)科技, 2005, (3): 28-30, 33-33.

[7]唐啟義, 馮明光. DPS數(shù)據(jù)處理系統(tǒng)——實驗設(shè)計、統(tǒng)計分析及數(shù)據(jù)挖掘. 北京: 科學出版社, 2007: 504-528.

[8]丁巖欽. 昆蟲數(shù)學生態(tài)學. 北京: 科學出版社, 1994: 22-69.

[9]Wang Z Q. Application of Geostatistics in Ecology. Beijing: Science Press, 1999: 10-72.

[10]Hohn M E, Liebhold A M, Gribko L S. Geostatistical model for forecasting spatial dynamics of defoliation caused by the gypsy moth (Lepidoptera: Lymantriidae). Environmental Entomology, 1993, 22(5): 1066-1075.

[11]李建慶, 楊忠岐, 張雅林, 梅增霞. 楊樹上云斑天牛種群的空間格局及抽樣技術(shù). 昆蟲學報, 2009, 52(8): 860-866.

[12]梅愛華. 云斑天牛幼蟲空間分布型及抽樣技術(shù). 昆蟲知識, 1997, 34(2): 94-95.

[13]張世權(quán), 楊寶祥, 鄭麗芳. 云斑天牛空間分布型與種群密度估計的研究. 河北林學院學報, 1992, 7(3): 210-213.

[14]吳開明, 張建強, 代方銀, 許恩遠, 唐萬成, 邵春明, 鞠先奎. 云斑天牛危害桑樹及生物學特性研究. 蠶業(yè)科學, 1995, 21(1): 53-54.

[15]Farley R A, Fitter A H. Temporal and spatial variation in soil resources in a deciduous woodland. Journal of Ecology, 1999, 87(4): 668-696.

[16]Liebhold A M, Rossi R E, Kemp W P. Geostatistics and geographic information systems in applied insect ecology. Annual Review of Entomology, 1993, 38(1): 303-327.

[17]孫鵬舉, 駱有慶, 姚東華, 劉春亮, 宗世祥. 柳毒蛾幼蟲和蛹的空間分布的地統(tǒng)計學分析. 植物保護學報, 2015, 42(4): 604-611.

[18]蔣杰賢, 萬年峰, 季香云. 基于地統(tǒng)計學和GIS的化學農(nóng)藥脅迫下桃樹桃蚜與草間鉆頭蛛種群空間格局. 中國生態(tài)農(nóng)業(yè)學報, 2015, 23(7): 906-913.

Geostatistical analysis on the spatial distribution ofBatoceralineolata(Coleoptera: Cerambycidae) populations in different types of ash tree forests

LI Jianqing1, MEI Zengxia1, YANG Zhongqi2,*

1ShandongProvinceKeyLaboratoryonEco-environmentsofYellowRiverDelta,BinzhouUniversity,Binzhou256603,China2KeyLaboratoryofForestProtectionofStateForestryAdministration,ResearchInstituteofForestEcology,EnvironmentandProtection,ChineseAcademyofForestry,Beijing100091,China

Batoceralineolatais an important woodorer pest in China, mainly attacking the trunk of ash treesFraxinusvelutina, a predominant forest species in the Yellow River Delta, causing severe damage. In the present study, the spatial distribution patterns of eggs, larvae, and adults or pupae of theB.lineolatapopulation were analyzed using geostatistical methods. The egg niche, frass hole, and emergence hole at four ash tree sample plots, including a city road greening forest, yard greening forest, highway forest, and patch forest, were investigated. The results showed that the semivariogram curve models of eggs, larvae, and adults (pupae) ofB.lineolatapopulations in different types of ash forests were all spherical, except the adults (pupae) in the patch forest plot were exponential, and the spatial distribution patterns were all aggregated. There were no obvious correlations between the aggregated behavior and living habitat ofB.lineolatapopulations. However, the geostatistical parameters of the spatial correlation ofB.lineolatapopulations in different types of ash forests were obviously different. The range (a), nugget (C0), sill (C0+C) of eggs, larvae, and adults (pupae) in the patch forest were lower than those of the other three types of forests. At the same time, the semivariogram model of adults (pupae) at the patch forest was exponential.

Batoceralineolata; spatial distribution; geostatistics; ash trees

國家“十一五”科技支撐計劃項目(2006BAD08A12)；山東省優(yōu)秀中青年科學家科研獎勵基金項目(BS2013NY014); 山東省自然科學基金項目(ZR2014CL031,ZR2015CL043); 山東高等學校科技計劃項目(J13LF10)

2015-06-29;

2015-10-22

Corresponding author.E-mail: yzhqi@caf.ac.cn

10.5846/stxb201506291334

李建慶,梅增霞,楊忠岐.不同林分白蠟樹云斑白條天牛種群空間格局地統(tǒng)計學分析.生態(tài)學報,2016,36(14):4540-4547.

Li J Q, Mei Z X, Yang Z Q.Geostatistical analysis on the spatial distribution ofBatoceralineolata(Coleoptera: Cerambycidae) populations in different types of ash tree forests.Acta Ecologica Sinica,2016,36(14):4540-4547.