余杰穎 任軻亮 薛文鵬 李添群 汪雪俊 耿坤

摘??? 要：為掌握貴陽市獼猴桃園桔小實蠅種群發(fā)生動態(tài)，找出影響種群動態(tài)變動的氣象因子，為桔小實蠅的監(jiān)測和防控提供科學(xué)指導(dǎo)。2017—2018年，在貴陽市主要獼猴桃基地獼猴桃園，采用Steiner誘捕器內(nèi)添加桔小實蠅誘劑對桔小實蠅種群進(jìn)行系統(tǒng)調(diào)查與監(jiān)測，并對影響種群變動的主要氣象因子進(jìn)行分析。結(jié)果表明，1—6月，貴陽市獼猴桃園桔小實蠅種群發(fā)生數(shù)量較低， 7月中旬開始上升，8月中旬至10月為種群盛發(fā)期，11月開始急劇下降，12月種群逐漸消退。桔小實蠅種群數(shù)量的日變化與日平均氣溫、日最低氣溫相關(guān)性極顯著，與日最高氣溫相關(guān)性顯著。日平均氣溫、日最低氣溫、日最高氣溫、日降雨量對種群的影響比日平均相對濕度更直接。氣溫是影響桔小實蠅種群變動的主要因素，受氣溫、降雨量以及其他因素的綜合影響，每年8月中旬至10月為貴陽市獼猴桃園桔小實蠅種群發(fā)生盛期，在8月下旬或9月中旬桔小實蠅種群數(shù)量達(dá)到最高峰。

關(guān)鍵詞：獼猴桃園;桔小實蠅;種群監(jiān)測;影響因子

中圖分類號：Q968.1;S436.611.2???????????? 文獻(xiàn)標(biāo)志碼：A?????????? DOI 編碼：10.3969/j.issn.1006-6500.2021.12.010

Population Dynamics and Influencing Factor of Bactrocera dorsalis in Kiwifruit Orchard in Guiyang City

YU Jieying1， REN Keliang1， XUE Wenpeng1， LI Tianqun2， WANG Xuejun2， GENG Kun1

（1. Plant Protection and Quarantine Station of Guiyang，? Guiyang， Guizhou 550081， China; 2. The Bureau of Agriculture and Rural of Xiuwen，? Guiyang，Guizhou 550200， China）

Abstract： To grasp the population dynamics of Bactrocera dorsalis in Guiyang， find out the meteorological factors that affect the population dynamics， and provide scientific guidance for the monitoring and prevention of Bactrocera dorsalis. From 2017 to 2018， in the kiwi orchard of the main kiwi fruit base in Guiyang City， the Steiner trap was used to add the Bactrocera dorsalis attractant to systematically investigate and monitor the population of Bactrocera dorsalis， and the main meteorological factors affecting population changes Perform analysis. The results showed that， From January to June， the population of Bactrocera dorsalis in kiwifruit orchards in Guiyang City was relatively low. It began to increase in mid-July. The population peaked from mid-August to October. The population began to decline sharply in November and gradually decreased in December. Subside. The daily variation of Bactrocera dorsalis population has a very significant correlation with daily average temperature and daily minimum temperature， and significant correlation with daily maximum temperature. The daily average temperature， daily minimum temperature， daily maximum temperature， and daily rainfall have a more direct impact on the population than the daily average relative humidity. Temperature is the main factor affecting the population changes of Bactrocera dorsalis. Due to the comprehensive influence of temperature， rainfall and other factors， mid-August to October each year is the peak period for the population of Bactrocera dorsalis in kiwifruit orchards in Guiyang. The population of Bactrocera dorsalis reached its highest peak in late August or mid-September.

Key words：kiwifruit orchard;Bactrocera dorsalis;population monitoring;influencing factor

收稿日期：2021-08-17

基金項目：2021年貴陽市植物疫病防控項目（筑財農(nóng)[2021]10號）

作者簡介：余杰穎（1983—），男，貴州甕安人，高級農(nóng)藝師，碩士，主要從事植物保護(hù)技術(shù)方面研究。

通訊作者簡介：耿坤（1976—），女，山東兗州人，研究員，博士，主要從事植物保護(hù)技術(shù)方面研究。

