全 文 ：中国生态农业学报 2015年 7月 第 23卷 第 7期
Chinese Journal of Eco-Agriculture, Jul. 2015, 23(7): 906913


* 上海市农业领军人才计划项目、上海市自然科学基金青年项目(12ZR1449100)、上海市农委青年人才成长计划项目(沪农青字(2014)
第 1-22号)和上海市青年科技启明星计划项目(15QB1402800)资助
** 通讯作者: 万年峰, 研究方向为昆虫生态学及农业生态学。E-mail: fnianwan_2004@163.com
蒋杰贤, 主要从事昆虫生态学及生物防治研究。E-mail: jiangjiexian@163.com
收稿日期: 20150114 接受日期: 20150507
http://www.ecoagri.ac.cn
DOI: 10.13930/j.cnki.cjea.150066
基于地统计学和 GIS的化学农药胁迫下桃树
桃蚜与草间钻头蛛种群空间格局*
蒋杰贤 万年峰** 季香云
(上海市农业科学院生态环境保护研究所/上海市设施园艺技术重点实验室 上海 201403)
摘 要 化学农药过量施用不仅杀死了害虫和自然天敌, 也影响天敌与害虫之间的空间关系。桃园生态系统
中草间钻头蛛(Hylyphantes graminicola)是桃蚜(Myzus persicae)的重要捕食性天敌, 本文系统调查了化学农药
长期胁迫下不同时期(4月中旬至 9月上旬)桃树桃蚜和草间钻头蛛种群数量, 并运用地统计学方法和地理信息
系统(GIS)分析了桃树桃蚜和草间钻头蛛种群的空间结构, 并采用基于高斯模型、指数模型、球型模型和圆型
模型的普通克立格插值法模拟了其种群空间分布格局。结果表明, 桃蚜和草间钻头蛛种群在 10次调查期间均
呈随机分布格局; 桃蚜和草间钻头蛛空间分布距离相关性较弱, 其变动范围均为 6.863 0~43.174 1 m; 桃蚜和草
间钻头蛛的空间结构比例分别为 0.788 8~0.983 9和 0.811 6~0.980 6, 块金值分别为 0.254 2~4.896 3和 0.218 4~0.749 9,
偏基台值分别为 0.010 5~0.250 0和 0.004 8~0.075 7。本研究结果显示, 长期使用化学农药致使桃树桃蚜和草间
钻头蛛种群的空间结构比例值均大于 0.7500, 这表明化学农药胁迫下桃蚜和草间钻头蛛种群之间的跟随效应
不明显, 草间钻头蛛对桃蚜的捕食作用不强。
关键词 化学农药胁迫 桃蚜 草间钻头蛛 空间格局 跟随效应 捕食作用
中图分类号: Q968.1; S436 文献标识码: A 文章编号: 1671-3990(2015)07-0906-08
Geostatistics and GIS analysis of spatial patterns of Myzus persicae and
Hylyphantes graminicola in peach orchards under chemical pesticide stress
JIANG Jiexian, WAN Nianfeng, JI Xiangyun
(Shanghai Key Laboratory of Protected Horticultural Technology / Eco-environment Protection Research Institute, Shanghai
Academy of Agricultural Sciences, Shanghai 201403, China)
Abstract The overuse of chemical pesticides not only kills insect pests and natural enemies, but also affects the spatial relationship
between insect pests and their natural enemies. Hylyphantes graminicola is a key predator of Myzus persicae in peach orchard
eosystems. In this paper, we systematically investigated M. persicae and its predator H. graminicola populations in peach orchards
under chemical pesticide stress at different times (from mid April to early September), used geostatistics and geographic information
system (GIS) to analyze the spatial structure, and used ordinary Kriging interpretation with Gaussian, Exponential, Spherical and
Circular models to simulate the spatial distribution of the two species.The aim of the study was to understand the spacial distribution
of M. persicae and its predator H. graminicola under long-term applicaiton of chemical pestcides and provide the theoretical support
for ecological control of peach garden pests. The results suggested that both H. graminicola and M. persicae had random spatial
arrangement within 10 iterative times of investigations. The proportions of spatical sturcture [C0/(C+ C0)] of populations of H.
graminicola and M. persicae were 0.788 80.983 9 and 0.811 60.980 6 indicating weak spacial relathship bewteen two populaitons
under long-term chemical pestcide stress. The nugget values and partial sills of H. Graminicola, M. persicae were respectivel 0.254 2
4.896 3, 0.218 40.749 9 and 0.010 50.250 0, 0.004 80.075 7, respectively, indicating random spatial arrangement, too. The
第 7期 蒋杰贤等: 基于地统计学和 GIS的化学农药胁迫下桃树桃蚜与草间钻头蛛种群空间格局 907


