环境科学  2020, Vol. 41 Issue (5): 2425-2433   PDF    
施用生物炭对塿土土壤微生物代谢养分限制和碳利用效率的影响
王强1, 耿增超1,2, 许晨阳1,2, 郭靖宇1, 李倩倩1, 刘莉丽3, 赵汉红4, 杜旭光5     
1. 西北农林科技大学资源环境学院, 农业部西北植物营养与农业环境重点实验室, 杨凌 712100;
2. 农业部农业环境重点实验室, 北京 100081;
3. 西北农林科技大学林学院, 杨凌 712100;
4. 汉中市勉县农村能源工作站, 汉中 724200;
5. 汉中市勉县农业技术推广站, 汉中 724200
摘要: 土壤微生物代谢对土壤养分循环和生态系统的稳定至关重要.为明确施加生物炭对土壤微生物代谢养分限制和碳利用效率(carbon use efficiency,CUE)的长效影响机制,于2012年将果树树干、枝条生物炭(450℃、限氧条件下裂解)以不同用量(0、20、40、60和80 t·hm-2)施入塿土,与耕层土壤(0~20 cm)混匀,小麦玉米轮作7 a后,通过生态酶化学计量学对土壤微生物代谢养分限制特征进行了定量分析和比较.结果表明:①随生物炭施用量的增加,土壤含水量、有机碳、全氮、碳氮比、碳磷比和氮磷比显著提高,碳氮磷活性组分、微生物生物量碳氮磷和总磷未表现出明显的规律性,而5种胞外酶活性(β-1,4-葡萄糖苷酶、纤维素酶、亮氨酸氨基多肽酶、β-1,4-N-乙酰氨基葡萄糖苷酶和磷酸酶)显著降低.②所有处理土壤微生物均受磷限制;在施加生物炭各处理中,随施用量的增加微生物代谢碳和磷限制显著提高,微生物CUE显著降低;当生物炭施用量为20 t·hm-2时,碳限制(0.625±0.022)和磷限制(62.153°±0.892°)最低,微生物CUE(0.511±0.007)最高.③偏最小二乘路径建模分析表明,土壤碳氮磷及其元素化学计量比对磷限制产生了直接的极显著正效应(P < 0.01),碳限制与磷限制呈正相关关系(R2=0.242,P < 0.001),而碳磷限制又对CUE产生了极显著的负效应(P < 0.001).综上,过量施用生物炭使土壤元素化学计量失衡是导致土壤微生物代谢磷限制加剧的重要因素,继而诱导了微生物碳限制的增强和CUE的降低.当生物炭施加量为20 t·hm-2时,微生物代谢所受碳磷限制最低,且具有最高的微生物CUE,对于调节土壤微生物代谢、维持生态功能和减少微生物二氧化碳排放最优.
关键词: 生物炭(BC)      生态化学计量学      矢量分析      微生物养分限制      微生物碳利用效率     
Effects of Biochar Application on Soil Microbial Nutrient Limitations and Carbon Use Efficiency in Lou Soil
WANG Qiang1 , GENG Zeng-chao1,2 , XU Chen-yang1,2 , GUO Jing-yu1 , LI Qian-qian1 , LIU Li-li3 , ZHAO Han-hong4 , DU Xu-guang5     
1. Key Laboratory of Plant Nutrition and Agri-environment in Northwest China, Ministry of Agriculture, College of Natural Resources and Environment, Northwest A&F University, Yangling 712100, China;
2. Key Laboratory for Agricultural Environment, Ministry of Agriculture, Beijing 100081, China;
3. College of Forestry, Northwest A&F University Yangling 712100, China;
4. Rural Energy Workstation in Mianxian County, Hanzhong 724200, China;
5. Agrotechnical Promotion Station in Mianxian Country, Hanzhong 724200, China
Abstract: Soil microbial metabolism is vital for nutrient cycling and stability of an ecosystem. To elucidate the long-term effects of biochar application on nutrient limitations and carbon use efficiency (CUE) of soil microbial metabolisms, biochars pyrolyzed at 450℃ from trunks and branches of fruit trees under an oxygen-limited condition were mixed with the top Lou soils (0-20 cm) with application amounts of 0, 20, 40, 60, and 80 t·hm-2 in 2012. Corn-wheat rotation was carried out afterwards for seven years. The nutrient limitations of soil microbial metabolisms were analyzed quantitatively through ecoenzymatic stoichiometry in 2019. The results indicated that:① With an increase in the biochar application amount, soil moisture, organic carbon, total nitrogen, C:N, C:P, and N:P significantly increased, whereas there were no clear patterns for the active components of carbon, nitrogen, and phosphorus, microbial biomass carbon, nitrogen, phosphorus and total phosphorus. In contrast, the activities of five extracellular enzymes (β-1, 4-glucosidase, cellobiohydrolase, leucine aminopeptidase, β-1, 4-N-acetylglucosaminidase, and phosphatase) were significantly reduced. ② The soil microorganisms suffered from the phosphorus limitation under all treatments in this study. In the treatments of biochar application, the carbon and phosphorus limitations of microbial metabolisms increased significantly with increasing application amount, whereas the microbial CUE decreased significantly. When the application amount was 20 t·hm-2, the carbon limitation (0.625±0.022) and phosphorus limitation (62.153°±0.892°) were lowest, and the microorganism CUE (0.511±0.007) was highest. ③ Partial least-squares path modeling (PLS-PM) showed that soil carbon, nitrogen, phosphorus, and their stoichiometry had a very direct positive effect on phosphorus limitation (P < 0.01), and there was a positive correlation between carbon limitation and phosphorus limitation (R2=0.242, P < 0.001); in contrast, the carbon and phosphorus limitations had a very significant negative effect on CUE (P < 0.001). It was revealed that the excessive application of biochar had caused a soil element stoichiometry imbalance, which deteriorated the phosphorus limitation of the soil microbial metabolism and further led to carbon limitation and reduction of CUE. When the biochar application amount was 20 t·hm-2, C and P limitations were lowest, and microbial CUE was highest. Therefore, 20 t·hm-2 was optimal for regulating soil microbial metabolism, maintaining ecological functions, and reducing carbon dioxide emission produced by microbial metabolism.
Key words: biochar (BC)      ecological stoichiometry      vector analysis      microbial nutrient limitation      carbon use efficiency     

