长期定位试验区始建于1979年, 位于哈尔滨市道外区民主镇黑龙江现代农业示范区(东经126° 51', 北纬45° 50')。 该区域地势平坦, 海拔151 m, 属松花江二级阶地, 成土母质为洪积黄土状粘土, 黑土层厚度为50 cm。 气候属中温带大陆性季风气候, 冬季寒冷干燥, 夏季炎热多雨, 年均气温3.5 ℃, 昼夜温差较大, 年均无霜期135 d, 年均降雨量533 mm。 1979年定位试验种植小麦, 采用小麦-大豆-玉米轮作制, 每年一季, 设置24个施肥处理, 每个处理3次重复, 无灌溉设施。 初始耕层(0~20 cm)土壤有机碳含量15.5 g· kg^-1, 全氮含量1.47 g· kg^-1, 全磷含量1.07 g· kg^-1, 全钾含量25.16 g· kg^-1, 有效氮含量151 mg· kg^-1, 速效磷含量51 mg· kg^-1, 速效钾含量200 mg· kg^-1, pH值7.2。

试验选择2017年10月收获后长期定位试验中不施肥(CK)、 常量氮磷钾(NPK)、 常量有机肥+常量氮磷钾(MNPK)、 常量有机肥(M)等4个处理。 所施用化肥为尿素(N 46%)、 重过磷酸钙(P₂O₅ 46%)、 磷酸二铵(N 18%、 P₂O₅ 46%)和硫酸钾(K₂O 50%), 化肥均为秋收后施用。 有机肥为发酵后纯马粪, 每次施用前测定养分和含水率, 常量施用水平按氮量75 kg N· hm^-2, 每轮作周期玉米季收获后秋施。 常量氮磷钾肥料水平, 在轮作周期的小麦季和玉米季为N 150 kg· hm^-2, P₂O₅ 75 kg· hm^-2和K₂O 75 kg· hm^-2, 大豆季为N 75 kg· hm^-2、 P₂O₅ 150 kg· hm^-2和K₂O 75 kg· hm^-2。

在2015年秋季玉米收获后施用纯马粪, 2016年和2017年分别种植作物小麦和大豆, 选择2017年夏季大豆开花结荚期进行田间的根际和非根际土壤样品取样, 此时大豆处于生殖生长期, 根际效应最为强烈。 在各处理小区随机选10株植株, 收获地上部, 挖取整个根系带回实验室后, 轻柔土壤将松散土壤抖掉去除, 采用毛刷收集附着在根系表面的土壤为根际土壤, 对应选取行间土壤为非根际土壤。 过2 mm筛混匀样品, 用镊子去除根系、 植物残体以及其他杂质, 土壤样品装入聚乙烯塑料袋中。 土样自然风干用于土壤基础理化性状、 有机碳结构测定。

(1)土壤有机碳测定, 采用重铬酸钾外加热法测定土壤有机碳含量。

(2)不同活性有机碳测定

不同活性有机碳测定采用高锰酸钾溶液氧化-比色法。 33 mmol· L^-1 KMnO₄测定值为高活性有机碳含量, 33和167 mmol· L^-1测定值的差值为中活性有机碳含量, 167和333 mmol· L^-1测定值的差值为低活性有机碳, 低、 中、 高活性碳比例是指分别占根际(非根际)土壤总有机碳的比例。

(3)核磁共振测定土壤碳结构

