亚热带地区3种喀斯特林分凋落物组成对其持水性能的影响
2022-02-19伟侯满福刘雨婷贺露炎吴兴澳
唐 伟侯满福刘雨婷贺露炎吴兴澳
(1.广西师范大学 环境与资源学院,广西 桂林541004;2.珍稀濒危动植物生态与环境保护教育部重点实验室/广西漓江流域景观资源保育与可持续利用重点实验室 广西师范大学,广西 桂林541004)
森林具有改善水质、保持水土、涵养水源的功能,是陆地生态系统服务的重要组成部分[1]。凋落物层是森林生态系统实现涵养水源,保持水土功能的核心,它通过截留降水,改变土壤物理性质,防止雨滴击溅侵蚀、减少土壤的蒸发等实现对降水的调蓄[2-3],对全球变化引发的极端降水具有缓冲作用。中国亚热带喀斯特区面积广泛,是珠江、长江上游的生态脆弱区,研究其森林的生态水文功能对下游生态安全具有重要意义。很多研究者通过浸泡法[2-4]和人工模拟降水[5-6]研究凋落物的持水性能,揭示了不同类型[7]、不同林龄[8]、不同林分密度[9]和不同海拔[10]林分凋落物层的持水性能。对于凋落物层的持水性能,目前认为阔叶林普遍高于针叶林、成熟林优于过熟林、中龄林和幼龄林、高林分密度林强于过高和过低林分密度林。然而,凋落物水文生态功能的变化本质上由凋落物的质和量决定,其中凋落量与森林类型密切相关[11],凋落物的质则会受物种组成和凋落物组成的影响,且会随时间、群落类型而变化[12]。但绝大多数已有研究仅关注凋落量的影响,只有少量研究对凋落物组成的影响进行了一定的探讨[13-14],对凋落物组成如何影响其水文功能的了解十分缺乏。为此,本研究选择亚热带地区3种喀斯特代表性林分类型,通过对凋落物层的分类并设置不同凋落物类型组成比例的分组试验,以验证凋落物组成对其持水性能的影响规律,为丰富森林生态水文研究和加深凋落物层与水文功能关系的认识提供科学依据。
1 材料与方法
1.1 研究区概况
研究区位于广西壮族自治区桂林市郊,属典型的喀斯特峰丛洼地和峰林平原地貌。该区属亚热带湿润季风气候,多年平均气温17.8~19.1℃,多年平均降雨量1 814~1 941 mm,雨季主要集中在4—7月。该区出露基岩主要为纯质石灰岩,岩溶发育强烈,土壤主要为棕色石灰土。该区植被以灌丛为主,局部区域分布森林,主要是以青冈(Cyclobalanopsis glauca)、岩樟(Cinnamomumsaxatile)、密 花 树(Myrsineseguinii)等常绿树种及光蜡树(Fraxinus griffithii)、化香树(Platycaryastrobilacea)、黄连木(Pistaciachinensis)、光皮梾木(Cornuswilsoniana)、圆叶乌桕(Triadicarotundifolia)等落叶树种组成的混交林。
1.2 研究方法
1.2.1 样地设置样品采集 选择研究区的青冈林、青冈+光蜡树林、化香树+密花树林3种代表性的林分,分别在林内设置20 m×20 m的调查样方(3个重复),开展森林群落调查[15]。各样地概况详见表1。在凋落高峰刚过的初冬期(2019年11月初)于各林分内随机布置7个0.5 m×0.5 m的小样方,将地表凋落物按未分解层和半分解层分别全部收集带回实验室,并进一步将未分解层凋落物分为枝、叶和其他,85℃烘干至恒重,称量干重并计算单位面积凋落物蓄积量。
表1 试验样地基本情况
1.2.2 试验设计 将每种林分类型的半分解、未分解枝和未分解叶3种凋落物成分(其他成分量过少,未纳入分组试验)按比例混合成总重40 g的处理组,每种林分类型21个分组进行浸水试验。表2为浸水试验各林分类型凋落物浸水试验的分组及其凋落物成分组成结果。
表2 各林分类型凋落物浸水试验的分组及其凋落物成分组成结果
1.2.3 凋落物持水率与吸水速率的测定 采用室内浸水法[8]测定凋落物的持水率及吸水速率。将每组凋落物分别装入规格为25 cm×25 cm、孔径为1 mm的尼龙网袋中放入盛有清水的容器中浸泡,分别在浸水1/12,1/4,1/2,1,2,4,6,8,10,12,16和24 h时捞出平置沥干,至凋落物不滴水时迅速称重(精确到0.01 g)。凋落物的持水量、持水率和吸水速率计算公式为:
式中:Wt为吸水湿重(g);W0为干重(g);Q为持水率(g/g);C为持水量(t·hm-2);W为凋落物蓄积量(t·hm-2);V为吸水速率〔g/(g·h)〕;t为浸水时间(h);当t为24:00时,样本持水率及持水量达到饱和[8]。
根据凋落物分组试验中单一组分的最大持水率,结合混合分组的凋落物组成情况,拟合计算其最大持水率,并与实测的最大持水率进行回归分析。具体计算公式为:
式中:Qn为拟合的最大持水率;Q24半为半分解最大持水率;Q24枝为枝最大持水率;Q24叶为叶最大持水率。
混合分组的实际与拟合最大持水率关系为:青冈林:
青冈+光蜡树林:
化香树+密花树林:
式中:X为拟合最大持水率;Y为实际最大持水率。
再根据公式(2)—(3)计算各林分类型凋落物层最大持水量(未分解的其他组分以枝和叶的平均值计)和最大吸水速率(浸水1/12 h)。
1.2.4 数据分析 运用SPSS 19.0进行单因素方差分析,对凋落物的持水率和吸水速率与凋落物各组分的相关关系采用Pearson相关性分析。
2 结果与分析
2.1 不同林分类型凋落物层组成特征
3种林分类型的凋落物层总蓄积量在8.49~9.37 t/hm2间,大小顺序为:青冈林>青冈+光蜡树林>化香树+密花树林,但差异不显著(p>0.05)。凋落物蓄积量中未分解层均大于半分解层,前者占比在53.75%~58.99%之间,其中又以未分解层枝占比较高(26.04%~33.23%),其次为未分解层叶(16.10%~31.61%),其他凋落物所占比例最小(<5%)。不同林分类型的凋落物层总蓄积量、半分解层、未分解层枝及其他凋落物量均无显著差异(p>0.05),仅青冈+光蜡树林未分解层叶显著低于化香树+密花树林(p<0.05),同时也较明显低于青冈林(为后者的约64%),但由于组内变异较大而使其均值的差异未达显著水平(表3)。
表3 不同林分类型凋落物蓄积量及其构成 t/hm2
2.2 不同分组凋落物持水过程
3种林分类型不同分组凋落物的持水率和吸水速率随浸水时间的增加均表现出相似的变化规律。其中持水率与浸水时间呈显著对数函数关系(R2≥0.718),吸水速率对浸水时间的响应可由幂函数描述(R2≥0.998)。凋落物在烘干状态下浸水初期(2 h内)持水率迅速增加(图1—2),2~12 h持水率缓慢增加,12~24 h持水率的变化趋于平缓逐渐饱和;吸水速率在浸水的0~1/12 h最快,1/12~4 h急剧下降,10~24 h维持低位并趋于稳定(图3—4)。
图1 不同林分类型不同分组凋落物持水率与浸水时间的关系
在各浸水时段不同林分类型0∶10∶0(纯枝)的持水率普遍最低(仅青冈林稍高于2∶8∶0)且明显低于其他分组,0∶0∶10(纯叶)的持水率普遍最高但不显著(图2)。此外化香树+密花树林的2∶8∶0和青冈林的0∶8∶2的持水率与纯枝分组相近,其他分组间多彼此交叉,处于连续变化当中。化香树+密花树林的0∶0∶10持水率最大,青冈+光蜡树林的0∶10∶0最小,前者在各浸水时段的持水率比后者多1.07~1.22(g/g)。不同分组凋落物达到最大持水率所需时间存在差异,青冈林的2∶6∶2,青冈+光蜡树林的6∶4∶0,化香树+密花树林的2∶4∶4等需时最短(16 h)。总的来看,从1/12 h到24 h,不同分组凋落物的持水率增长了0.44~0.89(g/g),以青冈+光蜡树林的0∶2∶8的增长值最大,是其他分组的1.04~2.04倍。
图2 3种林分类型单一组分凋落物持水率与浸水时间的关系
由图3可以看出,不同分组凋落物间吸水速率相差较大,但整个浸水时段仍以化香树+密花树林的0∶0∶10的吸水速率最快,青冈+光蜡树林的0∶10∶0最慢。在浸水初期(1 h内),3种林分类型的4∶0∶6吸水速率在各自混合组成分组中最快。不同分组凋落物持水率和吸水速率与浸水时间的关系详见表4(以化香树+密花树林为例)。
图3 不同林分类型不同分组凋落物吸水速率与浸水时间的关系
表4 化香树+密花树的不同分组凋落物持水率和吸水速率与浸水时间的关系
图4 3种林分类型单一组分凋落物吸水速率与浸水时间的关系
2.3 凋落物组成与持水率和吸水速率的关系
凋落物的持水率和吸水速率均随未分解叶凋落物量的增加而增加,随未分解枝量的增加而减小(图1—4)。
Pearson相关分析表明,在全浸水时段3种林分类型不同分组凋落物的持水率和吸水速率均与未分解叶的量呈极显著正相关(n=21,p<0.01),与未分解枝的量呈极显著负相关(n=21,p<0.01),与半分解的量无显著关系(n=21,p>0.05)(表5,以最大持水率和最大吸水速率为例)。
表5 不同林分类型不同分组凋落物最大持水率和最大吸水速率与凋落物各组分的相关系数
2.4 不同林分类型凋落物的最大持水量、最大持水率与最大吸水速率
拟合得到的不同混合分组凋落物的最大持水率总体小于实际最大持水率;当半分解量一定时,随着未分解叶量增加,凋落物最大持水率呈现上升趋势;当未分解枝量一定时,含同等量的未分解叶凋落物的最大持水率大于含同等量的半分解凋落物(图5—7)。
2.4 不同林分类型凋落物的最大持水量、最大持水率与最大吸水速率
拟合得到的不同混合分组凋落物的最大持水率总体小于实际最大持水率;当半分解量一定时,随着未分解叶量增加,凋落物最大持水率呈现上升趋势;当未分解枝量一定时,含同等量的未分解叶凋落物的最大持水率大于含同等量的半分解凋落物(图5—7)。
图5 青冈林拟合的不同分组凋落物最大持水率与实际对比
图6 青冈+光蜡树林拟合的不同分组凋落物最大持水率与实际对比
根据线性方程得到的化香树+密花树林凋落物层最大持水量和最大吸水速率均最小,且其半分解、未分解枝和未分解叶的最大吸水速率显著小于其他2种林分(p<0.05),但未分解叶的最大持水量显著高于青冈+光蜡树林的未分解叶(p<0.05);青冈林凋落物层的最大持水量和最大吸水速率均最大,且最大吸水速率显著高于化香树+密花树林(p<0.05)(表6)。
表6 不同林分类型凋落物层最大持水量和最大吸水速率
图7 化香树+密花树林拟合的不同分组凋落物最大持水率与实际对比
3 讨论
3.1 凋落物组成对凋落物持水性能的影响
凋落物的持水率和吸水速率与浸水时间呈对数函数和幂函数关系,即在浸水前2 h内凋落物持水率迅速增加,之后趋于平缓,直至饱和。这主要是因为浸水初期凋落物死细胞间或不同凋落物组分表面水势差较大,吸水迅速;随着时间的推移,凋落物内外水势差降低,吸水速率放缓[16]。本研究中不同分组凋落物持水率和吸水速率与浸水时间的关系与前人研究结果一致[2-3,7-9],表明凋落物组成差异不影响其持水过程函数关系的性质。
本研究显示凋落物的持水率和吸水速率与未分解叶的量呈极显著正相关,与未分解枝的量呈极显著负相关,与半分解层的量无显著关系。这与凋落物表面对水分子的吸附和凋落物内部空隙对水分的储存有关[17]。未分解叶具有较大的比表面积[18]和叶片间的结构空隙,有利于水分子的吸附和持留。相反,未分解枝的比表面积和可持水的结构空隙要小得多,其比例的增加会降低凋落物持水率和吸水速率。这与刘忠玲等[13]的研究结果类似。半分解凋落物的比表面积和结构空隙与分解程度有关,持水率和吸水速率受组成比例及其分解程度的共同影响,从而未表现出相关性。因此,当凋落物层总蓄积量差别不大时,未分解叶和枝的组成比值将可以作为判断其持水性能的有效指标,其值越大,持水性能越好。
本研究中不同林分类型间大致以纯叶的持水性能最好,纯枝的持水性能最差。这与王金悦等[2]和刘玉国等[4]半分解高于未分解的研究结果不同。这可能是因为本研究将持水性能低的未分解枝分离,未分解叶优越的持水性能得以突显。不同的组成凋落物持水性能在林分类型间也存在差异,如青冈林2∶8∶0持水率最低,青冈林0∶8∶2和化香树+密花树林2∶8∶0亦显著低于其他组成,达到最大持水率所需的时间3种林分类型的凋落物组成也各不相同。
此外,根据单一组分凋落物持水特性拟合得到的不同混合组成凋落物最大持水率总体小于实际最大持水率,两者符合线性函数关系(R2≥0.681)。这可能是凋落物各组分的混合会增大凋落物彼此间的空隙,如枝可起到混合态的骨架作用,增加了对水分的吸持。这与王忠禹等[17]对刺槐凋落物的研究结果类似。根据二者的关系方程,可以得到一种估算凋落物持水性能的较简便方法。
3.2 不同林分类型凋落物水文效应
凋落物的水文效应与凋落物蓄积量成正相关[8-9],而蓄积量由凋落物的输入动态与凋落物分解动态共同决定。群落的物种组成对凋落量和分解速率具有显著的控制作用[19],很大程度决定了凋落物蓄积量。林冠层的透光度也会对凋落物分解产生影响,从而影响蓄积量。本研究中3种林分凋落物层的总蓄积量及最大持水量均无显著差异,可能与采样时间在凋落高峰刚过时、受分解的影响较弱有关。但3种林分类型凋落物层的水文效应的大小顺序与郁闭度有较好的对应关系,一定程度上支持了前人[20]的研究结果。此外,凋落物组成特别是未分解枝和未分解叶对持水性能有显著影响,由此使得不同林分类型凋落物层的水文效应也出现差别。本研究中青冈林凋落物层的最大吸水速率显著高于化香树+密花树林,主要是青冈林的未分解枝高吸水速率的贡献(表6)。但是,不同林分类型的物种组成会影响凋落物的性质,从而使凋落物组成对其持水性能的影响会因林分类型而异。因此,应对更广泛的林分类型进行凋落物持水性能对比研究和规律总结,并在评估凋落物层的水文效应时,综合考虑凋落物量及其组成情况[13]。
3.3 展望
本研究选取了一个时间节点的凋落物层,且部分林分类型果的凋落量较少,未开展不同比例凋落果的持水性能对比研究。但果的性质和分解速率与枝、叶不同,对于果凋落量较大的林分应对其水文效应开展进一步研究。浸水试验的结果只反映了饱和条件下的吸水效应,自然降水的强度、过程复杂多变,开展的原位观测或者模拟试验可以更清晰地揭示凋落物对降水的调控规律。此外,凋落量和降水均存在明显的季节动态,开展不同季节的对比并与降水的季节变化相结合,有助于更好地理解凋落物层的水文效应及其季节动态。
4 结论
(1)凋落物持水率和吸水速率与未分解叶的量呈极显著正相关,与未分解枝的量呈极显著负相关,二者的比值是反映凋落物层水文效应的有效指标。
(2)凋落物的持水率与浸水时间呈对数函数关系,吸水速率与浸水时间呈幂函数关系,不受林分类型及其凋落物组成比例的影响。
(3)亚热带喀斯特林分凋落物持水性能由蓄积量及其组成共同决定,科学评价凋落物层的水文效应需综合考虑蓄积量及其具体组成的影响。