生态化学计量比作为生态系统的组成部分, 是限制生物地球化学循环的重要因素, 在评估生态系统稳定性、碳循环和微生物代谢限制上得到了广泛应用[1, 2].生态化学计量学认为土壤中元素化学计量比(土壤碳氮比、碳磷比和氮磷比)强烈影响生态系统中土壤微生物代谢碳(C)、氮(N)和磷(P)限制, 而且原因可能是由于养分的不平衡供应[3].Delgado-Baquerizo等[4]的研究表明, 细菌多样性及结构组成主要是由土壤元素化学计量的变化(C :N :P)所驱动.因此, 土壤元素化学计量学通过影响元素平衡和微生物群落结构, 强烈调节微生物群落的代谢, 可能是导致微生物代谢受限的重要因素.

生物炭(biochar, BC)是由生物质材料(如农业和林业废料)在限氧条件下热解而成[5].起初由于其在土壤中的长期持久性, 被用来减少土壤中温室气体的排放[5];近年来被广泛应用于改善土壤理化特性、吸附农田土壤有机无机污染物以提高土壤质量[6].有研究表明, 施用生物炭可以改善土壤孔隙、团聚体和养分含量等理化性质, 促进农田土壤增汇减排, 改善水热特性, 对维持农田土壤生态系统平衡有着不可忽视的作用[7, 8].陕西关中平原是我国主要的粮食产区, 土壤类型以塿土为主, 该土壤有机质含量偏低, 且由于化肥的长期使用导致农田土壤质量下降, 故农业上常使用生物炭来改善土壤质量[9], 但目前生物炭对塿土土壤微生物代谢活动的影响尚未可知.因此, 本研究基于陕西关中塿土的田间定位试验, 分析了果树树干、枝条制成的生物炭施加7a后对土壤微生物代谢养分限制和碳利用效率(carbon use efficiency, CUE)的影响, 从生态酶化学计量角度探究生物炭的施加对土壤微生物代谢的长效影响机制, 对塿土最适生物炭添加量进行评估, 并通过全面和系统的土壤理化指标来完善生物炭添加对塿土土壤质量影响的数据支撑和科学依据.

1 材料方法 1.1 试验区概况

试验地位于陕西省杨凌区西北农林科技大学试验田(34°16′ N, 108°04′ E), 海拔458.60 m, 属温带大陆性季风气候, 年均日照时数2 196 h, 年均气温11~13℃, 年均降雨量500~700 mm.土壤类型为褐土类, 塿土亚类, 红油土属, 黄土母质, 系统分类为土垫旱耕人为土, 属石灰性土壤.试验地耕层(0~20 cm)土壤基础理化特性(平均值)如下:pH为7.33, 容重为1.30 g ·cm-3, 有机质(SOM)为14.99 g ·kg-1, 全氮0.71 g ·kg-1, 总磷0.43 g ·kg-1, 全钾19.64 g ·kg-1, 硝态氮18.20 mg ·kg-1, 铵态氮15.90 mg ·kg-1, 速效磷12.38 mg ·kg-1, 速效钾193.00 mg ·kg-1.

1.2 供试材料及试验设计

本研究所用生物炭由废弃果树树干、枝条在裂解炉450℃、限氧条件下热裂解所得(陕西亿鑫生物能源科技开发有限公司), 磨细过1 mm筛, 备用.供试生物炭pH为10.43(土水比为1 :10), 比表面积为86.70 m2 ·g-1, 碳、氢、氧和氮的质量分数分别为72.38%、2.62%、23.81%、1.19%, 灰分的质量分数为13.98%, 硝态氮和铵态氮含量分别为0.52 mg ·kg-1和1.86 mg ·kg-1.

