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

（1遼寧省綠色農業(yè)技術中心，沈陽110034）（2沈陽農業(yè)大學植物保護學院）

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

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

1 材料與方法

1.1 試驗地點

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

1.2 調查方法

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

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

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

1.4 種群聚集度分析方法

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

2 結果與分析

2.1 調查結果

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

表1 2016年榛黃達癭蚊蟲癭數(shù)動態(tài) 個

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

2018年調查結果表明，榛黃達癭蚊蟲癭增量最多的日期是6月12日，較上次調查時間（6月3日）增加109個蟲癭（表3）。

表2 2017年榛黃達癭蚊蟲癭數(shù)動態(tài) 個

表3 2018年榛黃達癭蚊蟲癭數(shù)動態(tài) 個

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

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

2.3 防治適期的推測

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

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

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

2.4 種群空間分布型的時間動態(tài)分析

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

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

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

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

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

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

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

2.5 垂直分布特性分析

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

3 結論

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

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