本试验土壤样品使用Bruker Avance 400固体核磁共振仪器。 应用¹³C CP/TOSS技术来辨别有机质特定的官能团。 本试验采用4 mm双共振魔角旋转(MAS)探头, ¹³C-NMR的频率为100 MHz。 应用¹³C CP/TOSS试验转速为5 kHz, 选好延迟0.8 s, ¹H 90° 脉宽4 μ s, 交叉极化接触时间为1 ms。 在检测之前, 应用四脉冲总旋转边带抑制和双脉冲相位调制解耦优化分辨率。 本研究中所有NMR图谱中每一官能团在全谱中所占比例均通过积分获取, 其中甲氧基积分后所占比例需再除以0.65, 用以校正因偶极去相时所引起信号的损失。 波谱被分为7个化学位移区域^{[10, 11]}(见表1), 烷基碳(δ =0~45)/烷氧基碳(δ =45~110)和芳香碳(δ =110~160)/总碳(δ =0~160)比值, 可以表征土壤有机质的分解程度和有机质质量^{[12, 13]}, 脂肪碳/芳香碳和疏水性碳/亲水性碳分别表示土壤有机碳的脂化度及疏水性, 其中脂肪碳/芳香碳=[烷基碳(δ =0~45)+甲氧基碳(δ =45~60)+碳水化合物类碳(δ =60~94)+双氧烷基碳(δ =94~110)]/[芳香碳(δ =110~142)+酚芳基碳(δ =142~160)], 疏水性碳/亲水性碳=[烷基碳(δ =0~45)+芳香碳(δ =110~142)+酚芳基碳(δ =142~160)]/[甲氧基碳(δ =45~60)+碳水化合物类碳(δ =60~94)+双氧烷基碳(δ =94~110)+羰基碳(δ =160~220)]^[14]。

表1

Table 1

表1(Table 1)

表1 13C NMR CP/TOSS分析波谱中化学位移及其对应成分 Table 1 Signal assignment in the 13C NMR CP/TOSS spectra of soil 化学位移 对应官能团(大分子化合物) 文献 0~45 烷基碳(脂肪类) [10, 11] 45~60 甲氧基碳(木质素和氨基酸) [11, 13] 60~94 碳水化合物类碳(纤维素、 半纤维素和木质素) [11, 13] 94~110 双氧烷基碳(纤维素和半纤维素) [11, 13] 110~142 芳香基碳(单宁和木质素) [10, 11] 142~160 O-芳香基碳(多酚、 单宁和木质素) [10, 11] 160~220 羰基碳(酰胺、 羰基酸和酯类化合物) [10, 11]

表1 13C NMR CP/TOSS分析波谱中化学位移及其对应成分 Table 1 Signal assignment in the 13C NMR CP/TOSS spectra of soil

统计分析采用WPS 2016软件对试验数据进行处理和作图, 采用SPSS 19软件对数据进行方差分析(ANOVA), 不同处理之间数据多重比较采用Duncan新复极差法检验(p< 0.05)。

图1显示, 施肥(M, NPK和MNPK)处理大豆季根际土壤有机碳含量显著高于不施肥处理, 且施肥处理间根际土壤有机碳含量差异不明显, 大小为MNPK> NPK> M。 在非根际土壤中, 不同施肥处理有机碳含量变化规律与根际土壤一致, 以CK处理有机碳含量最低(14.75 g· kg^-1), 施肥处理显著增加了非根际土壤有机碳含量。

图1

Fig.1

Figure Option ViewDownloadNew Window

图1 黑土长期不同施肥处理对根际和非根际土壤有机碳影响 注: 数据柱上标识的不同小写字母表示根际土壤之间差异显著, 数据柱上标识的不同大写字母表示非根际土壤之间差异显著(p< 0.05); * 表示根际与非根际土壤之间差异显著(p< 0.05), 下同Fig.1 Effects of long-term different fertilization on SOC in rhizosphere and bulk soils Note: Different lowercase letters above the bars indicate significant difference among rhizosphere soils, and different capital letters above the barsindicatesignificant difference among bulk soils at the 0.05 level; * indicate significant difference between rhizosphere and bulk soils at 0.05 level, respectively. The same below

所有处理中根际土壤有机碳含量均高于非根际, CK和M处理的根际土壤有机碳比非根际分别高1.2和1.1 g· kg^-1, 均未达到显著水平, 而NPK和MNPK处理均达到显著水平, MNPK处理的根际效应最显著, 比非根际增加了14.9%。

从表2可知, 在根际土壤中不同处理的低活性有机碳含量为MNPK> NPK> CK> M, 但处理间差异不明显。 施肥处理显著影响非根际土壤低活性有机碳含量, 大小为MNPK> CK> NPK> M, 其中MNPK处理的非根际土壤低活性有机碳含量显著大于M处理。 MNPK和M处理的根际土壤低活性有机碳高于非根际土壤, 而CK和NPK处理根际土壤低活性有机碳含量低于非根际, 但均未达到显著水平(表2)。

表2

Table 2

表2(Table 2)

