姜彩鴿，王國珍*，韓江霞，樊仲慶

(1.寧夏農(nóng)林科學(xué)院植物保護研究所，銀川 750002； 2.寧夏小任果業(yè)發(fā)展有限公司，永寧 750100)

設(shè)施葡萄上薊馬種群消長規(guī)律和空間分布

姜彩鴿1，王國珍1*，韓江霞2，樊仲慶1

(1.寧夏農(nóng)林科學(xué)院植物保護研究所，銀川 750002； 2.寧夏小任果業(yè)發(fā)展有限公司，永寧 750100)

為摸清寧夏設(shè)施葡萄上薊馬種群消長規(guī)律和空間分布，于2011年進(jìn)行了調(diào)查研究。從各項聚集度指標(biāo)、空間分布型適合度卡方(χ2)檢驗及單因素隨機區(qū)組設(shè)計分析看，薊馬在設(shè)施葡萄上均呈聚集分布，棚內(nèi)北、南間差異極顯著，蟲口分布具有北側(cè)聚集的特點。以蟲態(tài)來說，在整個葡萄生育期除3月下旬一次大量孵化時期以若蟲為害為主外，其他時期均以成蟲為害為主，尤其在種群穩(wěn)定階段，3月17日左右成蟲蟲口數(shù)量達(dá)到高峰，若蟲也在3月24日左右急劇增加達(dá)到高峰；4月初之后若蟲消退，以少量成蟲為害。從蟲口量比例看，薊馬種群的空間分布，呈現(xiàn)出在大棚北側(cè)、中上層居多的特點。

設(shè)施葡萄； 薊馬； 種群消長； 空間分布型

葡萄薊馬又稱煙薊馬(ThripstabaciLindeman)，是纓翅目(Thysanoptera)薊馬科(Thripidae)昆蟲。中國記錄薊馬種類有336種[1]，主要以成蟲和若蟲的銼吸式口器吸取寄主幼嫩組織的汁液為生，進(jìn)食時會造成葉片、花蕾及果實直接受害。據(jù)報道，為害果樹的薊馬有40多種[2]，且在我國已有廣泛的分布。

寧夏地處西北中部，因具有充足的光熱資源，近年來設(shè)施葡萄的發(fā)展規(guī)模逐年擴大，尤其是促早栽培葡萄，薊馬的危害呈現(xiàn)逐年加重的趨勢，而此蟲的發(fā)生直接關(guān)系到葡萄的產(chǎn)量、品質(zhì)和產(chǎn)品的經(jīng)濟價值。國內(nèi)有關(guān)薊馬對葡萄的為害狀況及其時空分布方面的研究資料較少，而寧夏尚屬空白，寧夏設(shè)施(日光溫室)葡萄薊馬防治技術(shù)可借鑒的資料不多，因此在初步了解此害蟲對設(shè)施葡萄為害的基礎(chǔ)上，2011年制定了詳盡的方案研究其種群的發(fā)展動態(tài)及空間分布，對制定科學(xué)的防治策略，明確防治關(guān)鍵期，有效控制薊馬對寧夏設(shè)施(日光溫室)葡萄的為害，具有重要意義。

1 材料與方法

1.1 調(diào)查品種、地點和葡萄棚管理

在寧夏永寧縣選擇小任果業(yè)設(shè)施促早栽培葡萄有代表性的25號溫棚，品種為‘紅地球’，樹齡5年，單籬架栽培，共有葡萄植株64行。此棚采用無公害的病蟲害防治措施，農(nóng)藥均在推薦濃度下使用，管理良好。葡萄薊馬多年發(fā)生偏重。

1.2 調(diào)查方法

2011年3月3日到5月19日，每隔7 d調(diào)查一次，共調(diào)查12次。

