干旱事件对华北落叶松和油松生态弹性的影响差异.pdf

第 46 卷第 4 期2023 年 7 月河 北 农 业 大 学 学 报JOURNAL OF HEBEI AGRICULTURAL UNIVERSITYVol.46 No.4Jul.2 0 2 3干旱事件对华北落叶松和油松生态弹性的影响差异张博奕，李 雪，张先亮（河北农业大学 林学院，河北 保定 071000）摘要：为了探究干旱事件对人工林主要树种树木生长的影响，本研究以华北落叶松（Larix principis-rupprechtii）以及油松（Pinus tabulaeformis）人工林为研究对象，分析了 2 个树种径向生长对不同气候因子的响应以及单年、连年干旱事件影响下的生态弹性差异。结果表明：不同区域树种径向生长对气候因子和标准化降水蒸散指数(Standardized precipitation evaportranspiration index,SPEI)的响应存在差异，且与 SPEI 的相关性明显强于其他因子，随着干旱程度的增加华北落叶松的抵抗力稳定性呈现下降趋势，而油松较为稳定。单年干旱事件影响下的种间稳定性差异增大，连年干旱事件影响下的种间稳定性差异先增大后减小，在连年干旱事件的影响下，华北落叶松抵抗力明显下降，因此注重连年干旱下华北落叶松人工林的经营与管理，以保持其树木生长的稳定。关 键 词：华北落叶松；油松；径向生长；气候因子；干旱中图分类号：S791.22 开放科学（资源服务）标识码（OSID）:文献标志码：AInfluences of drought events on ecological resilience of Larix principis-rupprechtii and Pinus tabulaeformisZHANGBoyi,LIXue,ZHANGXianliang（College of Forestry,Hebei Agricultural University,Baoding 071000,China）Abstract:This study used Larix principis-rupprechtii and Pinus tabulaeformis plantation to explore the effects of drought events on the growth of main tree species in planted forest.The radial growth of two tree species were analyzed in response to different climatic factors and the differences of ecological elasticity were evaluated under the influence of single and successive drought events.The results showed that there were differences in the response of different tree species to regional climatic factors.Drought is the main limiting factor for the radial growth of trees in the study area.With the increase of drought intensity,the resistance stability of Larix principis-rupprechtii decreased while Pinus tabulaeformis showed stable resistance,and the difference of interspecific stability increased under the influence of annual drought events.The difference of interspecific stability under the influence of successive drought events increased at first and then decreased,and the resistance of Larix principis-rupprechtii decreased obviously under the influence of successive drought events.Therefore,we should pay attention to the management of Larix principis-rupprechtii plantation under continuous drought in order to maintain the stability of its tree growth.Keywords:Larix principis-rupprechtii,Pinus tabulaeformis,radial growth,climatic factors,drought收稿日期：2022-09-26基金项目：河北省中央引导地方科技发展项目（226Z6801G）.第一作者：张博奕（1995），男，河北石家庄人，硕士研究生，主要从事树轮生态学研究.E-mail：通信作者：张先亮（1985），男，山东菏泽人,博士，教授，主要从事树轮生态学研究.E-mail：本刊网址：http:/文章编号：1000-1573(2023)04-0065-09DOI：10.13320/ki.jauh.2023.006066第 46 卷河 北 农 业 大 学 学 报近年来我国北方干旱的频率和强度都有不断增强的趋势1。随着暖干化的加剧，干旱胁迫引起树木死亡率的增加2，进而威胁着森林的稳定性。由于人工林易受干旱、极端气候的影响，干旱加剧会对人工林树木的径向生长产生更严重的影响3。在闽东南地区，随着干旱胁迫增强，马尾松(Pinus massoniana)遭受干旱的影响加强，部分树木出现死亡4。连续的干旱使得内蒙古大兴安岭的山杨(Poplar)更新出现困难，面临着消失的风险5。人工林各个树种之间对干旱的响应存在明显差异6，云杉(Spruce)比冷杉(Abies)对干旱扰动的适应能力差，干旱增强时云杉的生产力可能会大幅度下降7。强烈的极端干旱事件导致黄土高原东南部的油松和白皮松(Pinus bungeand)径向生长出现抑制8。干旱程度加剧使树木抗旱性降低9，西班牙干旱地区的阿勒颇松(Pinus halepensis)比在湿润地区的树木具有更差的抵抗力，且该研究区树木的抵抗力随干旱程度的增加而下降10。生态弹性是树木抵抗干扰及受干扰后恢复到原初状态的能力，干旱事件会影响到树木的生态弹 性11。连年干旱比单年干旱对树木生长造成更严重的影响。马尾松在连续 2 次的干旱事件中表现出比杉木（Cunninghamia lanceolata)更强的恢复能力12。在连续 3 次极端干旱事件之后，伊比利亚半岛东北部 3 种针叶树黑松(Pinus thunbergii)、地中海白松(Pinus halepensis)及欧洲赤松(Pinus sylvestris)对干旱的脆弱性增加13，我国华北地区面对干旱日益严重的现状14，其生态稳定性面临更严重的威胁。华北落叶松和油松作为华北地区人工林的主要造林树种，都属于气候敏感型树种15。华北落叶松、油松人工林在长期干旱频发的影响下出现了林分退化、自然更新差、“小老树”、树梢干枯和林下植被稀疏等现象，导致区域内生态环境进行着恶行循环，进一步加速了森林资源的枯竭16。本研究为明确不同人工林主要优势树种树木径向生长与气候因子的关系，量化径向生长应对干旱事件能力，比较分析不同树种应对干旱的生态弹性的差异，以期探索研究区内华北落叶松与油松生态弹性对单年、连年干旱事件的响应差异，为预测不同树种在经历干旱后的动态变化提供科学依据。材料与方法1.1 研究区概况研究区属温带大陆性季风气候，夏天炎热多雨，春秋多风地干，冬季寒冷少雪，采样点位于华北地区(36 26 42 40 N，111 24 118 14 E)，年均温在5.312.5，7月最高温可以达到27.232.4，1 月最低气温-7.9 -17.8，年均降水除大同市广灵县 288 mm 外，其他区域均在 352 373 mm，降水主要集中于 78 月份，海拔位于 1 320 1 400 m 之 间。主 要 树 种 有 华 北 落 叶 松(Larix principis-rupprechtii)，油松(Pinus tabuliformis)，樟子松(Pinus sylvestris)，云杉(Picea asperata Mast.)，白桦(Betula platyphylla)，蒙古栎(Quercus mongolica)等。1.2 树轮数据在研究区选择 4 个采样点，分别为临汾市霍州市（112 E,36.6 N），吕 梁 市 交 城 县（111.6 E,37.7 N），大 同 市 广 灵 县（114.3 E,39.6 N），承德市围场县（117.4 E,42.3 N），共采集华北落叶松和油松 186 棵，样芯 371 个，其中华北落叶松200 芯（100 棵）；油松 171 芯（86 棵）。按照树木年轮样本处理的基本程序对样本进行晾干、固定、打磨、定年后用精度为 0.001 mm 的 LinTab6 年轮分析测量仪测量样芯，将测完的数据用 COFECHA 程序进行检验，对未通过的样芯重新定年和测量，直到所有的样芯都通过交叉定年检验。使用 ARSTAN程序运用线性函数或负指数函数拟合方法去除年龄趋势得到树轮宽度标准化年表（STD）。1.3 气象数据选取距采样点最近气象站（霍州市气象站、交城市气象站、广灵市气象站、围场县气象站）的月值气象数据（19512020 年），包括月平均最高温，月平均最低温，月均温，月总降水量，气象数据来源于中国气象数据网（http:/）。使用 PDSI 作为地域干旱的代表，其原理为对水分偏差进行适当的加权可生成可用于时空比较的指数17计算公式如下：Xj=Zj3+0.897Xj-1（1）Z=dK（2）SPEI 由 Vicente-Serrano 等人在 2010 年提出，Beguera 等人在 2014 年改进,是在 SPI 的理论基础上,把潜在蒸散发对干旱的影响也加以考虑 SPEI 除67第 4 期了继承 SPI 空间一致性和多时间尺度的优点,又考虑到了潜在蒸散发对干旱的影响,从而能够更好地发挥对干旱的监测和评估作用18。SPEI 引入了潜在蒸散发包含了温度信息，具有多尺度性的特点19。PETi=16(H10Ti)m（3）Di=Pi-PETi（4）F(X)=1+(x-)-1（5）=(1+1)(1-1)(0-21)（6）=61-0-6221-0（7）=0-(1+1)(1-1)（8）SPEI=-1+d1+d22+d33c0+c1+c22（9）=-2lnP （10）式中：PETi为第 i 月潜在蒸散量，mm；Pi为月降水量，mm；Ti为月平均温度，；H 为年热量常 数；m 为 常 数，m=0.49+0.179H-0.000 077 1H2+0.000 000 675H3；为尺度函数，为形状函数，为 origin 参数，P 为累计概率密度；F(x)为 log-logistic概率分布的累计函数；为阶乘函数；0、1、2 为数据序列 Di 的概率加权矩；c0、c1、c2、d1、d2、d3均为常数；为距离加权矩。若 P0.5,则 P=F(x)；若 P0.5,则 P=1-F(x)。全球的格点的 SPEI 和 PDSI 数据集被很多研究使用，具有很高的可靠度，因此本研究中各采样区标准化降水蒸散发指数（SPEI）和帕默尔干旱指数（PDSI）数据序列为 19012020 年分辨率为 0.50*0.50 的格点数据下载自 KNMI 网站（http:/climexp.knmi.nl）。1.4 数据处理1.4.1 相关性计算 STD 年表保留有更多的低频信号20，故基于各采样点 STD 年表与月平均最高温，月平均最低温，月均温，月总降水量，SPEI 以及生长季组合月进行相关性分析。多月份组合对树木生长的影响比单个月份更有气候代表性，因此我们计算了多月份组合与年表的相关性。为了研究组合月树木径向生长对 SPEI 的响应随时间的变化，用年轮数据作步长为 20 年的滑动相关。1.4.2 干旱事件划分 根据参考气象干旱等级对SPEI干湿等级的划分21：-0.5SPEI为正常，-1 SPEI-0.5 为轻度干旱，-1.5 SPEI-1 为中度干旱，-2 SPEI-1.5 为重度干旱，SPEI-2为特别干旱。以干旱发生作为划分即 SPEI -0.5，将各个干旱事件和对应的年表状态标出。1.4.3 生态弹性计算 生态弹性包括抵抗力和恢复力：抵抗力为系统吸收外来干扰并保持其原有状态的能力，恢复力为生态系统从干扰中恢复能力。抵抗力和恢复力的计算公式12如下：Resistance(Rt)=PreDrDr（11）Recovery(Rc)=DrPostDr（12）式中：Resistance(Rt)为抵抗力，Recovery(Rc)为恢复力，Dr=干旱年份的平均树轮宽度，PreDr=每次干旱事件发生时前 3 年的平均树轮宽度，PostDr=每次干旱事件发生时后 3 年的平均树轮宽度。Rt1说明干旱胁迫导致树木生长下降，Rt=1 说明树木完全抵抗干旱胁迫，Rt1 说明干旱胁迫导致树木生长上升，Rc1 说明干旱胁迫导致树木生长水平进一步下降，Rc=1 说明树木完全干旱后生长水平不变，Rc1 说明干旱事件后树木生长得到一定程度恢复。利用 STD 年表计算出所有干旱事件影响下各树种的抵抗力和恢复力，接着将干旱事件分为单年、连年干旱事件使用每根芯的年轮宽度计算各树种的生态弹性，并用 SPSS 对种间抵抗力和恢复力进行显著性分析并在图中进行显著性标注。2 结果与分析2.1 年表的基本统计特征华北落叶松和油松年表统计特征如表 1 所示，可以看出所有采样点样本代表性（EPS）均大于 0.85，年表质量高，可信度高；且信噪比（SNR）都比较高，说明样本中包含大量的环境信息。华北落叶松的平均敏感度（0.266、0.394、0.542、0.222）普遍大于油松（0.126、0.238、0.339、0.199），说明华北落叶松的径向生长的年际变化相对更大。吕梁市交城县的油松（0.427）一阶序列自相关（AC1）高于落叶松（0.019），说明此地的油松的径向生长更容易受到前一年气候因素的影响，围场县塞罕坝华北落叶松、临汾市霍州市油松同样更容易受到前一年气候因素的影响。张博奕，等：干旱事件对华北落叶松和油松生态弹性的影响差异68第 46 卷河 北 农 业 大 学 学 报表 采样点统计特征Table 1 Statistic characteristics of standard chronologies样地名称 Site name样地编号 Site ID树芯/棵数Core/Tree时段 Time span年龄 Length 样本代表性EPS标准差 SD平均敏感度MS信噪比SNR一阶序列自相关 AC1海拔/m Altitude霍州市L144/2219912018280.9770.2490.26643.4410.1861 340霍州市Y140/2019582018610.9110.1930.12610.2080.4901 340交城县L248/2419812018380.9860.3110.39470.1330.0191 320交城县Y242/2119822018370.9890.2790.23888.8050.4271 320广灵县L368/3419862018330.9910.4000.542112.114-0.0221 370广灵县Y346/2319492018700.9830.2800.33959.4870.0031 370围场县L440/2019642020570.9720.2250.22234.3330.3741 400围场县Y443/2219682020530.9510.2130.19919.3450.3441 4002.2 气候年轮宽度的相关关系霍州市地区油松年表与上年 9 月和当年 5 月的降水达到显著正相关，与上年 12 月月均温和上年 12月、当年 5 月月平均最高温成显著负相关关系；交城县华北落叶松与各个气候因子的响应均未达到显著水平（图 1）。广灵县油松年表与上年 11 月月均温成显著正相关关系与当年 5 月月均温成显著负相关关系与当年 2 月、5 月月均最高温成显著负相关关系；围场县华北落叶松年表与当年 5 月降水成显著正相关关系，除此之外，其他月份都没有显著性相关，可见这些对气候因子响应不敏感的树种，气候因子没有显著影响到径向生长。霍州市华北落叶松的年轮宽度与 68 月的月平均最低温呈显著正相关关系，与降水和 SPEI 多呈负相关关系，但均未达到显著。交城县的油松，广灵县华北落叶松以及围场县油松的年表与全年 SPEI 都呈正相关关系且在 68 月达到显著正相关，与其他气候因子的响应并不突出，说明SPEI 是研究区树木径向生长的主要限制因子。2.3 树轮宽度年表与干旱指数的滑动相关分析干旱对霍州市华北落叶松人工林和油松人工林68 月径向生长的影响在减弱（图 2），对吕梁市交城县油松人工林 68 月径向生长的影响有增强的趋势且滑动年份都达到了显著水平（P0.05），大同市广灵县的华北落叶松和围场县塞罕坝油松的滑动响应都比另一树种响应明显，但随时间变化大同市广灵县的华北落叶松 68 月树木径向生长对SPEI 的响应系数呈“中间高两头低”的变化趋势，围场县塞罕坝油松 68 月树木径向生长对 SPEI 的响应系数呈“中间低两头高”的变化趋势，整体为下降趋势。2.4 单年及连年干旱事件临 汾 市 霍 州 市 地 区 单 年 干 旱 事 件：1992、1995、1997、2009、2011 年（图 3）。吕梁市交城县地区单年干旱事件：1987、1992、1995、2009、2011、2015 年。大同市广灵县地区单年干旱事件：1989、1992、2018 年。承德市围场县塞罕坝地区单年干旱事件：1983、1989、2000、2002 年。临汾市霍州市地区连年干旱事件：19992002 年，20052007 年，20142015 年。吕梁市交城县地区连年干旱事件：19972001 年，20052007 年。大同市广灵县地区连年干旱事件：19992002 年，20052007 年，20092011 年，20142015 年。承 德 市围场县塞罕坝地区连年干旱事件：20042011 年，20142015 年，20172018 年。总体上看，连续干旱事件发生在 1997 年之后，四个区域公共干旱事件是 2000 年，20052007 年，2009 年，2011 年，2015 年。2.5 生态稳定性分析各区域树种生态稳定性存在差异（图 4）。霍州市华北落叶松抵抗力均值（r=1.05）比油松（r=0.99）略高，可见华北落叶松抵抗力略强；交城县油松的抵抗力均值（r=0.96）比华北落叶松（r=0.87）高约 10%，此区域油松抵抗力更强。围场县华北落叶松抵抗力均值（r=1.04）与油松抵抗力均值（r=1.05）基本相同，但华北落叶松抵抗力中位数高于油松，相比之下华北落叶松抵抗力更强一些。广灵县华北落叶松的抵抗力均值（r=0.72）低于油松抵抗力均值（r=0.99）约 37%，华北落叶松的抵抗力稳定性明显低于油松。整体可以看出，随地区干旱程度的升高，华北落叶松的抵抗力明显下降，下降约 40%，恢复力并无明显变化，相比之下油松的抵抗力和恢复力较为稳定且在 1 值附近，可见油松面对干旱事件时相比华北落叶松表现出更强的生态 稳定性。69第 4 期Y4L4Y3L3Y2L2Y1L1Y4L4Y3L3Y2L2Y1L1Y4L4Y3L3Y2L2Y1L1Y4L4Y3L3Y2L2Y1L1Y4L4Y3L3Y2L2Y1L1abcde1.00.50.0-0.5-1.0月份MonthP9P11P10P1201020304050607080906-08注：L1：霍州市华北落叶松，Y1：霍州市油松，L2：交城县华北落叶松，Y2：交城县油松，L3：广灵县华北落叶松，Y3：广灵县油松，L4：围场县塞罕坝华北落叶松，Y4：围场县塞罕坝油松。a：月平均最高温，b：月平均最低温，c：月均温，d：月总降水量，e：标准化降水蒸散指数，表示在 P 0.05 水平显著相关，下同。图 1 各采样点 STD 年表与气候因子相关性分析Fig.1 Correlation between STD chronology and climatic factors in each sampling site张博奕，等：干旱事件对华北落叶松和油松生态弹性的影响差异70第 46 卷河 北 农 业 大 学 学 报0.888 00.666 00.444 00.222 00.000 0Y1L1Y2L2Y3L3Y4L4198220011983200219842003198520041986200519872006198820071989200819902009199120101992201119932012199420131995201419962015199720161998201719992018注：相关系数大于 0.444 达到显著水平（P0.05），横坐标为滑动窗口图 2 组合月年表与 SPEI 相关性滑动相关Fig.2 Sliding correlation between chronology and seasonally SPEISPEI1.51.00.50.0543210-1-219921996199419982000200820022010200420122016200620142018落叶松油松SPEI11.51.00.50.0-0.5-1.043210-1-219921984199619881994198619981990200020082002201020042012201620062014201821.51.00.50.0-0.5-1.043210-11992198619961990199419881998 20002008200220102004201220162006201420181.51.00.50.0-0.543210-1-21992198419821996198819941986199819902000200820022010200420122016200620142018年份Year年表Tree-ring standard34注：1:霍州市，2：交城县，3：广灵县，4：围场县塞罕坝，单竖贯穿式虚线：单年干旱事件，双竖夹击型虚线：连年干旱事件。图 3 各地区单年、连年干旱事件标出Fig.3 Single and successive drought events in various regions are marked71第 4 期25%75%1.5IQR 内的范围中位线均值异常值25%75%1.5IQR 内的范围中位线均值异常值1.81.61.41.21.00.80.60.40.22.01.51.00.5抵抗力Resistance恢复力RecoveryL1L2L4L3Y1Y2Y4Y3L1L2L4L3Y1Y2Y4Y3bbaaaaaaaaaaabbb样地Site样地Site注：不同小写字母表示同一组内多重比较结果差异显著，下同。图 4 种间抵抗力和恢复力稳定性分析Fig.4 Stability analysis of resistance and recovery for different speceis2.6 单年、连年干旱事件影响下的种间差异将干旱事件分为单年干旱事件和连年干旱事件分析各地树种的生态弹性（图 5），随着干旱程度增强，单年干旱事件影响下华北落叶松抵抗力呈现波动略微下降趋势，所有单年干旱事件抵抗力均值下降约 23.7%，连年干旱事件影响下华北落叶松抵抗力出现明显的持续下降，所有连年干旱事件抵抗力均值下降约 73.4%。而油松所有干旱事件的抵抗力和恢复力的均值都较为平稳且在1值附近及以上，整体上可以看出无论是在单年干旱事件还是连年干旱事件影响下，研究区油松生态稳定性高于华北落叶松。2.52.01.51.00.50.02.52.01.51.00.50.0抵抗力Resistance抵抗力Resistance恢复力RecoveryL1L2L4L3Y1Y2Y4Y3L1L2L4L3Y1Y2Y4Y3样地Site样地Site25%75%1.5IQR 内的范围中位线均值数据25%75%1.5IQR 内的范围中位线均值数据25%75%1.5IQR 内的范围中位线均值数据25%75%1.5IQR 内的范围中位线均值数据2.52.01.51.00.50.0L1L2L4L3Y1Y2Y4Y3样地SiteDLbbbbbbbbbaaaaaaaaaaaaaaaaaaaaaaaDLL1L2L4L3Y1Y2Y4Y3样地Site3.02.52.01.51.00.50.0恢复力Recovery注：D：单年干旱事件，L：连年干旱事件。图5 单年、连年干旱事件稳定性分析Fig.5 Stability analysis of drought events in single year and multiple years张博奕，等：干旱事件对华北落叶松和油松生态弹性的影响差异72第 46 卷河 北 农 业 大 学 学 报3 讨论3.1 不同树种径向生长对气候因子响应分析前一年 9 月大同市广灵县地区华北落叶松径向生长与月均最高温成显著负相关关系与降水成显著正相关关系，这在生理学上被称为“时滞效应”22。上年 9 月充足的降雨量提高了土壤含水率，营养物质得到积累，为第 2 年春季树木的生长提供有利条件，此时树木已经基本停止生长温度的升高增加水分以及营养物质的消耗23。不同地区对树木径向生长的限制因子也不尽相同24，气候变暖后树木对单一因子的敏感性降低，局部地区会由单一因子或多因子引发的干旱胁迫成为研究区树种径向生长的主要影响因子25。干旱是影响研究区树木生长的重要限制因子，且其生长受多因素共同制约。生长季的开始受控于最低温26，春季升温有利于形成层细胞分裂的提前27，而根系及土壤解冻相对较慢，造成早春生理干旱，随着温度的升高，超过了树木本身正常需求的阈值范围，导致树叶气孔关闭光合速率降低土壤中的水分蒸腾速率升高，植物蒸腾速率和呼吸作用加剧，生理干旱逐渐积累，干旱胁迫也会越来越严重28。3.2 干旱事件对生态弹性影响分析干旱对森林的影响有所增加29，吕梁市交城县，围场县塞罕坝，大同市广灵县 3 个区域干旱程度依次增强，受 SPEI 显著正相关树种的抵抗力稳定性低于另一树种，苏佳佳也指出对干旱越敏感的森林表现出越弱的抵抗力30且连续的干旱事件之后，更容易受到干旱的影响31。相比单年干旱事件的影响对华北落叶松抵抗力稳定性的影响，连年干旱影响是导致华北落叶松抵抗力稳定性下降的重要原因。连年干旱年会产生“累积效应”32，增加了华北落叶松对干旱的脆弱性，而华北落叶松的恢复力并未出现明显的上升趋势，在连年干旱事件中干旱程度的增加，树木会消耗更多储存的干物质量从而降低抵抗力，后通过光合作用进行恢复33。干旱事件影响下种间的生态弹性差异与不同树种应对干旱采取的生理生态策略有关，华北落叶松更优先满足根系生长调节而后地上生长，干旱胁迫加重，对地上的部分的生长影响就更大34，油松的木质部栓塞脆弱性小且由分子小有横壁的管胞组成，应对干旱胁迫时不易发生气穴，保持水分运输的稳定35，所以，在单年干旱影响下，随着干旱程度的增加，种间生态稳定性差异就会越来越明显。连年干旱事件影响下抵抗力和恢复力的种间差距先增加后降低且恢复力的种差下降极为明显几乎降至0，可能在连年干旱达到中度干旱（-1.5SPEI-1）的时候，人为对人工林经营管理减弱了其稳定性的种间差异。4 结论研究区内华北落叶松、油松人工林对气候因子的响应情况以及抵御干旱事件的能力存在差异，干旱是影响华北落叶松、油松人工林生长季径向生长的主要因子。单年干旱事件中，华北落叶松、油松人工林的种间稳定性差异随地区干旱程度的增大而增大，干旱程度较大的区域，华北落叶松人工林的抵抗力有明显下降趋势而油松人工林的生态弹性保持得相对稳定；连年干旱事件成为华北落叶松人工林抵抗力下降的重要因素。因此随着暖干化日益严重，在今后人工林的经营管理中，应更加关注华北落叶松受连年干旱事件的影响以及两树种对有效水的利用，以更好的应对干旱事件，保持长期的生态适应性。参考文献：1 马柱国,符淙斌,杨庆,等.关于我国北方干旱化及其转折性发生J.大气科学,2018,42(4):951-961.2 谢安,孙永罡,白人海.中国东北近 50 年干旱发展及对全球气候变暖的响应J.地理学报,2003,58(S1):75-82.3 Bai M W，Dong Z P,Chen D,et al.Different responses of the radial growth of the planted and natural forests to climate change in humid subtropical ChinaJ.Geografiska Annaler:Series A,Physical Geography,2020,102(3):235-246.4 李颖俊,方克艳,白毛伟,等.闽东南地区马尾松古树对气候变化和虫灾的生态弹性J.应用生态学报,2021,32(10):3539-3547.5 刘畅.内蒙古大兴安岭不同森林类型树木死亡研究D.呼和浩特:内蒙古农业大学,2016.6 Lyu S N,Wang X C,Zhang Y D,et al.Different responses of Korean pine(Pinus koraiensis)and Mongolia oak(Quercus mongolica)growth to 73第 4 期recent climate warming in northeast ChinaJ.Dendrochronologia,2017,45:113-122.7 Seidl R,Thom D,Kautz M,et al.Forest disturbances under climate changeJ.Nature Climate Change,2017,7(6):395-402.8 张启,闫明,梁寒雪.山西省长治市过去 150 年森林的生长抑制和释放历史J.生态学报,2017,37(9):3115-3123.9 李晓琴,张凌楠,曾小敏,等.黄土高原中部针叶树与 灌木径向生长对气候的响应差异J.生态学报,2020,40(16):5685-5697.10 Gazol A,Ribas M,Gutirrez E,et al.Aleppo pine forests from across Spain show drought-induced growth decline and partial recoveryJ.Agricultural and Forest Meteorology,2017,232:186-194.11 张晓英,王飞,铁牛.兴安落叶松林径向生长与气候因子的关系J.林业调查规划,2021,46(6):106-113.12 曹新光,胡红兵,李颖俊,等.亚热带人工和天然马尾松、杉木林生长对干旱的生态弹性差异J.应用生态学报,2021,32(10):3531-3538.13 Serra-Maluquer X,Mencuccini M,Martnez-Vilalta J.Changes in tree resistance,recovery and resilience across three successive extreme droughts in the northeast Iberian PeninsulaJ.Oecologia,2018,187(1):343-354.14 杨思遥,孟丹,李小娟,等.华北地区 2001 2014 年植被变化对 SPEI 气象干旱指数多尺度的响应.生态学报,2018,38(3):1028-1039.15 Chen Z J,Zhang X L,He X Y,et al.Response of radial growth to warming and CO2 enrichment in southern Northeast China:a case of Pinus tabulaeformisJ.Climatic Change,2015,130(4):559571.16 熊千志,杜恩在,薛峰,等.塞罕坝地区人工针叶林径向生长对水热条件的响应J.生态学报,2022,42(13):5371-5380.17 Palmer W C.Meteorological droughtJ.Weather Bureau Research Paper,1965,30:45-58.18 Vicente-Serrano S M,Begueria S,Lopez-Moreno J I,et al.A multiscalar drought index sensitive to global warming:the standardized precipitation evapotran-spiration indexJ.Journal of Climate,2010(23):1696-1718.19 王小雪,王恒,张俊飞,等.塞罕坝林区华北落叶松径向生长对气候变化的响应J.林业与生态科学,2022,37(2):192-197.20 Eriksson M.Methods of dendrochronology:applications in the environmental sciencesJ.Forest Science,1991,37(2):734-735.21 中国气象局政策法规司.气象标准汇编,气象干旱等级(GB/T20481 2006)M.北京:气象出版社,2006:33-34.22 于大炮,王顺忠,唐立娜,等.长白山北坡落叶松年轮年表及其与气候变化的关系J.应用生态学报,2005,16(1):14-20.23 王亚军,高尚玉,马玉贞,等.基于树轮重建的宁夏河东沙地西部公元 1899 年来年降水量的变化J.干旱区地理,2010,33(3):377-384.24 杨镒如,张茗珊,张凌楠,等.秦岭中西部油松径向生长对气候因子的响应差异研究J.生态学报,2022,42(4):1474-1486.25 郑勤莹.帽儿山不同坡位水曲柳和红松径向生长及其对气候变化的响应D.哈尔滨:东北林业大学,2021.26 Li X X,Liang E Y,Griar J,et al.Age dependence of xylogenesis and its climatic sensitivity in Smith fir on the south-eastern Tibetan PlateauJ.Tree Physiology,2013,33(1):48-56.27 Deslauriers A,Rossi S,Anfodillo T,et al.Cambial phenology,wood formation and temperature thresholds in two contrasting years at high altitude