文献参数: 张爽, 杨佳, 刘清华, 等.小黑杨与紫穗槐人工林对退化典型黑土土壤呼吸的影响[J].水土保持通报,2026,46(4):286-296. Citation:Zhang Shuang, Yang Jia, Liu Qinghua, et al. Effects of Populus × xiaohei and Amorpha fruticosa plantations on soil respiration in degraded typical black soils [J]. Bulletin of Soil and Water Conservation,2026,46(4):286-296.
中国东北黑土区是全球四大黑土区之一,其土壤有机碳储量占全国耕地总量的25%以上,是中国重要的商品粮生产基地,对于保障中国乃至全球的粮食安全具有重要的战略意义
[1]。然而,由于长期高强度农业利用和不合理开垦,该区域正面临严重的土壤退化问题。研究表明,在过去30 a间,东北黑土区表层土壤有机碳含量下降了约20%~30%
[2],且黑土呈现出变薄、变硬、变瘦的综合退化特征。土壤退化不仅会威胁粮食产量,还会导致土壤结构恶化和生物活性降低,深刻影响陆地生态系统的碳循环过程
[3]。土壤呼吸是土壤碳库向大气输出的主要途径,是陆地生态系统碳循环的关键环节,其速率是反映土壤生物活性强弱和碳循环的关键指标,对生态系统功能的恢复至关重要
[4]。土壤呼吸受植物根系活动、土壤微生物群落及有机质供给等生物因素,以及土壤温度、含水量等非生物因素的共同调控
[5]。在退化黑土的修复中,植被恢复是调控土壤碳循环的核心策略,在黑土修复和维持土壤碳平衡方面起着至关重要的作用。植被可以通过改变凋落物输入和根系活动,调节土壤有机碳的累积与分解过程,从而影响土壤呼吸和土壤碳库动态
[6-7]。同时,植被冠层及地表覆盖还可通过改变遮阴程度和土壤结构,间接影响土壤温度和含水量,对土壤碳排放过程产生进一步调节作用。
已有研究表明,土壤温度和土壤体积含水量可以解释土壤呼吸的大部分变异,其通过影响根系代谢活动和土壤微生物过程而对土壤呼吸产生重要作用。例如,张克胜的研究表明,土壤温度与土壤呼吸速率呈显著的指数关系
[8]。也有研究表明,土壤呼吸与土壤体积含水量的关系通常不显著,但土壤温度与土壤体积含水量的交互作用可以显著影响土壤呼吸
[9]。在紫穗槐纯林和北方杨树人工林中土壤呼吸随土壤温度也呈现明显的季节变化,土壤温度是解释其变异的主要因子之一,土壤含水量多表现为次要或限制性因子。因此,准确监测退化黑土土壤呼吸及其对温度、含水量的响应,对于评估修复措施对土壤碳循环过程的影响至关重要。
尽管已有研究关注人工林对土壤呼吸的影响,但在东北退化黑土区沿退化梯度,系统比较以小黑杨和紫穗槐为代表的乔木型与灌木型人工幼龄林对土壤呼吸及其温度响应的研究仍相对有限。因此,探究黑土区不同退化程度下人工幼龄林土壤呼吸及其影响机制,以评估人工幼龄林对退化黑土土壤碳循环过程和土壤质量改善的作用,揭示退化黑土的质量改善效果以及土壤碳循环过程具有重要的理论和实践意义。
小黑杨(
Populus ×
xiaohei)为东北黑土区主要防护林乔木树种,具有深根系和快速生长的特点,能显著改变土壤微环境;紫穗槐(
Amorpha fruticosa)为灌木,对土壤要求不高,可改善土壤结构,提高土壤肥力,还可通过根瘤固氮作用影响土壤氮循环
[10]。这种差异可能导致小黑杨与紫穗槐人工林在相同退化程度黑土中的土壤呼吸调控效应并不相同,且在不同退化程度黑土中对土壤呼吸的调控效应也有所不同。这限制了我们从土壤碳循环管理的角度,为不同退化程度黑土修复准确选择造林植物种的能力。因此,本研究以轻度、中度和重度退化黑土为对象,采用野外原位连续测定的方法,对比分析种植小黑杨和紫穗槐人工幼龄林黑土的土壤呼吸速率差异,重点探讨: ①对分别种植小黑杨和紫穗槐人工林的不同退化程度黑土土壤呼吸速率的影响; ②土壤温度和湿度对退化黑土土壤呼吸速率的影响; ③植被类型、土壤退化程度对土壤呼吸速率的潜在影响机制。研究结果旨在更好地理解乔木和灌木对退化黑土土壤呼吸、土壤碳排放与土壤碳循环过程的影响,从而为基于人工林的退化黑土碳管理及修复造林植物种的选择提供科学依据,为实现黑土的可持续利用提供理论和决策支撑。
1 材料与方法
1.1 研究区概况
研究区位于黑龙江省齐齐哈尔市克山农场(125°08′—125°37′E,48°12′—48°23′N)。克山农场为典型国有农场,自建场以来长期以机械化粮食作物种植为主,是松嫩平原重要商品粮基地之一。该农场地处讷河市与克山县交界处,位于松嫩平原腹地,地势平坦,略有起伏
[11]。克山农场的地貌特征为漫川漫岗,平均海拔315 m,坡长而缓,坡长多为200~1 000 m,坡度平均3°~5°,土壤类型以黏化湿润均腐土为主。农场的主要土地利用类型包括耕地和林地:耕地主要用于农作物种植,林地则涵盖了人工经济林和农田防护林等。该地区夏季高温多雨,冬季寒冷干燥。其中,年均降雨约为600 mm,且主要集中在6—8月;年平均气温约0.9 ℃。克山农场所属地区为典型黑土区,土质肥沃,黑钙土占土地面积的98%,适宜农作物生长。
1.2 试验地设计
2022年5月在克山农场科技园的一处农田开展试验,在距离农田地边5 m之处,选定3处长60 m,宽7 m的农田,基于 “分层剥离和交错回填”方法设置退化样地。根据土壤有机质含量差异模拟不同退化程度土壤
[12-13]。分别收集垂直剖面的3层土层,即0—20,20—40,40—60 cm进行回填。回填0—20 cm土壤作为轻度退化土壤(lightly degraded, LD),20—40 cm土壤为中度退化土壤(moderately degraded, MD),40—60 cm土壤为重度退化土壤(severely degraded, SD)。退化程度以模拟土壤剖面构建方式并结合 土壤有机质(SOM),土壤有机碳(SOC),全氮(TN)等理化指标的显著差异来界定(
表1),从而形成由 LD向 SD逐渐加重的退化梯度。2023年5月,在每个退化黑土样地上设置大小相同(4 m×7 m)的区域分别实施小黑杨(两根一干苗,树高约为1.3 m)的种植和紫穗槐的行播,每个小区种植7行,3列,共91株苗,种植行、列距均为0.5 m,并设置与试验样地大小相同的空白对照样地(3种退化程度×3种处理×3个重复=27块小区),均为4 m×7 m。试验样地的基本信息见
表1。
1.3 土壤样品采集及理化性质测定
于植物生长旺季,2024年9月10日,在种植小黑杨和紫穗槐的各小区以及空白对照小区内,使用铲子在每个小区沿对角线按5点取样法采集0—20 cm表层土壤,将5个点的土样混合为一个混合样,依次装入密封袋并做好标记带回实验室。土壤样品经过自然风干后,进行土壤理化性质测定。土壤pH值采用电位法测定;土壤容重采用环刀法测定;土壤有机碳含量采用重铬酸钾法测定;土壤碱解氮通过碱解扩散法测定;土壤全碳通过碳氮分析仪测定;土壤全氮质量分数通过凯氏法测定;土壤全磷质量分数采用硫酸-高氯酸消解钼锑抗比色法测定;土壤全钾含量测定采用原子吸收分光光度法;土壤速效磷采用氟化铵浸提-钼锑抗比色法测定;土壤速效钾采用乙酸铵浸提-原子吸收仪测定
[14]。
1.4 土壤呼吸速率测定
首先在每个种植小区和空白小区内设置土壤呼吸环,具体为在每个小区的对角线交叉位置设置1个内径20 cm,高12 cm的PVC呼吸环作为固定样点,并在整个观测期内保持位置不变。将PVC呼吸环底端朝下垂直于地表插入土壤中约9 cm,保留环顶部距离土壤表面3 cm。随后,在2024年5月至10月分别对土壤呼吸速率进行测定。测定时间为每月10日,若遇下雨则根据降雨量大小延迟2~3 d测定土壤呼吸速率。在测量前一天,贴地面清除土壤呼吸环内的杂草。在测定日的7:00—11:00,使用便携式土壤碳通量观测系统(LI-870)对土壤呼吸进行测定(测定时间3 min,包括30 s的预吹风,120 s的观测时间和30 s的换气)。同时,使用仪器配备的温湿度传感器测量呼吸环附近5 cm深处的土壤温度(T)和土壤体积含水量(W)。
1.5 数据处理
首先建立土壤呼吸速率(以C计)与土壤温湿度之间的单变量模型。根据已有研究结果,使用
公式(1)拟合
[15]土壤呼吸速率与土壤5 cm土壤温度间的关系:
式中:Rs 表示土壤呼吸速率〔μmol/(m²·s),以C计〕; R0表示温度为0 ℃时的土壤呼吸速率〔μmol/(m²·s),以C计〕;T表示5 cm土壤温度(℃);a和b为拟合参数。
在
公式(1)的基础上,使用
公式(2)[16]推导土壤呼吸温度敏感性系数:
式中:Q10表示土壤呼吸温度敏感性系数; a, b为公式(1)中的拟合参数。
使用公式(3)探究土壤呼吸速率与土壤体积含水量间的关系:
式中:Rs 表示土壤呼吸速率〔μmol/(m²·s),以C计〕; W表示土壤体积含水量(cm3/cm3); a和b为拟合参数。
使用公式(4)拟合土壤呼吸速率与5 cm土壤温度和土壤体积含水量间的多元回归模型:
式中:Rs 表示土壤呼吸速率〔μmol/(m²·s)〕; T表示5 cm土壤温度(℃); W表示土壤体积含水量(cm3/cm3); a, b,c和d为拟合参数。
1.6 数据分析
使用多因素方差分析探究月份、退化程度、植被类型及其交互作用对土壤呼吸速率、5 cm土壤温度和土壤体积含水量的影响。采用单因素方差分析研究不同土地退化程度间平均土壤呼吸速率的差异,如果方差分析显著,采用Duncan法进行多重比较。采用一元非线性回归分析探究土壤呼吸速率和5 cm土壤温度间的关系;采用一元线性回归分析探究土壤呼吸速率和土壤体积含水量间的关系;采用多元回归分析探究土壤呼吸速率与5 cm土壤温度和土壤体积含水量间的协同关系。采用Pearson相关分析探究土壤性质对土壤呼吸速率的影响。采用RDA分析探究土壤性质对土壤呼吸速率及Q10的多元影响。采用偏最小二乘法路径模型(PLS-PM)探究植被类型、土壤因素对土壤呼吸速率的直接和间接影响。研究的数据分析均在R 4.0.1中完成,所有显著性水平均设为α=0.05。
2 结果与分析
2.1 不同退化程度黑土土壤呼吸速率、5 cm土壤温度和土壤体积含水量变化
从
表2中可以看出,月份、土壤退化程度、植被类型及其交互作用对土壤呼吸速率均有显著影响(
p<0.001)。相比之下,月份、土壤退化程度、植被类型、月份与土壤退化程度交互作用及土壤退化程度与植被类型交互作用对5 cm土壤温度有显著影响(
p<0.05);月份、植被类型及月份与植被类型交互作用对土壤体积含水量有显著影响(
p<0.001)。从
图1中可以看出,种植小黑杨和紫穗槐及对照组的黑土土壤呼吸速率、5 cm土壤温度和土壤体积含水量在不同退化程度黑土中均呈现明显的季节动态。其中,小黑杨和紫穗槐的土壤呼吸速率均在7,8月达到峰值,随后逐月下降。整体上,不同退化程度和植被措施下5 cm土壤温度均表现为先增加后降低的趋势,且均在7月份达到峰值。相比之下,种植小黑杨和紫穗槐的黑土土壤体积含水量在不同退化程度黑土中呈现先降低后增加再降低的趋势。在各退化土壤中,紫穗槐种植地的平均土壤呼吸速率2.31 µmol/(m²·s)(以CO
2计,下同)总体上高于小黑杨栽植地的平均土壤呼吸速率1.96 µmol/(m²·s)。从
图2中可以看出,小黑杨、紫穗槐及对照组的平均土壤呼吸速率均随着退化程度的增加逐渐降低。在LD土壤中(
图2),紫穗槐的平均土壤呼吸速率显著高于小黑杨和对照组(
p<0.05),且小黑杨和对照组间差异不显著。在MD和SD土壤中,小黑杨和紫穗槐平均土壤呼吸速率差异不显著,但均显著高于对照组(
p<0.05)。
2.2 土壤呼吸速率与5 cm土壤温度和土壤体积含水量的关系
从
图3和
表3中可以看出,土壤呼吸速率与5 cm土壤温度存在显著的指数关系。整体上,小黑杨、紫穗槐和对照组的土壤呼吸速率均随土壤温度的升高而增大。其中,MD下对照组的
R2最大,为0.656,SD下对照组的
R2最低,为0.143。由此可见,MD下对照组的土壤温度对土壤呼吸速率的贡献最大,SD下对照组的贡献最小。此外,不同处理下,土壤呼吸的
Q10值范围为1.162~1.804。其中,种植紫穗槐的MD土壤呼吸速率对土壤温度最敏感(
Q10为1.804)。
相比之下,小黑杨、紫穗槐和对照组的土壤呼吸速率与土壤体积含水量的关系均不显著。从
表4中可以看出,土壤呼吸速率与5 cm土壤温度和体积含水率的双因素复合模型显著提高了模型的预测精度。其中MD下小黑杨的R
2最大,为0.863,SD下紫穗槐的
R2最低,为0.314。
2.3 土壤呼吸与各因子的关系
采用Pearson相关分析探究土壤理化性质与土壤呼吸速率(
Rs )的关系(
表5)。结果表明,
Rs 与5 cm土壤温度、碱解氮呈极显著正相关(
p<0.01),与全钾、容重呈极显著负相关(
p<0.01)。研究选择土壤理化性质(全碳、有机质、全氮、速效磷、全磷、碱解氮、5 cm土壤温度、体积含水量、全钾和速效钾)和土壤呼吸速率及
Q10进行冗余分析(RDA)(
图4)。结果表明,RDA第一轴和第二轴对土壤呼吸速率及
Q10总方差的解释度分别为43.17%和12.20%。其中,土壤呼吸速率与全氮、速效磷、全磷、碱解氮、5 cm土壤温度、体积含水量显著正相关;
Q10与全碳、有机质、全钾、5 cm土壤温度、体积含水量显著正相关。偏最小二乘法路径模型(见
图5),结果表明,植被类型、土壤退化程度和5 cm土壤温度均可以显著直接影响土壤呼吸。此外,植被类型和土壤退化程度可以通过改变土壤养分间接影响土壤呼吸。
3 讨 论
3.1 不同退化程度黑土土壤呼吸的变化规律
种植小黑杨和紫穗槐人工林土壤呼吸都具有明显的月变化规律,且最大值均出现在7月和8月。在同一退化程度下,种植小黑杨和紫穗槐人工林的峰值月份一致,但季节变化幅度和平均水平存在差异。这一结果与王旭等
[17]观测的山杨白桦次生林、原始阔叶红松林土壤呼吸季节变化规律基本一致。这可能是因为适中的土壤温度加快了微生物活动与根系代谢,进而导致呼吸速率骤增。相比之下,9月后呼吸速率下降不仅与温度降低有关,还可能与植物生长后期光合产物分配减少,根系分泌物输入下降相关
[18]。土壤呼吸包括土壤中植物根系、微生物和土壤动物等组分的呼吸,其中,根系呼吸和微生物呼吸占据主导地位
[19]。研究中种植小黑杨和紫穗槐的不同退化程度黑土土壤呼吸速率均随退化程度加剧逐渐降低,表现为:LD>MD>SD。这可能是因为LD土壤中保留了较高的活性碳库和相对较好的物理结构,支持旺盛的微生物活动。
而MD和SD土壤中由于有机质含量的显著降低和容重的增加可能限制了底物可用性与氧气扩散,导致呼吸速率降低。种植小黑杨与紫穗槐后,LD土壤呼吸速率较对照提升40%~60%,但SD土壤呼吸速率仅提升10%~15%,这表明植被修复对轻度退化土壤的碳循环激活效应更显著。整体上,种植紫穗槐的土壤的呼吸速率高于种植小黑杨的呼吸速率,这可能是因为其土壤高氮含量(C/N≈13∶1)加速了微生物矿化过程
[20],且紫穗槐种植地凋落物较多,凋落物输入增加了土壤中的养分供给,促进了土壤生物对养分吸收的有效性,进而刺激土壤微生物的活动,提高了土壤的呼吸速率
[21]。而小黑杨土壤高木质素特性(C/N≈22∶1)可能抑制了短期碳释放。
3.2 5 cm土壤温度和体积含水量对土壤呼吸速率的影响
研究结果表明,不同退化程度黑土土壤呼吸速率与土壤温度均呈现指数正相关(
图3,
表3),土壤温度可以解释土壤呼吸变异的14.3%~65.6%,表明土壤温度是影响土壤呼吸的重要因素,这一结果得到了大量研究的支持
[22]。值得注意的是,研究中温度对种植小黑杨和紫穗槐的人工幼龄林不同退化程度黑土土壤呼吸的影响存在差异。从温度指数模型拟合结果看(
表3),小黑杨人工幼龄林在LD, MD和SD土壤中的
R²分别为0.643,0.554,0.497,整体高于紫穗槐人工幼龄林对应的0.439,0.432,0.217。这说明在各退化程度下小黑杨人工幼龄林土壤呼吸对温度变化的响应较为稳定,而紫穗槐幼龄在重度退化黑土上的土壤呼吸对温度的解释力明显较弱。此外,紫穗槐幼龄在MD土壤中的
Q10值达到1.804,高于小黑杨的1.448,这说明中度退化土壤中紫穗槐人工林土壤呼吸对温度升高的敏感性更强。LD土壤呼吸与土壤温度的相关性比SD土壤呼吸与土壤温度的相关性更大。相比之下,土壤体积含水量对土壤呼吸的影响较为复杂。目前,对于土壤含水量与土壤呼吸的关系得出的结论存在很大的差异。研究结果中土壤含水量与土壤呼吸速率的总体拟合效果较差,对土壤呼吸速率的影响不显著,这与郝龙飞等
[23]在东北地区白桦次生林及水曲柳、红松、长白落叶松人工林、全权等
[24]在辽东栎林、华北落叶松林和油松林的结果一致,上述研究均表明土壤温度是土壤呼吸变异的主要控制因子,而土壤水分的影响相对较小。相比之下,闫宝龙等
[25]以大兴安岭南段蒙古栎林、林缘草地、华北落叶松人工林、山杏灌丛和山荆子林5种植被类型为对象,发现土壤呼吸与土壤体积含水量和温度均呈显著正相关,表明在该区域水分是限制土壤呼吸的重要环境因子之一。这可能是因为不同地区生长环境不一样。闫宝龙研究区域的土壤具有良好的持水性,土壤含水量位于微生物活性敏感区间,而本研究的退化黑土土壤孔隙度低导致有效持水能力下降,且由于研究期间土壤含水量变化较小,水分不是植物及微生物活动的胁迫因子
[16],土壤含水量对土壤呼吸速率的影响会被其他因素的影响所掩盖。目前也有多数研究者认为单独研究土壤含水量与土壤呼吸的关系很困难,但土壤含水量与土壤温度之间有交互作用,并且他们协同作用对土壤呼吸的影响大于单独对土壤呼吸的影响 。本研究显示,土壤呼吸速率与5 cm土壤温度和土壤含水量的交互作用显著相关(
p<0.05)(
表4),相较于土壤温度与土壤呼吸速率间的单因子拟合模型,加入土壤体积含水量后的多元回归模型对土壤呼吸速率解释能力有一定提高,达到31.4%~86.3%,这与马涛等
[26]研究结果一致,土壤呼吸变化都不能很好地用土壤体积含水量来解释,但是土壤体积含水量可以增强土壤温度对土壤呼吸的影响。
Q10值的范围在1.162~1.804之间,这与谭炯锐等
[27]在北京人工杨树林中得到的
Q10值1.63接近,小黑杨和紫穗槐在SD中
Q10值最低,原因是
Q10值大小取决于有机碳的含量和组分,SD土壤中有机碳含量最低,进而导致
Q10值降低。在各退化黑土土壤中紫穗槐的
Q10值均高于小黑杨,唐夫凯等
[28]发现,生长季乔木的
Q10值高于灌木。这与本研究结果不同,这可能是因为紫穗槐属豆科固氮类植物,其根部共生的根瘤菌可以更好地将大气中的氮固定于土壤中,进而为微生物提供优质底物,增强微生物活性,促进易分解碳代谢,进而导致种植紫穗槐黑土土壤的
Q10值显著升高。
3.3 影响土壤呼吸速率因素的相关关系
偏最小二乘法路径模型结果表明,植被类型、土壤退化程度和土壤5 cm处温度不仅对土壤呼吸速率具有直接影响,植被类型和土壤退化程度还通过改变土壤养分产生显著的间接影响。具体而言,Pearson相关性分析(
表5)及RDA(
图4)的结果共同表明,土壤呼吸速率与土壤温度和碱解氮呈现极显著的正相关关系(
p<0.01),而与全钾、容重呈现极显著负相关关系(
p<0.01)。表明土壤养分有效性是调控呼吸过程的关键中间变量。紫穗槐作为一种固氮植物,可以提高土壤氮素含量,为微生物提供易分解底物,进而增强了微生物活性和根系代谢。相关性分析中碱解氮与呼吸速率的极显著正相关(
p<0.01),这可能是本研究中紫穗槐土壤呼吸速率显著高于小黑杨和对照组的原因。这一结果与郭明的研究一致,表明高氮环境可促进碳的快速矿化
[29],这与前文推测一致。相比之下,土壤退化程度通过改变土壤结构和土壤养分间接影响呼吸速率。土壤退化导致土壤团聚体破坏,孔隙度降低,容重与呼吸速率的显著负相关进一步从物理结构上解释了退化对呼吸的抑制机理,同时造成土壤有机质和养分含量下降,进而使得土壤呼吸速率随土壤退化程度增加均显著降低。这与赵慧丽等
[30]研究的结果一致,均强调退化程度抑制了土壤呼吸。最后需要客观指出的是,本研究仍存在一定的局限性。需在未来的工作中加以拓展。由于当前观测期尚处于植被恢复早期(小黑杨为两年生扦插苗,紫穗槐为种子直播形成的幼龄灌丛),本研究结论主要反映的是人工幼龄林对退化黑土土壤呼吸的短期调控响应。本阶段未能同步获取精细的林分生物量、冠层结构和根系分布等群落特征指标。因此现有结果实质上评估的是这两种造林物种作为整体群落在退化初期的综合宏观效应,长期的群落演替过程及其对土壤碳循环的深远影响仍需后续持续跟踪。未来可在自然退化黑土背景下,结合更全面的退化诊断指标体系,增加群落和根系特征、土壤水热养分及微生物的长期联合监测,并通过增加呼吸测定点和引入控制模拟试验,进一步验证和拓展本研究的结论。
4 结 论
(1) 植被措施(采用小黑杨和紫穗槐修复黑土)均能显著提升退化土壤的呼吸速率,但这种提升效果随退化程度加剧而显著减弱。
(2) 土壤呼吸呈现明显的季节动态,峰值出现在7—8月高温期,且其变化主要受5 cm土壤温度驱动;土壤体积含水量未对土壤呼吸速率表现出显著影响,土壤含水量通过与温度的协同作用共同影响土壤呼吸。
(3) Pearson相关分析、RDA分析以及偏最小二乘法路径模型表明,植被类型、退化程度和5 cm土壤温度不仅直接影响土壤呼吸,还通过影响土壤养分(如碱解氮、全氮、速效磷)产生显著的间接影响。
(4) 研究从碳循环管理的角度揭示了小黑杨和紫穗槐在退化黑土修复中的差异作用机制,在退化黑土的植被修复实践中,轻度退化区域应优先选择种植紫穗槐等具有固氮能力的灌木,以更有效地激活土壤碳循环。
国家自然科学基金项目“小黑杨根系对不同退化程度黑土土壤微生物群落恢复演替的影响机制”(32471697)
黑龙江省博士后基金(LBH-Z25057┫与国家重点研发计划┣2021YFD150070506)