免费文献传递   相关文献

Geostatistics analysis on spatial patterns of Myzus persicae and Erigonidium graminicola in plum orchard

李园桃蚜和草间小黑蛛种群空间格局的地学统计学研究



全 文 :李园桃蚜和草间小黑蛛种群空间格局的地学
统计学研究*
丁程成  邹运鼎* *  毕守东  高彩球  刘小林  曹传旺  孟庆雷  李昌根
(安徽农业大学微生物防治省重点实验室, 合肥 230036)
摘要  利用地学统计学的原理和方法分析了李园桃蚜种群和草间小黑蛛种群的空间结构和时间相关
性.结果表明,桃蚜种群 6 月 27日和 11 月 22日此两次的半变异函数拟合模型为直线,表明此两次两种空
间分布为均匀分布或随机分布, 其余 9 次半变异函数拟合模型均为球形, 表明此 9 次李园桃蚜种群的空间
分布是聚集分布, 变程 RSD( m)为 4173903~ 347837. 草间小黑蛛空间格局 5 月21 日、5 月 31 日及 10 月
19 日、11 月 22 日 4 次半变异函数拟合模型为直线, 表明此 4 次种群的空间分布为均匀分布或随机分布,
其余 8 次半变异函数拟合模型均为球形, 表明此 8 次李园草间小黑蛛种群的空间分布是聚集分布, 变程
RSD( m)为 2 3169~ 42 00711, 表明草间小黑蛛种群对桃蚜种群在空间位置上具有追随关系,对其种群数
量的控制有较大作用.
关键词  桃蚜  草间小黑蛛  空间格局  地学统计学
文章编号  1001- 9332( 2005) 07- 1308- 05 中图分类号  Q968  文献标识码  A
Geostatistics analysis on spatial patterns of Myzus persicae and Erigonidium graminicola in plum orchard.
DING Chengcheng , ZOU Yunding , BI Shoudong, GAO Caiqiu, LIU Xiaolin, CAO Cuanw ang , MENG Qing lei, L I
Changgen ( A nhui Prov incial K ey Laboratory f or M icrobial Control , A nhui Agricultural University , H ef ei
230036, China) . Chin . J . A pp l. Ecol . , 2005, 16( 7) : 1308~ 1312.
Investig ations on the spatial construction and distribution o f Myzus persicae and Er igonidium graminicola in a
plum orchard were conducted from March 2003 to November 2003. The results indicated that t he semivariogram
of Myzus p er sicae could be described by spher ical model, except on June 27 and November 22, w hich should be
described by lined model, and that of Er igonidium graminicola could be descr ibed by spher ical model, except on
May 21, May 31, October 19 and November 22, which should be descr ibed by lined model. It could be concluded
that the amount and spatial distr ibut ion of Er igonidium graminicola w as closely related to those o f Myzus p er
sicae.
Key words  Myzus persicae, Er igonidium graminicola, Spatial pattern, Geostat istics.
* 国家农业综合开发省级科技项目 ( NFZ20011511)和安徽省教育
厅基金资助项目(2002KJ101) .
* * 通讯联系人.
2004- 07- 26收稿, 2004- 12- 06接受.
1  引   言
空间格局是种群的重要属性之一, 研究空间格
局可以揭示种群的空间结构、种群下结构的状况及
物种间的空间联系, 它是抽样技术的理论基础,过去
的研究方法都是依赖于样本频次分布或均值、方
差关 系及 其 相 应 的 判 定准 则 确 定 空 间 格
局[ 2, 4, 6, 13, 17, 19] , 用地学统计学方法研究种群空间格
局,考虑了样点的位置和彼此距离,直接测定空间结
构相关性和依赖性,克服了过去空间格局传统研究
方法的缺陷,用于研究有一定随机性和结构性的各
种变量和空间分布规律是一种较为理想的方
法[ 1, 4, 8, 10, 11, 14, 16, 19, 21, 23, 24] .本文应用地学统计学原
理和方法研究李树上桃蚜( Myz us p ersicae)及其天
敌草间小黑蛛( Erigonidium gr am inicola )的空间分
布,以揭示其空间分布型及其空间关系.
2  材料与方法
2 1 调查方法
在肥东共青林场选取样地, 南北为行.采用平行跳跃法,
株行距离为 4 m ! 4 m,共抽取样树 30 株, 分3 行每行 10 株;
每树在东、南、西、北 4个方位和上、中、下 3 层各取一代表性
枝条, 以一年生为主;调查每个枝条从梢部向内 30 cm 长度
内桃蚜和草间小黑蛛的种群数量, 同时对样地情况作详细记
录. 调查工作自 2003 年 3 月至 11 月. 具体安排: 2003 年 3
月, 1 次/月, 4~ 6 月, 2 次/月; 7 月~ 11 月, 1 次/月, 共调查
12 次, 其间果园按正常管理措施管理, 但不施化学农药.
2 2 分析方法
根据区域化变量的理论[ 3, 7, 12, 15] , 在空间上, 昆虫种群
应 用 生 态 学 报  2005 年 7 月  第 16 卷  第 7 期                               
CHINESE JOURNAL OF APPLIED ECOLOGY, Jul. 2005, 16( 7)∀1308~ 1312
数量是区域化变量.因此,可用区域化变量理论和方法进行
研究.本文通过计算桃蚜和草间小黑蛛种群的实验半变异函
数、拟合半变异函数模型、分析半变异函数结构来描述其空
间格局和空间相关关系.对于观察的数据系列 z ( x i ) , I = 1,
2, 3, #, n, 样本半变异函数 R * ( h) = 用公式计算: R * ( h )
= [ 1/ 2N ( h ) ] ∃ [ z ( x i) - z ( x i+ h) ] 2. 其中, N ( h ) 是被 h
分割的数据对( x i , x i + h ) 的对数, z ( x i ) 和 z ( x i + h ) 分别
是在点 x i和点 x i+ h 处样本的测量值, h 是分割两样点的距
离.
半变异函数有 3 个重要参数,即基台值、变程和块金值.
变程指在变异函数的值达到平衡时的空间间隔距离, 基台值
指达到平衡时的变异函数值,块金值是指变异曲线被延伸至
间隔间距为零时的截距.这 3 个参数可反映昆虫种群空间格
局或空间相关类型,给出该空间相关的范围, 几种常用的理
论模型如球型、指数型、高斯型和线型都可拟合实验半变异
函数. 球型半变异函数说明所研究的种群是聚集分布, 它的
空间结构是当样点间隔距离达到变程之前时,样点的空间依
赖性随样点间距离增大而逐渐降低. 指数型与球型模型类
似,但其基台值是渐近线,非水平直线型的半变异函数表明
种群是中等程度的聚集分布, 其空间依赖范围超过研究尺
度.如果是随机分布, 则 R * ( h )随距离无一定规律性变化,
完全随机或均匀的数据, R * ( h)呈水平直线或稍有斜率, 块
金值等于基台值, 表明抽样尺度下没有空间相关性.
3  结果与分析
31  桃蚜种群数量的半变异函数和空间格局
3月 22日调查的桃蚜种群数量为零,将从 2003
年 4~ 11月的 11次调查结果计算出的桃蚜种群半
变异函数值 R * ( h)列于表 1.对桃蚜空间种群结构
进行分析,结果显示, 6 月 27日和 11 月 22日此两
次的半变异函数拟合模型为直线, 表明此两次种群
的空间分布为均匀分布或随机分布, 其余 9次半变
异函数拟合模型均为球形, 表明此 9 次李园桃蚜种
群的空间分布是聚集分布, 如表 2所示.
  桃蚜不同时期 11 个半变异函数拟合模型具有
