刘 宁,丁昌俊,李 波,丁 密,苏晓华,黄秦军
(中国林业科学研究院林业研究所,北京 100091)
杨树(Populusspp.)为杨柳科杨属(Populus)植物的总称,属于速生阔叶树种,其用途广泛,既是北半球重要的速生工业、民用材原料树种,又是进行水土保持、防风固沙和生态绿化的主要速生树种[1]。杨树速生、高产、易繁殖、木材易加工、适应性和用途广泛,是温带地区种植的首选树种,根据第九次森林资源结果,在我国其人工林面积已达757.07万hm2,占全国人工林总面积的13.25%,木材产量占全国木材总产量近30%,对国家林业建设和地域经济发展的贡献巨大,并取得了显著的社会效益和生态效益[2-6]。黑杨派(Section aigeiros)是杨树中的高生产潜力树种,具有极高的育种价值和经济价值,目前,大部分人工栽培的杨树都源于黑杨派[7-8]。大量实践表明,黑杨派中美洲黑杨(Populus deltoidesMarsh.)及其与欧洲黑杨(P. nigraL.)的杂种后代欧美杨(P.×euramericana)已成为我国杨树规模和商业化栽培的主体[9],对于人工用材林发展具有不可替代的作用。
由于生长周期长,林木的产出、适应性、木材性质等因素在育种时必须要综合考虑,其中,对能够影响产量和木材性质的因素尤为重视[10]。密度是人工林合理结构的数量基础,对林木生长和木材质量具有重要影响,人工林中个体间的空间竞争随林龄的增长而变大,合理的密度能提供适当的环境条件使林木达到最佳生长,过大的林分密度对林木的生长产生抑制作用[11-12]。研究表明,密度不仅影响林冠特征、林下植被和林木生长,而且还影响木材质量[13-14]。西南桦(Betula alnoidesHamilt)生长受到密度影响,其中,密度与胸径生长呈显著负相关,与林分蓄积呈正相关,与单株材积呈负相关,且随林龄增长,不同密度林分间蓄积差异逐渐缩小[15]。造林密度对6年生尾叶桉(Eucalyptus urophyllaBlake)胸径、立木单株材积、冠幅及枝下高的影响达到极显著水平,对蓄积量、木材纤维宽度和导管分子长度的影响达到显著水平,对树高、木材气干密度和木纤维长度有一定的影响[16]。
鉴于国家木材资源将长期紧缺现状[14],选育和栽培适宜的速生高产优质欧美杨品种,以提高我国尤其是北方地区杨树人工林的生产力和木材质量非常重要。本文以2个国家审定欧美杨良种渤丰1号杨(P.×euramericana‘Bofeng 1’;良种编号:国S-SC-PE-002-2011)、渤丰3号杨(P.×euramericana‘Bofeng 3 hao’;良种编号:国S-SV-PE-003-2017)和欧美杨新品种京2杨(P.×euramericana(Dode)Guiner c1.‘J2’;品种权号:20 050 033)为试材[17-18],通过分析比较不同密度对不同林龄生长和木材性质的影响,确定合理造林密度,为欧美杨新品种的选育、北方杨树用材林特别是纸浆材林产量提高和木材品质优化提供有价值参考。
试验开展于辽宁省锦州凌海市原种场(41°17′ N,121°36′ E),年降水量400~600 mm,海拔17 m。2007年春季埋根造林,造林所用根为1年生苗木的根。试验采用4个密度处理(2 m×3 m、2 m×4 m、2 m×5 m和2 m×6 m),3个重复,林地四周栽植2行隔离行。实验材料为3个欧美杨品种:渤丰1号杨、渤丰3号杨和京2杨。
2007至2012年和2014年每年秋季使用测树围尺测量胸径,其中,使用2014年春季调查的胸径代替2013年秋的胸径数据;统计不同年份的保存率;木材性质的测定由东北林业大学在2017年完成。
使用Excel、SPSS25以及R3.6.0等统计软件对试验数据进行分析。对保存率进行整理并绘制连年变化图;对胸径进行方差分析;对3品种的胸径进行多重比较;对于木材解剖性质和木材化学成分,取3个品种的均值作为欧美杨人工林整体的木材性质进行多重比较。
2.1.1 保存率的连年变化 3个品种的保存率变化模式具有明显的差异(图1)。总体上看,京2杨的4个栽植密度均呈逐年下降的趋势,第4年开始逐年较大幅度下降,各造林密度间差异不明显;渤丰1号杨2 m×6 m造林密度的保存率始终高于87%,而2 m×5 m造林密度的保存率从第4年开始逐年小幅下降,较大造林密度(2 m×3 m和2 m×4 m)的保存率是从造林开始逐年小幅下降至第3年,第3至6年保存率稳定在70%左右,在第7年有小幅下降,第7、8年的保存率维持在55%左右;渤丰3号杨较小造林密度(2 m×5 m和2 m×6 m)的保存率变化较稳定,至第8年两密度保存率在75%左右,2 m×4 m造林密度的保存率大幅下降至第2年,后维持在54%~63%,而2 m×3 m造林密度的保存率,在第1至3年大幅下降,在第6至8年均小幅下降,最终保存率为41%。
2.1.2 胸径的方差分析及多重比较结果 对连续测量的胸径数据进行双因素方差分析(表1),结果表明:造林密度和品种对胸径生长均有极显著的影响(p<0.01);造林密度与品种的互作效应在第2至6年均达显著水平(p<0.05),其中,第3、4年达极显著(p<0.01)水平。
图1 保存率变化Fig. 1 Changes in preserving rates
表1 欧美杨胸径的双因素方差分析Table 1 Two-Way ANOVA of DBH of Populus×euramericana
对不同年份不同造林造林密度下各品种的胸径进行多重比较(表2),结果表明:不同年份间,各造林密度下,品种胸径的表现存在明显差异。2 m×3 m造林密度下,除第2年外,渤丰3号杨胸径最大且除第3年外显著(p<0.05)高于其他品种,第2年京2杨的胸径最大且显著高于其他品种;2 m×4 m造林密度下渤丰3号杨的胸径在第2至8年最大且除第2、第7年外显著高于其他品种;2 m×5 m造林密度下,京2杨在第2、3年胸径最大且显著高于其他品种,第4年各品种无显著差异,第5至8年渤丰3号杨胸径最大且显著高于京2杨;2 m×6 m造林密度下,第2年京2杨胸径最大且显著优于其他品种,第3年渤丰1号杨胸径最小且显著劣于其他品种,第4至8年渤丰3号杨胸径最大且显著高于其他品种。第8年试验结束时,2 m×6 m造林密度下渤丰3号杨的胸径最大,达25.01 cm;渤丰1号杨次之,达22.96 cm;京2杨最差,仅21.34 cm。
2.2.1 对木材物理性质的影响 由表3可以看出:除第3年外,欧美杨在同一林龄不同栽植密度下木材基本密度差异不显著;第3年2 m×4 m造林密度的欧美杨木材基本密度显著大于2 m×3 m造林密度,其余造林密度差异不显著。整体上看,欧美杨木材基本密度每年的总均值在第4、5年下降,除此之外呈逐年增长的趋势;各栽植密度中,随林龄的增长,欧美杨的木材基本密度存在增长的趋势,各造林密度下欧美杨第7、8年的木材基本密度均高于第1、2年。2 m×3 m、2 m×4 m、2 m×5 m和2 m×6 m造林密度下欧美杨木材基本密度首先出现显著性差异的年份分别为第4年、第3年、第7年和第7年。较小造林密度下,欧美杨的木材基本密度呈增长趋势,较大造林密度下欧美杨的木材基本密度,整体上是逐年增长的趋势,但其在第4至6年存在木材基本密度下降的现象;至第8年,各造林密度下欧美杨的基本密度为0.33~0.34 g·cm-3。
随着林龄的增大,欧美杨的微纤丝角均呈先增加后减小的趋势,至第8年,各造林密度下欧美杨的微纤丝角为15.00~15.65;只有第3年和第5年各造林密度间微纤丝角存在差异,在第3年,2 m×4 m和2 m×6 m造林密度的欧美杨微纤丝角显著大于2 m×3 m造林密度;在第5年,2 m×5 m和2 m×6 m的欧美杨微纤丝角显著大于2 m×3 m造林密度。只有2 m×4 m和2 m×6 m造林密度下欧美杨微纤丝角的逐年变化存在差异,2 m×3 m和2 m×6 m密度处理的欧美杨微纤丝角随林龄的增长差异不显著,2 m×3 m、2 m×4 m、2 m×5 m和2 m×6 m造林密度下欧美杨微纤丝角开始明显下降年份分别为第3年、第4年、第6年和第6年(表3)。
表2 3品种不同林龄的胸径及多重比较Table 2 Multiple comparisons of 3 cultivars' DBH at different forest ages cm
欧美杨木材纤维的长、宽及长宽比在同一林龄时,各造林密度之间均差异不显著(表4)。至第8年,欧美杨木材的长宽比为49.65~52.83。随着林龄的增长,各造林密度下的欧美杨长宽比和纤维长均有显著增加,而纤维宽则没有显著的变化。长宽比总均值的排序为:2 m×3 m>2 m×4 m>2 m×5 m>2 m×6 m;纤维长总均值的排序为:2 m×6 m>2 m×3 m>2 m×5 m>2 m×4 m;纤维宽总均值的排序为:2 m×6 m>2 m×5 m>2 m×3 m>2 m×4 m。
2.2.2 对木材解剖性质的影响 欧美杨木材的壁腔比差异不显著,第6年时,2 m×3 m、2 m×4 m和2 m×5 m 密度的欧美杨壁腔比较低,2 m×6 m密度壁腔比相对稳定,其余年份壁腔比大部分在0.3以上(表5)。欧美杨木材的木纤维组织比量在同一林龄不同栽植密度间均差异不显著,除2 m×3 m造林密度外,其余造林密度第8年欧美杨的木纤维组织比量均显著小于造林初期;欧美杨木纤维组织比量较高,均在50%以上,木纤维组织比量木纤维组织比量总体上随林龄的增长呈减小的趋势(表5)。
2.2.3 对木材化学成分的影响 欧美杨木材的各化学成分在不同栽植密度下均差异不显著(表6)。欧美杨木材的纤维素含量为40.54%~42.76%,综纤维素含量均大于76%,欧美杨的酸溶木素含量均高于3.25%,总木素的含量在22%左右;2 m×3 m造林密度的酸不溶木素和总木素含量最高,2 m×5 m造林密度的水分和综纤维素含量最高,2 m×6 m造林密度的纤维素和酸溶木素含量最高。
表3 不同造林密度下的基本密度和微纤丝角Table 3 Basic density and microfibrillar angle of different planting densities
造林密度和欧美杨的品种特性与保存率的变化关系密切,研究其变化规律对欧美杨品种选择及生产的科学规划具有实际意义。不同造林密度下,3品种每年的保存率变化趋势差异显著,京2杨保存率在第1至8年保持中等幅度的持续下降,其第8年保存率最低;渤丰1号杨与渤丰3号杨保存率变化相似,均为较小造林密度的保存率优于较大造林密度,较大造林密度在第3年和第6年附近出现下降。整体看,2 m×6 m造林密度下渤丰1号杨的保存率最高,始终高于87%。保存率出现差异可能是因为渤丰3号杨和渤丰1号杨均是来自荷兰的欧洲黑杨145号(P.nigra)与河南的美洲黑杨65号 杨(P. deltoidescl. ‘55/65’ ×P.deltoidescl.‘2KEN8’)的杂交种,而京2杨属聚合型杂种欧美杨,但其遗传成分中欧洲黑杨仅占1/8,而美洲黑杨占7/8[17-18]。由遗传结构差异而带来的杂种优势[7,19-20],在本研究中极为明显:渤丰3号杨和渤丰1号杨既有欧洲黑杨的抗逆性,也获得了美洲黑杨的速生性,能适应当地冬季寒冷等环境条件而且有较好的生长表现;而京2杨主要成分是美洲黑杨,继承欧洲黑杨抗寒能力较弱,加之其大冠幅的特点,使京2杨生长处于逆境状态,这可能是其保存率低的重要原因。据报道[13-14,21],造林密度会影响杨树冠层特性、生物量、树体和叶片形态以及光合作用。本研究结果表明,欧美杨品种随着林龄的增长,其对生存空间的需求增加,当初始造林密度满足不了其生存需求时,空间上的竞争使保存率下降以增加其生存空间。人工林的造林密度与林分后期密度关系紧密,林木的生长受此影响,过高的造林密度会使林分中可利用的空间减少,从而加剧林木间的竞争,林木的死亡率上升、生长受到抑制,必需适当调整林分密度,才能促进林木生长发育,达到速生丰产的效果[11-12,15,22]。
表4 不同造林密度下欧美杨的纤维结构Table 4 The fiber structure of Populus×euramericana with different planting densities
造林密度和品种对欧美杨人工林生产力具有显著的影响。研究表明,密度主要与胸径生长相关显著[15-16]。本研究中,品种和造林密度对黑杨派品种胸径生长均有极显著(p<0.01)的效应;品种与造林密度的互作在第2至6年均有显著效应,而在第7和第8年互作效应不显著。说明品种对不同造林密度的响应不一,根据胸径表现筛选品种、造林密度有利于人工林的营造栽培。据悉,渤丰1号杨和渤丰3号杨是针对东北环渤海湾地区设计杂交的欧美杨造林品种,其综合表现优于108杨[18]。随林龄增长,首先郁闭的是造林密度较大的人工林,杨树胸径的生长会较早受到竞争而被抑制[23];而品种与造林密度的互作效应可以看做基因型与环境的互作效应,试验结果表明互作效应随林龄的增长会消失,这与相关的研究结果一致[24-25]。第8年试验结束时2 m×4 m和2 m×6 m造林密度下欧美杨各品种的胸径表现最好,但2 m×4 m造林密度下各品种的保存率较低。综合各品种保存率和胸径的表现,应当选择2 m×6 m造林密度来营造人工林。
表5 不同造林密度下的纤维形态Table 5 The fiber morphology of different planting densities
表6 不同密度处理下欧美杨的化学成分Table 6 Chemical composition of Populus×euramericana with different planting densities
造林密度对欧美杨品种木材物理和化学性质的效应不明显,第8年试验结束时,4种造林密度下欧美杨木材的解剖性质和化学成分均差异不显著,这可能是因为保存率的连年变化导致了生存空间的增加。鲍甫成等[26]认为,阔叶树木材密度会因林分中郁闭度的增加而减少。基本密度是木材在用作纤维原料时的重要性质,密度为0.4~0.6 g·cm-3的木材坚硬适中且易于制浆;木纤维长度不仅与木材的利用关系紧密,而且也是评价木浆的一项重要因素。纤维长度的增加会提高纸张的耐撕度、抗拉强度、耐破度和耐折度等性质;壁腔比对纸张的性能产生影响,壁腔比越小,纸张强度越大;上等造纸用材的标准是壁腔比小于1。高木纤维组织比量也是选择造纸材树种时的考虑因素[27]。纸浆材的木材纤维应细且长,长宽比大于30,且愈大愈好[28]。在试验的第8年,欧美杨平均木材基本密度为0.34 g·cm-3,平均微纤丝角为15.32,平均长宽比达51.15,平均纤维长达1 199.65 μm,平均纤维宽达23.75 μm,平均壁腔比为0.31,平均木纤维组织比量达56.58%。是易于制浆,并可以用于生产优质性能纸张的上等造纸用材。
适宜造林密度的确定应充分考虑培育材种和品种特性。本研究中,欧美杨品种人工林的培育主要是纸浆用材,因此,木材的化学成分倍受关注。结果表明,造林密度愈小,欧美杨木材纤维素含量愈高,这一趋势与姜岳忠等[28]、柴修武等[29]、曹福亮[30]对不同造林密度杨树的木材化学性质的研究结果相同。本研究中,不同密度下欧美杨的综纤维素含量间差异不显著,均大于76%;4种密度处理下的综纤维素含量排序为:2 m×3 m<2 m×6 m<2 m×4 m<2 m×5 m。总体上看,造林密度愈小,综纤维素的含量愈高。欧美杨的酸溶木素含量均高于3.25%,可以使用酸法制纸浆,总木素的含量在22%左右。酸不溶木素的含量均低于19.97%,这一结果与雷金选等[31]研究的8种杨树的材性结果略有差异。欧美杨木材的各化学成分在不同密度下均差异不显著(表7),2 m×4 m和2 m×6 m密度欧美杨的总木素含量和水分含量较低,2 m×5 m和2 m×6 m密度的纤维素含量较高,2 m×4 m和2 m×5 m密度的综纤维素较高。相比之下,造林密度较小的欧美杨木材化学成分更适宜纸浆的生产。
(1)较大的造林密度对欧美杨3个品种的保存率和胸径产生了抑制作用;较小造林密度下欧美杨的保存率更高。3个欧美杨派品种的适应性存在差异,渤丰1号杨和渤丰3号杨优于京2杨。(2)品种、造林密度、品种与造林密度的互作对欧美杨胸径的生长均有显著效应,品种与造林密度的互作效应对胸径的影响随时间增长而变弱。(3)造林密度对欧美杨的木材性质的影响不显著。随林龄的增大,木材基本密度、纤维长、纤维长宽比显著增大,而欧美杨木材的微纤丝角、木纤维组织比量显著减小,纤维宽的变异较小,其木材性质均符合纸浆材要求。(4)在大凌河流域平原及类似环境下,渤丰1号杨和渤丰3号杨更适合生产,且造林密度最少不能小于2 m×6 m。本试验中,渤丰3号杨有发展大径材的潜力,而高保存率的渤丰1号杨适合纸浆材人工林的营造。