桔小實蠅（Bactrocera dorsalis）又名東方果實蠅，屬雙翅目（Diptera）實蠅科（Tephritidae）寡毛實蠅亞科（Dacinae）離腹寡毛實蠅屬（Bactrocera）[1-2]。1911年首次發(fā)現(xiàn)于臺灣的柑桔上，目前廣泛分布于全球83個國家，被許多國家列為重要檢疫性害蟲。其可危害世界范圍內(nèi)的70科305種植物[3]，主要以幼蟲在果實內(nèi)取食，引起組織腐爛、落果，造成嚴(yán)重的經(jīng)濟(jì)損失。不同種植區(qū)域生態(tài)環(huán)境存在一定的差異，桔小實蠅的發(fā)生規(guī)律也存在一定差異。目前，國內(nèi)許多學(xué)者已對不同生境下桔小實蠅種群動態(tài)及影響因子進(jìn)行了大量研究，2014年林明光等[4]對海南番石榴園桔小實蠅種群動態(tài)與綜合防治進(jìn)行了研究報道;賀海雄等[5]2016年對黔南州都勻市柑桔園實蠅發(fā)生種群及消長規(guī)律進(jìn)行了研究報道;2019年朱俐遐等[6]對攀枝花地區(qū)杧果園桔小實蠅成蟲種群動態(tài)規(guī)律進(jìn)行了研究報道;周雙云等[7]2020年對廣西南寧市蓮霧上桔小實蠅的發(fā)生規(guī)律及產(chǎn)卵選擇性進(jìn)行了研究報道;黃慧欣等[8]對廣西南寧市火龍果產(chǎn)區(qū)桔小實蠅成蟲消長規(guī)律研究進(jìn)行了研究報道;2021年巫輔民等[9]對廣西柿主產(chǎn)區(qū)桔小實蠅的轉(zhuǎn)主為害規(guī)律進(jìn)行了研究報道，均為當(dāng)?shù)亟坌嵪壍谋O(jiān)測與防控提供了重要信息和依據(jù)。

近年來，貴陽市大力開展產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，獼猴桃作為特色果業(yè)得到了大力發(fā)展，截止2020年底，已發(fā)展獼猴桃產(chǎn)業(yè)1.11萬hm-2，對貴陽市經(jīng)濟(jì)社會發(fā)展，果農(nóng)增收致富發(fā)揮了至關(guān)重要的作用。隨著獼猴桃種植面積的不斷擴(kuò)大，桔小實蠅逐漸上升為主要害蟲，危害日趨嚴(yán)重，嚴(yán)重影響獼猴桃產(chǎn)業(yè)可持續(xù)健康發(fā)展。本研究于2017—2018年對貴陽市主要獼猴桃基地獼猴桃園桔小實蠅種群發(fā)生動態(tài)進(jìn)行系統(tǒng)監(jiān)測，旨在摸清桔小實蠅種群發(fā)生動態(tài)及影響種群變動的主要氣象因子，為貴陽市獼猴桃園桔小實蠅的監(jiān)測及防控提供科學(xué)依據(jù)。

1 材料和方法

1.1 材料

誘捕器：Steiner誘捕器（簡稱S），誘劑：桔小實蠅性誘劑Methyleugenol（簡稱Me），均由貴州省出入境檢驗檢疫局提供。

1.2 方法

2017—2018年，在貴陽市主要獼猴桃基地修文縣龍場鎮(zhèn)新寨村牛角土利用Steiner誘捕器與桔小實蠅性誘劑（Me）誘集桔小實蠅。將長約4 cm和直徑1 cm 搓實的脫脂棉誘芯固定在誘捕器螺絲孔的細(xì)鐵絲上， 在誘捕器內(nèi)放置適量脫脂棉誘芯，用細(xì)滴管將混有8%馬拉松的Me誘劑3～4 mL加入棉芯中，蓋上誘捕器蓋子。誘捕器懸掛于獼猴桃枝條下方10 cm及距離地面1 m。

固定專人檢查和記錄監(jiān)測結(jié)果，每7 d收集記錄誘捕器內(nèi)桔小實蠅數(shù)量。依棉芯的濕潤程度，15 d加1次誘劑，每次2 mL，棉芯使用2個月后更新。

1.3 數(shù)據(jù)來源及分析

氣象資料由貴陽市氣象局提供，主要包括日平均氣溫（X1）、日最高氣溫（X2）、日最低氣溫（X3）、日降水量（X4）、日平均相對濕度（X5）等5個氣候因子。對氣象因子先進(jìn)行主成分分析，然后再將主要氣象因子與桔小實蠅成蟲數(shù)量（Y）進(jìn)行相關(guān)分析、聚類分析。數(shù)據(jù)分析采用SPSS 19.0 進(jìn)行。

2 結(jié)果與分析

2.1 貴陽市獼猴桃園桔小實蠅成蟲種群年消長動態(tài)

桔小實蠅成蟲雌雄性比大致為1∶1[10]，因此，雄成蟲誘集量可作為桔小實蠅成蟲種群數(shù)量的估計參數(shù)。