生物炭于2012年4月一次性均匀撒施(此后不再施入生物炭)并浅翻, 与耕层土壤均匀混合.根据用量(t ·hm-2)设5个水平:0(B0)、20(B20)、40(B40)、60(B60)和80(B80), 采用随机区组试验设计, 3次重复, 小区面积为1.35 m×3.50 m, 共计15个小区.各小区除生物炭用量不同外, 氮、磷和钾肥(分别为尿素、磷酸二铵和硫酸钾)均作基肥于每次种植前施入, 其用量折合成纯养分用量分别为225 kg ·hm-2(以N计)、180 kg ·hm-2(以P2O5计)和150 kg ·hm-2(以K2O计).

2012年4月至2019年6月在供试土壤上进行玉米-小麦轮作, 最后一季小麦收获后, 于2019年6月5日用土钻采集耕层(0~20 cm)土壤, 除去杂质后过2 mm筛, 一部分于4℃保存用于测定土壤胞外酶活性及微生物生物量碳、氮和磷, 另一部分自然风干用于测定土壤基本理化性质, 同时采集土壤容重样品, 每个小区重复3次.

1.3 测定方法

土壤含水量采用105℃烘干法测定;土壤容重采用环刀法测定;土壤pH采用pH计测定(土水比为1 :2.5);土壤有机碳(SOC)采用重铬酸钾氧化-外加热法测定;全氮(TN)采用半微量开氏法测定(K12型全自动凯氏定氮仪);土壤总磷(TP)采用HClO4+H2SO4消煮-钼锑抗比色法测定;可溶性有机碳(DOC)采用去离子水浸提-TOC-VCPH自动分析仪(岛津, 日本)测定;硝、铵态氮采用1 mol ·L-1 KCl浸提-连续流动分析仪法测定;土壤速效磷采用0.5 mol ·L-1 NaHCO3浸提-钼锑抗比色法测定.

采用氯仿熏蒸法[10]测定微生物生物量碳、氮和磷(MBC、MBN和MBP).简言之, 将新鲜土壤样品分成两个子样品, 一组直接用0.5 mol ·L-1 K2SO4提取碳和氮, 用0.5 mol ·L-1 NaHCO3提取磷, 另一组用氯仿熏蒸后提取碳、氮和磷.对于微生物生物量碳、氮和磷, 转换系数分别为0.45、0.54和0.40.

土壤胞外酶活性测定采用96微孔酶标板荧光分析法[11, 12], 所用底物如表 1所示, 采用多功能酶标仪(Scientific Fluoroskan Ascent FL, Thermo)在激发波长365 nm、发射波长450 nm条件下测定.

表 1 土壤胞外酶缩写及所用底物 Table 1 Abbreviations and substrates of soil enzymes assayed

1.4 胞外酶活性矢量分析和微生物碳利用效率计算方法

为量化土壤微生物代谢C、N和P的限制, 对生态酶活性相对比例进行矢量分析[3], 根据胞外酶活性比例[(BG+CBH)/(BG+CBH+NAG+LAP)和(BG+CBH)/(BG+CBH+AP)]计算了所有处理的向量长度(碳限制)和角度(小于45°为氮限制, 大于45°为磷限制)[13].

向量长度为xy平方值之和的平方根, 其中x代表C与P获取酶的相对活性, y代表C与N获取酶的相对活性.

向量角度为原点到点(x, y)的直线的反正切.

使用生物地球化学平衡模型计算微生物碳利用效率(CUE)[14].土壤胞外酶活性(EEA)通过转换以nmol ·(g ·h)-1为单位(以SOM计, 下同).EEAC :N计算为(BG+CBH)/(NAG+LAP), EEAC :P计算为(BG+CBH)/AP. SOC :TN和SOC :TP的摩尔比被用作LC :NLC :P的估计值.微生物生物量C :N和C :P(BC :NBC :P)也按摩尔比计算:

式中, KC :NKC :P是基于C、N和P可用性的CUE的半饱和常数.该模型假定当可吸收营养物供应的比例与微生物生物量化学计量相匹配且生长效率与N和P供给相对于C的几何平均值成正比时, 增长率最大.假设所有模型的KC :NKC :P为0.5, CUEmax为0.6[14].

1.5 数据处理

采用Microsoft Excel 2016进行数据处理;SPSS 19.0进行不同处理间差异显著性分析, 显著性用One-way ANOVA(单因素方差分析)进行检验(P<0.05), 采用Duncan多重比较分析组间差异, 土壤各指标间的相关性采用Pearson相关系数表示;在R 3.6.0中利用“Plspm”包构造偏最小二乘路径模型(partial least squares path modelling, PLS-PM)确定土壤特性控制土壤微生物养分限制的可能路径;文中箱线图和散点图采用Origin 2019b作图.图表中数据均为平均值±标准差.

2 结果与分析 2.1 土壤理化性质、胞外酶活性及其化学计量分析 2.1.1 各处理土壤理化特性和元素化学计量比

生物炭施加7 a后, 随生物炭施加量的增加土壤含水量随之增加, 容重减小.土壤DOC、硝、铵态氮、速效磷、MBC、MBN和MBP并没有表现出一定的规律性(表 2), DOC在B40处理时含量最高, 为(265.83±3.03)mg ·kg-1;硝态氮在B60处理时含量达最高, 为(11.04±1.15)mg ·kg-1;铵态氮在B80处理时含量最高, 为(2.13±0.16)mg ·kg-1;MBC、MBN在B0时含量最高, 分别为(177.92±6.67)mg ·kg-1和(7.37±0.21)mg ·kg-1;MBP在B60时含量最高, 为(7.23±0.63)mg ·kg-1.

