大兴安岭低质落叶松林补植改造后土壤肥力的综合评价1)

2018-07-13 06:00韩贵杰唐亚森曲杭峰董希斌
东北林业大学学报 2018年6期
关键词:低质大兴安岭土壤肥力

韩贵杰 唐亚森         曲杭峰 董希斌

(大兴安岭新林林业局,新林,165023)   (森林持续经营与环境微生物工程黑龙江省重点实验室(东北林业大学))

土壤是森林生态系统的必要组成部分,它对植物的生长提供肥力,为生物的繁衍提供物质基础,土壤的改变会使森林生态系统发生改变,与此同时森林生态系统的变化又会对土壤造成影响[1]。土壤肥力的质量直接影响到林木的分布、生长和产量,因此,土壤肥力的研究对我国的森林可持续经营和林业生态安全具有重要的意义[2-5]。林地土壤肥力变化的研究,是当前国内土壤学领域的前沿课题[7-9],有关低质林土壤肥力的研究,已见诸多报道。吕海龙等[7]对小兴安岭进行带状和块状改造,发现马永顺林场的针阔混交林块状改造优于带状改造,且25 m2块状改造更有利于更新苗木的生长。张泱等[10]对不同采伐强度后小兴安岭低质林土壤肥力进行灰色聚类评价,发现当采伐强度为55%时的土壤肥力得到明显的改善。宋启亮等[11]对大兴安岭低质阔叶混交林进行不同改造模式的综合评价,得出10 m带状改造和100 m2块状改造明显优于其他改造模式。曾翔亮等[12]研究了带状改造和块状改造对大兴安岭低质林土壤肥力的影响,得出块状改造比带状改造更能有效提高林地土壤肥力,其中20 m×20 m块状改造方式最有利于土壤肥力的积累,更加适合大兴安岭蒙古栎低质林进行生态改造。毛波等[13]研究了诱导改造对大兴安岭低质山杨林土壤肥力的影响,分别种植落叶松,樟子松,西伯利亚红松,云杉,发现诱导改造后所有林地的土壤肥力均有所提升,其中樟子松林的诱导改造效果最佳。陈百灵等[14]通过对小兴安岭马永顺林场的针阔混交林进行研究,发现在改造后的第7年土壤肥力情况最佳,使土壤肥力情况明显改善,且所有顺山带改造以8 m改造带的改造效果最好,所有横山带改造以6 m改造带的改造效果最好。

目前国内外学者采用灰色关联度法、层次分析法、主成分分析法等对土壤肥力质量进行综合评价[15-17],主要从不同带宽带状改造、不同面积块状改造和不同采伐强度改造等方面研究对土壤肥力的影响,但关于补植改造对土壤肥力综合评价研究较少[18-20]。灰色关联分析能够很好的处理信息不完整的问题,尤其是对样本小且规律不明显的评价指标上,通过数据分析,得到与理想对象接近的最优方案[21-23]。本文以大兴安岭蒙古栎低质林为研究对象,对试验样地进行实地勘察,选取落叶松作为补植苗木,采用块状补植补造法进行改造,利用灰色关联分析法对补植改造后样地的土壤肥力进行综合评价分析,得出最佳的土壤肥力改造模式,即最佳的低质林补植改造密度,可为大兴安岭低质林补植改造提供参考和依据。

1 研究区概况

研究区设在大兴安岭新林林业局宏图林场,其在黑龙江大兴安岭山脉的东南坡,地理坐标为东经57°56′41.8″~124°5′78.2″,北纬51°36′37.15″~57°45′32″,属于低山丘陵地带。海拔高度介于713~745 m,林下土壤主要是暗棕壤和棕色针叶林土,土壤厚度15~30 cm,地势平缓,立地条件较好,坡度多在15°以下,平均坡度5°左右;无霜期介于85~130 d,年平均降水量大约494.8 mm左右,降雨多集中于夏季的6—8月,属温寒带大陆性季风气候,冬天气候干燥寒冷且漫长,夏季短暂,年平均气温-1.4 ℃,最髙气温37.3 ℃,最低气温-45.4 ℃。试验区内林分林分类型主要为落叶松(Larixdahurica),林分郁闭度0.4,补植苗木为樟子松,平均苗高1.2 m,灌木类型以杜鹃为主,盖度为7%,草木以越桔为主,盖度为55%。

2 研究方法

2.1 样地设置

