薄壳山核桃青皮分解及养分释放规律

侯 婷，袁 军，周文君，李 俊，周乃富

（1.中南林业科技大学 林学院，湖南 长沙 410004；2.湖南省水资源研究和利用合作中心，湖南 长沙 410031）

【研究意义】植物残体是植物在生长发育过程中新陈代谢的产物，其分解是生态系统养分循环的重要环节，也是维持土壤肥力的基础。在分解过程中逐渐把养分归还给土壤，对增强土壤结构稳定性，减少水土流失，改善微生物活性等方面发挥重要作用[1-3]。由于不同植物残体的化学组成不同，影响分解过程的因素又复杂多样，与其本身的理化性质以及外界条件等因素密切相关，因此不同分解物分解过程的研究结果存在差异[4-7]。【前人研究进展】薄壳山核桃（Carya illinoensis）又名长山核桃、美国山核桃、碧根果，原产美国和墨西哥，于1900 年左右引入我国，是国内外广泛种植的一种优良的果材兼用树种[8]。通过上百年的引种栽培及示范推广，目前薄壳山核桃在我国生产发展较为迅速[9]。近年来，薄壳山核桃种植面积和区域不断增加，据统计，目前薄壳山核桃种植面积约20 万hm2，盛产期薄壳山核桃青皮产量将达18 万t，同时栽培繁育手段的发展也促进了薄壳山核桃的品质和产量，随之而来的是薄壳山核桃外果皮废弃物的问题。薄壳山核桃外果皮现如今虽利用程度不高，但通过前人的研究表明，它可以被用于多个方面。Novak 等[10]研究了薄壳山核桃的外果皮在700 ℃下热解，可以形成生物炭，能有效提高土壤中的碳含量，并且对土壤的聚合、渗透和持水能力有一定的改善作用。Vaghett等[11-12]采用动力学模型研究了薄壳山核桃外壳作为生物吸收剂吸收溶液中Cr3+、Fe3+、Zn2+、Cu2+、Mn2+、Pb2+等的作用，发现其对这些有毒重属物质均具有较好的吸收作用。梁琼[13]研究表明，从薄壳山核桃外果皮中提取出的酚类化合物对植物炭疽病菌、腐烂病菌等具有良好的抑制作用。【本研究切入点】薄壳山核桃外果皮又叫做青龙衣、青皮[14]，通常情况下薄壳山核桃废弃青皮通过焚烧或堆弃的方式处理，处理方法极不科学，造成资源浪费。不仅造成环境的巨大压力，对人体的健康问题也产生一定程度的威胁[15]。目前有大量文献研究山核桃的生长发育规律，而对薄壳山核桃青皮分解的研究却鲜见报道。【拟解决的关键问题】本研究以粉碎后的薄壳山核桃青皮为试验材料，采用尼龙网袋原位分解法研究青皮的干物质分解和养分释放规律，分析不同还田位置对青皮分解的影响。如能充分利用薄壳山核桃青皮分解规律，对薄壳山核桃青皮的再利用不仅能变废为宝，更利于生态环境的保护，以期为缓解环境压力和能源危机提供参考资料。

1 材料与方法

1.1 样地概况

试验地点位于中南林业科技大学春晖山（28°8′N，112°59′E，海拔110 m）。该地为典型的亚热带湿润季风气候，年平均气温17.2 ℃，冬季平均气温4.6 ℃，无霜期年平均279.3 d，平均总降水量1 422.4 mm，冬春季节多雨，夏秋季节高温少雨，日照年平均1 726 h。土壤类型为第四纪酸性红壤，原生植被为亚热带常绿阔叶林。

1.2 试验材料

试验材料为薄壳山核桃‘金华’青皮，于2020 年12 月采集于云南省林业科学院漾濞核桃研究院。采集后首先将薄壳山核桃进行剥壳处理，将剥离的青皮放入高速粉碎机进行粉碎，得到的粉碎物再使用20 目筛网进行过筛，之后放入大小11 cm×17 cm，孔径0.178 mm 大小的尼龙网分解袋中，每袋准确称取30.0 g。表1 为薄壳山核桃青皮初始养分含量。

表1 薄壳山核桃青皮的初始化学组成Tab.1 Initial chemical composition of Carya illinoensis husks