表2 黑土长期不同施肥处理对根际和非根际土壤活性有机碳影响 Table 2 Effects of long-term different fertilization on labile organic carbon in rhizosphere and bulk soils 处理 低活性有机碳 /(g· kg-1) 中活性有机碳 /(g· kg-1) 高活性有机碳 /(g· kg-1) 根际 CK 3.94a 3.84a 0.56a M 3.91a 3.68a 0.52a NPK 3.98a 3.82a 0.54a MNPK 4.46a 3.90a 0.54a 非根际 CK 4.18AB 3.70A 0.59A M 3.88B 3.69A 0.59A* NPK 4.09AB 3.64A 0.59A* MNPK 4.21A 3.59A 0.60A* 注: * 表示根际与非根际土壤之间差异性显著(p< 0.05), 下同 Note: * Indicate significant difference between rhizosphere and bulk soils at 0.05 level, The same below

表2 黑土长期不同施肥处理对根际和非根际土壤活性有机碳影响 Table 2 Effects of long-term different fertilization on labile organic carbon in rhizosphere and bulk soils

不同施肥处理对大豆季土壤中、 高活性有机碳含量影响不明显, 均未达到显著水平。 所有施肥处理的根际土壤中活性有机碳含量均高于非根际, 而根际土壤高活性有机碳含量均低于非根际, 其中M, NPK和MNPK处理达到显著水平。

从表3可知, 土壤活性有机碳仅占总有机碳含量的42.9%~57.7%, 在活性有机碳中又以低活性有机碳为主, 占有机碳比例为18.9%~30.3%, 其次是中活性有机碳(16.9%~29.8%), 最低是高活性有机碳(2.67%~3.91%)。 根际土壤低活性有机碳/总有机碳比值在不同施肥间差异不明显。 在非根际土壤中, CK处理的非根际土壤低活性有机碳/总有机碳比值最高, 与M和NPK处理相比均达到显著水平。 施肥处理均表现为根际土壤低活性有机碳/总有机碳比值低于非根际的趋势。

表3

Table 3

表3(Table 3)

表3 黑土长期不同施肥处理对根际和非根际土壤活性有机碳/总有机碳含量比值影响 Table 3 Effects of long-term different fertilization on labile organic carbon/organic carbon ratio in rhizosphere and bulk soils 处理 低活性 有机碳/ 总有机碳 比值 中活性 有机碳/ 总有机碳 比值 高活性 有机碳/ 总有机碳 比值 根际 CK 0.25a 0.24a 0.035a M 0.21a 0.19b 0.028b NPK 0.22a 0.21ab 0.030b MNPK 0.23a 0.20ab 0.028b 非根际 CK 0.28A 0.25A 0.04A M 0.24B* 0.22B 0.036B* NPK 0.24B* 0.22B 0.035B* MNPK 0.26AB* 0.22B 0.037AB*

表3 黑土长期不同施肥处理对根际和非根际土壤活性有机碳/总有机碳含量比值影响 Table 3 Effects of long-term different fertilization on labile organic carbon/organic carbon ratio in rhizosphere and bulk soils

CK处理的根际土壤中活性有机碳/总有机碳比值最高。 与CK处理相比, M处理显著降低根际土壤中活性有机碳/总有机碳比值。 在非根际土壤中, CK处理的中活性有机碳/总有机碳比值也最高, 均显著高于M, NPK和MNPK处理。 无论在根际还是非根际土壤中, 施肥处理(M, NPK和MNPK)间中活性有机碳/总有机碳比值差异不显著(表3)。

在根际土壤中, 与CK处理相比, 其他处理均显著降低高活性有机碳/总有机碳含量比值, M, NPK和MNPK处理间差异不明显。 非根际土壤中, CK处理的高活性有机碳/总有机碳含量比值显著高于M和NPK处理, 施肥(M, NPK和MNPK)处理间差异不显著。 所有处理中根际高活性有机碳/总有机碳含量比值均显著低于非根际(表3)。

