槐豆木虱若虫空间分布型及抽样技术研究

槐豆木虱若虫空间分布型及抽样技术研究
吴建霞;沈平;吴建华;朱惠英;王新东;常承秀
【摘要】应用6种分布型指数法分析判定了槐豆木虱[ Cyamophila willieti (Wu)]若虫在国槐(Sophora japonica Linnaeus)上的空间分布型,利用Taylor(1961)幂法则和Iwao(1971)回归方程分析聚集原因,结果表明,槐豆木虱若虫在国槐上呈聚集分布,公共kc值为5.159 0,且符合负二项分布；其种群聚集原因是由昆虫本身行为习性所致.在此基础上,采用Iwao的方法确定了槐豆木虱若虫的田间理论抽样数和序贯抽样表.%The spatial distribution patterns of Cyamophila willieti nymphs in Sophora japonica garden were studied using six spatial distribution pattern analysis methods, and the assembling reasons were analyzed with Taylor power law and Iwao regression equation. The results indicated that the aggregated distribution of the C. Willieti nymphs in field fitted the negative binomial series (kc =5. 159 0); the reason for aggregation might be due to the insect natural habit. Based on the distribution pattern, the theoretical sampling number and sequential sampling table were calculated by Iwao's method.
【期刊名称】《植物保护》
【年(卷),期】2012(038)004
【总页数】4页(P47-49,67)
【关键词】槐豆木虱;空间分布型;抽样技术
【作者】吴建霞;沈平;吴建华;朱惠英;王新东;常承秀
【作者单位】甘肃省临夏市城郊镇林业工作站,临夏731100;甘肃省临夏回族自治
州森林病虫害防治检疫站,临夏731100;甘肃省临夏回族自治州农业学校,临夏731100;甘肃省临夏回族自治州森林病虫害防治检疫站,临夏731100;甘肃省临夏
回族自治州森林病虫害防治检疫站,临夏731100;甘肃省临夏回族自治州森林病虫
害防治检疫站,临夏731100
【正文语种】中文
【中图分类】Q968.1
槐豆木虱［Cyamophila willieti（Wu）]属半翅目（Hemiptera），木虱科（Psyllidae），豆木虱属（Cyamophila）［1]，是国槐（Sophora japonica Linnaeus）上的主要害虫之一，其次为害龙爪槐（Sophora japonicaf.pendula Hort.）、金枝国槐（Sophora japonica‘Golden Stem’）。

目前，国内对槐豆木虱形态学特征、生物学特性、发生与环境因子的关系、防治、天敌等，已有诸多的研究报道［2－12]，但就槐豆木虱若虫空间分布型及抽样技术研究迄今尚未见
报道。

本文报道槐豆木虱若虫在国槐上的空间分布结构，在此基础上，确定出了合理而简便的抽样方法。

1 材料与方法
1.1 调查地点概况
调查样地位于甘肃省临夏市北滨河东路，海拔1 860m，2011年年均气温7.8℃，年日照时数2 171h，无霜期144d以上，年降雨量378mm；该路段国槐单排栽植，东西走向，树龄10年，树高4～5m，株距7m，路面铺设彩砖，树盘面积
0.67m2，调查前未进行药剂防治。

1.2 取样方法
2010年4月下旬－5月上旬，在第1代若虫期，采用线路调查法，在道路北侧每间隔50m设1样段，共设12组样段；样段内每隔1株取1样株，每样段调查10株；每样株从东、西、南、北4个方位，随机选1年生的样枝（10cm），分别统计不同方位样枝上的若虫数量。

1.3 分布型测定
根据调查数据，计算均值m、方差S2、平均拥挤度。

依据聚块性指数／m、扩散系数c、负二项分布k值、扩散型指数Iδ、丛生指数I、聚集度指数Ca参数标准判定空间分布型［13－14]。

1.4 聚集原因检验
利用Taylor（1961）幂法则S2＝amb，即lgS2＝lga＋blgm［14]和Iwao （1971）的－m 回归模型＝α＋βm［15]分析测定种群的分布格局和内部结构；应用Blackith（1961）种群聚集均数（λ＝mγ／2k）检验聚集原因［13]。

1.5 理论抽样数的确定
依据已知Iwao（1971）的－m线性回归方程，应用Iwao的理论抽样数模型：N ＝（t／D）2［（α＋1）／m＋β－1]［14]来确定理论抽样数。

1.6 序贯抽样的确定
根据Iwao的方法［14]，设临界防治指标为m0，则：最大抽样数
1.7 数据统计分析
采用Excel 2003进行数据统计，参数均采用最小二乘法估计，并利用相关系数r 进行相关性检验。

2 结果与分析
2.1 分布型检验
依据槐豆木虱若虫11组样方的均值m、方差S2和平均拥挤度，计算出各样方的聚块性指数／m、扩散系数C、扩散型指数Iδ均大于1，丛生指数I、聚集度指数
Ca大于0（表1），根据各参数空间分布型判定标准［13－14]，得出槐豆木虱若虫在田间呈聚集分布。

