土壤呼吸是土壤碳库向大气输出碳的主要途径，也是大气CO₂重要的源，在陆地生态系统碳循环过程中具有重要地位^[1]。土壤呼吸主要包括微生物的异养呼吸和植物根系的自养呼吸，两者约占土壤呼吸的95%^[2]。因此，明确土壤呼吸组分对土壤呼吸的贡献率，对准确评估土壤碳过程具有重要意义^[3]。近年来，国内外学者对森林和草地生态系统土壤呼吸组分进行了大量研究^[4-5]，而针对旱作农田风沙土土壤呼吸组分的研究相对较少。风沙土主要分布在中国干旱、半干旱地区，是新疆分布面积最大的土类，其有机碳含量相对较低，是陆地生态系统潜在的碳汇^[6-8]。提高农田风沙土有机碳含量和碳汇功能，对农业可持续发展和缓解气候变化具有重要意义。

生物质炭是指生物质在缺氧或厌氧条件下经高温裂解后生成的一种高度芳香化难熔性固态富碳物质^[9]，其高度生物化学稳定性在土壤固碳减排方面发挥重要作用，对全球碳循环具有重要的影响^[10]。生物质炭在土壤改良中得到了较多的应用，施用生物质炭显著提升了表层土壤有机碳、无机碳含量及碳储量^[11]，降低了土壤有机碳矿化速率和累积矿化量^[12]，表现出明显的固碳作用。生物质炭进入土壤后随着时间的推移发生氧化和表面结构破裂，官能团种类、数量、芳香性和极性等都会发生改变^[13]，从而影响生物质炭的固碳减排能力^[14]。作为外源有机物，生物质炭施入土壤，在短期内会对土壤原有有机质产生激发效应^[15]。即使生物质炭施入土壤的时间相同，其对土壤呼吸的影响仍存在促进^[16-17]和无影响^[18]的差异。此外，生物质炭对土壤呼吸组分的影响也存在较大的不确定性。生物质炭添加可促进土壤自养呼吸，降低土壤异养呼吸^[19-20]；但也有研究发现生物质炭会显著降低土壤自养呼吸^[21]，这可能与土壤类型、生物质炭性质以及水热因子等有关^[22-24]。目前，生物质炭在农田生态系统中得到较为广泛的应用，而有关生物质炭对风沙土土壤呼吸及组分影响的研究较少。

本试验利用前期已建立的定位试验平台，拟探明单次施用不同量生物质炭5年后土壤呼吸及组分日变化特征、玉米不同生育时期土壤呼吸以及影响土壤呼吸的水热因子和土壤呼吸温度敏感系数变化特征，以期为干旱区绿洲农田风沙土培肥提供基础数据。

1 材料与方法 1.1 研究区概况

试验于新疆石河子垦区兵团农八师121团炮台土壤改良试验站进行，试验站位于新疆准噶尔盆地绿洲边缘（84°58′~86°24′E，43°26′~45°20′N）。研究区属内陆干旱半荒漠气候，冬季长而严寒，夏季短而炎热，年平均气温7.5 ℃，日照2 525 h，无霜期169 d，年降雨量225 mm，年蒸发量1 250 mm。供试土壤类型为风沙土，土壤砂粒、粉粒和黏粒含量分别为53.2%、27.2%和19.6%。试验前土壤0~20 cm土层基本理化性质：有机质1.38 g·kg^-1，全氮0.76 g·kg^-1，碱解氮7.4 mg·kg^-1，有效磷4.6 mg·kg^-1，速效钾97 mg·kg^-1，pH 8.51。试验所用生物质炭是由河南三利新能源有限公司提供的小麦秸秆炭，在450 ℃限氧条件下炭化4~8 h，粉碎后过2 mm筛制成，基本理化性质：有机碳670 g·kg^-1，有效磷82.2 mg·kg^-1，速效钾1 590 mg·kg^-1，pH 9.90。

1.2 试验设计

