李宗善,陈维梁,焦 磊,李兆林, 2,王 聪, 2,高光耀,2,孙婧雅,白应飞,温润泉

1 中国科学院生态环境研究中心 城市与区域生态国家重点实验室,北京 100085

2 陕西黄土高原地球关键带国家野外科学观测研究站,西安 710061

3 陕西师范大学 地理科学与旅游学院,西安 710119

4 延安市退耕还林工程管理办公室,延安 716000

5 延安市宝塔区南泥湾国有生态林场,延安 716000

刺槐(Robiniapseudoacacia)由于生长速率快,且具有较强的环境适应性和抗干旱胁迫能力,在黄土高原地区得到了广泛的种植,是黄土高原生态恢复建设的重要人工林类型之一。然而,刺槐林这种速生、高耗水的林型在属于干旱半干旱区的黄土高原广泛种植和推广造成了其生态适应性与当地气候条件出现了脱节,造成土壤水分过度消耗而出现土壤干层,导致土壤供水能力大大下降,局地生存条件恶劣,黄土高原延安以北地区的人工刺槐林已出现大面积生长衰退问题[1]。因而对黄土高原人工林刺槐林对环境因子的生态适应性研究具有较为重要的生态学意义,目前树木蒸腾耗水特征[2]、树木生态型和生活型谱特征[3]、植物化学计量参量[4─5]、植物功能形状[6]和生长适宜生态区划[7]等方法对该地区人工槐林的生态适应性均有一定的研究。

树木年轮准确地记录了树木年度生长量信息,具有气候信号明确、数据连续性强、样本地分布广泛等优点[8—9],已经被广泛地应用于树木生长对气候变化响应和适应研究当中。利用树木年轮材料对黄土高原人工刺槐林生长与气候条件的生态适应性研究也逐渐引起人们的关注,胡万银[10]利用山西中部地区刺槐和旱柳的树木年轮数据,来统计物种生态反应数值对其生态环境适应性进行了研究,表明生态反应值与树木生长情况呈负相关,可以作为树木对生态环境适应能力的评价指标;而李颖[11]则利用山西比较常见的五个外引树种, 法桐、新疆杨、樟子松、刺槐、水杉的树木年轮数据构建树木年轮生态敏感度,对这些物种树木生长情况对环境因子的适应性情况进行了评价。然而基于树木年轮材料,从黄土高原区域尺度上对人工刺槐生长对气候条件的适应特征研究还未引起人们的关注。

本研究在黄土高原沿降雨梯度,在区域尺度获取多点人工刺槐林树轮数据,并探讨了人工刺槐林树木生长对干旱胁迫响应特征及其空间表现规律,对科学认识当前气候条件下黄土高原人工刺槐林生态适应性有重要意义,也将为在气候变化背景下该地区人工刺槐林生态恢复建设提供科学依据。

1 材料与方法

1.1 研究区概况

黄土高原位于黄河中游(34°—40° N, 101°—113° E),位于我国中北部地区,总面积约64万km2,该区域自北而南跨越我国温带和暖温带,从东南向西北包括半湿润、半干旱和干旱地带,年均温度为 年均降雨量为250—600 mm、海拔范围为800—3000 m(图1)。黄土高原是世界上面积最大的黄土堆积区,具有连续的第四纪黄土堆积,区内高塬沟壑区和丘陵沟壑区黄土厚度大多在100—300 m[12—13]。 黄土高原大多数区域存在严重的土壤侵蚀问题,是世界上水土流失最为严重的区域之一[14]。

图1 黄土高原人工刺槐林树木年轮采样点及气象站空间分布图Fig.1 Location map of tree-ring sample sites and meteorological stations for the black locust plantations in the Loess Plateau of ChinaPDSI格点,毗邻年轮采样点的PDSI干旱指数网格点

刺槐林是黄土高原半干旱半湿润地区植被恢复和重建的主要人工林类型之一,作为一种典型的中生树种,既喜湿润肥沃的土壤又耐干旱瘠薄,有较强的适应性和抗逆性;在20世纪70年代末到80年代初,陕北黄土高原从南到北进行了大面积的栽植,对改善这一地区的生态环境、防治水土流失,乃至调节黄河流域的水文状况发挥了重要作用[1, 6]。