1.3 试验设计

在春晖山选取一块较为平整的试验地，向下挖掘30 cm，将24个内径20 cm，高30 cm的pvc管随机放入试验地，确保在试验地中分布均匀，且与地表平行，每个管子之间间隔20 cm，回填原土。将之前粉碎过筛处理的薄壳山核桃青皮分为24个尼龙网袋，样地划分为2个处理组，第1个处理组将12个网袋直接放置在pvc 管土柱的上层，即地表处理，用网纱覆盖，橡皮筋固定。第2 个处理组将另外12 个网袋放入pvc管子中部10 cm，即覆土处理。

2020 年12 月21 日，将样品置于样地。试验地中尼龙网袋的收集始于2021 年3 月21 日，为第1 次收集，此后每90 d收集1次。试验开始后分别于第3、6、9、12个月共4次进行对角线收集，取样时抖落尼龙网袋上附着的土壤和杂物，每次每处理各取3个重复，共计6袋，带回实验室。

1.4 指标测定

将取回的分解网袋装于纸质信封袋中置于60 ℃烘箱中烘干至恒重，用精度到万分之一的分析天平称量网袋中剩余青皮的质量，计算青皮干物质残留率。青皮总有机碳（TOC）含量采用重铬酸钾氧化-分光光度法测定。将青皮用浓硫酸-双氧水消解至澄清后用于其他元素含量的测定，氮（N）、磷（P）含量采用全自动间断式化学分析仪测定，钾（K）含量采用火焰光度计测定。

采用Olson负指数衰减模型[16]对青皮分解过程进行拟合，公式为：

式（1）中，y为青皮干物质月残留率（%）；k为分解系数；a为拟合参数；t为分解时间（月）。

青皮干物质残留率计算公式[17]为：

式（2）中，Re为青皮分解干质量残留率（%），Mt为t时刻的青皮剩余干质量（g），M0为青皮初始干质量（g）。

青皮元素残留率计算公式[18]为：

式（3）中，Re为元素残留率（%），Wt为t时刻养分元素含量（%），Mt为t时刻的青皮剩余干质量（g）；W0为初始养分元素含量（%），M0为青皮初始干质量（g）。

采用SPSS 26.0、MATLAB 2020b和Excel 2016软件进行数据统计分析，采用Origin 9.8软件绘图。

2 结果与分析

2.1 不同处理中薄壳山核桃青皮干物质残留率的动态变化

2个处理中薄壳山核桃青皮分解速率不同，地表处理的青皮在分解到第3 到第6 个月时分解速率最快，覆土处理的青皮在分解到第6到第9个月时分解最快，分别为6.49%和12.14%。地表处理的青皮只有在分解第3 到第6 个月时分解速率大于覆土处理部分，其他时间覆土处理的分解速率均大于地上部分（图1）。在分解到第6个月后，两者的干物质残留率产生显著差异（P＜0.05），在分解到第9 个月后，两者产生极显著差异（P＜0.001）。分解到第12个月时，地表处理和覆土处理的薄壳山核桃青皮干物质残留率分别为87.09%和66.25%。

图1 薄壳山核桃青皮干物质残留率变化Fig.1 Changes of remaining rates of dry matter in Carya illinoensis husks

利用Olson 负指数衰减模型对青皮干物质残留率进行动态拟合，R2表明拟合方程可以较好地表示薄壳山核桃青皮的分解速率（表2）。覆土处理的青皮k值大于地表处理的青皮，分解系数k越大，说明分解速率越快，分解50%和95%所需时间越短。根据方程预测，薄壳山核桃青皮两种处理下分解50%分别需要53.05 个月和19.74 个月，分解95%分别需要228.3 个月和85.72 个月。

表2 薄壳山核桃青皮的Olson分解模型Tab.2 Olson decomposition models of Carya illinoensis husks

2.2 不同处理中薄壳山核桃青皮各元素含量及残留率的动态变化

2.2.1 TOC 含量及其残留率的动态变化青皮TOC 含量随着分解时间的增加呈先上升后下降趋势（图2）。在分解的前3 个月，TOC 含量迅速上升，分解到第3～6 个月期间迅速下降，6 个月后含量逐渐降低，两组处理的青皮TOC 含量在整个分解期无显著差异。青皮TOC 残留率在分解的12个月期间先增加后减少（图3）。其中，分解到第3 到第6 个月降低幅度最大。两组处理的青皮TOC 残留率在整个分解期无显著差异。