2017、2018年桔小實蠅成蟲誘集數(shù)量年變動情況如圖1 所示。圖中的誘集量為每月誘集量的總和。從圖1 可以看出，貴陽市獼猴桃園桔小實蠅成蟲種群變動在2017、2018年基本一致，年度間無差異（df=9，F(xiàn)=1.165，P=0.405）。1—6月，每月誘集成蟲數(shù)僅為0頭或1頭，表明該時期無桔小實蠅成蟲活動或活動量較少;從7月中旬開始，誘集到的成蟲數(shù)量開始上升，表明其種群開始活躍; 2017年8、9、10月誘集量分別占全年誘集總量13.28%，65.56%，17.29%，3個月誘集量占全年誘集總量的96.13%;2018年8、9誘集量分別占全年誘集總量77.39%，12.17%，2個月誘集量占全年誘集總量的89.57%，表明8月至10月為桔小實蠅成蟲發(fā)生盛期;11月種群數(shù)量開始急劇下降，12月未誘集到成蟲，種群逐漸消退。

2.2 桔小實蠅成蟲種群變動影響氣象因子分析

2.2.1 主成分分析 將2017—2018年各監(jiān)測日的氣象因子（Xi）進(jìn)行主成分分析，結(jié)果見表1。從表中可看出，第1個主成分的特征根為3.037，方差貢獻(xiàn)率為60.747%，代表了全部性狀信息的60.747%，是最重要的主成分;第2個主成分的特征根為1.205，方差貢獻(xiàn)率為24.097%， 累計代表了全部性狀信息的84.844%，其他主成分的貢獻(xiàn)率之和為15.156%。前2個主成分的累積方差貢獻(xiàn)率超過84%，說明前2個主成分把影響貴陽市獼猴桃園桔小實蠅種群數(shù)量變動的氣象因子84%以上的性狀信息反映出來。由此，可應(yīng)用前2個主成分的氣象因子分析其對獼猴桃園桔小實蠅種群變動的影響。

從表2可以看出，在第1主成分中，反映氣溫的指標(biāo)即日平均氣溫（X1）、日最高氣溫（X2）、日最低氣溫（X3）等性狀具有較大的正系數(shù)值，說明第1主成分主要反映了影響獼猴桃園桔小實蠅種群數(shù)量變動的氣象因子中氣溫的特性，第1主成分可看作為氣溫指標(biāo);第2主成分中，日降雨量（X4）、日平均相對濕度（X5）的系數(shù)值較大，說明第2主成分主要反映了影響獼猴桃園桔小實蠅種群數(shù)量變動的氣象因子中降雨量的特性，因此，第2主成分可看作為降水指標(biāo)。

2.2.2 相關(guān)性分析 將主成分分析中代表氣溫、降水的氣象因子變量（Xi），與2017、2018年各監(jiān)測日誘集監(jiān)測的桔小實蠅成蟲數(shù)量（Y）一起進(jìn)行相關(guān)性分析。

結(jié)果表明（表3），桔小實蠅種群數(shù)量的日變化與日平均氣溫（X1）、日最低氣溫（X3）相關(guān)性極顯著，與日最高氣溫（X2）相關(guān)性顯著。

2.2.3 聚類分析 將2017—2018年對種群動態(tài)影響較大的氣象因子即日平均氣溫（X1）、日最高氣溫（X2）、日最低氣溫（X3）、日降雨量（X4）及日平均相對濕度（X5）等5個氣象因子進(jìn)行聚類分析，結(jié)果表明（圖2），5個因子聚成3類，日平均氣溫（X1）、日最低氣溫（X3）及日最高氣溫（X2）先聚在一起，說明三者聯(lián)系最緊密;日降雨量（X4）、日平均相對濕度（X5）和桔小實蠅成蟲誘集量（Y）單獨存在。之后3小組先聚成兩大類群后再聚在一起，日平均氣溫（X1）、日最低氣溫（X3）、日最高氣溫（X2）與日降雨量（X4）聚在一起，日平均相對濕度（X5）單獨存在，因此日平均氣溫（X1）、日最低氣溫（X3）、日最高氣溫（X2）、日降雨量（X4）對種群的影響比日平均相對濕度（X5）更直接，分析結(jié)果和主成分分析及相關(guān)性分析結(jié)果相符。

3 結(jié)論與討論

