云南省耕地以坡耕地为主,极易发生水土流失,同时不合理的利用和耕作管理,导致土壤侵蚀情况进一步加剧[1],最终导致山区农业土壤退化和植物生产力下降[2]。合理的作物种植模式可以显著提高光热水肥的利用效率及土壤团聚体的稳定性[3],有效降低土壤侵蚀。间作往往是一种有效的土地利用和可持续农业实践,可以改善土壤微生态环境,促进土壤养分的高效利用[4]。已有大量的研究表明,与单作相比,间作能显著增强土壤团聚体的稳定性和提高土壤的酶活性[5-8],但这些研究主要针对整体土壤。土壤酶是生态系统中养分循环的重要动力,其活性可以作为土壤质量的指标[7, 9]。研究特定条件下,不同粒径团聚体中的酶活性具有重要意义。
土壤团聚体作为土壤结构的基本单位[10],决定着土壤养分的积累与转化,是调节土壤水、肥、气、热的关键场所,在保持土壤孔隙度、持水量、渗透性及减少土壤侵蚀等方面起着至关重要的作用[11-12]。团聚体中的土壤养分在提高土壤肥力的同时也改变了土壤团聚体的粒径分布,从而影响了团聚体的稳定性[13]。而由于团聚体中养分、通气和水分等条件的不同,造成了微生物和酶参与的养分转化特性的差异[14]。因此,有必要明确不同粒径团聚体中的酶活性和养分含量。
目前,有不少学者对不同种植模式下土壤酶活性和养分的变化做了大量的研究,但这些研究主要以整体土壤为研究对象,而针对不同粒径团聚体的研究较少。有研究表明,不同粒径团聚体的形成环境和胶结物质不同,团聚体内部的物质组成和稳定性也有差异[10],从而导致不同粒径土壤团聚体中酶活性的差异。由此,本研究以坡耕地红壤为研究对象,研究坡耕地土壤团聚体中酶活性和养分含量变化特征及对不同种植模式的响应,以期为提升坡耕地土壤肥力、保障土壤健康、减少水土流失以及促进土壤资源可持续利用提供理论依据。
1 材料与方法 1.1 供试材料田间定位试验始于2013年,期间处理和施肥均保持一致。本试验于2021年5—10月在云南农业大学后山试验小区(25°8′18″N,102°45′58″E)开展。该区的年均降雨量为960 mm,年均气温为14.2 ℃,平均海拔为1 950 m,坡向为北偏西57°,是云南省典型的坡耕地红壤水土流失区。供试土壤类型为山原红壤,基本理化性质为:有机质39.25 g/kg,全氮0.75 g/kg,有效磷9.52 mg/kg,碱解氮55.25 mg/kg,速效钾155 mg/kg,pH 6.10。供试大豆品种为滇豆7号,玉米品种为乔单6号。
1.2 试验设计试验共设置4个处理,包括大豆单作(DD)、玉米单作(MM)、大豆玉米间作(MD)、裸地(CK),每个处理设置3个重复,共12个小区。每个小区面积均为4 m×7 m,各小区间隔40 cm。玉米单作采用宽窄行种植,宽行80 cm,窄行40 cm,株距25 cm。大豆单作采用等行距种植,行距60 cm,株距25 cm;间作采用2︰2种植模式,玉米与玉米间行距40 cm,大豆与大豆间行距40 cm,玉米与大豆间行距50 cm,株距均为30 cm。各小区的播种方式均为穴播,玉米每穴播种2粒,大豆每穴播种4粒,后期玉米每穴留1株,大豆每穴留2株。每年播种前用锄头整地,轻微翻动表土,各小区一致。
单作玉米每公顷施N 250 kg、P2O5 120 kg、K2O 120 kg,其中氮肥分3次施入(基肥35%,小喇叭口期追30%,大喇叭口期追35%)。大豆每公顷施N 120 kg、P2O5 240 kg、K2O 180 kg,全作为基肥一次性施入。间作条件下玉米、大豆施肥量与各自单作相同。施肥方式为穴施,各小区一致。
1.3 样品采集与分析于2021年8月15日采集玉米抽雄期的土壤样品,每个小区采用五点取样法采集0 ~ 20 cm土层非根际土,混合成1个土壤样品,玉米和大豆单作处理分别采集玉米和玉米、大豆和大豆之间的行间土,玉米大豆间作采集玉米和大豆之间的行间土。将采集后的土壤样品放入硬质塑料盒内,尽量避免挤压,以保持原状土壤结构带回实验室风干。
采用干筛法进行土壤团聚体分级[15]:称取等量风干样品置于套筛顶部振荡,左右振幅为10 cm,频率为150次/min,振荡时间为2 min,随后称重,计算机械性团聚体各粒级质量分数。筛分的2 ~ 1、1 ~ 0.25 mm及 < 0.25 mm粒径土壤团聚体,测定其中的酶活性和养分含量。其中,脲酶和蔗糖酶活性按照《土壤酶及其研究法》[16]进行测定,酸性磷酸酶活性采用土壤酸性磷酸酶试剂盒(上海优选生物科技有限公司)检测,酶活性单位用μg/(g.24h) 表示。
土壤养分含量根据《土壤农化分析》测定[17]。
1.4 数据处理与分析试验数据采用Excel 2016进行预处理,运用SPSS17.0进行方差分析,其中,采用双因素方差分析分析不同粒级、不同种植模式对土壤养分与酶活性的影响及其交互作用,单因素方差分析分析同一粒径不同种植模式下的酶活性和养分含量及同一种植模式不同粒径下的酶活性和养分含量的差异显著性,事后多重比较采用Duncan法(P < 0.05)。采用Origin 2018作图,Canoco 5进行冗余分析(RDA),并生成RDA图。
2 结果与分析 2.1 土壤团聚体粒径分布由表 1可知,不同种植模式处理土壤均表现为2 ~ 1 mm和1 ~ 0.25 mm粒径团聚体含量最高,二者占比总团聚体的67.7% ~ 71.8%,且显著性差异仅体现在2 ~ 1 mm和 < 0.25 mm粒级。其中,MD处理显著提升了2 ~ 1 mm粒径团聚体质量分数,分别较MM、DD、CK处理提高16.2%、13.3% 和22.2%。而MD处理显著降低了 < 0.25 mm粒径团聚体质量分数,分别较MM、DD、CK处理降低23.5%、14.2%、21.4%。
由表 2可知,种植模式对土壤各养分含量均有显著影响,粒径对团聚体中养分含量的显著影响体现在有效磷和速效钾上。种植模式和粒径的交互作用对有效磷、碱解氮、速效钾含量具有显著影响。
由图 1可知,从团聚体粒径看,有机质含量除CK处理外,其他处理均表现为2 ~ 1 mm粒径团聚体中最高,且差异性仅体现在2 ~ 1 mm和1 ~ 0.25 mm粒径团聚体中。有效磷含量除DD处理外,其他处理均表现为2 ~ 1 mm粒径团聚体中最高,且差异性仅体现在2 ~ 1 mm粒径团聚体中。速效钾和碱解氮含量,各处理无明显规律,全氮含量各处理均无显著差异。
从不同种植模式看,MD处理显著提升了2 ~ 1 mm粒径团聚体中的有机质含量,分别较MM、DD、CK提升5.0%、11.9%、29.2%,且MD处理与CK处理相比,显著提高了1 ~ 0.25 mm粒径团聚体中的有机质含量。有效磷含量在2 ~ 1 mm粒径团聚体中MD处理与DD和CK处理相比,分别显著提高51.1%、57.5%(P < 0.05)。而速效钾含量表现为MD处理显著降低2 ~ 1、1 ~ 0.25、< 0.25 mm粒径含量,分别较CK处理降低5.1%、22.4%、4.2%,较DD处理降低7.5%、33.1%、11.2%,较MM处理降低3.0%、21.2%、7.1%。此外MD处理显著提高了1 ~ 0.25 mm和 < 0.25 mm粒径团聚体中中的碱解氮含量,分别较MM、DD、CK处理提高12.5%、33.3%、33.3% 和10.4%、21.8%、37.6%。全氮含量在各粒径团聚体中均表现为MD处理显著高于CK处理。
2.3 土壤团聚体中的酶活性由表 3可知,粒径和种植模式对团聚体中酶活性均有显著影响,且二者的交互作用对脲酶和酸性磷酸酶具有显著影响。
图 2显示,从团聚体粒径看,各处理的脲酶活性均表现为2 ~ 1 mm和1 ~ 0.25 mm粒径团聚体中最高。蔗糖酶活性除MM处理外,各处理均无显著性差异。酸性磷酸酶活性表现为MD处理在 < 0.25 mm粒径中最高,CK和DD处理在2 ~ 1 mm粒径团聚体中最高。
从不同种植模式(图 2)看,MD处理显著提高了2 ~ 1 mm和1 ~ 0.25 mm粒径团聚体中的脲酶含量,分别较CK、DD、MM处理提高61%、18.8%、14.5% 和65%、17.9%、13.8%,而MD与MM处理相比,显著降低了 < 0.25 mm粒径团聚体中的脲酶活性。蔗糖酶活性表现为MD处理较CK、DD处理显著提高2 ~ 1、1 ~ 0.25、< 0.25 mm粒径团聚体中的酶活性,分别提高63.8%、57.1%,42.9%、40% 和62.5%、44.4%。MD处理显著提高1 ~ 0.25 mm和 < 0.25 mm粒径团聚体中的酸性磷酸酶活性,分别较CK、DD、MM处理提高7.1%、13.2%、6.7% 和80.4%、44.3%、74.1%;而与CK、DD处理相比,MD处理显著降低了2 ~ 1 mm粒径团聚体中的磷酸酶活性,降幅为16.5%、32.7%。
2.4 环境因子与土壤团聚体中酶活性的相关性如表 4所示,团聚体中的脲酶、蔗糖酶活性与有机质、有效磷、碱解氮、全氮含量极显著正相关(P < 0.01)。冗余分析(图 3)表明,不同粒径团聚体中影响酶活性的养分各不相同。2 ~ 1 mm粒径团聚体中的脲酶和蔗糖酶活性与该粒径中的有机质、有效磷、碱解氮含量显著正相关,其中有效磷对该粒径团聚体中酶活性的影响最大。1 ~ 0.25 mm粒径团聚体中的酶活性与该粒径中的有机质、碱解氮含量显著正相关,其中有机质对该粒径团聚体中的酶活性影响最大,< 0.25 mm粒径团聚体中的酶活性与有效磷、碱解氮含量极显著正相关,且碱解氮的影响最大。
土壤团聚体是土壤结构的重要组成部分[18],是土壤环境中微生物的营养载体[19]。土壤团聚体影响养分有效性、酶活性及土壤生物地球化学过程,如生物活性、植物生长和土壤抗侵蚀能力[20]。本研究发现,不同种植模式下均显著提高2 ~ 1 mm大团聚体和1 ~ 0.25 mm小团聚体的含量,这与乌达木等[21]研究结果相似。这主要是因为铁铝氧化物是南方红壤地区微团聚体的胶结物质,70% 以上的土壤有机碳存在于粉–黏团聚体中,而大团聚体的主要胶结物质是有机质,有机质将微团聚体胶结成大团聚体[9]。并且作物种植后新鲜有机物质的输入,作物根系、微生物及他们的代谢产物增多[22],土壤中有机胶结物质的量也增多,从而促进了大粒径团聚体的形成。此外,在本研究中,与单作相比,间作显著提高了2 ~ 1 mm粒径团聚体的质量分数,这与王婷等[3]、孙涛等[5]研究结果一致。这是因为根系生物量和根长密度对团聚体的形成和稳定起着重要作用,间作模式显著增加了玉米和大豆的总根长、总根面积,固土能力增强,将微团聚体缠绕固结形成大团聚体[3]。另外,由于两根系的交互作用,促使根系分泌物增加和微生物活力增强,植物根系和真菌菌丝分泌多糖及其他有机复合物促进了土壤颗粒和微团聚体的黏结,这有利于微团聚体转化为大颗粒团聚体[22]。
土壤有机碳、氮、磷等可以提升胶结物质胶结能力,从而降低团聚体的破坏程度,促进团聚体的形成和稳定[23]。本研究发现,各处理下的有机质和有效磷含量在大粒径团聚体中最高,全氮含量在各粒径团聚体间无显著差异。葛楠楠等[24]研究指出,养分在大团聚体中含量较高,其含量与团聚体的数量和稳定性显著正相关,本研究结果与之相似,说明大粒径团聚体相较于微团聚体能够保持更高的营养水平;相反,马建辉[13]等研究发现,有机碳、全磷、全氮等含量在微团聚体中较高,这种差异可能是由于作物种类、土壤类型及耕作方式不同造成的。说明团聚体中的土壤养分在提高土壤肥力的同时也改变了土壤团聚体的粒径分布,从而影响团聚体的稳定性。此外,种植模式对土壤养分含量也有显著影响。本研究中,MD处理显著提高了2 ~ 1 mm粒径团聚体中的有机质和有效磷含量,及1 ~ 0.25 mm和 < 0.25 mm粒径团聚体中的碱解氮含量。王鹏等[6]研究发现,桔梗大葱间作提高了土壤有效氮、有效磷含量。李孝梅等[25]研究发现,不同间作模式下均显著提高大粒径团聚体中的有机碳、氮含量。孟炎奇[26]研究发现,枣树和苜蓿间作显著提高各粒级团聚体中的有机碳、氮、磷的含量。本研究结果与以上研究相似,这可能是因为作物生长过程中以根系分泌的方式向土壤中输送碳、氮等有机化合物,从而影响土壤养分含量。由于不同种植模式根系分泌物数量和种类的差异造成土壤养分含量不同,豆科与禾本科间作根系分泌物的种类和数量显著大于单作,且根系密度高,活性强,使土壤养分转化效率和土壤有效养分含量提高[6]。此外,本研究还发现,间作处理下的速效钾含量较单作低,这可能是因为土壤养分比例失调造成的。已有研究表明,土壤养分比例失调易产生离子拮抗,抑制植物对离子的吸收[6],连作造成土壤有效氮磷钾比例失调[27]。本试验地大豆、玉米单作地连续种植多年,连作造成土壤养分失调,导致钾吸收障碍,这可能是单作土壤速效钾含量高于间作的主要原因。
土壤酶是土壤生态系统中养分循环的重要驱动力,也是评价土壤肥力的指标之一[28]。不同粒径团聚体中的酶活性存在显著差异。本研究发现,同一种植模式下脲酶活性均以大粒径团聚体最高,且显著高于小团聚体,这与曹良元等[19]研究结果相似。这可能是因为大团聚体结构和孔隙分布有利于微生物生长,而土壤团聚体粒径越大,有机物相对越多,这为酶促反应提供了更多易利用的基质,从而增强了酶活性[29]。张帅等[30]研究表明,大团聚体中的酶活性显著高于小团聚体,这是因为土壤微生物代谢活动影响着酶活性,分布在大团聚体中的微生物比小团聚体丰富,所以大团聚体中的酶活性比小团聚体高[31]。这也是导致大多数酶分布在大团聚体中的重要原因。另外,由微团聚体颗粒与有机碳的胶结、根系与菌丝的缠绕作用而逐渐形成大团聚体,其中包裹着较多的颗粒有机碳,为酶提供了充足的基质[32]。本研究中还发现,不同粒径团聚体中的蔗糖酶活性差异不显著,MD处理下的酸性磷酸酶活性在小粒径团聚体中最高,说明不同种类酶在不同类型土壤、不同级别团聚体中的活性存在差异,土壤利用方式、土壤类型、物质组成及生物群落等均可能对其产生影响[32]。此外,种植模式对酶活性也有显著影响,在本研究中与单作相比,间作显著提高了2 ~ 1 mm和1 ~ 0.25 mm粒径团聚体中的脲酶活性、2 ~ 1 mm粒径团聚体中的蔗糖酶活性和 < 0.25 mm粒径团聚体中的酸性磷酸酶活性。同样地,沈雪峰等[33]研究发现,花生与甘蔗间作显著提高了脲酶和磷酸酶活性。这是因为间作中两种作物根系之间相互交叉,根系分泌物的成分和含量发生改变,为根际微生物提供纤维素、糖类等养分来改变根际微环境,促使微生物代谢活性增强,间接导致酶活性增强[6-7]。
已有大量研究表明,土壤酶活性与土壤养分具有显著的相关性[8, 28]。本研究发现,有效磷、有机质、碱解氮是影响酶活性的主要因子。代会会等[8]在大豆与番茄间作中发现,土壤脲酶、蔗糖酶、磷酸酶活性与土壤氮、磷、钾及有机质含量之间有着显著或极显著的相关性;Ma等[34]研究表明,土壤养分有效性与土壤酶活性高度相关,这可能是它们的正反馈调节作用导致的。说明土壤酶活性的分布与土壤中主要养分的分布密切相关,可以代表土壤肥力水平和作为衡量土壤肥力的生物指标。
4 结论1) 在田间试验条件下,间作模式显著提高了坡耕地红壤2 ~ 1 mm大粒径团聚体的含量及该粒径团聚体中的有机质、有效磷、碱解氮含量;显著提高了2 ~ 1 mm和1 ~ 0.25 mm粒径团聚体中的脲酶活性和2 ~ 1 mm粒径团聚体中的蔗糖酶活性及 < 0.25 mm粒径团聚体中的酸性磷酸酶活性。
2) 冗余分析表明,团聚体中的有机质、碱解氮及有效磷含量是酶活性变化的主导因子,而且有机质和有效磷含量对各粒径团聚体中酶活性的影响最大。
[1] |
张立芸, 段青松, 范茂攀, 等. 玉米和大豆根系对滇中地区坡耕地红黏土抗剪强度的影响[J]. 土壤学报, 2022, 59(6): 1527-1539 (0) |
[2] |
Xiao L, Huang Y M, Zhao J F, et al. Effects of planting structure on soil water-stable aggregates, microbial biomass and enzyme activity in a catchment of Loess Plateau Terraces, China[J]. Applied Soil Ecology, 2021, 159: 103819 DOI:10.1016/j.apsoil.2020.103819 (0) |
[3] |
王婷, 李永梅, 王自林, 等. 间作对玉米根系分泌物及团聚体稳定性的影响[J]. 水土保持学报, 2018, 32(3): 185-190 (0) |
[4] |
Fu Z D, Zhou L, Chen P, et al. Effects of maize-soybean relay intercropping on crop nutrient uptake and soil bacterial community[J]. Journal of Integrative Agriculture, 2019, 18(9): 2006-2018 DOI:10.1016/S2095-3119(18)62114-8 (0) |
[5] |
孙涛, 冯晓敏, 赵财, 等. 西北绿洲区间作模式对土壤团聚体组成及其有机碳含量的影响[J]. 农业资源与环境学报, 2021, 38(5): 874-881 DOI:10.13254/j.jare.2020.0571 (0) |
[6] |
王鹏, 祝丽香, 陈香香, 等. 桔梗与大葱间作对土壤养分、微生物区系和酶活性的影响[J]. 植物营养与肥料学报, 2018, 24(3): 668-675 (0) |
[7] |
赵雅姣, 刘晓静, 吴勇, 等. 西北半干旱区紫花苜蓿-小黑麦间作对根际土壤养分和细菌群落的影响[J]. 应用生态学报, 2020, 31(5): 1645-1652 (0) |
[8] |
代会会, 胡雪峰, 曹明阳, 等. 豆科间作对番茄产量、土壤养分及酶活性的影响[J]. 土壤学报, 2015, 52(4): 911-918 (0) |
[9] |
Gong X W, Liu C J, Li J, et al. Responses of rhizosphere soil properties, enzyme activities and microbial diversity to intercropping patterns on the Loess Plateau of China[J]. Soil and Tillage Research, 2019, 195: 104355 DOI:10.1016/j.still.2019.104355 (0) |
[10] |
刘亚龙, 王萍, 汪景宽. 土壤团聚体的形成和稳定机制: 研究进展与展望[J/OL]. 土壤学报: 1-18[2023-05-24].
(0) |
[11] |
Liu D D, Ju W L, Jin X L, et al. Associated soil aggregate nutrients and controlling factors on aggregate stability in semiarid grassland under different grazing prohibition timeframes[J]. Science of the Total Environment, 2021, 777: 146104 (0) |
[12] |
赵桂茹, 安曈昕, 欧阳铖人, 等. 青贮玉米氮投入对坡耕地土壤侵蚀及水稳性团聚体的影响[J]. 水土保持学报, 2021, 35(5): 72-79 (0) |
[13] |
马建辉, 叶旭红, 韩昌东, 等. 膜下滴灌不同灌水控制下限对设施土壤团聚体养分、酶活性和球囊霉素含量的影响[J]. 应用生态学报, 2018, 29(8): 2713-2720 (0) |
[14] |
牛文静, 李恋卿, 潘根兴, 等. 太湖地区水稻土不同粒级团聚体中酶活性对长期施肥的响应[J]. 应用生态学报, 2009, 20(9): 2181-2186 (0) |
[15] |
Bach E M, Williams R J, Hargreaves S K, et al. Greatest soil microbial diversity found in micro-habitats[J]. Soil Biology and Biochemistry, 2018, 118: 217-226 (0) |
[16] |
关松荫. 土壤酶及其研究法[M].
农业出版社, 北京, 1986
(0) |
[17] |
鲍士旦. 土壤农化分析(3版)[M].
中国农业出版社, 北京, 1981
(0) |
[18] |
宋依依, 曹阳, 段鑫盈, 等. 秸秆还田深度对土壤团聚体组成及有机碳含量的影响[J]. 土壤, 2022, 54(2): 344-350 (0) |
[19] |
曹良元, 张磊, 蒋先军, 等. 长期垄作免耕对不同大小土壤团聚体中几种氮素形态分布的影响[J]. 植物营养与肥料学报, 2009, 15(4): 824-830 (0) |
[20] |
Kumar A, Dorodnikov M, Splettstößer T, et al. Effects of maize roots on aggregate stability and enzyme activities in soil[J]. Geoderma, 2017, 306: 50-57 (0) |
[21] |
乌达木, 范茂攀, 赵吉霞, 等. 不同种植模式下坡耕地红壤团聚体有机碳矿化特征[J]. 农业环境科学学报, 2021, 40(7): 1519-1528 (0) |
[22] |
张道勇, 李会科, 郭宏, 等. 间作白三叶对苹果/白三叶复合系统土壤微生物量碳、氮及酶活性的影响[J]. 水土保持研究, 2015, 22(5): 39-45 (0) |
[23] |
赵永超, 李晓鹏, 闫一凡, 等. 激发式秸秆还田对麦季潮土团聚体中酶活性的影响[J]. 土壤, 2018, 50(3): 498-507 (0) |
[24] |
葛楠楠, 石芸, 杨宪龙, 等. 黄土高原不同土壤质地农田土壤碳、氮、磷及团聚体分布特征[J]. 应用生态学报, 2017, 28(5): 1626-1632 (0) |
[25] |
李孝梅, 李永梅, 乌达木, 等. 玉米间作大豆、萝卜对红壤不同粒径水稳性团聚体碳氮分布的影响[J]. 中国土壤与肥料, 2022(1): 104-111 (0) |
[26] |
孟炎奇. 不同间距对枣树苜蓿间作系统土壤团聚体及养分的影响[D]. 阿拉尔: 塔里木大学, 2021.
(0) |
[27] |
祝丽香, 霍学慧, 孙洪信, 等. 桔梗连作对土壤理化性状和生物学性状的影响[J]. 水土保持学报, 2013, 27(6): 177-181 (0) |
[28] |
范祎玮, 信秀丽, 钟新月, 等. 玉米豆科覆盖作物间作对潮土酶活性和真菌群落特征的影响[J]. 土壤, 2021, 53(6): 1236-1243 (0) |
[29] |
马瑞萍, 安韶山, 党廷辉, 等. 黄土高原不同植物群落土壤团聚体中有机碳和酶活性研究[J]. 土壤学报, 2014, 51(1): 104-113 (0) |
[30] |
张帅, 成宇阳, 吴行, 等. 生物炭施用下潮土团聚体微生物量碳氮和酶活性的分布特征[J]. 植物营养与肥料学报, 2021, 27(3): 369-379 (0) |
[31] |
Liao H, Zhang Y C, Wang K, et al. Complexity of bacterial and fungal network increases with soil aggregate size in an agricultural Inceptisol[J]. Applied Soil Ecology, 2020, 154: 103640 (0) |
[32] |
钟晓兰, 李江涛, 李小嘉, 等. 模拟氮沉降增加条件下土壤团聚体对酶活性的影响[J]. 生态学报, 2015, 35(5): 1422-1433 (0) |
[33] |
沈雪峰, 方越, 董朝霞, 等. 甘蔗/花生间作对土壤微生物和土壤酶活性的影响[J]. 作物杂志, 2014(5): 55-58 (0) |
[34] |
Ma Y H, Fu S L, Zhang X P, et al. Intercropping improves soil nutrient availability, soil enzyme activity and tea quantity and quality[J]. Applied Soil Ecology, 2017, 119: 171-178 (0) |