分解速率与初始质量的关系
2019-03-26王多锋李国林戚建莉柴春山
魏 强, 凌 雷, 王多锋, 李国林, 戚建莉, 柴春山
(1.甘肃省林业科学研究院, 甘肃 兰州 730020; 2.甘肃兴隆山国家级自然保护区管理局, 甘肃 兰州 730117)
森林凋落物是森林植物生长发育过程中新陈代谢的产物,作为连接植物与土壤的“纽带”和养分的基本载体[1],是森林生态系统自肥的主要补给者,在能量流动和物质循环、维持土壤肥力、改善土壤理化性质、促进森林生态系统生产力、增强森林水源涵养与水土保持功能等方面起着重要作用;同时,凋落物还是森林土壤动物和微生物的能量和物质来源,对促进和维持整个森林生态系统平衡起着重要作用[2]。因此,森林凋落物是森林生态学、森林土壤学、森林水文学、生物地球化学及环境化学等学科的重要研究内容之一[3]。
凋落物分解是将生物大分子降解为无机物小分子,大部分转化为CO2和水的过程,最后只剩下少量的腐殖质进入土壤层,为地上植物生长与发育提供所需养分[4]。凋落物分解一般包括物理、化学和生物学过程,是在降水淋溶、动物咀嚼和啃食、土壤干湿交替和冻融、生物代谢等相互作用下共同完成的[5]。凋落物分解是生态系统物质循环和能量流动的重要组成部分,在维持土壤肥力、促进系统物质循环、土壤有机质形成以及群落演替等方面具有不可替代的作用和地位[6]。
凋落叶分解速率(k)是衡量凋落叶分解快慢和养分循环效率的关键指标,目前主要采用凋落物网袋(litter bag)法进行研究。凋落叶质量可表征凋落叶的相对可分解性[5],并依赖于其易分解成分(N,P,水溶性糖和氨基酸等)和难分解或较难分解的有机成分(木质素、纤维素、半纤维素和多酚类物质等)的含量及相对比例[4]。在区域气候尺度上,凋落叶质量是控制凋落叶分解快慢的关键因子[7]。探讨凋落叶初始质量与分解速率的关系,对于了解不同生态系统的物质循环速率具有重要意义[4]。
采用初始质量指标预测分解速率是凋落叶分解研究的重要内容。研究表明,凋落叶分解速率与凋落叶初始质量关系密切[8-9];碳含量及碳/氮是影响分解速率的关键因子,两者均与分解速率呈负相关[10];Dyer等[11]和陈法霖等[12]认为木质素含量是预测凋落叶分解速率的关键指标,且与分解速率呈负相关;也有人认为[13]分解速率与初始碳含量及碳/氮的关系不大,而木质素/氮是影响分解速率的主要因子;彭少麟等[14]认为碳/氮比和木质素浓度是制约凋落物分解速率最为重要的因素。此外,初始磷、半纤维素、酚类、钙、钾、锰和硅等含量对分解速率的影响也在许多研究中被证实[15]。可以看出,不同研究区域得到的初始质量指标与凋落叶分解速率的关系差别较大,研究还存在很多不确定性。
兴隆山森林位于黄土高原最西端与祁连山东延余脉的交汇处,是以黄土和石质山地为主要成土母质、在降水量不足条件下发育而成的隐域森林,为半干旱区保留较完整的一块森林,其中有大面积原始青杄(Piceawilsonii)林分布,而青杄纯林是该区域森林自然演替过程中的亚顶极群落。兴隆山森林是兰州市天然的生态屏障,对于兰州市及周边地区在保持水土、涵养水源、防风固沙、调节气候、改善空气质量等方面具有不可替代的作用。因此,对该区域森林生态系统的稳定性、生态功能、物质循环与能量流动等方面的研究显得尤为重要。目前,对兴隆山森林凋落物累积量、凋落物凋落量等进行了相关研究[16-17],但对凋落叶分解的研究还未见报道,特别是主要树种凋落叶分解速率大小及分解速率与凋落叶初始质量的关系如何还不清楚。因此,本文拟采用凋落物分解袋法,以典型山杨(Populusdavidiana)—白桦(Betulaplatyphylla)—青杄针阔混交林、青杄中龄林和青杄近熟林3种森林的青杄、山杨和白桦的凋落叶为研究对象,开展野外和室内试验,测定凋落叶分解速率和叶初始质量,阐明凋落叶各初始质量指标与分解速率的关系,以期为兴隆山森林生态系统物质循环的深入研究提供依据。
1 研究地区与研究方法
1.1 研究区概况
甘肃兴隆山国家级自然保护区位于兰州市东南约45 km处(103°50′—104°10′E,35°38′—35°58′N),属于祁连山的东延余脉,东西长37 km,南北宽17 km,总面积33 301 hm2,海拔1 800~3 670 m,是黄土高原最西端的一座绿色岩岛。属于温带半湿润半干旱气候类型,在海拔<1 800 m区域,年均降水量350 mm;海拔在1 800~2 600 m区域,年降水量为400~600 mm;海拔在2 600~3 000 m区域,年降水量为600~800 mm。地貌以石质山地和山间谷地为主要特征,土壤由高山草甸土、亚高山草甸土、灰褐土、栗钙土、黄绵土、新积土组成。保护区内天然林主要有山杨—白桦林、山杨—辽东栎(Quercuswutaishanica)林、山杨—白桦—青杄林、青杄纯林等[18]。
山杨—白桦—青杄林分布在阴坡,建群种为青杄、山杨、白桦和红桦(B.albosinensis),灌木层树种有陕甘花楸(Sorbuskoehneana)、红毛五加(Acanthopanaxgiraldii)、细枝绣线菊(Spiraeamyrtilloides)、藤山柳(Clematoclethralasioclada)、华西箭竹(Fargesianitida)等;青杄重要值达53.0%,成为群落的主要建群种,已进入主林冠层;白桦、红桦等重要值较小[18]。青杄中龄林分布在阴坡,建群种为青杄,乔木层中还混有少量山杨、白桦和皂柳(Salixwallichiana)等,灌木层树种有华西箭竹、红毛五加、陕甘花楸、细枝绣线菊和蒙古荚蒾(Viburnummongolicum)等;青杄重要值达78.2%,在林冠层中占绝对优势;山杨、白桦等重要值很小,趋于消失阶段[18]。青杄近熟林分布在阴坡,建群种仅为青杄,灌木层树种有华西箭竹、红毛五加、陕甘花楸等;青杄重要值达100.0%,初步形成青杄纯林;灌木层以华西箭竹为主要树种[18]。林分主要特征见表1。
表1 甘肃兴隆山自然保护区林分概况
注:林分组成中各树种前数字表示株数所占成数。
1.2 研究方法
1.2.1 样地设计 在山杨—白桦—青杄针阔混交林、青杄中龄林和青杄近熟林3种森林中,选择固定样地,总共样地11个。山杨—白桦—青杄林4块样地,分别位于海拔2 400,2 450,2 500,2 500 m处;青杄中龄林4块样地,分别位于海拔2 400,2 450,2 500,2 550 m;青杄近熟林3块样地,分别位于海拔2 300,2 400,2 500 m。样地均分布在阴坡,大小为30 m×30 m的正方形。2012年6月,在每个样地各安放5个长宽高1.00 m×1.00 m×0.20 m的凋落物收集器,用于凋落物年凋落量的动态研究和凋落叶分解试验研究。
1.2.2 试验设计 2013年9至10月,在兴隆山森林凋落物凋落的高峰期,收集各样地凋落物收集器内的凋落物,带回试验室,挑选分离出青杄针叶、山杨叶和白桦叶可做为凋落叶分解试验的材料。其中,白桦叶在山杨—白桦—青杄林海拔2 400,2 450 m样地内收集,山杨叶在山杨—白桦—青杄林海拔2 500,2 550样地内收集。将各样地所挑选分离的青杄中龄林针叶、青杄近熟林针叶、山杨叶和白桦叶分别进行均匀混合,等样品自然风干后,每种凋落叶分别称取5份,每份10 g,于65℃烘箱烘干至恒重,测定各凋落叶中的含水量,由此推算各凋落叶样品的含水率及初始干重。
采用尼龙网制成大小20 cm×20 cm、网孔大小为0.15 mm(100目)的凋落袋,进行凋落叶分解试验研究。称取自然风干凋落叶,每种凋落叶每袋称取相当于烘干重的10 g左右,青杄中龄林针叶108袋,青杄近熟林针叶36袋,山杨叶36袋,白桦叶36袋,总共216袋。2013年11月5日,将青杄中龄林针叶、青杄近熟林针叶、山杨叶和白桦叶放置在位于海拔2 450 m的青杄中龄林中,用于研究不同森林主要树种凋落叶在同一立地条件下的分解速率,每种凋落叶各放36袋;将其余青杄中龄林针叶放置在位于海拔2 550 m和2 650 m的青杄中龄林中,用于研究同一森林类型青杄针叶在不同海拔梯度下的分解速率,每种凋落叶各放36袋。凋落物分解袋平铺于林地土壤层表面,使其紧贴土壤层,各凋落袋间隔>5 cm。分别于2014年、2015年、2016年的5月5日、8月5日和11月5日采回凋落叶分解袋,每次每种凋落叶每个试验点各采集4袋,采集时去除凋落物袋表面的土块、碎石和琐屑等杂物;带回实验室,用水冲洗,去除表面杂质和土壤以及袋内由降水冲刷而进入的土壤,65 ℃下烘至恒量,测定凋落叶残留率,用3 a的凋落叶分解观测数据来计算各树种凋落叶的分解速率。
1.2.3 测定项目与方法 将烘干的凋落叶样品研磨,分析其全碳、全氮、全磷、全钾、木质素和综纤维素含量。全碳含量采用红外吸收光谱测定,仪器为CS2008碳硫分析仪(北京纳克公司);全氮含量采用半微量凯氏定氮法测定,仪器为K9860凯氏定氮仪;全磷含量采用硫酸—高氯酸消煮—钼锑抗分光光度法,仪器为紫外可见分光光度计;全钾含量采用硫酸—高氯酸消煮—原子吸收光谱法测定,仪器为M6原子吸收光谱仪(美国热电公司);木质素含量采用索氏提取器法测定;综纤维素含量采用索氏提取器法测定。
1.3 数据处理与分析
凋落叶残留率计算公式为:
y=Wt/W0×100%[4,19]
式中:y——凋落叶经时间t后的残留率(%);W0——凋落叶初始干质量(g);Wt——凋落叶经时间t后的残留干物质(g)。
采用修正的Olson指数衰减模型[20]模拟凋落叶分解过程,计算凋落叶分解速率(k)和半分解期(t0.5)及分解95%所需时间(t0.95):y=ae-kt
式中:a——修正系数;t——分解时间(a)。凋落叶50%的分解期(50%分解)计算式为t0.5=ln(0.5/a)/(-k);凋落叶95%的分解期(95%分解)计算式为t0.95=ln(0.05/a)/(-k)。
数据统计与分析采用SAS 9.0和Excel 2007软件进行。
2 结果与分析
2.1 凋落叶初始质量
甘肃省兴隆山3种乔木树种凋落叶初始质量差异较大(表2),碳、氮、磷、钾含量分别为362.000~504.000 g/kg,9.600~10.500 g/kg,0.740~1.250 g/kg,2.050~3.800 g/kg,纤维素含量为516.950~577.225 g/kg,木质素含量为388.430~420.715 g/kg;碳/氮为38.850~52.506,木质素/氮为37.077~43.829,碳/磷为403.846~673.096。对于碳、钾、碳/氮含量来说,阔叶明显高于针叶;对于氮、磷、纤维素、木质素、木质素/氮、碳/磷含量,针叶和阔叶差别不明显,无明显变化规律。对于青杄树种来说,不同龄林间除了氮和磷之外,其他成分含量差别较大。
表2 兴隆山森林3种主要树种凋落叶初始值
2.2 凋落叶分解速率
由表3可以看出,在相同条件下(即凋落叶分解所处立地条件和林分因子、凋落物分解袋网孔大小均一致),兴隆山森林3种主要树种的凋落叶分解速率(k),50%分解时间(t0.5)和95%分解时间(t0.95)差异比较明显。青杄中龄林针叶、青杄近熟林针叶、山杨叶和白桦叶分解系数k值变化范围为0.11~0.16,其中最低的是山杨叶(0.11)和白桦叶(0.11),其次是青杄近熟林针叶(0.13),k值最高的是青杄中龄林针叶(0.16)。研究表明,在凋落叶分解相同条件下,阔叶分解速率均较慢,青杄针叶分解快于阔叶,特别是青杄中龄林针叶明显要快于青杄近熟林针叶。说明同一树种林龄不同,凋落叶分解速率也不相同,龄林较小的针叶分解速率较快。3种树种凋落叶t0.5为4.45~6.71 a,t0.95为19.08~28.57 a。在甘肃兴隆山林区,在分解袋网孔大小为0.15 mm的条件下,青杄针叶分解速率明显要快于山杨和白桦,凋落叶全部分解至少需要19~24 a以上的时间;山杨和白桦叶全部分解所需时间更长,可达29 a以上,这一结论与目前研究结果不相吻合。按目前研究结果来看,阔叶分解速率明显要快于针叶,产生这一结果的原因很可能是凋落叶分解主场效应和分解袋孔径较小所致。对于青杄中龄林针叶来说,海拔2 450 m处k值为0.16,而2 550,2 650 m处k值均为0.15;t0.5为4.45 ~4.65 a,t0.95为19.08 ~19.98 a。3种海拔梯度间青杄针叶分解速率差别不明显,说明在100 m 海拔梯度下,海拔因子对青杄中龄林针叶分解速率影响不明显。
表3 兴隆山森林3种主要树种凋落叶分解的指数模型及参数
总体看来,在分解网孔大小为0.15 mm分解袋所处林分和立地因子(即3种树种凋落叶均放置在青杄中龄林中进行试验)均为一致的条件下,针叶分解速率明显要快于阔叶,如青杄中龄林针叶较快,青杄近熟林针叶此之,而山杨—白桦—青杄林的山杨和白桦叶最慢。
2.3 凋落叶初始质量与分解速率的关系
2.3.1 凋落叶初始质量与分解速率的回归分析 由图1可以看出,凋落叶年分解速率(k)与凋落叶初始氮含量呈显著线性正相关,而与木质素含量、碳/氮值、木质素/氮值和钾含量呈显著线性负相关,特别是与木质素含量、氮含量和木质素/氮值,相关系数较高,均达0.700 0以上,说明木质素含量、氮含量和木质素/氮值对兴隆山3种树种凋落叶分解速率影响较大,在一定程度上决定着凋落叶的分解速率。凋落叶年分解速率(k)与凋落叶初始纤维素含量、磷含量和碳/磷值线性回归方程不显著,说明纤维素含量、磷含量和碳/磷对兴隆山3种主要树种凋落叶分解速率影响不明显。
由此看出,凋落叶初始质量在一定程度上决定着凋落叶的分解速率,凋落叶初始氮含量越高,分解速率越高;木质素含量越高,分解速率越低;木质素/氮比值越高,凋落叶分解速率越快;另外,钾含量也是影响凋叶分解速率的一个因素。
图1 兴隆山森林3种主要树种凋落叶分解速率与凋落叶初始质量的线性回归分析
2.3.2 凋落叶初始质量各组分对分解速率的相对贡献 对兴隆山3种主要树种凋落叶分解速率与凋落叶初始质量影响因子(木素含量、全氮、全磷、全钾、全碳、碳/氮、碳/磷和木质素/氮)进行主成分分析,结果见表4。由表4看出,第1主成分的特征根为5.481 5,方差贡献率为60.91%,在主成分函数表达式中全钾含量、木质素含量、木质素/氮值对主成分有较大的权重贡献,特别是全钾贡献最大。第2主成分的特征根为2.551 9,方差累积贡献率达89.26%,在主成分函数表达式中碳/磷值、纤维素含量和全磷含量对主成分有较大的权重贡献,特别是碳/磷贡献最大。第3主成分的特征根为0.767 3,方差累积贡献率达97.79%,在主成分函数表达式中全碳含量、氮含量、碳/氮值对主成分有较大的权重贡献,特别是全磷含量贡献最大。由主成分分析可以看出,全钾含量、木质素含量、木质素/氮、碳/磷、纤维素含量是影响兴隆山青杄、山杨和白桦凋落叶分解速率的重要指标。
表4 兴隆山森林3种主要树种凋落叶初始质量主成分分析
3 讨 论
3.1 兴隆山森林主要树种凋落叶的分解速率
Olson模型拟合凋落叶失重过程得到的分解速率(k)是一个恒值,可以揭示凋落叶分解过程的基本特征[4]。本研究发现兴隆山森林3种主要树种凋落叶的分解速率,其中青杄k值较低,仅为0.13~0.16,这一结果与川西亚高山岷江冷杉k相当[4];山杨、白桦k值更低,仅为0.11,远低于其他研究区的研究结果,如长白山山杨k值为0.34~0.48,白桦k值为0.41~0.44[13,21]。对于山杨、白桦凋落叶来说,分解速率很低,其原因之一可能是凋落物分解袋网孔太小所致。本研究中所采用尼龙网孔大小为0.15 mm,并且表面和贴地面大小一致,均为0.15 mm,而李雪峰等研究所用网孔大小为贴地面0.5 mm、表面为1.0 mm[13,21]。其二可能是凋落物所处林分环境不同所致,即凋落物分解存在主场效应[22],凋落物具有在其生长的栖息地(主场)比在别的生境(客场)分解更快的特征。本研究是将山杨—白桦—青杄林的山杨和白桦叶、青杄近熟林的青杄针叶、青杄中龄林的青杄针叶均放入青杄中龄林内,进行山杨、白桦、青杄中龄林针叶和青杄近熟林针叶4种凋落叶在相同条件下分解速率的研究,青杄中龄林和山杨—白桦—青杄林由于林下光照条件及地表凋落物的不同,造成了2种森林类型土壤化学特性和土壤微生物的不同,进而影响到地表凋落物分解速率的不同,由此造成了山杨、白桦叶分解速率较慢这一结果。李雪峰等[13,21]是将山杨和白桦叶放在山杨—白桦林内进行其分解速率的研究,他们的研究结果反映了山杨、白桦叶在其自身森林中的分解速率。另外,青杄中龄林的针叶分解速率高于青杄近熟林的针叶,这有可能也是凋落物分解存在主场效应原因所致。因此,凋落物分解袋孔径大小、凋落物分解所处林分环境状况(如土壤微生物和光照等)将很可能影响到凋落叶的分解速率,具体如何影响还有待于进一步深入研究。
3.2 凋落叶初始质量与分解速率的关系
兴隆山森林主要树种凋落叶初始质量在一定程度上决定了凋落叶的分解速率,凋落叶分解速率随初始氮含量增高而增大,随凋落叶初始木质素含量、木质素/氮、全钾含量增高而减小,其余指标对凋落叶分解速率无显著作用。凋落叶分解速率随氮含量增高而增大,这与唐仕珊等[9]、Zhang等[8]和杨林等[4]对国内外森林凋落叶分解研究结果相一致。氮的有效量常被作为确定分解速率的控制因子,因为氮决定微生物生物量的增长和微生物矿化的有机碳周转[4]。在本研究中,凋落叶分解速率随初始木质素含量增高而减小,符合凋落叶分解的一般理论[11],因木质素结构复杂,它是凋落叶中难分解的主要成分。
目前普遍认为碳/氮值和木质素/氮值最能反映凋落叶的分解速率[23],也有研究发现初始碳含量及碳/氮是影响分解速率的关键因子[10]。还有研究发现,碳/氮与分解速率关系不大,木质素/氮才是影响凋落叶分解速率的主要因子[13]。可见,不同研究区域相同树种或同一研究区域不同树种凋落叶各初始质量指标与分解速率的关系不尽相同。在本研究中,木质素含量、木质素/氮、碳/磷、纤维素含量、钾含量是影响兴隆山森林主要树种凋落叶分解速率的重要指标,而木质素含量、木质素/氮是控制兴隆山青杄、山杨和白桦凋落叶分解速率的关键因子,凋落叶初始木质素/氮越高,分解速率越低。这一结果与Constantinides等[24]、Scott等[25]和李雪峰等[13]的研究相一致。本研究也发现,除了木质素含量和木质素/氮外,钾含量也是影响凋落叶分解速率的一个主要因子,这与唐仕珊等[9]研究结果相似。
3.3 凋落叶分解速率与分解袋孔径大小的关系
在凋落物分解速率研究中,除了凋落叶自身质量因素外,凋落物分解袋网孔大小对凋落叶分解速率有着重要影响。本研究发现,山杨、白桦凋落叶分解速率较慢,k值仅为0.11,并且低于青杄针叶。在正常情况下,阔叶分解速率远大于针叶,因阔叶中难分解物质较少,但本研究结果却低于其他相同类型区山杨、白桦凋落叶的分解速率[13,21]。笔者认为,导致山杨、白桦叶分解速率低的一个主要原因是凋落叶分解袋网孔较小所致。本研究所用尼龙网孔径为0.15 mm,而其他地区为0.5 mm或1.0 mm[13,21]。李艳红等[26]采用不同孔径凋落袋(6目、30目和260目)对四川省乐山市苏稽镇不同比例巨桉(Eucalyptusgrandis)与台湾桤木(Alnusformosana)混合凋落物分解过程进行了研究,结果表明孔径6目的分解速率最大,30目的次之,260目的最小。当网孔大小为6目时可允许各种土壤动物的作用;网孔30目时则基本排除了大型土壤动物的影响,允许中小型土壤动物自由出入;网孔260目时则尽可能排除土壤动物的作用[26]。从以上分析看出,由于凋落物袋网孔大小的不同,将直接影响土壤动物能否参与分解袋中凋落物的分解,进而影响到凋落物分解速率的大小。因此,在凋落物分解研究过程中,凋落物分解袋网孔大小对研究结果影响较大,不同树种凋落物分解袋网孔大小应不尽相同。
4 结 论
研究表明,甘肃兴隆山青杄凋落叶分解较慢,分解速率为0.15~0.16,半分解期为4.45~4.65 a,95%分解期为19.08~19.98 a,一个完整的分解周期至少要20 a以上的时间,该地区海拔较高,全年温度偏低,分解者活性不高可能是其重要原因。线性回归和主成分分析表明,木质素含量、氮含量、木质素/氮、全钾含量是控制兴隆山森林乔木树种凋落叶分解速率的关键因子,凋落叶初始木质素/氮越高,分解速率越低。除了凋落叶基质质量之外,凋落物分解主场效应和分解袋孔径大小也是影响凋落叶分解速率快慢的重要因子,而这两个因子是不能被忽略的。