油茶-花生间作土壤酶活性与养分的关系

莫 晶 ，闫文德 ，刘曙光 ，3，吴小红

（1. 中南林业科技大学，湖南 长沙 410004；2. 南方林业生态应用技术国家工程实验室，湖南 长沙 410004；3. 美国地质勘探局 地球资源观测与科学研究中心，南达科他州 苏福尔斯 SD57198）

油茶-花生间作土壤酶活性与养分的关系

莫 晶1，2，闫文德1，2，刘曙光1，2，3，吴小红1，2

（1. 中南林业科技大学，湖南 长沙 410004；2. 南方林业生态应用技术国家工程实验室，湖南 长沙 410004；3. 美国地质勘探局 地球资源观测与科学研究中心，南达科他州 苏福尔斯 SD57198）

本研究比较了丘陵地区油茶单作、花生单作及油茶-花生间作三种种植模式下五种土壤酶（蔗糖酶、蛋白酶、脲酶、酸性磷酸酶、过氧化氢酶）活性的差异，并采用相关分析、通径分析对土壤酶活性和土壤养分的相关性进行了研究，旨在探讨油茶种植模式对土壤酶活性的影响，以期为油茶生产中培肥和改良土壤提供参考。结果表明，不同种植模式对土壤酶活性影响明显。与油茶单作相比，油茶-花生间作显著提高0～10 cm土层土壤蔗糖酶、蛋白酶、脲酶及过氧化氢酶的活性，以及10～20 cm土层蔗糖酶活性；与花生单作相比，油茶-花生间作显著提高了0～10 cm土层土壤蛋白酶活性。相关性分析表明土壤酶活性与土壤养分的相关性较高，特别是脲酶和酸性磷酸酶。通径分析结果表明，土壤速效磷和速效钾对蔗糖酶产生显著的直接影响，有机质和速效磷是影响蛋白酶的主导因子，土壤全氮是影响脲酶的重要因素，土壤全氮和速效磷是直接影响酸性磷酸酶的重要因子。综合看来，与油茶单作相比，油茶-花生间作模式可以改善土壤质量，提高土壤酶活性，对提高土壤肥力有着重要意义。

油茶；间作；土壤酶活性；土壤养分

油茶Camellia oleifera又名茶籽树，为山茶科Theaceae山茶属Camellia常绿小乔木或灌木，是我国特有的木本食用油料树种，主要分布在我国南方地区[1]。由于油茶种植前几年树冠扩展较慢，林地内有一定的空间，因此，在油茶新造林前期可以利用林地内的空间种植作物来提高土地资源的利用率，增加经济收益。研究表明，在油茶林下种植农作物不仅可以提高土地的经济效益，还有助于提高土壤肥力和土壤质量，改善油茶生长环境，促进油茶生长结实[2]。

油茶林下间作的农作物种类繁多，包括花生、大豆、绿豆、红薯、豇豆、旱稻、绿肥、油菜、辣椒等[3-4]，现有的文献主要集中于研究不同油茶农作物间作模式对油茶生长生理特性和土壤理化性质的影响，而有关油茶与农作物间作对土壤酶活性的影响以及土壤酶活性与土壤养分关系的研究还少有报道。土壤养分是土壤中含有植物生长代谢所需要的营养物质，是反映土壤肥力的基础指标[5]。土壤酶是具有催化活性的蛋白质，具有生物催化剂之称，它能促进有机体中的化学反应，是土壤生物学活性的总体现[6]。土壤酶既参与包括土壤生物化学过程在内的自然界物质循环，又是植物营养元素的活性库[7-8]，在生态系统物质循环和能量流动中发挥着重要的作用，是土壤肥力评价的一个重要指标。因此，本研究通过分析油茶单作、花生单作、油茶-花生间作三种种植模式的土壤酶活性差异，并采用相关分析和通径分析研究土壤酶活性与土壤养分之间的关系，探讨不同种植模式对土壤酶活性的影响，旨在为油茶复合经营土壤管理提供科学的理论依据。

1 材料与方法

1.1 试验样地概况

试验地位于湖南省长沙市林业科学院试验林场（113°01′30″E，28°06′40″N），试验地属低丘岗地地貌，海拔50～108 m，属亚热带季风湿润气候。年平均气温为17.5℃，极端最高温度42.7℃，极端最低温度-10.4℃。年平均降水量1 378 mm，年均相对湿度81%，年日照1 814.8 h，全年无霜期275 d。