采用階層抽樣方法調(diào)查，在棚內(nèi)由東到西選樹勢一致的6行葡萄樹，每行由北向南均勻分4個方位各選1株(即由北向南所調(diào)查的4個方位均設(shè)6個重復(fù))，共24株葡萄，每株分上(1.5 m以上)、中(1～1.5 m)、下(1m以下)3層，每層掛牌調(diào)查具代表性的2個枝條和全部花序，輕輕拍打至藍(lán)板上，記錄薊馬所在部位的蟲口數(shù)量及蟲態(tài)。同時調(diào)查受害狀、葡萄果粒出現(xiàn)害狀的時間、受害程度(果穗數(shù)、果粒總數(shù)、受害穗數(shù)、受害粒數(shù))。

1.3 分析方法

利用DPS數(shù)據(jù)處理系統(tǒng)進(jìn)行種群空間分布型聚集度指標(biāo)測定[3-8]。

計算公式如下:

將調(diào)查獲得的數(shù)據(jù)列成頻次分布表，進(jìn)行頻次分布檢驗,求得各種空間分布型的理論頻次,與實測頻次相比較，利用DPS數(shù)據(jù)處理系統(tǒng)進(jìn)行卡方檢驗，確定其所屬分布型[3,9-10]。

利用DPS數(shù)據(jù)處理系統(tǒng)進(jìn)行單因素隨機區(qū)組設(shè)計分析[4,11]。

計算每次調(diào)查葡萄棚的有蟲株數(shù)、成若蟲蟲口數(shù)及各個方位的有蟲株率和蟲口量比例，用Excel軟件繪制折線及柱形動態(tài)圖。其中有蟲株數(shù)、成若蟲蟲口數(shù)分別用以衡量薊馬隨時間變化的發(fā)生動態(tài)；有蟲株率是指某一方位(如溫棚南側(cè))發(fā)生薊馬的植株數(shù)量占調(diào)查植株總數(shù)的比例，蟲口量比例是指某一方位(如葡萄植株上部)的薊馬發(fā)生量占調(diào)查的薊馬發(fā)生總量的比例，兩者分別用來衡量薊馬在發(fā)生率和發(fā)生量上的空間分布結(jié)構(gòu)[12-18]。

2 結(jié)果分析

2.1 薊馬在設(shè)施葡萄上的空間分布型

2.1.1 薊馬空間分布的主要參數(shù)

取3月17日落花期薊馬種群發(fā)展階段所調(diào)查的數(shù)據(jù)進(jìn)行分析，計算其空間分布的主要參數(shù)，結(jié)果見表1。

表1設(shè)施葡萄薊馬聚集度的測定結(jié)果1)

Table1Aggregationindicesofgrapethripsonfacilitiesgrapes

方位Position擁擠度M?I指標(biāo)M?/M指標(biāo)CA指標(biāo)擴散系數(shù)CK指標(biāo)空間分布型Spatialdistributionpattern131．977214．69151．84990．849915．69151．1766聚集分布Aggregateddistribution211．93044．21611．54650．54655．21611．8297聚集分布Aggregateddistribution36．23472．09181．50490．50493．09181．9805聚集分布Aggregateddistribution42．01800．87511．76560．76561．87511．3061聚集分布Aggregateddistribution

1) 方位1、2、3、4分別表示由北向南依次4個方位的各重復(fù)植株，對其上蟲口數(shù)量進(jìn)行調(diào)查后計算各項指標(biāo)。 Position 1,2,3and 4 indicate the repeated plants in 4 directions from north to south,respectively.Each indicator is calculated based on the number of its population.

經(jīng)對取樣行由北向南依次4個方位的各重復(fù)植株上蟲口數(shù)量的調(diào)查結(jié)果計算得出，其M*值、M*/M值、C值均大于1，I值、CA值、K值均大于0，空間分布均表現(xiàn)為聚集分布。

2.1.2 Iwao M*-M回歸模型

運用4組調(diào)查數(shù)據(jù)可以建立薊馬密度M和平均擁擠度M*間的回歸方程為：M*=-1.232 73+1.885 08M。

相關(guān)系數(shù)r=0.9966，說明線性相關(guān)極為顯著，線性方程可靠。在該方程中，α=-1.232 73<0，說明其個體之間相互排斥；β=1.885 08>1，說明其為聚集分布。

2.1.3 Taylor冪法則

運用4組調(diào)查數(shù)據(jù)同樣可以建立很有意義的處理平均數(shù)與方差的對數(shù)值之間的冪法則回歸方程：