较高的决定系数, 其变化范围 05607~ 09363, 表
明实际半变异函数值与理论曲线的拟合度极高; 4
月 6日、6月 8日、8 月 27 日和 9 月 21 日的模型及
其拟合曲线如图 1A 所示. 不同时期桃蚜种群半变
异函数的变程 RSD( m)即为样点间的空间依赖性距
离,其变动范围为 4173903~ 3478371, 如表 2所
示.
表 1  不同时期桃蚜的半变异函数
Table 1 Observed semivariogram R * ( h) for Myzus persicae at different time
时间
Date
距离 h Distances( m)
4 8 12 16 20 24 28 32 36
46 426E02 147E02 161E02 179E02 002 227 E02 263 E02 003125 0
427 9610682 4713898 409962 4357767 496772 56480 6328832 2880703 7144654
521 2994372 1040838 1124306 1263625 141542 159625 1862132 2211438 03846154
531 1877545 6492002 7119961 7875318 8828758 1003268 1161678 1379491 0
68 01808511 02647059 01935484 01964286 006 227E02 526E02 00625 385E02
627 2179787 2629412 1643548 1291071 686 1125 1044737 1128125 0
719 1856383 2117647 1108065 1905357 119 4704545 4657895 575 06538461
827 837234 1229412 982258 1216071 712 05681818 05526316 003125 0
921 2202128 2367647 1564516 1232143 164 1522727 09210526 059375 01923077
1019 851E02 01176471 00483871 536E02 002 682E02 789E02 0 0
1122 9114893 1212941 4312903 46875 4478 6265909 7231579 0125 01538462
表 2  不同时期桃蚜的理论半变异函数模型 R* ( h)、参数值、决定系数 R 2
Table 2 Models, parameters, R 2 value and spatial patterns for Myzus persicae
时间
Date
模型
Models
a b c d 变程
RSD( m)
决定系数
R 2
空间格局
Spat ical pat tern
46 球形1) - 100E05 00007 - 00124 00802 3478371 08519 聚集3)
427 球形 - 20032 11947 - 21441 160091 2607902 09363 聚集
521 球形 - 08246 49868 - 87726 56538 2735173 08516 聚集
531 球形 - 51538 31172 - 54863 353964 2734689 0852 聚集
68 球形 400E05 - 00021 00265 0129 8257886 08392 聚集
627 直线2) - 06078 25186 07309 R/ U
719 球形 00009 - 00579 04363 18231 4173903 0809 聚集
827 球形 00025 - 01609 24823 08746 1008354 09227 聚集
921 球形 - 00001 00071 - 01607 28769 1872253 08896 聚集
1019 球形 - 100E05 00007 - 00137 01481 3270231 05607 聚集
1122 直线 - 27222 10805 05749 R/ U
注: R * / U 指随机分布或均匀分布 R ( h) = ah3+ bh2+ ch+ d. 1) Sphere; 2)L inear; 3) Aggregat ion.下同 Th e same below .
13097 期          丁程成等:李园桃蚜和草间小黑蛛种群空间格局的地学统计学研究           
图 1  桃蚜( A)及草间小黑蛛( B)半变异函数曲线
Fig. 1 Model of semivariogram of Myz us persicae ( A) and Er igonidium graminicola ( B) .
a) 46; b) 68; c) 827; d)921.
表 3  不同时期的草间小黑蛛半变异函数
Table 3 Observed semivariogram R * ( h) for Erigonidium graminicola at different time
时间
Date
距离 h Distances( m)
4 8 12 16 20 24 28 32 36
322 319E02 441E02 323E02 179E02 002 227E02 0 0 0
46 319E02 441E02 323E02 179E02 002 227E02 0 0 0
427 04893617 07794118 0516129 03035714 026 02954545 03157895 009375 769E02
521 156383 2445588 4935484 8428572 85 7318182 8684211 20625 03461539
531 1934043 2122059 1145161 1498214 752 9340909 6473684 38125 01923077
68 1105319 1554412 816129 9714286 922 5227273 2789474 16875 769E02
627 1119149 1517647 6370968 5410714 57 9522727 8342105 38125 3846154
719 45 6882353 3225806 2160714 148 2886364 1578947 10625 06538461
827 9808511 8426471 3274194 3375 314 3795455 3236842 303125 3076923
921 4340425 5 4451613 3857143 474 35 3263158 234375 1961538
1019 1051064 8205882 4951613 8267858 528 8113636 5394737 6 7769231
1122 745E02 735E02 00483871 536E02 004 227E02 263E02 0 0
32  草间小黑蛛种群数量半变异函数和空间格局
不同时期的草间小黑蛛半变异函数值列于表
3,以此为实际值对理想半变异函数最优模型进行拟
合,如图 1B所示.草间小黑蛛空间种群结构进行分
析,结果显示, 5月 21日、5月 31日及 10月 19日、
11月 22日 4次的半变异函数拟合模型为直线, 表
明此 4次种群的空间分布为均匀分布或随机分布,
其余 8 次半变异函数拟合模型均为球形, 表明此 8
次李园草间小黑蛛种群的空间分布是聚集分布(表
4) .
1310                    应  用  生  态  学  报                   16 卷
表 4  不同时期草间小黑蛛的理论半变异函数模型 R * ( h)、参数值、决定系数 R 2
Table 4 Models, parameters, R 2 value and spatial patterns for Erigonidium graminicola at di fferent time
时间
Date
模型
Models
a b c d 变程
RSD( m)
决定系数
R 2
空间格局
Spat ial pattern
322 球形1) 300E06 - 00002 00022 00293 6430366 08531 聚集3)
46 球形 300E06 - 00002 00022 00293 6430366 08531 聚集
427 球形 200E05 - 00013 00057 05948 23161 07389 聚集
521 直线2) - 05082 19093 0582 R/U
531 直线 - 06017 22516 08879 R/U
68 球形 700E04 - 00486 05462 1061 6544767 08752 聚集
627 球形 - 110e03 00724 - 15587 18457 2494001 05349 聚集
719 球形 200E05 00017 - 02487 64815 4200711 07003 聚集
827 球形 - 800E04 0061 - 15257 1551 2862581 09026 聚集
921 球形 500E05 - 00068 01137 41167 9317905 08882 聚集
1019 直线 - 00702 85702 01708 R/U
1122 直线 - 00025 00871 09462 R/U
  草间小黑蛛不同时期 12 个半变异函数拟合模
