劉雪鋒
(桂陽縣林業局, 湖南 桂陽 424000)
性信息素引誘監測預報馬尾松毛蟲試驗
劉雪鋒
(桂陽縣林業局, 湖南 桂陽 424000)
采用人工合成的性信息素橡膠誘芯和船形誘捕器對馬尾松毛蟲成蟲進行林間誘捕和幼蟲蟲口密度監測試驗,結果表明:幼蟲低蟲口密度、成蟲低誘捕量時(幼蟲蟲口密度100條/株以下,成蟲誘捕量在65只以下),成蟲誘捕量與幼蟲蟲口密度有較明顯的相關性,可以用模擬的冪函數1n(T)=2.278×1n(p)-2.521進行監測預報;幼蟲高密度、成蟲高誘捕量情況下,成蟲誘捕量與幼蟲蟲口密度相關性不明顯。
馬尾松毛蟲; 成蟲; 幼蟲; 性信息素; 監測預報
馬尾松毛蟲(Dendrolimuspunctatus)屬鱗翅目枯葉蛾科,在淮河以南廣大地區均有分布,是馬尾松林最主要的歷史性害蟲。該害蟲以幼蟲蠶食馬尾松、濕地松、火炬松等松樹針葉,影響松樹生長,危害嚴重時可造成松樹成片死亡。人體接觸其毒毛易引起皮炎、關節腫痛等,影響人體健康[1]。監測預報是馬尾松毛蟲防治中最基礎的工作,以前采用的調查枝條幼蟲、卵塊等監測預報方法落后,其影響因素多,人工勞動強度大,測報準確率低。性信息素引誘監測是一項先進、操作簡便、靈敏度高的監測實用技術[2]。為了提高林間監測水平,積累數據,建立有效適用的馬尾松毛蟲監測預報數學模型,我們在桂陽縣馬尾松毛蟲發生區開展了性信息素引誘監測預報馬尾松毛蟲試驗。
試驗地選在桂陽縣馬尾松毛蟲常災區的流峰鎮、樟木鄉與偶災區的仁義鎮。 該3個鄉鎮位于桂陽縣中北部,以丘陵地貌為主,平均海拔350m,土壤為紅壤,氣候屬中亞熱帶濕潤季風氣候,年平均氣溫15.3~18.2℃,無霜期274d,日照時數1758h,年降雨量1300~1600mm。常災區為馬尾松純林,集中連片,林下植被稀少,有鐵芒萁、白茅等,生物多樣性不豐富,森林生態環境穩定性差,松毛蟲災害暴發周期短(3~4年);偶災區也為馬尾松純林,但伴生油茶、楓香、杜英等其它闊葉樹,林下灌木、草本種類較多,有拔葜、檵木、杜鵑、烏藥、雜竹、白茅、鐵芒萁、野莓等,蓋度大,生物多樣性豐富,森林生態系統較穩定,松毛蟲災害暴發周期長(8~10年)。
2.1試驗材料
誘芯采用北京中捷四方生物科技有限公司生產的、含有人工合成的馬尾松性信息素的新型袖口式復合橡膠塞,其藥效持續期在30d以上。誘捕器為紙塑板做成的雙層船形結構,下板放膠片,誘芯用鐵絲固定于距下板膠片2cm的空中。
2.2試驗方法
2.2.1 誘捕器的掛放地點與懸掛方法 試驗區共設置30個樣地,每1個樣地掛放誘捕器2個,共掛放誘捕器60個,其中常災區的流峰鎮設11個樣地,樟木鄉設13個樣地;偶災區仁義鎮設6個樣地。誘捕器懸掛處要通風、相對空曠,間距應大于100m,高2m,誘捕器保持水平[3-4]。
2.2.2 成蟲誘捕效果觀測 在成蟲羽化前按要求正確懸掛誘捕器,成蟲發生期每3d觀測1次誘捕器,并記錄誘捕成蟲數量和天氣狀況,同時清空誘捕到的成蟲。要注意檢查膠片和誘芯的誘捕效果,并及時更換膠片和誘芯[5]。
2.2.3 幼蟲蟲口密度調查 幼蟲蟲口密度調查在3~4齡期進行,調查時間要相對固定,第1代調查時間6月15—30日,第2代8月25至9月10日,越冬代11月1—15日,樣地幼蟲蟲口密度調查采用樣枝調查法,每1個樣地取樣枝100枝,記錄樣枝蟲口數。第1次幼蟲調查時記錄每株樣木的枝盤數、小枝數、枝條數[6]。
2.3數據處理方法
樣地每木平均蟲口密度計算公式如下。
Nt=Q×S×T×M/W
成蟲誘捕量與幼蟲蟲口密度的相關性分析和模擬曲線圖利用SPSS13.0進行數據處理與分析。
3.1成蟲誘捕量與幼蟲蟲口密度的相關性
成蟲誘捕量與幼蟲蟲口密度的相關性見表1。從表1可以看出,平均成蟲誘捕量從2007年越冬代4.5只較快遞增至2007年第2代41只,達最高峰,然后又較快遞減至2008年第1代27.5只;平均幼蟲蟲口密度相應的從2007年第1代開始緩慢增長,經越冬代快速遞增后,在2008年第1代達到猖獗暴發,蟲口密度最高,達105條/株,至2008年第2代又顯著下降至10條/株。充分說明在一定范圍內成蟲誘捕量與幼蟲蟲口密度是呈正相關的,即成蟲誘捕量多,下代幼蟲蟲口密度大,成蟲誘捕量超過一定數量后,下代幼蟲蟲口密度反而下降。這與我們實際調查時蟲口密度由平衡狀態的低蟲口密度在誘因影響下逐漸上升,經增殖階段蟲口迅速達到極高密度,暴發成災后,因食物銳減和天敵的增加,蟲口密度迅速下降和雄蟲比例增加的自然規律相符。