试验始于2014年，采用随机区组设计，根据秸秆炭化量（2.625 t·hm^-2·a^-1）共设5个处理：单施化肥（T0），分别施加15.75（T1）、31.50（T2）、63.00 t·hm^-2（T3）和126.00 t·hm^-2（T4）生物质炭。每个处理设3次重复，共计18个小区，小区面积为4.6 m×7.0 m，地膜宽2.1 m，每个小区2行地膜，每行地膜3行滴灌带，6行玉米，小区间隔1 m，四周均有保护行。供试生物质炭于2014年播种前一次性撒施，并与表层0~20 cm土壤均匀混合。试验所用化肥施用量一致：尿素660 kg·hm^-2，过磷酸钙225 kg·hm^-2，磷酸二氢钾165 kg·hm^-2，硫酸钾15 kg·hm^-2，其中过磷酸钙和硫酸钾作为基肥一次性施入，其余氮肥和磷酸二氢钾随滴灌分次施入。在此期间，各处理农艺措施均保持一致。2019年5月16日播种，供试玉米品种为新玉53号，灌溉方式为滴灌，覆膜处理滴灌管置于膜下，种植体系为一年一茬，田间管理措施与当地保持一致。

1.3 土壤呼吸、土壤温度和含水率以及土壤容重的测定

土壤呼吸速率采用LI-8100土壤碳通量测量系统（LI-COR，Lincoln，美国）进行测定。每个小区内沿对角线设2个内径为20 cm、高10 cm的PVC管底座，将底座嵌入土壤，顶端距离地面3 cm，其中一个设置在玉米株间，株距为25 cm；另一个设置在垄间裸地，垄间宽度为50 cm。底座具体埋设位置如图 1所示。采用根系排除法测定土壤呼吸组分^[25]，玉米株间测定的土壤呼吸为土壤总呼吸（R_t），垄间裸地测定的土壤呼吸为土壤异养呼吸，用土壤总呼吸减去相应的土壤异养呼吸，即为玉米根系的自养呼吸。试验期间及时清理裸地及底座内所有地上植物，为减少嵌入底座短期内对土壤呼吸的扰动，在初次测量土壤呼吸速率前24 h固定好PVC管底座，固定好的底座留在试验小区内直至玉米生育期结束。于2019年6月22—24日（拔节期）、7月28—30日（抽穗期）、8月20—22日（灌浆期）和10月5—7日（成熟期）对各小区土壤呼吸速率进行监测。测定时段为07：00—22：00，每3 h监测一次，每个PVC管测量时间为120 s。

土壤温度采用地温计测定，监测前24 h在底座附近设置地温计，分别对5、10、15 cm处土壤地温进行监测。以底座为圆心，在半径50 cm的范围内，采集0~20 cm土壤用于含水率的测定（烘干法）。采用环刀法对土壤0~20 cm容重进行测定。

1.4 数据处理与统计分析

将土壤呼吸速率数值经对数转换后，采用线性方程模拟土壤呼吸与土壤温度之间的关系，见公式（1）:

(1)

式中：R_t为土壤总呼吸速率，μmol·m^-2·s^-1；T为土壤温 度，℃；a、b为拟合参数。

采用Jian等^[26]方法计算不同土层温度变化敏感性对土壤呼吸速率的响应，Q₁₀表示土壤温度升高10 ℃时，土壤呼吸速率变化的倍数，见公式（2）：

(2)

式中：R₁和R₂分别为对应温度T₁和T₂（℃）时的土壤呼吸速率，μmol·m^-2·s^-1。

采用Microsoft Excel 2016、SPSS 24.0和Origin2018分别进行数据处理、统计分析和绘图。采用LSD最小显著差数法检验不同处理间差异显著性（α= 0.05）；用Pearson相关系数表示土壤呼吸与水热因子之间的相关性；采用主成分分析法对土壤呼吸速率影响因子进行提取，并计算各因子的影响权重。

2 结果与分析 2.1 土壤呼吸日动态变化