1.2 试验设计

供试油茶林为2010年所造，油茶株行距为4 m×3 m，种植品种为湘林4号。试验设置油茶单作、花生单作和油茶-花生间作3种种植模式，每种种植模式设置3个重复。油茶单作模式设置：选择长势较一致的连片油茶3行，每行3株为一个小区，重复3次，油茶株间不间作花生；花生单作模式设置：在油茶林旁边空地按照穴距行距0.25 m×0.1 m种植花生，每个小区48 m2，重复3次；油茶-花生间作模式设置：选择长势较一致的连片油茶3行，每行3株为一个小区，重复3次，于2016年4月在油茶株间间作花生。花生与油茶植株的间距为1 m，花生穴距0.25 m，株距0.1 m，花生品种为市场购买的适合湖南地区种植的品种。三种种植模式均使用“旺亿农”有机肥进行一次性施入，锄草、打药等其他田间管理措施均一致。

1.3 土壤样品采集及测定方法

1.3.1 土壤样品采集

土壤样品于2016年7月采集。分别在油茶单作、花生单作及油茶-花生间作模式内以蛇形布置5个采样点，采用土钻分层采集0～10 cm和10～20 cm土壤样品后，将各样点相同层次土壤混合，取样时避开施肥区。现场采用四分法弃去多余的样品，保留1 kg土壤样品装于自封袋中，写好标签，带回实验室。将取回的土样摊放在实验室，除去土壤中可见的动、植物残体及石块等，让其自然风干。风干后土样分为两份，一份过1 mm筛混匀后置于常温保存，用于土壤酶活性的测定。另一份土样过0.25 mm筛后供土壤有机质、全氮、全磷、全钾、速效氮、速效磷及速效钾的测定。

1.3.2 测定方法

土壤酶活性测定：土壤酶活性均采用比色法进行测定。蔗糖酶活性的测定采用3，5-二硝基水杨酸比色法[9]，单位以1 g土壤恒温培养（37℃）24小时内分解蔗糖释放出的葡萄糖毫克数表示；蛋白酶活性的测定采用茚三酮比色法[9]，单位以1 g土壤在37℃恒温培养箱中培养24 h后的氨基氮的毫克数表示；酸性磷酸酶活性的测定采用磷酸苯二钠比色法[9]，单位以1 g土壤释放的酚的毫克数表示；脲酶活性的测定采用苯酚钠-次氯酸钠比色法[9]，单位以1 g土壤在37℃恒温培养箱中培养24 h释放出的铵态氮毫克数表示；过氧化氢酶活性的采用高锰酸钾滴定法[9]，单位以20 min内1 g土壤分解的过氧化氢的毫克数表示。

土壤养分测定：土壤有机质采用水合热重铬酸钾氧化法；全氮用半微量凯氏定氮法；全磷用王水酸溶钼锑抗比色法；全钾用王水酸溶火焰光度计法；速效氮用蒸馏法；速效磷用双酸浸提分光光度比色法；速效钾用醋酸铵浸提火焰光度计法。

1.4 数据分析

采用 R（3.3.1）软件对实验数据进行统计分析。不同种植模式间的差异显著性采用单因素方差分析（One-way Anova）检验（P＜0.05表示存在显著性差异），不同土层之间土壤酶活性差异显著性采用Tukey’s HSD 检验（P＜0.05表示存在显著性差异）；相关性分析采用Pearson统计方法进行分析；采用R软件中的sem包进行通径分析。

2 结果与分析

2.1 不同种植模式对土壤酶活性的影响

