暖温带北缘3种栎树径向生长对气候要素的响应差异及预测

王奇 孙婧依 刘建锋 常二梅 冯健 张晓放

摘要：[目的]比较同一气候区不同栎类树种生长对气候要素的响应差异，特别是在分布边缘不同树种的种内种间差别。[方法]利用树木年轮学方法，分析了我国暖温带北缘3种优势落叶栎树（蒙古栎、栓皮栎和麻栎）径向生长对气候要素的响应差异；同时利用混合线性效应模型构建了3个树种径向生长与气候的关系，并基于未来气候情景（RCP2.6和RCP8.5）预测了其未来生长趋势。[结果]同一地区栎类树种年表具有较高的相似性，而不同地区间差异较大；在径向生长对气候响应方面，同一地区的不同栎类树种和不同地区同一栎类树种也存在一定差别。其中北京地区的3种栎类树种的径向生长与当年6月气温呈负相关关系，但在辽宁仙人洞地区却与气温无显著相关关系；北京地区的麻栎、河北地区的蒙古栎、辽宁大黑山地区的栓皮栎的径向生长均受降水条件限制。模型预测显示，至本世纪末，北京地区和辽宁大黑山地区栎类树种的径向生长呈现下降趋势；而辽宁仙人洞地区栎类树种径向生长呈增加趋势。[结论]我国暖温带北缘不同地区的栎类树种对气候要素的响应存在显著差别．主要为取样地区的气候所支配，下一步尚需结合林分及立地因子，来厘清各类环境要素的相对贡献度。

关键词：气候变化；树本年轮；栎类树种；径向生长

中图分类号：S718.43 文献标识码：A 文章编号：1001-1498（2023）03-0062-09

气候是影响全球植物分布的关键驱动因素。随着全球气候变化的不断加剧，植物的生长及其分布范围，特别是在其分布边缘将受到更大的影响。有研究发现，相对于植物分布范围的变化，植物的生长对气候变化更加敏感。因此，探索不同植物生长动态对气候变化的敏感性，可为预测未来植物分布动态提供重要参考。森林是陆地生态系统碳汇的主体，研究某一区域森林树种的生长动态及其对气候变化的响应差异，有助于理解气候变化背景下森林组成、结构及碳汇功能的变化动态，为适应性经营管理措施的制定提供科学依据。

已有研究表明，树种分布北缘的群体生长主要受低温影响，但在其分布南缘则受到高温导致的干旱缺水的限制。这些研究大多基于阿尔卑斯山脉以北的高纬度地区或以南的地中海地区。然而，在暖温带气候区（如暖温带湿润—半干旱地区），树种的生长是否也存在上述现象尚缺乏研究。此外，不同分布范围的树种的生长可能会对气候变化做出不同的响应。如Huang研究加拿大东部的4种北方树种时发现，毛白杨（Populus tomentosa Carr.）的径向生长主要受水分条件的影响，而白桦（Betula platyphylla Suk）的生长受1月的温度和生长季降水影响，云杉（Picea asperata Mast.）和班克松（Pinus banksiana Lamb.）的径向生长则与当年冬季和春季或整个生长季的温度呈正相关关系。然而，Tardif等发现，白栎（Quercus alba L）和红栎（Q.rubraL）的生长在其北部范围内对气候的响应并没有显著差异。因此，不同树种在分布范围边缘对气候的生长反应是复杂的，目前尚缺乏对同域分布但处于不同树种各自不同分布边缘的群体对气候变化响应的比较研究。气候暖化背景下森林的响应与适应是当前生态学研究的重要方向，一方面，气候变暖可以加速土壤有机质的分解速率，从而增加植物对土壤矿物质营养的吸收，特别是对最常见的限制性营养素N的吸收；另一方面，温室气体特别是大气CO2浓度上升，加上环境温度的升高，提高了植物光合同化物质的生产和水分利用效率，从而增加树木生物量的积累。有研究发现，气候暖化可能导致北半球温带树种在其分布北缘呈现增长趋势。然而，在我国暖温带树种分布的边缘群体是否会有类似趋势尚缺乏研究。