本研究通過對貴陽市主要獼猴桃基地獼猴桃園桔小實蠅進(jìn)行系統(tǒng)誘集監(jiān)測，結(jié)果表明桔小實蠅成蟲種群變動趨勢在所調(diào)查的兩年中基本一致，1—6月其種群發(fā)生數(shù)量較低，種群數(shù)量從7月開始上升，8月至10月為種群盛發(fā)期，11月種群數(shù)量開始急劇下降，12月種群逐漸消退。

在氣溫、降雨量等因素的綜合影響下，桔小實蠅成蟲種群數(shù)量發(fā)生變動。氣溫與降雨量之間，氣溫是影響貴陽市獼猴桃基地獼猴桃園桔小實蠅成蟲種群變動的最主要氣候因子。貴陽市1—6月，氣溫大部分在桔小實蠅種群發(fā)生最適溫度以下，種群增長較慢，7月份以后，溫度逐漸上升，種群數(shù)量開始上升，但誘集量很少，可能與7月前獼猴桃園桔小實蠅種群基數(shù)較低有關(guān)，也可能與其他氣候因子、寄主食物等因素有一定關(guān)系; 8月至10月桔小實蠅種群維持較高的水平，種群發(fā)生最高峰2017年出現(xiàn)在9月中旬，2018年出現(xiàn)在8月下旬，兩年發(fā)生最高峰期的不一致，可能是由于其他因子及當(dāng)年獼猴桃桃園內(nèi)及周邊種植的寄主植物有一定的關(guān)系，11月至12月氣溫驟降，低溫抑制了種群的增長，桔小實蠅種群活動受限，種群逐漸消退。

降雨量雖然對獼猴桃園桔小實蠅種群影響沒有氣溫影響直接，但降雨量直接影響獼猴桃園土壤濕度、平均相對濕度等要素，間接影響桔小實蠅成蟲交配、產(chǎn)卵、幼蟲化蛹及蛹羽化。其他因子對獼猴桃園桔小實蠅種群的影響，還有待進(jìn)一步研究。

受氣溫、降雨量、寄主植物以及其他因素的綜合影響，每年8月至10月為貴陽市獼猴桃基地獼猴桃園桔小實蠅種群發(fā)生盛期，在8月下旬或9月中旬桔小實蠅種群數(shù)量達(dá)到最高峰。在實際生產(chǎn)中應(yīng)做好種群發(fā)生動態(tài)監(jiān)測，在種群發(fā)生初期及盛期，采取必要的防治措施加以控制。

參考文獻(xiàn)：

[1] 劉奎， 唐良德， 邱海燕， 等. 海南儋州地區(qū)橘小實蠅成蟲數(shù)量動態(tài)[J]. 生物安全學(xué)報， 2016， 25（1）： 23-26.

[2] 雷邦海， 鐘建明， 潘文桃. 凱里地區(qū)果蔬實蠅種類監(jiān)測初報[J]. 西南農(nóng)業(yè)學(xué)報， 2009， 22（4）： 963-965.

[3] 王在凌， 徐婧， 張潤志. 中國重要檢疫性實蠅的全球分布和入侵情況[J]. 生物安全學(xué)報， 2020， 29（3）： 164-169.

[4] 林明光， 張艷， 曾玲， 等. 海南番石榴園桔小實蠅種群動態(tài)與綜合防治研究[J]. 中國南方果樹， 2014， 43（6）： 60-63.

[5] 賀海雄， 梁家燕， 張季， 等. 黔南州都勻市柑桔園實蠅發(fā)生種群及消長規(guī)律[J]. 中國南方果樹， 2016， 45（4）： 35-36.

[6] 朱俐遐， 李桂珍， 堯美英， 等. 攀枝花地區(qū)杧果園桔小實蠅成蟲種群動態(tài)規(guī)律[J]. 中國南方果樹， 2019， 48（2）： 30-33.

[7] 周雙云， 劉增亮， 龍興， 等. 桔小實蠅在蓮霧上的發(fā)生規(guī)律及產(chǎn)卵選擇性研究[J]. 中國南方果樹， 2020， 49（1）： 65-69.

[8] 黃慧欣， 何昌瑋， 李媛， 等. 南寧市火龍果產(chǎn)區(qū)桔小實蠅成蟲消長規(guī)律研究[J]. 中國南方果樹， 2020， 49（2）： 65-70.

[9] 巫輔民， 王小云， 鄭霞林， 等. 廣西柿主產(chǎn)區(qū)桔小實蠅的轉(zhuǎn)主為害規(guī)律研究[J]. 中國南方果樹， 2021， 50（1）： 134-140.

[10] 和萬忠， 孫兵召， 李翠菊， 等. 云南河口縣桔小實蠅生物學(xué)特性及防治[J]. 昆蟲知識， 2002， 39（1）： 50-52.