表 2 各处理土壤理化性质1) Table 2 Soil physicochemical properties under different biochar application amounts

与空白相比(表 2), 施加生物炭各处理显著提高了土壤有机碳、全氮、C :N、C :P和N :P, 并且在B60或B80时土壤有机碳、全氮、C :N、C :P和N :P达到最大, 分别为(29.02±0.86)g ·kg-1、(1.32±0.02)g ·kg-1、22.34±0.64、34.01±1.44和1.52±0.04, 在B0时SOC、TN、C :N、C :P和N :P均最小;土壤TP随生物炭施加量的增加未表现出明显的规律性.

2.1.2 各处理土壤胞外酶活性及其化学计量比

施加生物炭各处理均显著降低了5种胞外酶活性(表 3), 且随生物炭施加量的增加酶活性逐渐降低, 在B80处理下最小, BG、CBH、NAG、LAP和AP活性分别为(0.99×103±38.48)、(0.18×103±10.05)、(0.11×103±6.59)、(0.47×103±24.65)和(2.58×103±116.25)nmol ·(g ·h)-1;(BG+CBH)/(NAG+LAP)在B80(2.00±0.08)中显著高于其他处理(P < 0.05), (BG+CBH)/AP在B0(0.47±0.02)和B80(0.45±0.02)处理中显著高于其他处理, 而(NAG+LAP)/AP在B20(0.33±0.02)最大, 显著高于其他处理.

表 3 各处理土壤胞外酶活性及酶活性比值1) Table 3 Soil extracellular enzyme activities and enzyme activity ratios under different biochar application amounts

2.2 生态酶化学计量矢量分析

不同生物炭施加量各处理间生态酶化学计量特征存在差异[图 1(a)], 这表明生物炭强烈影响微生物养分的供应.如图 1(a)所示, 所有数据点均高于虚线(1 :1线), 这表明土壤微生物群落受到磷的严重限制.本研究通过使用胞外酶活性相对比例计算矢量的长度和角度来量化微生物碳、氮和磷的限制[图 1(b)图 1(c)].向量长度(从0.587~0.746)和角度(从60.987°~65.739°)差异显著(P < 0.05).B20处理矢量长度和角度最低, 分别为(0.625±0.023)和(62.153°±0.892°), 表明B20时微生物受到的碳限制和磷限制最低.B80的矢量长度(0.736±0.007)最高, 表明相对碳限制最高;同样, B80的矢量角度(64.904°±0.690°)最高, 表明相对磷限制最强.而且矢量角度与矢量长度呈正相关关系[图 1(d)], 随着生物炭施加量的增加, 土壤微生物代谢所受碳限制和磷限制趋于严重.

图 1 碳氮获取酶与碳磷获取酶相对比例的生态酶化学计量、矢量长度和角度的变化及其关系 Fig. 1 Ecoenzymatic stoichiometry of relative proportions of C to N acquisition versus C to P acquisition, the variations in vector length, angle, and their relationships

2.3 各处理土壤微生物CUE

土壤微生物CUE在施加生物炭各处理间差异显著(图 2), 均值为0.469(F=173.429, P < 0.001), 在B20和B60时达到最高和最低值, 分别为0.511±0.007和0.412±0.011. B20和B40与空白相比, 微生物CUE差异不显著, 但随着生物炭施加量的增加, 微生物CUE显著降低.

图 2 不同处理下土壤微生物群落的碳利用效率(CUE) Fig. 2 Box plots showing CUE of microbial communities under different biochar application amounts

2.4 酶活性和生态酶化学计量特征与土壤特性的相关性分析

相关性分析表明(表 4), 5种土壤胞外酶活性与含水量、pH、SOC、TN、TP、C :N、C :P、N :P、MBC、MBN、硝态氮和速效磷表现出显著或极显著相关关系, DOC与BG、铵态氮与LAP呈显著相关关系, C :N-酶和N :P-酶与含水量、容重、SOC、TN、TP、C :N、C :P、N :P、铵态氮和速效磷表现出显著或极显著相关关系, C :N-酶和C :P-酶与MBN和DOC呈极显著相关关系.

表 4 土壤理化性质和微生物生物量与微生物代谢特征相关性分析1) Table 4 Correlations between soil physicochemical properties and microbial biomass/microbial metabolic characteristics

将矢量长度、矢量角度和CUE与土壤物理化学指标进行相关性分析表明(表 4), 矢量长度(碳限制)、矢量角度(磷限制)和CUE均与土壤含水量、容重、SOC、TP、C :N、C :P、N :P和速效磷呈极显著相关关系;且矢量角度和矢量长度均与CUE呈极显著负相关关系;DOC与矢量长度(碳限制)呈极显著负相关.

2.5 土壤特性对微生物代谢碳磷限制以及CUE的影响

PLS-PM分析确定了土壤物理化学特性和微生物生物量对微生物代谢碳和磷限制以及CUE的直接和间接影响.