于2016年6月经过野外实地勘察设计,在大兴安岭新林林业局宏图林场设立试验样地,改造样地编号记为BZ1~BZ6,对照样地为CK,其中7块样地的大小为50 m×50 m。在改造样地补植樟子松(Pinussylvestris),通过块状补植法对试验样地进行改造,分别在6个不同的改造样地补植不同密度的樟子松(Pinussylvestris)。BZ1~BZ6样地的补植密度分别为300、400、500、600、700、800株·hm-2,CK作为对照样地,不进行补植作业,在改造试验区相邻处,选取林分和立地条件接近的保留地作为对照样地CK,且改造样地分两排横山排列。大兴安岭低质落叶松林补植改造完成后,须对补植树种进行基本的管理和维护,并对其进行科学合理的抚育,当年的抚育工作主要有浇水、扩穴、扶正、培土、踏实、除草。

2.2 土壤肥力指标的收集与测定

2017年9月进行野外取样,在每块补植改造样地和对照样地上,按照“S”型混合采样法进行取样,每个样地选择5个土壤取样点,在每个取样点均取厚度为0~10 cm的土壤样本,按照四分法混合取样,每个土壤样本为1 kg。将土壤样本带回实验室,在室内土壤经过自然风干、研磨过筛,然后分析化学性质,土壤化学性质的测定方法如表1所示。土壤的物理性质采用环刀法进行测量,环刀容积为100 cm3。

表1 森林土壤化学性质测定方法

2.3 土壤肥力的综合评价方法

将理化性质实验所得的数据录入到Excel2010中,进行标准化处理后,运用灰色关联度分析法进行行综合评价,求出各个指标的权重,然后运用SPSS20.0计算出低质林补植改造后土壤肥力的灰色关联度,灰色关联度值越大,说明越适合大兴安岭低质落叶松林的补植改造。

3 结果与分析

3.1 补植改造对土壤肥力的影响

大兴安岭低质落叶松林补植改造后各个样地土壤物理指标的描述性统计分析见表2。BZ1、BZ2、BZ3、BZ4、BZ5、BZ6改造样地的土壤含水率均高于对照样地;随着补植密度的增大,各个改造样地土壤含水率无明显变化趋势。在各个改造样地中,BZ3改造样地土壤含水率91.46%最大,高于对照样地38.18%;BZ6改造样地土壤含水率61.03%最小,高于对照样地7.75%。

由表2可知:BZ1、BZ2、BZ3、BZ4、BZ5、BZ6样地,土壤密度均低于对照样地,差异性显著;随着补植密度的增大,各个改造样地土壤密度无明显变化趋势。在各个改造样地中,BZ3改造样地的土壤密度0.69 g·cm-3最大,低于对照样地0.12 g·cm-3。

由表2可知:BZ1、BZ2、BZ3、BZ4、BZ5、BZ6样地,最大持水量均高于对照样地,差异性显著(P<0.05);随着补植密度的增大,各个改造样地最大持水量无明显变化趋势。在各个改造样地中,BZ1改造样地的最大持水量143.67%最大,高于对照样地47.37%;BZ2改造样地的最大持水量106.23%最小,高于对照样地9.93%。

表2 补植改造后各样地土壤含水率、土壤密度和最大持水量实测值

注:表中数据为“平均值±标准差”;同列不同小写字母表示差异显著(P<0.05)。

由表3可知:BZ1、BZ2、BZ3、BZ4、BZ5、BZ6改造样地的非毛管孔隙度均高于对照样地,差异性显著;随着补植密度的增大,各个改造样地非毛管孔隙度呈现先增大后减小的趋势。在各个改造样地中,BZ4改造样地的非毛管孔隙度33.16%最大,高于对照样地22.38%;BZ1改造样地的非毛管孔隙度13.56%最小,高于对照样地2.78%;其中BZ1与BZ6改造样地的非毛管孔隙度相差不大。

由表3可知:BZ1、BZ3、BZ6改造样地的毛管孔隙度高于对照样地,经方差分析与对照样地差异性显著(P<0.05);其他改造样地的毛管孔隙度与对照样地,经方差分析与对照样地差异性不显著(P≥0.05)。随着补植密度的增大,各个改造样地毛管孔隙度无明显变化趋势。在各个改造样地中,BZ3改造样地的毛管孔隙度58.23%最大,高于对照样地2.67%;BZ5改造样地的毛管孔隙度53.78%最小,低于对照样地1.78%。

由表3可知:BZ2、BZ3、BZ4、BZ5的改造样地总孔隙度大于对照样地,经方差分析与对照样地差异性显著(P<0.05)。随着补植密度的增大,各个改造样地总孔隙度呈现先增大后减小的趋势。在各个改造样地中,BZ4改造样地的总孔隙度83.62%最大,高于对照样地19.59%;BZ1改造样地的总孔隙度67.28%最小,高于对照样地3.25%。