表1 馬尾松毛蟲成蟲誘捕量與幼蟲蟲口密度Tab1 ThetrappingnumberofadultandthepopulationdensityoflarvaofD.punctatus樣地號誘捕器號誘捕量(只)蟲口密度(條/株)2007年2008年2007年2008年越冬代第1代第2代越冬代第1代第1代第2代越冬代第1代第2代11~2 2 7 642926 510220 16 623~4 1712632154537662135~6 0145080331532179287047~8 4812640317916297059~1019102514056283956611~120034292524151211713~140362363511161302815~162043383001011710917~1841424252210129811019~20112127201002531121~222271446047101278166261223~2411261223121098451325~262036272354929861427~2841342042510176371529~30393814271202162311631~322316253723234141733~34532632320213430

續表1 馬尾松毛蟲成蟲誘捕量與幼蟲蟲口密度ContinuedTab1 ThetrappingnumberofadultandthepopulationdensityoflarvaofD.punctatus樣地號誘捕器號誘捕量(只)蟲口密度(條/株)2007年2008年2007年2008年越冬代第1代第2代越冬代第1代第1代第2代越冬代第1代第2代1835~36456333241602831937~38283230352203932039~4011045393210118352141~423162332391102752243~44011424230215722345~463658312412122922447~483313302802011302549~5016520181257231012651~521041715214012012753~5414431201322131432855~561361526162700422957~58169312515764943059~60129152219104110平均456441030727589533191053103
3.2成蟲誘捕量與幼蟲蟲口密度相關性模擬曲線圖估計
將2007—2008年馬尾松毛蟲各代成蟲誘捕量和下1代幼蟲蟲口密度疊加,利用SPSS15.0,選用線性估計和二次方估計2種方法進行曲線圖模擬估計,結果見圖1。
從圖1中可看出,成蟲誘捕量在65只以下,幼蟲蟲口密度100條/株以下,實測數與估計的曲線離散度低,成蟲誘捕量與幼蟲蟲口密度相關性較明顯;而成蟲誘捕量超過70只,幼蟲蟲口密度在100條/株以上時,實測數偏離估計的曲線較遠,成蟲誘捕量與幼蟲蟲口密度相關性不明顯。
3.3幼蟲不同蟲口密度下成蟲誘捕量與幼蟲蟲口密度的相關性模擬曲線估計

圖1 成蟲誘捕量與幼蟲蟲口密度模擬曲線圖Fig.1 The simulation curve of adult trapping number and larval population density
將幼蟲低蟲口密度分為0~5、5~10、10~15、15~20、20~25、25~30、30~40、40~50條/株等8個梯度,高蟲口密度分為100~150、150~200、200~250、250~300、300~350條/株等5個梯度,對樣地的成蟲誘捕量取算數平均值后,利用SPSS15.0以線性估計、二次方估計、冪函數估計和直屬估計等4種方法進行曲線估計,并計算幼蟲低蟲口密度下成蟲誘捕量與幼蟲蟲口密度的相關性曲線估計參數。模擬和計算結果見圖2、圖3和表2、表3。
由圖2、圖3和表2、表3可以看出,在幼蟲低蟲口密度時,模擬曲線與處理數據點緊密,方差小,成蟲誘捕量和幼蟲蟲口密度之間存在著較明顯的相關性;

圖2 幼蟲低蟲口密度下成蟲誘捕量與幼蟲 蟲口密度的相關性曲線模擬圖Fig.2 The correlation simulation curve of adult trapping number and larval population density with low po- pulation density of larva

