蓟马在大葱田间分布型及其抽样技术

李平

摘要：采用空间分布型检验、聚集强度指标检验和线性回归方法调查研究了甘肃省武威市凉州区大葱田间蓟马空间分布型及其抽样技术。结果表明，在甘肃省金昌武威地区，大葱田间蓟马空间分布型呈聚集分布，蓟马聚集受环境影响较大。建立了若虫与成虫的理论抽样模型。

关键词：蓟马;大葱田;空间分布型;理论抽样模型

中图分类号：S436.33      文献标志码：A      文章编号：1001-1463（2021）10-0018-05

doi：10.3969/j.issn.1001-1463.2021.10.005

Spatial Distribution Pattern and Sampling Technology of Thrips in the Green Chinese Onion Field

LI Ping

（Wuwei Agricultural and Technology Extension Center， Wuwei Gansu 733000， China）

Abstract：The spatial distribution pattern and sampling technology of thrips in the green Chinese onion field in the Jinwu area of Gansu Province were investigated by using the methods of spatial distribution pattern test， aggregation intensity index test and linear regression. The results showed that the spatial distribution of thrips was aggregated. The aggregation was influenced by the environment. The theoretical sampling model of the immature and adults thrips was established.

Key words：Thrips; Green Chinese onion field; Spatial distribution pattern; Theoretical sampling model

甘肃河西灌区是大葱育苗及主要栽培地区之一，加强大葱田间病虫预测预报是保障大葱产业发展的科学措施之一[1 - 4 ]。近年来，在甘肃金昌武威地区大葱田间调查发现，蓟马是为害大葱最主要的害虫，其若虫、成虫取食葱的叶片、叶心，直接影响了大葱的生产和农户的收益。因此，我们采用空间分布型检验、聚集强度指标检验和线性回归方法调查研究了武威市凉州区大葱田间蓟马种群动态、空间分布及其抽样方法，旨在为金昌武威地区大葱田蓟马科学防治和预测预报提供参考。

1   材料与方法

1.1   调查地点和方法

调查地点为甘肃省武威市凉州区金羊镇松涛村。当地平均海拔1 506 m，年均降水160 mm，土壤类型为薄层灌漠土，耕层土壤有机质含量17.1 g/kg。指示大葱品种为新选章丘大葱，种植密度270 000～3 000 000株/hm2，种植面积900 m2。于2020年4月8日至6月2日在大葱田间进行Z字形5点标记，每点连续调查20株，每株调查全部叶片及叶心，分别逐次统计100个植株上的若虫与成虫总量，视天气状况每周调查1次，制作χ2频次表。

1.2   空间分布型检验

1.2.1   聚集度指标检验   采用扩散系数C、Cassie的CA指数、Lloyd聚集指数M*/、David&Moore丛生指数I以及种群聚集均数λ检验空间分布型[5 - 9 ]。

1.2.2   线性回归方程   检验将平均拥挤度M*與平均密度值做Iwao回归M*=α+β。α为基本扩散指数，β为密度扩散系数。当α > 0，个体间相互吸引，分布的基本成分是个体群;当α=0，分布的基本成分是单个个体;当α < 0，个体间相互排斥。当β=1时，随机分布;当β < 1时，均匀分布;当 β > 1时，聚集分布。将方差S2与平均密度取对数值后做Taylor回归lg（S2）=lga+blg（）。当b=1时，空间分布为随机分布;当b>1时，空间分布为聚集分布;当b趋近于0时，空间分布为均匀分布。

1.3   理论抽样模型和序贯抽样模型

Iwao理论抽样模型n=t2/D2[（α+1）/+ β-1]，n为最适抽样数或理论抽样数，为平均密度，D为允许误差，t为置信度分布值，α、β同Iwao回归模型参数。

Iwao序贯抽样模型T（1、 2）=nm0±t，加号计算可得到害虫密度的上限值T1，减号计算可得到害虫密度的下限值T2。n即抽样数，m0为防治指标，t为置信度分布值，一般取95%置信区间即t=1.96;α、β同Iwao理论抽样模型参数。田间调查时，若累计查得害虫数量大于上限值T1，说明害虫密度高于防治指标，需要进行防治;若累计查得害虫数量小于下限值T2，说明害虫密度低于防治指标，不需要防治;若累计查得害虫数量处于上下限值之间，需继续取样调查。

