遼寧省榛黃達(dá)癭蚊種群動(dòng)態(tài)和防治適期研究

李維根，姜 策，陳 靜，馬海軍，孟 威，高 萍，王洪平

（1遼寧省綠色農(nóng)業(yè)技術(shù)中心，沈陽110034）（2沈陽農(nóng)業(yè)大學(xué)植物保護(hù)學(xué)院）

榛是四大堅(jiān)果之一，因其獨(dú)特的風(fēng)味和豐富的營養(yǎng)深受廣大消費(fèi)者喜愛。榛種類豐富，有平榛、歐榛以及平歐雜交榛，其中平榛主要栽培在我國東北地區(qū)。遼寧省作為榛的主要產(chǎn)區(qū)之一，以平榛居多，栽培面積近13.3萬hm2，主要分布在遼寧省東部和北部山區(qū)，平均667 m2產(chǎn)量50～60 kg，按目前市場售價(jià)，正常情況下平均667 m2收益約1 000元[1]。多數(shù)平榛處于半野生狀態(tài)，病蟲害發(fā)生較重，榛黃達(dá)癭蚊（Dasinura corylifalva）為其常發(fā)害蟲之一。該蟲通過幼蟲取食葉片，形成蟲癭，爭奪葉片養(yǎng)分，影響葉片功能，減少光合產(chǎn)物，導(dǎo)致產(chǎn)品質(zhì)量下降；也可為害果實(shí)，引起榛果品質(zhì)低下。危害嚴(yán)重時(shí)，被害株率可高達(dá)98%，榛果被害率26%～87%[2]。

目前，生產(chǎn)中防治榛黃達(dá)癭蚊仍以化學(xué)防治為主，且尚無明確的防治適期和防治指標(biāo)，造成盲目用藥和濫用藥，不利于榛產(chǎn)業(yè)的健康發(fā)展[2-5]。為此，我們于2016—2018年進(jìn)行了相關(guān)的試驗(yàn)，現(xiàn)將試驗(yàn)結(jié)果總結(jié)如下。

1 材料與方法

1.1 試驗(yàn)地點(diǎn)

試驗(yàn)地設(shè)在遼寧省撫順市東洲區(qū)章黨鎮(zhèn)上二村上二北溝，地理位置為東經(jīng)123°57′29.43′′，北緯41°59′21.12′′，海拔243 m。試驗(yàn)地年平均氣溫為6.6 ℃，≥10 ℃年活動(dòng)積溫2 700～3 200 ℃，無霜期為130～150 d，平均年降水量為760～790 mm。試驗(yàn)調(diào)查區(qū)總面積約21 000 m2。

1.2 調(diào)查方法

2016—2018年每年從5月10日左右開始調(diào)查，到7月10日左右結(jié)束，每隔7 d左右到田間調(diào)查榛黃達(dá)癭蚊的發(fā)生情況。采用對(duì)角線法5點(diǎn)取樣，每點(diǎn)隨機(jī)調(diào)查300個(gè)枝條，隨機(jī)分成4個(gè)大約相等的樣方，總計(jì)每次調(diào)查1 500個(gè)枝條，采取有放回抽樣，下次抽樣仍有抽到的概率。調(diào)查時(shí)，按樣方、樣點(diǎn)記錄調(diào)查到的有效蟲癭數(shù)量。有效蟲癭是指新鮮有生長趨勢(shì)的蟲癭，干癟、枯死的蟲癭不計(jì)算。

1.3 數(shù)據(jù)處理與分析方法

因每次調(diào)查均為隨機(jī)區(qū)組方式調(diào)查，調(diào)查到的蟲癭都保留在枝條上，下次調(diào)查的就是累計(jì)結(jié)果，因此，需要減去上次調(diào)查到的枝條蟲癭數(shù)，才能體現(xiàn)蟲癭的增量。

1.4 種群聚集度分析方法

Kunou CA指標(biāo)[6]：CA=（s2-m）/m2；叢生指標(biāo)[7]：I=s2/m-1；聚塊性指標(biāo)[8]：m*/m，m*=m+s2/m-1；平均擁擠度[9-10]：m*=α+βm。式中，s2為樣本平方和，m*為平均擁擠度，m為樣本平均值，α 為截距，β 為斜率。

