森林健康是森林生态系统健康与恢复的简称, 是指森林在发挥必要的生态服务功能的同时, 保持其自身良性存在和更新的状态^[1]。近年来由于全球气候变化所带来的高温和干旱等极端气候事件的频发, 使得世界上大多数国家的森林生态系统都存在不同程度的退化^[2-4], 主要表现为树木在生长发育过程中的生理机能下降、生产力降低、生长发育滞缓或死亡等现象。在我国, 森林覆盖率已由建国初期的12.5%上升至2013年的21.63%^[5]。然而与快速恢复性增长的森林覆盖率呈鲜明反差的是, 许多现存的森林处于破碎化或次生演替状态。中国森林退化问题突出, 森林健康状况堪忧, 影响了森林发挥必要的生态服务功能^[1], 使得森林健康的研究十分迫切。

目前国内外学者对森林健康的研究提出了许多评价方法^[6-8], 但是由于对森林群落缺乏长期观测, 人们难以得知其健康历史。树木年轮记录了树木个体的生长历史, 具有精确定年、分辨率高、连续性强和复本量好等特点^[9], 在反映过去森林衰退方面有着突出的贡献^[10-12]。研究表明, 树木在生长过程中会受到不同范围生物和非生物多元化组合的影响, 可能出现树冠枯死、低生长或死亡等不健康生长的现象^{[4, 13-14]}。由于限制因子发生的顺序、时间和强度的不同, 其影响可能是持续或者短暂的, 森林健康在空间上的表现较为复杂^{[2, 15-17]}。因此, 过去的森林不健康事件在发生程度和持续时间上有何特征是一个亟待解决的生态学问题。对森林健康历史时空特征以及其驱动因子的深入研究很大程度上丰富了历史时期森林生态学特征的研究。

忻州地处山西高原中北部, 位于吕梁山生态脆弱区。针对上述问题, 我们该选择地区保存较好的4个森林为研究对象, 试图利用树木年轮学评价该地区森林健康状况, 旨在研究以下两个具体问题：1)忻州地区森林健康状况在过去100年中发生了怎样的变化？2)研究时间段内多次森林不健康事件在空间上有无同步性？研究结果所得森林健康时空信息也可为林业部门制定合理的经营管理方案提供数据依据。

1 研究材料与方法 1.1 研究区域概况

本研究区位于山西省忻州市地区(38°08′—39°40′N, 110°53′—113°58′E), 西临黄河, 东止太行山。该区地形西高东低, 逐步倾斜, 北、西、南三面环山, 东部开阔平坦, 为忻定盆地的主体部分。境内流径滹沱河、云中河及牧马河等重要河流。由于受温带大陆性季风气候的影响, 该地区夏季多东南风, 冬季多西北风, 春温高于秋温。年平均气温为8.5℃左右, 年平均降水量为405 mm, 降水集中于7—9月份^[18]。

1.2 树轮样品采集与气象资料

忻州地区自然植被资源丰富, 全区森林面积4.07×105 hm², 森林覆盖率为16.4%。本研究选取忻州地区植被覆盖率较高的4个县作为采样点。其中宁武梅洞(38°71′N, 111°96′E；海拔1772 m；青扦Picea wilsonii)和大石洞林场(38°92′N, 112°08′E；海拔1752 m；青扦)属于吕梁山脉, 坡度5°—10°。繁峙沙河宽滩林场(39°06′N, 113°47′E；海拔2184 m；臭冷衫Abies nephrolepis)属于五台山区域, 坡度为0°—5°。岢岚县中寨林场(38°58′N, 111°43′E；海拔1745 m；油松Pinus tabuliformis)属于黑茶山区域, 坡度为15°—20°(图 1)。4个采样点区域天然次生林植被丰富, 主要有白扦林、青扦林、华北落叶松林、白桦林、山杨林、辽东栎林、油松林等。野外样品采集于2012年10月, 我们选择区域内年龄较老的树木, 使用生长锥在胸高处(1.3 m)采集树芯样本。取样方向为沿山坡的等高线方向, 每株立木采集一根树芯样本, 共采集了86棵树的树芯样本。