http://www.ecoagri.ac.cn
correlation of spatial distribution distance of two species was relatively weak, the distance ranges for the two species was 6.863 0
43.174 1 m. Though the model parameters of semivariograms for M. persicae and H. graminicola at different times were changed
greatly due to population density, temperature, and peach growth, the spatial patterns of the two populations were random under
long-term chemical pesticide stress. Our study demonstrated that under long-term and excessive use of chemical pesticides in
peach orchards, the proportions of H. graminicola and M. persicae spatial structure was greater than 0.750 0, and their spatical
patterns were random. This suggested that there was no obvious effect of H. graminicola on M. persicae and that H.
graminicola had minimal predatory effect on M. persicae.
Keywords Chemical pesticide stress; Myzus persicae; Hylyphantes graminicola; Spatial pattern; Following effect; Predatory
effect
(Received Jan. 14, 2015; accepted May 7, 2015)
地统计学(Geostatistics)和地理信息系统(Geographic
Information System, GIS)技术已广泛运用于生态农
业中[12], 为分析昆虫种群时空动态提供了新方法、
新工具[34]。较传统空间格局研究方法, 地统计学方
法的优点是考虑了样点的位置和彼此距离, 直接测
定了空间结构相关性和依赖性, 可用于研究有一定
随机性和结构性的各种变量的空间分布规律[5]; GIS
技术的优势是可以进行昆虫空间显示与查询、影
响种群动态的因子分析、可视化种群时空动态模
拟 [67]。运用地统计学和GIS技术可以充分利用野外
调查所提供的各种信息及稀疏或无规律的空间数据,
通过空间插值方法形成种群在时间序列上的空间分
布格局[7]。研究表明, 自然状态下李子园草间钻头蛛
(Hylyphantes graminicola, 原名草间小黑蛛)种群对
桃蚜(Myzus persicae)种群在空间位置上具有追随关
系 , 对其种群数量的控制有较大作用 [8]; 桃园生草
模式下桃树草间钻头蛛对桃蚜在空间位置上具有追
随关系, 对其种群数量的控制有较大作用, 草间钻
头蛛是桃蚜的优势天敌[9]。
截至目前 , 联合运用地统计学和GIS研究昆虫
种群空间结构已有报道 [910], 而研究化学农药长期
胁迫对害虫及其天敌种群空间格局的报道鲜见。为
此, 本文应用该方法研究了长期使用化学农药模式
下不同时期桃树桃蚜与其捕食性天敌草间钻头蛛种
群空间结构空间相关性, 以期深入了解其在化学农
药胁迫下桃园中的具体空间分布及其动态, 为桃园
生态修复控害机制研究奠定理论基础。
1 研究区域与研究方法
1.1 研究区概况
试验在上海市浦东新区新场镇果园村水蜜桃基
地(31.03°N, 121.41°E, 海拔 4.3 m)进行。桃园地处
长江三角洲冲积平原前沿, 属北亚热带南缘东亚季
风盛行地区, 距离长江约 20 km。为防治有害生物
危害, 桃园周年性化学农药有效成分使用量为 9.98~
14.80 kg·hm2, 主要为杀灭菊酯、氯氰菊酯、吡虫啉、
乙酰甲胺磷、敌敌畏、灭幼脲、多菌灵、百菌清等
化学农药。其他有害生物管理措施包括: 冬季人工
修剪受侵染的桃枝、6 月中下旬果实套袋、桃休眠
期喷施波美度 4~5 的石硫合剂。该种施药模式长达
9年以上。试验区面积 2 hm2, 水蜜桃品种‘新凤蜜露’,
8年树龄为主, 树高 2~2.5 m、株行距为 4 m× 4 m。
1.2 调查方法
在桃园以桃树南北为行, 采用4 m4 m网络化
取样调查方法 [1113], 共抽取样树 60 株, 分5行每行
12株[9]; 在每棵桃树树冠东、南、西、北4个方位和
上、中、下3层各取1根代表性枝条; 调查每根枝条
从梢部向内0.30 m长度内桃蚜和草间钻头蛛的种群
数量, 同时对样地情况作详细记录[1415]。调查时间
为4月中旬至9月上旬, 每隔15~20 d调查1次, 共调
查10次。调查期间, 桃园按常规管理措施进行农事
操作, 但调查当天桃园不施药。
1.3 分析方法
1.3.1 空间结构分析
根据区域化变量理论[1617], 计算草间钻头蛛和
桃蚜种群的试验半变异函数(Semivariogram), 拟合
其理论半变异函数, 分析半变异函数结构, 来描述
草间钻头蛛和桃蚜种群的空间结构和空间相互依赖
关系。对于观察的数据系列Z(xi) (i=1, 2, 3, ⋯, n), 样
本半变异函数值 (h)可用下式计算:
( )
1
1( ) [ ( ) ( )]
2 ( )
N h
i i
i
h Z x Z x h
N h


   (1)