表3 补植改造后各样地非毛管孔隙度、毛管孔隙度和总孔隙度实测值 %

注:表中数据为“平均值±标准差”;同列不同小写字母表示差异显著(P<0.05)。

大兴安岭低质落叶松林补植改造后各个样地土壤有机质质量分数的描述性统计分析见表4。BZ1改造样地的有机质质量分数小于对照样地,经方差分析与对照样地差异性显著(P<0.05);BZ3、BZ4、BZ5改造样地的有机质质量分数大于对照样地,经方差分析与对照样地差异性显著(P<0.05)。随着补植密度的增大,各个改造样地有机质质量分数呈现先增大后减小的趋势。在各个改造样地中,BZ5改造样地的有机质质量分数32.87 g·kg-1最大,比对照样地高18.49 g·kg-1;BZ1改造样地的有机质质量分数11.49 g·kg-1最小,比对照样地低2.89 g·kg-1;其中BZ6改造样地的有机质质量分数14.78 g·kg-1与对照样地相差不大。

由表4可知:BZ5改造样地的pH值大于对照样地,经方差分析发现差异性不显著(P≥0.05);其他改造样地的pH值小于对照样地(BZ4改造样地除外),经方差分析与对照样地差异性不显著(P≥0.05)。随着补植密度的增大,各个改造样地pH值大体上呈现先增大后减小的趋势。在各个改造样地中,BZ5改造样地的pH值6.62最大,高于对照样地0.36;BZ3改造样地的pH值5.57最小,低于对照样地0.69;BZ4改造样地的pH值6.26与对照样地大小相等。

表4 补植改造后各样地有机质质量分数和pH实测值

注:表中数据为“平均值±标准差”;同列不同小写字母表示差异显著(P<0.05)。

大兴安岭低质落叶松林补植改造后各个样地全氮质量分数的描述性统计分析见表5。BZ1、BZ3、BZ4、BZ5、BZ6改造样地的全氮质量分数小于对照样地,经方差分析与对照样地差异性显著(P<0.05);BZ2改造样地的全氮质量分数小于对照样地,经方差分析与对照样地差异性不显著(P≥0.05)。随着补植密度的增大,各个改造样地全氮质量分数无明显变化趋势。在各个改造样地中,BZ2改造样地的全氮质量分数21.56 g·kg-1最大,低于对照样地1.37 g·kg-1;BZ4改造样地的全氮质量分数9.78 g·kg-1最小,低于对照样地12.15 g·kg-1。

由表5可知:BZ1、BZ2、BZ3、BZ4、BZ5、BZ6改造样地的全磷质量分数均小于对照样地,经方差分析与对照样地差异性显著(P<0.05)。随着补植密度的增大,各个改造样地全磷质量分数无明显变化趋势。在各个改造样地中,BZ2改造样地的全磷质量分数9.72 g·kg-1最大,低于对照样地1.74 g·kg-1;BZ4改造样地的全磷质量分数3.32 g·kg-1最小,低于对照样地6.78 g·kg-1。

由表5可知:BZ1、BZ2、BZ3、BZ4、BZ5、BZ6改造样地的全钾质量分数小于对照样地,随着补植密度的增大,各个改造样地全钾质量分数无明显变化趋势。在各个改造样地中,BZ2改造样地的全钾质量分数7.96 g·kg-1最大,低于对照样地0.84 g·kg-1;BZ1改造样地的全钾质量分数6.22 g·kg-1最小,低于对照样地2.58 g·kg-1。

表5 补植改造后各样地全氮质量分数、全磷质量分数和全钾质量分数实测值 g·kg-1

注:表中数据为“平均值±标准差”;同列不同小写字母表示差异显著(P<0.05)。

大兴安低质落叶松林补植改造后各个样地速效氮质量分数的描述性统计分析见表6。BZ1、BZ2、BZ3、BZ6改造样地的速效氮质量分数小于对照样地,经方差分析与对照样地差异性显著(P<0.05);其他改造样地的速效氮质量分数大于对照样地,经方差分析与对照样地差异性显著(P<0.05)。随着补植密度的增大,各个改造样地速效氮质量分数无明显变化趋势。在各个改造样地中,BZ5改造样地的速效氮质量分数456.58 mg·kg-1最大,高于对照样地248.30 mg·kg-1;BZ2改造样地的速效氮质量分数194.69 mg·kg-1最小,低于对照样地13.59 mg·kg-1。