玉米拔节期、抽穗期、灌浆期和成熟期，不同生物质炭处理的土壤呼吸速率日动态变化趋势基本一致，均呈先增后降的单峰曲线（图 2a）。T4、T3、T0、T2、T1的曲线峰值分别为3.24、2.28、2.22、2.12、2.06 μmol· m^-2·s^-1，且均出现在13：00左右；各处理土壤总呼吸速率平均值分别为2.74、2.09、2.04、1.89、1.88 μmol·m^-2· s^-1，T4处理显著高于其他处理（P＜0.05）。T1处理的土壤呼吸速率在19：00时降至1.69 μmol·m^-2·s^-1，随后又增加；其余处理的土壤呼吸速率最低值均出现在22：00。不同生物质炭处理土壤异养呼吸总体呈先增后降的趋势（图 2b），在13：00达到峰值。各处理土壤自养呼吸总体变化趋势不明显（图 2c），除T4外，其他各处理波动幅度较小。

2.2 不同生育时期土壤总呼吸变化

施用生物质炭的处理，土壤呼吸速率随着施炭量的增加而增加（图 3）。在玉米拔节期、抽穗期、灌浆期，T4处理土壤呼吸速率较T0显著增加了24.13%、26.22%和34.69%（P＜0.05），较T1处理显著增加了50.32%、32.95%和74.57%。拔节期、抽穗期、灌浆期和成熟期T1处理土壤呼吸速率均低于T0处理，说明单次施用生物质炭处理5年后可降低风沙土土壤呼吸速率。

玉米生育期内，不同处理间土壤呼吸速率表现为T4>T3>T0>T2>T1（图 4）。T1处理土壤总呼吸速率较T0处理降低了7.10%（P＜0.05），T4处理的土壤总呼吸速率较T1处理则显著增加了27.04%（P＜0.05）。不同处理间的土壤异养呼吸具有显著差异（P＜0.05），表现为T4>T0>T3>T2>T1，T0处理的异养呼吸占土壤总呼吸的比例为76.78%，显著高于其他处理（P＜0.05）；T1处理异养呼吸占土壤总呼吸的比例为59.89%，显著低于其他处理（P＜0.05）。各处理土壤异养呼吸对土壤总呼吸的贡献率为59.89%~ 76.78%，土壤自养呼吸对土壤总呼吸的贡献率为23.22%~40.11%。

2.3 土壤呼吸与水热因子的相关性及温度敏感系数

不同处理土壤呼吸速率与不同土层土壤温度显著相关（图 5）。除T4处理外，土壤呼吸速率与15 cm土层温度的相关性高于5 cm土层和10 cm土层。温度敏感系数Q₁₀随生物质炭施用量的增加呈增加趋势，T3、T4处理的Q10值显著高于T0、T1处理（图 6）。相关性分析表明，不同处理土壤呼吸与土壤含水率无显著相关性（P>0.05）。

2.4 土壤因子对土壤呼吸速率的贡献率分析

采用主成分分析法对5、10、15 cm土层温度及含水率和容重进行主成分提取，明确其对土壤呼吸速率影响的贡献率。由表 1可知，所选取的水热因子的贡献率达到了79.98%，主成分1和主成分2分别为55.83%和24.15%（表 1）。通过计算可知5、10、15 cm土层温度对土壤呼吸速率影响的权重最大。

3 讨论 3.1 施用生物质炭5年后风沙土土壤呼吸特征

施用生物质炭5年后风沙土土壤呼吸速率存在明显的日变化规律，呈先增后降的单峰曲线，由于新疆夏季昼夜温差较大，在7：00、19：00和22：00时根系呼吸较弱导致土壤呼吸速率较低。在13：00时土壤呼吸速率达到峰值，此时土壤酶活性较高，促进了植物根系的代谢速率^[27-28]。由此可知，生物质炭在施入土壤后，不同时刻作用效果存在显著差异，峰值过后土壤呼吸随着时间的推移缓慢降低，这与Smith等^[29]的研究结果相似。在玉米生育期测得不同处理土壤呼吸速率为T4>T3>T2>T1，随着生育期的推进则呈逐渐降低趋势，生物质炭施用量在一定阈值内可抑制土壤呼吸，达到减排的效果；但超过一定阈值，随着施用量的增加反而会促进土壤呼吸。成熟期的土壤呼吸显著低于其他3个时期，可能是由于玉米生长后期死根数量增加以及受氮肥供应水平的限制，减缓了作物的生长，从而降低土壤呼吸^[30]。