式中: N(h)是被h分割的数据对(xi, xi+h)的对数, Z(xi)
和Z(xi+h)分别是在点xi和xi+h处样本的测量值, h是
分割两样点的距离。
块金值(C0, nugget)、变程(, range)和基台值
(C0+C, sill)是半变异函数的3个重要参数 , 可反映
昆虫种群空间格局及空间相互依赖关系, 同时还给
出种群间的空间相关范围[1617]。变程是指使半变异
908 中国生态农业学报 2015 第 23卷


http://www.ecoagri.ac.cn
函数达到平衡 (基台值)时的h值 , 说明该昆虫种群
在此范围内存在空间相关性 ; 基台值指达到平衡
时的半变异函数值 ; 块金值指变异函数曲线被延
伸至间隔距离为零时的截距 , 表示随机部分的空
间异质性。
在空间上昆虫种群数量是区域化变量, 因而可
用区域化变量理论和方法进行研究。通过计算桃蚜
和草间钻头蛛种群的半变异函数曲线和选择适合的
半变异函数模型, 分析其种群的空间格局及空间相
关关系。在ArcGIS 9.3软件平台的Geostatistical
Analysis 模块上进行空间结构分析。
1.3.2 构建最优半变异函数模型
根据半变异图的性质和特征选择模型。ArcGIS
软件中常见的模拟模型为球型模型、指数模型、高
斯模型和圆型模型[1819], 具体选择时主要根据模拟
误差最小原则, 即平均预测误差(mean error)尽可能
接近0、均方根误差(root mean square error)尽可能
小、平均克立格标准差(average standard deviation)
尽可能小、无偏估计(mean standardized)尽可能接近0、
一致性估计(root mean square standardized)接近1[20]。
球型、指数型半变异函数表明该种群空间分布属
聚集分布; 非水平直线型的半变异函数表明该种群是
中等程度的聚集分布, 其空间依赖范围超过研究尺
度。如果是随机分布, 则 (h)随距离无一定规律性变化,
完全随机或均匀的数据,  (h)呈水平直线或稍有斜率,
表明在抽样尺度下没有空间相关性或空间依赖性[17]。
揭示变量的空间相关性程度的参数是空间结构
比例(即块金值与基台值之比)。若比值<0.25, 表明
变量具有强烈的空间相关性; 若比值>0.75, 则空间
相关性很弱; 若比值介于0.25~0.75, 表明变量具有
中等的空间相关性[2122]。
1.3.3 空间分布模拟
在用空间统计学的理论和方法解决未知点性状
值的预测或估计问题时, 用Kriging插值法定量地分
析所采集样本性状值的空间分布特征, 再用半变异
函数确定权重后进行拟合[7]。插值的生成在ArcMap
9.3软件平台的Spatial Analyst模块支持下进行。
2 结果与分析
2.1 化学农药胁迫下桃树桃蚜和草间钻头蛛种群
空间分布的半变异函数和最优模型的构建
在 ArcMap 软件平台下 , 用普通克立格空间插
值法分别模拟上述4个半变异函数模型 , 分别得出
不同时期桃蚜和草间钻头蛛种群空间分布的最优半
变异函数模型(表 1)。
表 1 用普通克立格空间插值法计算的桃蚜和草间钻头蛛种群空间分布的最优半变异函数模型的误差
Table 1 Error of the optimal semivariogram models with ordinary Kriging interpolation for the spatial distribution of Myzus
persicae and Hylyphantes graminicola populations
取样日期(月-日)
Sampling date
(month-day)
物种
Species
模型类型
Model type
平均预
测误差
Mean error
均方根误差
Root mean
square error
平均预测标准差
Average standard
deviation
无偏估计
Mean
standardized
一致性估计
Root mean square
standardized
桃蚜 M. persicae 圆型 Circular 0.000 4 1.074 0 1.127 0 0.000 7 0.954 7 04-13
草间钻头蛛 H. graminicola 指数 Exponential 0.013 0 0.520 6 0.511 4 0.024 1 1.019 0
桃蚜 M. persicae 球型 Spherical 0.019 5 2.355 0 2.295 0 0.009 0 1.026 0 04-29
草间钻头蛛 H. graminicola 指数 Exponential 0.008 9 0.523 9 0.537 9 0.014 5 0.974 8
桃蚜 M. persicae 圆型 Circular 0.000 5 0.850 0 0.828 1 0.000 1 1.027 0 05-12
草间钻头蛛 H. graminicola 指数 Exponential 0.002 9 0.513 7 0.535 8 0.004 2 0.958 6
桃蚜 M. persicae 高斯 Gaussian 0.008 9 0.729 6 0.768 8 0.010 3 0.951 7 05-28
草间钻头蛛 H. graminicola 高斯 Gaussian 0.003 5 0.563 9 0.589 1 0.005 1 0.959 9
桃蚜 M. persicae 高斯 Gaussian 0.035 2 1.306 0 1.245 0 0.028 0 1.048 0 06-14
草间钻头蛛 H. graminicola 指数 Exponential 0.017 8 0.653 0 0.611 9 0.026 5 1.065 0
桃蚜 M. persicae 高斯 Gaussian 0.024 1 0.962 6 0.922 3 0.024 8 1.043 0 06-30
草间钻头蛛 H. graminicola 圆型 Circular 0.001 6 0.674 3 0.698 6 0.002 5 0.967 7
桃蚜 M. persicae 高斯 Gaussian 0.006 4 0.543 4 0.525 3 0.011 4 1.034 0 07-15
草间钻头蛛 H. graminicola 指数 Exponential 0.021 1 0.596 6 0.575 7 0.035 4 1.036 0
桃蚜 M. persicae 高斯 Gaussian 0.008 8 0.538 0 0.522 8 0.016 4 1.028 0 07-30
草间钻头蛛 H. graminicola 球型 Spherical 0.011 4 0.604 3 0.651 3 0.016 7 0.932 8
桃蚜 M. persicae 指数 Exponential 0.024 5 0.654 9 0.607 3 0.038 5 1.076 0 08-16
草间钻头蛛 H. graminicola 圆型 Circular 0.008 0 0.807 6 0.797 6 0.009 6 1.011 0
桃蚜 M. persicae 圆型 Circular 0.011 2 1.177 0 1.199 0 0.007 4 0.982 4 09-03
草间钻头蛛 H. graminicola 球型 Spherical 0.016 8 0.938 4 0.943 5 0.015 9 0.995 1
第 7期 蒋杰贤等: 基于地统计学和 GIS的化学农药胁迫下桃树桃蚜与草间钻头蛛种群空间格局 909