由表6可知:BZ1、BZ2、BZ3、BZ4、BZ5改造样地的速效磷质量分数大于对照样地,经方差分析与对照样地差异性显著(P<0.05);BZ6改造样地的速效磷质量分数大于对照样地,经方差分析与对照样地差异性不显著(P≥0.05)。随着补植密度的增大,各个改造样地速效磷质量分数大体上呈现先增大后减小的趋势。在各个改造样地中,BZ3改造样地的速效磷质量分数194.55 mg·kg-1最大,高于对照样地140.33 mg·kg-1;BZ6改造样地的速效磷质量分数56.45 mg·kg-1最小,高于对照样地2.23 mg·kg-1。

由表6可知:BZ3、BZ4、BZ5改造样地的速效钾质量分数大于对照样地,经方差分析与对照样地差异性显著(P<0.05);BZ1、BZ6改造样地的速效钾质量分数小于对照样地,经方差分析与对照样地差异性不显著(P≥0.05);BZ2改造样地的速效钾质量分数大于对照样地,经方差分析与对照样地差异性不显著(P≥0.05)。随着补植密度的增大,各个改造样地速效钾质量分数呈现先增大后减小的趋势。在各个改造样地中,BZ3改造样地的速效钾质量分数42.47 mg·kg-1最大,高于对照样地24.90 mg·kg-1;BZ1改造样地的速效钾质量分数15.48 mg·kg-1最小,低于对照样地2.09 mg·kg-1。

表6 补植改造后各样地速效氮质量分数、速效磷质量分数和速效钾质量分数实测值 mg·kg-1

注:表中数据为“平均值±标准差”;同列不同小写字母表示差异显著(P<0.05)。

3.2 补植改造后土壤肥力的综合评价

3.2.1确定决策矩阵

设有n种补植改造模式,每种改造模式有m个评价指标,本研究中m=14,n=7,因此,决策矩阵X为:

3.2.2标准化决策矩阵

为了消除量纲对评价结果的影响,需要对决策矩阵X进行标准化处理,得到初始化决策矩阵X′:

本研究中采用极大化法和极小化法对土壤肥力的各指标进行标准化评价。正向评价指标采用最大化法,逆向评价指标采用最小化法。

3.2.3确定灰色关联判断矩阵

S={si}m×1,i=1、2、…、m。

式中:si为X′中第i行的最大值。由初始化后的决策矩阵X′可知,本研究中的理想对象矩阵S为:

计算理想对象2矩阵与其初始化决策矩阵对应的绝对差值,形成如下绝对差值矩阵:

i=1、2、…、n,j=1、2、…、m。

式中:rij为灰色关联系数,rij反映的是其他待决策的序列与相对于最佳方案序列的灰色关联度。λ为分辨系数,其取值范围为0~1,其值只影响各生态改造模式灰色关联度的大小,而不会影响各生态改造模式灰色关联度的排列顺序,一般取0.5。通过计算出不同密度补植改造后土壤肥力的评价指标,得到灰色关联度判断矩阵R如下:

3.2.4确定评价指标权重

在土壤肥力的评价指标中,不同的指标对土壤肥力的综合评价结果的影响不尽相同,需要对每个指标赋予不同的权重。使用变异系数法计算出大兴安岭加格达奇林业局低质落叶松林在不同密度补植改造模式的土壤肥力指标权重,得到各评价指标的权重矩阵W:

3.2.5计算灰色关联度

已知灰色关联系数和指标权重之后,根据如下公式进行计算,可得到补植改造后各个样地土壤肥力的关联度bj。

表7 补植改造后土壤肥力的灰色关联度

在灰色评价理论体系中,理想对象S的质量是评价体系中最高的,不同密度的补植改造后土壤肥力的灰色关联度越大说明与理想对象S越接近,说明其补植效果越有利于大兴安岭低质落叶松林土壤肥力的积累。

通过灰色关联分析法,对大兴安岭低质落叶松林补植改造初期的土壤肥力进行综合评价分析,灰色关联度从从大到小依次为:BZ5(0.726)、BZ4(0.671)、BZ3(0.667)、BZ2(0.580)、BZ6(0.570)、BZ1(0.503)、CK(0.472)。可知BZ5(0.726)改造样地的灰色关联度最高,表明BZ5改造样地补植改造效果最理想,700株·hm-2改造方式最适宜大兴安岭低质落叶松林的补植改造。

4 结论与讨论