土壤蔗糖酶、蛋白酶、脲酶、酸性磷酸酶和过氧化氢酶广泛存在于土壤中，在土壤有机化合物和动植物残体等的分解转化方面发挥着重要的作用。由图1可知，油茶单作、花生单作和油茶-花生间作3种种植模式5种土壤酶活性表现出不同程度差异。对于0～10 cm土壤，油茶-花生间作模式中土壤蔗糖酶、蛋白酶、脲酶及过氧化氢酶活性分别比油茶单作模式提高了92.00%、85.68%、71.07%、16.68%，且差异达显著性差异水平（P＜0.05），但油茶-花生间作模式中土壤酸性磷酸酶活性比油茶单作模式降低了百分比14.49%。对于10～20 cm土壤，与油茶单作相比，油茶-花生间作模式提高了蔗糖酶、蛋白酶和过氧化氢酶的活性，但降低了脲酶和酸性磷酸酶的活性。与花生单作模式相比，在0～10 cm土层，油茶-花生间作模式中土壤蛋白酶、脲酶、酸性磷酸酶及过氧化氢酶活性分别提高了77.81%、37.25%、25.08%、5.07%，且对于蛋白酶达显著性差异水平（P＜0.05），但蔗糖酶的活性降低了16.93%。在10～20 cm土层，油茶-花生间作模式中的蔗糖酶、蛋白酶、酸性磷酸酶活性高于花生单作模式，但脲酶和过氧化氢酶的活性有所降低。此外，分析结果还表明，在3种种植模式中，土壤酶活性均随着土壤深度的增加而呈现出垂直递减的趋势，除花生单作模式的蔗糖酶外，其他模式的土壤酶活性在0～10 cm及10～20 cm并未表现出显著差异。

图1 不同种植模式下土壤的酶活性Fig.1 Soil enzymes activities under different cropping patterns

2.2 土壤酶活性与土壤养分之间的相关分析

2.2.1 相关分析

将不同种植模式土壤酶活性与土壤养分含量进行Pearson相关分析得到表1。可以看出，土壤蔗糖酶活性与速效氮存在显著负相关关系（P＜0.05），与其他土壤养分相关性弱；蛋白酶活性与有机质和速效磷存在极显著正相关（P＜0.01）和显著正相关关系（P＜0.05），与全氮、全磷、全钾、速效氮和速效钾相关性弱；脲酶活性与有机质和全氮存在极显著正相关关系（P＜0.01），与速效磷和速效钾呈显著正相关（P＜0.05），与其他土壤养分指标相关性不显著；酸性磷酸酶活性与有机质、全氮、速效磷及速效钾均存在极显著正相关关系（P＜0.01），与全磷存在显著正相关关系（P＜0.05），但与全钾和速效氮相关性弱；过氧化氢酶活性与全氮和速效钾存在显著正相关关系（P＜0.05），与其他土壤养分相关性弱。这说明土壤酶活性和土壤养分含量之间存在相互影响、相互作用的关系。

表1 土壤酶活性与土壤养分的相关性†Table 1 Coefficient between soil enzyme activities and the soil nutrient contents

2.2.2 通径分析

通径分析是标准化的多元线性分析，可以通过对自变量与因变量之间表面直接相关性的分解，来研究自变量对因变量的直接重要性和间接重要性。经数据标准化处理后的土壤蔗糖酶活性（Y1）、蛋白酶活性（Y2）、脲酶活性（Y3）、酸性磷酸酶活性（Y4）、过氧化氢酶活性（Y5）与土壤养分因子的多元线性回归方程为：

式中：X1、X2、X3、X4、X5、X6、X7分别是标准化的土壤有机质、全氮、全磷、全钾、速效氮、速效磷、速效钾。方程（1）～（5）中的系数为直接通径系数，直接通径系数乘以各土壤养分之间的相关系数即得到表2中的间接通径系数。直接通径系数表示各土壤养分因子对酶活性的直接影响，间接通径系数则是土壤养分因子通过其他养分因子间接对酶活性产生的影响。

由表2可以看出，土壤速效磷、速效钾对土壤蔗糖酶活性的直接通径系数较大，分别为-0.766 9、0.624 0，两者对土壤蔗糖酶活性分别有较强烈的负效应和正效应，土壤全磷、全钾对土壤蔗糖酶活性的直接影响力虽然很小（直接通径系数分别为-0.062 7、0.118 9），但二者通过速效磷对蔗糖酶活性的间接通径系数达到了-0.636 5和-0.521 5。可见，土壤全磷、全钾对蔗糖酶的影响主要体现在间接影响上。土壤有机质及速效磷对土壤蛋白酶的直接通径系数较大，分别为0.601 4和0.637 3，表明二者对蛋白酶活性有较强的直接效应，土壤全氮对蛋白酶的影响主要体现在通过土壤有机质的间接影响上。土壤全氮对脲酶的影响主要是直接影响（直接通径系数为0.444 0），其他养分因子对脲酶的直接、间接通径系数都较小，对脲酶活性的影响也较小。土壤全氮和速效磷对土壤酸性磷酸酶活性的通径系数较大，分别为0.827 7、0.631 7，全磷对土壤酸性磷酸酶活性的影响主要表现在通过速效磷的间接影响上，说明全氮和速效磷是影响酸性磷酸酶的主要因素。相对于其他四种酶，土壤养分因子对过氧化氢酶的直接和间接通径系数都相对较小，表明土壤养分因子对过氧化氢酶活性的影响较小。

