施肥是提高作物产量的有效手段,当施肥量超出作物养分需求量时,增产效果不但下降,而且会导致环境污染加剧、土壤微生物多样性下降、土壤肥力退化等方面的风险[1-3]。土壤肥力维持及其可持续发展在一定程度上受到土壤微生物活动的影响。细菌是土壤微生物的重要组成,其群落结构变化显著影响土壤养分循环。因此,调控土壤细菌群落结构可有效改善农田土壤肥力状况[4]。施肥措施的改变,会同时改变农田土壤理化性质和微生物群落结构。相关研究表明,长期施用化肥的农田土壤细菌数量和多样性都表现出显著降低趋势,群落结构也发生显著变化[5]。近年来,我国大力发展生态循环农业,配施有机肥是减肥增效、培肥土壤、防控农田面源污染的有效措施[6-7]。有机肥与化肥合理配施在促进作物生长的同时,还改善土壤微生物群落结构,增加有益细菌数量,对粮食高产高效和土壤质量提升起积极作用[7-8]。前人在不同种类有机肥或施用量对耕层土壤细菌群落结构变化影响方面开展了大量研究,但是在化肥减量条件下配施有机肥对不同土层土壤影响方面研究较少。本研究选用江苏里下河地区水稻高产优质创建过程中常用的两种有机肥种类(菜籽饼肥、蚯蚓粪肥)开展田间试验,采用Illumina高通量测序技术比较了不同种类有机肥替代部分化肥的施肥措施对稻季不同土层土壤环境因子、细菌群落结构特征变化的影响,以为施用有机肥替代部分化肥,实现减肥增效、培肥土壤提供科学依据。
1 材料与方法 1.1 试验材料试验地位于江苏省扬州市扬州大学农学院试验农场(32°30’ N,119°25’ E),属于亚热带季风性湿润气候,年均气温为14.5℃,年均降水量为1 000 mm,年均日照数2 300 h。试验地土壤质地为砂壤土,其有机质、碱解氮、有效磷、速效钾含量分别为23.1 g/kg、113.2 mg/kg、44.5 mg/kg、71.5 mg/kg。前茬作物为小麦,水稻供试品种为江苏沿江稻区广泛应用的武运粳24号。
试验中施用的有机肥包括菜籽饼肥(N 50.2 g/kg、P2O5 11.4 g/kg、有机质816.2 g/kg)、蚯蚓粪肥(N 26.3 g/kg、P2O5 28.3 g/kg、有机质272.7 g/kg)。菜籽饼肥购买自扬州农资公司,蚯蚓粪肥取自于扬州大学的农牧场,由培育蚯蚓消解牛粪后得到。
1.2 试验设计试验设置3个处理(表 1):单施化肥(CF)、菜籽饼氮替代化肥氮30%(RC)、蚯蚓粪氮替代化肥氮30%(WT)。CF处理为当地推荐施氮量N 270 kg/hm2,氮肥为尿素(含N 464 g/kg),分基肥、分蘖肥、促花肥施用,比例为5∶2∶3;施磷量P2O5 90 kg/hm2,磷肥为过磷酸钙(含P2O5 135 g/kg),作为基肥一次性施入;施钾量K2O 120 kg/hm2,钾肥为氯化钾(含K2O 630 g/kg),分基肥和拔节肥两次施用,比例为6∶4。有机肥随基肥一次性施入。以CF处理为参照,其他各处理中氮磷钾养分不足的部分,以尿素、过磷酸钙和氯化钾补充。各处理随机区组排列,重复3次,小区面积为30 m2,为了防止肥水串灌,在每个小区之间筑起田埂,并覆盖塑料薄膜。水稻于2019年6月6日播种,6月23日移栽,株行距为13.3 cm × 30 cm,10月31日收获。其他的田间管理措施同当地常规方法。
在水稻收获前1天,用5 cm直径的土钻采取0 ~ 5、5 ~ 10和10 ~ 15 cm土层土壤,按“S”型取15个点的土样,每个土层土样混为1个样品,剔除石砾和植物残体等杂质,装入灭菌袋暂存于冷冻保温箱中,尽快带回实验室。将每份土壤样品分为2份,一份风干研磨后,采用常规土壤农化分析方法测定土壤有机质(SOM)、pH、电导率(EC)、全氮、全磷、碱解氮、有效磷等土壤理化性质[9]。
另一份于–80℃保存,用于土壤DNA提取。土壤DNA采用Power soilTM DNA试剂盒(MoBio Laboratories, Solana Beach, CA)进行提取,每个样本3个重复。在土壤DNA提取完成后,采用16S RNA基因V3 ~ V4区通用引物338F(5′-ACTCCTACGGG AGGCAGCAGF-3′),806R(5′-GGACTACHVGGGTW TCTAATR-3′)进行PCR扩增,后续文库构建、Miseq测序均在上海美吉生物医药有限公司(Shanghai Majorbio Bio-pharm Technology Co., Ltd)的Illumina MiSeq PE300测序平台上进行[10]。
1.4 数据处理与分析采用随机区组统计多因素方差分析数据,Duncan法进行多重比较,分析处理、土层以及处理和土层交互作用对土壤理化性质和细菌多样性的影响。采用Microsoft Excel 2010和SPSS 19.0等软件进行统计分析和绘图。
在美吉平台生物云(http://www.i-sanger.com)进行高通量数据的生物信息学分析。序列拼接采用Mothur软件完成。根据序列的相似度,将有效序列聚类为OTU(97% 相似度)。采用细菌16S rRNA对比库Silva(Release132 http://www.arb-silva.de)对各OTU进行物种注释。采用基于Bray-Curtis距离矩阵的非度量多维尺度分析(non-metric multidimensional scaling,NMDS) 法分析各处理土壤细菌群落的结构差异。根据样品门水平上物种注释及丰富度绘制群落柱图,并计算细菌多样性获得相关指数。不同样本间门水平上的物种丰度差异性采用单因素方差法分析比较,获得样本间显著性差异物种。采用student’s T test方法比较各样品的细菌群落多样性。环境因子、样品细菌群落组成分布,及主要OTU之间关系的相关关系用冗余分析(RDA)比较。
2 结果与分析 2.1 减施化肥配施有机肥对稻田土壤环境因子的影响土壤有机质含量随土层加深而逐渐下降(表 2),0 ~ 5、5 ~ 10、10 ~ 15 cm土层土壤有机质含量分别平均为28.77、23.49、21.98 g/kg;与CF处理相比,RC和WT处理各土层土壤有机质含量都有显著增加,RC处理增加幅度比WT处理高。0 ~ 5、5 ~ 10、10 ~ 15 cm土层土壤pH分别平均为6.82、6.96、7.08;各处理中以RC处理土壤pH最低,CF和WT处理间差异不显著。土壤EC随着土层加深而有下降趋势,各处理中以RC处理最高,CF和WT处理间差异不显著。随着土层的加深,土壤全氮和碱解氮含量都有下降的趋势,而土壤全磷和有效磷含量变化趋势不明显,分别平均为2.16 g/kg和95.43 mg/kg。各处理中RC和WT处理的全氮、全磷、碱解氮、有效磷、速效钾含量均显著高于CF处理;0 ~ 5、5 ~ 10 cm土层土壤WT处理的碱解氮、有效磷、速效钾含量高于RC处理。
总体而言,配施不同有机肥的单一效应对收获期土壤有机质、pH、EC、全磷和速效钾含量有显著性影响;而不同土层的单一效应对pH和有效磷无显著影响,对其他各环境因子都有显著影响;有机肥和土层交互效应与土层的单一效应对土壤环境因子的影响几乎相同,可见,配施不同有机肥对土壤环境因子的影响受土层的影响作用不可忽视。
2.2 减施化肥配施有机肥对稻田土壤细菌群落结构的影响各处理土壤样品按照97% 相似性对非重复序列(不含单序列)进行OTU聚类,OTU总数是7 715个(表 3)。由Alpha多样性分析可知,与CF处理相比,RC和WT处理都能提高不同土层土壤细菌丰富度Sobs指数,以RC处理最大。Shannon指数考虑了样品的丰富度与均匀度,是对群落结构更加综合的反映。RC和WT处理都能提高不同土层土壤细菌多样性Shannon指数,其中RC处理最大。
配施不同有机肥的单一效应对收获期土壤细菌丰富度Sobs指数无显著影响,却显著影响多样性Shannon指数;不同土层的单一效应极显著影响土壤细菌丰富度Sobs指数和多样性Shannon指数。可见,土层的作用在配施有机肥对土壤细菌多样性的影响中起到重要作用。
基于Bray-Curtis距离不同处理土壤细菌非度量多维尺度分析(NMDS)表明(图 1),在0 ~ 5 cm土层,RC处理显著不同与CF和WT处理,CF和WT处理之间无显著差异(P = 0.734);在5 ~ 10 cm土层,各处理之间差异显著(P = 0.001);在10 ~15 cm土层,各处理之间的细菌群落组成呈现重合。可见,各处理对土壤细菌群落组成产生显著影响主要在5 ~ 10 cm土层。
各处理土壤细菌群落结构多样性较为丰富,变形菌门(Proteobacteria)、酸杆菌门(Acidobacteria)和绿弯菌门(Chloroflexi)是门水平上位于前3位的优势菌群,平均分别占细菌总数量的36.24%、18.30% 和13.87%(图 2)。整体而言,相对丰度大于1% 的菌门在各处理间有所差异,但是并不显著。
比较不同处理在各细菌门水平上的差异可以看出(图 3),0 ~ 5 cm土层,各处理仅在稀有菌(占比小于1%)单糖菌门(Saccharibacteria)、脱铁杆菌门(Deferribacteres)和护微菌门(Tectomicrobia)存在显著性差异;5 ~ 10 cm土层,受配施有机肥影响的菌门较多,且除放线菌门(Actinobacteria)和匿杆菌门(Latescibacteria)之外均为稀有菌门;10 ~ 15 cm土层,小基因组菌总门(Microgenomates)、降氨酸菌门(Aminicenantes)、BRC1菌门各处理间呈现显著性差异。
RDA分析结果表明,0 ~ 5 cm土层的第一、二主轴分别解释了土壤细菌群落变异的50.91% 和13.10%,两轴共解释了64.01% 的细菌群落总变异,5 ~ 10、10 ~ 15 cm土层前两轴共解释了77.40% 和58.75% 的细菌群落总变异。根据箭头长度可以看出,细菌群落分布受不同土层土壤环境因子影响作用有所差异。Mantel检验的结果表明,pH(R2 = 0.5138,P = 0.03),全磷(TP)(R2 = 0.4973,P = 0.029)是影响0 ~ 5 cm土层土壤细菌群落结构主要的土壤环境因子;EC(R2 = 0.5135,P = 0.035)、TP(R2 = 0.5721,P = 0.026)和pH(R2 = 0.4987,P = 0.031)是影响5 ~ 10 cm土层土壤细菌群落结构主要的土壤环境因子;全氮(TN)(R2 = 0.5382,P = 0.029)和pH(R2 = 0.5027,P = 0.027)是影响10 ~ 15 cm土层土壤细菌群落结构主要的土壤环境因子。
3 讨论 3.1 减施化肥配施有机肥对土壤环境因子的影响本研究表明,与单施化肥相比,减施化肥配施有机肥处理显著提高各土层土壤有机质含量,配施菜籽饼肥增幅较大(表 2),其主要原因可能是在等氮量条件下配施菜籽饼肥处理的有机质投入量较大。配施菜籽饼肥显著降低土壤pH,这可能是因为菜籽饼肥为偏酸性有机肥(pH 5.5 ~ 6.0),且有研究同样表明施用菜籽饼肥后土壤pH有下降的趋势[11]。配施蚯蚓粪肥对土壤pH则无显著影响,这可能是因为蚯蚓粪肥为中性有机肥(pH 6.8 ~ 7.1)[12]。与单施化肥相比,配施菜籽饼肥可显著增加各土层土壤电导率,而配施蚯蚓粪肥则无显著变化。其主要原因可能是配施菜籽饼肥释放出较多营养物质,从而增加溶液中的盐分含量,使土壤溶液中离子的浓度增大;而蚯蚓粪肥含大量腐殖质,阳离子交换量较大、吸附能力强,使土壤溶液中离子浓度的变化不明显[13]。配施菜籽饼肥和蚯蚓粪肥都显著增加各土层土壤全氮、全磷、碱解氮、有效磷和速效钾含量。与化肥养分释放速率相比,有机肥养分释放速率较慢,从而有利于土壤养分库的累积,特别是提高收获期土壤氮磷钾养分含量[14]。另外,有机肥的输入改善土壤理化及微生物环境,促进化学养分的固定,减少养分损失[15]。大量研究表明,配施有机肥可以显著提高土壤养分含量[14],本研究结果与之一致。另外,本研究结果表明配施蚯蚓粪肥的土壤碱解氮、有效磷、速效钾养分含量显著高于配施菜籽饼肥处理,其主要原因可能是蚯蚓粪肥更有利于土壤微生物环境的改善,有利于增强土壤团粒结构,从而促进对土壤速效养分的吸附与固持[12]。
土壤环境因子在不同土层间差异显著(表 2)。0 ~ 5 cm土层土壤有机质含量最高,其主要原因可能是,在常规耕作方式下,施用的有机肥主要集中在表层土壤中;5 ~ 10 cm土层土壤有机质含量居中,而10 ~ 15 cm土层土壤有机质含量最低,可能是该土层处于耕作下层,植株根系生长残体及其有机分泌物相对较少,因此有机质含量较低[16-17]。土壤EC、全氮、碱解氮和速效钾含量都随土层加深而呈现下降趋势,其主要原因可能是因为较深层土壤有机质含量下降,其有机质分解而释放的养分物质也随之减少[18-19];另一方面,植物根系从土壤中吸收养分,然后再通过枯枝落叶等形式将部分养分归还于表层土壤并提供了丰富的养分源,而在深层土壤造成各养分逐渐减少[20]。10 ~ 15 cm土层土壤全磷含量显著高于表层土壤,其主要原因可能是本研究中磷肥作为基肥一次性深施10 cm以下,且磷素在土壤中的移动性很小[20-21],因此收获期10 ~ 15 cm土层磷含量较高。
3.2 减施化肥配施有机肥对土壤微生物群落结构的影响施肥种类和用量都会对土壤微生物多样性产生影响。研究发现施用有机肥改变了土壤的细菌组成,其Sobs和Shannon指数都显著增加[22]。本研究同样表明,与单施化肥处理相比,配施有机肥的土壤细菌Sobs和Shannon指数都有所增加。各处理中变形菌门、酸杆菌门和绿弯菌门是土壤细菌群落中的优势菌门(图 2),这与其他研究结果基本一致[23]。不同土层深度也是影响土壤细菌多样性和群落组成差异的主要原因(表 2,图 1)。5 ~ 10 cm土层土壤细菌丰富度和多样性指数较高(表 2),其主要原因可能是该土层是植株根系生长的旺盛区域,根系残体及其有机分泌物相对较多,促进细菌繁殖[24];0 ~ 5 cm土层土壤有机质和养分含量较高,但是其细菌丰富度和多样性并非最高,可能主要是因为水稻收获前土壤长期处于厌氧环境,排水收获时表层土壤厌氧环境的消失会显著改变厌氧微生物生长环境,从而降低其丰富度和多样性。10 ~ 15 cm土层土壤细菌丰富度和多样性指数较低的主要原因可能是因为该土层处于耕作下层,根系残体及其有机分泌物等有机质含量较低[16-17]。本研究表明施用不同有机肥下不同土层土壤细菌群落结构有显著差异(图 1)。与单施化肥处理相比,5 ~ 10 cm土层土壤中放线菌门的相对丰度显著增加(图 3),其主要原因可能是:一方面,有机肥的投入,增加了有机物料及养分的释放,促进了土壤中原有放线菌的快速繁殖;另一方面,蚯蚓粪自身包含大量微生物,其中富含放线菌。这可能也是蚯蚓粪肥处理中放线菌门相对丰度较大的原因之一。配施菜籽饼肥降低了Latescibacteria菌门的相对丰度,这与前人研究结果一致[25]。整体而言,各处理之间存在显著差异的菌门主要集中在相对丰度小于1% 的菌门群落中。0 ~ 5 cm土层土壤中,与单施化肥处理相比,配施有机肥处理的单糖菌门(Saccharibacteria)相对丰度显著增加;脱铁杆菌门(Deferribacteres)、护微菌门(Tectomicrobia)相对丰度显著下降,这些与其他研究结果一致[25]。然而,目前大多研究主要集中在相对丰度较大的菌门,对这些配施有机肥后差异显著的稀有菌门的深入研究,可能有助于进一步了解配施有机肥对细菌群落结构的影响机理。
3.3 土壤细菌群落结构的影响因子有研究表明,土壤理化性质对土壤细菌群落结构具有重要作用,其中pH是影响土壤细菌群落变化的关键因子[26]。配施菜籽饼肥处理土壤pH显著小于单施化肥和配施蚯蚓粪肥处理(表 1),不同施肥处理下0 ~ 5 cm土层土壤细菌群落与pH有显著的相关关系(P < 0.05,图 5),可见pH是影响细菌群落结构的主要驱动因子之一。5 ~ 10 cm土层土壤细菌群落与EC有显著的负相关关系,其主要原因可能是土壤EC越高,养分(SOM、TN、NH4+、NO3-)累积越多,不利于土壤细菌群落结构组成的变化[27]。其他研究同样认为秸秆还田后土壤EC增大是影响细菌群落组成变化的环境因素之一[28]。10 ~ 15 cm土层土壤细菌群落与全氮有显著的负相关关系(P < 0.05,图 5),其主要原因可能也是因为土壤养分累积越多,不利于土壤细菌群落结构组成的变化。通过RDA分析发现,各土层中有机质对细菌群落结构影响并没有达到显著水平,其原因可能是配施有机肥主要通过改变土壤理化性状来调控细菌群落组成。有研究同样表明有机质是影响微生物的重要因素,但是土壤理化性状与土壤细菌群落具有显著相关性[28-29]。因此,在实施有机肥配施过程中,不仅要考虑有机肥替代化肥实现减肥增效的作用,而且要综合考虑有机肥施入土壤后对土壤理化性质的改变如何影响土壤微生物群落结构,从而改善土壤微生物环境,促进作物稳产高产。
1) 与单施化肥相比,减施化肥配施有机肥显著提高各土层土壤有机质含量和EC,配施菜籽饼肥增幅较大;显著增加各土层土壤氮磷钾养分含量,配施蚯蚓粪肥的土壤速效氮磷钾养分含量较高。
2) 减施化肥配施有机肥可提高土壤细菌Sobs指数和Shannon指数;各施肥处理中前3位优势菌门为变形菌门(Proteobacteria)、酸杆菌门(Acidobacteria)、绿弯菌门(Chloroflexi),分别占细菌数量的36.24%、18.30%、13.87%;各处理间存在显著差异的菌门主要集中在相对丰度小于1% 的稀有菌门。
3) 0 ~ 5 cm土层土壤细菌群落结构主要受pH、全磷影响;5 ~ 10 cm土层土壤细菌群落结构主要受EC、全磷、pH影响;10 ~ 15 cm土层土壤细菌群落结构主要受全氮、pH影响。
[1] |
汪洋, 杨殿林, 王丽丽, 等. 农田管理措施对土壤有机碳周转及微生物的影响[J]. 农业资源与环境学报, 2020, 37(3): 340-352 (0) |
[2] |
巨晓棠, 张翀. 论合理施氮的原则和指标[J]. 土壤学报, 2021, 58(1): 1-13 (0) |
[3] |
Wang G L, Ye Y L, Chen X P, et al. Determining the optimal nitrogen rate for summer maize in China by integrating agronomic, economic, and environmental aspects[J]. Biogeosciences, 2014, 11(11): 3031-3041 DOI:10.5194/bg-11-3031-2014 (0) |
[4] |
朱利霞. 不同调控措施对旱作农田土壤碳氮及微生物学特性的影响[D]. 杨凌: 西北农林科技大学, 2018.
(0) |
[5] |
王兴龙, 朱敏, 杨帆, 等. 配施有机肥减氮对川中丘区土壤微生物量与酶活性的影响[J]. 水土保持学报, 2017, 31(3): 271-276 (0) |
[6] |
王桂良, 张家宏, 王守红, 等. 沼液替代化肥氮对冬小麦产量、品质及生长发育的影响[J]. 农业资源与环境学报, 2018, 35(5): 467-475 (0) |
[7] |
李静. 稻麦轮作下有机无机肥料配施对作物生长及土壤微生物特性的影响研究[D]. 南京: 南京农业大学, 2015.
(0) |
[8] |
柳玲玲, 芶久兰, 何佳芳, 等. 生物有机肥对连作马铃薯及土壤生化性状的影响[J]. 土壤, 2017, 49(4): 706-711 (0) |
[9] |
鲍士旦. 土壤农化分析[M].
中国农业出版社, 北京, 1981
(0) |
[10] |
张俊华, 贾萍萍, 刘吉利, 等. 宁夏养鸡场粪污及周边土壤重金属和细菌群落特征研究[J]. 农业环境科学学报, 2020, 39(8): 1692-1705 (0) |
[11] |
祝红蕾, 储大勇, 刘华. 菜籽饼粕有机肥的应用[J]. 安徽化工, 2013, 39(5): 13-14, 24 (0) |
[12] |
张立功. 有机果园永恒的"造粪机"和产物——蚯蚓和蚯蚓粪[J]. 果农之友, 2019(1): 40-43 (0) |
[13] |
冯腾腾. 蚯蚓粪和秸秆生物炭对大棚多年连作黄瓜的调控效果[D]. 扬州: 扬州大学, 2017.
(0) |
[14] |
石纹碹, 刘世亮, 赵颖, 等. 猪粪有机肥施用对潮土速效养分含量及团聚体分布的影响[J]. 农业资源与环境学报, 2017, 34(5): 431-438 (0) |
[15] |
Jiang B S, Shen J L, Sun M H, et al. Soil phosphorus availability and rice phosphorus uptake in paddy fields under various agronomic practices[J]. Pedosphere, 2021, 31(1): 103-115 (0) |
[16] |
张大伟, 刘建, 王波, 等. 连续两年秸秆还田与不同耕作方式对直播稻田土壤理化性质的影响[J]. 江西农业学报, 2009, 21(8): 53-56 (0) |
[17] |
崔荣美, 李儒, 韩清芳, 等. 不同有机肥培肥对旱作农田土壤团聚体的影响[J]. 西北农林科技大学学报(自然科学版), 2011, 39(11): 124-132 (0) |
[18] |
李凤博, 牛永志, 高文玲, 等. 耕作方式和秸秆还田对直播稻田土壤理化性质及其产量的影响[J]. 土壤通报, 2008, 39(3): 549-552 (0) |
[19] |
程文龙, 韩上, 武际, 等. 连续秸秆还田替代钾肥对作物产量及土壤钾素平衡的影响[J]. 中国土壤与肥料, 2019(5): 72-78 (0) |
[20] |
佀国涵, 彭成林, 徐祥玉, 等. 稻虾共作模式对涝渍稻田土壤理化性状的影响[J]. 中国生态农业学报, 2017, 25(1): 61-68 (0) |
[21] |
吴明, 邵学新, 胡锋, 等. 围垦对杭州湾南岸滨海湿地土壤养分分布的影响[J]. 土壤, 2008, 40(5): 760-764 (0) |
[22] |
陆海飞, 郑金伟, 余喜初, 等. 长期无机有机肥配施对红壤性水稻土微生物群落多样性及酶活性的影响[J]. 植物营养与肥料学报, 2015, 21(3): 632-643 (0) |
[23] |
Guo J J, Liu W B, Zhu C, et al. Bacterial rather than fungal community composition is associated with microbial activities and nutrient-use efficiencies in a paddy soil with short-term organic amendments[J]. Plant and Soil, 2018, 424(1/2): 335-349 (0) |
[24] |
魏雪娇, 吴建富, 方加海. 水稻根系生长与产量形成的关系及根系生长调控途径研究进展[J]. 安徽农业科学, 2007, 35(36): 11785-11786, 11788 (0) |
[25] |
刁生鹏. 增施有机肥对旱作农田土壤环境及玉米产量的影响[D]. 呼和浩特: 内蒙古农业大学, 2019.
(0) |
[26] |
Xun W B, Huang T, Zhao J, et al. Environmental conditions rather than microbial inoculum composition determine the bacterial composition, microbial biomass and enzymatic activity of reconstructed soil microbial communities[J]. Soil Biology and Biochemistry, 2015, 90: 10-18 (0) |
[27] |
赵凤艳. 蚓粪改善设施番茄连作土壤质量和番茄品质的机理[D]. 沈阳: 沈阳农业大学, 2018.
(0) |
[28] |
刘玲玲. 秸秆还田与施氮对土壤养分、微生物及水稻产量的影响[D]. 扬州: 扬州大学, 2018.
(0) |
[29] |
张军, 刘菁, 陈长青. 有机稻鸭共作对土壤理化性状和细菌群落空间分布的影响[J]. 生态学杂志, 2020, 39(3): 822-829 (0) |