http://www.ecoagri.ac.cn
由表 1可以看出: 1)桃蚜种群空间分布在 4月 13
日的圆型函数模型平均预测误差和无偏估计最接近
0, 而草间钻头蛛在 6 月 30 日的圆型函数模型平均
预测误差和无偏估计最接近 0; 2) 桃蚜种群空间分
布在 4 月 29 日的球型函数模型一致性估计最接近
1, 而草间钻头蛛在 8 月 16 日的圆型函数模型一致
性估计最接近 1; 3)桃蚜种群空间分布的圆型函数模
型均方根误差最小值(0.538 0)在 7月 30日, 而草间
钻头蛛种群空间分布的球型模型均方根误差最小值
(0.513 7)在 5月 12日; 4)桃蚜和草间钻头蛛种群空间
分布的圆型函数模型平均预测标准差最小值分别在
7月 30日(0.522 8)、4月 13日(0.511 4)。
2.2 化学农药胁迫下桃树桃蚜和草间钻头蛛种群
的空间结构
化学农药长期大量施用后桃树桃蚜和草间钻头
蛛种群的空间结构比例值均大于0.7500(表2), 这表
明化学农药胁迫模式下桃蚜和草间钻头蛛种群的空
间相关性很弱, 化学农药长期胁迫致使草间钻头蛛
种群对桃蚜种群的空间依赖性不强、草间钻头蛛种
群对桃蚜种群的捕食效应不明显。10次调查数据分
析结果显示, 桃蚜和草间钻头蛛种群空间格局均表
现为随机分布(表 2、图 1、图 2)。
表 2 不同时期桃蚜和草间钻头蛛种群数量的半变异函数模型参数及空间格局
Table 2 Model parameters of semivariograms and spatial patterns for Myzus persicae and Hylyphantes graminicola at different
sampling times
日期(月-日)
Date
(month-day)
物种
Species
模型类型
Model type
变程()
Range
(m)
块金值
(C0)
Nugget
偏基台值
(C)
Partial sill
空间结构比例
[C0/(C+ C0)]
Proportion of
spatial structure
空间格局
Spatial pattern
桃蚜 M. persicae 圆型 Circular 6.863 0 0.919 3 0.246 2 0.788 8 随机分布 Random distribution04-13
草间钻头蛛 H. graminicola 指数
Exponential
43.174 1 0.243 2 0.004 8 0.980 6 随机分布 Random distribution
桃蚜 M. persicae 球型 Spherical 43.174 1 4.896 3 0.182 1 0.964 1 随机分布 Random distribution04-29
草间钻头蛛 H. graminicola 指数
Exponential
6.863 0 0.232 7 0.021 6 0.915 1 随机分布 Random distribution
桃蚜 M. persicae 圆型 Circular 32.859 7 0.639 7 0.010 5 0.983 9 随机分布 Random distribution05-12
草间钻头蛛 H. graminicola 指数
Exponential
6.863 0 0.218 4 0.036 1 0.858 2 随机分布 Random distribution
桃蚜 M. persicae 高斯 Gaussian 6.863 0 0.479 1 0.048 5 0.908 1 随机分布 Random distribution05-28
草间钻头蛛 H. graminicola 高斯 Gaussian 6.863 0 0.280 0 0.030 4 0.902 1 随机分布 Random distribution
桃蚜 M. persicae 高斯 Gaussian 43.174 1 1.440 6 0.250 0 0.852 1 随机分布 Random distribution06-14
草间钻头蛛 H. graminicola 指数
Exponential
43.174 1 0.326 2 0.075 7 0.811 6 随机分布 Random distribution
桃蚜 M. persicae 高斯 Gaussian 43.174 1 0.792 9 0.084 5 0.903 7 随机分布 Random distribution06-30
草间钻头蛛 H. graminicola 圆型 Circular 6.863 0 0.382 9 0.053 8 0.876 8 随机分布 Random distribution
桃蚜 M. persicae 高斯 Gaussian 43.174 1 0.257 0 0.031 9 0.889 6 随机分布 Random distribution07-15
草间钻头蛛 H. graminicola 指数
Exponential
43.174 1 0.300 2 0.030 8 0.906 9 随机分布 Random distribution
桃蚜 M. persicae 高斯 Gaussian 43.174 1 0.254 2 0.041 6 0.859 4 随机分布 Random distribution07-30
草间钻头蛛 H. graminicola 球型 Spherical 6.863 0 0.336 4 0.038 0 0.898 5 随机分布 Random distribution
桃蚜 M. persicae 指数 Exponential 43.174 1 0.335 8 0.028 9 0.920 8 随机分布 Random distribution08-16
草间钻头蛛 H. graminicola 圆型 Circular 6.863 0 0.498 8 0.070 5 0.876 2 随机分布 Random distribution
桃蚜 M. persicae 圆型 Circular 6.863 0 1.178 8 0.089 9 0.929 1 随机分布 Random distribution09-03
草间钻头蛛 H. graminicola 球型 Spherical 6.863 0 0.744 9 0.032 4 0.958 3 随机分布 Random distribution