大兴安岭低质落叶松林植改造后,BZ4、BZ5改造样地的有机质质量分数较大,其中BZ5改造样地的有机质质量分数最大,比对照样地高18.44 g·kg-1;其中BZ4、BZ5改造样地的pH值较大,以BZ5改造样地的pH值最大,比对照样地高0.36,BZ4改造样地的pH值与对照样地大小相等;各个改造样地的全氮质量分数小于对照样地(BZ2改造样地除外),经方差分析与对照样地差异性显著(P<0.05),BZ2改造样地的全氮质量分数最大,低于对照样地1.37 g·kg-1。所有改造样地的全磷质量分数均小于对照样地,经方差分析与对照样地差异性显著(P<0.05),其中BZ2改造样地的全磷质量分数最大,低于对照样地1.74 g·kg-1;所有改造样地的全钾质量分数均小于对照样地,经方差分析与对照样地差异性显著(P<0.05),BZ2改造样地的全钾质量分数最大。所有改造样地的速效氮质量分数,BZ5改造样地的速效氮质量分数456.58 mg·kg-1最大,高于对照样地248.30 mg·kg-1;所有改造样地的速效磷质量分数均大于对照样地,其中BZ3、BZ5改造样地的速效磷质量分数较大,以BZ3改造样地的速效磷质量分数最大,比对照样地高140.33 mg·kg-1;BZ3改造样地的速效钾质量分数42.47 mg·kg-1最大,高于对照样地24.90 mg·kg-1;随着补植密度的增大,各个改造样地速效钾质量分数呈现先增大后减小的趋势。

大兴安岭低质落叶松林补植改造后,各个改造样地的土壤含水率均高于对照样地,经方差分析与对照样地差异性显著(P<0.05),其中BZ3、BZ5改造样地的土壤含水率较大,以BZ3改造样地的土壤含水率最大,比对照样地高38.18%;各个改造样地的土壤容重均低于对照样地,其中BZ2改造样地土壤容重最小,低于对照样地0.45 g·cm-3;各个改造样地的最大持水量均高于对照样地,BZ1改造样地的最大持水量最大,高于对照样地60.25%;各个改造样地的非毛管孔隙度均高于对照样地,其中BZ4、BZ5改造样地的非毛管孔隙较大,以BZ4改造样地的非毛管孔隙最大,比对照样地高22.38%;BZ1、BZ3、BZ6改造样地的毛管孔隙度高于对照样地,经方差分析与对照样地差异性显著(P<0.05),其中BZ3改造样地的毛管孔隙度最大,高于对照样地2.67%;BZ2、BZ3、BZ4、BZ5的改造样地总孔隙度大于对照样地,其中BZ4、BZ5改造样地的总孔隙度较大,以BZ4改造样地的总孔隙度最大,比对照样地高19.59%。

通过灰色关联分析法,对大兴安岭低质落叶松林补植改造初期的土壤肥力进行综合评价分析,灰色关联度从从大到小依次为:BZ5(0.726)、BZ4(0.671)、BZ4(0.667)、BZ2(0.580)、BZ6(0.570)、BZ1(0.503)、CK(0.472)。可知BZ5(0.726)改造样地的灰色关联度最高,表明BZ5改造样地补植改造效果最理想,700株·hm-2改造方式更有利于土壤肥力的积累,当大于此值时,评价土壤肥力的灰色关联值呈现下降趋势。随着补植密度的增大,土壤肥力的灰色关联度先升高后降低,原因是对低质林进行补植改造后,随着林分密度的增大,补植改造样地的森林微气候发生了改变,有助于土壤肥力的积累,当补植密度过大时,补植样地土壤肥力的积累无法完整地提供植被对肥力的吸收时,造成土壤肥力的综合质量下降。

大兴安岭低质落叶松林改造后土壤肥力的综合评价除了与补植密度有关,还和立地条件、经济和社会等因素有关,这几个方面的问题还有待探讨。同时,本文只对大兴安岭低质落叶松林补植改造两年后的土壤肥力进行了研究,但土壤肥力的监控是一个漫长的过程,如果想更加全面的研究大兴安岭低质落叶松林改造后的土壤肥力,还需要连续并长期的观测与分析。

猜你喜欢
低质大兴安岭土壤肥力
大兴安岭不会忘记你
想离离不了,想过过不好——低质婚姻面面观
走出“低质婚姻”的雾霾
大兴安岭的铁道兵
2架增雨飞机为大兴安岭送雨
拉萨市土壤肥力变化趋势与改良策略
大兴安岭四季(四首)
不同类型低质林生物多样性的综合评价
江西省赣州市低质低效林成因与特征分析
安吉白茶园土壤肥力现状分析