最大抽样数模型Nmax=t2/d2[（α+1）m0+（β-1）m02）]，d即允许误差D，m0、t、α、β同Iwao序贯抽样模型参数。当田间调查到最大抽样数时，若累计查得害虫数量仍在上下限之间，则根据该点最靠近的界限值判断是否需要防治。

1.4   数据统与分析

采用Excel 2003和DPS17.10软件处理试验数据。

2   结果与分析

2.1   蓟马的种群动态

由图1可见，2020年4月8日至6月2日，蓟马发生的总量随时间推移或气候变化呈2个发生高峰。4月上旬总体气温偏低（平均气温9.9 ℃），蓟马总虫量增长缓慢;4月中旬气温开始上升（平均气温12.2 ℃），总虫量也开始逐渐升高;4月下旬受春季晚霜冻影响（其中4月20 — 25日平均气温9.9 ℃），总虫量下降明显。5月气温大幅上升，其中5月上旬平均气温17.0 ℃，蓟马总虫量也大幅升高;5月中下旬平均气温17.5 ℃，但降水偏多，蓟马总虫量在5月中旬初期达到最高，至5月下旬初期蓟马总虫量降至最低，然后缓慢增加。

2.2   空间分布型检验

从表1可知，2020年4月8日至6月2日蓟马的χ2值均有小于该自由度下奈曼分布、P-E分布、负二项分布P0.05时的χ2值，表示4月8日至6月2日蓟马的实际分布与或奈曼分布模型、P-E分布模型、负二项分布模型显著相符。其中5月23日、28日蓟马的χ2值也小于该自由度下泊松分布P0.05时的χ2值，表示5月23日、28日蓟马的实际分布也与泊松分布模型显著相符。奈曼分布、P-E分布和负二项分布都是聚集分布，因此可以得出2020年4月8日至6月2日蓟马的空间分布型均显著呈聚集分布。

从表2可知，2020年4月8日至6月2日蓟马分布的扩散系数C >1，Cassie指数CA> 0，Lloyd聚集指数M*/> 1，丛生指数I > 0，表示上述蓟马的空间分布型均呈聚集分布;种群聚集均数λ均小于2，表示蓟马聚集是受环境因素决定。种群聚集均数λ和平均密度的值回归模型极显著，方程式是λ=0.919 2-0.078 4（R2=0.931 4），经检验，F=95.00>F0.01，表示蓟马的聚集程度与平均密度极显著正相关。

2.3   线性回归检验

2.3.1   Iwao回归平均拥挤度   M*和平均密度的值回归模型极显著，方程式为M*=0.177 4+ 1.584 8（R2=0.896 7），经检验，F=60.78>F0.01。式中基本扩散指数α=0.177 4>0，表示蓟马个体间相互吸引，分布的基本成分是个体群;密度扩散系数β=1.584 8>1，表示薊马的空间分布型呈聚集分布。

2.3.2   Taylor回归方差   S2和平均密度的对数值的回归模型极显著，方程式为lg（S2）= 0.191 2+1.105 5 lg（）（R2=0.992 7），经检验，F=948.77 > F0.01。式中b=1.105 5 > 1，表示蓟马在大葱田间的分布型呈聚集分布。

2.4   抽样模型

根据Iwao回归模型和Iwao理论抽样模型，一般取95%置信度（即t=1.96），可得出大葱田间蓟马的最适抽样模型为n=3.841 6/D2（1.177 4/+0.584 8）。一般允许误差D可取0.1、0.2和0.3，可得出相应蓟马密度（例如平均密度=1、2、3、4、5、6、7、8、9、10头/株）的最适抽样数（表3）。若蓟马密度相同，则抽样数量随着允许误差的增大而减少;若允许误差相同，则抽样数量随着蓟马密度的增加而递减。

根据Iwao序贯抽样模型，例如选择大葱田间蓟马的防治指标是5头/株，即m0=5.0;取95%置信区间即t=1.96，可得出蓟马的序贯抽样模型方程为T（1、 2）=5n±8.876。例如取调查株数n=10、20、30、40、50、60、70、80、90、100时，可得到相应蓟马数量的序贯抽样表（表4）。在田间调查中，若累计调查的蓟马总量大于表中上限值T1，表示蓟马发生高于防治指标，需要开展防治;若累计调查的蓟马总量小于表中下限值T2，表示蓟马发生低于防治指标，不需要防治;若调查蓟马总量处于T1和T2，仍需继续取样调查。

