杨良辰, 张春茹

(1.河北省涉县职业技术教育中心 河北 涉县056400; 2.河北木兰围场国有林场管理局, 河北 围场 068450)

枯落物层和土壤层对森林截留降水、减缓降雨速度和涵养水分具有很重要的作用[1]。木兰围场位于沿坝地区，森林多以针叶为主，其中华北落叶松(Larixprincipis-rupprechtii)、油松(Pinustabulaeformis)为主要组成部分，阔叶则以山杨(Populusdavidiana)、白桦(Betulaplatyphylla)为主[2-3]。森林是生态系统的重要组成部分，其很重要的一个功能是水源涵养，减少水土流失，不同森林组成在水源涵养方面的能力也是不一样的，具有一定的差异[4-6]。该地区的研究主要集中林分的结构、生长特性、物种多样性、森林各层次的水源涵养能力，但是对林分水源涵养综合能力进行评价的研究则相对较少。本文利用熵权法将沿坝地区3种典型森林类型的枯落物层与土壤层水源涵养能力进行量化，能够更加直观地进行评价和比较，以便于为该地区的森林生态效益评价提供一定的理论依据[7-8]。

1 研究区概况

围场县位于承德境内，地理坐标为东经116°32′—118°14′，北纬41°35′—42°40′。与内蒙古接壤，处于华北山石区的沿坝地区，是山区与坝上的过渡地带，海拔750～1 829 m，年平均温度较低，为-1.5～4.8℃,年无霜期90～125 d，年均降水时间分布很不均匀，集中性比较强，降雨量为380～560 mm，土壤的种类也较多，总共有143个土种，植物种类资源也是非常丰富，野生的种子植物达到了793种，蕨类植物22种。此次的试验地点在围场县木兰林管局的北沟林场。

2 研究方法

2.1 枯落物层数据测定

2.1.1 枯落物生物量 2016年8月在沿坝地区的北沟林场，选择华北落叶松纯林、油松纯林、针阔混交林3种林分类型设置标准的样地，样地的大小为30 m×30 m，并对标准地的基本情况进行了调查(表1)，其中郁闭度的测定是通过在样地内机械地设置100个样点，在各样点位置上进行抬头垂直昂视，判断该样点是否被树冠覆盖，统计被覆盖的样点数，该点数与样点数的比值则是林分的郁闭度。而枯落物生物量的测定则是在选取的标准地里选取样地的4个角和中心位置设置5块枯落物小样方，大小为50 cm×50 cm，枯落物测定分未分解层和半分解层，以下枯落物的数据都是均值。

表1 3种林分类型标准地概况

注：Ⅰ代表华北落叶松纯林，简称落叶松林；Ⅱ代表油松纯林；Ⅲ代表针阔混交林,下表同。

2.1.2 枯落物持水量和吸水速率 利用室内浸泡法对3种林分的枯落物持水量和持水速率进行测定，先测量枯落物的厚度，根据枯落物分解的程度来进行区分，划分为未分解层与半分解层，然后要快速对其鲜重进行称量，此后带回实验室进行烘干，对其重量进行再次称重，最后将需要测定的枯落物在水中进行浸泡，在0.5,1,2,4,6,8,10,24 h对枯落物的重量变化要进行测定，从而得出枯落物的持水量、持水速率和最大持水率测定枯落物的持水量[9-11]。利用下式进行计算：

C=(m1-m2)/m2×100%

(1)

S=(m3-m2)/m2×100%

(2)

Wm=(Rm-R0)M

(3)

式中：C为枯落物自然含水量(%)；m1为样品鲜质量(g)；m2为样品烘干质量(g)；S为饱和持水率，m3为样品浸水24 h后的质量(g)；Wm为枯落物的最大拦蓄量(t/hm2)；Rm为最大持水率(%)；R0为平均自然含水量(%)；M为枯落物蓄积量(t/hm2)。

2.1.3 枯落物有效拦蓄量 枯落物对降雨的实际拦蓄量的测定，可以通过有效拦蓄量(modified interception)来进行计算,即:

W=(0.85Rm-Ro)M

