免费文献传递   相关文献

Degradation-detoxification behavior of methamidophos in phaiozem and burozem rhizosphere

黑土和棕壤中甲胺磷的根际降解脱毒模拟研究



全 文 :黑土和棕壤中甲胺磷的根际降解脱毒模拟研究 3
于 颖1 ,2  周启星1 3 3
(1 中国科学院沈阳应用生态研究所 陆地生态过程重点实验室 ,沈阳 110016 ;
2 大连海事大学环境科学与工程学院 ,大连 116026)
【摘要】 利用塑料根际盒研究了甲胺磷在黑土和棕壤大豆根际和非根际环境中的降解脱毒行为. 结果表
明 ,甲胺磷虽是急性毒性较高的农药 ,但在土壤环境中能很快降解 ,并且同等条件下 ,甲胺磷在黑土中的残
留量普遍低于棕壤. 在无大豆种植情况下 (对照处理) ,培养试验第 2 天 ,棕壤甲胺磷残留量约为 33 % ,黑
土只有 26 %. 在大豆根际圈中 ,甲胺磷的降解明显加快 ,尤其是在黑土中. 第 9 天 ,根际盒中层黑土和棕壤
的农药残留分别比无植物对照低 8715 %和 7610 %. 甲胺磷的土壤降解过程符合一级动力学方程 ,降解半
衰期为 2 d 左右.
关键词  黑土  棕壤  甲胺磷  根际圈  降解脱毒
文章编号  1001 - 9332 (2005) 09 - 1761 - 04  中图分类号  X13115  文献标识码  A
Degradation2detoxif ication behavior of methamidophos in phaiozem and burozem rhizosphere. YU Ying1 ,2 ,
ZHOU Qixing1 (1 Key L aboratory of Terrest rial Ecological Process , Institute of A pplied Ecology , Chinese Acade2
my of Sciences , S henyang 110016 , China ; 2 College of Envi ronmental Science and Engineering , Dalian M ar2
itime U niversity , Dalian 116026 , China) . 2Chin. J . A ppl . Ecol . ,2005 ,16 (9) :1761~1764.
With plastic rhizobox system ,this paper studied the degradation behavior of methamidophos in the rhizosphere
and non2rhizosphere of phaiozem and burozem planted with soybean. The results suggested that methamidophos
could be rapidly decomposed in soil environment despite its higher acute toxicity , and its residual amount in
phaiozem was generally below the level in burozem under same condition. In the absence of soybean (control) ,the
residual amount of methamidophos at the 2nd day of incubation was about 33 % in burozem ,whereas only about
26 % in phaiozem. An accelerated degradation of methamidophos in soybean rhizospheric soil was observed ,espe2
cially in phaiozem. At the 9th day of incubation , the insecticide residue in the middle rhizobox phaiozem and
burozem with soybean was decreased by 8715 % and 7610 % ,respectively ,compared with that of the control. The
degradation process of methamidophos in soil environment followed the first2order equation ,and its half2life was
about 2 days.
Key words  Phaiozem , Burozem , Methamisophos , Rhizosphere , Degradation2detoxification.3 国家杰出青年科学基金项目 (20225722) 、国家自然科学基金重点
项目 ( 20337010 ) 和国家重点基础研究发展规划资助项目
(2004CB418503) .3 3 通讯联系人. E2mail :zhouqixing2003 @yahoo. com
2004 - 10 - 12 收稿 ,2005 - 03 - 26 接受.
1  引   言
甲胺磷 (Methamidophos ,化学名 O , S - 二甲基
氨基硫代磷酸酯) 是一种广谱、高效的有机磷杀虫
剂 ,我国年需求量保持在万吨以上[6 ] . 由于长期使
用 ,我国许多地区土壤中存在着一定的残留污染 ,其
含量范围达 1317~3318 mg·kg - 1 [21 ,22 ] . 同时 ,甲胺
磷易溶于水、急性毒性较高 ,对土壤生态系统生命组
分的毒害效应和对地表水体、地下水污染的威胁不
容忽视[3 ,7 ,16 ,18 ] . 目前 ,尽管有些省份已明令禁止销
售和使用甲胺磷 ,但由于甲胺磷价格便宜、药效高 ,
在东北地区仍然普遍应用[19 ] .
根际圈以植物根系为中心 ,聚集了大量的生命
物质及其分泌物 ,构成了极为独特的“生态单元”. 鉴
于其独特的物理、化学和生物学活性[14 ] ,根际圈在
污染土壤生态脱毒与修复中的作用日益受到关
注[23 ,24 ] . 有研究表明 ,农药在根际中的分解速率明
显加快. 对硫磷在水稻根际中降解 2216 %时 ,非根
际仅降解 515 %[10 ] . 涕灭威在玉米、绿豆和豇豆根
际环境中的降解速率常数分别是其对照的 115~
210 倍[13 ] . 当处理浓度为 110 mg·kg - 1时 ,丁草胺在
小麦非根际土与根际土的半衰期分别为 1810 d 和
717 d ,降解速率常数分别为 010385 和 010902 [20 ] .
尽管对农药的根际降解动态已有一些文献报道 ,但
是与现有偏重重金属根际效应的研究[8 ,15 ]相比 ,这
方面的资料非常缺乏 ,仍需从根际环境的自身脱毒
角度进行深入分析 ,以便明确农药根际降解脱毒过
程的主要影响因子.
应 用 生 态 学 报  2005 年 9 月  第 16 卷  第 9 期                               
CHIN ESE JOURNAL OF APPL IED ECOLO GY , Sep. 2005 ,16 (9)∶1761~1764
大豆 ( Glyci ne m ax )是我国东北的主栽作物 ,对
于该地区农业的可持续发展具有重要意义. 但是 ,大
豆根际是否具有较强的降解甲胺磷的生态行为则没
有相关报道. 为此 ,本研究以大豆为供试植物 ,研究
甲胺磷在大豆根际中的降解脱毒 ,为利用根际环境
降低有机污染物毒性 ,减少污染物的食物链传递和
富集提供理论依据.
2  材料与方法
211  供试材料
21111 土壤  供试黑土和棕壤是东北地区的典型土壤 ,分别
采自中国科学院海伦生态实验站 (47°26′N ,126°38′E) 和沈
阳十里河生态试验站 (41°31′N ,123°241′E) 未被污染的休耕
地.黑土和棕壤的 p H 分别为 6148 和 6109 ,有机质分别为
3198 %和 2131 % ,CEC 分别为 32192 和 19123 cmol·kg - 1 , <
01001 mm 黏粒组成分别为 29102 %和 17110 %. 土壤风干后
过 210 mm 筛 ,按 25 %含水量均匀拌入 40 %甲胺磷乳油 ,使
土壤含纯甲胺磷的浓度为 20 mg·kg - 1 (风干土) .
21112 作物  东北大豆 ,品种为铁丰 30 号 ,购自沈阳农业大
学.
21113 根际盒装置  由 3 个大小不一的塑料方盒相叠组成
(图 1) ,规格分别为 :上盒 910 cm ×910 cm ×317 cm ;中盒
913 cm ×913 cm ×315 cm ;下盒 1010 cm ×1010 cm ×610
cm. 上盒底部镂空 ,用 400 目尼龙纱网与中盒分隔 ,内装相
当于 200 g 干土重的黑土和棕壤. 中盒和下盒内装同样土
壤 ,分别相当于干土重 100 g 和 200 g. 为避免中盒与下盒内
土壤水分状况相差悬殊 ,中盒底部有两条平行的细透水缝.
图 1  根际盒装置示意图
Fig. 1 Rhizobox constructed profile.
212  研究方法
21211 根际培养试验  催芽后 6 粒大豆种子等距离播种于
上盒内 ,每盒加入 25 ml 蒸馏水并称重记录 ,作为其后补充
水分的参照. 根际盒放入恒温光照生化培养箱内 25 ℃培养 ,
并在间歇光照条件 (约 3 000 lx 光照 14 h ,黑暗 10 h)下培养
2、5、9、14、20 和 27 d 后取样分析土壤甲胺磷残留量. 试验共
有对照和种植大豆两种处理 ,每处理 3 次重复. 对照除上盒
不种大豆外 ,其余条件完全相同.
21212 甲胺磷土壤残留测定  称相当于干土重 10 g 的鲜土
于具塞三角瓶中 ,加 4 g NaCl 和 60 ml 蒸馏水 ,振荡 1 h 后将
土壤溶液转移到玻璃离心管中 ,离心 ,过滤入另一三角瓶 ,再
向具塞三角瓶中加入 20 ml 蒸馏水清洗土壤残渣 ,过滤 ,三
角瓶中合并两次提取液 ,从中吸取 40 ml 置于 250 ml 分液漏
斗内 ,加 10 g 无水硫酸铵 (饱和) ,顺次加入 20 ml 乙酸乙酯、
20 ml 乙酸乙酯/ 甲醇混合液 (乙酸乙酯与甲醇的比例为 5∶
1) 、15 ml 乙酸乙酯萃取 ,静置分层后 ,过盛装无水硫酸钠漏
斗 (无水硫酸钠经 400 ℃马福炉 3 h 烘过) ,合并有机相 ,乙
酸乙酯定容 50 ml 并在旋转浓缩仪上浓缩至 1 ml ,利用气相
色谱 ( HP6890)测定[11 ] . 气相色谱条件 :毛细管柱长 30 m ,内
径 0132 mm ,内装 HP2210 , 190912613 (50 % Trifluoropropyl
ME Siloxame) ;温度 :柱温 130 ℃,进样口 220 ℃,氮磷检测
器 280 ℃;气体流速 :氮气 30 ml·min - 1 、氢气 3 ml·min - 1 、空
气 60 ml·min - 1 ;进样量 1μl ,保留时间为 213 min. 样品回收
率在 70 %以上.
3  结果与讨论
311  甲胺磷的降解脱毒动态
由图 2 可知 ,甲胺磷是急性毒性较高的杀虫剂 ,
土壤降解十分迅速. 在无大豆种植情况下 (对照处
理) ,培养试验开始第 2 天 ,棕壤甲胺磷残留量约为
33 %、黑土甲胺磷残留量约 26 % ;至第 14 天 ,根际
盒上层黑土已无甲胺磷检出.
种植大豆对甲胺磷降解影响显著. 在根际盒上
层土壤中 ,第 3 次取样 (第 9 天) 是值得注意的转折
点 ,在此之前大豆生长对甲胺磷降解影响不大. 在第
图 2  甲胺磷在黑土 (A)和棕壤 (B)中的残留动态 (25 ℃)
Fig. 2 Methamidophos residue in phaiozem ( A) and burozem soils (B)
(25 ℃) .
a)根际盒上层土壤 Upper soil in plastic rhizobox ;b) 根际盒中层土壤
Middle soil in plastic rhizobox ;c) 根际盒下层土壤 Lower soil in plastic
rhizobox. 1)对照 Control ;2)种植大豆 Planted soybean.
2671 应  用  生  态  学  报                   16 卷
9 天时 ,甲胺磷在 2 种土壤中的残留水平明显低于
对照处理 ,黑土低 014 mg·kg - 1 、棕壤低 015 mg·
kg - 1 (图 2a) .
根际盒中层土壤是根际效应最明显的层次. 随
着大豆根量的增加 ,甲胺磷的根际降解效应逐渐明
显. 尽管培养试验的第 2 天是大豆膨胀萌芽时期 ,从
理论上讲 ,此时对照与种植大豆处理之间甲胺磷残
留量比较接近 ,而根际盒中层土壤甲胺磷的残留量
仍表现出一定的差异 ,在种植大豆的黑土和棕壤根
际盒中层 ,甲胺磷残留量分别比其对照处理低 015
mg·kg - 1 (916 %) 和 018 mg·kg - 1 (1112 %) . 第 5 天
时 ,根际盒中层黑土和棕壤甲胺磷残留量分别比对
照处理低 2413 %和 1413 % ;至第 9 天 ,种植大豆黑
土根际盒中层甲胺磷残留量不足 011 mg·kg - 1 、棕
壤不足 013 mg ·kg - 1 , 分别较对照低 8116 %和
7610 %(图 2b) .
根际盒下层是基本不受大豆根系影响的土壤 ,
相当于非根际土 ,但在本试验中 ,甲胺磷在根际盒下
层土壤中的降解仍表现出一定的类“根际效应”. 由
图 2c 可见 ,第 2 天和第 5 天 ,种植大豆的根际盒土
壤中 ,下层甲胺磷残留量与对照处理相当 ,至第 9 天
方显现出明显差距 ,种植大豆黑土和棕壤根际盒下
层甲胺磷残留量分别比对照低 7818 %和 4719 %. 在
本试验末期 (第 27 天) ,无大豆种植土壤根际盒下层
仍有微量甲胺磷检出 ,而有大豆种植情况下 ,根际盒
下层黑土和棕壤分别至第 14 天、第 20 天已检测不
出甲胺磷.
在根际盒设计中 ,尽管大豆生长于上层土壤 ,但
由于大豆播种较深 ,根系向下伸展的过程中距离尼
龙网越来越近 ,培养中后期已在尼龙网处形成少量
根垫. 因此 ,与尼龙网紧密接触的中层土壤成为根际
效应最为明显的层次. 下层土壤虽然距离大豆根系
较远 ,但甲胺磷在该层土壤中的降解仍被加速 ,可能
是甲胺磷是水溶性农药 ,易于随每日补充的水分 ,沿
根际盒中盒底部的两条细透水缝向根际盒下层淋
溶 ,由于甲胺磷在对照中层土壤中的残留量高于同
层种植大豆处理 ,因此淋溶进入对照下层土壤的农
药量较多 ,造成种植大豆的下层土壤似有根际效应
的假象.
312  甲胺磷的降解脱毒模式
一般情况下 ,农药在土壤中的降解可以用一级
动力学方程来描述[2 ] :
C = C0 e - kt (1)
式中 , C 为某时刻农药浓度 , C0 为初始浓度 , k 为降
解速率常数 , t 为时间. 农药降解半衰期可由式 t015
= 01693/ k 算得. 在本试验中 ,甲胺磷在各土壤样品
中的降解符合一级动力学方程 ,相关系数 r 达到了
显著水平 (表 1) .
表 1  一级动力学方程拟合表达式、显著性分析及甲胺磷半衰期
Table 1 Fitting representation , test of goodness of f it of f irst2order e2
quation and half2life of methamidophos
根际盒
Rhizobox
土壤样品
Soil sample
表达式
Representation
r 半衰期 t015
Half2life (d)
上 层 Ⅰ C = 10146e - 01332 t 01996 3 3 (n = 3) 1) 211
Upper Ⅱ C = 20138e - 01569 t 01981 3 (n = 3) 112
Ⅲ C = 12165e - 01307 t 01994 3 3 3 (n = 4) 213
Ⅳ C = 14154e - 01376 t 01999 3 3 3 (n = 4) 118
中 层 Ⅰ C = 9194e - 01313 t 01995 3 3 3 (n = 4) 212
Middle Ⅱ C = 17161e - 01566 t 01992 3 3 (n = 3) 112
Ⅲ C = 12101e - 01296 t 01988 3 3 (n = 4) 213
Ⅳ C = 10168e - 01357 t 01988 3 3 (n = 4) 119
下 层 Ⅰ C = 8131e - 01273 t 01998 3 3 3 (n = 6) 215
Lower Ⅱ C = 17136e - 0151 t 01995 3 3 (n = 3) 114
Ⅲ C = 10154e - 01257 t 01998 3 3 3 (n = 6) 217
Ⅳ C = 9128e - 01264 t 01988 3 3 (n = 4) 216
1)可检测出的时间点数 The detectable time. Ⅰ. 黑土2对照 Phaiozem2control ; Ⅱ. 黑土2大豆
Phaiozem2soybean ; Ⅲ. 棕壤2对照Burozem2control ; Ⅳ. 棕壤2大豆Burozem2soybean. 3 P < 0102 ;3 3 P < 0101 ; 3 3 3 P < 010011
  由于甲胺磷降解较快 ,一些样品至第 4 次取样