表1 槐豆木虱若虫的聚集度指标样段号均值m 方差S2 平均拥挤度聚块性指数m ＊／m扩散性系数C负二项分布K值k扩散型指数Iδ丛生指数I聚集度指数Ca 聚集均数λγ 1 20.175 74.194 4 22.852 5 1.132 7 3.677 5 7.534 91.134 1
2.677 5 0.132 7 19.290 0 14.408 7 2 21.125 82.309 4 24.021 3 1.137 1
3.896
3 7.293 8 1.138 5 2.896 3 0.137 1 20.167 8 13.926 6 3 16.275 51.049
4 18.411 7 1.131 3 3.136 7 7.617 0 1.133 0 2.136 7 0.131 3 15.568 6 14.572 7 4 10.87
5 48.459 4 14.331 0 1.317 8 4.45
6 0 3.146
7 1.320
8 3.456 0
0.317 8 9.747 3 5.640 8 5 7.675 26.969 4 10.188 9 1.327 5 3.513 9 3.053 0
1.331 9
2.513 9 0.327 5 6.855 2 5.453 8 6 1
3.175 68.794 4 17.396 6 1.320 4 5.221 6 3.120 9 1.322 9
4.221 6 0.320 4 11.797 7
5.589 2 7 28.025 152.174 4 32.455 0 1.158 1 5.430 0
6.326 3 1.159 1 4.430 0 0.158 1 26.563 1
11.992 5 8 31.750 162.437 5 35.866 1 1.129 6 5.116 1 7.713 5 1.130 5
4.116 1 0.129 6 30.388 7 14.765 6 9 28.975 143.024 4 32.911 1 1.135 8
4.936 1 7.361 3 1.136 8 3.936 1 0.135 8 27.674 1 14.061 6 10 9.350 44.727
5 13.133 7 1.404 7 4.783 7 2.471 1 1.408 4 3.783 7 0.404 7 8.122 8 4.293
6 11 21.125 91.459 4 24.454 4 1.15
7 6 4.329 4 6.3449 1.159 0 3.329 4
0.157 6 20.026 3 12.029 8
Waters（1959）认为，k值与虫口密度无关，k值愈小，种群的聚集度愈大，若k值趋于∞时（一般在8以上时），则逼近Poisson分布［17]。

表1中的k值均依据矩法求得，且2＜k＜8，因此槐豆木虱若虫的分布型不符合Poisson分布。

利用Bliss和Owen（1958）提出的用矩法来估公共kc 值［18]，计算得公共kc 值为5.159 0，kc 值的χ2 值为16.202 1，且df为10，经χ2（χ2＝18.387）检
验，即χ20.05＝18.387＞χ2＝16.202 1，其差异不显著（p＞0.05），因此认为
槐豆木虱若虫分布符合负二项分布。

2.2 聚集原因分析
利用 Taylor（1961）的方法［15]，以lgS2 和lgm进行线性回归（相关系数R＝0.958 6，当df＝11－2＝9，查相关系数检验表得α0.05＝0.602 1，即｜R｜＝0.958 6＞α0.05＝0.602 1，相关系数R在α＝0.05水平上显著）得方程：
lgS2 ＝0.479 3＋1.128 5lgm式中：lga＝0.479 3＞0，b＝1.128 5＞1，S2＝3.015m1.128 5，说明种群在一切密度下均是聚集分布，且具密度依赖性，聚集强度随种群密度的升高而增强。

利用Iwao（1971）的－m 线性回归模型［16]，计算出槐豆木虱若虫在国槐上的线性回归方程为＝2.680 1＋1.038 4m。

11组抽样数据的统计结果（相关系数R
＝0.616 0，当df＝11－2＝9，查相关系数检验表得α0.05＝0.602 1，即｜R｜
＝0.616 0＞α0.05＝0.602 1，相关系数R在α＝0.05水平上显著）为α＝2.680 1＞0，说明槐豆木虱若虫个体间相互吸引，存在个体群；β＝1.038 4＞1，表明基本成分呈聚集分布。

应用Blackith（1961）种群聚集均数（λ）进行检验［13]，公式为：
λ＝mγ／2k
式中：k为负二项分布的指数，γ是自由度等于2k，以0.5概率值对应处的χ2值，应用内插法计算得出。

根据表1，当λ＜2时，昆虫种群的聚集可能是由某些环境
作用引起的，当λ≥2时，其聚集是昆虫本身行为和环境因素中的任一因子引起的。

经计算，11组样方的总λ值（λ＝17.836 5＞2），说明槐豆木虱若虫聚集是由昆虫本身行为习性所造成的。

2.3 抽样技术
2.3.1 理论抽样数技术
在已知Iwao（1971）的－m 线性回归方程＝2.680 1＋1.038 4 m 的前提下，应用Iwao的理论抽样数模型［14]，确定出理论抽样数公式：
N ＝（t／D）2（3.680 1／m＋0.038 4）取t＝1.96（置信概率为0.95），相对误差（D＝d′／m）取0.1、0.2和0.3，平均若虫数（m）在5、10、15、20……80的条件下，得出槐豆木虱若虫的理论抽样数（N）（表2）。