栎类是我国面积和蓄积最大的树种组，在木材生产、固碳增汇、土壤保持、水源涵养及生物多样性维持等方面发挥着重要的功能。在栎类树种中，拴皮栎（Q.variabills Bl.）、麻栎（Q.acutissima Carruth.）和蒙古栎（Q.mongolica Fisch）是我国分布最广、数量最多的3种落叶阔叶树种，其范围分别为（19°～42°N、97°～140°E）、（19°～42°N、97°-140°E）和（34°-53°N、112°-134°E）。也是我国暖温带森林的主要成分。我国暖温带北缘是栓皮栎和麻栎天然分布的北缘，蒙古栎分布的南缘。有研究表明，该区域气候暖干化的趋势明显。因此，本研究以生长在暖温带北缘的栓皮栎、麻栎、蒙古栎为研究对象，通过树木年轮学方法，分析径向生长与气候因子的关系；同时，考虑未来气候变化情景（RCP2.6；RCP8.5），对该区域栎类树种的未来生长展开预测，以期为本区域制定应对气候变化的森林管理对策提供依据。

1材料与方法

1.1研究区域概况

在我国暖温带北缘（华北地区北部和东北地区南部）选择栎树树种分布的6个典型样点，分别为西部的北京市昌平、密云、怀柔和河北省棋盘山地区；东部的辽宁省大黑山和仙人洞地区（图1）。研究区域属典型暖温带季风气候，夏季炎热多雨，冬季寒冷干燥。但3个地区的水热条件存在明显差异，如东部辽宁地区温暖湿润，降水量丰富，年平均气温为10.2℃，年降水量644.7 mm；其中仙人洞和大黑山年内降水季节性也存在一定差异，如仙人洞降水集中在5-7月，大黑山降水集中在6-8月。西部的北京地区气温较高，雨水相对较少，年平均气温为9.8℃，年降水量490.3 mm；河北棋盘山地区气候干燥寒冷，雨量相对最少；年平均气温仅为2.3℃，年降水量415.1 mm（图2）。

1.2年轮样品采集与处理

于2018年8-9月展开野外取样工作。在每个样点，选择人为干扰较小且栓皮栎、麻栎和蒙古栎为优势树种的林分，建立20m x 20m的临时样地，选择每个样地的冠层优势个体，按照标准树木年代学方法进行采样。对于每棵样树，使用5.15 mm树木生长锥在1.3 m处按十字交叉取2个树芯，并记录胸径（离地面1.3m）以及样地基本信息（纬度、经度、海拔、坡度、坡向等）（表1）。将样品储存在定制塑料吸管中，风干后用白胶安装木芯样品，用200目至1 000目砂纸逐级打磨，直到年轮边界清晰可判读，然后用扫描仪（中晶i800plus）扫描成TIF格式图片（分辨率1 200 dpi），再用WinDENDRO（Regent Instruments，Canada）测量（精度0.001 mm）。所有轮宽序列通过目视检查和COFECHA进行交叉定年验证，最后共获得329个样芯（表1）。

1.3气候数据

过去时期（1930-2018）的气候数据来自Climate Research Unit（CRU TS 4.05）地表逐月气候格点数据集（https：//crudata.uea.ac.uk/cru/data/hrg/，网格距为0.05°×0.05°）。根据样地经纬度进行提取，气候要素包括月均气温、月降水量和月尺度帕默尔干旱指数（PDSI）。未来时期（2018-2100）的气候数据来自亚太地区历史和未来气候数据集ClimateAP v2.30（https：//asiapacifle forestry.ubc.ca/research-approaches/climate-modeling/）。本研究选择了2种未来温室气体排放情景，即CanESM2-RCP2.6和CanESM2-RCP8.5。未来气候要素包括年平均温度（MAT）、年降水量（MAP）、Hargreaves参考蒸发散（EREF）和哈格里夫斯气候水分亏缺（CMD）。

1.4统计分析

由于当年的气象因子会影响树木次年或后几年的径向生长，存在一定的“滞后效应”，本研究用上年6月至当年9月的气象数据与当年的轮宽指数（RWI）进行相关分析。为了确定研究区域栎类树种生长趋势是否趋同或趋异，根据RWI进行了主成分分析（PCA）和聚类分析。未来气候变化下的树木径向生长预测使用线性混合效应模型（LMM）。以上分析均利用R软件（V4.13）完成，使用R中树轮气候分析程序包“treeclim”（V2.0.5.1）中的“dcc”函数来进行树轮年表与气候要素的相关性分析，使用“corrplot”（V 0.92）程序包绘制上述分析的相关关系图；使用“Vegan”（V2.6-4）包来进行主成分分析，使用“factoextra”（v1.0.7）包来进行聚类分析，使用“Ime4”（V1.1-31）包进行混合效应模型的拟合。

2结果与分析

2.1研究区气候状况及年表特征