表2 幼蟲低蟲口密度下成蟲誘捕量與幼蟲蟲口密度的相關性參數估計表Tab2 Thecorrelationparameterestimationofadulttrappingnumberandlarvalpopulationdensitywithlowpopulationdensityoflarva參數線性估計二次方估計冪函數估計指數函數參數估計R0864087709080854R20746076908240729經調整的R20710069207990690方差2425250102170269

圖3 幼蟲高蟲口密度下成蟲誘捕量與幼蟲蟲口密度 的相關性曲線模擬圖Fig.3 The correlation simulation curve of adult trapping number and larval population density with high population density of larva

表3 幼蟲高蟲口密度下成蟲誘捕量與幼蟲蟲口密度的相關性參數估計表Tab3 Thecorrelationparameterestimationofadulttrappingnumberandlarvalpopulationdensitywithhighpopulationdensityoflarva參數線性估計二次方估計冪函數估計指數函數參數估計R0688069507150699R20473048305110488經調整的R20297003403480317方差503796110405980612
在幼蟲高蟲口密度時,模擬曲線與處理數據點離散度大,方差也大,成蟲誘捕量和幼蟲蟲口密度的相關性不顯著。幼蟲低蟲口密度時,利用冪函數曲線估計的方差最小,成蟲誘捕量和幼蟲蟲口密度的相關性最高。幼蟲低蟲口密度時冪函數曲線估計公式如下。
ln(T)=lnb0+b1ln(p)
經計算,模擬的函數方程為:
ln(T)=2.278×ln(p)-2.521
(1)馬尾松毛蟲成蟲性信息素誘捕量與幼蟲蟲口密度存在著制約關系,即當幼蟲蟲口密度很大時,成蟲誘捕量反而下降;成蟲誘捕量極大時,幼蟲蟲口密度又迅速下降。在幼蟲低蟲口密度、成蟲低誘捕量(幼蟲蟲口密度100條/株以下,成蟲誘捕量在65只以下)時,成蟲誘捕量和幼蟲蟲口密度之間呈冪函數相關關系,我們可以利用幼蟲低蟲口密度時建立的數學模型ln(T)=2.278×ln(p)-2.521對馬尾松毛蟲的種群動態進行較準確監測預測。在幼蟲高蟲口密度、成蟲高誘捕量時,成蟲誘捕量和幼蟲蟲口密度之間的相關關系不明顯,當馬尾松毛蟲大暴發時,利用誘捕器監測的數據不準確。
(2) 成蟲誘捕量與幼蟲蟲口密度相關性的準確度,必須要充分保證誘芯質量,性信息素的揮發要持續均勻穩定,并要及時排除人為因素干擾。本試驗2007年第2代2號樣地成蟲誘捕量極少,人為過度修枝和取柴是一個重要原因。同時,還要綜合考慮針葉保存率、防治率、發生類型等影響因素。
[1] 蕭剛柔.中國森林昆蟲(第2版)[M].北京:中國林業出版社,1992:76-77.
[2] 陳昌杰. 松毛蟲綜合管理[M]. 北京:中國林業出版社,1990:5-18.
[3] 侯陶謙.中國主要害蟲綜合防治[M]. 北京:科學出版社, 1979:370-400.
[4] 張愛兵,譚聲江,高偉,等.應用性外激素監測馬尾松毛蟲發生數量的初步研究[J]. 昆蟲知識,2001,38(3):211-214.
[5] 蔣杰賢,嚴巍.城市綠地有害生物預警及控制[M].上海:上??茖W技術出版社,2006 :5-18.
[6] 徐明慧. 園林植物病蟲害防治[M]. 北京:中國林業出版社,1993:347-349.
MonitoringandforecastingofDendrolimuspunctatusbysexpheromonelure
LIU Xuefeng
(Forestry Bureau of Guiyang County, Guiyang 424000, China)
The test on trapping adultDendrolimuspunctatusby synthetic sex pheromone lure and monitoring of larval population density in forest was conducted. The result showed that, with low population density of larva (below 100 per plant) and low trapping number of adult (below 65), the trapping number of adult had significantly correlation with the population density of larva, which could be simulated by power function equation 1n(T)=2.278×1n(p)-2.521.With high population density of larva and high trapping number of adult, the trapping number of adult had no obvious correlation with the population density of larva.
Dendrolimuspunctatus; adult; larva; sex pheromone; monitoring and forecasting
2010 — 12 — 01
2011 — 03 — 24
S 763.712.48
A
1003 — 5710(2011)02 — 0017 — 04
10. 3969/j. issn. 1003 — 5710. 2011. 02. 005
(責任編輯:唐效蓉)