3 讨 论

3.1 不同种植模式下土壤酶活性差异分析

土壤酶直接参与了土壤营养元素的有效化过程，在一定程度上反映了土壤养分转化的动态，能够在较短的时间内反映出土壤状况的变化[10]。蔗糖酶来自植物根系和微生物，能催化蔗糖水解成葡萄糖和蔗糖，对土壤的碳、氮循环起到重要的作用[11]。蛋白酶可以水解蛋白质为肽，最终形成氨基酸，加速土壤的氮素循环。脲酶是氮素循环的关键酶类，能酶促尿素的水解，水解形成的NH3是植物的氮源之一，其活性可反映土壤的供氮能力和水平[12]。磷酸酶是土壤中最活跃的酶类之一，其能酶促有机磷化合物的水解，可为植物生长提供所需的速效磷。过氧化氢酶能促进过氧化氢对各种化合物的氧化，减轻土壤中过氧化氢的毒害作用，亦可反映土壤中总的呼吸强度。

表2 土壤养分对5种酶活性的通径系数†Table 2 Path coefficients between the soil nutrients affecting activities of the five enzymes

合理的间作有利于油茶生长、改善土壤的理化性质[13]、提高土壤酶活性及土壤微生物的数量和多样性[3]。贾树海[14]及田亚玲[15]等人的研究发现农林间作可提高0～20 cm土壤蔗糖酶、脲酶、蛋白酶、酸性磷酸酶及过氧化氢酶的活性。本研究发现，除0～20 cm土层脲酶外，油茶-花生间作模式土壤蔗糖酶、蛋白酶、脲酶及过氧化氢酶的活性均显著高于油茶单作模式（P＜0.05），与前人的研究结果一致[14-15]。油茶-花生间作模式脲酶活性显著增加可能与豆科植物花生有关，花生根系根瘤菌具有较好的固氮作用，提高了土壤氮素的含量，而脲酶活性与土壤氮素转化密切相关。此外，花生根系的穿插作用会改变土壤的孔隙度和土壤的通气状况，根系在生长过程中向土壤释放的分泌物也会影响土壤微环境，从而对蔗糖酶、蛋白酶与过氧化氢酶活性产生影响。张萌萌[16]等人对桑树/苜蓿间作的研究表明间作提高了土壤酸性磷酸酶的活性，而本研究中油茶-花生间作对土壤酸性磷酸酶活性的影响不显著，这可能与作物种类以及作物生长阶段有关[17]。研究显示，酸性磷酸酶的活性受植物腐化分解的影响较大，而本研究采样时期正值花生旺盛生长阶段，土壤中可供分解的花生植物残体少，因而对土壤酸性磷酸酶活性的影响较小。

3.2 土壤酶活性与土壤养分的相关性分析

土壤酶活性与土壤养分之间关系的研究一直是研究的热点问题，早期许多学者运用简单的相关分析对两者的关系做了相应的研究，运用相关系数说明两者之间的关系。本研究对土壤酶活性与土壤有机质、全氮、全磷等养分因子的相关关系分析结果表明，土壤酶活性与养分有着密切的关系，土壤的肥力水平很大程度上受土壤酶活性的影响。但相关分析结果只反映了土壤养分与土壤酶活性关系的大小，并不能揭示土壤养分对土壤酶活性的直接和间接影响程度[18]。因此进一步采用通径分析全面的考察了土壤养分因子与酶活性之间的直接关系与间接关系，使二者之间的影响过程更为明晰。

