文章信息
- 朱惠英, 沈平, 吴建华, 王新东, 常承秀, 田炜
- ZHU Huiying, SHEN Ping, WU Jianhua, WANG Xindong, CHANG Chengxiu, TIAN Wei
- 斑膜合垫盲蝽若虫在国槐上的空间分布型及抽样技术
- Spatial distribution pattern and sampling technique for Orthotylus (O.) sophorae nymphs on Sophora japonica
- 生态学报, 2014, 34(4): 832-836
- Acta Ecologica Sinica, 2014, 34(4): 832-836
- http://dx.doi.org/10.5846/stxb201209231341
-
文章历史
- 收稿日期:2012-9-23
- 修订日期:2013-3-15
2. 甘肃省临夏回族自治州农业学校, 临夏 731100
2. Agricultural School of Linxia Hui Nationality Autonomous Region, Gansu Province, Linxia 731100, China
斑膜合垫盲蝽Orthotylus (O.) sophorae Josifov,属半翅目Hemiptera,盲蝽科Mlridae,合垫盲蝽属 Orthotylus[1]。在国内分布于四川省、河南省、甘肃省、陕西省、湖北省、天津市,近年来在临夏地区国槐 Sophora japonica Linnaeus 上发生危害,1年发生1代,6月下旬成虫开始交尾产卵,块产,以卵在国槐主干树皮缝内及2—3年生分枝处丝织的巢内越冬;翌年4月下旬开始孵化。以成虫、若虫聚集于新梢叶背刺吸为害,叶面呈现失绿斑点,致使叶片皱缩干枯,叶柄下垂提前脱落。Josifov[1]根据朝鲜的标本形态描述该种,刘国卿[2]报道斑膜合垫盲蝽为中国新记录种;目前,国内外尚无有关斑膜合垫盲蝽生态学方面的研究报道。2011年对斑膜合垫盲蝽若虫在国槐上的空间分布型进行了调查,在此基础上,提出了操作简便、快捷的田间抽样设计方案,确定出能够在允许误差和可信区间范围内对总体参数进行估计的样本数据,为获取准确、系统的虫情调查资料提供了依据。
1 材料与方法 1.1 试验地概况调查样地位于兰郎公路临夏市西川段(安多园—后杨村口),该路段距离10km,海拔1917.2m,2011年年均气温7.8℃,年日照时数2171.3h,无霜期144d以上,年降水量378mm;栽植地土质坚硬、贫瘠,顺其生长,管理粗放,整体树势衰弱,东西走向,单排栽植,树龄20n,树高5—6m,株距5m,调查前未进行药剂防治。
1.2 取样方法2011年6月上旬,采用线路调查法,在公路北侧每间隔50m设1样段,共设12组样段;样段内每隔1株取1样株,每样段调查10株;每样株从东、西、南、北4个方位,每方位上、中、下分层随机选1年生枝条为样枝,每样枝调查长度为10cm,分别统计不同方位样枝上的若虫和叶片数。
1.3 分布型测定根据调查数据,计算均值(若虫数/叶)(m)、方差(s2)、平均拥挤度 (m*)。依据Lloyd聚块性指数(m*/m)、Beall扩散系数(c)、Waters负二项分布k值、森下正明扩散型指数(Iδ)、David和Moore丛生指数(I)、Cassie聚集度指数(Ca)参数标准判定空间分布型[3, 4]。
1.4 聚集原因检验利用Taylor幂法则[5]和Iwao的m*-m回归模型[6],测定种群的分布格局和内部结构;以Blackith种群聚集均数λ分析聚集原因[3]。
1.5 抽样方法根据Iwao最适理论抽样数模型[4]计算最适抽样数;应用Iwao的方法[4]确定序贯抽样表;依据Gerrard等的零频率模型[9]估计种群密度和理论抽样[10]。
1.6 数据统计分析采用Excel 2003进行数据统计,参数均采用最小二乘法估计,并利用相关系数r进行相关性检验。
2 结果与分析 2.1 分布型检验依据斑膜合垫盲蝽若虫12组样段的均值(m)、方差(s2)和平均拥挤度(m*),计算出各样段的聚块性指数(m*/m)、扩散系数(c)、扩散型指数(Iδ)均大于1,丛生指数(I)、聚集度指数(Ca)大于0(表 1),根据各参数空间分布型判定标准[3, 4],得出斑膜合垫盲蝽若虫在国槐上呈聚集分布。
Waters认为,k值与虫口密度无关,k值愈小,种群的聚集度愈大,若k值趋于∞时(一般在8以上时),则逼近Poisson分布[7]。表 1中的k值均依据矩法求得,且0<k<8,因此斑膜合垫盲蝽若虫的分布型不符合Poisson分布。
样段号 | 均值
m | 方差
s2 | 平均
拥挤度 m* | 聚块性
指数 m*/m | 扩散性
系数 c | 负二项分
布K值 k | 扩散型
指数 Iδ | 丛生指数
I | 聚集度
指数 Ca | 聚集
均数 λ | γ |
1 | 0.0634 | 0.0714 | 0.1899 | 2.9953 | 1.1265 | 0.5012 | 2.9131 | 0.1265 | 1.9953 | 0.0289 | 0.4571 |
2 | 0.1614 | 0.2066 | 0.4417 | 2.7369 | 1.2803 | 0.5758 | 2.7631 | 0.2803 | 1.7369 | 0.0830 | 0.5921 |
3 | 0.1957 | 0.2534 | 0.4905 | 2.5070 | 1.2949 | 0.6636 | 2.5135 | 0.2949 | 1.5070 | 0.1107 | 0.7511 |
4 | 0.0619 | 0.0681 | 0.1619 | 2.6156 | 1.1000 | 0.6190 | 2.6280 | 0.1000 | 1.6156 | 0.0335 | 0.6703 |
5 | 0.2213 | 0.2982 | 0.5688 | 2.5702 | 1.3475 | 0.6369 | 2.5754 | 0.3475 | 1.5702 | 0.1221 | 0.7027 |
6 | 0.2209 | 0.3033 | 0.5941 | 2.6897 | 1.3732 | 0.5918 | 2.6948 | 0.3732 | 1.6997 | 0.1159 | 0.6212 |
7 | 0.2242 | 0.3001 | 0.5625 | 2.5086 | 1.3383 | 0.6629 | 2.5115 | 0.3383 | 1.5086 | 0.1268 | 0.7498 |
8 | 0.0727 | 0.0801 | 0.1741 | 2.3958 | 1.1015 | 0.7164 | 2.4053 | 0.1015 | 1.3958 | 0.0430 | 0.8468 |
9 | 0.1161 | 0.1456 | 0.3702 | 3.1885 | 1.2541 | 0.4569 | 3.1970 | 0.2541 | 2.1885 | 0.0479 | 0.3770 |
10 | 0.2614 | 0.3775 | 0.7057 | 2.7000 | 1.4443 | 0.5882 | 2.7030 | 0.4443 | 1.7000 | 0.1366 | 0.6147 |
11 | 0.2266 | 0.3005 | 0.5523 | 2.4370 | 1.3257 | 0.6959 | 2.4399 | 0.3257 | 1.4370 | 0.1318 | 0.8095 |
12 | 0.2308 | 0.3161 | 0.6003 | 2.6014 | 1.3696 | 0.6245 | 2.6049 | 0.3696 | 1.6014 | 0.1257 | 0.6803 |
利用Bliss和Owen提出的用矩法来估公共kc值[8],计算得公共kc值为0.6169,kc值的χ2值为5.1098,且df为11,经χ2(χ20.05=19.675)检验,即χ20.05=19.675>χ2=5.1098,其差异不显著(P>0.05),因此认为斑膜合垫盲蝽若虫分布符合负二项分布。
2.2 聚集原因分析利用Taylor幂法则s2=αmb,即lgs2=lgα+ blgm[5],以lgs2和lgm进行线性回归(当df=10时,即 r =0.9995>α0.05=0.5760)得出方差与平均密度的回归方程:
式中,lgα=0.2249>0,b=1.1622>1,s2=1.0274m1.1622,说明种群在一切密度下均是聚集分布,且具密度依赖性,聚集强度随种群密度的升高而增强。
利用Iwao的m*-m线性回归模型m*=α+βm[6],得出斑膜合垫盲蝽若虫在国槐上的线性回归方程为m*=0.0224+2.5011m。12组抽样数据的统计结果(当df=10时,即 r =0.9020>α0.05=0.5760)为α=0.0224>0,说明斑膜合垫盲蝽若虫个体间相互吸引,存在个体群;β=2.5011>1,表明基本成分呈聚集分布。测定结果与聚集度指标检验相一致。
应用Blackith种群聚集均数λ=mγ/2k分析聚集原因[3]。
式中,m为平均虫口密度,k为负二项分布的指数,γ是自由度等于2k,以0.5概率值对应处的χ2值,应用内插法得出。经计算,12组样段的总λ=0.0922,因此根据Blackith对昆虫的聚集原因分析可知:λ=0.0922<2,说明斑膜合垫盲蝽若虫聚集可能是由某些环境作用引起的。
2.3 抽样技术 2.3.1 理论抽样数技术应用Iwao最适理论抽样数模型N=(t/D)2[(α+1)/m+β-1][4],确定出最适抽样数公式:N=(t/D)2(1.0224/m+1.5011)
取t=1.96(置信概率为0.95),相对误差(D=d′/m)为0.1、0.2和0.3时,平均若虫数(m)在0.5、0.10、0.15、0.20、…、0.80的条件下,得出斑膜合垫盲蝽若虫的最适抽样数(N)(表 2)。
相对误差
Relative variation (D) | 不同若虫密度下的最适抽样数Optimal sampling number at different density of nymphs | ||||||||||
0.05 | 0.10 | 0.15 | 0.20 | 0.25 | 0.30 | 0.40 | 0.50 | 0.60 | 0.70 | 0.80 | |
0.1 | 8432 | 4504 | 3195 | 2541 | 2148 | 1886 | 1559 | 1362 | 1231 | 1137 | 1068 |
0.2 | 2108 | 1126 | 799 | 635 | 536 | 472 | 389 | 341 | 307 | 284 | 266 |
0.3 | 937 | 501 | 355 | 282 | 239 | 209 | 173 | 151 | 137 | 126 | 118 |
由表 2可知,在相对误差(D)固定时,斑膜合垫盲蝽若虫的最适抽样数(N)随平均若虫数(m)的增大而下降,只是下降幅度逐渐缩小;在相同的平均若虫数(m)下,随着相对误差(D)的增大,所需抽样数(N)依次减少。
2.3.2 序贯抽样技术根据Iwao的方法[4],设临界防治密度为m0,则:, 取t=1.96(置信概率为0.95),已知α=0.0224,β=2.5011,初步拟定防治指标m0为0.3头/叶,得:
用上述两式,计算出斑膜合垫盲蝽若虫序贯抽样表(表 3)供防治参考。
调查株数(n)
Tree number | 累计若虫数/头
The accumulative amount of nymphs | |
上限(T′(n))Upper | 下限(T″(n))Lower | |
50 | 24 | 6 |
60 | 28 | 8 |
70 | 32 | 10 |
80 | 36 | 12 |
90 | 39 | 15 |
100 | 43 | 17 |
110 | 47 | 19 |
120 | 50 | 22 |
130 | 54 | 24 |
140 | 57 | 27 |
150 | 61 | 29 |
160 | 65 | 32 |
170 | 68 | 34 |
180 | 72 | 37 |
最大抽样数:nmax=(t/d)2(1.0224m0+1.5011m20)
d为m=m0时估计密度所允许的置信限,当d=0.1时,最大抽样数为170株;d=0.2时,最大抽样数为42株。
2.3.3 零频率估计法根据Gerrard等的零频率回归模型:x=α(-lnP0)b[9],得出膜合垫盲蝽若虫的平均密度(x)与零频率(P0)回归模型:
本试验总样本的零样本频率P0=0.8715,求得x=0.1923,而实测值m为0.1714,相对误差(D= m-x /x)仅0.1089,可使用该模型估计膜合垫盲蝽若虫的平均密度。
根据Gerrard的方法对种群进行理论抽样(n)[10],其公式为:
式中,β为模型中-lnP0的指数,P0为零频率,D抽样精度。已知P0=0.8715,D=0.2时,理论抽样数(n)为240.9,即需要调查241片叶是否有虫,求得P0值后代入Gerrard模型x=1.7457(-lnP0)1.1119,即可得出每叶的虫口密度。
将田间按理论抽样数(n)调查取得的P0值代入已知Gerrard模型进行验证,得出理论值x与实测值m的比较结果(表 4)。
样段1
Sample 1 | 样段2 | 样段3 | 样段4 | 样段5 | 样段6 | 样段7 | 样段8 | 样段9 | 样段10 | 样段11 | 样段12 | |
1 )相对误差(精度) D= m-x /x | ||||||||||||
P0 | 0.9417 | 0.8637 | 0.8477 | 0.9421 | 0.8286 | 0.8296 | 0.8257 | 0.9332 | 0.9016 | 0.8108 | 0.8236 | 0.8248 |
m | 0.0634 | 0.1614 | 0.1957 | 0.0619 | 0.2213 | 0.2209 | 0.2242 | 0.0727 | 0.1161 | 0.2614 | 0.2266 | 0.2308 |
x | 0.0765 | 0.2063 | 0.2358 | 0.0759 | 0.2723 | 0.2703 | 0.2779 | 0.0894 | 0.1404 | 0.3074 | 0.2820 | 0.2796 |
D | 0.1717 | 0.2176 | 0.1703 | 0.1850 | 0.1873 | 0.1828 | 0.1930 | 0.1871 | 0.1729 | 0.1497 | 0.1964 | 0.1747 |
由表 4可见,该模型的抽样精度除样段2略大于预定D值0.2,其余11个样段均小于预定D值0.2。因此利用Gerrard模型x=1.7457(-lnP0)1.1119,根据0—1抽样所得P0值,可有效估计国槐叶片上膜合垫盲蝽若虫的种群密度。
3 讨论本文就斑膜合垫盲蝽若虫的12组样段的调查结果进行了空间分布型检验分析,结果表明,斑膜合垫盲蝽若虫在田间呈聚集分布;其公共kc值为0.6169,符合负二项分布;聚集原因是由某些环境作用引起的;这主要是由于斑膜合垫盲蝽若虫对国槐为害部位(新梢叶背)的选择性所致,其基本上受国槐新梢长势的影响,趋向于阳面树冠外围的新梢。
在虫情调查时,以叶片为抽样单位,在应用本文抽样技术研究确定的模型时,采用随机抽样,因为模型建立的基础是随机抽样。利用本研究确定的斑膜合垫盲蝽若虫理论抽样数,根据精确度与工作量,在允许误差和可信区间范围内,即在试验研究时,相对误差要求严格取D=0.1,在大田调查时相对误差可放宽到D=0.2或0.3,确定出所需最适抽样株数。用序贯抽样表进行抽样时,当调查的株数,累计若虫数超过上限时,就必须进行防治,若累计若虫数小于下限时,则不需要进行防治;累计若虫数在上限和下限之间时,应继续增加抽样株数。当不易做出判断时,应需确定最大抽样数目。Gerrard等的零频率回归模型,被广泛应用于一些体小的有害生物(如蚜虫、木虱、螨类等)种群密度的估计[11];在应用时首先需满足理论抽样数以保证抽样的精度,在调查过程中不需要检查所抽取样本的虫口密度,只需记录有虫样数和无虫样数,根据建立的零频率模型,即可估计该环境寄主上的虫口密度;用零频率法估计昆虫田间种群密度,由于避免了直接计数的烦琐过程,因而具有操作简便和效率高的优点[12]。因此,在保证一定的虫情调查准确度情况下,利用合理简便的抽样设计方案进行调查就可取得较准确的虫情资料,为斑膜合垫盲蝽制定有效的防治措施提供参考依据。
致谢: 承蒙南开大学刘国卿教授鉴定斑膜合垫盲蝽种名并提供相关文献资料;中国农业大学彩万志教授和本站张山林站长帮助写作,特此致谢。
[1] | Josifov M.Drei neue Orthylus-Arten aus Korea (Heteroptera, Miridae). Reichenbachia, 1976, 16(11):143-146. |
[2] | Liu G Q.A study on systematics of the subfamily orthotylinae from China. Tianjin: Nankai University, 1998. |
[3] | Ding Y Q.Principle and Application of Insect Mathematical Population Ecology. Beijing: Science Press, 1980, 84-124. |
[4] | Xu R M. Population Ecology of Insect., Beijing Normal University Press, 1988, 8-60. |
[5] | Taylor L R.Aggregation,variance and the mean. Nature, 1961, 189:732-735. |
[6] | Iwao S, E Kuno.An approach to the analysis of aggregation pattern in biological populations//Patil G P, Pielou E C,Waters W E. Statistical Ecology. The Pennsylvania State University Press, Pennsylvania, 1971,1:461-513. |
[7] | Waters W E.A quantitative measure of aggregation in insects. Journal of Economic Entomology, 1959, 52:1180-1184. |
[8] | Bliss C I, Owen A R G.Negative binomial distributions with a common K. Biometrika, 1958, 45:37-58. |
[9] | Gerrard D J, H C Chiang.Density estimation of com rootworm egg populations based upon frequency of occurrence.Ecology,1970, 51 (2):237-245. |
[10] | Ding Y Q. Mathematical Ecology of Insects.Beijing. Beijing: Science Press, 1994,126-134. |
[11] | Wang X G,Liu S S.Development of binomial sampling models for predicting insect population density. Journal of Biomathematics,1996, 11(3):45-51. |
[12] | Cai D C,Li W,Xiao D M.Estimation of the larval population density of Liriomyza sativae on tomato leaves by zero frequency method.Plant Protection, 2011, 37(2):159-161. |
[2] | 刘国卿.中国合垫盲蝽亚科系统学研究.天津:南开大学, 1998. |
[3] | 丁岩钦.昆虫种群数学生态学原理与应用.北京:科学出版社, 1980, 84-124. |
[4] | 徐汝梅.昆虫种群生态学.北京:北京师范大学出版社, 1988, 8-60. |
[10] | 丁岩钦.昆虫数学生态学.北京:科学出版社, 1994, 126-134. |
[11] | 汪信庚,刘树生.昆虫种群密度的二项抽样估计模型研究进展.生物数学学报,1996, 11(3):45-51. |
[12] | 蔡笃程,李炜,肖冬梅.番茄叶片上美洲斑潜蝇幼虫种群密度的零频率估计法.植物保护,2011, 37(2):159-161. |