表 2 表明, 在调查区域内, 桃蚜和草间钻头蛛
空间分布距离相关性比较弱 , 其变动范围均为
6.863 0~43.174 1 m; 桃蚜和草间钻头蛛的空间结构
比例分别为 0.788 8~0.983 9、0.811 6~0.980 6, 其种
群均一直表现为随机分布。
桃蚜发生初期(4 月 13 日), 种群密度较低, 块
金值(C0)和偏基台值 (C)较低 , 其空间分布属随机
分布, 种群空间相关范围为 6.863 0 m(图 1a), 即桃
蚜在直径为 6.863 0 m的区域中能相互影响。进入
4月中下旬, 桃蚜数量急速上升, 4月 29日, 其块金
值(C0)和偏基台值(C)分别增大至 4.896 3、0.182 1,
但仍属随机分布, 变程扩大至 43.174 1 m (图 1b)。
5—6 月份是桃蚜发生高峰, 块金值和偏基台值分别
为 0.479 1~1.440 6、0.010 5~0.250 0, 种群相关范围
变化较大(图 1c~f)。7—9 月气温较高, 桃蚜种群密
度骤降, 块金值和偏基台值处在较低水平, 变程为
6.863 0~43.174 1 m(图 1g~f), 桃蚜种群均呈随机分
布格局。
在桃生长早期(5月中旬前), 草间钻头蛛种群密
度较低, 块金值和偏基台值相对比较小, 而种群空
间相关范围为 6.863 0~43.174 1 m(图 2a~c)。5—7月
份是桃生长关键阶段, 桃园化学农药使用较多, 桃
910 中国生态农业学报 2015 第 23卷


http://www.ecoagri.ac.cn

图 1 化学农药胁迫下不同时间桃树桃蚜种群的半变异函数曲线
Fig. 1 Semivariograms for Myzus persicae in peach orchards under the stresses of chemical pesticides at different sampling times
图中小写字母后括号内为调查日期(月-日)。下同。 In the figures, data in the blankets following small letters are sampling date
showing in month-day. The same below.