lgv= 0.14037+1.755 19 lgM

相關(guān)系數(shù)r=0.9904，說明線性相關(guān)極為顯著，線性方程可靠。在該方程中，b=1.755 19>1，也說明其為聚集分布型。

2.2 空間分布型適合度卡方(χ2)檢驗

將3月17日落花期薊馬種群發(fā)展階段所調(diào)查的4組數(shù)據(jù)列成頻次分布表，進(jìn)行卡方檢驗，計算結(jié)果如下：

二項分布頻次檢驗:χ2=6609.9038，p=0.0000；

波松分布頻次檢驗: χ2=3027.6069，p=0.0000；

負(fù)二項分布頻次檢驗：

矩法估計參數(shù)：K=0.609 7，χ2=1.274 0，p=0.865 8；

樣本中零頻率較小，不適合用零頻率法估計參數(shù)K；

最大或然法估計參數(shù)：K=0.6907，χ2=1.115 3，p=0.8918；

Neyman分布頻次檢驗:χ2=5.2864，p=0.152 0；

波松-二項分布頻次檢驗:N=40，P=0.353952 1，χ2=175.6164，p=0.0000；

復(fù)合波松分布頻次檢驗：MM=1.7918，A=0.8709，χ2=23.445 7，p=0.0000；

由以上結(jié)果可看出，各種空間分布型中均以負(fù)二項分布的χ2最小；在負(fù)二項分布中，又以最大或然法估計參數(shù)的擬合結(jié)果最好，因此最適合用負(fù)二項分布來表達(dá)其空間分布型(服從負(fù)二項分布的生物種群在空間結(jié)構(gòu)上呈聚集分布，特點是呈疏密相間的極不均勻狀或嵌紋狀)，即薊馬在設(shè)施葡萄上呈聚集分布。

2.3 單因素隨機區(qū)組設(shè)計處理間差異分析

運用3月17日落花期薊馬4組調(diào)查數(shù)據(jù)進(jìn)行單因素隨機區(qū)組設(shè)計分析，得差異顯著結(jié)果(表2)如下。

表2新復(fù)極差檢驗結(jié)果1)

Table2Newmultiplerangetestresults

方位Position蟲口數(shù)量/頭Population1(20．00±1．15)aA2(8．67±1．28)abAB3(4．33±1．58)bAB4(0．67±0．33)bB

1) 方位1、2、3、4分別表示由北向南依次4個方位的各重復(fù)植株。表中數(shù)據(jù)為平均數(shù)±標(biāo)準(zhǔn)誤。同列數(shù)據(jù)后不同大、小寫字母表示經(jīng)Duncan氏新復(fù)極差法檢驗差異顯著(P<0.01,P<0.05)。Position 1,2,3and 4 indicate the repeated plants in 4 directions from north to south,respectively.The data in the table are mean±SE.Different lowercase and uppercase letters within the same column indicate significant difference atP<0.05 andP<0.01levels by Duncan’s new multiple range test,respectively.

從表2結(jié)果可看出，相鄰的方位間無顯著差異。在5%顯著水平上，方位1與方位2間無顯著差異，而與方位3、4有顯著差異。在1%極顯著水平上，方位1與方位2、3無顯著差異，而與方位4有顯著差異?？偟膩碚f，不管是在何種顯著水平下，蟲口分布由北向南呈遞減趨勢；此外，方位1與方位4，即北、南植株蟲口數(shù)量間差異極顯著，蟲口分布具有北側(cè)聚集的特點。

2.4 薊馬種群空間分布動態(tài)

2.4.1 種群消長規(guī)律

將開花期至著色期調(diào)查的數(shù)據(jù)匯總分析，得種群消長規(guī)律結(jié)果(圖1～圖3)。