施加生物炭改变了土壤养分供应及元素计量比等土壤理化特性, PLS-PM分析表明, 相对于土壤碳氮磷含量及其计量比(-0.045), 活性组分含量(-0.191)对于矢量长度(碳限制)的影响较大[图 3(a)];而对于矢量角度(磷限制)而言, 土壤碳氮磷含量及其计量比(0.763)对矢量角度产生了直接的极显著正效应[图 3(b)].

偏最小二乘路径建模(PLS-PM)解释了土壤物理化学特性和微生物生物量影响土壤微生物代谢养分限制和CUE的主要途径;实线和虚线箭头分别表示因果关系的正负(P < 0.05);箭头上的数字表示标准化路径系数, 星号表示显著性;R2表示由模型解释的因变量的方差; 基础特性包括含水量、容重和pH;碳氮磷及其计量比包括SOC、TN、TP、C :N、C :P和N :P;活性组分包括DOC、硝、铵态氮和速效磷;微生物生物量包括MBC、MBN和MBP;碳磷限制包括矢量角度(磷限制)和矢量长度(碳限制);*表示显著性水平0.05;**表示显著性水平0.01;***表示显著性水平0.001 图 3 偏最小二乘路径建模分析 Fig. 3 Partial least-squares path modelling (PLS-PM)

此外, PLS-PM分析还表明[图 3(c)], 微生物代谢碳磷限制(-0.326)和土壤碳氮磷含量及其计量比(-0.410)对微生物CUE产生了极显著的负效应, 而且活性组分含量(0.630)对微生物代谢养分限制产生了极显著正效应.

3 讨论 3.1 生物炭对土壤微生物代谢活动的影响

生物炭的施加均显著降低了5种胞外酶的活性, 不施加生物炭时胞外酶活性最高, 表明其微生物C、N和P的代谢水平最高.生物炭对微生物代谢活性的影响可能是由于生物炭改变了土壤理化特性和养分状况[7, 8, 9], 进而影响微生物对养分的利用, 造成土壤微生物的生物量和群落结构发生变化[15], 继而对胞外酶活性产生影响.

生物炭对土壤碳氮磷含量及其元素计量比产生了不同程度的影响(表 2), 而且相关性分析显示SOC、TN和TP含量及其化学计量比与胞外酶活性和生态酶化学计量存在极显著相关关系(表 4), 表明土壤碳氮磷与胞外酶活性和生态酶化学计量的变化存在紧密联系.Sinsabaugh等[16, 17]的研究表明, 土壤胞外酶活性和生态酶化学计量通常受土壤养分状况的影响, 当养分利用率低时, 微生物可以分泌更多的相关酶来满足养分需求[16], 并且土壤养分含量会通过改变可利用的底物和C、N、P的化学计量而对胞外酶活性和生态酶化学计量产生重大影响[18].因此, 笔者认为生物炭的添加通过改变土壤营养供应状况极大地影响了土壤微生物的代谢活动和养分需求, 从而对土壤胞外酶活性和生态酶化学计量比产生影响, 甚至可能会影响农田生态系统的稳定性.

3.2 不同生物炭施用量条件下微生物的碳磷限制

长期的生态观测和试验表明, 磷是许多生态系统的主要限制元素[19].本研究中微生物群落同样受到磷的限制.在5个生物炭处理(所有矢量角度>45°)下, 微生物群落均受到磷的限制[图 1(a)].磷限制分别在B20和B80时达到最低和最高值[图 1(c)].

生物炭施加7 a后仍显著提高了土壤碳氮含量, 而土壤总磷甚至有所降低(表 2), 这就会对土壤养分计量比(C :N、C :P和N :P)产生很大影响, PLS-PM分析表明, 土壤碳氮磷及其计量比对矢量角度(磷限制)产生了极显著的正效应, 活性碳氮磷组分对矢量角度产生了显著负效应[图 3(b)], 说明生物炭的施加一方面通过提高土壤碳氮含量及元素计量比加剧了土壤微生物代谢的磷限制, 另一方面通过调节速效养分含量缓解了磷的限制.元素化学计量学可以强烈影响微生物群落的结构和活性[20, 21].Delgado-Baquerizo等[4]的研究表明, 土壤元素化学计量的变化(C :N :P)驱动了细菌多样性和组成, 而且土壤元素化学计量也可以控制丛枝菌根与根源非菌根真菌的竞争能力[22].因此, 笔者认为过量生物炭的施加在一定程度上加剧了磷的限制, 主要是由于生物炭改变了土壤碳氮磷含量及其计量比所造成的.

本研究中微生物群落同样受到碳的限制[图 1(a)], 分别在B20和B80时所受限制最低和最高[图 1(b)].不同生物炭施加量对土壤胞外酶活性和微生物生物量有强烈影响(表 2表 3), 这表明生物炭通过影响微生物群落结构和土壤养分供应而间接影响微生物代谢养分限制.施加生物炭各处理SOC含量的变化会对微生物碳源的供应产生影响[15].