树草间钻头蛛种群数量始终维持在较低水平, 块金
值和偏基台值也一直较低, 分别为 0.218 4~0.336 4、
0.030 4~0.075 7, 而变程为 6.863 0~43.174 1 m(图
2c~h)。8月初桃收获完毕, 桃园停止施药, 草间钻头
蛛种群数量有上升的趋势, 块金值和偏基台值都增
大, 分别为 0.498 8~0.744 9、0.032 4~0.070 5, 变程
始终维持在较小范围内(6.863 0 m), 而空间结构比
例仍大于 0.75, 其空间格局均属随机分布。
第 7期 蒋杰贤等: 基于地统计学和 GIS的化学农药胁迫下桃树桃蚜与草间钻头蛛种群空间格局 911


http://www.ecoagri.ac.cn

图 2 化学农药胁迫下不同时间桃树草间钻头蛛种群的半变异函数曲线
Fig. 2 Semivariograms for Hylyphantes graminicola in peach orchards under the stresses of chemical pesticides at different
sampling times
3 讨论
联合运用地统计学方法和GIS技术对不同时期
桃蚜及其捕食性天敌草间钻头蛛空间分布进行模拟,
从时空两个角度直观地分析害虫及其天敌的动态变
化, 既可定性判断、定量分析、可视化查询桃蚜及
其草间钻头蛛混合种群在不同时期的发生程度和变
化, 又可确定同一时期桃蚜及其草间钻头蛛在空间
上分布位置范围及其差异, 同时还提供了空间连续
的数量分布和详细的地理环境信息[23]。对桃蚜及其
912 中国生态农业学报 2015 第 23卷