本研究选用距离采样点较近位置的原平气象站(38°44′N, 112°43′E；海拔828.2 m, 数据自1954年起测定)和五寨气象站(38°55′N, 111°49′E；海拔1401 m, 数据自1957年起测定)的气候数据(图 1)。气候因子包括月总降水量(Monthly total precipitation, Pm)以及月平均气温(Monthly mean temperature, Tm)。由于两地同属暖温带季风气候区域, 气象站数据差异性较小, 因此将两个气象站的气候数据平均值作为该区域的气候数据^[19], 分别分析不同样点与区域气候的响应关系。干旱指数Palmer drought severity index(PDSI)数据来自荷兰皇家气象研究所的数据共享网(http://climexp.knmi.nl/), 区域范围为38.59°—39.06°N, 111.43°—113.47°E。

1.3 树轮数据的分析方法

样本在自然状态下干燥后, 按照基本程序进行固定、磨光处理和交叉定年^[20]。在测量精度为0.001 mm的LINTAB树轮宽度测量仪(RINNTECH, Heidelberg, German)上测量每一年轮的宽度, 并利用国际树木年轮库的软件程序COFECHA^[21]进行定年质量控制。更正因缺轮、伪轮和测量误差等带来的错误, 最终使树芯样本上的每一个树轮都对应其实际生长年份。应用ARSTAN^[22]程序, 用负指数函数或直线拟合每个测量序列的生长趋势, 用加权平均法将去趋势序列合并成树轮宽度指数序列。研究计算了宁武县两个采样点年表间的滑动相关系数矩阵, 将梅洞和大石洞的样本合并。同时设定了分析健康历史的起始年, 采用子样本信号强度SSS(Sub-sample Signal Strength)＞85%的样本量作为分析的起始点, 高于该样本量的健康序列认为是可靠的。

森林中树木个体生长的健康状况由于树木自身或局部生长小环境的原因会产生差别, 在不利环境条件下, 一些树木具有较强的恢复力而能够及时响应环境的变化, 而一些树木受到外界影响后不能及时恢复到健康的状况, 从而导致了树木的健康与不健康的分异。不健康树木较健康树木年轮生长较慢, 而且不健康树木的低生长是一个持续的过程。结合树木年轮数据, 定义树木的不健康时段表现为树轮指数连续5 a低于正常值。基于此定义, 从下面3个步骤挑选出不健康时段的树木：1)挑出树木至少持续5 a低生长片段(树轮指数小于0.9)；2)挑出同时期树轮指数平均值低于0.7, 且期间指数最小值低于0.5的树为非健康树；3)统计各点每年中出现非健康树的样本量所占百分比。将挑出的非健康树年份中的树轮值平均作为非健康树年表, 同时将非健康树轮值剔除平均剩余值作为健康树年表。通过对比各采样点结果, 找出该区域时间和空间上非健康历史的特征。

为了研究该地区森林生长在对当地气候变化的响应特征, 研究利用DendroClim2002程序^[23]分别分析各点的树轮宽度标准年表与逐月气候因子的相关关系, 逐月气候因子包括前1年10月到当年9月的气候指标(月总降水量和月平均气温、PDSI)。

2 研究结果 2.1 树轮样本及交叉定年结果

经过对树芯样本测量和交叉定年后, 得到了共86个年龄准确可用于分析的树轮数据。各年表的统计特征见表 1, 其中该地区年龄最大的树为223 a, 年龄最小的树为25 a。各采样点序列间平均相关系数都达到了0.5以上, 表明该区域内各点树木个体间的轮宽变化均较为一致, 对外界环境变化有着一致的响应。一阶自相关系数也都达到了0.6以上, 表明上一年气候状况对当年树轮生长有着持续性影响。样本总体代表性都高于0.85, 各样点含着大量的共同信号。

2.2 树轮和气候因子的响应关系

各样点标准年表与气候响应分析结果如图 2所示。从图中可以看出, 研究区域树轮宽度指数与月平均气温和月总降水量的相关结果有较好的一致性。各样点树轮主要响应于生长季中期的温度与生长季末期的降水, 其中, 以与当年5—6月月平均温度的负相关关系最为显着(r=-0.51—-0.26, P＜0.05)；在树木径向生长与降水因子的关系中, 主要与7月份的降水显着的正相关(r=0.29—0.43, P＜0.05)。从树轮指数与PDSI相关的分布图可以看出, 4个采样点树轮与该区域PDSI呈显着的正相关关系(r=0.22—0.55, P＜0.05)。