根據(jù)圖1所表現(xiàn)的變化幅度，可將設(shè)施葡萄薊馬種群消長動態(tài)劃分為3個階段：(1)2011年3月初開花期蟲口基數(shù)已經(jīng)很高，蟲口數(shù)量101頭/株，植株受害率已達(dá)87.50%；至3月中旬落花期蟲口密度穩(wěn)步增加，為種群發(fā)展階段，蟲口數(shù)量發(fā)展到202頭/株，增長了1倍，植株受害率達(dá)91.67%。(2)3月中旬落花期至3月底幼果期，薊馬對設(shè)施栽培‘紅地球’的為害呈穩(wěn)定狀態(tài)，植株受害率保持在91.67%，蟲口數(shù)量達(dá)到高峰395頭/株。(3)4月初幼果期至5月下旬著色期，蟲口數(shù)量受溫度波動及修剪施藥等管理措施影響大幅度下降，為種群衰退階段，蟲口數(shù)量下降至1頭/株，植株受害率4.17%，薊馬的為害減輕很多。

圖1 薊馬蟲口數(shù)量及為害植株動態(tài)變化Fig.1 Dynamics of the number of thrips population and injured plants

圖2 不同聚集部位蟲口數(shù)量動態(tài)變化Fig.2 Dynamics of population quantity in different aggregation sites

圖3 果穗/果粒受害率動態(tài)變化Fig.3 The dynamics of damage rate of ears and fruits

以蟲態(tài)來說，在整個葡萄生育期除3月下旬一次大量孵化時期以若蟲為害為主外，其他時期均以成蟲為害為主。尤其在種群穩(wěn)定階段，3月17日左右成蟲蟲口數(shù)量達(dá)到高峰，蟲口密度達(dá)到195頭/株；若蟲也在3月24日左右急劇增加達(dá)到高峰，最高蟲口密度達(dá)到了259頭/株。4月初之后若蟲消退，以少量成蟲為害。

由圖2可見，3月末即坐果期前薊馬多在花序或果穗上活動為害，坐果期后因新生幼嫩枝條增多及疏果管理措施影響而轉(zhuǎn)移到葉片為害。其為害動態(tài)與圖1相符。

由圖3可看出，薊馬對葡萄整個生育期的為害是一個不斷疊加的過程，其果穗受害率及果粒受害率不斷加大，至著色期套袋時分別已達(dá)到92%和14.76%，直接影響了鮮食葡萄的品質(zhì)和價值。

2.4.2 空間動態(tài)特點

從蟲口量比例看，隨著時間的推移，薊馬種群顯現(xiàn)出由北到南(圖4)、從中層到上下層(圖5)轉(zhuǎn)移的趨勢。在整個葡萄生育期，薊馬主要在大棚北側(cè)偏涼處活動，只有在降溫涼爽氣候下才稍向南側(cè)轉(zhuǎn)移，但仍以北側(cè)為多(圖4)；而在垂直方向上，5月初著色初期之前多分布在嫩梢嫩葉較多的上層和花果較集中的中層取食為害，隨著時間推移，5月上旬由于植株中層果穗套袋和下層修剪枝條留出通風(fēng)帶及葉片老化，薊馬就轉(zhuǎn)移至嫩梢較集中的上層為害(圖5)。

圖4 溫棚水平方向蟲口量比例動態(tài)Fig.4 Dynamics of the proportion of population quantity in horizontal direction in greenhouses

圖5 溫棚垂直方向蟲口量比例動態(tài)Fig.5 Dynamics of the proportion of population quantity in vertical direction in greenhouses

從有蟲株率看，薊馬空間分布動態(tài)也顯現(xiàn)了大致相同的趨勢(圖6、7)。但從整個發(fā)生時期看，不管是蟲口量比例還是有蟲株率，大棚北側(cè)和植株中上層的發(fā)生程度均較嚴(yán)重。

圖6 溫棚水平方向有蟲株率動態(tài)Fig.6 Dynamics of infestation rate in horizontal direction in greenhouses

圖7 溫棚垂直方向有蟲株率動態(tài)Fig.7 Dynamics of infestation rate in vertical direction in greenhouses

3 討論

3.1 設(shè)施葡萄薊馬來源