在本研究中, B80的SOC含量显著高于其他处理(表 2), 理论上讲, 应该具有更多的碳源来供应微生物代谢, 微生物群落所受碳限制也会随之减小.但本研究的结果表明随着生物炭施加量的持续增加土壤微生物代谢受到的碳限制愈加严重, 这一方面可能是由于生物炭具有碳含量高、芳香化结构复杂和固有的化学惰性等特征[23], 属惰性碳库, 在土壤环境中具有较高的化学和微生物稳定性[24], 此外, 生物炭的施用可能改变了土壤原有有机质的组成, 形成比较稳定的土壤有机质[25], 而且本研究得出DOC的含量并没有随生物炭施加量的增加而持续增加(表 2), 对于微生物代谢碳的限制, DOC的影响最大, 因为它的含量决定了土壤碳的有效性.进一步分析表明, DOC与矢量长度(碳限制)呈极显著负相关(表 4), 随着DOC含量的增加微生物代谢所受碳限制减弱, 而且铵态氮和速效磷与矢量长度呈极显著正相关关系, 也就是说随着铵态氮和速效磷含量的增加, 根据养分平衡原理土壤微生物代谢所受碳限制增强.PLS-PM分析也表明土壤碳氮磷及其计量比和活性组分均对微生物代谢碳限制产生了负效应[图 3(a)], 且速效养分对碳限制的影响较大.

另一方面, 最可能的解释就是微生物磷限制加剧了微生物代谢的碳限制, 线性回归分析也显示微生物代谢碳限制和磷限制呈正相关[图 1(d)].前人的研究表明, 如果磷成为相对受限的元素, 则微生物群落可能会消耗更多的碳和氮来生产和排泄与磷代谢相关的酶[26].土壤酶活性与土壤理化性质之间的强相关性(表 4)还表明, 微生物可以通过在酶促合成过程中优化碳氮磷的分配来获得相对稀缺的资源(磷)[20].另外, 微生物所需要的大部分磷来自土壤有机质的分解[27], 因此, 微生物对磷的需求可以诱导碳相关酶的分泌和微生物代谢的碳限制, 因为微生物分泌更多的碳摄取酶来水解有机物, 并减轻土壤磷的限制.因此, 土壤磷的限制在一定程度上控制着微生物的碳代谢.

Liang等[28]的研究表明生物炭可增强土壤生物活性和团聚性, 改善土壤结构, 提高土壤持水特性, 而土壤含水量和容重也是影响生物炭分解和微生物活动的重要因素, 从而影响微生物的代谢和养分吸收.

因此, 笔者认为适量的生物炭添加有助于缓解塿土土壤微生物代谢所受的养分限制, 但过量施加生物炭会影响微生物代谢养分供应水平继而加剧这种限制.

3.3 微生物代谢碳和磷限制对微生物CUE的影响

土壤微生物代谢所受碳和磷限制会对生态平衡产生一定影响, 这种限制可以加速土壤有机质的分解, 提供有效的碳和磷, 刺激土壤碳的释放[20].地球系统模型还表明, 氮和磷限制在一定程度上减少了土壤碳汇, 使得一些区域性碳汇变成了净碳源[29].许多生态系统都已证明了微生物碳代谢所受的磷限制, 该限制增加了溶解性有机物的微生物矿化作用, 并增加了土壤CO2的外流[30].因此, 微生物的代谢限制极大地影响了土壤碳的固存.

为了定量评估微生物代谢限制对土壤碳循环的生态影响, 本研究引入了关键性指标土壤微生物CUE, 作为将土壤微生物代谢与有机物分解联系起来的定量指标.

土壤微生物CUE在施加不同量的生物炭过程中变化显著, 且随施加量的增加在B40和B60之间发生剧烈降低.土壤碳氮磷及其化学计量比可以通过改变微生物生物量和微生物代谢限制来显著影响微生物的CUE.而且由于分解者细胞需要保持碳和营养成分的平衡以及微生物生物量的稳态, 元素化学计量(尤其是碳)的波动将会对土壤微生物的生长和呼吸产生影响[31, 32], 这会进一步改变CUE.

本研究将微生物代谢所受碳磷限制与微生物CUE联系起来发现, 微生物CUE与微生物碳磷限制呈极显著负相关(表 4), 这表明它们具有高度依赖性.通过PLS-PM进一步分析还表明[图 3(c)], 微生物碳磷限制对微生物CUE具有强烈的负效应.随着微生物代谢碳限制的增加, CUE的下降可能归因于微生物从生长到维持呼吸的转变, 以及随着有机质分解的加快, 增加了碳对微生物产生酶的投入[31].微生物代谢限制是代谢的重要特征, 总体而言, 它通过以下方式影响微生物CUE, 可直接利用的土壤养分直接参与微生物的生长和代谢, 当可利用的营养不足时, 例如本研究区域的碳和磷限制, 微生物必须首先分解复杂化合物[33], 这就需要大量酶促反应降解化合物, 从而降低了转化为新生物质的效率.另外, 各种底物需要不同的代谢途径才能被微生物完全吸收[34], 这可能导致单位碳的转化率更低.因此, 该过程将大大抵消固定在土壤中的一些碳, 使其通过呼吸作用以二氧化碳形式从土壤中大量释放.

4 结论

(1) 本研究通过生态酶化学计量学分析表明, 研究区各处理土壤微生物均受磷的限制;

(2) 过量施用生物炭使土壤元素化学计量失衡是导致土壤微生物代谢磷限制加剧的重要因素, 继而诱导了微生物碳限制的增强和CUE的降低;