型具有较高的决定系数, 其变化范围 01708 ~
09462, 4月 6日、6月 8日、8月 27 日和 9月 21日
的模型及其拟合曲线如图 1B所示. 不同时期草间
小黑蛛种群半变异函数的变程 RSD( m)即为样点间
的空间依赖性距离, 其变动范围为 23169 ~
4200711(表 4) .
4  讨   论
  通过对本研究 12次调查结果分析,李园节肢动
物群落中天敌共有 1890只,其中 100只以上的天敌
有瓢虫科 124只,蟹蛛科 176只,肖蛸科 101只, 微
蛛科 1236只, 微蛛科主要是草间小黑蛛,占调查的
天敌总数的 6539%. 李园节肢动物群落中害虫共
有6912 只, 其中 100只以上的害虫有刺蛾科 108
只,毒蛾科 208只, 叶甲科 441 只, 蚜科 5703只, 蚜
科主要是桃蚜, 占调查的害虫总数的 8251%. 王昌
贵等[ 20]、周琼等[ 22]、蒋金妹等[ 5]研究表明, 草间小
黑蛛可捕食松大蚜、棉蚜、槐蚜和桃蚜等同翅目蚜科
害虫.
利用昆虫种群空间格局的经典方法难以区别不
同空间格局的差别[ 8, 9] , 地学统计学方法分析昆虫
空间格局可以避免上述问题.
对不同时期桃蚜及草间小黑蛛种群的地学统计
学分析可以看出,除 6月 27日和 11月 22日外, 其
余时间桃蚜种群的空间格局均呈聚集分布, 其半变
异函数皆为球型曲线, 变程在 4173908~ 347873
cm 之间, 这是因为一定空间和时间中环境因子或自
然资源的供应最适宜其生活, 同时社会行为的结果
有利于其生存和繁殖, 个体聚集在一起又保持相当
距离,以维持彼此食物和环境资源的平衡. 除 5 月
21日、5月 31日、10月 19日和 11月 22日外, 草间
小黑蛛种群亦呈聚集分布,其理论半变异函数曲线
图与桃蚜种群理论半变异曲线图形状相同,变程在
23169~ 420071 cm 之间,不同时期草间小黑蛛种
群变程的变化趋势与桃蚜种群变程的变化趋势一
致.通过地统计学分析,调查的 11次数据,草间小黑
蛛有 7次为聚集格局, 桃蚜有 9次为聚集格局, 草间
小黑蛛和桃蚜同时都为聚集格局的有 6次,占桃蚜
聚集格局总次数的 67%, 占草间小黑蛛聚集格局总
次数的 86% , 害虫聚集格局是一种最佳的防御天
敌,自我生存保护的空间格局,天敌的聚集格局是攻
击捕获食饵,提高攻击效率的最佳空间格局.这充分
说明草间小黑蛛种群对桃蚜种群在空间位置上具有
较强的追随关系, 是确定草间小黑蛛是桃蚜的优势
种天敌的证据之一.实践证明,地学统计学能为探索
昆虫种群的形成和为何形成其空间分布以及种间竞
争机制提供强有力的工具.
参考文献
1  Bi SD(毕守东) , Zou YD(邹运鼎) , Geng JG(耿继光) , et al .
2000. Geostat iscal analysis on spatial pat terns of Aphi s gossypi i and
Prop ylaca japonica . Chin J A ppl Ecol (应用生态学报 ) , 11 ( 3 ) :
421~ 424( in Chinese)
2  Hassell MP, May RM. 1974. Aggregat ion of predators and insect
parasites and it s effect on stabilit y. J A nim Ecol , 43: 567~ 594
3 H ou JR ( 侯景儒 ) , Huang JX ( 黄竞先 ) . 1985. Theory and
Method of Geostat ist ics. Beijing:Ch ina Geology Press. ( in Chinese)
4  Hurlbert SH. 1990. Spat ial distribut ion of the montane unicorn.
Oikos , 58: 257~ 271
5  Jiang JM (蒋金妹) , Luo MH ( 罗梅浩 ) , Guo XR (郭线茹 ) .
2004. The seasonal characters and structure of the arthropod com
munity in aubergin e and capsicum field in Zhengzhou suburbs. J
Henan Agri c Univ (河南农业大学学报 ) , 38( 1) : 54~ 58( in Ch i
nese)
6  Jumars PA, Thist le D, Jones ML. 1977. Detect ing twodimensional
spat ial st ructure in biological data. Oecologia, 28: 109~ 123
7  Lecoust re L, Farget te D, et al . 1989. Analysis and mapping of the
spat ial spread of Af rican cassava mosaic virus using geostat ist ics
and kriging technique. Phytop athology , 79: 913~ 920
8  Liebhold AM, Rossi RE, Kemp WP. 1993. Geostatist ics an d geo
graphic in form at ion systems in applied insect ecology. Annu Rev
En t , 38: 303~ 327
13117 期          丁程成等:李园桃蚜和草间小黑蛛种群空间格局的地学统计学研究           
9  Liebhold A M, Zhang X, et al . 1991. Geostat ist ical analysis of gyp
symoth ( Lepidoptera: Lymant riidae) egg mass populat ions. E nv i
ron Ent , 20: 1407~ 1417
10  Li YC (李友常) , Xia NB (夏乃斌) , Tu QH (屠泉洪 ) , et al .
1997. A geostatist icae analysis on spatial patterns of A noplophora
g labr iper nis in popears. Acta Ecol S in (生态学报) , 17( 4) : 393~
401( in Chinese)
11  Li L(李  磊) , Zou YD(邹运鼎 ) , Bi SD(毕守东 ) , et al . 2004.
Geostatist ic analysis on spatial patterns of A phis gosspi i and Er igo
nidium graminicola . Chin J Ap pl E col (应用生态学报) , 15 ( 6) :
1043~ 1046 ( in Chinese)
12  Rossl RE, Mulla DJ , Journel AG, et al . 1992. Geostatist ical tools for
modeling an d interpret ing ecological spatial dependence. Ecol
Monogr , 62( 2) : 277~ 314
13  Saw yer AJ. 1989. Inconstancy of Tabor% s sb: Simulated sam pling
w ith dif ferent quadrat sizes and spat ial dist ribut ion . Res Pop ul
Ecol , 31: 11~ 24
14  S chotzko DJ, O% keeffe LE . 1989. Geostatical descript ion of the spa
t ial dist ribut ion of Lygus hesper us( Heteroptera: M iridae) in lentils.
J Econ Ent , 82( 5) : 1277~ 1288
15  Schotzko DJ, O% keffe LE . 1990. Ef fect of sample placement on the
geostat ist ical analysis of the spatial distribut ion of L ygus hesper us
( Heteroptera: M iridae) in Lent ils. J Econ Ent , 83( 5) : 1888~ 1900
16  Shi GS(石根生) , Li DM( 李典谟 ) . 1997. S patial dist ribution of
massonpine caterpillar pupae and its parasit ic natural enemy guilds
for the overw intering generat ion in the different pine forest s. Acta
Ecol Sin (生态学报) , 17( 4) : 386~ 392 ( in Chinese)
17  Southwood TRE. 1978. Ecological Methods. London: C hapman and
Hall.
18  Taylor LR. 1984. Assessing and interpret ing the spat ial distribu
tions of insect populat ion. A nnu Rev E nt , 29: 321~ 357
19  Vieia SR, Hat if ield JL, et al . 1989. Geostat ical theory and applica
tion to variabilit y of some agronomical propert ies. Hilgardia, 51: 1
~ 75
20  Wang CG( 王昌贵 ) , Zhen HM (郑华美) , Li SY( 李素英) , et
al . 2003. The study of reproduct ion habit s and appet ite observat ion
of Erigonidium g raminicola. Shand ong For Sci Technol (山东林
业科技) , 5: 39( in Chinese)
21  Williams IIIL, Schotzko DJ, McCaff rey JP. 1992. Geostat ist ical de
script ion of the spat ial distribut ion of L imoniu s calif or nicus
( Coleoptera: Elateridae) w irew orms in the northw estern United
States, w ith comments on sampling. En vi ron Ent , 21 ( 5) : 983~
995
22  Zhou Q(周  琼) , Liang GW(梁广文) , Zeng L(曾  玲) . 2004.
Populat ion dynamics of aphids and their natural enemies on crucifer
vegetables in Shenzhen Suburbs. Chin J Biol Contr ol (中国生物防
治) , 20( 1) : 70~ 71( in Chinese)
23  Zou YD(邹运鼎) , Bi SD( 毕守东) , Wang XS(王祥胜) , et al .
2001. Geostat iscal analysis on spat ial paterns of Macrosiph um g ra
narium and A phid ius . Chin J Ap pl Ecol (应用生态学报) , 12( 6) :
887~ 891( in Chinese)
24  Zou YD(邹运鼎) , Bi SD(毕守东) , Zhou XZ(周夏芝) , et al .
2002. Geostat ist ical analysis on spatial patterns of Erythroneu ra su
dra and Er igonidium graminicola . Chin J App l Ecol (应用生态
学报) , 13( 12) : 1645~ 1648( in Chinese)
作者简介  丁程成, 女, 1980 年生, 硕士, 博士研究生. 主要
从事昆虫生态学研究, 发表论文 10 篇.
欢迎订阅 2006年&应用生态学报∋
&应用生态学报∋( 1990年创刊)是经国家科委批准、科学出版社出版的国内外发行的综合性学术刊物.
本刊宗旨是坚持理论联系实际的办刊方向,结合科研、教学、生产实际,报导生态科学诸领域在应用基础研究
方面具有创新的研究成果,交流基础研究和应用研究的最新信息,促进生态学研究为国民经济建设服务.
本刊专门登载有关应用生态学(主要包括森林生态学、农业生态学、草地牧业生态学、渔业生态学、自然
资源生态学、全球生态学、污染生态学、化学生态学、生态工程学和恢复生态学等)的综合性论文、创造性研究
报告和研究简报等.
本刊读者对象主要是从事生态学、地学、林学、农学和环境科学研究、教学、生产的科技工作者,有关专业
的大学生及经济管理和决策部门的工作者.
本刊与数十家相关学报级期刊建立了长期交换关系, &中国科学引文索引∋、&中国生物学文摘∋、美国
&生物学文摘∋( BA)、美国&化学文摘∋( CA)、英国&生态学文摘∋( EA)、日本&科学技术文献速报∋( CBST)和俄
罗斯&文摘杂志 ∋(   )等十几种权威检索刊物均收录本刊的论文摘要(中英文) ,并被认定为&中国核心期刊
(遴选)数据库∋和&中国科学引文数据库∋来源期刊.本刊的整体质量与水平已达到新的高度, 1992年荣获全
国优秀科技期刊三等奖和中国科学院优秀期刊二等奖, 1996年荣获中国科学院优秀期刊三等奖. 2000年荣
获中国科学院优秀期刊二等奖. 2001年入选中国期刊方阵双效期刊.
本刊为月刊, A4开本, 192页,每月 18日出版,期定价 45. 00元,全国各地邮政局(所)均可订阅, 邮发代
号 8- 98.错过订期也可直接向本刊编辑部邮购,个人订阅优惠 30%.地址: 110016 辽宁省沈阳市文化 72号
&应用生态学报∋编辑部. 电话: ( 024) 83970393, Email: cjae@ iae. ac. cn
1312                    应  用  生  态  学  报                   16 卷