2 結(jié)果與分析

2.1 調(diào)查結(jié)果

2016年調(diào)查結(jié)果表明，榛黃達(dá)癭蚊蟲癭增量最多的日期是6月19日，較上次調(diào)查（6月10日）增加154個(gè)蟲癭，此后蟲癭數(shù)量仍然很多，但增量是下降的（表1）。

表1 2016年榛黃達(dá)癭蚊蟲癭數(shù)動(dòng)態(tài) 個(gè)

2017年調(diào)查結(jié)果表明，榛黃達(dá)癭蚊蟲癭增量最多的日期是6月20日，比上次調(diào)查時(shí)間（6月9日）增加84個(gè)蟲癭，此后數(shù)量開始下降（表2）。

2018年調(diào)查結(jié)果表明，榛黃達(dá)癭蚊蟲癭增量最多的日期是6月12日，較上次調(diào)查時(shí)間（6月3日）增加109個(gè)蟲癭（表3）。

表2 2017年榛黃達(dá)癭蚊蟲癭數(shù)動(dòng)態(tài) 個(gè)

表3 2018年榛黃達(dá)癭蚊蟲癭數(shù)動(dòng)態(tài) 個(gè)

2.2 種群動(dòng)態(tài)數(shù)據(jù)分析

由圖1可知，榛黃達(dá)癭蚊蟲癭增量變化情況為：2016年6月19日達(dá)到高峰，然后蟲癭增量下降，6月26日增量為負(fù)值，說明6月26日后沒有新的蟲癭產(chǎn)生，而原有的蟲癭逐漸干癟、枯死。2017年6月20日達(dá)到高峰，6月末以后不再增加新蟲癭，原有蟲癭逐漸干癟、枯死。2018年6月12日有1個(gè)顯著的高峰，6月下旬蟲癭增量不再增加。連續(xù)3年的試驗(yàn)結(jié)果表明，榛黃達(dá)癭蚊1年完成1代，每年危害造成的蟲癭數(shù)量變化呈單峰型曲線，2018年稍早些，2016年和2017年相近。這些與以往研究結(jié)果時(shí)間表述方面有所差異[11]。

2.3 防治適期的推測

由圖2可知，榛黃達(dá)癭蚊蟲癭數(shù)量變化情況為：2016年6月19日達(dá)到高峰，2017年6月27日達(dá)到高峰，2018年6月25日達(dá)到高峰。一般按蟲癭率達(dá)到10%進(jìn)行防治。從連續(xù)3年得到的數(shù)據(jù)分析，2016年達(dá)到10%蟲癭率的時(shí)間約為6月10日，2017年約為6月9日，2018年約為6月7日，說明年度間差異較小，基本可以按照特定時(shí)間進(jìn)行藥劑防治?？紤]到卵的發(fā)育歷期需要5～7 d，因此防治成蟲產(chǎn)卵盛期需要再提前7 d左右。

圖1 2016—2018年榛黃達(dá)癭蚊蟲癭增量變化情況

圖2 2016—2018年榛黃達(dá)癭蚊蟲癭累積量動(dòng)態(tài)

2.4 種群空間分布型的時(shí)間動(dòng)態(tài)分析

對(duì)2016—2018年調(diào)查結(jié)果按樣方數(shù)據(jù)進(jìn)行分析，計(jì)算每次調(diào)查各樣方的平均擁擠度和平均密度值，分析二者的線性關(guān)系，計(jì)算出Iwao公式中的α值和β 值，計(jì)算出決定系數(shù)R2，然后計(jì)算出每次調(diào)查的平均擁擠度以及平均種群密度，計(jì)算出聚塊性指標(biāo)，然后計(jì)算每次調(diào)查的叢生指標(biāo)和CA指標(biāo)。

2016年5月13日、20日調(diào)查，蟲癭數(shù)量極少，所以無法獲得以樣方為單位的平均擁擠度截距值和斜率值（α 值和β 值），也無法計(jì)算決定系數(shù)R2，后面的結(jié)果表明，種群密度較低時(shí)，個(gè)體間排斥，但呈聚集分布，單個(gè)個(gè)體存在，而種群密度較高時(shí)，則以個(gè)體群為存在單位。叢生指標(biāo)、CA指標(biāo)和聚塊性指標(biāo)結(jié)果表明，5月27日開始調(diào)查的數(shù)量結(jié)果，都是聚集性分布（表 4）。

