全 文 ：第 49 卷 第 10 期
2 0 1 3 年 10 月
林 业 科 学
SCIENTIA SILVAE SINICAE
Vol. 49，No. 10
Oct．，2 0 1 3
doi:10．11707 / j．1001-7488．20131016
收稿日期: 2013 － 01 － 08; 修回日期: 2013 － 03 － 20。
基金项目: 林业公益性行业科研专项 ( 201004003) ; “十二五”国家科技支撑计划项目 ( 2012BAD19B07 ) ; 国家自然科学基金项目
(31270693)。
* 骆有庆为通讯作者。
基于 GS +的大黄枯叶蛾卵的空间分布关系*
章一巧1 刘永华1 宗世祥1 陆鹏飞1 齐连珍2 骆有庆1
(1．北京林业大学省部共建森林培育与保护教育部重点实验室 北京 100083; 2．陕西省吴起县林业站 吴起 717600)
摘 要: 基于 GS +软件，应用地统计学方法对大黄枯叶蛾卵的空间分布关系进行研究。结果表明: 63. 2%的大
黄枯叶蛾卵块被产在雌茧上，且卵在沙棘上的平均高度为 66. 5 cm，卵距最近雌茧的平均距离为 22. 9 cm; 沙棘林
中大黄枯叶蛾卵的空间依赖范围为 95. 782 m，空间变异值为 0. 680，最佳拟合模型为高斯模型，呈聚集分布; 卵块
的空间分布模拟结果同雌蛹最为接近，两者在空间分布上存在一定的跟随性。大黄枯叶蛾种群密度深受 8 月下旬
至 9 月中旬气温与降水的影响。
关键词: 大黄枯叶蛾; 卵; 空间分布; 地统计学; 气象因子
中图分类号: S718. 7 文献标识码: A 文章编号: 1001 － 7488(2013)10 － 0100 － 06
Spatial Distribution of Trabala vishnou gigantina Eggs with GS +
Zhang Yiqiao1 Liu Yonghua1 Zong Shixiang1 Lu Pengfei1 Qi Lianzhen2 Luo Youqing1
(1 ． Key Laboratory for Silviculture and Conservation of Ministry of Education，Beijing Forestry University Beijing 100083;
2 ． Wuqi County Forestry Station，Shaanxi Province Wuqi 717600)
Abstract: In recent years，Trabala vishnou gigantina ( Lepidoptera: Lasiocampidae) had caused enormous economic
losses in the northwest drought region of China，where Hippophae rhamnoides forests were distributed． As one of the most
important defoliators in that area，larvae of T． vishnou gigantina feeding on tender leaves primarily caused damage to the
crown layers of the trees． Using Software GS + and geostatistics methods，we investigated the spatial distribution of eggs of
T． vishnou gigantina in pure H． rhamnoides forest． Results showed that adult moths preferred to spawn on the female
pupae，with an incidence of 63. 2% ． Eggs gathered mainly at heights of 66. 5 cm on the tree stems． Average distance
between eggs and nearby female pupae was 22. 9 cm． Analysis of isotropic variograms indicated that the distance of spatial
dependence of eggs was 95. 782 m，while the intensity of local spatial continuity was 0. 680． Furthermore，semivariogram
analysis revealed that the Gaussian model was fitted best to eggs distribution，implying an aggregation distribution． The
maps derived from the Kriging interpolation exhibited that the spatial distribution of T． vishnou gigantina eggs was similar
to that of the female pupae． Taking climatic factors into consideration，we found that oviposition of T． vishnou gigantina
was affected by the temperature and rainfall from late-August to mid-September． This research provided a theoretical basis
for controlling T． vishnou gigantina，particularly in predicting their potential distribution areas and future possible
population density．
Key words: Trabala vishnou gigantina; eggs; spatial distribution pattern; geostatistics; climatic factors
大黄枯叶蛾(Trabala vishnou gigantina)是一种
食叶害虫，又名栎黄枯叶蛾、黄绿枯叶蛾，属鳞翅目
(Lepidoptera)枯叶蛾科( Lasiocampidae)黄枯叶蛾属
(Trabala)，国内主要分布在北京、山西、内蒙古、河
南、陕西、甘肃等地区，主要寄主有栎 (Quercus)、板
栗(Castanea mollissima)、核桃( Juglans regia)、苹果
(Malus pumila )、海 棠 ( Malus spectabilis )、沙 棘
(Hippophae rhamnoides)等(刘友樵等，2006)。大黄
枯叶蛾幼虫具有群体危害的特点，常聚集在寄主叶
片背面取食叶肉，轻者造成叶片孔洞和缺刻，严重时
仅残留叶柄; 2008 年沙棘林上大发生时单株虫口密
度达 30 ～ 58 头，有虫株率 100% (吴宗凯，2006)。
该虫的危害导致树势衰弱，果实产量降低，严重威胁
着西北地区的植被恢复与重建。
第 10 期 章一巧等: 基于 GS +的大黄枯叶蛾卵的空间分布关系
空间分布的研究有助于揭示种群动态结构，提
高抽样技术水平和制定虫害防治策略 (郑汉业等，
1993; 丁岩钦，1994; Grego et al．，2006)。自 20 世
纪 80 年代末开始，地统计学(Geostatistics)作为空间
分布研究重要手段之一，在昆虫生态学研究领域已
被广泛应用于定量区域化变量的空间相关性、空间
插值或估值以及分析种群的时空规律性上(Kemp，
1987; Schotzko et al．，1989)。Kemp 等(1987)运用
地统计学探讨蝗虫(Locust)在 110 m 分布范围中的
空间相关性关系; Liehold 等 ( 1991; 1995; 1996;
1998)对舞毒蛾( Lymantria dispar)卵块不同尺度下
半方差图的分析，应用 Kriging 插值法得到密度阀值
概率图从而建立预报模型; 毕守东等(2000)在研究
棉蚜 ( Aphis gossypii) 及天敌龟纹瓢虫 ( Propylaea
japonica)间的种群关系时发现两者在数量和空间上
追随关系明显; 王玉玲(2004)利用地统计学研究发
现栗角斑蚜(Mysoxallis karicola)在板栗林拟合模型
为球型，呈聚集分布; 宗世祥等(2006)运用生物学
方法和 GS 对宁夏地区沙棘木蠹蛾 ( Holcocerus
hippophaecolus)幼虫和蛹的空间特性进行探索并确
定最适样方大小; 季荣等(2006)在 GIS 和 GS 平台
下分析东亚飞蝗(Locusta migratoria manilensis)受环
境影响后卵块的空间异质性和格局分布; 汪文俊等
(2012)应用回归分析、地统计学和聚集强度指数法
对合肥地区李肖叶甲(Ceoporus variabilis)成虫数量
动态及三维空间格局进行研究。可见，地统计学已
经成为生态学领域的新宠。它不仅考虑了种群同环
境的空间相关性和依赖性，避开了传统生物学统计
方法随机选取的理论体系，而且还利用半变异函数
揭示区域化变量的内在联系，结合样点的空间位置
和方向，对样点中稀疏的或无规律的数据进行最优
估值，最大限度地揭示了变量的空间分布规律 (孙
洪泉，1990; 张仁铎，2005)。
笔者针对不同坡向和疏密度的沙棘纯林上大黄
枯叶蛾蛹的空间分布已经作了详细研究(章一巧
等，2012)，但到目前为止，关于大黄枯叶蛾卵的空
间分布特征尚未见报道。鉴于此，本文基于 GS +软
件，利用地统计学方法，研究大黄枯叶蛾卵的空间分
布关系，并结合气象因子(温度、降水量、相对湿度
等)深入分析环境对种群密度的影响，以期进一步
揭示其动态变化规律和扩散趋势，为大黄枯叶蛾虫
害的发生和预警奠定更为全面的理论基础。
1 材料与方法
1. 1 试验地概况
试验地设在陕西省吴起县 (107°38—108°32
N，36°33—37°24E)。吴起是我国沙棘种植大县，
成林沙棘 8. 3 万 hm2，以人工沙棘纯林为主，间有少
量山 杏 ( Armeniaca sibirica )、山 桃 ( Armygdalus
davidiana)(方向军，2006; 吴宗凯，2006)。该县属
黄土高原梁状丘陵沟壑区，海拔 1 233 ～ 1 809 m。
境内为半干旱温带大陆性季风气候，年平均气温
7. 8 ℃，年降水量 478. 4 mm，无霜期 146 天。
于吴起县吴仓堡乡刘堡村场崾岘，选取一块人
工沙棘纯林为样地，大小为 50 m × 40 m，海拔自西
北向东南呈下降趋势。样地内沙棘均为 10 年生，平
均地径 4. 23 cm，树高 170. 09 cm，平均虫口密度为
每株 8. 3 头，有虫株率 11. 1%。
1. 2 研究方法
1. 2. 1 调查方法 采取逐株调查的方法，统计每株
沙棘上的雌雄蛹数、卵块数，同时记录卵块距最近雌
茧的距离，并将每株沙棘的相对空间位置绘制于坐
标纸上。
1. 2. 2 分析方法 1)半方差函数 半方差函数是
区域化变量 Z( xi)和 Z( xi + h)的增量平方的数学期
望，即区域化变量增量的方法。基于地统计学的区
域化变量理论，使用 GS + 8. 0 软件中的半方差函数
分析大黄枯叶蛾卵的空间分布。其公式如下:
γ (h) = 12N(h)Σ
N( h)
i = 1
［Z( xi) － Z( xi + j)］
2。
式中: γ(h)为相对距离 h 的的半方差函数值，N(h)
是被 h 分割的数据对( xi，xi + h)的对数，Z ( xi )和
Z( xi + h)分别是点 xi和 xi + h 处样本的测量值，h 为
两样点间的距离。
半方差函数有 3 个重要参数，即块金常数
(Nugget)、基台值( Sill)和变程(Range)。块金常数
是指变异曲线被延伸至间隔距离为零时的截距，反
映了区域化变量内部随机性的可能程度。基台值是
指达到平衡时的变异函数值，它反映了变量变化的
幅度。变程是指在半方差函数达到平衡或最高时的
间隔距离，它反映了区域化变量的影响范围或空间
依赖性的大小(侯景儒等，1990)。这 3 个参数决定
了半方差函数的形状、结构，最终反映出种群的空间
分布特征。
2)变异函数理论模型的拟合 计算所得的半
方差函数值需应用一定的模型进行拟合，常用来拟
合空间变异曲线的理论模型主要包括线性模型、球
状模型、高斯模型和指数模型。除线性模型揭示种
群随机或均匀分布外，后 3 种模型均表明种群聚集
状态分布(周国法等，1998)。在拟合过程中还需要
结合实际调查情况综合判断，选出最优拟合模型。
101
林 业 科 学 49 卷
根据模拟误差最小原则，依次考虑的因素为决定系
数 R2 大小、残差 RSS 大小、变程 A 和块金值 C0的大
小(宗世祥，2006; 陈强等，2007; 潘杰等，2011;
Ramírez-Davila et al．，2012)。
3)空间分布图的生成 Kriging 法是一种对估
计观测样点间内插值的地统计学方法，是由以 xi为
中心块段的平均值影响范围内的 n 个有效样本值
Z( xi) 的线性组合得到，其公式如下:
Z * ( xo) = ∑
n
i = 1
λ i Z( xi)。
式中: Z * ( xo) 为待估点的值，λ i 是与样点观察
Z( xi) 有关的加权系数，用来表示各个样点 Z( xi)
对 Z * ( xo) 的贡献。
当获得某个变量的半变异函数的拟合模型后，
可对研究区域上未取得样点的区域化变量进行最小
误差估计。借助 Kriging 法对卵、蛹在样地中的密度
进行插值分析，在 Surfer 8 软件中生成空间分布图。
1. 2. 3 气象资料的获取及分析 从陕西吴起气象
局获取 2010，2011，2012 年每年 8 月中旬至 9 月下
旬的气象资料，借助 SPSS 18. 0 软件进行统计分析。
2 结果与分析
2. 1 大黄枯叶蛾卵块的调查情况
大黄枯叶蛾卵多见于茧壳、树干和枝条上，常数
十粒双行相间排列，其上覆有深褐色长毛。解剖 30
头刚羽化的雌蛾发现，雌虫平均怀卵量为 302. 5 粒，
单雌最大怀卵量 360 粒，平均产卵量 297 粒。野外
观察发现: 由于雌虫腹部怀有大量虫卵，使得腹部
较重不利飞行，故羽化交配后雌蛾通常将卵产在雌
茧附近，有未产完的卵时再分批产在同株或者相邻
树上的枝条等处。调查可知: 沙棘样地内有卵株率
达 9. 05%，雌茧上有卵率则为 16. 20%，雄茧上未见
有卵附着; 63. 2% 的卵产在雌茧上，36. 8% 的卵产
在沙棘上; 雌茧上卵粒平均数为 66. 3 个，树干枝条
上卵粒平均数为 56. 0 个; 卵在沙棘上的平均高度
为 66. 5 cm (图 1 )，距最近雌茧的平均距离为
22. 9 cm。
2. 2 大黄枯叶蛾卵的半变异函数及空间分布分析
利用 GS +软件对沙棘样地内卵的空间分布进
行处理，分析得出卵的最优拟合模型为高斯模型
(图 2)，决定系数 R2 为 0. 916，说明卵的数量在沙棘
样地中呈现较明显的空间聚集状态，这与雌蛾飞行
能力弱，交配后将多数卵产在同株或相邻树上密切
相关。卵的空间依赖范围为 95. 782 m，空间变异值
为 0. 680，表明在相应尺度范围内的总变异有
图 1 大黄枯叶蛾卵在沙棘树上不同高度的分布频次
Fig． 1 Height frequency of T． vishnou gigantina eggs in
H． rhamnoides
68. 0%是由空间自相关现象产生的。卵的块金常数
为 277. 000，数值较大的原因可能是人为误差的结
果，即大黄枯叶蛾的初产卵为黑褐色，后逐渐变成铅
灰色，且卵上覆有深褐色长毛，同沙棘黑灰色的枝干
颜色较接近，调查时较难被发现，未被统计在内。
图 2 大黄枯叶蛾卵的变异函数
Fig． 2 Variogram for T． vishnou gigantina eggs
2. 3 大黄枯叶蛾卵与蛹的空间关系
由大黄枯叶蛾卵数量的等值线与矢量叠加图
(图 3a)可以看出:卵在沙棘样地中呈斑块状零星分
布，每个斑块所代表的卵数由聚集中心逐渐向四周
减少，反映出卵的聚集强度和空间尺度，也说明了它
的扩散方向。
为进一步探究卵块数量与雌雄蛹数量之间的空
间分布关系，对比三者的等值线与矢量叠加图 (图
3b，3c)发现: 卵块聚集的区域正是雌雄蛹大量密集
的位置，并且卵块的空间模拟结果同雌蛹最为接近，
两者在空间分布上存在跟随性，分析的结果与大部
分卵被产在雌茧上这一特性相吻合; 大黄枯叶蛾蛹
在沙棘上的平均高度为 66. 7 cm (章一 巧等，
2012)，卵的平均高度为 66. 5 cm，两者近似相等也
证实了该推论。
201
第 10 期 章一巧等: 基于 GS +的大黄枯叶蛾卵的空间分布关系
图 3 大黄枯叶蛾卵( a)、雌(b)雄( c)蛹的等值线与矢量叠加
Fig． 3 Overlay map of isoline and vector for T． vishnou gigantina eggs( a) and pupae( female b and male c)
2. 4 大黄枯叶蛾卵数与气候因子的关系
观察气象因子中的相对湿度、温度及降水量发
现: 2011—2012 年的气温在 8 月中旬至 9 月中旬呈
下降趋势，9 月中旬至 9 月下旬略有回升，而 2010
年 9 月上旬至 9 月下旬期间的温度显著高于 2011，
2012 年同期，最大温差达 7. 5 ℃ ; 降水量在 2010 年
8 月下旬至 9 月下旬明显低于 2011，2012 年同期，
在 9 月上旬差异最为显著，累计毫米数超过 30 mm;
由于相对湿度受降水量影响，变化趋势同降水量近
乎一致。
野外调查结果表明 63. 2%的大黄枯叶蛾卵块产
在雌茧上，故雌茧有卵率可作为衡量成虫产卵情况的
一个指标。2010 年 9 月中旬连降 1 周大雨，温度骤降
7. 5 ℃，此时恰逢大黄枯叶蛾羽化盛期，使得蛹无法
满足羽化的有效积温，羽化期延长，部分羽化出的成
虫还因遭受雨水冲刷造成机械损伤。成虫产卵量明
显下降，使得雌茧上有卵率从 2010 年的 40. 0%骤减
至 2011 年的 16. 2%，平均虫口密度也从每株 8. 3 头
降至每株 0. 35 头，继而影响了卵的空间分布格局。
大黄枯叶蛾虫口密度的变化趋势同 8 月下旬至 9 月
中旬的温度下降曲线和降水量上升曲线相吻合，这说
明其种群密度深受 8 月下旬至 9 月中旬这段羽化交
配产卵期间内气温与降水的影响(图 4)。
图 4 2010—2012 年陕西吴起主要气象因子情况
Fig． 4 Main factors of climate in Wuqi，Shaanxi，2010 － 2012
MA:8 月中旬 Mid Aug． ; LA:8 月下旬 Last Aug． ; ES:9 月上旬
Early Sep． ; MS:9 月中旬 Mid Sep． ; LS:9 月下旬 Last Sep．
301
林 业 科 学 49 卷
3 讨论
了解和掌握害虫的空间分布特性和区域内发生
动态是进行害虫宏观管理的前提和基础。地统计学
作为一种分析昆虫种群动态变化规律和迁移扩散的
有效工具，通过半方差函数分析预测模型，揭示出大
黄枯叶蛾在调查样地内卵聚集分布、蛹随机分布的
空间分布类型; Kriging 插值法作为生物量的整体估
计，有效地解决了空间异质性问题，揭示出大黄枯叶
蛾卵与雌雄蛹的跟随关系，进而反映了成虫就近产
卵的 习 性。陆 永 跃 等 ( 2002 ) 在 研 究 棉 铃 虫
(Helicoverpa armigera)卵的空间分布也同样论证了
卵的空间分布依赖于成虫产卵这一生物学特性。
值得注意的是，在地统计学下探讨空间分布特
性时不能抛却生态学基本理论，要结合昆虫自身的
生物学特性进行全面研究。随着信息技术的发展及
各学科间的交叉，生态研究中还可充分利用地理信
息系统、遥感、全球定位系统、专家系统等多种现代
技术对空间数据强大的处理能力 (李哲等，2004;
郝立武，2011)，结合种内和种间对生物空间分布的
影响，融合时间因素建立种群时空动态模型，综合气
象、地理等生境因子来模拟害虫区域扩散模型，为昆
虫生态学研究、有害生物的监测和决策开辟新的
道路。
目前，国内外已有不少学者对降雨在昆虫种群
密度及空间分布方面开展了研究: 季荣等(2006)发
现产卵盛期遭遇大降雨会影响半方差函数模型的拟
合和东亚飞蝗卵块的空间异质性程度; 赵静等
(2010 ) 研 究 表 明 大 雨 严 重 影 响 了 烟 盲 蝽
(Nesidiocoris tenuis ) 空间分布的理论模型选择;
Marcelo 等(2011)揭示出频繁的降雨和高湿度的环
境条件导致咖啡果小蠹(Hypothenemus hampei)种群
密度的迅速下降; 汪文俊等(2012)也报道了李肖叶
甲成虫的羽化出土时期与雨季的迟早密切相关。影
响昆虫种群动态的因子很多，包括研究区域的地理
位置、海拔、气候条件、天敌、栽培管理等 (张蓉等，
2007)，而探究气象因子等环境变量对昆虫产卵的
影响能直接或间接地揭示种群密度和结构的演变
(赵成章等，2012)。据吴起 2012 年温度和降水量
的气象资料分析，预测 2013 年大黄枯叶蛾在该地的
虫口密度将会有一定程度的增加，当地林业部门应
做好虫害发生的防控准备。
参 考 文 献
毕守东，邹运鼎 ． 2000． 棉蚜及天敌龟纹瓢虫空间格局的地统计学
研究 ． 应用生态学报，11 (3) : 421 － 424．
陈 强，吴伟坚，张振飞，等 ． 2007． 越北腹露蝗若虫空间格局的地
学统计学分析 ． 应用生态学报，18 (2) : 467 － 470．
丁岩钦 ． 1994． 昆虫数学生态学 ． 北京: 科学出版社，22 － 69．
方向军 ． 2006． 吴起县沙棘林分配置模式 ． 国际沙棘研究与开发，4
(3) : 40 － 44．
郝立武 ． 2011． 地理信息系统和地质统计学在昆虫生态学中的应
用 ． 农业工程，1 (3) : 96 － 99．
侯景儒，黄竞先 ． 1990． 地质统计学的理论与方法 ． 北京: 地质出版
社，69 － 78．
季 荣，谢宝瑜，李 哲，等 ． 2006． 基于 GIS 和 GS 的东亚飞蝗卵块
空间格局的研究 ． 昆虫学报，49 (3) : 410 － 415．
李 哲，季 荣，谢宝瑜，等 ． 2004． 论昆虫空间生态学研究 ． 昆虫知
识，40 (1) : 1 － 6．
刘友樵，武春生 ． 2006． 中国动物志 ． 昆虫纲 ．第 47 卷 ．鳞翅目 ． 枯叶
蛾科 ． 北京: 科学出版社，338 － 339．
陆永跃，粱广文 ． 2002． 棉铃虫卵空间分布的地统计学分析 ． 华中农
业大学学报，21 (1) : 13 － 17．
潘 杰，王 涛，宗世祥，等 ． 2011． 红脂大小蠹种群空间格局地统
计学分析及抽样技术 ． 生态学报，31 (1) : 195 － 202．
孙洪泉 ． 1990． 地质统计学及其应用 ． 北京: 中国矿业大学出版社，
38 － 39．
汪文俊，林雪飞，邹运鼎，等 ． 2012． 李肖叶甲成虫数量及三维空间
格局动态 ． 生态学报，32 (14) : 4544 － 4552．
王玉玲，方国飞，汪礼鹏 ． 2004． 栗角斑蚜与红点唇瓢虫种群空间格
局的地统计学分析 ． 安徽农业大学学报，31 (3) : 330 － 334．
吴宗凯 ． 2006． 吴起县沙棘种植区划 ． 国际沙棘研究与开发，4 (3) :
12 － 14．
张仁铎 ． 2005． 空间变异理论及应用 ． 北京: 科学出版社，13．
张 蓉，冷允法，朱猛蒙，等 ． 2007． 基于地统计学和 GIS 的苜蓿斑蚜
种群空间结构分析和分布模拟 ． 应用生态学报，18 (11) : 2580
－ 2585．
章一巧，宗世祥，刘永华，等 ． 2012． 陕西省大黄枯叶蛾蛹的空间分
布 ． 生态学报，32 (20) : 6308 － 6317．
赵成章，李丽丽，王大为，等 ． 2012． 黑河上游天然草地蝗虫空间异
质性与空间分布 ． 生态学报，32 (13) : 4166 － 4172．
赵 静，赵 鑫，王玉军，等 ． 2010． 烟盲蝽及其天敌蜘蛛空间格局的
地统计学分析 ． 生态学报，30 (15) : 4196 － 4205．
郑汉业，夏乃斌 ． 1993． 森林昆虫生态学 ． 北京: 中国林业出版社，
120 － 167．
周国法，徐汝梅 ． 1998． 生物地理统计学———生物种群时空分析的
方法及其应用 ． 北京: 科学出版社，42 － 65．
宗世祥，骆有庆，许志春，等 ． 2006． 沙棘木蠹蛾蛹的空间分布 ． 生态
学报，26 (10) : 3232 － 3237．
宗世祥 ． 2006． 沙棘木蠹蛾生物生态学特性的研究 ． 北京: 北京林
业大学博士学位论文 ．
Grego C R，Vieira S R，Lourenco A L． 2006． Spatial distribution of
Pseudaletia sequax Franclemlont in triticale under no-till
management． Scientia Agricola，63(4) : 321 － 327．
Kemp W P． 1987． Probabiling of outbreak for rangeland grasshoppers
(Orthoptera: Acrididae ) in Montana: application of Markwvian
principle． Environmental Entomology，80(6) : 1100 － 1105．
Liebhold A M，Elkinton J S，Zhou C，et al． 1995． Regional correlation of
401
第 10 期 章一巧等: 基于 GS +的大黄枯叶蛾卵的空间分布关系
gypsy moth ( Lepidoptera: Lymantriidae) defoliation with counts of
egg masses，pupae，and male moths． Environmental Entomology，24
(2) : 193 － 203．
Liebhold A M，Luzader E，Reardon R，et al． 1996． Use of a geographic
information system to evaluate regional treatment effects in a gypsy
moth ( Lepidoptera: Lymantriidae ) management program．
Environmental Entomology，89(5) : 1192 － 1203．
Liebhold A M，Luzader E，Reardon R． 1998． Forecasting gypsy moth
( Lepidoptera: Lymantriidae ) defoloation with a geographic
information system． Environmental Entomology，91(2) : 464 － 472．
Liebhold A M，Zhang X，Hong M E，et al． 1991． Geostastical analysis of
gypsy moth ( Lepidoptera: Lymantriidae ) egg mass populations．
Environmental Entomology，20(5) : 1407 － 1417．
Marcelo C A，Fabio M S，Jair C M，et al． 2011． Geostastical analysis of
the berry borer and leaf miner in a coffee agroecosystem． Precision
Agriculture，12(1) : 18 － 31．
Ramírez-Davila J F，Porcayo-Camargo E，Sánchez-Pale J R，et al． 2012．
Spatial behavior comparision of Bactericera cockerelli Sulc．
(Hemiptera: Triozidae) in Mexico． Neotrop Entomol，41: 9 － 16．
Schotzko D J，O’Keeffe L E． 1989． Geostastical description of the spatial
distribution of Lygus hesperus ( Heteroptera: Miridae) in lentils．
Economic Entomology，82(5) : 1277 － 1288．
(责任编辑 朱乾坤)
501