0 引言
纤维细胞是树木次生木质部重要组成部分,在植物体结构中起机械支撑作用
[1],林木的水分储量受纤维性状的影响
[2],植物可以对木质部中的纤维细胞壁、纤维细胞管腔和轴向薄壁组织进行权衡提高环境适应性
[3],从而保证植物在恶劣环境下的生长发育。研究纤维细胞解剖特征在不同气候条件下的变化规律,为解析树木对气候变化的响应机制提供新支持。
温度是影响树木生长的重要因子
[4-5],低温事件会影响树木正常的生理生态活动。已有研究发现低温冻害显著降低油桐(
Vernicia fordii)幼苗的光合效率
[6]。随着自然低温胁迫的加强,楸树(
Catalpabungei)通过增加细胞膜流动性和提高细胞液浓度等抵御机制提高自身抗寒性
[7]。树木长时间生长在寒冷环境中,会逐步演化出御寒机制以减轻低温胁迫
[8]。
水曲柳((
Fraxinus mandshurica)是木犀科(
Oleaceae)白蜡树属(
Fraxinus)阔叶乔木,主要分布在我国东北地区
[9],是珍贵的三大硬阔之一
[10],具有重要的经济价值。作为东北地区阔叶红松林的重要伴生树种,水曲柳在增加生态系统固碳能力
[11]、促进群落演替和维护生态系统稳定
[12]方面有重要作用。近年来,木材需求不断上升,水曲柳大径级材数量急剧减少,天然林资源几近枯竭
[13],营建优质人工林是缓解木材供需矛盾的重要途径。东北地区气候寒冷,目前该地区水曲柳对低温的响应与适应策略还知之甚少。了解低温事件对水曲柳木材纤维解剖特征的影响,对于东北地区水曲柳选育具有重要意义。本研究以小兴安岭地区水曲柳子代测定林为研究对象,采用树轮年代学和木材解剖学方法,研究水曲柳纤维解剖特征与主要气候因子之间关系,明确低温对木材纤维解剖特征的影响,为该地区耐寒水曲柳选育提供科学依据。
1 研究区概况和研究方法
1.1 研究区概况
研究地点位于黑龙江省伊春市带岭区(47°10′—47°14′N,128°53′—128°55′E),小兴安岭南麓,海拔280~1 090 m,属于大陆性季风气候,地势西北高、东南低。年均气温1.4 ℃,年均降水量661.5 mm,无霜期115 d左右。地带性土壤为暗棕壤,垂直地带性明显。主要植被类型为针阔混交林,树种有红松(
Pinus koraiensis)、红皮云杉(
Picea koraiensis)、冷杉(
Abies nephrolepis)、水曲柳、胡桃楸(
Juglans mandshurica)、黄菠萝 (
Phellodendron amurense)、山杨(
Populus davidiana)和紫椴(
Tilia amurensis)等
[14]。
1.2 样本采集及测定
于2023年4月在黑龙江省小兴安岭地区水曲柳子代测定林进行样品采集工作。水曲柳子代测定林于2002年采种育苗、2004年营建,采用完全随机区组设计,5次重复。各区组中的每个家系为一个小区,每个小区10株,双行排列,共计52个家系。根据国际树轮库的取样标准,每个家系选取20株以上生长良好的样树,使用内径为5.15 mm的生长锥,在每株样木胸高处(1.3 m)采集2根树芯
[15]。确保钻过髓心,树芯编号带回试验室后进行固定、风干、打磨。交叉定年后读取年轮宽度,制作年轮宽度标准年表。每个家系筛选出与主序列相关性最高的3根来自不同采样树的南向样芯,共156株树芯进行木质部纤维解剖特征的测定及分析。
样芯在75 ℃恒温水浴锅中软化,用GSL-1滑走式切片机切片
[16]。用番红乙醇试剂对切片进行染色,树胶封片,完全干燥后用XSP-6X6A光学显微镜在20倍物镜下观察并拍照,如
图1所示。用Adobe Photoshop CC 2019裁剪处理图像,搭载Image-Pro Plus 6.0软件,使用Roxas软件测定样芯年轮宽度(ring width,RW,mm)、木纤维细胞数量(fiber cell number,FN)、平均木纤维细胞面积(mean fiber cell area,MFA,μm
2)、木纤维细胞总面积(total fiber cell area,TFA,μm
2)、木纤维细胞密度(fiber density,FD,个/mm
2)、木纤维细胞占比(mean percentage of fiber cell area within xylem,RFTA, %)、木纤维细胞壁厚度(overall mean thickness of all fiber cell walls,CWTall,μm)、总体平均木纤维细胞壁增强指数(square of the total-to-bound ratio,(T/B)
2)数据。所采集的树芯能够有效判读的年份为2008—2022年。
1.3 气候数据
由CRU(
https://crudata.uea.ac.uk/cru/data/hrg/)获取月平均气温(
T)、月平均最低气温(
Tmin)、月平均最高气温(
Tmax)和月降水量(
P),分辨率为0.5°×0.5°。从中国气象局国家气象科学数据中心(
http://data.cma.cn/)获得月相对湿度(relative humidity,RH)。小兴安岭地区夏季高温多雨、冬季寒冷干燥,1962—2022年月平均气温为-22.20~20.91 ℃,月平均最低气温为-28.56~15.65 ℃,月平均最高气温为-15.89~26.24 ℃,月降水量在4.16~158.95 mm,如
图2(a)所示。年降水量有缓慢增高的趋势,如
图2(b)所示。年平均气温有显著增高的趋势,如
图2(c)所示。年相对湿度有下降的趋势,如
图2(d)所示。
1.4 数据分析
采用样条函数去除树木生长中与年龄有关的生长趋势并进行标准化处理,将年轮宽度与主要纤维解剖参数进行相关分析。鉴于树木生长对气候因子的响应存在时间滞后效应
[17],为确保气候数据与树木生长周期的精准匹配,本研究选取2007—2022年间每年的上一年10月至当年9月的气候数据作为分析基数;同时纤维解剖特征原始数据标准化后进行主成分分析(PCA),最终采用Pearson相关分析构建主成分与气候数据的关系。低温年采用相对(百分位)阈值法
[18]和累计距平法
[19]确定,通常极端气温事件百分位数值多取90%(10%)或95%(5%),本研究时间区间较短,因此以85%(15%)来确定百分位数值。利用IBM SPSS 27软件和R语言进行统计分析,利用Origin Pro 2021、R 4.2.3、Adobe Illustrator CC 2017和Excel软件绘制图表。
2 结果与分析
2.1 水曲柳年轮宽度和纤维解剖特征
小兴安岭地区水曲柳子代测定林在20 a内年轮宽度和木质部纤维解剖特征表现出不同的变化趋势,如
图3所示。年轮宽度在6~20 a表现出先升后降的变化趋势,在13 a(2015年)达到最高值,为4.18 mm,较6 a(2008年)时的年轮宽度增长了49.76%。木纤维细胞数量和木纤维细胞总面积变化趋势与年轮宽度相似,在11、13 a时达到最大值,分别达到了2.294×10
3个和9.05×10
5 mm
2,较6 a时分别上升48.17%和116.32%。木纤维细胞密度与总体平均木纤维细胞壁增强指数均呈现逐年下降的趋势,20 a(2022年)时木纤维细胞密度和总体平均木纤维细胞壁增强指数分别为3.87×10
-3个/mm
2和0.9,比6 a时分别下降了30.27%和20.35%。平均木纤维细胞面积和木纤维细胞壁厚度则呈现逐年增加的趋势,平均木纤维细胞面积年际增长速率较快,而木纤维细胞壁厚度的年际增长较为缓慢,20 a时平均木纤维细胞面积和木纤维细胞壁厚度较6 a时分别增加了66.26%和3.57%。
由
图4可知,水曲柳年轮宽度与木纤维细胞数量、木纤维细胞密度和木纤维细胞总面积极其显著正相关(
R=0.74~0.99,
P<0.001)。木纤维细胞数量与木纤维细胞密度和木纤维细胞总面积极其显著正相关(
R=0.82~0.98,
P<0.001)。木纤维细胞密度与木纤维细胞总面积极其显著正相关(
R=0.79,
P<0.001),与木纤维细胞占比极显著正相关(
R=0.73,
P<0.01)。平均木纤维细胞面积与总体平均木纤维细胞壁增强指数显著负相关(
R=-0.63,
P<0.05)。木纤维细胞总面积与木纤维细胞占比显著正相关(
R=0.56,
P<0.05)。木纤维细胞占比与总体平均木纤维细胞壁增强指数极显著负相关(
R=-0.72,
P<0.01)。木纤维细胞壁厚度与总体平均木纤维细胞壁增强指数极其显著正相关(
R=0.90,
P<0.001)。
2.2 气候因子对水曲柳纤维解剖特征的影响
水曲柳纤维解剖参数可以概括为2个具有代表性的主成分,如
图5所示。第1主成分(PC1)和第2主成分(PC2)的贡献率分别为52.83%和33.27%,2个主成分总共解释了水曲柳子代测定林木质部纤维解剖参数变异86.10%的总方差,表明2个主成分已成功提取足够的信息,能够充分反映原始变量的特征。
PC1主要由木纤维细胞总面积、木纤维细胞数量、木纤维细胞密度和年轮宽度构成。4个参数因子载荷量绝对值在0.87以上,均与PC1负相关,PC1是构成水曲柳木纤维细胞数量和径向生长量的复合因子。PC2主要由总体细胞壁增强指数、平均木纤维细胞面积和木纤维细胞壁厚构成。总体平均木纤维细胞壁增强指数和木纤维细胞壁厚度与PC2负相关,平均木纤维细胞面积与PC2正相关,总体平均木纤维细胞壁增强指数和木纤维细胞壁厚度与平均木纤维细胞面积负相关,PC2是构成水曲柳木纤维细胞面积和厚度的复合因子。
由水曲柳纤维解剖参数主成分与主要气候因子相关性如
图6所示。由
图6可知,PC1与上一年10、12月以及当年1、5、6、9月的温度因子负相关,与2、3、7月的温度因子正相关。PC1与8月平均温度正相关,与8月最低温度显著正相关(
P<0.05)。PC2与上一年10、11月以及当年1、2、5月温度因子负相关,与6、7月的温度因子正相关,其中与7月平均温度和最低温度显著正相关(
P<0.05)。
PC1与3、4月和6—9月的降水量负相关,与5月份降水量正相关。PC1与3、8月相对湿度正相关,与5—7月和9月相对湿度负相关。PC2与上一年11月和6—8月降水量负相关,与6、8月相对湿度显著负相关(P<0.05),与4—5月和9月的降水量和相对湿度正相关。小兴安岭地区水曲柳子代测定林木纤维解剖特征受温度影响较大,尤其是生长季最低温。
2.3 低温事件对水曲柳纤维解剖特征的影响
由2008—2022年小兴安岭地区年平均气温、年平均最高气温和年平均最低气温距平值可知,2010年和2012年发生了低温事件,使水曲柳的径向生长减小,如
图7所示。
由
图8可知,在低温事件下,水曲柳年轮宽度、木纤维细胞数量、平均木纤维细胞面积、木纤维细胞总面积、木纤维细胞占比极显著低于非低温年(
P<0.01),分别下降了24.6%、16.9%、18.3%、31.5%、8.6%。低温年木纤维细胞密度和总体平均木纤维细胞壁增强指数显著高于非低温年(
P<0.05),分别升高11.6%和9.2%,其中,木纤维细胞密度极显著高于非低温年(
P<0.01)。水曲柳通过减少木纤维细胞数量,提高木纤维细胞密度这一相对保守的生存策略来适应低温胁迫。
3 讨论
3.1 水曲柳年轮宽度与木纤维解剖特征的关系
年轮是木质部解剖特征的宏观表现,能够反映出树木生长的速度,其特征的变化是对自身生长状态和外界环境的综合响应
[20]。小兴安岭地区水曲柳子代测定林在幼龄林阶段,水曲柳的年轮宽度与纤维细胞数量、木纤维细胞总面积年际变化趋势一致,且存在极显著的正相关关系,表明纤维细胞数量和木纤维细胞总面积是影响年轮宽度的重要因子。幼龄林阶段水曲柳年轮宽度年际变化呈现先升后降的趋势,与杨树的年轮宽度径向变异规律相似
[21-22],符合林木径向生长速率沿半径梯度遵循典型的单峰模式
[23],是受到林木个体生长发育趋势控制的表现。在发育初期林木生长旺盛,形成层细胞分裂活跃,能够增殖形成较多的纤维细胞。此外,幼龄林木生长受环境因子的限制较小
[24],能够满足细胞生长发育的物质需求,加速细胞发育成熟,进而提高木纤维细胞的总面积,使得年轮宽度在生长发育初期有较大的年际增长。林木的生长发育过程中,植物体叶片和根系的功能逐渐成熟
[25],能够吸收更多的水分和矿质养分
[26],合成更多的光合同化物并用于细胞的发育,平均纤维细胞的面积逐渐增大,木纤维密度减小,但细胞数量减少和平均面积增大存在差异,导致木纤维占比的减少。在林木的发育过程中,纤维细胞厚度不断增加,与平均纤维细胞面积表现出一致的变化趋势,但与年轮宽度和平均细胞面积的相关性较弱,对其影响较小。
3.2 气候因子对小兴安岭地区水曲柳纤维解剖特征的影响
林木的径向生长是一个复杂的过程,其中木质部细胞经历了细胞分裂、增大、细胞壁增厚、木质化和程序性细胞死亡5个阶段
[27],受气候因子的影响
[28],其中气温和水分状况对林木生长和木质部解剖特征的影响尤为显著
[29-30]。小兴安岭地区水曲柳径向生长和纤维解剖特征受生长季气温和水分状况影响。以往的研究表明,气温是影响高寒地区的树木生长的主要气候因子
[31-32],对木材解剖特征的影响存在一定的迟滞性
[33],上一年10月的温度与年轮宽度呈现负相关关系,可能是低温增加了树体对生长季积累养分的消耗,进而对来年的木材径向生长产生不利影响。8月的月均最低气温与年轮宽度和纤维细胞数量显著正相关,最低温度的升高能够使光合酶等植物生理活动更长时间处于活跃状况,有利于营养物质的积累。初春时期的降水与水曲柳径向生长负相关,这是由于形成层活性对温度比较敏感
[34],林木在此阶段逐渐解除休眠
[35],降水增加容易对植物产生冻害,会抑制纤维细胞的分裂和发育,最终造成林木当年径向生长量的减少。
低温是限制植物生长的重要因素
[36],严重时可以导致林木的死亡。小兴安岭地区的低温事件对水曲柳子代测定林木纤维解剖特征存在显著影响,抑制了水曲柳的径向生长,使其较非低温年下降了32.6%,纤维细胞数量和平均纤维细胞面积减小,木纤维细胞密度有所增加。低温对树木根、茎和叶的发育存在消极的影响,既会直接损伤树木组织与器官,还会改变环境状况,间接影响树木的生长发育
[37]。为适应低温环境,植物会消耗储藏的营养物质用于维持正常的生命活动,进而阻碍了林木的生长
[38]。已有研究表明东北地区水曲柳对极端低温具有较高的抵抗力,能够适应一定的低温环境
[5],这可能与木质部解剖特征的改变具有一定的关系。低温导致形成层细胞分裂能力下降
[39],木纤维细胞数量减少,同时平均纤维细胞面积的下降提高了木纤维细胞密度。木纤维密度是木材密度的重要组成部分,木材密度与植物对抗重力、风力等环境力的支撑作用有关
[40],高的木材密度与植物维持高存活状态有关
[41],较高的木纤维密度一定程度上有利于林木在低温环境下的生存。减少纤维细胞数量、增强木材密度是水曲柳抵抗低温事件的保守生存策略。
4 结论
运用木材解剖学方法探究小兴安岭地区水曲柳子代测定林对低温事件响应的研究中发现,年轮宽度是由木质部纤维解剖特征的动态变化调控的,其中木纤维细胞数量、木纤维细胞总面积与年轮宽度年际变化趋势一致,显著正相关,是年轮宽度的决定因子。在幼龄林阶段,水曲柳木纤维细胞密度和总体平均木纤维细胞壁增强指数逐年下降,平均木纤维细胞面积和木纤维细胞壁厚度逐年增加。低温事件使水曲柳径向生长受到抑制,使年轮宽度、木纤维细胞数量、平均木纤维细胞面积、木纤维细胞总面积和木纤维细胞占比显著下降32.6%、20.3%、22.4%、45.9%、9.4%,而木纤维细胞密度升高11.6%。低温事件发生时,水曲柳会采取减少木纤维细胞数量、提高木纤维细胞密度这一相对保守的生存策略,研究结果为小兴安岭地区耐寒水曲柳选育提供科学依据。
国家重点研发计划项目(2021YFD2200303)