http://www.ecoagri.ac.cn
草间钻头蛛空间分布的模拟并不只是了解其动态信
息 , 更重要的是揭示害虫时空动态的内在机制 [6],
为害虫的控制和预警奠定基础, 同时也为保护和利
用天敌提供依据。
草间钻头蛛是桃蚜的优势天敌, 自然状态下两
物种间具有较好的空间依赖关系和捕食被捕食关
系[89]。丁程成等[8]利用地学统计学的原理和方法分
析了李园桃蚜种群和草间钻头蛛种群的空间结构和
时间相关性, 结果表明: 绝大多数情况下桃蚜和草
间钻头蛛种群的空间分布均为聚集分布, 草间钻头
蛛种群对桃蚜种群在空间位置上具有依赖关系。害
虫的聚集格局是一种最佳的防御天敌、实现自我保
护的空间格局 , 天敌的聚集格局则是攻击捕获食
饵、提高攻击效率的最佳空间格局[8]。研究证实, 若
害虫与其天敌均存在较明显的聚集分布格局, 则害
虫与其天敌之间的跟随效应明显, 反之则无跟随效
应[2324]。有研究表明, 不施用化学农药李园(石榴园),
桃蚜(棉蚜)和草间钻头蛛种群均呈明显的空间聚集
格局, 草间钻头蛛种群对桃蚜(棉蚜)种群在空间位
置上具有追随关系[8,24]。
长期大量施用化学农药不仅杀死了害虫和自然
天敌, 也可能影响害虫与天敌之间的空间分布。本
研究表明, 4月中旬至9月上旬10次调查中桃园桃蚜
和草间钻头蛛种群的空间结构比例值均大于0.75,
桃蚜和草间钻头蛛种群始终呈随机分布, 这表明化
学农药胁迫下桃蚜和草间钻头蛛种群之间的跟随效
应不明显、草间钻头蛛对桃蚜的捕食作用不强。然
而, 前期研究表明: 桃园生草模式下桃蚜和草间钻
头蛛种群的空间结构比例值均小于0.75, 桃蚜和草
间钻头蛛种群均呈聚集分布格局, 草间钻头蛛对桃
蚜在空间位置上具有追随关系, 对其种群数量的控
制有较大作用[9]。这可能缘于桃园生草能够修复化
学农药胁迫下桃蚜和草间钻头蛛种群之间的空间相
关关系。
桃园生态系统具有多种害虫和天敌, 研究其空
间分布对科学管理桃树害虫和保护天敌具有重要意
义。本文仅报道了化学农药长期胁迫下桃树桃蚜和
草间钻头蛛种群的空间结构关系, 其他害虫与天敌
之间的空间分布有待进一步探讨。影响昆虫种群空
间格局及其动态的因子不仅仅是化学农药, 可能还
包括区域的地理位置(地势、海拔、地貌等)[7]、气象
(温度、降雨等)[25]和其他生态干扰因素(果园生草)
等 [9], 这些因子也许会直接或间接地影响昆虫种群
数量的变化。因此, 其他因子影响害虫及其天敌时
空分布及其动态变化的研究还有待进一步探讨。
参考文献
[1] 张婵婵, 张瑞芳, 张建恒, 等. 高阳县农田土壤速效养分空
间变异特征研究 [J]. 中国生态农业学报 , 2013, 21(6):
758–764
Zhang C C, Zhang R F, Zhang J H, et al. Spatial variability of
available nutrients contents in cropland soils in Gaoyang
County of Hebei Province, China[J]. Chinese Journal of
Eco-Agriculture, 2013, 21(6): 758–764
[2] 石淑芹, 曹祺文, 李正国, 等. 气候与社会经济因素对土壤
有机质影响的空间异质性分析——以黑龙江省中部地区为
例[J]. 中国生态农业学报, 2014, 22(9): 1102–1112
Shi S Q, Cao Q W, Li Z G, et al. Influence of spatial
heterogeneity of climatic and socio-economic factors on soil
organic matter — A case study of the central Heilongjiang
Province, China[J]. Chinese Journal of Eco-Agriculture, 2014,
22(9): 1102–1112
[3] Foster J R, Townsend P A, Mladenoff D J. Spatial dynamics
of a gypsy moth defoliation outbreak and dependence on
habitat characteristics[J]. Landscape Ecology, 2013, 28(7):
1307–1320
[4] Rijal J P, Brewster C C, Bergh J C. Spatial distribution of
grape root borer (Lepidoptera: Sesiidae) infestations in
Virginia vineyards and implications for sampling[J]. Environ-
mental Entomology, 2014, 43(3): 716–728
[5] Liebhold A M, Rossi R E, Kemp W P. Geostatistics and
geographic information systems in applied insect ecology[J].
Annual Review of Entomology, 1993, 38: 303–327
[6] 袁哲明, 柏连阳, 王奎武, 等. 二化螟种群密度的克力格估
值及其模拟抽样[J]. 应用生态学报, 2004, 15(7): 1166–1170
Yuan Z M, Bai L Y, Wang K W, et al. Krigle estimation and its
simulated sampling of Chilo suppressalis population
density[J]. Chinese Journal of Applied Ecology, 2004, 15(7):
1166–1170
[7] 张蓉, 冷允法, 朱猛蒙, 等. 基于地统计学和 GIS的苜蓿斑
蚜种群空间结构分析和分布模拟[J]. 应用生态学报, 2007,
18(11): 2580–2585
Zhang R, Leng Y F, Zhu M M, et al. Spatial structure analysis
and distribution simulation of Therioaphis trifolii population
based on geostatistics and GIS[J]. Chinese Journal of Applied
Ecology, 2007, 18(11): 2580–2585
[8] 丁程成, 邹运鼎, 毕守东, 等. 李园桃蚜和草间小黑蛛种群
空间格局的地学统计学研究[J]. 应用生态学报, 2005, 16(7):
1308–1312
Ding C C, Zou Y D, Bi S D, et al. Geostatistics analysis on
spatial patterns of Myzus persicae and Erigonidium
graminicola in plum orchard[J]. Chinese Journal of Applied
Ecology, 2005, 16(7): 1308–1312
[9] 万年峰, 季香云, 蒋杰贤, 等. 桃园生草复合生态系统中桃
蚜与草间小黑蛛种群的地统计学分析[J]. 应用昆虫学报 ,
2013, 50(2): 488–495
Wan N F, Ji X Y, Jiang J X, et al. Geostatistic analysis for
Myzus persicae and Erigonidium graminicola population in
the compound ecosystem of peach orchards with ground cover
getation[J]. Chinese Journal of Applied Entomology, 2013,
第 7期 蒋杰贤等: 基于地统计学和 GIS的化学农药胁迫下桃树桃蚜与草间钻头蛛种群空间格局 913