根据1930年至2018年的气候数据（图3），研究区域各样点的年平均温度增幅为0.015℃。年平均温度由高到低依次为辽宁仙人洞（8.76～11.90℃）、辽宁大黑山（7.41-11.70℃）、北京（8.38～11.55℃）、河北棋盘山（0.60～4.25℃）。年总降水量由高到低依次为辽宁大黑山（466.32～1147.03 mm）、辽宁仙人洞（378.41～1 008.51mm）、北京（243.77-735.82 mm）、河北棋盘山（293.74～579.22 mm）。年降水量在1950年至2018年明显有下降趋势，降幅由高至低依次为辽宁仙人洞、辽宁大黑山、北京、河北棋盘山。

为更好保留年表的低频气候信息，本研究构建了标准化年表，最终建立9个年表（表1，图4）。由表1可知，在3个树种中，北京怀柔地区麻栎年表最长（1931-2018）．北京栓皮栎年表最短（1978-2018）。栓皮栎的样本总体代表性最大，北京、河北地区树木年表相对于辽宁地区有更高的样本总体代表性。

2.2树轮宽度指数的主成分分析及聚类分析

不同地区栎类树种年表PCA显示，第一主成分（PC1）的方差累积贡献率为38.56%；第二主成分贡献率为14.43%。可以看出，不同栎树的RWI在同一地区具有更高的相似性。除河北棋盘山地区蒙古栎生长与其他地区差异较大外，北京地区栓皮栎、麻栎、蒙古栎生长相似，辽宁仙人洞地区麻栎、栓皮栎生长相似，大黑山地区麻栎、栓皮栎和仙人洞地区蒙古栎生长相似（图5A）。聚类分析结果进一步验证了PCA分析的结果，如辽宁大黑山地区栎类树种为一类，辽宁仙人洞地区和北京地区栎类树种聚为一类，而棋盘山地区栎类树种则单独为一类（图5B）。

2.3树轮宽度对气候的响应

各树种对主要气候变量的相关分析如图6所示。在本研究区域的西部（华北地区北部），北京地区栓皮栎和麻栎的RWI与冬季（上年12月）温度呈负相关关系；而北京地区和河北棋盘山地区的蒙古栎的RWI则分别与上年12月和当年2月温度表现出正相关关系；华北北部的3种栎类树种均与夏季（5-7月）温度呈负相关关系。北京地区栓皮栎和麻栎的RWI对逐月降水量的敏感性（正相关关系）要高于本地区的蒙古栎，河北棋盘山蒙古栎的RWI与月降水量也表现出较高的敏感性，且全年均呈显著正相关关系。北京地区的蒙古栎RWI与冬季PDSI（上年11月至当年1月）呈现负相关关系，北京地区的栓皮栎、麻栎及河北棋盘山的蒙古栎RWI与PDSI多表现出正相关关系。

在研究区东部（东北地区南部），辽宁大黑山地区的栓皮栎和麻栎的RWI较之仙人洞地区3个栎类树种对春季（当年2-4月）温度的敏感性更高，均表现为显著的负相关关系，而仙人洞地区的栎类树种则与春季温度无显著相关关系；但上年11-12月温度对仙人洞地区栓皮栎RWI的正相关显著。东北地区南部的所有栎类树种对各月降水量均不敏感，但与干旱指数PDSI多表现为正相关关系。

2.4不同地区栎类树种生长预测

基于CanESM2-RCP2.6和CanESM2RCP8.5两种气候变化情景，通过构建混合效应模型，对不同地区栎类树种生长的模拟预测时发现（图7），在RCP2.6情景下，研究区域的栎类树种生长状况较为稳定；而在RCP8.5预测情景下，仅北京地区栎类树种生长状况较为稳定，辽宁地区栎类树种在2050年左右开始出现较大变化，其中仙人洞地区栎类树种径向生长加快，但大黑山地区的径向生长开始逐渐减弱。

3讨论

识别森林群落中不同树种生长对气候响应的差异，是评价气候变化下森林动态变化的重要前提。本研究研究了我国暖温带北缘栓皮栎、麻栎、蒙古栎径向生长对气候的响应差异，发现生长在同一地区的不同树种存在较高的相似性，但在不同地区差别较大。这与王兆鹏对福建柏（Fokienia hodginstiDunn.）、铁杉（Tsuga chinensis Franch）、冷杉（Abies beshanzuensis Mast.）、马尾松（Pinus massoniana Lamb.）的研究结果相一致。如罗霄山南部4个针叶树种的径向生长与单月和季节性气候因子的相关结果大致相同；气温发生突变后，4个树种树轮宽度指数的变化趋势也相似。聚类分析中，河北棋盘山栎树（蒙古栎）年表与其他区域差异较大，被最先单独分离出来，可能的原因是由于河北棋盘山地区纬度较高，年降水量、年均温明显低于其他地区。