2.3 健康树与非健康树的径向生长历史

各采样点健康和非健康树轮指数序列如图 3所示。结果表明, 岢岚地区健康树1900—1930年间, 生长量呈现明显的下降趋势, 此时非健康树树轮指数也较低, 其中1910—1930年期间均值仅为0.43；在1960年之后健康树和非健康树也表现出低生长的趋势, 一直持续到2005年。繁峙地区非健康树在1920—1940年间树轮指数均值为0.55；健康树和非健康树在1990—2010年期间表现出明显的持久下降趋势。宁武地区非健康树在1910—1940年期间, 树轮指数均值为0.51, 且在1926年极低；在1970—1980年和2000年以后树轮指数较低, 均值分别为0.58和0.54存在稳定下降趋势。

研究统计了各采样点非健康树所占百分比情况(图 4), 发现样本量百分比峰值年较为一致。在过去1个世纪的公共时期内, 岢岚地区最为显着的不健康事件发生在1913—1936年期间, 持续了25年, 通过定义判断的非健康树木量达总体的40%以上。而宁武地区的非健康树样本量从1915年就开始持续增多, 1923年达到峰值, 期间非健康树所占比例达到30%以上。超过30%样本量的时期持续了18 a。繁峙地区为1925—1942年, 持续了18 a。70年代的非健康事件仅在岢岚和宁武地区有所体现, 都是持续了5 a, 宁武从1973年开始发生, 岢岚从1983年开始。20世纪末, 该区域的不健康事件在3个研究区域都有体现, 都是从1999年开始。其中沙河地区持续时间最长为14 a, 一直持续到采样时间；岢岚地区11 a；宁武地区持续时间较短, 为6 a。

3 讨论 3.1 树木径向生长与气候的关系

忻州市地处中纬度地带内陆, 受大陆性季风气候以及北部地区的风蚀, 使得该地区的年蒸发量远超过于降水量。在干旱及半干旱地区, 降水量是树木生长的主要限制因子之一^[24]。本文中树轮宽度主要与生长季中期5—6月的月平均温度以及7月总降水量呈显着的相关关系。树轮宽度与降水和温度的这种相关关系与戴君虎等^[25]在山西忻州五台山地区、李颖俊等^[26]在芦芽山地区、蔡秋芳等^[27]在吕梁山地区的研究结果相似。忻州地区5—6月份降水量不足, 较高的温度必然会加速树木的蒸腾量, 使树木体内的水分散失过快, 同时也会加速土壤水分的蒸发量, 形成水分胁迫限制了树木的生理代谢活动, 从而降低了树木的生长量。7月份降水充足, 同时较高的气温和地温也使该地区冻土融化, 有效的提供树木径向细胞分裂和伸长所需的必要水分, 促进了树木的生长。树轮与PDSI的相关系数结果要高于降水或温度(图 4), 也表明该地区的树木主要受到水分作用的限制, 更多地体现在土壤水分的有效供应方面。Cai等^[28]在山西太行山中部地区和Sun等^[29]在汾河流域的研究也得到了类似的结果。

3.2 森林健康的时空特征

由于微环境、限制因子、树木个体耐受力、群落竞争等状况的影响, 树木生长个体间存在着差异^[30]。本文利用生长指数的差异性挑选出特定时期健康树和非健康树, 得到了该区域在时间和空间尺度上的森林健康信息, 并且非健康树所占百分比的统计进一步补充了不健康事件发生的程度。各采样点非健康树样本情况(图 4)显示该区域出现高于30%阈值的不健康事件在空间上表现出明显的同步性, 主要集中在1910—1940、1970—1987和1990—2012年期间。各地区不同时期的不健康事件持续时间和强度不尽相同, 如岢岚地区在20世纪30年代不健康事件持续时间最长强度也最大。20世纪70年代时期, 岢岚和宁武非健康树比例都较高, 而该时期繁峙地区非健康仅占少数比例。20世纪末, 繁峙地区不健康事件持续时间最长。梁尔源等^[31]对中国中西部地区树轮研究发现, 极端干旱事件导致20世纪20年代和20世纪30年代初期大范围的树木生长下降现象。Liu等^[10]探究了亚洲中部半干旱地区森林衰退现象, 也发现20世纪70年代末80年代初和21世纪初期树轮表现出显着的生长下降。我们认为, 大尺度的气候事件导致森林不健康事件发生的空间格局具有很强的一致性。不同林分特征的森林对气候变化的抗干扰能力和恢复能力有显着的差异^[32], 并且不同树种对气候变化的敏感度不同^[30], 这可能是影响不健康事件发生强度和持续时间不一致的重要因素。