表2 槐豆木虱若虫理论抽样数相对误差（D）不同若虫密度下的抽样量／株5 10 15 20 25 30 40 50 60 70 80 0.1 252 133 94 74 62 54 44 38 35 32 30 0.2 74 39 27 21 18 15 13 11 10 9 8 0.3 33 17 12 9 8 7 6 5 4 4 4
由表2可知，在相对误差（D）固定时，槐豆木虱若虫的理论抽样数（N）随平均若虫数（m）的增大而下降，只是下降幅度逐渐缩小，当平均若虫数（m）上升到一定数量时，抽样数（N）也趋于稳定；在相同的平均若虫数（m）下，随着相对误差（D）的增大，所需抽样数（N）依次减少。

2.3.2 序贯抽样技术
根据Iwao的方法［14]，取t＝1.96（置信概率为0.95），已知α＝2.680 1，β
＝1.038 4，设初步拟定防治指标m0为10头／株，代入序贯抽样模型，得下式：
用上述两式，计算出槐豆木虱若虫序贯抽样表（表3）供防治参考。

表3 槐豆木虱若虫序贯抽样表
最大抽样数：nmax＝（t／d）2（3.680 1 m0＋0.038 4 m20）
d为m＝m0时估计密度所允许的置信限，当d＝1时，最大抽样数为156株；d
＝2时，最大抽样数为39株。

3 讨论
本文就槐豆木虱若虫的11组样段调查结果进行了空间分布型检验分析，结果表明，槐豆木虱若虫在田间呈聚集分布；其公共kc值为5.159 0，符合负二项分布；聚
集原因是由昆虫本身行为习性所造成的；这与槐豆木虱若虫喜聚集在嫩芽、嫩叶、嫩梢及始花期花序为害习性相符。

在虫情调查时，以样枝为抽样单位，采用随机抽样，根据精确度与工作量，在试验研究时，相对误差要求严格取D＝0.1，在大田调查时相对误差可放宽到D＝0.2
或0.3。

用序贯抽样表进行抽样时，当调查的株数，累计若虫数超过上限时，就必须进行防治，若累计若虫数小于下限时，则不需要进行防治；累计若虫数在上限和下限之间时，应继续增加抽样株数。

当不易做出判断时，需确定最大抽样数目。

因此，在保证一定的虫情调查准确度情况下，只要合理的利用抽样技术进行调查就可取得较准确的虫情资料。

参考文献
［1]李法圣.中国木虱志（昆虫纲：半翅目）［M].北京：科学出版社，2011：695－699.
［2]王维翊，王维中.槐豆木虱研究初报［J].辽宁林业科技，1996，22（2）：38
－39.
［3]王本辉，饶晓明，沈彦刚.槐木虱的发生规律与防治措施［J].甘肃农业科技，2001，38（10）：34－35.
［4]杨友兰，王红武，吕小虎，等.槐豆木虱生物学特性及防治［J].昆虫知识，2002，39（6）：433－436.
［5]贺伟生，菅广林，赵丽，等.龙爪槐木虱生物学特性的初步观察［J].内蒙古林
业科技，2005，31（3）：32－33.
［6]强中兰，马长林.国槐木虱的发生与防治实验［J].甘肃林业科技，2007，32（1）：60－61.
［7]沈平，常承秀，张永强，等.槐豆木虱形态特征及发生规律［J].甘肃林业科技，2008，33（1）：30－33.
［8]吴建霞，朱惠英，沈平，等.临夏市国槐行道树栽植现状、主要虫害及防治对
策［J].甘肃林业科技，2008，33（2）：35－36.
［9]沈平，常承秀，张永强，等.槐豆木虱发生与主要环境因子关系的调查研究［J].甘肃林业科技，2008，33（3）：40－42.
［10]常承秀，沈平，张永强，等.八种药剂防治槐豆木虱药效试验［J].昆虫知识，2009，46（3）：472－474.
［11]沈平，常承秀，朱惠英，等.异色瓢虫对槐豆木虱的捕食作用［J].植物保护，2009，35（1）：66－69.
［12]张永强，沈平，常承秀，等.临夏地区异色瓢虫生物学特性研究［J].甘肃林业科技，2010，35（1）：71－73.
［13]丁岩钦.昆虫种群数学生态学原理与应用［M].北京：科学出版社，1980：84－124.
［14]徐汝梅.昆虫种群生态学［M].北京：北京师范大学出版社，1988：8－60. ［15]Taylor L R.Aggregation，variance and the mean［J].Nature，1961，189：732－735.
［16]Iwao S，Kuno E.An approach to the analysis of aggregation pattern in biological populations［C]∥ Patil G P，Pielou E C，Waters W E.Statistical Ecology，Pennsylvania：The Pennsylvania State University Press，1971，1：461－513.
［17]Waters W E.A quantitative measure of aggregation in insects
［J].Journal of Economic Entomology，1959，52：1180－1184.
［18]Bliss C I，Owen A R G.Negative binomial distributions with a common K［J].Biometrika，1958，45：37－58.。