http://www.ecoagri.ac.cn
50(2): 488–495
[10] Lasmar O, Zanetti R, dos Santos A, et al. Use of geostatistics
to determine the spatial distribution and infestation rate of
leaf-cutting ant nests (Hymenoptera: Formicidae) in
eucalyptus plantations[J]. Neotropical Entomology, 2012,
41(4): 324–332
[11] Wan N F, Gu X J, Ji X Y, et al. Ecological engineering of
ground cover vegetation enhances the diversity and stability
of peach orchard canopy arthropod communities[J].
Ecological Engineering, 2014, 70: 175–182
[12] Wan N F, Ji X Y, Gu X J, et al. Ecological engineering of
ground cover vegetation promotes biocontrol services in
peach orchards[J]. Ecological Engineering, 2014, 64: 62–65
[13] Wan N F, Ji X Y, Jiang J X. Testing the enemies hypothesis in
peach orchards in two different geographic areas in eastern
China: The role of ground cover vegetation[J]. PLoS ONE,
2014, 9(6): e99850
[14] 蒋杰贤, 万年峰, 季香云, 等. 桃园生草对桃树节肢动物群
落多样性与稳定性的影响[J]. 应用生态学报, 2011, 22(9):
2303–2308
Jiang J X, Wan N F, Ji X Y, et al. Diversity and stability of
arthropod community in peach orchard under effects of
ground cover vegetation[J]. Chinese Journal of Applied
Ecology, 2011, 22(9): 2303–2308
[15] 万年峰, 季香云, 蒋杰贤, 等. 桃园生草对桃树上主要害虫
及天敌生态位的影响[J]. 生态学杂志, 2011, 30(1): 30–39
Wan N F, Ji X Y, Jiang J X, et al. Effects of ground cover on
the niches of main insect pests and their natural enemies in
peach orchard[J]. Chinese Journal of Ecology, 2011, 30(1):
30–39
[16] Schotzko D J, O’keffe L E. Effect of sample placement on the
geostatistical analysis of the spatial distribution of Lygus
hsperus (Heteroptera: Miridae) in Lentils[J]. Journal of
Economic Entomology, 1990, 83(5): 1888–1900
[17] 王政权. 地统计学及其在生态学中的应用[M]. 北京: 科学
出版社, 1999
Wang Z Q. Statistics and Its Application in Ecology[M].
Beijing: Science Press, 1999
[18] Cliff A D, Ord J K. Spatial Processes: Models and Applica-
tions[M]. London: Pion Ltd, 1981
[19] 张蓉, 朱猛蒙, 马建华, 等. 基于地理信息系统的耕地苜蓿
斑蚜种群发生的适宜生境[J]. 应用生态学报, 2009, 20(8):
1998–2004
Zhang R, Zhu M M, Ma J H, et al. Suitable habitat of
Therioaphis trifolii population on arable lands of Southern
Ningxia based on GIS[J]. Chinese Journal of Applied Ecology,
2009, 20(8): 1998–2004
[20] 王正军, 李典谟, 谢宝瑜. 基于 GIS和 GS的棉铃虫卵空间
分布与动态分析[J]. 昆虫学报, 2004, 47(1): 33–40
Wang Z J, Li D M, Xie B Y. Analysis on spatial distribution
and dynamics of Helicoverpa armigera (Hübner) eggs, based
on GIS and GS[J]. Acta Entomologica Sinica, 2004, 47(1):
33–40
[21] Cambardella C A, Moorman A T, Parkin T B, et al. Field-scale
variability of soil properties in central Iowa soils[J]. Soil
Science Society of America Journal, 1994, 58(5): 1501–1511
[22] 郭旭东, 傅伯杰, 陈利顶, 等. 河北省遵化平原土壤养分的
时空变异特征——函数与 Kriging 插值分析[J]. 地理学报,
2000, 55(5): 555–566
Guo X D, Fu B J, Chen L D, et al. The spatio-temporal
variability of soil nutrients in Zunhua Plain of Hebei Povince:
Semivariogram and Kriging analysis[J]. Acta Geographica
Sinica, 2000, 55(5): 555–566
[23] 苗进, 武予清, 郁振兴, 等. 麦红吸浆虫及其卵寄生蜂混合
种群空间格局[J]. 应用生态学报, 2011, 22(3): 779–784
Miao J, Wu Y Q, Yu Z X, et al. Spatial pattern of Sitodiplosis
mosellana and its egg parasitoids mixed population[J].
Chinese Journal of Applied Ecology, 2011, 22(3): 779–784
[24] 李磊, 邹运鼎, 毕守东, 等. 棉蚜和草间小黑蛛种群空间格
局的地统计学研究[J]. 应用生态学报, 2004, 15(6): 1043–1046
Li L, Zou Y D, Bi S D, et al. Geostatistic analysis on spatial
patterns of Aphis gossypii and Hylyphantes graminicola[J].
Chinese Journal of Applied Ecology, 2004, 15(6): 1043–1046
[25] 周强, 张润杰, 古德祥, 等. 大尺度下褐飞虱种群空间结构
初步分析[J]. 应用生态学报, 2001, 12(2): 249–252
Zhou Q, Zhang R J, Gu D X, et al. Spatial structure and
pattern of Nilaparvata lugens population in large-scale[J].
Chinese Journal of Applied Ecology, 2001, 12(2): 249–252