1.2 年轮样品采集及处理

本文所用人工刺槐林树轮资料是2015年6月至2016年9月在黄土高原南北样带上获取的,一共采集了6个年轮样点,从南到北依次为永寿、宜君、富县、延安、绥德、神木,采样点年均降雨量变化范围为420—600 mm,海拔高度变化范围为 1180—1350 m(图1,表1)。选择人工刺槐林分布典型、受人为活动干扰较少的地段采集年轮样品;为获得较长的树木年轮样本,在树轮样品采集点选择年龄较老的树木个体获取样本。样芯采集通常在胸高部位,在垂直于山坡方向上用生长锥钻取树轮样本。取到的样芯放置在塑料样管内,并在样管上用油性笔标注代码。每个采样点保障采集30—40刺槐植株,每个植株采集一根树木年轮样芯。

表1 黄土高原人工刺槐林树木年轮样点统计表Table 1 The basic information of tree-ring chronologies for black locust plantations in the Loess Plateau

样品预处理基本按照Stokes和Smiley[15]的方法进行。首先将样品放置于平坦处晾干,然后将其粘在特制的木槽中;将粘牢的样芯依次用颗粒由粗到细(280—800目)的砂纸进行打磨。

1.3 树轮年表构建

样本经预处理后, 用骨架示意图进行交叉定年[8],然后使用LINTAB5.0年轮分析仪对树木年轮宽度逐年进行测量,该测量系统精度可达0.01 mm。 利用COFECHA程序[16]对定年和测量结果进行了检验,保证定年与测量的准确性。对于那些效果较差、与主序列之间的相关系数达不到99%置信区的样芯予以剔除,最终206根刺槐树芯用于年表的构建。

经过交叉定年的年轮序列,用ARSTAN程序[17]进行去趋势和标准化,这一过程消除了树木生长中与年龄增长相关联的生长趋势及部分树木之间的非一致性扰动,排除了其中的非气候信号。首先对每个年轮宽度序列进行负指数曲线拟合,以去除与树木年龄有关的生长趋势;如果不成功,则用任意斜率的线性回归进行直线拟合,并以取树轮宽序列和与其对应的生长趋势拟合曲线的比值,对所得到的去趋势年轮序列进行双权韧性平均处理,合成了平均值为l、无量纲的树轮宽度标准年表。

1.4 气象资料