通径分析的研究结果表明，不同土壤养分对酶活性的影响有所差异。土壤速效磷对蛋白酶和酸性磷酸酶存在着直接的正效应，对蔗糖酶却产生了负效应。这是因为速效磷是土壤中可被植物、动物、微生物吸收的磷组分，而蛋白酶和酸性磷酸酶酶促反应的底物主要就是来源于土壤植物的根系、残体，土壤动物及其残骸和微生物的分泌物。全氮是影响脲酶和酸性磷酸酶的重要因素，对酶活性存在较大的直接正效应，因为脲酶是一种促进氮有机物的水解酶，与土壤氮素供应和转化情况密切相关。前人的研究也表明脲酶及酸性磷酸酶活性的变化与土壤氮素状况有关[19-20]。袁霞等对八角林地土壤酶活性的研究表明土壤有机质与蛋白酶活性显著相关，与本文研究结果一致。土壤速效钾是影响蔗糖酶的主要因素，可能是由于本地区的土壤中钾素的亏缺有关[21]，吴雪[22]等对阿拉尔垦区土壤理化因子与酶活性关系的研究也表明土壤速效钾与蔗糖酶存在显著相关关系。

本文对3种模式的土壤酶活性及其与养分关系进行了分析，得出油茶-花生间作模式对提高土壤酶活性有明显的效果。然而不同作物及不同生长阶段具有不同的生物学特性及生境条件，今后的研究应增加农作物的种类，并增加研究时间，以期筛选出效果更佳的间作模式，更好的为油茶生产中培肥和改良土壤提供参考依据。

4 结 论

本文得到的主要结论如下：

（1）油茶-花生间作模式0～10 cm土层的土壤蔗糖酶、蛋白酶、脲酶及过氧化氢酶的活性显著（P＜0.05）高于油茶单作模式，在10～20 cm土层，油茶-花生间作模式的蔗糖酶的活性显著高于油茶单作模式（P＜0.05），与花生单作相比，油茶-花生间作显著（P＜0.05）提高了0～10 cm土层土壤蛋白酶的活性。

（2）三种种植模式的土壤多数酶与土壤有机质、全氮、速效磷、速效钾呈极显著或显著正相关，但与全钾没有显著相关性。

（3）土壤速效磷和速效钾对蔗糖酶产生显著的直接影响，有机质、速效磷是影响蛋白酶的主导因子，土壤全氮是影响脲酶的重要因素，土壤全氮及速效磷可直接影响酸性磷酸酶。

[1]陈隆升, 杨小胡, 李志刚, 等. 间种迷迭香对油茶幼林生长及病虫害的影响[J]. 中南林业科技大学学报, 2016,36(5):38-40.

[2]曾 静, 彭 秀, 李秀珍. 油茶林地间作研究[J]. 绿色科技,2013(2): 51-52.

[3]滕维超. 油茶—农作物间作系统生理生态及经济效益评价[D].南京:南京林业大学, 2013.

[4]郑 灏, 杨志坚,冯金玲,等.不同林下套种模式对油茶幼林根区土壤化学特性及微生物的影响[J].福建农林大学学报:自然版, 2015, 44(2): 147-153.

[5]王宪奎, 李建贵, 刘隋赟昊, 等. 不同施肥措施对灰枣园土壤速效养分含量的影响[J]. 经济林研究, 2016, 34(2):35-40.

[6]Abdul Kadir Salam, Arata Katayama, Makoto Kimura. Activities of some soil enzymes in different land use systems after deforestation in hilly areas of West Lampung, South Sumatra,Indonesia[J]. Soil Science and Plant Nutrition,1998,44(1):93-103.

[7]Badiane N N Y, Chotte J L, Pate E, et al. Use of soil enzyme activities to monitor soil quality in natural and improved fallows in semi-arid tropical regions[J]. Applied Soil Ecology, 2001,18(3):229-238.

[8]Acostamartinez V, Zobeck T M, Gill T E, et al. Enzyme activities and microbial community structure in semiarid agricultural soils[J]. Biology and Fertility of Soils, 2003, 38(4):216-227.

[9]关松荫. 土壤酶及其研究法[M]. 北京: 农业出版社, 1986.

[10]杨 瑞, 刘 帅, 王紫泉, 等. 秦岭山脉典型林分土壤酶活性与土壤养分关系的探讨[J]. 土壤学报, 2016,53(4):1037-1046.

[11]李 娟, 林位夫, 周立军. 成龄胶园间作不同姜科作物对土壤养分与土壤酶的影响[J]. 中国农学通报,2014,30(31):192-198.

[12]Wei L, Hao H L, Wu W, et al. Transgenic Bt rice does not affect enzyme activities and microbial composition in the rhizosphere during crop development[J]. Soil Biology & Biochemistry, 2008,40(2): 475-486.