葡萄薊馬在整個葡萄生育期前后世代重疊，多以成蟲和若蟲在葡萄園附近的雜草、蔬菜和死株上越冬，少數(shù)以蛹越冬[19]。第二年在葡萄初花期，開始發(fā)現(xiàn)有薊馬為害。2011年3月初大棚內(nèi)已發(fā)現(xiàn)有薊馬活動，有蟲株率達(dá)87.50%左右，危害極大極早。永寧地區(qū)2月下旬溫度已在10℃以上，此溫度下薊馬可成功繁殖，其設(shè)施葡萄于2月下旬開始開花，因此應(yīng)在冬季及扣棚前及時有效地清理棚內(nèi)外的枯枝殘葉等越冬寄主，減少蟲口基數(shù)，做好防治措施。

3.2 薊馬種群在大棚內(nèi)空間分布的特點

設(shè)施葡萄上薊馬的空間分布型總體上屬于聚集分布，呈現(xiàn)出在大棚北側(cè)、中上層居多的特點。原因可能是南側(cè)靠近大棚拱腳，薄膜低，空間狹窄空氣流通差，同時溫度回升快、偏高；北側(cè)空間大，空氣流通好，偏涼；而薊馬適宜在25 ℃以下活動，因此多聚集于北側(cè)。植株中上層多果穗和嫩梢，薊馬喜食，多活動于此；而植株下層一般為防治病蟲害而修剪枝條作通風(fēng)帶，留下少量老葉不利薊馬取食。利用大棚內(nèi)薊馬種群空間分布的特點，控制薊馬種群發(fā)展值得深入研究。

3.3 為害程度調(diào)查及癥狀鑒定

薊馬對葡萄的為害，最主要表現(xiàn)在對果粒的為害程度上，影響其外觀、品質(zhì)及產(chǎn)量，今后應(yīng)完善薊馬對葡萄果粒為害程度調(diào)查和受害癥狀鑒定的研究，精確掌握蟲口密度、發(fā)展動態(tài)，有利于做好預(yù)測預(yù)報和防治措施。

3.4 設(shè)施葡萄薊馬的防治

3.4.1 農(nóng)業(yè)防治

避免在保護地種植油桃、蔥等替代寄主植物，以免為薊馬越冬、越夏創(chuàng)造有利條件，同時入冬后及時清除落葉和雜草，減少越冬蟲源。

3.4.2 生物防治

在田間調(diào)查中，發(fā)現(xiàn)有絨螨，且當(dāng)其數(shù)量較多時，薊馬蟲口數(shù)量有所減少，但這種絨螨是否是其天敵還有待于進(jìn)一步研究鑒定。與薊馬對應(yīng)的天敵種群消長動態(tài)、保護和提高棚內(nèi)天敵數(shù)量進(jìn)行有效防治是值得關(guān)注和研究的方向。

3.4.3 藥劑防治

據(jù)調(diào)查，薊馬為害葡萄以果實為主，因此噴藥時期應(yīng)在開花前1～2 d或初花期進(jìn)行，以降低蟲口基數(shù)；此外，還應(yīng)對薊馬活動較集中的大棚北側(cè)、樹冠中上層部位進(jìn)行重點防治。

[1] Zhang W Q,Tong X L.Checklist of thrips (Insecta:Thysanoptera) from China[J].Pure and Applied Zoology,1993,4:409-443.

[2] 韓運發(fā).中國經(jīng)濟昆蟲志(第五十五冊纓翅目)[M].北京:科學(xué)出版社,1997:45-287.

[3] 丁巖欽.昆蟲種群生態(tài)學(xué)原理與應(yīng)用[M].北京:科學(xué)出版社,1994:84-124.

[4] 徐汝梅.昆蟲種群生態(tài)學(xué)[M].北京: 北京師范大學(xué)出版社,1985:7-13.

[5] 鄭伯平,鄭長英,顧松東,等.白三葉草生長期西花薊馬的空間分布型變化[J].環(huán)境昆蟲學(xué)報,2010,32(1):1-5.

[6] Arzone A,Alma A,Rapetti S.Frankliniellaoccidentalis(Perg.) (Thysanoptera: Thripidae)nuovo fitomizo delleserre in Italia[J].Informatore Fitopatologico,1989,39(10):43-48.

