第4l卷第7期 东北林业大学学报 Vo1.41 No.7 2013年7月 JOURNAL OF NORTHEAST FORESTRY UNIVERSITY Ju1.2013 不同诱导改造后大兴安岭蒙古栎低质林 土壤养分的灰色关联评价 ) 曾翔亮 董希斌 高 明 (东北林业大学,哈尔滨,150040) 摘要对大兴安岭蒙古栎低质林进行不同面积的带状和块状皆伐后,种植樟子松、西伯利亚红松、落叶松进 行诱导改造.采用灰色关联分析法对各诱导改造样地和对照样地的土壤养分建立综合评价模型。结果显示:和对 照样地相比.诱导改造后各样地的土壤养分在总体上均有所上升;块状改造比带状改造更能有效提高林地土壤养 分.其中20 mx20 m林窗改造方式最有利于土壤养分的积累,最适宜大兴安岭蒙古栎低质林的改造;采用灰色关联 分析法对土壤养分建立评价模型,不但简单可靠,而且模型具有通用性,稍加修改便可对其他区域的土壤养分进行 比较评价。 关键词 蒙古栎低质林;诱导改造;土壤养分;灰色关联分析法 分类号s714.8:s756.5 Grey Correlation Evaluation of Soil Nutrients to Quercus mongolica Low-Quality Forest in Daxing’an Mountains after Diferent Induced Transformation/Zeng Xiangliang,Dong Xibin,Gao Ming(Northeast Forestry University,Har- bin 15oo40)//Journal of Northeast Forestry University.一2013,41(7).一48~52 An evaluation was conducted to do the clear—cutting in Quercus mongolica low.quality forest of Daxing’8.1"1 Mountains with different area of the stirp and block and perform the induced transformation by planting Mongdica Litv,P/nus koraien— sis and Larix gmelinii.The soil nutirent evaluation model in the induced transformation zone and the check zoon were as- tablished by using the grey correlation analysis.Compared with the check zoon.the soil nutrient in the induced transforma- tion zone can rise in genera1.The block transformation can improve the land soil nutrient of low-quality forest effectively rather than the stirp transformation and the 20 mx20 m forest gap is the best way to accumulate soil nutrient.Therefore, the block transformation is the most suitable way to transforiB the Ouerc ̄mongolica lOW.quality forest in Daxing’an Mountains.It is simple and reliable to establish the soil nutirent evaluation model by using grey correlation analysis.The model Can compare and evaluate soil nutrients of other areas with minimal modifications. Keywords Quercus mongolica lOW.quality forest:Induced transformation;Soil nutirent;Grey correlation analysis method 土壤作为植物生长的重要介质.不仅为植物的 壤养分的评价主要侧重于定性评价和单因素评价, 生长提供了物理支撑,而且为植物生长提供了水分、 导致评价结果粗糙,可用信息量少,而传统的定量化 微生物和必要的养分,是生态系统中物质和能量循 描述则存在较大的主观随意性,已经无法满足人们 环的重要场所…。因此,土壤养分不但直接影响植 对土壤研究的需要[5]。灰色关联分析法是采用关 物的生长发育,而且也对森林群落内植物种类的分 联度来量化研究系统内各指标的相互联系、相互影 布格局具有重要影响,但由于受到人为因素和自然 响与相互作用的一种方法,若两指标参数矩阵构成 因素的共同作用.不同地区的土壤养分具有高度的 的空间几何曲线越接近.则关联度越大[6 ]。由于 空间变异性_2]。周莉等[3]对长白山阔叶红松混交 灰色关联分析法简单、高效,因此,在水质评价、大气 林的不同采伐年限迹地的土壤养分质量分数进行研 质量评价等环境质量评价领域[8 ]得到广泛应用, 究,发现采伐年限对迹地土壤养分质量分数具有重 但在土壤养分评价方面却鲜见报道。 要影响。张泱等[4]对小兴安岭低质林不同采伐强 大兴安岭林区的总经营面积达8.35×10 hm , 度改造后土壤养分的动态变化进行研究,发现改造 森林资源丰富,是我国重要的林业基地之一,木材贮 后的样地土壤养分高于对照样地.且在改造后初期 量约占全国的一半。但由于20世纪在强烈的自然和 土壤养分上升,然后下降。但目前国内外学者对土 人为因素共同干扰破坏下.林分结构成分大面积缺 失,生物多样性大幅下降,形成了大片质量低下、林 1)“十二五”农村领域国家科技计划课题(2Ol2BAD22B02O2— 相残败、林下土壤侵蚀严重的低质林[1 。文中以大 3)、林业公益性行业科研专项(201004043)。 第一作者简介:曾翔亮,男,1989年7月生,东北林业大学工程 兴安岭蒙古栎(Querc ̄mongolica)低质林为研究对 技术学院,硕士研究生。 通信作者:董希斌.东北林业大学工程技术学院,教授。E—mail: 象.探讨不同诱导改造方式对林地土壤养分的影响, xibindong@sina.COrn。 并尝试利用灰色关联分析法对各样地的土壤养分进 收稿日期:2013年5月2日。 责任编辑:任俐。 行综合评价,旨在为蒙古栎低质林今后的更新和培 第7期 曾翔亮等:不同诱导改造后大兴安岭蒙古栎低质林土壤养分的灰色关联评价 49 育提供参考依据。 1试验区概况 试验区设在黑龙江省大兴安岭林区的加格达奇 林业局(5O。34 9.17”~50。34 32.00”N.124。23 47.8”~ 124。24 35.1”E)所属翠峰林场174林班。该地区属 于温带性气候,春秋气候分明,冬季漫长并且寒 冷干燥,夏季时间较短,年平均气温为一1.2 oC,大于 l0℃的年积温只有1 806 oC。年平均降水量为494.8 cm,主要集中在6—8月份,无霜期为85~130 d;土 壤以暗棕壤和棕色针叶林土为主,土壤平均厚度为 22 cm;该地区海拔为370~420 m,地势平缓,坡向东 南,坡度多在15。以下;试验区林分类型为蒙古栎低 质林,乔木平均树高在9 m左右,郁闭度约为0.3, 下层木主要为胡枝子,草本主要为苍术。 2材料与方法 2.1样地设置 在加格达奇林业局翠峰林场174林班的蒙古栎 低质林进行带状(图1)和块状皆伐(图2),并通过 更新造林进行诱导改造。其中,带状改造试验区包 括4条顺山带,面积依次为6 mx200 m(S1)、10 m× 200 m(S2)、14 mx200 m(S3)和18 mx200 m(S4), 每条皆伐带沿山坡平均分成A、B、C 3段进行更新 造林.造林树种分别为樟子松、西伯利亚红松、落叶 松,栽植苗木时遵循以下原则:与相邻保留林带距离 1 m.株行距为2 mx1.5 m。 图1带状试验区设置图 G1 G2 G3 G4 G5 G6 图2块状试验区设置图 块状改造试验区共3组,每组由6个沿横坡方 向排列且面积不同的林窗试验区组成。6个试验区 的面积分别为:5 mx5 m(G1)、10 mX10 m(G2)、15 mx15 m(G3)、20 mx20 m(G4)、25 mx25 m(G5)、30 mx30 m(G6)。在3组林窗改造试验区内分别栽植 西伯利亚红松、樟子松和落叶松。栽植苗木时遵循以 下原则:与原有林分边缘间隔1 m左右。株行距配置 为1.5 mx1.5 m。同时,在蒙古栎低质林试验区中 设置一个未采伐的样地作为对照。 2.2试验方法 测定时间选在2012年的6月底(改造后第3 a),在对照样地和带状改造试验区的每条采伐带沿 山坡上、中、下按“S”型分别选取5个采样点。在块 状改造试验区随机选取3个采样点.每个样点取土 壤剖面为0~10 cm的土壤1 kg带回实验室.鲜土在 实验室做自然风干处理,然后去除杂质并研磨过筛, 用于分析土壤养分质量分数,分析方法见表1… 。 表1森林土壤分析方法 2.3数据处理 将测得的数据录入Excel 2010,进行基本的处 理后,再导入Matlab7.0中对矩阵进行计算。 2.4评价方法 2.4.1确定决策矩阵 由n个样地的m个评价指标实测值组成的集 合即为决策矩阵 。 2.4.2初始化决策矩阵 由于各评价指标的量纲和量纲单位均有所不 同.为了消除其对评价结果造成的影响,应该对各评 价指标进行无量纲化处理[1 ,从而得到初始化决策 矩阵 。 2.4.3计算灰色关联系数 理想对象矩阵 为: S={S } ( =1,2,…,m)。 式中:s 为初始化后的决策矩阵 中第i行的最大 值。 已知初始化后的决策矩阵 和理想对象矩阵S 后,就可以利用公式(1)计算出各自的关联系数r minminl¥i--X l+Amaxmaxl¥i- ̄f I — l亡 Si-习 I +Am—axm ax 8l i-l 。 ( ) 式中:A为分辨系数,其取值范围为0~1,其值只影 东北林业大学学报 第41卷 响样地关联度的大小,而不会影响各样地关联度的 的土壤养分指标的实测平均值见表2。从表2中可 排列顺序,一般取0.5; 一 表示s 与 在点 = 知,各改造带的土壤有机质质量分数与对照样地 1,2,…,n的绝对差;minminIs I为因素 1,2,…, (CK)相比,均有不同程度地下降,总体趋势是随着 m在点 =1,2,…,n的最小绝对差,也称二级最小 带宽增加,有机质质量分数下降(S4带除外);土壤 差;maxmaxI Si- ̄ l为因素i=1,2,…,m在点 =1, 全磷的质量分数和有机质类似,总体趋势也是随着 2,…,n的最大绝对差,也称二级最大差。 带宽增加,质量分数下降,且除S1带略高于对照样 2.4.4确定评价指标权重 地(CK)外,其他改造带的全磷质量分数均有不同程 由于土壤养分受各指标的影响程度不一样,因 度地下降;而土壤全氮正好相反,各诱导改造带全氮 此,需要对不同的影响指标赋予不同的权重。 的质量分数均比对照样地高(CK),且随着带宽增 2.4.5计算灰色关联度 加,土壤全氮质量分数升高(s2带除外);各诱导改 已知灰色关联系数和指标权重后,即可根据公 造带土壤全钾和土壤全氮类似。总体趋势是随着带 式(2)计算出各样地的关联度b : 宽增加,质量分数升高,但S1带、S4带与对照样地 相比,土壤全钾质量分数有所降低;在碱解氮和有效 = (Wi ̄rq)(,:1’2'… 。 (2) 磷方面.各诱导改造带的质量分数均不同程度地高 于对照样地,但和带宽没有明显的相关性;而速效钾 3结果与分析 却正好相反.其质量分数与带宽呈显著正相关,但和 3.1传统定量分析 对照样地(CK)相比,均有不同程度地下降。 蒙古栎低质林各诱导改造带、林窗和对照样地 表2样地土壤养分实测值 在蒙古栎低质林的林窗改造中,土壤有机质的 对照样地相比,面积较小的林窗(G1、G2、G3)有所 质量分数为15.47~4O.68 g・kg~,其中G1、G3、G5 升高,面积较大的林窗(G4、G5、G6)有所降低;在速 林窗的有机质质量分数与对照样地相比有所降低, 效钾方面,和对照样地相比,除G5林窗的质量分数 其他林窗则有不同程度地升高;在土壤全氮方面,除 有所降低外,其他林窗均有不同程度地升高。 G4林窗低于对照样地外,其他林窗均有不同程度地 3.2灰色关联分析评价 升高,其中最高的是G3林窗(13.75 g・kg );各林 3.2.1确定决策矩阵 窗改造中,全磷的总体趋势是随着林窗面积增加,其 选择合适的土壤养分评价指标是土壤养分评价 质量分数下降,其中除G2、G3林窗高于对照样地 的基础和关键,在遵循生产性、主导性和稳定性三原 外,其他林窗的全磷质量分数均有所下降;全钾的质 则的基础上l ,同时考虑大兴安岭的土壤特征并参 量分数除G4林窗低于对照样地外,其他林窗均有 考有关专家的指导建议.文中选取对土壤养分和植 所升高,且总体上和林窗面积呈负相关;各林窗改造 物生长影响较大且差异显著的土壤有机质、全氮、全 中,碱解氮的质量分数为428.07 595.57 mg・ 磷、全钾、碱解氮、有效磷和速效钾7个土壤养分指 ,均远远高于对照样地和诱导改造带;各林窗改 标作为评价指标[1 。由表2可知,本研究中m=7, 造有效磷的质量分数为l3.45 18.67 mg・kg~.和 rt=11。决策矩阵 为: 第7期 21.46 8.60 1.20 X= 24.33 曾翔亮等:不同诱导改造后大兴安岭蒙古栎低质林土壤养分的灰色关联评价 19.31 8.94 1.21 18.62 8.25 1.06 18.57 9.40 1.13 20.08 10.54 1.11 16.62 31.22 19.68 40.68 15.47 25.30 10.11 1.04 29.60 10.31 1.62 37.75 13.75 1.45 29.81 6.88 1.07 20.76 l2.75 0.90 25.50 10.82 1.06 25.49 23.39 26.86 29.43 22.99 291.12 316.40 366.03 341.21 341.21 428.07 390.84 465.29 595.57 521.12 446.68 18.42 46.81 24.97 32.11 26.08 37.46 20.74 42.91 23.95 45.03 18.67 52.45 18.67 51.95 21.03 51.26 13.45 59.03 15.68 35.70 15.99 49.61 3.2.2初始化决策矩阵 式中: 表示第 个样地的第i种土壤养分指标初 始化后的值; 表示第 个样地的第i种土壤养分指 标的实测值; 表示 在第i种土壤养分指标上的 最大值,即决策矩阵 中第i行的最大值。 (3) 目前,无量纲化处理的方法有很多 14-15],由于 本研究中的指标均为效益型指标,因此,采用公式 (3)对决策矩阵进行初始化处理: = /xm( =1,2,…,e;rj’=1,2,…,n) 因此,初始化后的决策矩阵 为: O.53 0.47 0.46 0.46 0.49 0.41 0.77 0.48 1.00 0.38 0.62 0.63 0.65 0.6O 0.68 0.77 0.74 0.75 1.00 0.50 O.93 0.79 0.74 0.75 0.65 0.70 0.69 0.64 1.00 0.90 0.66 O.56 0.66 = 0.64 0.62 0.7l 0.78 0.61 0.78 1.00 0.79 0.55 0.68 0.68 0.49 0.53 0.61 0.57 0.57 0.72 0.66 0.78 1.0o 0.88 0.75 0.71 O.96 1.0o 0.80 0.92 0.72 O.72 0.81 0.52 0.60 0.61 0.79 0.54 0.63 O.73 0.76 0.89 0.88 0.87 1.oo 0.60 0.84 3.2.3确定灰色关联判断矩阵 根据公式(1)进行计算,可得到灰色关联判断 矩阵R: 本研究中的理想对象矩阵 为: S =[1 1 1 1 1 1 1]。 0.40 0.37 0.36 0.36 0.38 0.34 0.57 0.38 1.00 0.33 0.45 0.45 0.47 0.44 0.49 0.57 0.54 0.55 1.00 0.38 0.81 0.59 0.55 0.55 0.47 0.51 0.50 0.46 1.00 0.75 0.48 0.41 0.47 = 0.47 0.45 0.52 O.58 0.44 0.59 1.00 0.60 0.41 0.49 0.49 0.55 0.38 0.40 0.45 0.42 0.42 O.52 0.47 0.59 1.00 0.710.51 0.88 1.00 0.60 0.79 0.52 0.52 0.62 0.39 o.44 0.44 O.60 0.40 0.46 0.53 0.57 0.74 0.72 0.70 1.OO 0.44 0.66 3.2.4确定评价指标权重 3.2.5确定灰色关联评价矩阵 目前,常用的权重确定方法有层次分析法、熵值 法和变异系数法等 1 ,本研究采用变异系数法来确 定各土壤养分指标的权重。采用变异系数法可以直 接利用土壤养分各影响因子所包含的信息,通过计 算得到各土壤养分影响因子的权重,是一种客观赋 根据公式(2)进行计算,可得到各个样地关联 度b .组成的集合即为灰色关联评价矩阵曰: CK S1 s2 s3 S4 Gl G2 G3 G4 G5 G6 口=[O.47 0.49 0.52 O.48 0.51 0.51 0.67 0.63 0.70 0.51 0.52]。 根据灰色系统理论中关联度分析原则,由于 权的方法[1 。各土壤养分指标的变异系数计算公 式如下: = “理想对象”的质量是系统中质量“最高”的,如果被 评价样地土壤养分的关联度越大.则其与“理想对 ( 1,2,…,m)。 (4) 象”越接近.表明其土壤养分越高[1 。通过灰色关 联评价矩阵可知:各诱导改造样地的关联度为0.48~ 0.70,均大于对照样地(0.47),表明各诱导改造样地 的土壤养分均优于对照样地;在各诱导改造带中,s2 带的关联度最高,但在所有的改造样地中,关联度最 高的是G4林窗(0.70),说明20 ex20 rm林窗改造方 式中: 是第i种土壤养分指标的变异系数,即标准 差系数; 是第i种土壤养分指标的标准差; 是第 i种土壤养分指标的平均值。 各土壤养分指标的权重计算公式为: (i=1,2,…,m)。 (5) 式比其他改造方式更有利于土壤养分的积累。 将表2中的数据按公式(4)和(5)可计算出各 土壤养分指标的权重系数,因此,各指标的权重矩阵 W为: '',=[0.22 0.13 0.12 0.12 0.15 0.14 0.12]。 4结论与讨论 各诱导改造带的土壤有机质、全磷质量分数与 带宽表现为负相关,且和对照样地相比,各诱导改造 52 东北林业大学学报 第41卷 带的有机质和全磷质量分数均有所降低;各诱导改 造带的土壤全氮、全钾、速效钾质量分数与带宽表现 用了已有的评价指标信息,避免了定性评价和传统 定量评价主观随意性大、结果粗糙的缺点,评价的结 果与实际情况是比较相符的,说明采用灰色关联分 为正相关,但与对照样地相比,各诱导改造带的土壤 全氮质量分数有不同程度地升高,速效钾质量分数有 所降低,全钾质量分数则表现为有升有降;在碱解氮 质量分数和有效磷质量分数方面,各诱导改造带均比 对照样地高,但和带宽没有明显的相关性。在林窗诱 导改造中.各改造样地的土壤有机质质量分数的变化 析法对大兴安岭蒙古栎低质林的土壤养分进行综合 评价是科学合理的,并且所建的模型具有通用性,稍 加修改便可对其他区域的土壤养分进行比较评价。 参 考 文 献 『1]徐建明,张甘霖,谢正苗,等.土壤质量指标与评价[M].北京: 科学出版社,2010. 规律不明显,但土壤全氮、全钾、碱解氮、速效钾质量 分数和对照样地相比,均有不同程度地升高,其中全 钾质量分数随着林窗面积增加而下降;全磷质量分数 的总体趋势是随着林窗面积增加而下降;各林窗改造 样地有效磷的质量分数和对照样地相比,面积较小的 林窗有所升高。面积较大的林窗有所降低。 土壤养分是林木生长的物质基础,二者相互作 用.相互影响:土壤养分状况的好坏直接影响着林木 生长速度及可持续性,而林木生长又反过来影响土 壤的养分状况[1 。研究采伐后林地土壤养分的分 布特征.对于了解土壤养分循环、森林土壤肥力和森 林生态系统的更新演替规律具有十分重要的意 义[20 ‘。。本研究采用灰色关联分析法对大兴安岭 不同诱导改造后蒙古栎低质林的土壤养分进行多指 标综合评价,结果显示:诱导改造后各样地的土壤养 分在总体上均有所上升,这主要是因为皆伐改造后 初期的林地透光率增加。土壤温度高于未采伐的对 照样地,森林微气候发生改变,林地上残留的大量采 伐剩余物变得易于分解,从而导致林地土壤养分增 加I2 ;林窗改造比带状改造更能有效提高林地土壤 养分,这和“林窗的全氮质量分数均高于改造带”的 试验结果一致,这可能是因为林窗改造后的森林微 气候最适宜采伐剩余物的分解,并且和林窗改造样 地不易发生水土流失有关:在所有的诱导改造方式 中,20 mx20 m林窗改造方式最有利于土壤养分的 积累,从土壤养分的角度看,最适宜大兴安岭蒙古栎 低质林的改造。 对于“部分信息已知,部分信息未知”的“小样 本”、“贫信息”的不确定性系统。灰色系统理论能够 精确的描述和理解_2 ,而在土壤养分评价中.只有 部分影响因子的信息已知,而其他影响因子的信息 是未知的,并且土壤养分的评价样本一般都比较少, 可见,灰色关联分析法在土壤养分评价上是非常合 适的。土壤养分的指标较多,单独从某个指标的角 度来评价各样地的土壤养分高低,是不合理的.本研 究采用灰色关联分析法和变异系数法对各样地的土 壤养分建立的综合评价模型,简单可靠,充分开发利 [2] 徐剑波,宋立生,彭磊,等.土壤养分空间估测方法研究综述 『J].生态环境学报,2011,20(8):1379—1386. [3] 周莉,代力民,谷会岩,等.长白山阔叶红松林采伐迹地土壤养 分含量动态研究[J].应用生态学报,2004,15(10):1771— 1775. [4] 张泱,宋启亮,董希斌.不同采伐强度改造对小兴安岭低质林 土壤理化性质的影响[J].东北林业大学学报,2011,39(11): 22—24. 『51姚荣江,杨劲松,陈小兵,等.苏北海涂围垦区耕层土壤养分分 级及其模糊综合评价[J].中国土壤与肥料,2009(4):16—20. [6]刘碉,薛喜成,.灰色关联分析法在土壤重金属污染评价 中的应用[J].安全与环境工程,2009,16(1):l5一l7. [7] 刘继明,宋启亮,李芝茹,等.大兴安岭白桦低质林生态功能评 价指标的灰色关联聚类分析[J].东北林业大学学报,2012,40 (8):112-l15. [8]翟国静.灰色关联分析在水质评价中的应用[J].水电能源科 学,1996,14(3):183—187. [9] 张磊,张波.灰色关联分析模型对大气环境质量预测的应用研 究[J].上海环境科学,1996,15(8):18—2O. [10]李勇,刘继明,董希斌,等.大兴安岭林区低质林成因及改造 方式[J].东北林业大学学报,2012,40(8):105—107. 『11] 国家林业局.LY/T 1210—1275—1999 森林土壤分析方法 [s].北京:中国标准出版社,2000. [12] 叶宗裕.关于多指标综合评价中指标正向化和无量纲化方法 的选择[J].浙江统计,2003(4):24—25. [13]李洪斌,王翠红,廖超林,等.城市边缘带土壤肥力特征及质 量评价:以长沙市东郊为例[J].农业现代化研究,2008,29 (2):219-222. [14] 武月荣,苏根成.灰色关联度法在蒙古高原北部草原土壤质 量评价中的应用[J].内蒙古师范大学学报:自然科学汉文 版,2008,37(6):775—779. [15]樊红艳,刘学录.基于综合评价法的各种无量纲化方法的比 较和优选:以兰州市永登县的土地开发为例f J].湖南农业科 学,2010(17):163—166. [16]戴西超,张庆春.综合评价中权重系数确定方法的比较研究 [J].煤炭经济研究,2003(11):37. [17] 赵宏,马立彦,贾青.基于变异系数法的灰色关联分析模型及 其应用[J].黑龙江水利科技,2007,35(2):26—27. [18]张磊,苏芳莉,郭成久,等.灰色关联分析在不同生态修复模 式土壤质量评价中的应用[J].沈阳农业大学学报,2009,40 (6):703—707. [19]许明祥,刘国彬.黄土丘陵区刺槐人工林土壤养分特征及演 变[J].植物营养与肥料学报,2004,10(1):40—46. [2O]周莉,代力民,谷会岩,等.长白山阔叶红松林采伐迹地土壤 养分含量动态研究[J].应用生态学报,2004,15(1O):1771— 1775. [21]程先富,史学正,于东升,等.丘陵山区林地土壤养分状况研 究[J].水土保持学报,2003,17(2):28—3O. [22]刘美爽,纪浩,董希斌.诱导改造对大兴安岭低质林土壤理化 性质的影响[J].林业科学,2012,48(5):67—71. [23]刘思峰,谢乃明.灰色系统理论及其应用[M].北京:科学出 版社,2010.