(第 14 天)已经无检出 ,因此各一级动力学方程的时
间区段不一致. 有些方程的 r 值为显著 ,但计算值
与实测值误差较大 , C0 值发生较大的漂移. 上中下
层黑土2对照和上中下层棕壤2对照处理的计算 C0
值与施加初始浓度 20 mg·kg - 1相去甚远 ,最为接近
的处理是种植大豆的黑土上层土壤 ,计算 C0 值为
29138 mg·kg - 1 (表 1) . 产生 C0 值漂移的方法学原
理在于 :一级动力学方程对数转换后 ,用最小二乘原
理进行目标函数最小化计算时 ,采用无约束最小化 ,
导致估算 C0 值漂移的现象普遍存在. 有时也与降
解条件、降解机制的变化及试验误差有关. 但不能因
为存在较大的 C0 值漂移就认定农药降解不符合一
级动力学反应[17 ] .
通过计算一级动力学方程中的降解常数 k 和
半衰期 t015可以发现甲胺磷在各土壤样品中降解存
在一定的规律性. 甲胺磷在黑土中的降解速率高于
棕壤 ;种植大豆使甲胺磷的降解速率明显加快 ,尤其
是对于黑土. 在无大豆种植情况下 (对照处理) ,甲胺
磷在上中下层黑土和棕壤中降解半衰期的平均值分
别为 213 d 和 214 d ;在有大豆种植情况下 ,甲胺磷
在上中下层黑土和棕壤中降解半衰期的平均值分别
为 113 d 和 211 d.
有报道表明 ,甲胺磷在未灭菌的土壤中比在灭
菌的土壤中降解快 ,在富含有机质的土壤中比在较
贫瘠的土壤中降解快 ,证明微生物对于土壤中的甲
36719 期             于  颖等 :黑土和棕壤中甲胺磷的根际降解脱毒模拟研究            
胺磷降解起重要作用[1 ] . 事实上 ,甲胺磷降解菌在
自然界分布广泛 ,且类型多样. 曲霉 ( Aspergill us) 、
青霉 ( Penicilli um ) 、酵母 ( S accharom yces) 等真菌以
及假单胞菌 ( Pseudomonas) 、芽孢杆菌 ( B acill us) 等
被发现可降解甲胺磷[4 ] . 从甲胺磷的酯键结构考
虑 ,推测微生物降解甲胺磷主要为水解作用[1 ,5 ,12 ] .
本试验中 ,大豆根际对甲胺磷降解有明显促进效应 ,
说明微生物因素对该降解过程扮演重要角色. 此外 ,
一些土壤组分可通过酯键吸附2催化模式水解有机
磷农药. 土壤中存在着较多的氧化物 (如 O3 、H2O2
以及有机质等) ,可在体系中产生较多的 - OH ,使有
机磷农药快速彻底的水解[9 ] . 甲胺磷在黑土中的消
失速率较棕壤快 ,表明有机质的作用不可低估 ,也进
一步佐证土壤水解机制可能是影响甲胺磷降解的重
要因子之一.
4  结   论
在大豆根际环境中 ,甲胺磷的降解速率显著加
快 ,并且在黑土中表现出更为明显的趋势. 同时 ,甲
胺磷在黑土中的降解快于棕壤 ,主要在于土壤有机
质可能通过酯键吸附2催化水解机制促进甲胺磷降
解脱毒. 甲胺磷在土壤中的降解行为可以很好地用
一级动力学方程来描述 ,相关系数 r 达到显著水
平.从一级动力学方程式中计算甲胺磷土壤降解半
衰期为 2 d 左右.
参考文献
1  Cao Z2F (曹志方) , Wang Y2S (王银善) . 1996. Biodegradation of
methamidophos by microbes. A dv Envi ron Sci (环境科学进展) ,4
(6) :32~35 (in Chinese)
2  Cockerham L G ,Shane BS. 1994. Basic Environmental Toxicology.
Boca Raton :CRC Press.
3  de Castro VL ,Chiorato SH ,Pinto NF. 2000. Relevance of develop2
mental testing of exposure to methamidophos during gestation to its
toxicology evaluation. Toxicol Letters ,118 :93~102
4  Fang X2H (方晓航) , Qiu R2L (仇荣亮) . 2003. The degradation
transformation of organophosphorus pesticides in soil environment .
Envi ron Sci Technol (环境科学与技术) ,26 (2) :57~59 (in Chi2
nese)
5  Heydorn LN , Wong CY ,Srinivas R , et al . 2003. The isobaric ions
CH3O2P = O·+ and CH3O2P2NH +2 and their neutral counterparts :
A tandem mass spectrometry and CBS2QB3 computational study.
Int J M ass S pect rom ,225 :11~23
6  Hua X2M (华小梅) ,Shan Z2J (单正军) . 1996. The production and
application of pesticides and factor analysis of their pollution in envi2
ronment in China. A dv Envi ron Sci (环境科学进展) ,4 (2) :33~
45 (in Chinese)
7  Liang J2D(梁继东) ,Zhou Q2X(周启星) . 2003. Single and binary2
combined toxicity of methamidophos ,acetochlor and Cu on earth2
worm Esisenia f oeli de. Chin J A ppl Ecol (应用生态学报) , 14
(4) :593~596 (in Chinese)
8  Mc Grath SP , Shen ZG , Zhao FJ . 1997. Heavy metal uptake and
chemical changes in the rhizosphere of Thlaspi caerulescens and
Thlaspi ochroleucum grown in contaminated soils. Plant Soil ,
188 :153~159.
9  Meng Z2Q (孟紫强) . 2000. Environmental Toxicology. Beijing :
China Environmental Science Press. (in Chinese)
10  Reddy BR ,Sethanathan N. 1983. Mineralization of parathion in the
rice rhizosphere. A ppl Envi ron Microbial ,45 :826~829
11  Shi L2L (石利利) ,Lin Y2S (林玉锁) , Xu Y2G(徐亦钢) , et al .
2001. Detection of pesticide residues in soil and crops in south of
Jiangsu Province. A gro2Envi ron Proc (农业环境保护) ,20 (3) :158
~159 (in Chinese)
12  Singh A K , White T , Spassova D , et al . 1998. Physicochemical ,
molecular2orbital and electronic properties of acephate and
methamidophos. Com p Biochem Physiol ,119C:107~117
13  Sun HW ,Xu J , Yang SH , et al . 2004. Plant uptake of aldicarb from
contaminated soil and its enhanced degradation in the rhizosphere.
Chemosphere ,54 :569~574
14  Wang M2E(王美娥) ,Zhou Q2X(周启星) ,Zhang L2H (张利华) .
2003. Chemical behavior and ecological effects of pollutants acting
on root2soil interface. Chin J A ppl Ecol (应用生态学报) ,14 (11) :
2067~2071 (in Chinese)
15  Wang ZW , Shan XQ , Zhang SZ. 2002. Comparison between frac2
tionation and bioavailability of trace elements in rhizosphere and
bulk soils. Chemosphere ,46 :1163~1171
16  Wu YL ,Chang PA ,Li M , et al . 2003. Effect of tri2o2cresyl phos2
phate and methamidophos on 45Ca uptake by brain synaptosomes in
hens. Pest Biochem Physiol ,77 :18~23
17  Xiang F(向 锋) . 1998. Discussion on some issues on parameter of
first2order chemical kinetics reaction equation. R ural Eco2Envi ron
(农村生态环境) ,14 (3) :20~25 (in Chinese)
18  Yen J H ,Lin KH , Wang YS. 2000. Potential of the insecticides a2
cephate and methamidophos to contaminate groundwater. Ecotoxi2
col Envi ron S af ety ,45 :79~86
19  Yu Y ,Zhou QX. 2003. Effects of methamidophos on sorption2des2
orption behavior of copper in soils. B ull Envi ron Contam Toxicol ,
71 :979~987
20  Yu YL ,Chen YX ,Luo YM , et al . 2003. Rapid degradation of bu2
tachlor in wheat rhizosphere soil. Chemosphere ,50 :771~774
21  Zhang HW , Zhang QR , Zhou QX , et al . 2003. Binary2joint effects
of acetochlor ,methamidophos ,and copper on soil microbial popula2
tion. B ull Envi ron Contam Toxicol ,71 :746~754
22  Zhou Q2X (周启星) , Huang G2H (黄国宏) . 2001. Environmental
Biogeochemistry and Global Environmental Changes. Beijing : Sci2
ence Press. (in Chinese)
23  Zhou Q2X(周启星) ,Song Y2F (宋玉芳) . 2001. Technological im2
plications of phytoremediation and its application in environment
protection. J S af ety Envi ron (安全与环境学报) ,1 (3) :48~53 (in
Chinese)
24  Zhou Q2X (周启星) , Song Y2F (宋玉芳) . 2004. Remediation of
Contaminated Soils : Principles and Methods. Beijing :Science Press.
(in Chinese)
作者简介  于  颖 ,女 ,1975 年生 ,博士. 主要从事污染物环
境行为与生态毒理学研究 ,参编专著 1 部 ,发表论文 10 篇.
E2mail :yuying @newmail. dlmu. edu. cn 或 yyuhelen @hotmail.
com
4671 应  用  生  态  学  报                   16 卷