从图2可知, 土壤有机碳NMR谱各峰值差异比较明显, 在δ =0~45烷基碳范围内, 有两个主要吸收峰为δ =21~23和δ =30~31, 分别表示与芳香环相连的CH₃和脂肪族多肽链中无定型(CH₂)_n; δ =59~60代表甲氧基碳(— OCH₃)吸收峰, 与木质素及木质素替代物有关; δ =71~72表示碳水化合物中仲醇(— CR₂— OH)吸收峰; δ =103~105为双氧烷基碳(— C— Di— O)吸收峰; δ =45~100范围内的碳主要来自纤维素、 半纤维素及木质素的烷氧基碳; δ =127~130表示被羧基(— COOH)或者羧甲基(— CH₂COO— )取代的芳香基碳; δ =150~155表示与O, N等取代基相连的芳香基碳; δ =110~160范围峰主要来自多酚、 软木质、 单宁和木质素等的芳香碳; δ =173~174代表羧酸、 酞胺COO/N— C=O和酯类吸收峰。 从NMR图谱官能团相对含量可知(表4), 碳水化合物类碳所占比例最大, 其次是烷基碳、 芳香基碳及羰基碳, 其中烷氧基碳(δ =45~110)占到38.3%~43.9%。

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 不同处理的土壤有机碳化学组成的NMR谱图Fig.2 13C NMR of peat profile samples in rhizosphere and bulk soils under long-term different fertilization

表4

Table 4

表4(Table 4)

表4 不同施肥下根际与非根际土壤不同官能团相对含量 Table 4 Assignments at different chemical shift regions and their relative proportions in rhizosphere and bulk soils under long-term different fertilization Treatments (Alk C 0~45 Met C 45~60 Cad C 60~94 Di-O-alk C 94~110 O-alk C 45~110 Ary C 110~142 O-Ary C 142~160 Aro C 110~160 Car C)% 160~220 Alk C/ O-alk C Aro C/ Tot C Ali C/ Aro C Hb C/ Hl C RS CK 18.15 8.86 27.96 6.03 42.85 17.53 4.59 22.12 16.89 0.42 0.34 2.76 0.67 M 20.41 9.09 27.11 5.85 42.05 17.22 5.50 22.72 14.82 0.49 0.35 2.75 0.76 NPK 20.93 8.60 26.71 6.04 41.35 17.12 6.02 23.14 14.58 0.51 0.36 2.69 0.79 MNPK 20.78 9.21 28.38 5.96 43.55 15.40 4.73 20.13 15.53 0.48 0.32 3.20 0.69 BS CK 19.54 8.30 25.32 5.52 39.14 19.54 5.88 25.42 15.92 0.50 0.39 2.31 0.82 M 18.24 8.02 23.53 5.14 36.69 22.33 6.94 29.27 15.79 0.50 0.44 1.88 0.91 NPK 19.38 7.65 25.93 5.70 39.28 19.82 6.04 25.86 15.48 0.49 0.40 2.27 0.83 MNPK 22.14 9.40 26.67 6.06 42.13 18.10 5.29 23.39 12.35 0.53 0.36 2.75 0.84 Note: Alk C: Alkyl C, Met C: Methoxyl C, Cad C: Carbohydrate C, Di-O-alk C: Di-O-alkyl C, O-alk C: O-alkyl C, AryC: Aryl C, O-Ary C: O-ArylC, Aro C: Aromatic C, Car C: Carbonyl C, Alk C/O-alk C: Alkyl C/O-alkyl C, Aro C/Tot C: Aromatic C/Total C, Ali C/Aro C: Aliphatic C/Aromatic C, Hb C/Hl C: Hydrophobic C/Hydrophilic C

表4 不同施肥下根际与非根际土壤不同官能团相对含量 Table 4 Assignments at different chemical shift regions and their relative proportions in rhizosphere and bulk soils under long-term different fertilization

在根际土壤中, 不同施肥处理能够明显影响土壤有机碳不同官能团相对含量, 与CK处理相比, 施肥处理增加烷基碳、 烷基碳/烷氧基碳和疏水性碳/亲水性碳比值, 施肥处理的烷基碳相对含量大小为NPK> MNPK> M, 烷基碳/烷氧基碳比值为NPK> M> MNPK, 疏水性碳/亲水性碳比值为NPK> M> MNPK。 与CK处理相比, M处理降低烷氧基碳、 芳香基碳和羰基碳。 与NPK处理相比, 增加有机肥后MNPK处理升高甲氧基碳、 碳水化合物类碳和脂肪碳/芳香碳比值, 降低烷基碳、 芳香基碳、 烷基碳/烷氧基碳、 芳香碳/总碳和疏水性碳/亲水性碳比值(表4)。