(4)

式中:W为有效拦蓄量(t/hm2);Rm为最大持水率(%);Ro为平均自然含水率(%);M为枯落物蓄积量(t/hm2)。

2.2 土壤层物理性质和土壤入渗过程测定

按照不同的土层深度取样，即土壤剖面法，采用土壤烘干法和环刀法测定需要的土壤特性。其中土壤持水力S=10000hp,式中：S为土持水力(t/hm2)；h为土壤层厚度(m)；p为非毛管孔隙度(%)。

采用双环法测定水分入渗土壤的过程，初渗率=最初入渗时段内渗透量/入渗时间；平均渗透速率=达到稳渗时的渗透总量/达到稳渗时的时间。利用Kostiakov模型来对水入渗土壤的过程来拟合，在f=at-b中：f为t时间时的瞬时入渗速率；t为入渗时间；a，b则是指试验资料拟合的参数。

2.3 水源涵养综合评价

本研究为了能够将枯落物层和土壤层的水文效应更加直观地进行比较评价，利用熵权法对林分各因子进行量化，在一个标准系统下进行计算，对3种林分类型的水源涵养能力进行综合评价，主熵的概念源于热力学，是对系统状态不确定性的一种度量，通过熵值法得到各个指标的信息熵，信息熵越小，信息的无序度越低，其信息的效用值越大，指标的权重越大[12-13]。

(1) 依据以下公式计算出各个指标的熵值，记为Hi，各指标具体熵值，计算熵值公式为：

(5)

式中:Hi为熵值；k为常数1/lnn；fij为在第j个评价指标下第i个林分Hi的贡献度；i表示林分类型数；j表示评价指标数。

(2) 依据以下公式计算各个指标的权重，记为wi，计算权重公式为：

(6)

式中：wi为权重；Hi为熵值；i为林分类型数。

3 结果与分析

3.1 枯落物生物量

枯落物的存量多少是由很多因素决定的，决定枯落物存量的是林分枯落物的进入量和分解量，而林分枯落物的进入量与林分树种组成、郁闭度和林下植被生长情况都是有关的。由表2可以看出，针阔混交林的枯落物生物量最大，为17.19 t/hm2,其次是油松纯林，大小为16.9 t/hm2，最小的是落叶松纯林，为15.86 t/hm2。此外，3种林分的未分解层和半分解层的组成比例也是有差异的，其中落叶松纯林和油松纯林的未分解层的比例都是大于半分解层，而针阔混交林的半分解层的比例是大于未分解层的，针阔混交林未分解层所占比例最小，占总储量的46.42%，落叶松纯林则最大，占总储量的54.54%，而油松居中为52.6%。

表2 3种林分类型枯落物生物量

3.2 枯落物水文效应

3.2.1 枯落物最大持水量 3种林分类型中针阔混交林的最大持水量最大，为37.59 t/hm2,其次为油松纯林和落叶松纯林，分别为27.78 t/hm2和26.04 t/hm2；3种林分类型的最大持水率的范围163.55%～216.55%，依次为针阔混交林分>油松纯林>落叶松纯林。

表3 3种林分类型枯落物最大持水量和最大持水率

3.2.2 枯落物持水过程 枯落物在浸泡持水过程中有一定的规律，在0.5 h内是吸水最快的阶段，随着时间的推移枯落物的持水量不断增加，最后趋于稳定，但是不同的林分枯落物的持水速率也是有差距的，枯落物的这一过程与降雨时枯落物的拦蓄过程有一定的相似性，在降雨的初期,枯落物的拦蓄能力比较强，但是随着枯落物的持水量增加，枯落物的吸持能力降低。枯落物的未分解和半分解层的持水量也有一定的差距，达到饱和的时间也不相同(表4)。对0.25～24 h之间3种类型林分枯落物各层持水量与浸泡时间的关系进行回归分析,得出该时间段内持水量与浸泡时间之间存在如下关系(表5):

Q=alnt+b

(7)

式中:Q为枯落物持水量(g/kg);t为浸泡时间(h);a为方程系数;b为方程常数项。