[7] Gerin C,Hance T,van Impe G.Impact of flowers on the demography of western flower thripsFrankliniellaoccidentalis(Pergande) (Thysanoptera: Thripidae)[J].Journal of Applied Entomology,1999,123:569-574.

[8] Hoddle M S,Robinson L,Morgan D.Attraction of thrips(Thysanoptera: Thripidae and Aeolothripidae)to colored sticky cards in a California avocado orchard[J].Crop Protection， 2002,21:383-388.

[9] Hulshof J,Ketoja E,Vanninen L.Life history characteristics ofFrankliniellaoccidentalison cucumber leaves with and without supplemental food[J].Entomologia Experimentalis et Applicata,2003,108(1):19-32.

[10]Kirk W D J,Terry L I.The spread of the western flower thripsFrankliniellaoccidentalis(Pergande)[J].Agricultural and Forest Entomology， 2003,5:301-310.

[11]Kirk W D J.Pollen feeding in thrips(Insecta:Thysanoptera)[J].Bangladesh Journal of Zoology,1984,20(4):107-117.

[12]侯茂林，文吉輝，盧偉.煙粉虱成蟲在日光溫室內(nèi)的分布和日活動規(guī)律[J].生態(tài)學(xué)報,2006,26(5):1431-1437.

[13]黃頂成，張潤志，王斌,等.保護地桃蚜種群時空動態(tài)與防治對策[J].植物保護,2009,35(6):108-113.

[14]Lewis T.Thrips,their biology,ecology and economic importance[M].London:Press Academic,1973:349.

[15]劉麗輝,張帆,吳珍泉.溫度對西花薊馬(Frankliniellaoccidentalis)生長發(fā)育和存活率的影響[J].生態(tài)學(xué)報,2008,28(10):4891-4895.

[16]路虹,宮亞軍,石寶才,等.西花薊馬在黃瓜和架豆上的空間分布型及理論抽樣數(shù)[J].昆蟲學(xué)報,2007,50(11):1187-1193.

[17]Trichilo P J,Leigh T F.Influence of resource quality on the reproductive fitness of flower thrips(Thysanoptera:Thripidae)[J].Annals of the Entomological Society of America,1988,81(1):64-70.

[18]Van Rijin P C J,Sabelis M W.Does alternative food always enhance biological control the effect of pollen on the interaction between western flower thrips and its predators[J].Bulletin OILB/SROP， 1993,16(8): 123-125.

[19]鄭玉梅,王福利.葡萄薊馬的發(fā)生與防治[J].北京農(nóng)業(yè),2003(9):20.

Populationfluctuationandspatialdistributionofthethripsongrapesinprotectivefacilities

Jiang Caige1, Wang Guozhen1, Han Jiangxia2, Fan Zhongqing1

(1.InstituteofPlantProtection,NingxiaAcademyofAgriculturalandForestrySciences,Yinchuan750002,China;2.NingxiaXiaorenFruitIndustryDevelopmentCo.,Ltd.,Yongning750100,China)

Population fluctuation and spatial distribution of the thrips on grapes in protective facilities was investigated in Ningxia in 2011.From various aggregation indices,chi-square test of spatial distribution pattern and one-factor randomized block design analysis,the thrips showed aggregated distribution on grapes in facilities,showing significant difference between the south and the north in the greenhouse,and insect distribution characteristics was aggregated in the northern side of the greenhouse.The adults caused damage throughout grape growing period,except that,in late March,a large number of hatched nymphs caused harm.The adult population reached its peak around March 17 and nymphs dramatically increased and reached its peak at March 24.Nymph population subsided in early April,and a small number of adults caused harm.According to the proportion of population quantity,the thrips were mainly distributed in the upper north side in the greenhouses.

grapes in facilities; thrips; population fluctuation; spatial distribution pattern

2013-03-02

：2013-06-08

寧夏自然科學(xué)基金項目(NZ1168)；現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系專項資金(CARS-30)

S 436.631.2

：ADOI：10.3969/j.issn.0529-1542.2014.02.027

* 通信作者 E-mail: wangguozhen00@sina.com