图2 青皮TOC含量的动态变化Fig.2 Dynamic changes of TOC content in Carya illinoensis husks

图3 青皮TOC残留率的动态变化Fig.3 Dynamic changes of TOC remaining rate in Carya illinoensis husks

2.2.2 N 含量及其残留率的动态变化随着分解时间的增加，地表处理的青皮N 含量呈“上升-下降-上升”的趋势，覆土处理的青皮N含量呈逐渐上升的趋势（图4）。地表处理的青皮在分解到第3个月时N含量高于覆土处理部分，其余时间覆土处理的青皮N 含量均高于地表处理。当分解到第9 个月时，两组处理N 含量存在显著差异（P＜0.05）。地表处理的青皮N 残留率呈“富集-释放-富集”的趋势，覆土处理的青皮N 残留率表现为富集（图5）。地表处理的青皮N 残留率波动幅度高于覆土处理。研究结果表明埋深明显加快了N元素的积累。

图4 青皮N含量的动态变化Fig.4 Dynamic changes of N content in Carya illinoensis husks

图5 青皮N残留率的动态变化Fig.5 Dynamic changes of N remaining rate in Carya illinoensis husks

2.2.3 P 含量及其残留率的动态变化随着分解时间的增加，青皮P 含量呈“降低-增加-降低-增加”的波动趋势（图6）。在分解到第3 个月时，两组处理青皮P 含量无显著差异，第6 到第9 个月覆土处理的青皮含量高于地表处理，分解到第9～12 个月时两者无显著差异。在分解到第3 个月和第12 个月时，地表处理青皮P残留率略高于覆土处理（图7），其余时间覆土处理的青皮P残留率均高于地表处理，2个处理青皮P残留率无显著差异（P＜0.05）

图6 青皮P含量的动态变化Fig.6 Dynamic changes of P content in Carya illinoensis husks

图7 青皮P残留率的动态变化Fig.7 Dynamic changes of P remaining rate in Carya illinoensis husks

2.2.4 K含量及其残留率的动态变化在分解的前3个月，2个处理中青皮K含量急剧下降，分解到第3～6 个月期间后K 含量缓慢增加，分解到第6～12 个月K 含量缓慢降低（图8），地表处理和覆土处理两者无显著差异。青皮K残留率的变化趋势和K含量相似（图9），在分解的前3个月急剧下降，地表处理和覆土处理分别释放了总量的98.8%和97.91%，随后缓慢增加后又缓慢降低，2 个处理中青皮K 残留率无显著差异，这与国内外研究中关于枯落物中K元素动态的研究结果一致。

图8 青皮K含量的动态变化Fig.8 Dynamic changes of K content in Carya illinoensis husks

图9 青皮K残留率的动态变化Fig.9 Dynamic changes of K remaining rate in Carya illinoensis husks

3 结论与讨论

3.1 讨论

植物残体分解包括粉碎、淋溶和分解代谢，是一个复杂的物理、化学和生物学过程，受气候条件、分解物类型、土壤微生物和动物等因素的影响[19]。薄壳山核桃青皮分解先慢后快，是外界环境条件和青皮成分变化共同作用的结果。分解前期青皮的质量损失主要是通过非结构性化合物的淋溶作用实现，试验开始时（12 月），薄壳山核桃青皮中易分解物质的含量高，但由于此时气温低，降水少，淋溶作用弱，微生物活性降低，不利于青皮分解，因此分解速度缓慢；从3 月开始，气温回升，降水量增加，雨热同期，良好的水热条件使微生物繁殖迅速，青皮分解速率明显增加。本研究表明，在薄壳山核桃青皮的分解过程中，覆土处理的青皮在整个分解期分解速率比地表处理的青皮快，主要是因为覆土处理的青皮N 含量高于放于地表处理，N 含量越高，微生物可以利用的氮源就越高，且地下水热条件好，更有利于增加微生物的生物量及活性[20]，且覆土处理青皮与土壤的接触面积大，因此有利于微生物参与分解，分解迅速。还有研究表明[21]，在5～10 cm 深度的土壤脲酶、多酚氧化酶和过氧化氢酶含量较高，有利于植物残体分解。