2017年調(diào)查結(jié)果與2016年相似。根據(jù)平均擁擠度的分析，5月26日以后調(diào)查，密度較高時(shí)，仍然是個(gè)體群為主，密度較低時(shí)，以單個(gè)個(gè)體為主。但5月26日以前，各項(xiàng)聚集度指標(biāo)都表明，均勻分布和隨機(jī)分布現(xiàn)象都可出現(xiàn)，這是由于種群密度低，調(diào)查誤差較大導(dǎo)致。而后面數(shù)次結(jié)果都是聚集性分布（表5）。

表4 2016年榛黃達(dá)癭蚊蟲癭數(shù)聚集度指標(biāo)

表5 2017年榛黃達(dá)癭蚊蟲癭數(shù)聚集度指標(biāo)

2018年試驗(yàn)結(jié)果與前2年差異不大。平均擁擠度參數(shù)表明，5月28日表現(xiàn)為個(gè)體間排斥，表現(xiàn)為聚集分布，再以后則都是個(gè)體間相互吸引聚集，表現(xiàn)為不均勻的聚集分布。聚塊性指標(biāo)則表明早期為均勻分布，后來則為聚集分布。叢生指標(biāo)與聚塊性指標(biāo)結(jié)果相同（表6）。

表6 2018年榛黃達(dá)癭蚊蟲癭數(shù)聚集度指標(biāo)

2.5 垂直分布特性分析

垂直分布特點(diǎn)是：蟲癭基本都在上部葉片或果上，中部及以下的部位基本沒有蟲癭。蟲癭在葉片上的分布極其不均勻。有時(shí)一片葉片上多達(dá)14個(gè)蟲癭，而很多葉片上卻沒有蟲癭，說明其分布聚集度變化很大。

3 結(jié)論

本試驗(yàn)結(jié)果表明，隨著時(shí)間的推移，榛黃達(dá)癭蚊蟲癭數(shù)量逐漸增加，至高峰后下降。種群數(shù)量變化的前期，空間分布型呈稀疏的泊松分布，后來都趨近于不均勻的核心分布或嵌紋分布。由于5月蟲癭數(shù)量較少，調(diào)查誤差會(huì)較大，導(dǎo)致計(jì)算結(jié)果誤差較大，甚至出現(xiàn)結(jié)果與事實(shí)不相符的情況。這就需要增加調(diào)查數(shù)量，降低相應(yīng)的誤差。

昆蟲種群空間分布型的聚集性指標(biāo)還有幾種，如Taylor冪法則等也是常用方法[12-13]，雖然本文沒有采用，但用這些方法分析出來的結(jié)果應(yīng)該與本文結(jié)論相似。本文通過定期抽樣調(diào)查幼蟲造成蟲癭的數(shù)量動(dòng)態(tài)，掌握幼蟲危害的數(shù)量動(dòng)態(tài)。這種方法對(duì)抽樣時(shí)間要求不嚴(yán)格，獲得的蟲癭數(shù)量和變化動(dòng)態(tài)數(shù)據(jù)都十分穩(wěn)定，便于觀察和計(jì)數(shù)。根據(jù)蟲癭數(shù)量的動(dòng)態(tài)變化，結(jié)合文獻(xiàn)記載榛黃達(dá)癭蚊在遼寧地區(qū)為1年發(fā)生1代等生物學(xué)特性，就可以確定合適的防治時(shí)期。最后獲得統(tǒng)計(jì)學(xué)意義上的幼蟲數(shù)量動(dòng)態(tài)，再根據(jù)卵的發(fā)育進(jìn)度，推算產(chǎn)卵前期的末期，作為成蟲防治的最后時(shí)期，也是最佳防治適期。通過3年試驗(yàn)，發(fā)現(xiàn)不同年份的發(fā)生量和發(fā)生時(shí)間以及增長數(shù)量都存在差別，但最佳防治時(shí)期卻相近，最有利于化學(xué)防治。