整体上，相对于蒙古栎，栓皮栎、麻栎对温度和降水具有更高的敏感性。本研究结果显示，北京地区栓皮栎和麻栎的RWI与冬季（上年12月）温度呈显著负相关关系；而北京地区和河北棋盘山地区的蒙古栎的RWI则分别与上年12月和当年2月温度表现出显著的正相关关系。可以认为，华北地区冬季增温会导致该地区的栓皮栎和麻栎生长下降，而蒙古栎生长增加；但夏季温度的升高会导致整个区域栎类树种生长的下降。华北北部的3个栎类树种均与夏季（5-7月）温度呈负相关关系，这可能是由于夏季高温会诱导植物叶片气孔关闭，进而导致光合同化物质生产的下降。在华北地区，同域分布的栓皮栎北缘群体和蒙古栎南缘群体的生长与当年冬春降水存在正相关吲，表明生长均受到降水的限制。北京地区麻栎和河北棋盘山地区蒙古栎RWI与全年降水呈显著正相关，辽宁大黑山栓皮栎RWI与当年4月降水及上年6-10月降水呈显著正相关，表明北京地区麻栎、河北棋盘山地区蒙古栎、辽宁大黑山地区栓皮栎生长均受到降水制约。可以认为，不同区域树木生长的限制因子有所不同，但在同一区域，不同树种生长的限制因素也存在差异。如本研究发现北京地区麻栎与河北地区蒙古栎受干旱制约，但在辽宁仙人洞地区3种栎树均与降水、温度无明显相关，表明该地区水热条件较为适宜栎类树种生长。有研究表明，较低的水分可利用性是干旱、半干旱地区森林生态系统发育和林木生长受到限制的主要原因。相对于华北北部，虽然东北南部的辽宁大黑山地区降水更丰沛，但该地区栓皮栎RWI与降水却表现出正相关关系，这与Bimal对红椿（Toona cillaataRoem.）的研究结果一致。这可能与降水的季节性差异有关，比如大黑山地区降水主要集中在6-9月，生长季前期降水较少，水热的不同步限制了树木的初生生长。

为了对未来树木生长作出评估，本研究采用了IPCC第五份评估报告4个温室气体浓度情景中的最理想的RCP2.6情景和最极端的RCP8.5情景下的气候数据进行预测，分别表示到2100年辐射强度为2.6 W·m-2和8.5 W·m-2。RCP2.6情景下，北京·河北地区增温2.6℃，降水量增加156.35 mm；辽宁仙人洞地区增温2.5℃，降水量增加117.22 mm；辽宁大黑山地区增温2.4℃，降水量增加175.54 mm。RCP8.5情景下，北京-河北地区增温6.6℃，降水量增加341.38 mm；辽宁仙人洞地区增温5.3℃，降水量增加318.74mm；辽宁大黑山地区增温5.3℃，降水量增加385.45 mm。根据模型预测，截至本世纪末，北京、辽宁大黑山地区栎类树种生长量将会持续减少，特别是在RCP8.5情景下，栎类树种生长下降愈加明显。这与李宗善对北京东灵山辽东栎（Q.wutaishanica）的预测一致。可能的解释是，北京地区处于半干旱地区，气候变暖导致蒸散量增加，树木的可利用水分降低，水分限制愈加严重。辽宁仙人洞地区栎类树种径向增长预计会增加，可能的原因是辽宁仙人洞地区有适宜的温度和相对丰沛的降水，RWI与气候因子的相关性分析中，显著相关的气候因子均为正相关，因此降水量和气温的增加导致了树木的加速生长。然而，树木生长适宜温度的上限是多少，超过适宜温度对树木的限制作用、立地条件及林分结构等因素也应该进行考虑。

4结论

本研究利用年轮生态学研究了我国暖温带北缘3个栎类树种对气候要素的响应差异，发现同一地区的栎类树种生长对气候响应相似；不同区域的同一栎树的生长对气候响应存在差别，可以认为在北京、河北地区，水分和低温是制约栎类树种生长的主要因素。模型预测还发现，至本世纪末，北京、辽宁大黑山地区栎类树种的生长可能发生衰退，而辽宁仙人洞地区栎类树种的生长则表现出增加趋势。然而，除气候要素外，树木生长还受到立地条件、林分结构与树木年龄等因素及各类因子间交互作用的影响，下一步需要结合更多环境因子展开更深入的研究，以提高模拟预测的可靠性。