(3) 当生物炭施加量为20 t ·hm-2时, 土壤微生物代谢所受碳和磷限制相对最低, 且CUE最高, 对于调节微生物代谢、维持生态功能和减少土壤微生物二氧化碳排放最优.

参考文献
[1] Cui Y X, Zhang Y L, Duan C J, et al. Ecoenzymatic stoichiometry reveals microbial phosphorus limitation decreases the nitrogen cycling potential of soils in semi-arid agricultural ecosystems[J]. Soil and Tillage Research, 2020, 197: 104463.
[2] Xu Z W, Yu G R, Zhang X Y, et al. Soil enzyme activity and stoichiometry in forest ecosystems along the north-south transect in eastern China (NSTEC)[J]. Soil Biology and Biochemistry, 2017, 104: 152-163.
[3] Cui Y X, Fang L C, Guo X B, et al. Natural grassland as the optimal pattern of vegetation restoration in arid and semi-arid regions:evidence from nutrient limitation of soil microbes[J]. Science of the Total Environment, 2019, 648: 388-397.
[4] Delgado-Baquerizo M, Reich P B, Khachane A N, et al. It is elemental:soil nutrient stoichiometry drives bacterial diversity[J]. Environment Microbiology, 2017, 19(3): 1176-1188.
[5] 祁乐, 高明, 郭晓敏, 等. 生物炭施用量对紫色水稻土温室气体排放的影响[J]. 环境科学, 2018, 39(5): 2351-2359.
Qi L, Gao M, Guo X M, et al. Effects of biochar application rates on greenhouse gas emissions in the purple paddy soil[J]. Environment Science, 2018, 39(5): 2351-2359.
[6] Xie T, Reddy K R, Wang C W, et al. Characteristics and applications of biochar for environmental remediation:a review[J]. Critical Reviews in Environmental Science and Technology, 2015, 45(9): 939-969.
[7] 王月玲, 耿增超, 王强, 等. 生物炭对塿土土壤温室气体及土壤理化性质的影响[J]. 环境科学, 2016, 37(9): 3634-3641.
Wang Y L, Geng Z C, Wang Q, et al. Influence of biochar on greenhouse gases emissions and physico-chemical properties of Loess soil[J]. Environment Science, 2016, 37(9): 3634-3641.
[8] 尚杰, 耿增超, 赵军, 等. 生物炭对塿土水热特性及团聚体稳定性的影响[J]. 应用生态学报, 2015, 26(7): 1969-1976.
Shang J, Geng Z C, Zhao J, et al. Effects of biochar on water thermal properties and aggregate stability of Lou soil[J]. Chinese Journal of Applied Ecology, 2015, 26(7): 1969-1976.
[9] 李倩倩, 许晨阳, 耿增超, 等. 生物炭对塿土土壤容重和团聚体的影响[J]. 环境科学, 2019, 40(7): 3388-3396.
Li Q Q, Xu C Y, Geng Z C, et al. Impact of biochar on soil bulk density and aggregates of Lou soil[J]. Environment Science, 2019, 40(7): 3388-3396.
[10] Tapia-Torres Y, Elser J J, Souza V, et al. Ecoenzymatic stoichiometry at the extremes:how microbes cope in an ultra-oligotrophic desert soil[J]. Soil Biology and Biochemistry, 2015, 87: 34-42.
[11] 高雨秋, 戴晓琴, 王建雷, 等. 亚热带人工林下植被根际土壤酶化学计量特征[J]. 植物生态学报, 2019, 43(3): 258-272.
Gao Y Q, Dai X Q, Wang J L, et al. Characteristics of soil enzymes stoichiometry in rhizosphere of understory vegetation in subtropical forest plantations[J]. Chinese Journal of Plant Ecology, 2019, 43(3): 258-272.
[12] Duan C J, Razavi B S, Shen G T, et al. Deciphering the rhizobium inoculation effect on spatial distribution of phosphatase activity in the rhizosphere of alfalfa under copper stress[J]. Soil Biology and Biochemistry, 2019, 137: 107574.
[13] Deng L, Peng C H, Huang C B, et al. Drivers of soil microbial metabolic limitation changes along a vegetation restoration gradient on the Loess Plateau, China[J]. Geoderma, 2019, 353: 188-200.
[14] Sinsabaugh R L, Manzoni S, Moorhead D L, et al. Carbon use efficiency of microbial communities:stoichiometry, methodology and modelling[J]. Ecology Letters, 2013, 16(7): 930-939.
[15] 周凤, 耿增超, 许晨阳, 等. 生物炭用量对塿土土微生物量及碳源代谢活性的影响[J]. 植物营养与肥料学报, 2019, 25(8): 1277-1289.
Zhou F, Geng Z C, Xu C Y, et al. Effect of biochar addition on soil microbial biomass and metabolic activities of carbon sources in Lou soil[J]. Journal of Plant Nutrition and Fertilizers, 2019, 25(8): 1277-1289.
[16] Sinsabaugh R L, Lauber C L, Weintraub M N, et al. Stoichiometry of soil enzyme activity at global scale[J]. Ecology Letters, 2008, 11(11): 1252-1264.
[17] Sinsabaugh R L, Hill B H, Shah J J F. Ecoenzymatic stoichiometry of microbial organic nutrient acquisition in soil and sediment[J]. Nature, 2009, 462(7274): 795-798.
[18] Kivlin S N, Treseder K K. Soil extracellular enzyme activities correspond with abiotic factors more than fungal community composition[J]. Biogeochemistry, 2014, 117(1): 23-37.
[19] 陈美领, 陈浩, 毛庆功, 等. 氮沉降对森林土壤磷循环的影响[J]. 生态学报, 2016, 36(16): 4965-4976.
Chen M L, Chen H, Mao Q G, et al. Effect of nitrogen deposition on the soil phosphorus cycle in forest ecosystems:a review[J]. Acta Ecologica Sinica, 2016, 36(16): 4965-4976.
[20] 周正虎, 王传宽. 微生物对分解底物碳氮磷化学计量的响应和调节机制[J]. 植物生态学报, 2016, 40(6): 620-630.
Zhou Z H, Wang C K. Responses and regulation mechanisms of microbial decomposers to substrate carbon, nitrogen, and phosphorus stoichiometry[J]. Chinese Journal of Plant Ecology, 2016, 40(6): 620-630.
[21] Cui Y X, Fang L C, Guo X B, et al. Ecoenzymatic stoichiometry and microbial nutrient limitation in rhizosphere soil in the arid area of the northern Loess Plateau, China[J]. Soil Biology and Biochemistry, 2018, 116: 11-21.
[22] Chagnon P L, Bradley R L. Evidence that soil nutrient stoichiometry controls the competitive abilities of arbuscular mycorrhizal vs. root-borne non-mycorrhizal fungi[J]. Fungal Ecology, 2013, 6(6): 557-560.
[23] 叶协锋, 周涵君, 于晓娜, 等. 热解温度对玉米秸秆炭产率及理化特性的影响[J]. 植物营养与肥料学报, 2017, 23(5): 1268-1275.
Ye X F, Zhou H J, Yu X N, et al. Physiochemical properties and yields of corn-stalk-biochar under different pyrolyzed temperatures[J]. Journal of Plant Nutrition and Fertilizer, 2017, 23(5): 1268-1275.
[24] 高诚祥, 刘玉学, 汪玉瑛, 等. 生物炭的稳定性及其对矿物改性的响应机制研究进展[J]. 应用生态学报, 2019, 30(9): 3245-3251.
Gao C X, Liu Y X, Wang Y Y, et al. Stability of biochar and the mechanisms underlying its response to mineral modification:a review[J]. Chinese Journal of Applied Ecology, 2019, 30(9): 3245-3251.
[25] 尚杰, 耿增超, 陈心想, 等. 施用生物炭对旱作农田土壤有机碳氮及其组分的影响[J]. 农业环境科学学报, 2015, 34(3): 509-517.
Shang J, Geng Z C, Chen X X, et al. Effects of biochar on soil organic carbon and nitrogen and their fractions in a rainfed farmland[J]. Journal of Agro-Environment Science, 2015, 34(3): 509-517.
[26] Marklein A R, Houlton B Z. Nitrogen inputs accelerate phosphorus cycling rates across a wide variety of terrestrial ecosystems[J]. New Phytologist, 2012, 193(3): 696-704.
[27] Ru J Y, Zhou Y Q, Hui D F, et al. Shifts of growing-season precipitation peaks decrease soil respiration in a semiarid grassland[J]. Global Change Biology, 2018, 24(3): 1001-1011.
[28] Liang B, Lehmann J, Sohi S P, et al. Black carbon affects the cycling of non-black carbon in soil[J]. Organic Geochemistry, 2010, 41(2): 206-213.
[29] Exbrayat J F, Pitman A J, Zhang Q, et al. Examining soil carbon uncertainty in a global model:response of microbial decomposition to temperature, moisture and nutrient limitation[J]. Biogeosciences, 2013, 10(11): 7095-7108.
[30] Cui Y X, Fang L C, Deng L, et al. Patterns of soil microbial nutrient limitations and their roles in the variation of soil organic carbon across a precipitation gradient in an arid and semi-arid region[J]. Science of the Total Environment, 2019, 658: 1440-1451.
[31] Manzoni S, Taylor P, Richter A, et al. Environmental and stoichiometric controls on microbial carbon-use efficiency in soils[J]. New Phytologist, 2012, 196(1): 79-91.
[32] Zhu Z K, Ge T D, Luo Y, et al. Microbial stoichiometric flexibility regulates rice straw mineralization and its priming effect in paddy soil[J]. Soil Biology and Biochemistry, 2018, 121: 67-76.
[33] Cui Y X, Fang L C, Gao X B, et al. Responses of soil microbial communities to nutrient limitation in the desert-grassland ecological transition zone[J]. Science of the Total Environment, 2018, 642: 45-55.
[34] Van Hees P A W, Jones D L, Finlay R, et al. The carbon we do not see-the impact of low molecular weight compounds on carbon dynamics and respiration in forest soils:a review[J]. Soil Biology and Biochemistry, 2005, 37(1): 1-13.