薄壳山核桃青皮分解伴随着养分元素的变化，这些元素的释放或富集在循环过程中起着重要作用。本研究表明，青皮分解过程中的TOC、N、P、K 变化特征不同，这是由于青皮分解及元素自身特性所决定的。薄壳山核桃青皮TOC 的主要成分是木质素和纤维素，TOC 含量随着分解时间的增加逐渐降低。在分解前期，TOC 的损失量小于干物质损失量，所以TOC 浓度表现为上升。分解到第3 个月时，随着木质素和纤维素等有机物的大量分解，TOC浓度开始下降[22]。由于薄壳山核桃青皮在分解过程中主要以C为主导的物质循环模式，易分解碳水化合物的分解以及可溶性物质的淋失都会引起青皮TOC 的释放。青皮在分解中N 的释放与富集主要取决于分解初期N 含量能否维持微生物活动的需要[23]。覆土处理的青皮在分解期间N 含量一直处于增加状态，地表处理的青皮N 含量先增加后降低再增加，这可能是由于土壤的淋溶作用导致有效态的养分向下淋失，有机质在分解过程中带有负电荷，可以吸附阳离子，覆土处理的青皮作为有机质在分解过程中可能吸附了上层土壤淋溶下的养分，2个处理的N 含量在分解第9个月时存在显著差异。覆土处理青皮N 含量显著高于地表处理，这是由于地下部分的水热条件好，微生物大量繁殖，青皮N 含量不能满足分解的需要，微生物得不到充足的养分来维持自身生长，因此分解者会在周围环境中固持N[24]，即青皮N含量增加以促进分解。P元素属于易淋溶元素，主要以硫酸根或化合物的形态存在[25]，在青皮分解初期，温度较低，微生物量少，P 可以供微生物的营养需要，所以P 会表现为释放。到分解第3个月到第8个月，温度上升，降雨增多，湿度较大，微生物活性增大，微生物对P的需求不足时，会出现累积[26]。分解到第9～12 个月，温度逐渐下降，微生物数量减少，青皮腐解P 元素释放较快。温度是P 元素的主要限制因子。K 元素在植物体内主要以离子状态存在，在分解过程中流动性很强，其浓度在分解开始时便急剧下降，后趋于平缓[27-28]。国内外研究中关于枯落物中K元素动态的研究结果基本为直接释放[29]，与本研究结果一致。综上所述，覆土加快了薄壳山核桃青皮的分解和N的累积，建议在薄壳山核桃青皮还田时采用深施覆土方式。

分解是一个漫长的过程，本研究中仅探究了薄壳山核桃青皮的分解规律及埋深对青皮分解的影响，而青皮分解过程对土壤微生物、养分的影响，氮沉降、磷沉降和酸沉降对青皮分解的影响以及薄壳山核桃青皮还田对大田作物的影响等尚需进一步研究。

3.2 结论

薄壳山核桃市场前景广阔，综合效益高，是国内外广泛种植的一种优良的果材兼用树种。为给薄壳山核桃青皮的资源化利用提供参考，本研究采用分解网袋法进行原位分解试验，分析不同还田位置处理下薄壳山核桃青皮的分解及养分释放规律，主要结论如下：

（1）不同还田位置处理下的薄壳山核桃青皮分解速率不同。覆土处理的青皮分解速率显著快于地表处理，分解速率随分解时间的增加呈增大趋势，分解时间越久，两组处理的分解速率差异就越显著。

（2）薄壳山核桃青皮的干物质残留率符合Olson 负指数衰减模型，前期缓慢，后期分解加快，地表处理和覆土处理的青皮半衰期分别为53.05 个月和19.74 个月，分解95%所需时间分别为228.3 个月和85.72个月。

（3）两组处理的青皮TOC 残留率均呈富集-释放趋势，N 残留率地表处理的青皮呈“富集-释放-富集”的趋势，覆土处理的青皮表现为富集趋势，P 残留率均呈释放-富集-释放趋势，K 残留率均呈释放趋势，薄壳山核桃青皮养分元素残留率与环境条件和微生物活性密切相关。