非根际土壤中, MNPK处理相比于CK处理增加了烷基碳、 烷氧基碳、 烷基碳/烷氧基碳比值和脂肪碳/芳香碳比值相对含量。 与CK处理相比, M处理升高芳香基碳、 O-芳香基碳、 芳香碳/总碳和疏水性碳/亲水性碳比值。 与NPK处理相比, MNPK处理的烷基碳、 甲氧基碳、 碳水化物类碳、 双氧烷基碳、 烷基碳/烷氧基碳、 脂肪碳/芳香碳比值增加, 其他官能团降低或者持平。

不同施肥处理的根际与非根际土壤¹³C-NMR图谱中吸收峰相似, 只是吸收峰的相对强度存在明显差异, 表明不同施肥处理未能改变根际与非根际土壤中有机碳的碳骨架, 仅是改变了有机碳中官能团的数量。 与CK处理的非根际土壤相比, 根际土壤增加甲氧基碳、 碳水化物类碳、 双氧烷基碳、 羰基碳、 脂肪碳/芳香碳比值, 而降低其他官能团相对含量及烷基碳/烷氧基碳、 芳香碳/总碳和疏水性碳/亲水性碳比值。 与根际土壤中CK处理相比, 单纯施用M和NPK处理能够增加根际土壤有机碳中烷基碳相对含量, NPK处理降低羰基碳及增加烷基碳/烷氧基碳比值, M处理降低根际土壤羰基碳; 而MNPK处理调节根际土壤有机碳中甲氧基碳和双氧烷基碳, 导致其相对含量降低。

由于植物向根系转运的沉积碳和根系脱落物在根际微域的持续输入, 导致根际微生物大量繁殖, 土壤有机碳转化快, 与非根际土壤存在明显差异。 有研究显示, 土壤碳素的内循环主要发生在植物根际微域中, 在土壤、 微生物和动植物体内不断快速迁移和转化。 本研究发现, 不论是否施肥, 大豆季根际土壤有机碳含量均高于非根际土壤, 长期施肥处理能够显著增加外源碳输入量, 增加根际与非根际土壤中总有机碳和低活性有机碳的累积, 其中以NPK配施有机肥处理提升效果最好。 ¹³C-核磁共振结果表明, 相比于CK处理, NPK配施有机肥处理明显增加根际土壤烷基碳(δ =0~45)、 烷氧基碳(δ =45~110)比例以及烷基碳/烷氧基碳比值, 降低芳香基碳(δ =110~160)和芳香碳/总碳比值, 在非根际土壤中尤其显著。 研究证明烷氧基碳(δ =45~110)升高能够改善土壤物理结构、 腐殖酸合成以及生物活性, 促进土壤团聚体形成和增加土壤颗粒结构稳定性^[15], 烷氧基碳的转化效率与土壤粘粒含量呈正相关^[16], 粘粒又促进团聚体的形成和稳定, 降低微生物对团聚体中碳素的降解, 同时, 烷基碳比例、 烷基碳/烷氧基碳比值可以用来显示有机碳的分解程度以及稳定性, 烷基碳比例和烷基碳/烷氧基碳比值越大, 表明有机质分解越彻底且有机碳越稳定^{[12, 13]}, 有机无机配施通过改变土壤有机碳组分比例和结构起到固碳的作用。 与CK处理相比, 单施化肥(NPK)处理无论在根际和非根际土壤均增加芳香基碳(δ =110~160)比例和芳香碳/总碳比值, 而烷基碳(δ =0~45)比例、 烷氧基碳(δ =45~110)和烷基碳/烷氧基碳比值在根际和非根际土壤存在差异, 在根际土壤烷基碳(δ =0~45)比例和烷基碳/烷氧基碳比值增加, 烷氧基碳(δ =45~110)比例降低, 非根际土壤中结果相反。 表明NPK处理能够增加根际有机碳难降解组分烷基碳(δ =0~45)和芳香基碳(δ =110~160)的比例, 降低根际土壤的团聚化, 增加根际土壤微生物对有机碳的降解。