表4 3种林分类型枯落物层持水过程 g/kg

3.2.3 枯落物吸水速率 通过试验表明，3种类型林分的枯落物的吸水速率具有一定的规律，在0.25～1 h内，枯落物的吸水速率是最大的时候，随后速率出现了较大的下降，在6 h左右时枯落物的下降速率趋于平缓，随着时间推移，枯落物的吸水速率趋于稳定，这种情况是因为随着浸泡的时间不断增长，枯落物的持水量达到了最大持水量，即枯落物的持水量趋于饱和，持水速率也会趋于稳定。

3.2.4 枯落物有效拦蓄量 林分类型下不同层次枯落物的拦蓄能力见表6，枯落物的累积量的大小排序为：Ⅲ>Ⅱ>Ⅰ，针阔混交林的枯落物累积量最大为17.19 t/hm2，落叶松纯林的最小为15.86 t/hm2，油松居中为16.9 t/hm2；最大持水量的大小排序为：Ⅲ>Ⅱ>Ⅰ，针阔混交林最大为37.59 t/hm2，落叶松纯林的最小为26.04 t/hm2，油松居中为27.78 t/hm2；有效拦蓄量大小排序为：Ⅲ>Ⅱ>Ⅰ，有效拦蓄量针阔混交林的枯落物累积量最大为27.65 t/hm2，落叶松纯林的最小为20.24 t/hm2，油松居中为21.71 t/hm2。最大持水率、最大拦蓄率和有效拦蓄量的大小排序都是一样的，其中针阔混交林的最大，油松其次，落叶松纯林最小。

表5 不同林分枯落物层持水量、持水率与浸泡时间关系

图1 3种林分类型枯落物持水量、持水速率与浸泡时间关系曲线

3.3 土壤层的水文效应

3.3.1 土壤层的物理性质 由表7可以看出，3种林分类型土壤容重的总体差距不大，其中落叶松纯林的最大，为1.31 g/cm3，油松纯林的次之，为1.21 g/cm3，针阔混交林的最小，为1.19 g/cm3；总孔隙度的排序为：针阔混交林(54.15%)>油松纯林(49.31%)>落叶松纯林(41.87%)；而非毛管孔隙度与土壤持水力密切相关，其大小排序为：针阔混交林(5.12%)>落叶松纯林(4.89%)>油松纯林(4.69%)。土壤持水能力大小排序为：针阔混交林(363.52 t/hm-2)>落叶松纯林(388.51 t/hm-2)>油松纯林(267.33 t/hm-2)。

表6 不同林分枯落物层有效拦蓄能力

表7 3种林分类型土壤的物理性质

3.3.2 土壤层的渗透性 土壤渗透性是土壤物理性质的主要特征之一，它代表了入渗的速度，而且还是林分水源涵养的指标之一。土壤的渗透性越好，林分的水源涵养能力越强，土壤的流失也会比较小，地表的径流也会比较少，从而减少林分土壤被侵蚀的程度。从表8可以看出，3种林分土壤层的初渗速率差距比较大，大小排序为：针阔混交林(24.7 mm/min)>落叶松纯林(22.1 mm/min)>油松纯林(18.9 mm/min)；随着时间的不断推移，土壤的入渗速率会趋于稳定，各林分的稳渗速率大小排序为：针阔混交林(1 mm/min)>落叶松纯林(0.9 mm/min)>油松纯林(0.7 mm/min)；呈现这种情况的主要原因为：针阔混交林的土壤层的非毛管孔隙数量多，对于土壤的水分渗入非常有利，渗透性也会比较强，而油松纯林土壤层的非毛管孔隙较少，不利于土壤层水分的渗入，渗透性也就比较差。3个林分达到稳定入渗的时间是不同的，针阔混交林时间最长为32 min，油松纯林最小为25 min，落叶松纯林居中为29 min。通过对入渗速率与入渗时间进行线性回归分析，发现二者之间存在较好的幂函数关系：

f=at-b(R>0.96)

式中：f为入渗速率(mm/min)；a，b为常数；t为入渗时间(min)。

表8 3种林分类型土壤渗透的速率及渗透数学模型