在95%置信度即t=1.96，根据最大抽样数模型，当允许误差d=0.1时，可得出Nmax=7 878.0，表示防治指标（5.0±0.1）头/株时，田间调查的最大抽样数是7 878个。当允许误差d=0.2时，可得出Nmax=1 969.5，表示防治指标（5.0±0.2）头/株时，田间调查的最大抽样数是1 970个。当允许误差d=0.3时，可得出Nmax=875.33，表示防治指标（5.0±0.3）头/株时，田间调查的最大抽样数是876个。实际应用中，在一定允许误差内调查到最大抽样数时，若累计查得的蓟马总量仍在T1和T2之间，则根据该数值靠近的边界值来决定是否开展防治。

3   结论与讨论

通过空间分布型检验和聚集强度指标检验表明，武威大葱田蓟马空间分布呈聚集分布，其聚集强度随虫口密度的增加而增大，该结论与蓟马在蔬菜、花卉田间的空间分布规律基本一致[10 - 15 ]。本研究发现，蓟马在大葱田间聚集受环境影响较大，而受蓟马本身特性影响较小。建立了大葱田间蓟马若虫与成虫的最适抽样模型n=3.841 6/D2（1.177 4/ +0.584 8）以及序贯抽样模型T（1、 2）=5n±8.876。在实际应用中，可根据蓟马田间序贯抽样表开展虫量调查，以决定是否开展防治。

参考文献：

[1] 曹   雪，朱晓蕾，平淑娥，等.  我国北方地区大葱主要害虫的发生及防治[J].  长江蔬菜，2019（17）：49-53.

[2] 赵干清，杨   峰，吴春峰.  葱蓟马的发生特点及综合防治技术[J].  河南农业， 2007（1）：15.

[3] A M SHELTON，J PLATE，M CHEN.  Advances in control of onion thrips（Thysanoptera：Thripidae） in cabbage[J].  Journal of Economic Entomology，2014（101）：438-443.

[4] 刘忠钱.  大葱蓟马的综合防治措施[J].  福建农业，2008（5）：24.

[5] 李   平，戴   伟.  藜科杂草在洋葱育苗田的空间分布型及其抽样技术[J].  甘肃农业科技，2021，52（4）：49-52.

[6] 李   平.  洋葱根腐病在育苗初期的空间分布型及抽样技术[J].  甘肃农业科技，2021，52（5）：26-29.

[7] 李   平，戴   伟.  潜叶蝇幼虫在二月兰的田间空间分布型及其抽样技术[J].  甘肃农业科技，2021，52（5）：53-57.

[8] 李   平.  冬小麦返青期地下害虫危害空间分布型及其抽样技术[J].  甘肃农业科技，2021，

52（7）：9-13.

[9] 李   平.  设施油白菜地藜科杂草的空间分布型及其抽样技术[J].  甘肃农业科技，2021，

52（7）：63-66.

[10] 郑伯平，郑长英，顾松东，等.  西花蓟马在月季上的空间分布和种群动态[J].  中国农学通报，2012（19）：194-198.

[11] 苑广迪，万岩然，何秉青，等.  叶用莴苣上西花蓟马的防治指标研究[J].  中国蔬菜， 2016（11）：37-39.

[12] 谢国栋，杨成德.   蓟马在临泽蔬菜上的发生与防治[J].  甘肃农业，2014（20）：54.

[13] JIANHUA MO，SCOTT MUNRO，ALAN BOULTON，et al. Within-plant distribution of onion thrips （Thysanoptera： Thripidae） in onions[J].   Journal of Economic Entomology，2008（101）：1331-1336.

[14] KIMBERLEY A STONER，A M SHELTON. Influence of variety on abundance and within-plant distribution of onion thrips （Thysanoptera：Thripidae） on cabbage[J].  Journal of Economic Entomology，1988（101）：1190-1195.

[15] 王厚振，華尧楠，牟吉元.  棉铃虫预测预报与综合治理[M].  北京：中国农业出版社，1999.

（本文责编：郑立龙）

收稿日期：2021 - 03 - 29;修订日期：2021 - 08 - 10

作者简介：李   平（1983 — ），男，陕西西安人，农艺师，硕士，主要从事植物保护研究和推广工作。联系电话：（0）13884093137。Email：274620558@qq.com。