[13]彭 秀, 曾 静, 李秀珍. 间作对油茶幼林生长及土壤养分的影响[J]. 四川林业科技, 2013, 34(4): 14-17.

[14]贾树海, 李 明, 邢兆凯, 等. 不同农林复合模式对土壤理化性质及酶活性的影响[J]. 土壤通报, 2014, (3): 648-653.

[15]田亚玲, 曹福亮. 银杏-茶间作模式对土壤养分和酶活性的影响[J]. 林业科技开发, 2012, 26(5): 41-45.

[16]张萌萌, 敖 红, 李 鑫, 等. 桑树/苜蓿间作对根际土壤酶活性和微生物群落多样性的影响[J]. 草地学报, 2015, 23(2):302-309.

[17]唐秀梅, 钟瑞春, 蒋 菁, 等. 木薯/花生间作对根际土壤微生态的影响[J]. 基因组学与应用生物学, 2015, 34(1): 117-124.

[18]杨敬天, 苏智先, 胡进耀, 等. 珙桐林土壤有机质与酶活性的通径分析[J]. 应用与环境生物学报, 2010, 16(2): 164-167.

[19]滕维超, 刘少轩, 曹福亮, 等. 油茶大豆间作对盆栽土壤化学和生物性质的影响[J]. 中南林业科技大学学报, 2013, 33(2):24-27.

[20]Soil Acid Phosphomonoesterase Activity to Simulated Nitrogen Deposition in Three Forests of Subtropical China[J]. Pedosphere,2012, 22(5): 698-706.

[21]罗汉东, 雷先高, 朱丛飞, 等. 不同施钾水平对油茶树体生长和林地养分含量的影响[J]. 经济林研究, 2016, 34(2): 1-6.

[22]吴 雪, 贡 璐, 冉启洋, 等. 阿拉尔垦区土壤理化因子与酶活性的通径分析[J]. 水土保持研究, 2013, 20(3): 48-54.

Soil enzyme activities and their relations with soil fertility in Camellia oleifera peanut intercropping

MO Jing1,2, YAN Wende1,2, LIU Shuguang1,2,3, WU xiaohong1,2
(1. Central South University of Forestry and Technology, Changsha 410004, Hunan, China; 2. National Engineering Lab for Applied Technology of Forestry and Ecology in Southern China, Changsha 410004, Hunan, China; 3. Earth Resource Observation and Science Center, United States Geological Survey, Sioux Falls SD57198, South Dakota)

In order to explore the effect of different cropping patterns of Camellia oleifera on enzyme activity, and to afford the reference of fertilization and improving of soil during C. oleifera production. This study compared the difference of fi ve soil enzyme activities(sucrose, protease, urease, acid phosphatase and catalase) under three cropping patterns: C. oleifera monoculture, peanut monoculture and C. oleifera peanut intercropping, and explored the relationships between soil enzyme activities and soil fertility using correlation analysis and path analysis. Compared with C. oleifera monoculture, C. oleifera peanut intercropping can significantly improve the activities of sucrose, protease, urease and catalase in 0-10cm soil layer and sucrose activity in 10-20 cm soil layer. Compared with peanut monoculture, C. oleifera peanut intercropping can signi fi cantly improve protease activity in 0-10 cm soil layer. Correlation analysis result indicated that the relationship between soil enzyme and soil fertility were signi fi cant, especially for urease and acid phosphatase.The result of path analysis showed that soil rapidly available phosphorus and rapid available potassium have a signi fi cant effect for sucrose, soil organic matter and soil rapidly available phosphorus were dominant factors in effecting protease, soil total nitrogen was important for urease and the soil total nitrogen and rapidly available phosphorus may directly affect the acid phosphatase. In conclusion, C.oleifera peanut can improve the quality of soil and enzyme activity comparing with C. oleifera monoculture, and show signi fi cant mean in promoting soil fertility.

Camellia oleifera; intercropping; soil enzyme activity; soil fertility

S714

A

1673-923X(2017)06-0089-07

10.14067/j.cnki.1673-923x.2017.06.015

2016-11-16

国家科技支撑计划（2015BAD07B00）

莫 晶，硕士研究生

闫文德，教授，博士生导师；E-mail：csfuyywd@hotmail.com

莫 晶，闫文德，刘曙光，等. 油茶-花生间作土壤酶活性与养分的关系[J].中南林业科技大学学报，2017, 37(6): 89-95.

[本文编校：吴 彬]