3.4 水源涵养能力综合评价

3.4.1 不同森林类型评价指标数据表 本研究采用熵权法，选择枯落物层和土壤层的指标来进行综合评价，其中枯落物层水文效应指标选择3个，包括枯落物蓄积量、有效拦蓄量和最大持水量；土壤层水文效应选择5个评价指标，包括土壤容重、毛管孔隙度、非毛管孔隙度、土壤持水力、初渗速率,总计8项指标，综合8项指标进行权重赋值来评价沿坝地区3种不同森林类型水源涵养能力(表9)。

表9 3种不同林分类型评价指标数据

注：A1为枯落物蓄积量，A2为枯落物有效拦蓄量，A3为枯落物最大持水量，A4为土壤容重，A5为土壤毛管孔隙度，A6为土壤非毛管孔隙度，A7为土壤持水力，A8初渗速率，下表同。

3.4.2 构造标准化指标矩阵并计算权重值 根据第j个属性下第i个方案Ai的贡献度，即先求出各列的和，然后用每行的数值比上列和，形成新的标准矩阵(表10)，再根据公式(5)计算出各指标的熵值，最后利用公式(6)计算出各林分类型的权重值(表11)，利用熵值法计算得出的权重为枯落物蓄积量占1%，枯落物有效拦蓄量占19%，枯落物最大持水量占28%，土壤容重占2%，土壤毛管孔隙度占14%，土壤非毛管孔隙度占1%，土壤持水力占19%，初渗速率占16%，故我们在对林分水源涵养进行综合评价的时候更多要考虑林分枯落物的最大持水量、枯落物有效拦蓄量、土壤持水力和初渗速率，3种林分指标的具体值与其权重的乘积之和为其综合评价值，这样得出3种林分类型综合评分排序为：针阔混交林(16.968 9)>油松纯林(12.854 4)>落叶松纯林(12.103 8)。

表10 各指标标准化矩阵

表11 各指标权重值

表12 水源涵养能力评价

4 结 论

(1) 枯落物层最大持水量：针阔混交林>油松纯林>落叶松纯林，3种林分中针阔混交林的持水潜力是最大，落叶松纯林相对是最小的；有效拦蓄量：针阔混交林>油松纯林>落叶松纯林，说明3种林分中真实拦蓄量最大的为针阔混交林，最次的落叶松纯林，出现这种差异主要是因为林分枯落物层的生物量不同，因此三种林分中最优的水源涵养林为针阔混交林；持水量与浸水时间的回归方程为Q=alnt+b(R2>0.97)，持水速率与浸水时间的回归方程为V=Ktn(R2>0.94)。枯落物持水测定采用了浸泡法，虽然在一定程度上反映了林分的真实拦蓄量，但是对于特定降雨情况的拦蓄能力还需要进一步的研究。

(2) 3种林分类型总孔隙度、非毛管孔隙度和土壤的初渗速率最高的均为混交林，高于其他两种纯林，说明混交林的地理环境更有利于林木的生长和涵养水源。非毛管孔隙度与土壤持水力密切相关，其大小排序与土壤持水能力相同，大小依次排序为针阔混交林、落叶松纯林、油松纯林，3种林分土壤层的初渗速率差距比较大，大小依次排序为针阔混交林、落叶松纯林、油松纯林；3个林分达到稳定入渗的时间是不同的，针阔混交林时间最长为32 min，油松纯林最小为25 min，落叶松纯林居中为29 min；入渗速率与入渗时间回归方程为：f=at-b(R>0.96)。

(3) 利用熵权法计算得出的权重大小排序为：枯落物最大持水量(28%)>枯落物有效拦蓄量(19%)=土壤持水力(19%)>初渗速率(16%)>土壤毛管孔隙度(14%)>土壤容重(2%)>枯落物蓄积量(1%)=土壤非毛管孔隙度(1%)，3种林分类型综合评分排序为：针阔混交林(16.968 9)>油松纯林(12.854 4)>落叶松纯林(12.103 8)。针阔混交林为最优的水源涵养林，其在保持水土、涵养水源方面功能最强。