3.3 干旱事件对森林健康的指示

Liu等^[10]认为气候变暖和干旱加速了亚洲中部半干旱地区森林树木生长的下降。Liang等^[11]探究了青藏高原东北部半干旱地区云杉林在全球变暖背景下的森林衰退现象, 也认为干旱事件导致了森林衰退, 快速的气候变暖对森林树木的生长有着非常重要的影响。本研究发现该区域三次较为严重的非健康事件爆发的高峰主要集中在1930s、1970s和2000s左右。根据研究区的气象记录, 忻州地区5—9月平均气温与同一时期降水量显着负相关, 即研究区气候具有高温少雨、低温多雨的组合模式。在中国北方地区, 高温少雨总是和干旱联系在一起。20世纪20—30年代前后由于东亚夏季风减弱而导致中国北方降水量明显减少, 而这次著名干旱事件在许多树轮研究中也得到印证^{[31, 33]}。李强等^[34]对山西宁武地区降水重建资料显示, 1966—1972年和1988—2003年为降水较少的时期(低于多年平均降水量413 mm), 其中1972年干旱尤为严重。蔡秋芳等^[27]对山西吕梁山气温变化的研究中发现, 1994—2002年的温度是过去近170年中温度最高的时段, 这一持续高温事件对此次干旱事件的爆发无疑贡献巨大。Williams等^[35]认为干旱会导致气孔关闭, 光合作用降低, 从而限制了形成层活动, 使得树木在生长季受到严重的水分胁迫。本研究发现的3次较大的不健康事件与历史时期严重干旱事件时间相吻合, 表明该地区降水减少或高温所带来的极端干旱是影响该区域森林健康生长的主要原因。

20世纪中叶以来的气候变暖是不容置疑的^[36]。忻州地区地形崎岖, 多为山地丘陵, 水分条件较差, 干旱事件发生的频率和干旱程度在未来气候变化背景下都可能会增加。高温对育苗造成的严重危害和干旱对造林成活率和保存率的严重影响等都是该地区林业发展的不利因素。要针对这些不利的限制性气候因子, 选择合理的造林季节, 避开高温和干旱时期因地适宜地选择耐旱树种进行栽培和经营, 实现林业的可持续发展。

4 结论

本文以忻州地区4个采样点森林为研究对象, 利用树木年轮生态学方法重建了自1900年以来忻州森林的健康历史, 分析了森林健康状况的时空变化特征, 并对产生不健康事件的原因进行深入探讨。研究结果表明：过去100年中该地区森林不健康事件主要集中发生在1910—1940、1970—1987和1990—2012年期间, 4个采样点在20世纪30年代和世纪末不健康事件在空间上较为一致, 20世纪70年代的不健康事件在宁武和岢岚地区表现较强, 而在繁峙地区仅有少量树木表现出非健康生长。在经历了3次较严重的不健康事件的影响, 树木能够度过环境不利时期依然存活, 表明该地区的森林具有较强的抵抗力或生态弹性。升温或降水减少所引起的极端干旱是造成该区域森林非健康生长的主要原因, 也是不健事件空间一致性的主要驱动力。研究结果有助于进一步认识区域森林健康在全球气候背景下时间和空间上的变化特征以及其驱动因子, 对科学地制定森林保护措施具有积极的现实意义。在以后的研究中仍需加强空间同步信号的鉴别能力, 在更广泛的区域和气候背景下开展森林健康的研究, 来提高我们对森林衰退的预防和治理能力。

致谢: 中国科学院植物研究所树木年轮实验室提供树轮样本, 张齐兵研究员在论文写作过程中给予建议, 邱红岩老师帮助实验工作, 特此致谢。