气象数据来自与刺槐年轮采样点距离较近的气象站,一共在黄土高原选取了5个气象站台站,从南到北依次有长武站(对应永寿采样点)、洛川站(对应宜君和富县采样点)、延安站(对应延安采样点)、绥德站(对应绥德采样点)、榆林站(对应神木采样点)。气 候 资 料气象数据由国家气象局(http://cdc.cma.gov.cn/)提供,时间区间为1951—2020年,气候要素包括月平均气温、月平均最高气温、 月平均最低气温、月降雨量。为了能够较好反映出采样点气候干旱程度,本研究选取了与采样点临近的CRU网格化PDSI干旱指数数据,共提取了6个临近样点的PDSI指数格点数据。CRU 干旱指数格点数据从CRU TS v. 3.22 全球气候数据库(http://www.cru.uea.ac.uk/cru/data/hrg/)[18]下载,时间跨度为1951—2020年,精度为0.5度(约为25 km)。

1.5 数据分析

采用相关分析统计树轮生长对于气候要素的响应,考虑到前一年气候状况对当年树木生长的影响[8],所用气候资料的时间跨度为前一年6月到当年9月,气候要素包括月平均气温、 月平均最高气温、 月平均最低气温、月降雨量、PDSI干旱指数,树轮年表与气候要素的响应关系均用DENDROCLIM2002 程序实现[19],相关统计显著性水平为P<0.05。

2 研究结果

2.1 研究区气候变化特征

从气象站过去70 年(1951—2020)气候数据来看(图2),黄土高原年均温度呈线性上升趋势,每10 年上升0.16—0.32℃,以延安采样点气温上升幅度最大;年总降雨量则大致有不断下降的趋势,每10 年下降16.1—47.5 mm,以延安采样点的降雨下降最为明显,另外神木采样点的年总降雨量略有上升趋势;PDSI干旱指数呈线性下降趋势,以分布在南部采样点(永寿、宜君、富县、延安)的干旱指数下降幅度较大(R=-0.34—-0.514,P<0.01),而分布在北部采样点(绥德、神木)的干旱指数下降幅度较小(R=-0.17—-0.268,P>0.05)。

图2 黄土高原人工刺槐林年轮采样点气候数据波动趋势Fig.2 Climatic changing trend of the meteorological station and CRU grid near tree-ring site of black locust plantations in the Loess Plateau of ChinaSM,神木;SD,绥德;YA,延安;FX,富县;YJ,宜君;YS,永寿;PDSI:PDSI干旱指数;R,相关系数;Slope,气候要素上升速率

从总体上看,黄土高原近代大致向气候暖干化方向发展,具体表现特征是年均温度迅速上升,而年总降雨量趋于下降,PDSI干旱指数也成明显下降趋势;另外黄土高原暖干化气候程度以处于中部延安地区最为明显,其次是延安以南地区,而延安以北地区暖干化气候特征较不明显。

2.2 树轮年表统计特征

本研究构建的刺槐树轮宽度标准年表(简称年表)时间长度为33—50年,以宜君刺槐采样点年表年龄最短,长度为33年(1985—2017年),以永寿刺槐年表年龄最长,长度为50年(1968—2017年)(图3)。从刺槐年表指数数十年来波动趋势来看,延安以北样点(绥德、神木)的年表指数有不断下降趋势,树木生长衰退现象比较明显;延安以南样点(延安、富县、宜君、永寿)的年表指数有不断上升趋势,树木生长无生长衰退现象。

图3 黄土高原人工刺槐林树轮标准年表Fig.3 Climatic changing trend of the meteorological station and CRU grid near tree-ring site of black locust plantations in the Loess Plateau of China

年表的统计量特征可以反映树木生长的一些基本特征以及树轮年表所含不同频率信息量的多少(表2)。结果表明,本研究采集的刺槐年表间具有较好的相关性(样芯间相关系数Rbar=0.257—0.557,P<0.01),表明不同样芯之间的轮宽变化有很好的一致性。各年表的一阶自相关系数为0.22—0.558,说明树木前一年的生长状况对当年生长的影响较小[8]。 年表的平均敏感度(MS)数值为0.363—0.446,这说明了树轮资料具有较强的高频波动信号。年表的信噪比(SNR)为3.411—40.827,达到了较高的水平;年表具有较高的样本群体表达信号(EPS),数值为0.773—0.976,这说明了本调查采集样本量所含的信号是能代表总体特征的[20]。以上年表各种统计特征表明了样本之间的树轮宽度变化有很好的一致性,记录了可靠区域气候信息,所以本次调查所用的树轮资料适于进行年轮气候学分析。

表2 黄土高原人工刺槐林树木年轮年表统计量Table 2 Main characteristic parameters of tree-ring chronologies developed from the black locust plantations of the Loess Plateau

2.3 树木生长对气候的响应特征

年表与月际气候数据分析结果表明(图4),刺槐年表与温度相关性以延安采样点最为明显,其中与去年夏季7月(R=-0.387,P<0.01)和秋季9月(R=-0.401,P<0.01),以及当年春季5—6月(R=-0.412—-0.498,P<0.01)和夏季7—8月(R=-0.392—-0.425,P<0.01)温度负相关关系达到了显著性水平;其他样点刺槐年表与温度的关系大致也呈相关关系,具体表现为绥德年表与春季3月温度(R=-0.42,P<0.01)和永寿年表与夏季6月(R=-0.44,P<0.01)温度的相关性。刺槐年表与降雨相关性也以延安采样点最为明显,其中与去年夏季7月(R=0.515,P<0.01)和秋季9月(R=0.489,P<0.01),以及当年冬季2月(R=0.512,P<0.01)降雨正相关关系达到了显著性水平;其他刺槐样点年表与降雨关系较弱,大多未达到显著性水平。刺槐年表与PDSI干旱指数相关性同样以延安采样点最为明显,其中与去年秋季9月到当年秋季9月的正相关水平达到了显著性水平(R=0.426—0.632,P<0.01),其他刺槐样点年表与PDSI干旱指数关系较弱,大多未达到显著性水平,仅永寿年表与当年春季3—5月PDSI干旱指数具有较强的相关关系(R=0.389—0.497,P<0.01)。年表与季节平均(三个月均值)气候要素相关分析结果表明(图5),年表与季节平均温度和PDSI干旱指数的相关结果基本与月际数据相一致,与季节平均降雨相关结果要略弱于月际数据,但是基本规律是相似的,进一步验证了年表与气候数据相关结果的可靠性。

图4 黄土高原刺槐人工林年表与月际气候要素的响应分析Fig.4 Correlation analysis of monthly climate data with tree-ring chronology of black locust plantations in the Loess Plateau of China TMP: 温度;PRE: 降雨;PDSI: PDSI干旱指数;*,达到95%显著水平;-代表上一年

图5 黄土高原刺槐人工林年表与季节平均(3个月)气候要素的响应分析Fig.5 Correlation analysis of seasonal (averaged from three months) climate data with tree-ring chronology of black locust plantations in the Loess Plateau of China

2.4 树木生长对气候响应的多元线性模型

根据刺槐年表与气候要素的相关分析结果,挑选出与刺槐年表相关显著的气候要素(当年5至8月份温度、当年2月份降雨、去年9月份降雨和年平均PDSI干旱指数)与刺槐年表进行多元线性回归和逐步回归分析,以便从整体水平揭示刺槐年际生长变化方差在多大程度上能被气候要素所解释。结果表明(表3),处于黄土高原中部延安刺槐年表与气候要素回归模型结果最为显著,气候要素总体上可解释刺槐年表年际波动的64.5%的方差,回归方差中各气候要素与年表的相关性均通过显著性检验;其次是绥德、富县和永寿刺槐年表,气候要素总体上可解释刺槐年表年际波动的20%左右的方差,而宜君和神木刺槐年表与气候要素的回归模型结果较差,总体解释方差量也较低。

表3 黄土高原南北样带人工刺槐林树轮年轮年表与气候要素的多元线性回归方程Table 3 Most supported ordinary least-squares regression models of climate on tree-ring chronologies of the black locust (Robinia pseudoacacia) plantations along the south-north transect in the Loess Plateau of China

3 讨论

基于黄土高原南北样带尺度的人工刺槐林树木年轮材料,本研究发现延安以北的刺槐样点(绥德、神木)在1990年代以来出现明显生长衰退趋势,延安以南的刺槐样点(延安、富县、宜君、永寿)在1990年代以来有明显生长趋势,无生长衰退现象。基于黄土高原刺槐林生长耗水和土壤水分动态关系研究表明,延安以北地区刺槐林林地土壤出现严重退化,深层土壤含水量产生严重亏缺,在1.6—2.0 m以下出现持续性干层,部分土壤退化严重地点3.0—5.0 m平均土壤含水量已接近土壤萎蔫含水量,严重影响到人工刺槐林的持续健康生长,树木生长出现明显衰退,最终导致大面积生态低效的刺槐“小老头树”林出现[21—22]。王力等[1]通过对陕北黄土高原南北 6 个人工刺槐林样地进行标准地调查,分析了人工刺槐林树高、胸径、材积等指标的生长状况与土壤含水量亏缺的关系,研究发现黄土高原北部刺槐林的平均树高(4.5 m)、胸径(5.8 cm)、材积(10.79 m3/hm2)水平要明显低于南部刺槐林(平均树高、胸径、材积分别为11.8 m、11.5 cm、1.62 m3/hm2),这表明黄土高原北部刺槐林树木生长受到明显的干旱因素限制,出现了明显的生长衰退现象,这与本研究的树木年轮学结果具有较好的对应关系。

从刺槐树轮年表与气候要素响应特征来看,处于黄土高原中部的延安地区刺槐年轮生长对气候响应敏感性最高,其年表与温度呈负相关关系,而与降雨和PDSI干旱指数呈正相关关系,体现出对干旱胁迫气候具有较强的响应特征。延安地区处于黄土高原森林植被向草原植被的过渡地带[1,7],年降雨量大约在500 mm左右,应该是处于黄土高原人工刺槐林适宜生长区域的北界,一般认为森林生长在其空间分布的边缘地带,即森林向灌丛或草地的过渡地带,对外界环境条件变化最为敏感,譬如森林分布海拔上限的树木生长对温度响应最为敏感,而森林分布海拔下限的树木生长对降雨响应最为敏感[23—24]。基于延安南泥湾地区天然辽东栎的树木年轮材料,买尔当·克依木等[25]发现该地区辽东栎径向生长对气候响应敏感性具有时间不稳定性特征,表现在与温度由正相关转为负相关,与降雨正相关关系则逐渐加强,这表明延安地区天然辽东栎林生长主要受到水分胁迫条件的影响,而这是与本研究相互吻合的。

延安以北地区属于典型草原植被带,年降雨量一般在450 mm左右,该地区的水分条件已不足以支撑人工刺槐林的正常发育,本研究树轮年表数据显示该地区刺槐树木生长速率较慢,树木生长对气候响应敏感性也较低,体现出刺槐生长对干旱气候条件体现出的一定的适应性特征。已有许多研究表明黄土高原北部干旱区刺槐林为了维持在干旱条件下的新陈代谢,并减少水分散失,在一系列生态生理特征上体现出对干旱气候的适应性,譬如比叶面积较小、叶片较厚、叶片栅栏组织发达、植物比根长和可溶性糖分含量较高等[26—28];刺槐对干旱气候的适应策略需要分配更多的物质给吸收和保持水分的组织和器官,从而对刺槐生长的物质投入较少,使得刺槐生长速率较慢且对干旱胁迫条件不具有敏感性特征[29—30]。 延安以南地区属于典型森林植被带,年降雨量一般在550 mm左右,该地区气候条件相对优越,土壤水分条件足以支撑人工刺槐林生长需求,本研究表明该地区刺槐树轮年表生长速率较快,但年表对气候要素响应敏感性较低,也无明显干旱胁迫信号,这主要是因为延安以南地区是人工刺槐林生长适宜地区,刺槐生长还未受到水分胁迫的限制性影响的原因[31—32]。

4 结论

本研究依照树木年轮气候学传统技术和方法,在黄土高原沿南北样带尺度(降雨梯度为180 mm量,从420 mm到600 mm)构建了6个地点(永寿、宜君、富县、延安、绥德、神木)人工刺槐林树轮年表宽度序列,并对这些刺槐年表的生长趋势和气候响应敏感性特征进行了分析。结果表明延安以北绥德和神木采样点(降雨量大约450 mm)属于典型草原植被带,气候条件较为干旱恶劣,刺槐林年轮指数近期趋于下降,树木有较为明显的生长衰退现象,刺槐年表中气候信号较弱,体现出树木生长对干旱胁迫气候一定的适应性特征;延安采样点(降雨量大约500 mm)位于黄土高原中部,属于森林植被向草原植被过渡地带,气候条件适中,刺槐树木生长速率有上升趋势,树轮年表对气候响应敏感性较高,树木生长对干旱胁迫响应最为敏感;延安以南富县、宜君和永寿采样点(降雨量大约550 mm)属于典型森林植被带,气候条件较为优越,刺槐树木生长速率有上升趋势,树轮年表对气候响应敏感性较低,这是由于干旱胁迫还未对刺槐生长产生限制影响的原因。

就黄土高原区域人工刺槐林经营和维护方面而言,黄土高原南部(延安以南地区)人工刺槐虽未发生生长衰退问题,刺槐树轮年表对气候响应特征也未体现出明显的干旱胁迫特征,但是由于刺槐种植密度过高,需要对人工刺槐林进行适当疏伐和群落结构改良措施,以免刺槐林在其生命周期中后期出现大面积生长衰退问题;黄土高原北部(延安以北地区)人工刺槐林已出现明显生生长衰退问题,对气候响应敏感性也较差,刺槐林生长处于明显亚健康状态,且具有明显小老头树特征,很难发挥其生态效益,在延安以北地区应逐步以耗水量较小的灌丛和草地恢复措施替代人工刺槐林。