土壤呼吸分为自养呼吸和异养呼吸，区分两者在土壤呼吸中所占比例对明确土壤呼吸机理具有重要意义^[31]。对于不同生态系统，根系呼吸对土壤总呼吸的贡献率在10%~90%之间波动^[32]，根系呼吸变异较大。本研究结果表明，不同处理的土壤自养呼吸差异不显著，对土壤总呼吸的贡献率为23.22%~40.11%。李亚森等^[33]的研究发现，根系的自养呼吸对总呼吸的贡献率为34.55%~41.56%，这一变化略高于本研究结果，主要原因可能是本研究中生物质炭在土壤中发生了老化现象^[34]，其内部的结构和性质发生了改变，影响了土壤养分转化，从而影响玉米根系呼吸。同时本研究发现土壤异养呼吸随着生物质炭施用量的增加而增加，这一结果与Li等^[20]得出的结论相反。本研究中生物质炭在土壤中风化时间长达5年，一方面新鲜生物质炭含有的易氧化有机碳已经被微生物矿化分解，降低了微生物分解过程中不稳定碳的利用率；另一方面生物质炭对本底有机碳的吸附防止了微生物及其胞外酶的进入^[15]，从而降低生物质炭对本底有机碳的分解。Zhou等^[35]研究发现，施用生物质炭对亚热带森林土壤的自养呼吸和异养呼吸没有显著影响，且生物质炭的不同施用量对土壤自养呼吸和异养呼吸也没有显著差异，这与本研究结果不一致。可能是由于本研究绿洲农田风沙土经历了长期的耕作模式，地表微环境以及土壤有机质的变化都能够改变土壤异养呼吸，同时土地管理措施以及水热因子的变化间接影响土壤自养呼吸^[36]。由此可见，生物质炭对土壤呼吸组分的影响，可能与土壤类型、作物种类、生物质炭的类型^[37]、施用量以及试验时间有一定关系。

3.2 施用生物质炭5年后土壤呼吸及水热因子的关系

土壤温度和含水率是影响土壤呼吸的关键因子。土壤温度的变化在很大程度上可以解释土壤呼吸日变化^[38]。本研究采用线性方程对土壤呼吸和土层5、10、15 cm温度进行拟合，发现土壤呼吸与土壤温度显著相关，这与张阳阳等^[39]研究结果相似。温度会影响土壤酶与微生物活性以及植物根系生长^[40]，因此与土壤微生物呼吸和植物根系呼吸密切相关。本试验测得的不同处理土壤Q₁₀介于1.16~1.23之间，这与Zheng等^[41]测定的农田生态系统Q₁₀范围（1.28~4.75）不同，可能是由于生物质炭进入土壤5年后，对碳排放的激发效应随着土壤中易氧化有机碳的矿化逐渐消失^[20]。同时，也有研究发现，不同环境下的Q₁₀值是变化的，在温度较低的冬春季Q₁₀值较高，而在夏季Q₁₀值较低^[42]。

水分对土壤呼吸的影响具有复杂性与不确定性^[43]，土壤含水率主要受降雨和农田灌溉的影响。本研究发现，不同处理土壤呼吸与土壤含水率之间无显著相关性，这可能是由于一方面干旱区降雨量小，风沙土保水性能差；另一方面覆膜滴灌的灌溉方式导致土壤含水率变化范围较小，不足以影响微生物活性以及群落结构^[44]，故未能体现出土壤含水量对土壤呼吸速率的作用，而土壤温度则是影响土壤呼吸的主要环境因子。

4 结论

（1）风沙土土壤呼吸与异养呼吸日动态均呈单峰曲线，峰值出现在13：00左右；施用126.00 t·hm^-2生物质炭处理日平均呼吸速率为2.74 μmol·m^-2·s^-1，显著高于其他处理（P＜0.05）；生物质炭显著改变了土壤异养呼吸，对自养呼吸无显著影响。

（2）土壤呼吸速率随生物质炭施用量的增加显著增加，而随玉米生育时期的推进逐渐降低。

（3）不同生物质炭处理土壤呼吸速率与土层5、10、15 cm温度显著相关，土壤呼吸温度敏感系数随生物质炭施用量的增加而提升，而土壤呼吸与土壤含水率无显著相关性。