健康树

第４０卷第２４期２０２０年１２月

生态学报ＡＣＴＡＥＣＯＬＯＧＩＣＡＳＩＮＩＣＡ

Ｖｏｌ．４０，Ｎｏ．２４

Ｄｅｃ．，２０２０

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

基金项目：京津风沙源人工林与土壤水分互馈机制（２０１７ＺＹ０２）；国家自然科学基金项目（４１８７７５３９）

收稿日期：２０２０⁃０３⁃０４； 网络出版日期：２０２０⁃１１⁃０５

∗通讯作者Ｃｏｒｒｅｓｐｏｎｄｉｎｇａｕｔｈｏｒ．Ｅ⁃ｍａｉｌ：ｘｉｎｚｈｏｎｇｂａｏ＠１２６．ｃｏｍ

ＤＯＩ：１０．５８４６／ｓｔｘｂ２０２００３０４０４０６

刘亚玲，信忠保，李宗善，买尔当·克依木．近４０年河北坝上地区杨树人工林径向生长对气候变化的响应差异．生态学报，２０２０，４０（２４）：９１０８⁃９１１９．

ＬｉｕＹＬ，ＸｉｎＺＢ，ＬｉＺＳ，ＭａｉｅｒｄａｎｇＫｅｙｉｍｕ．ＣｌｉｍａｔｅｅｆｆｅｃｔｏｎｔｈｅｒａｄｉａｌｇｒｏｗｔｈｏｆＰｏｐｕｌｕｓｓｉｍｏｎｉｉｉｎＮｏｒｔｈｗｅｓｔｏｆＨｅｂｅｉｆｏｒｌａｓｔｆｏｕｒｄｅｃａｄｅｓ．Ａｃｔａ

ＥｃｏｌｏｇｉｃａＳｉｎｉｃａ，２０２０，４０（２４）：９１０８⁃９１１９．

近４０年河北坝上地区杨树人工林径向生长对气候变

化的响应差异

刘亚玲

１，信忠保１，∗，李宗善２，买尔当·克依木２

１北京林业大学水土保持学院，北京 １０００８３

２中国科学院生态环境研究中心城市与区域生态国家重点实验室，北京 １０００８５

摘要：运用树木年轮气候学方法，研究近４０年河北坝上地区健康和衰退小叶杨人工林径向生长对气候响应敏感性差异，揭示健

康和衰退杨树生长与气候关系的时间变异规律。结果表明：（１）衰退杨树径向生长对温度、降水等气候因素响应较健康杨树敏

感。衰退杨树年表中的气候信号较强，与当年生长季（４、８—１０月）的气温因素呈显著负相关，与上一年休眠期（９月—当年１

月）和当年生长季（７月）的降水和相对湿度呈显著正相关。健康杨树年表中气候信号较弱，主要与上一年冬季（１２月）和当年

生长季（４月）的气温因素呈显著负相关，与上一年生长末期（８—１１月）降水和相对湿度呈显著正相关。（２）从各年表与帕默尔

干旱指数（ＰＤＳＩ）的响应强度来看，衰退杨树生长更易受夏季干旱胁迫影响。衰退杨树年表与上一年９月—当年３月、６—１０月

的ＰＤＳＩ呈显著正相关，而健康杨树径向生长与ＰＤＳＩ呈弱的正相关。（３）１９７５—２０１７年间，随气温升高，健康和衰退杨树生长

对温度的敏感性下降；健康杨树生长对降水和ＰＤＳＩ的敏感性较为稳定，适应能力强，而衰退杨树生长对降水和ＰＤＳＩ的敏感性

增强，适应能力变弱。综上所述，干旱胁迫是限制衰退杨树生长的主要因素，而健康杨树生长受气候影响较弱，能适应当地气候

条件。衰退杨树对气候变化响应较健康杨树明显，在气候变暖背景下，衰退杨树生长的气候限制因子由温度转变为水分，导致

河北坝上地区遭受干旱灾害时发生退化的趋势更加明显。

关键词：河北坝上；小叶杨；衰退；径向生长；气候响应

ＣｌｉｍａｔｅｅｆｆｅｃｔｏｎｔｈｅｒａｄｉａｌｇｒｏｗｔｈｏｆＰｏｐｕｌｕｓｓｉｍｏｎｉｉｉｎＮｏｒｔｈｗｅｓｔｏｆＨｅｂｅｉｆｏｒ

ｌａｓｔｆｏｕｒｄｅｃａｄｅｓ

ＬＩＵＹａｌｉｎｇ１，ＸＩＮＺｈｏｎｇｂａｏ１，∗，ＬＩＺｏｎｇｓｈａｎ２，ＭａｉｅｒｄａｎｇＫｅｙｉｍｕ２

１ＣｏｌｌｅｇｅｏｆＷａｔｅｒａｎｄＳｏｉｌＣｏｎｓｅｒｖａｔｉｏｎ，ＢｅｉｊｉｎｇＦｏｒｅｓｔｒｙＵｎｉｖｅｒｓｉｔｙ，Ｂｅｉｊｉｎｇ１０００８３，Ｃｈｉｎａ

２ＳｔａｔｅＫｅｙＬａｂｏｒａｔｏｒｙｏｆＵｒｂａｎａｎｄＲｅｇｉｏｎａｌＥｃｏｌｏｇｙ，ＲｅｓｅａｒｃｈＣｅｎｔｅｒｆｏｒＥｃｏ⁃ＥｎｖｉｒｏｎｍｅｎｔａｌＳｃｉｅｎｃｅｓ，ＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ

１０００８５，Ｃｈｉｎａ

Ａｂｓｔｒａｃｔ：Ｄｅｎｄｒｏｃｌｉｍａｔｉｃｍｅｔｈｏｄｓｈａｖｅｂｅｅｎｕｓｅｄｔｏｅｓｔａｂｌｉｓｈｓｔａｎｄａｒｄｃｈｒｏｎｏｌｏｇｉｅｓｆｏｒｔｈｅｈｅａｌｔｈｙａｎｄｔｈｅｄｅｃｌｉｎｉｎｇｔｒｅｅｓ

ｏｆｔｈｅＰｏｐｕｌｕｓｓｉｍｏｎｉｉａｒｔｉｆｉｃｉａｌｆｏｒｅｓｔｉｎｎｏｒｔｈｗｅｓｔｏｆＨｅｂｅｉＰｒｏｖｉｎｃｅｏｖｅｒｔｈｅｐａｓｔｆｏｕｒｄｅｃａｄｅｓ．Ｃｏｒｒｅｌａｔｉｏｎａｎａｌｙｓｅｓｗｅｒｅ

ｃｏｎｄｕｃｔｅｄｂｅｔｗｅｅｎｔｈｅｓｅｃｈｒｏｎｏｌｏｇｉｅｓａｎｄｌｏｃａｌｃｌｉｍａｔｉｃｆａｃｔｏｒｓｔｏｅｓｔａｂｌｉｓｈｔｈｅｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｓｏｆｔｈｅｇｒｏｗｔｈ⁃ｃｌｉｍａｔｅ

ｒｅｌａｔｉｏｎｓｈｉｐｆｏｒｔｈｅｔｈｅｈｅａｌｔｈｙａｎｄｔｈｅｄｅｃｌｉｎｉｎｇｔｒｅｅｓｏｆｔｈｅＰｏｐｕｌｕｓｓｉｍｏｎｉｉｆｏｒｅｓｔ．Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔ（１）ｔｈｅｔｒｅｅ⁃

ｒｉｎｇｇｒｏｗｔｈｏｆｔｈｅｄｅｃｌｉｎｉｎｇｔｒｅｅｓｗａｓｍｏｒｅｓｅｎｓｉｔｉｖｅｔｏｍｅｔｅｏｒｏｌｏｇｉｃａｌｆａｃｔｏｒｓｔｈａｎｔｈａｔｏｆｔｈｅｈｅａｌｔｈｙｔｒｅｅｓ．Ｔｈｅｓｅｒｉｅｓｏｆ

ｔｒｅｅ⁃ｒｉｎｇｗｉｄｔｈｓｏｆｔｈｅｄｅｃｌｉｎｉｎｇｔｒｅｅｓｅｘｈｉｂｉｔｅｄｓｉｇｎｉｆｉｃａｎｔｌｙｎｅｇａｔｉｖｅｃｏｒｒｅｌａｔｉｏｎｓｗｉｔｈｔｈｅｔｅｍｐｅｒａｔｕｒｅｆａｃｔｏｒｓｉｎＡｐｒｉｌａｎｄ

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

ｆｒｏｍＡｕｇｕｓｔｔｏＯｃｔｏｂｅｒｏｆｔｈｅｃｕｒｒｅｎｔｙｅａｒ；ｍｅａｎｗｈｉｌｅ，ｉｔｙｉｅｌｄｅｄｐｏｓｉｔｉｖｅｃｏｒｒｅｌａｔｉｏｎｓｗｉｔｈｔｏｔａｌｐｒｅｃｉｐｉｔａｔｉｏｎａｎｄｒｅｌａｔｉｖｅ

ｈｕｍｉｄｉｔｙｆｒｏｍＳｅｐｔｅｍｂｅｒｏｆｔｈｅｐｒｅｖｉｏｕｓｙｅａｒｔｏＪａｎｕａｒｙａｎｄｆｏｒＪｕｌｙｏｆｔｈｅｃｕｒｒｅｎｔｙｅａｒ．Ｆｏｒｔｈｅｈｅａｌｔｈｙｔｒｅｅｓ，ｔｈｅｓｅｒｉｅｓ

ｗａｓｓｉｇｎｉｆｉｃａｎｔｌｙｎｅｇａｔｉｖｅｌｙｃｏｒｒｅｌａｔｅｄｗｉｔｈｔｈｅｔｅｍｐｅｒａｔｕｒｅｆａｃｔｏｒｓｆｏｒＤｅｃｅｍｂｅｒｏｆｔｈｅｐｒｅｖｉｏｕｓｙｅａｒａｎｄＡｐｒｉｌｏｆｔｈｅｃｕｒｒｅｎｔ

ｙｅａｒａｎｄｐｏｓｉｔｉｖｅｌｙｃｏｒｒｅｌａｔｅｄｗｉｔｈｔｈｅｔｏｔａｌｐｒｅｃｉｐｉｔａｔｉｏｎａｎｄｒｅｌａｔｉｖｅｈｕｍｉｄｉｔｙｆｒｏｍＡｕｇｕｓｔｔｏＮｏｖｅｍｂｅｒｏｆｔｈｅｐｒｅｖｉｏｕｓ

ｙｅａｒ．（２）ＦｒｏｍｔｈｅｃｈｒｏｎｏｌｏｇｉｃａｌｒｅｓｐｏｎｓｅｔｏｔｈｅＰａｌｍｅｒＤｒｏｕｇｈｔＩｎｄｅｘ（ＰＤＳＩ），ｔｈｅｄｅｃｌｉｎｉｎｇｔｒｅｅｓｗｅｒｅｍｏｒｅｓｕｓｃｅｐｔｉｂｌｅ

ｔｏｓｕｍｍｅｒｄｒｏｕｇｈｔｓｔｒｅｓｓ．Ｔｈｅｔｒｅｅ⁃ｒｉｎｇｗｉｄｔｈｓｅｒｉｅｓｏｆｔｈｅｄｅｃｌｉｎｉｎｇｔｒｅｅｓｈａｄａｓｔｒｏｎｇｌｙｐｏｓｉｔｉｖｅｃｏｒｒｅｌａｔｉｏｎｓｗｉｔｈｔｈｅＰＤＳＩ

ｆｒｏｍＳｅｐｔｅｍｂｅｒｏｆｔｈｅｐｒｅｖｉｏｕｓｙｅａｒｔｏＭａｒｃｈｏｆｔｈｅｃｕｒｒｅｎｔｙｅａｒ，ｗｈｅｒｅａｓｆｒｏｍＪｕｎｅｔｏＯｃｔｏｂｅｒｏｆｔｈｅｃｕｒｒｅｎｔｙｅａｒ，ｗｈｉｌｅ

ｔｈｅｈｅａｌｔｈｙｔｒｅｅｓｗｅｒｅｗｅａｋｌｙａｆｆｅｃｔｅｄｂｙｄｒｏｕｇｈｔ．（３）Ａｓｔｈｅｔｅｍｐｅｒａｔｕｒｅｒｏｓｅｆｒｏｍ１９７５，ｔｈｅｔｅｍｐｅｒａｔｕｒｅｓｅｎｓｉｔｉｖｉｔｙｏｆｔｈｅ

ｈｅａｌｔｈｙａｎｄｔｈｅｄｅｃｌｉｎｉｎｇｔｒｅｅｓｄｅｃｒｅａｓｅｄ，ａｎｄｔｈｅＰＤＳＩａｎｄｔｏｔａｌｐｒｅｃｉｐｉｔａｔｉｏｎｓｅｎｓｉｔｉｖｉｔｙｏｆｔｈｅｄｅｃｌｉｎｉｎｇｔｒｅｅｓｉｎｃｒｅａｓｅｄ；

ｍｅａｎｗｈｉｌｅ，ｔｈｅｈｅａｌｔｈｙｔｒｅｅｓｗｅｒｅｓｔａｂｌｅ．Ｉｎｓｕｍｍａｒｙ，ｄｒｏｕｇｈｔｓｔｒｅｓｓｗａｓｔｈｅｍａｉｎｆａｃｔｏｒｌｉｍｉｔｉｎｇｔｈｅｇｒｏｗｔｈｏｆｔｈｅ

ｄｅｃｌｉｎｉｎｇｔｒｅｅｓ，ｗｈｅｒｅａｓｔｈｅｈｅａｌｔｈｙｔｒｅｅｓｗｅｒｅｌｅｓｓａｆｆｅｃｔｅｄｂｙｃｌｉｍａｔｅａｎｄｃｏｕｌｄａｄａｐｔｔｏｌｏｃａｌｃｌｉｍａｔｉｃｃｏｎｄｉｔｉｏｎｓ．Ｔｈｅ

ｄｅｃｌｉｎｉｎｇｔｒｅｅｓｗｅｒｅｍｏｒｅｓｅｎｓｉｔｉｖｅｔｏｃｌｉｍａｔｅｃｈａｎｇｅｔｈａｎｔｈｅｈｅａｌｔｈｙｔｒｅｅｓ．Ｗｉｔｈｔｈｅｗａｒｍｉｎｇａｎｄｄｒｙｉｎｇ，ｔｈｅｃｌｉｍａｔｉｃ

ｌｉｍｉｔｉｎｇｆａｃｔｏｒｏｆｔｈｅｄｅｃｌｉｎｉｎｇｔｒｅｅｓｗｉｌｌｃｈａｎｇｅｆｒｏｍｔｅｍｐｅｒａｔｕｒｅｔｏｍｏｉｓｔｕｒｅ，ｗｈｉｃｈｌｅａｄｓｔｏａｍｏｒｅｐｒｏｎｏｕｎｃｅｄｔｒｅｎｄｏｆ

ｄｅｇｒａｄａｔｉｏｎｄｕｒｉｎｇｄｒｏｕｇｈｔｓｔｒｅｓｓ．

ＫｅｙＷｏｒｄｓ：ＨｅｂｅｉＰｒｏｖｉｎｃｅ；Ｐｏｐｕｌｕｓｓｉｍｏｎｉｉ；ｄｅｃｌｉｎｉｎｇ；ｔｒｅｅｇｒｏｗｔｈ；ｃｌｉｍａｔｅｒｅｓｐｏｎｓｅ

我国西北、华北及东北西部地区，风沙灾害和水土流失严重，为减少自然灾害，改善其生态环境，于１９７８

年批复建设三北防护林体系建设工程

［１］。自工程实施以来，我国三北地区植被覆盖率显著提高，森林覆盖率

由５．０５％增加到１３．５７％，尤其华北地区植被覆盖率增长最快，人工林面积逐年增加

［２⁃４］。河北省张家口坝上

地区杨树防护林是三北工程华北地区的重要组成部分，是预防京津冀沙尘天气发生和生态安全的重要保障，

在防风固沙、降低风速、保持水土、涵养水源、调节小气候等方面发挥着至关重要作用

［５⁃７］。

近年来，该区域杨树防护林出现了明显衰退现象，退化杨树面积占杨树防护林总面积的４／５，其中濒临衰

退和枯死杨树面积达到３．３９

×

１０

４ｈｍ２，占杨树防护林总面积的１／３，主要表现为林木生长不良、大面积枯死和

病虫害频发，防护林防护功能明显下降

［８⁃９］。已有诸多学者从树木生理过程［９⁃１１］、土壤理化性质［１２］、气候变

化

［１３⁃１４］、病虫害［１２，１５］

以及人为因素

［１３］

等方面研究了该地区杨树防护林的衰退问题，但该区域杨树衰退的具

体原因未能明确。树木年轮可以记录其过去径向生长过程

［１６⁃１８］，为研究杨树对生长环境和气候变化的响应

提供了难得契机，可以通过比较健康和衰退杨树的径向生长过程了解其衰退原因。

张北地区属于半干旱地区，受全球气候变化影响显著，自２０世纪８０年代以来，该地区温度持续上升，降

水量减少，呈现暖干化趋势

［１９］，会对当地杨树人工林径向生长产生巨大影响。最新研究表明气候暖干化可能

促进部分地区树木的生长

［２０⁃２１］，也会抑制部分地区树木的生长［２１⁃２２］。气候变化对不同区域不同树种生长的

影响具有很大的区别。如在我国东北地区气候变暖引起土壤水分降低，干旱胁迫严重限制了红松、红皮云

杉

［２３］

的径向生长，西藏地区温度升高限制了拉萨河流域大果圆柏生长，而在西北天山地区，温度升高促进了

雪岭云杉林的生长

［２０］。在同一地区不同物种对气候变化响应也有所不同，如小兴安岭地区黄菠萝径向生长

随温度的升高而增加，而蒙古栎出现降低趋势

［２１］。气候变化对树木生长的影响会随区域和树种的不同而产

生较大差异，需要开展张北地区杨树生长对气候变化的响应研究，对全面了解气候变化背景下树木生长动态

具有重要意义。

利用树木年轮气候学方法，以该地区广泛种植的小叶杨（Ｐｏｐｕｌｕｓｓｉｍｏｎｉｉ）为研究对象，通过采集健康和衰

退杨树树芯，研究了健康和衰退杨树径向生长对气候因素响应的差异以及气候变化对杨树生长的影响。研究

有助于提升当前河北坝上地区三北防护林衰退过程与机制的认识，为全球气候变化背景下坝上地区过熟人工

林经营与改造提供科学参考。

９０１９

２４期 刘亚玲 等：近４０年河北坝上地区杨树人工林径向生长对气候变化的响应差异

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

１ 研究区与研究方法

１．１ 研究区概况

研究区二台镇林场位于河北省张家口市张北县境内，海拔１３７０—１３９０ｍ，面积为６．４ｋｍ２（图１）。地处内

蒙古南缘的坝上高寒区，属中温带大陆性季风气候，光照充足，昼夜温差大，干旱、多风、少水、无霜期短，年平

均气温为１—６℃，年降水量在４００ｍｍ左右，多集中于夏季（图２）。研究区年平均日照时数２８９７．８ｈ，年平均７级以上大风日数３０天左右。该地区土壤风蚀严重，主要地带性土壤为栗钙土，土壤ｐＨ为５．５２—７．３７。植

被以２０世纪７０年代种植的小叶杨（Ｐｏｐｕｌｕｓｓｉｍｏｎｉｉ）为主，造林时使用的是一年生扦插苗，初植密度偏大，后

期抚育管理措施（间伐、肥水等）缺乏，现存杨树树龄基本一致，密度在７６９—１８４９株／ｈｍ

２

之间

［９］。其它防护

林种还包括青杨（Ｐｏｐｕｌｕｓｃａｔｈａｙａｎａ）、北京杨（Ｐｏｐｕｌｕｓｂｅｉｊｉｎｇｅｎｓｉｓ）和樟子松（Ｐｉｎｕｓｓｙｌｖｅｓｔｒｉｓｖａｒ．ｍｏｎｇｏｌｉｃａ）等。

样地中的杨树为小叶杨，包含健康和衰退杨树。

图１ 张北地区树木年轮采样点及气象站分布图

Ｆｉｇ．１ Ｌｏｃａｔｉｏｎｍａｐｏｆｔｒｅｅ⁃ｒｉｎｇｓａｍｐｌｅｓｉｔｅａｎｄｍｅｔｅｏｒｏｌｏｇｉｃａｌｓｔａｔｉｏｎｉｎＮｏｒｔｈｗｅｓｔｅｒｎｏｆＨｅｂｅｉ，Ｃｈｉｎａ

１．２ 样品采集及年表建立

２０１８年７—８月，在张家口市张北县二台镇林场采集健康和衰退的杨树树芯样品。研究以枯梢率来划分

健康和衰退杨树，是指枯树梢长度占树高的百分比（％）。当杨树枯梢率低于２０％时就将其定义为健康杨树，

当枯梢率在２０％之上时将其定义为衰退杨树。用５．１５ｍｍ的生长锥在胸径（１．３ｍ）处沿东西方向钻取树芯。

为对杨树的伤害减至最小，每棵树只钻取一根树芯，置于塑料管内编号封存。共采集杨树样芯２６８根，其中，

健康杨树１９７根，衰退杨树７１根。

将树芯带到实验室后，根据国际上标准方法，晾干后用白乳胶固定在木槽内

［２４］。由于树木材质和人为原

因导致部分样芯断裂严重难以固定粘贴，最终固定了杨树树芯１５８根，其中健康杨树１１４根，衰退杨树４４根。

自然晾干后用４００目、８００目和１２００目的干砂纸依次打磨，直到在显微镜下可以清楚看到杨树的细胞或年轮

界线。利用精度为０．００１ｍｍ树木年轮分析仪ＷｉｎＤＥＮＤＲＯ

ＴＭ（Ｒｅｇｅｎｔｉｎｓｔｒｕｍｅｎｔｉｎｃ．Ｃａｎａｄａ）进行树轮宽度

０１１９

生 态 学 报 ４０卷

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

图２ １９７５—２０１７年张北地区温度降水季节分布

Ｆｉｇ．２ Ｓｅａｓｏｎａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆｔｅｍｐｅｒａｔｕｒｅａｎｄｐｒｅｃｉｐｉｔａｔｉｏｎｉｎ

ＮｏｒｔｈｗｅｓｔｅｒｎｏｆＨｅｂｅｉ，Ｃｈｉｎａ（１９７５—２０１７）

参数的获取。为保证测量的准确性，用ＣＯＦＥＣＨＡ程序

进行交叉定年检验和校准

［２５］。

为消除树龄和其他非气候因素导致的生长趋势，用Ｒ语言中ｄｐｌＲ扩展包对每个树轮宽度序列进行去趋势

（ｍｏｄｉｆｉｅｄＨｕｇｅｒｓｈｏｆｆｃｕｒｖｅ）处理［２６］，最后利用双权重平

均法合成健康与衰退杨树标准年表、差值年表和自回归

年表

［２７］。考虑到标准年表包含更多的低频信息并且剔

除了非气候信号和生长趋势

［１６，２８］，利用健康和衰退杨

树标准年表进行一系列研究。

１．３ 气象数据及统计分析

气象资料来源于中国气象数据网（ｈｔｔｐ：／／ｄａｔａ．ｃｍａ．ｃｎ／），从距离取样点最近的张北气象站（１１４．４２°Ｅ，

４１．０９°Ｎ，海拔为１３９３．３ｍ）获取月降水量（ＰＲＥ）、月平

均气温（ＴＭＰ）、月平均最低气温（ＴＭＮ）、月平均最高气

温（ＴＭＸ）和月平均相对湿度（ＲＨ）气象数据（１９７５—２０１７年）。选取不同年表与气象数据的最大公共时段进

行相关性分析，健康和衰退杨树年表与气象数据的最大公共时段分别为１９７５—２０１７年和１９８２—２０１７年。考

虑到树木生长对气候变化存在滞后性，选用上一年６月到当年１０月的气候数据进行树木年表与气候要素相

关性分析，该过程由ＤｅｎｄｒｏＣｌｉｍ２００２程序

［２９］

实现。

张北地区干旱少水，研究选用了全球气候数据库（ＣＲＵＴＳ３．２２ＧｌｏｂａｌＣｌｉｍａｔｅＤａｔａｂａｓｅ，ｈｔｔｐ：／／ｗｗｗ．ｃｒｕ．ｕｅａ．ａｃ．ｕｋ／）中的帕默尔干旱指数（Ｐａｌｍｅｒｄｒｏｕｇｈｔｓｅｖｅｒｉｔｙｉｎｄｅｘ，ＰＤＳＩ），来更好的反映杨树生长与水分的关

系

［３０］。从距采样点最近的ＣＲＵ网格点（１１４．５０°—１１５．００°Ｅ，４１．００°—４１．５０°Ｎ）选取１９７５—２０１７年ＰＤＳＩ数

据。分析单月和季节性气候因子对杨树径向生长的影响，并利用滑动相关分析揭示健康、衰退杨树生长与气

候关系的稳定性，滑动窗口为２０ａ。

２ 研究结果

２．１ 气候年际变化趋势

该地区１９７５以来年平均气温以０．４１℃／１０ａ（Ｐ＜０．０１）的速度上升（图３），年平均最低气温（０．５０℃／１０ａ，Ｐ＜０．０１）的上升速率明显高于年平均最高气温（０．３５℃／１０ａ，Ｐ＜０．０１）。年降水量和相对湿度呈现降低趋

势，但减少量不明显，出现小范围的波动现象（图３）。１９７５—２０１７年多年平均降水量为（３９１．５±７１．０）ｍｍ，１９７８年降水量最大为５４８．４ｍｍ，１９９７年降水量最少，仅有２４５．２ｍｍ。研究区ＰＤＳＩ下降趋势明显（

－

０．３８／

１０ａ，Ｐ＜０．０１），多年份的ＰＤＳＩ值小于

－

２，该地区出现明显的气候暖干化的特征（图３）。

２．２ 健康和衰退杨树年表统计特征比较

建立了健康和衰退杨树树轮宽度标准年表，年表的起始年份最低样本量均大于５根样芯（图４），建立的

两个年表能够反映研究区健康和衰退杨树的径向生长年际波动特征。年表统计特征表明，健康和衰退杨树年

表的标准差（ＳＤ）分别为１．５４４和１．９００，蕴含较多的气候信息（表１）。衰退杨树的一阶自相关系数（ＡＲ１）较

高（０．６０９），表明其受前一年气候的影响较大，存在“滞后效应”

［１６］。２个年表样芯间相关系数（Ｒｂｔ）分别为

０．３０６和０．４４５，表明衰退杨树不同样芯间树轮宽度变化有很好的一致性。两个年表的样本总体代表性（ＥＰＳ）

均大于０．８５，表明采集的样本量包含的信号基本能代表当地的总体特征

［３１］。健康和衰退杨树年表的信噪比

（ＳＮＲ）分别为３３．９９５和２０．６１７，达到了较高水平，表明包含丰富的气候信号。以上各年表统计特征说明本次

调查所用的树轮资料包含可靠的气候信息，适用于进行树轮气候学分析。

１１１９

２４期 刘亚玲 等：近４０年河北坝上地区杨树人工林径向生长对气候变化的响应差异

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

图３ １９７５—２０１７年张北各气象要素的变化趋势

Ｆｉｇ．３ ＣｌｉｍａｔｉｃｃｈａｎｇｉｎｇｔｒｅｎｄｏｆＮｏｒｔｈｗｅｓｔｅｒｎｏｆＨｅｂｅｉ，Ｃｈｉｎａｆｒｏｍ１９７５ｔｏ２０１７

表１ 张北坝上地区树轮采样点概况、标准年表主要特征参数及共同区间分析结果

Ｔａｂｌｅ１ Ｓｉｔｅｉｎｆｏｒｍａｔｉｏｎ，ｃｈｒｏｎｏｌｏｇｙｓｔａｔｉｓｔｉｃｓａｎｄｒｅｓｕｌｔｓｏｆｃｏｍｍｏｎｉｎｔｅｒｖａｌｓｐａｎａｎａｌｙｓｉｓｏｆｓｔａｎｄａｒｄｔｒｅｅ⁃ｒｉｎｇｃｈｒｏｎｏｌｏｇｙｆｒｏｍｔｈｅ

ＮｏｒｔｈｗｅｓｔｅｒｎｏｆＨｅｂｅｉ，Ｃｈｉｎａ

类型Ｔｙｐｅ地点Ｌｏｃａｔｉｏｎ海拔Ｅｌｅｖａｔｉｏｎ／ｍ

时间段Ｔｉｍｅ

ｌｅｎｇｔｈ

样芯数Ｎｕｍｂｅｒ

ｏｆｃｏｒｅｓ

平均树

轮宽度Ｍｅａｎ／ｍｍ

标准差ＳＤ

一阶

自相关ＡＲ１

样芯间相

关系数Ｒｂｔ

信噪比ＳＮＲ

样本群体

表达信号ＥＰＳ

健康杨树Ｔｈｅｈｅａｌｔｈｙｐｏｐｌａｒ

ｔｒｅｅｓ

１１４．８８°Ｅ，

４１．３４°Ｎ

１３７９．０１９７５—２０１７１０４２．５３３１．５４４０．５９００．３０６３３．９９５０．９７１

衰退杨树ＴｈｅＤｅｃｌｉｎｉｎｇｐｏｐｌａｒ

ｔｒｅｅｓ

１１４．８８°Ｅ，

４１．３４°Ｎ

１３７９．０１９８２—２０１７３５２．５７０１．９０００．６０９０．４４５２０．６１７０．９５４

气象站Ｍｅｔｅｏｒｏｌｏｇｉｃａｌ

ｓｔａｔｉｏｎ

１１４．４２°Ｅ，

４１．０９°Ｎ

１３９３．３１９５６—２０１８———————

ＣＲＵ格点

ＣＲＵｇｒｉｄ

１１４．５０°Ｅ，

４１．００°Ｎ

—１９０１—２０１８———————

Ｍｅａｎ：平均树轮宽度Ｍｅａｎｔｒｅｅ⁃ｒｉｎｇｗｉｄｔｈ（ｍｍ）；ＳＤ：标准差Ｓｔａｎｄａｒｄｄｅｖｉａｔｉｏｎ；ＡＲ１：一阶自相关Ｆｉｒｓｔ⁃ｏｒｄｅｒａｕｔｏｃｏｒｒｅｌａｔｉｏｎ；Ｒｂｔ：样芯间相关系数Ｍｅａｎｉｎｔｅｒ⁃

ｓｅｒｉｅｓｃｏｒｒｅｌａｔｉｏｎ；ＳＮＲ：信噪比Ｓｉｇｎａｌ⁃ｔｏ⁃ｎｏｉｓｅｒａｔｉｏ；ＥＰＳ：样本群体表达信号Ｅｘｐｒｅｓｓｐｏｐｕｌａｔｉｏｎｓｉｇｉｎａｌ

２．３ 气候因素对杨树径向生长的影响

２．３．１ 温度对杨树径向生长的影响

健康和衰退杨树径向生长与气温多呈负相关关系，衰退杨树年表表达的温度信号较健康杨树强，衰退杨

２１１９

生 态 学 报 ４０卷

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

图４ 张北地区健康和衰退杨树标准年表

Ｆｉｇ．４ Ｓｔａｎｄａｒｄｔｒｅｅ⁃ｒｉｎｇｃｈｒｏｎｏｌｏｇｙｏｆｈｅａｌｔｈｙａｎｄｄｅｃｌｉｎｉｎｇＰｏｐｕｌｕｓｓｉｍｏｎｉｉｆｏｒｅｓｔｉｎＮｏｒｔｈｗｅｓｔｏｆＨｅｂｅｉ，Ｃｈｉｎａ

绿色实线为健康杨树标准年表，红色实线为衰退杨树标准年表，蓝色虚线为样本量

树对气温响应敏感性强于健康杨树。健康杨树年轮宽度与上年１０月—当年１月和当年３—６月温度显著相

关（图５，Ｐ＜０．０５），衰退杨树年轮宽度与上年８—１０月、当年４—６月以及当年８—１０月温度显著相关（Ｐ＜０．０５）。健康杨树与上年１２月和当年４月温度指标呈显著负相关（Ｐ＜０．０５），其中与当年４月平均温度、平均

最高温度以及平均最低温度达到极显著相关水平（Ｐ＜０．０１）。衰退杨树与当年４月、７月、９月和１０月平均气

温，７月和１０月最高气温以及９月和１０月最低气温呈显著负相关（图５，Ｐ＜０．０５），其中与当年４月平均温度、

平均最高温度以及平均最低温度达到极显著相关水平（Ｐ＜０．０１）。两个年表与单月尺度４月份温度因素的相

关系数比季节性尺度４—６月温度因素高。

２．３．２ 降水和相对湿度对杨树径向生长的影响

健康和衰退杨树径向生长与降水和相对湿度多呈正相关关系，衰退杨树年表表达的降水和相对湿度信号

较健康杨树强。健康杨树年轮宽度与上年８—１１月降水以及上年８月—当年１月相对湿度呈显著正相关（图６，Ｐ＜０．０５）；衰退杨树与上年９月—当年１月、５—９月降水以及上年９月—当年１月、６—９月相对湿度呈

显著正相关（Ｐ＜０．０５）。衰退杨树与上年９月、１１月和当年７月降水以及上年９月和当年７月相对湿度呈显

著正相关（Ｐ＜０．０５），其中与当年７月降水和相对湿度达到极显著相关水平（Ｐ＜０．０１）；而健康杨树仅与上年９

月降水和相对湿度呈显著正相关（图６，Ｐ＜０．０５）。衰退杨树年表与７月单月尺度降水的相关系数比５—７月

季节尺度的相关系数高，但年表与上年１１月—当年１月的降水呈极显著正相关，与单月尺度降水相关性较

低。健康杨树年表与上年８—１０月降水以及上年８—１１月相对湿度呈极显著正相关，而与单月尺度降水和相

对湿度相关性仅９月达到显著水平。

２．３．３ 干旱对杨树径向生长的影响

帕默尔干旱指数（ＰＤＳＩ）基于水平衡方程的供需概念，根据降水和温度数据以及土壤含水量计算，是反映

土壤湿度的代用指标，在干旱区的相关研究中被广泛应用，其数值大小代表着土壤水分有效性状况

［３２］。衰退

杨树年表与多月份ＰＤＳＩ呈显著正相关，而健康杨树径向生长与ＰＤＳＩ相关性不显著（图７）。衰退杨树径向生

长与ＰＤＳＩ在多月份（上年９月—当年３月、当年６月—１０月）呈显著正相关（Ｐ＜０．０５），其中与上年１１月以及

当年夏季（６—８月）ＰＤＳＩ的相关性（ｒ

＝

０．４４１—０．５０９，Ｐ＜０．０１）达到极显著水平。

３１１９

２４期 刘亚玲 等：近４０年河北坝上地区杨树人工林径向生长对气候变化的响应差异

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

图５ 张北杨树年表与单月和三个月（季节）尺度温度要素的相关关系

Ｆｉｇ．５ ＣｏｒｒｅｌａｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔｓｏｆＰｏｐｕｌｕｓｓｉｍｏｎｉｉｃｈｒｏｎｏｌｏｇｉｅｓｂｅｔｗｅｅｎｔｅｍｐｅｒａｔｕｒｅｆａｃｔｏｒｓ

∗代表Ｐ＜０．０５，∗∗代表Ｐ＜０．０１，月份前的“

－

”表示上一年，如

－

６表示上一年６月份（单月尺度），

－

６

－－

８表示上一年６月份—８月份（季节

尺度）

２．４ 杨树径向生长对气候响应敏感性年际变化趋势

衰退杨树与温度因子的负相关关系较健康杨树显著，研究时段内健康和衰退杨树径向生长对温度因素的

敏感性呈下降趋势（图８）。健康和衰退杨树年表对４—６月平均温度的响应分别在１９８２—２００１年和１９７５—２０１４年间呈显著负相关，之后负相关性逐渐降低；两个年表与４月平均气温在１９７５—２０１６年间都呈显著负相

关，并且相关性年际间呈明显下降趋势。两个年表与４月和４—６月最高温度的响应有相同的变化趋势。健

康杨树年表对８—１０月平均温度的响应在１９８２—２０１７年间呈弱负相关，而衰退杨树在１９８６—２０１７年间与之

呈显著负相关（Ｐ＜０．０５），在２０１０年之后负相关性逐渐降低。

衰退杨树与降水和ＰＤＳＩ的正相关性较健康杨树强，研究时间段内衰退杨树径向生长对降水和ＰＤＳＩ的

敏感性逐步增强，而健康杨树径向生长对降水和ＰＤＳＩ的敏感性较为稳定，呈现弱正相关性（图８）。由健康和

衰退杨树年表对５—７月降水量和６—８月ＰＤＳＩ的滑动相关关系可见，衰退杨树年表在１９８２—２０１７年间始终

与降水量和ＰＤＳＩ呈正相关，且对降水和ＰＤＳＩ的敏感性逐渐增强；健康杨树年表在研究时间段内对两者的敏

感性较衰退杨树弱，对降水的敏感性在２０１４年时由负相关转为正相关，和ＰＤＳＩ的敏感性也在逐步加强但是

呈现弱相关性。由健康和衰退杨树年表对７月降水量和８月的ＰＤＳＩ的滑动相关关系可以看出，衰退杨树年

表在１９８２—２０１７年间与降水量和ＰＤＳＩ呈正相关性有上升趋势，与降水量在１９９２—２０１１年之后达到显著性

水平，与ＰＤＳＩ在１９９７—２００６年之后达到显著正相关；健康杨树对降水量和ＰＤＳＩ的响应敏感性在１９８２—

４１１９

生 态 学 报 ４０卷

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

图６ 张北杨树年表与单月和三个月（季节）尺度降水和相对湿度的相关关系

Ｆｉｇ．６ ＣｏｒｒｅｌａｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔｓｏｆＰｏｐｕｌｕｓｓｉｍｏｎｉｉｃｈｒｏｎｏｌｏｇｉｅｓｂｅｔｗｅｅｎｐｒｅｃｉｐｉｔａｔｉｏｎａｎｄｒｅｌａｔｉｖｅｈｕｍｉｄｉｔｙ

∗代表Ｐ＜０．０５，∗∗代表Ｐ＜０．０１，月份前的“

－

”表示上一年，如

－

６表示上一年６月份（单月尺度），

－

６

－－

８表示上一年６月份—８月份（季节

尺度）

图７ 张北杨树年表与单月和三个月（季节）尺度ＰＤＳＩ的相关关系

Ｆｉｇ．７ ＣｏｒｒｅｌａｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔｓｏｆＰｏｐｕｌｕｓｓｉｍｏｎｉｉｃｈｒｏｎｏｌｏｇｉｅｓｂｅｔｗｅｅｎＰＤＳＩ

ＰＤＳＩ：帕默尔干旱指数Ｐａｌｍｅｒｄｒｏｕｇｈｔｉｎｄｅｘ；∗代表Ｐ＜０．０５，∗∗代表Ｐ＜０．０１，月份前的“

－

”表示上一年，如

－

６表示上一年６月份（单月尺

度），

－

６

－－

８表示上一年６月份—８月份（季节尺度）

５１１９

２４期 刘亚玲 等：近４０年河北坝上地区杨树人工林径向生长对气候变化的响应差异

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

２０１７年间都较弱。

３ 讨论

３．１ 健康和衰退杨树径向生长与气候因子关系

河北坝上地区杨树人工林中衰退杨树主要受温度和降水条件的制约，而健康杨树对温度响应敏感，对水

分变化响应相对较弱。该结果在河北地区内也有体现，由于受夏季水热条件的限制，河北木兰围场油松与５、６月气温呈显著负相关，而与６月降水呈显著正相关［３３］；太行山南段侧柏径向生长受春末夏初水热条件的限

制

［３４］；冀北山地落叶松受夏季气温和降水两个因子共同影响［３５］；６００ｍｍ年降水量等值线是中国树木径向生

长对气候变化响应的重要空间边界，在该线的北部，树轮宽度通常与降水呈正相关，与温度呈负相关

［３６］，张北

地区位于６００ｍｍ年降水量等值线北面，说明结果较为准确。

衰退杨树年表与生长季前４月、生长季７月和生长季末９—１０月温度呈显著负相关关系，健康杨树与生

长季前４月温度呈显著负相关。该地区小叶杨生长发育关键时期主要集中在５—９月，一般５月开始发芽，６—８月为生长期，９月末进入落叶休眠期［９］，在生长季中多月份平均气温过高对小叶杨生长产生胁迫，从而

限制小叶杨径向生长。健康杨树对当年４月最低气温以及衰退杨树对４月、９月和１０月最低温度的响应较

敏感，４月份是小叶杨生长季前，当温度过低时，杨树生长受到限制，树木各种生理活动都会延迟

［３７］。同时，

受９月最低气温的影响，杨树提前进入休眠期，树木生长季变短，促使杨树形成窄树轮

［３８］。健康杨树与当年４

月最高气温以及衰退杨树与当年４、７和１０月最高气温呈显著负相关，生长季前高温会使土壤蒸发增加，导致

土壤的有效水分显著下降，进入生长季时没有充足的水分

［３９］，促使生长季缩短。在生长季期间，高温促使树

木蒸腾作用加强，呼吸作用加强，消耗的养分过多，在生长季中树木径向生长减缓

［３７，４０］。

杨树径向生长与降水存在“滞后效应”，上年８—１１月降水对衰退杨树的生长有促进作用，当生长季后期

降水充足时可用于补充根系层的土壤水分，是下年生长季初期杨树生长的主要水分来源

［４１］。张北雨热同期，

在夏季时温度基本稳定在２０℃上下，树木各项生理活动进入旺盛时期，在持续高温情况下树木水分蒸发快

速

［４２］，此时降水成为衰退杨树生长的主要限制因素。本地区７月份降水量最大，降水成为杨树快速生长的有

利因素，所以与７月降水量正相关关系最显著。

总之，温度和降水对衰退杨树生长的限制作用明显强于健康杨树，可能与杨树衰退后生理功能下降有关。

在高温干旱胁迫的影响下，杨树径向生长得不到充足的水分，杨树在低水势的条件下易形成导管的栓塞，它自

身会调节牺牲部分构件来维持正常生理需水，如发生枯树梢现象降低植物蒸腾作用强度

［４３］。发生林分枯梢

现象说明对干旱的适应策略已超出其正常生长范围，气候变化极易影响衰退杨树生长，对气象因子响应敏感。

当杨树遭遇干旱胁迫，它会调节气孔来减少水分蒸发，气孔关闭会阻碍植物与外界进行气体交换，ＣＯ

２

供应不

足，使植物体内碳水化合物储存不够，同时植物呼吸作用和新陈代谢消耗体内碳水化合物，而使本体处于“饥

饿”状态，会使杨树径向生长下降生成窄轮

［４４］。从生理过程方面理解杨树衰退过程还需要做出进一步研究。

３．２ 气候暖干化对健康和衰退杨树的影响

根据张北气象站数据，在本研究区域内，１９５６—２０１７年间年均降水量为（３９１．５±７１．０）ｍｍ，为半干旱区域。

近６０年来温度呈现出显著上升的趋势，增长速率为０．３６℃／１０ａ，多年份的ＰＤＳＩ值小于

－

２，达到了严重干旱

的程度

［４５］，这说明该地区温度快速升高引起干旱，气候呈现暖干化趋势。衰退杨树与ＰＤＳＩ在多月份呈显著

正相关，径向生长下降，说明其深受水分胁迫的影响，而健康杨树与ＰＤＳＩ的相关性较弱，该地区的健康杨树

对干旱具有较好的适应性。

在全球气候变暖的背景下，北半球树木生长对气候因素的敏感性随时间发生了明显变化

［４６］。树木生长

对干旱胁迫响应程度呈上升趋势，而对温度的响应出现下降趋势，致使干旱胁迫对树木生长产生了较大影

响

［４７］。在本研究中衰退杨树年表与气象数据的滑动相关分析表明，年表与春季温度负相关性有减弱的趋势，

与夏季降水和ＰＤＳＩ的正相关性有增强的趋势，可能是气温升高引起的干旱胁迫严重影响了衰退杨树的生

６１１９

生 态 学 报 ４０卷

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

图８ 健康和衰退杨树年表与主要月份气候要素的滑动相关分析

Ｆｉｇ．８ Ｍｏｖｉｎｇｃｏｒｒｅｌａｔｉｏｎｆｕｎｃｔｉｏｎａｎａｌｙｓｉｓｃｏｍｐｕｔｅｄｂｅｔｗｅｅｎｔｈｅｔｗｏｃｈｒｏｎｏｌｏｇｉｅｓａｎｄｔｈｅｓｅａｓｏｎａｌａｎｄｍｏｎｔｈｌｙｃｌｉｍａｔｅｄａｔａｆｏｒｔｈｅｍｏｓｔ

ｉｍｐｏｒｔａｎｔｍｏｎｔｈｓｄｕｒｉｎｇｔｈｅｐｅｒｉｏｄ１９７５—２０１７

滑动相关分析窗口设置为２０ａ；显著性相关水平为Ｐ＜０．０５；４表示当年４月份（单月尺度），４—６表示当年４月份—６月份（季节尺度）

长，研究结果与Ｂａｂｓｔ等

［４７］

的结论一致。小兴安岭蒙古栎生长与夏季最低温度负相关性有减弱的趋势，与夏

季相对湿度的正相关性有增强的趋势，其径向生长出现衰退趋势

［２１］，本研究结果与其一致。伊春地区红松与

红皮云杉随温度的上升，高温对其径向生长的抑制作用加强将不利于两种树木的生长，使其产生生长衰退现

７１１９

２４期 刘亚玲 等：近４０年河北坝上地区杨树人工林径向生长对气候变化的响应差异

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

象，与本研究结果一致

［２３］。有研究指出气温升高会促进长白山地区红松以及天山地区雪岭云杉的生长，本研

究与其结果相反，可能是气候区不同而导致的

［２０，４８］。健康杨树年表与春季和夏季温度的负相关性以及夏季

降水和ＰＤＳＩ的正相关性有增强的趋势，但对健康杨树径向生长影响不大，说明当地气候适宜健康杨树的

生长。

４ 结论

以河北坝上地区三北防护林工程中小叶杨农田防护林为研究对象，利用树木年轮学方法建立健康和衰退

杨树年轮宽度标准年表，探讨了健康杨树和退化杨树对温度、降水等气候因素响应的差异以及年际变化特征。

发现健康和衰退杨树径向生长对气候因素的响应存在差异。温度是二者生长的主要限制因子，降水的限制作

用相对较弱，衰退杨树较健康杨树对气候因素响应更加敏感。生长季前（４月）温度抑制健康和衰退杨树的生

长，生长季（７月）降水促进衰退杨树的生长。高温引起的干旱胁迫明显抑制衰退杨树的生长，而健康杨树受

其影响较小，健康杨树对张北地区气候具有较好的适应性。随着温度的持续升高，健康和衰退杨树与生长季

温度、降水的相关性存在明显差异，衰退杨树径向生长由温度胁迫转变为水分胁迫，衰退杨树生长将更加受限

制。衰退小叶杨对气候响应较健康小叶杨更加敏感，对干旱灾害敏感，适应能力下降，因此，干旱是小叶杨衰

退乃至死亡的重要驱动因素。

参考文献（Ｒｅｆｅｒｅｎｃｅｓ）：

［１］ 李建树．三北防护林体系工程发展现状及对策．世界林业研究，１９９４，（１）：６４⁃６９．

［２］ 申丽娜，孙艳玲，杨艳丽，景悦．基于ＮＤＶＩ的三北防护林工程区植被覆盖度变化图谱特征．环境科学与技术，２０１７，４０（４）：７０⁃７７．

［３］ ＤｕａｎＨＣ，ＹａｎＣＺ，ＴｓｕｎｅｋａｗａＡ，ＳｏｎｇＸ，ＬｉＳ，ＸｉｅＪＬ．Ａｓｓｅｓｓｉｎｇｖｅｇｅｔａｔｉｏｎｄｙｎａｍｉｃｓｉｎｔｈｅｔｈｒｅｅ

－

ｎｏｒｔｈｓｈｅｌｔｅｒｆｏｒｅｓｔｒｅｇｉｏｎｏｆｃｈｉｎａｕｓｉｎｇ

ＡＶＨＲＲＮＤＶＩｄａｔａ．ＥｎｖｉｒｏｎｍｅｎｔａｌＥａｒｔｈＳｃｉｅｎｃｅｓ，２０１１，６４（４）：１０１１⁃１０２０．

［４］ ＮｉｕＱＦ，ＸｉａｏＸＭ，ＺｈａｎｇＹ，ＱｉｎＹＷ，ＤａｎｇＸＨ，ＷａｎｇＪ，ＺｏｕＺＨ，ＤｏｕｇｈｔｙＲＢ，ＢｒａｎｄｔＭ，ＴｏｎｇＸＷ，ＨｏｒｉｏｎＳ，ＦｅｎｓｈｏｌｔＲ，ＣｈｅｎＣ，

ＭｙｎｅｎｉＲＢ，ＸｕＷＨ，ＤｉＧＺ，ＺｈｏｕＸＭ．Ｅｃｏｌｏｇｉｃａｌｅｎｇｉｎｅｅｒｉｎｇｐｒｏｊｅｃｔｓｉｎｃｒｅａｓｅｄｖｅｇｅｔａｔｉｏｎｃｏｖｅｒ，ｐｒｏｄｕｃｔｉｏｎ，ａｎｄｂｉｏｍａｓｓｉｎｓｅｍｉａｒｉｄａｎｄ

ｓｕｂｈｕｍｉｄＮｏｒｔｈｅｒｎＣｈｉｎａ．ＬａｎｄＤｅｇｒａｄａｔｉｏｎ＆Ｄｅｖｅｌｏｐｍｅｎｔ，２０１９，３０（１３）：１６２０⁃１６３１．

［５］ 郑春雅，许中旗，马长明，底凯，程月明，孙守家，闫腾飞．冀西北坝上地区退化防护林的土壤性质．水土保持学报，２０１６，３０（１）：

２０３⁃２０７．

［６］ ＬｉｕＷＨ，ＺｈｕＪＪ，ＪｉａＱＱ，ＺｈｅｎｇＸ，ＬｉＪＳ，ＬｏｕＸＤ，ＨｕＬＬ．ＣａｒｂｏｎｓｅｑｕｅｓｔｒａｔｉｏｎｅｆｆｅｃｔｓｏｆｓｈｒｕｂｌａｎｄｓｉｎＴｈｒｅｅ

－

Ｎｏｒｔｈｓｈｅｌｔｅｒｂｅｌｔｆｏｒｅｓｔ

ｒｅｇｉｏｎ，Ｃｈｉｎａ．ＣｈｉｎｅｓｅＧｅｏｇｒａｐｈｉｃａｌＳｃｉｅｎｃｅ，２０１４，２４（４）：４４４⁃４５３．

［７］ ＷａｎｇＸＭ，ＺｈａｎｇＣＸ，ＨａｓｉＥ，ＤｏｎｇＺＢ．ＨａｓｔｈｅＴｈｒｅｅＮｏｒｔｈｓＦｏｒｅｓｔＳｈｅｌｔｅｒｂｅｌｔＰｒｏｇｒａｍｓｏｌｖｅｄｔｈｅｄｅｓｅｒｔｉｆｉｃａｔｉｏｎａｎｄｄｕｓｔｓｔｏｒｍｐｒｏｂｌｅｍｓｉｎ

ａｒｉｄａｎｄｓｅｍｉａｒｉｄＣｈｉｎａ？ＪｏｕｒｎａｌｏｆＡｒｉｄＥｎｖｉｒｏｎｍｅｎｔｓ，２０１０，７４（１）：１３⁃２２．

［８］ 邢海富．坝上杨树防护林死亡原因浅析．河北林业科技，２０１５，（１）：７６⁃７７．

［９］ 苗博，孟平，张劲松，何方杰，孙守家．基于稳定同位素和热扩散技术的张北杨树水分关系差异．应用生态学报，２０１７，２８（７）：

２１１１⁃２１１８．

［１０］ 张欢，曹俊，王化冰，宋波，贾国栋，刘自强，余新晓，曾佳．张北地区退化杨树防护林的水分利用特征．应用生态学报，２０１８，２９（５）：

１３８１⁃１３８８．

［１１］ 杜海峰．河北坝上杨树退化的原因及改造技术求索．环境与发展，２０１９，３１（７）：１９５⁃１９５，１９７⁃１９７．

［１２］ 程月明，张愿，崔硕，马长明．张北县杨树种群衰退原因初探———以盘城营为例．林业资源管理，２０１５，（１）：１５２⁃１５５，１６７⁃１６７．

［１３］ 孙守家，李春友，何春霞，张劲松，孟平．基于树轮稳定碳同位素的张北杨树防护林退化原因解析．应用生态学报，２０１７，２８（７）：

２１１９⁃２１２７．

［１４］ 李晓航，刘晓龙．干旱导致防护林衰退死亡研究进展．防护林科技，２０１６，（１２）：５９⁃６０．

［１５］ 郑春雅，许中旗，马长明，孙守家，闫腾飞．冀西北坝上地区杨树防护林退化的影响因素．林业资源管理，２０１８，（１）：９⁃１５，１４７⁃１４７．

［１６］ ＦｒｉｔｔｓＨＣ．ＴｒｅｅＲｉｎｇｓａｎｄＣｌｉｍａｔｅ．ＮｅｗＹｏｒｋ：ＡｃａｄｅｍｉｃＰｒｅｓｓ，１９７６．

［１７］ 吴祥定．树木年轮与气候变化．北京：气象出版社，１９９０．

［１８］ 彭剑峰，刘玉振，王婷．神农山白皮松不同龄组年轮———气候关系及ＰＤＳＩ重建．生态学报，２０１４，３４（１３）：３５０９⁃３５１８．

［１９］ 杨思遥，孟丹，李小娟，吴新玲．华北地区２００１—２０１４年植被变化对ＳＰＥＩ气象干旱指数多尺度的响应．生态学报，２０１８，３８（３）：

１０２８⁃１０３９．

［２０］ ＱｉＺＨ，ＬｉｕＨＹ，ＷｕＸＣ，ＨａｏＱ．Ｃｌｉｍａｔｅ⁃ｄｒｉｖｅｎｓｐｅｅｄｕｐｏｆａｌｐｉｎｅｔｒｅｅｌｉｎｅｆｏｒｅｓｔｇｒｏｗｔｈｉｎｔｈｅＴｉａｎｓｈａｎＭｏｕｎｔａｉｎｓ，ＮｏｒｔｈｗｅｓｔｅｒｎＣｈｉｎａ．Ｇｌｏｂａｌ

ＣｈａｎｇｅＢｉｏｌｏｇｙ，２０１５，２１（２）：８１６⁃８２６．

［２１］ 韩金生，赵慧颖，朱良军，张远东，李宗善，王晓春．小兴安岭蒙古栎和黄菠萝径向生长对气候变化的响应比较．应用生态学报，２０１９，

８１１９

生 态 学 报 ４０卷

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

３０（７）：２２１８⁃２２３０．

［２２］ 石松林，靳甜甜，刘国华，王东波，王景升，李魁．气候变暖抑制西藏拉萨河大果圆柏树木生长．生态学报，２０１８，３８（２４）：８９６４⁃８９７２．

［２３］ 苑丹阳，赵慧颖，李宗善，朱良军，国淼，张远东，王晓春．伊春地区红松和红皮云杉径向生长对气候变化的响应．生态学报，２０２０，４０

（４）：１１５０⁃１１６０．

［２４］ ＳｔｏｋｅｓＭＡ，ＳｍｉｌｅｙＴＬ．ＡｎＩｎｔｒｏｄｕｃｔｉｏｎｔｏＴｒｅｅＲｉｎｇＤａｔｉｎｇ．Ｃｈｉｃａｇｏ：ＵｎｉｖｅｒｓｉｔｙｏｆＣｈｉｃａｇｏＰｒｅｓｓ，１９６８．

［２５］ ＨｏｌｍｅｓＲＬ．Ｃｏｍｐｕｔｅｒ⁃ａｓｓｉｓｔｅｄｑｕａｌｉｔｙｃｏｎｔｒｏｌｉｎｔｒｅｅ⁃ｒｉｎｇｄａｔｉｎｇａｎｄｍｅａｓｕｒｅｍｅｎｔ．Ｔｒｅｅ⁃ＲｉｎｇＢｕｌｌｅｔｉｎ，１９８３，４３：５１⁃６７．

［２６］ ＢｕｎｎＡＧ．ＡｄｅｎｄｒｏｃｈｒｏｎｏｌｏｇｙｐｒｏｇｒａｍｌｉｂｒａｒｙｉｎＲ（ｄｐｌＲ）．Ｄｅｎｄｒｏｃｈｒｏｎｏｌｏｇｉａ，２００８，２６（２）：１１５⁃１２４．

［２７］ 邵雪梅，吴祥定．华山树木年轮年表的建立．地理学报，１９９４，６１（２）：１７４⁃１８１．

［２８］ 吴祥定，邵雪梅．中国秦岭地区树木年轮密度对气候响应的初步分析．应用气象学报，１９９４，５（２）：２５３⁃２５６．

［２９］ ＢｉｏｎｄｉＦ，ＷａｉｋｕｌＫ．ＤＥＮＤＲＯＣＬＩＭ２００２：ａＣ

＋＋

ｐｒｏｇｒａｍｆｏｒｓｔａｔｉｓｔｉｃａｌｃａｌｉｂｒａｔｉｏｎｏｆｃｌｉｍａｔｅｓｉｇｎａｌｓｉｎｔｒｅｅ⁃ｒｉｎｇｃｈｒｏｎｏｌｏｇｉｅｓ．Ｃｏｍｐｕｔｅｒｓ＆

Ｇｅｏｓｃｉｅｎｃｅｓ，２００４，３０（３）：３０３⁃３１１．

［３０］ ＨａｒｒｉｓＩ，ＪｏｎｅｓＰＤ，ＯｓｂｏｒｎＴＪ，ＬｉｓｔｅｒＤＨ．Ｕｐｄａｔｅｄｈｉｇｈ⁃ｒｅｓｏｌｕｔｉｏｎｇｒｉｄｓｏｆｍｏｎｔｈｌｙｃｌｉｍａｔｉｃｏｂｓｅｒｖａｔｉｏｎｓ⁃ｔｈｅＣＲＵＴＳ３．１０Ｄａｔａｓｅｔ．Ｉｎｔｅｒｎａｔｉｏｎａｌ

ＪｏｕｒｎａｌｏｆＣｌｉｍａｔｏｌｏｇｙ，２０１４，３４（３）：６２３⁃６４２．

［３１］ ＷｉｇｌｅｙＴＭＬ，ＢｒｉｆｆａＫＲ，ＪｏｎｅｓＰＤ．Ｏｎｔｈｅａｖｅｒａｇｅｖａｌｕｅｏｆｃｏｒｒｅｌａｔｅｄｔｉｍｅｓｅｒｉｅｓ，ｗｉｔｈａｐｐｌｉｃａｔｉｏｎｓｉｎｄｅｎｄｒｏｃｌｉｍａｔｏｌｏｇｙａｎｄｈｙｄｒｏｍｅｔｅｏｒｏｌｏｇｙ．

ＪｏｕｒｎａｌｏｆＣｌｉｍａｔｅａｎｄＡｐｐｌｉｅｄＭｅｔｅｏｒｏｌｏｇｙ，１９８４，２３（２）：２０１⁃２１３．

［３２］ ＢｈｕｙａｎＵ，ＺａｎｇＣ，ＭｅｎｚｅｌＡ．ＤｉｆｆｅｒｅｎｔｒｅｓｐｏｎｓｅｓｏｆｍｕｌｔｉｓｐｅｃｉｅｓｔｒｅｅｒｉｎｇｇｒｏｗｔｈｔｏｖａｒｉｏｕｓｄｒｏｕｇｈｔｉｎｄｉｃｅｓａｃｒｏｓｓＥｕｒｏｐｅ．Ｄｅｎｄｒｏｃｈｒｏｎｏｌｏｇｉａ，

２０１７，４４：１⁃８．

［３３］ 崔明星，何兴元，陈玮，陈振举，周长虹，吴涛．河北木兰围场油松年轮生态学的初步研究．应用生态学报，２００８，１９（１１）：２３３９⁃２３４５．

［３４］ 王亚军，张永，邵雪梅．河北太行山南段树木年轮指示的１６７年来相对湿度变化研究．生态学报，２０１９，３９（１２）：４５７０⁃４５７８．

［３５］ 范慧涛，陈立标，陈忠震，谷建才．冀北山地油松、落叶松径向生长对气象因子的响应．中南林业科技大学学报，２０１９，３９（１）：５２⁃５７．

［３６］ ＬｉｕＹ，ＳｏｎｇＨＭ，ＳｕｎＣＦ，ＳｏｎｇＹ，ＣａｉＱＦ，ＬｉｕＲＳ，ＬｅｉＹ，ＬｉＱ．Ｔｈｅ６００⁃ｍｍｐｒｅｃｉｐｉｔａｔｉｏｎｉｓｏｌｉｎｅｄｉｓｔｉｎｇｕｉｓｈｅｓｔｒｅｅ⁃ｒｉｎｇ⁃ｗｉｄｔｈｒｅｓｐｏｎｓｅｓｔｏ

ｃｌｉｍａｔｅｉｎＣｈｉｎａ．ＮａｔｉｏｎａｌＳｃｉｅｎｃｅＲｅｖｉｅｗ，２０１９，６（２）：３５９⁃３６８．

［３７］ 贺敏，魏江生，石亮，周梅，赵鹏武．大兴安岭南段山杨径向生长和死亡对区域气候变化的响应．生态学杂志，２０１８，３７（１１）：

３２３７⁃３２４４．

［３８］ 张先亮，崔明星，马艳军，吴涛，陈振举，丁玮航．大兴安岭库都尔地区兴安落叶松年轮宽度年表及其与气候变化的关系．应用生态学

报，２０１０，２１（１０）：２５０１⁃２５０７．［３９］ 常永兴，陈振举，张先亮，白学平，赵学鹏，李俊霞，陆旭．气候变暖下大兴安岭落叶松径向生长对温度的响应．植物生态学报，２０１７，

４１（３）：２７９⁃２８９．

［４０］ 王婷，于丹，李江风，马克平．树木年轮宽度与气候变化关系研究进展．植物生态学报，２００３，２７（１）：２３⁃３３．

［４１］ 韦景树，李宗善，焦磊，陈维梁，伍星，王晓春，王帅．黄土高原羊圈沟小流域人工物种和自然物种径向生长对气候变化的响应差异．生

态学报，２０１８，３８（２２）：８０４０⁃８０５０．［４２］ 王童，孙玉军，乔晶晶．将乐林场马尾松树轮宽度对气候变化的响应．北京林业大学学报，２０１９，４１（９）：３０⁃３９．

［４３］ 韦景树，李宗善，冯晓玙，张园，陈维梁，伍星，焦磊，王晓春．黄土高原人工刺槐林生长衰退的生态生理机制．应用生态学报，２０１８，２９

（７）：２４３３⁃２４４４．

［４４］ ＭｃＤｏｗｅｌｌＮ，ＰｏｃｋｍａｎＷＴ，ＡｌｌｅｎＣＤ，ＢｒｅｓｈｅａｒｓＤＤ，ＣｏｂｂＮ，ＫｏｌｂＴ，ＰｌａｕｔＪ，ＳｐｅｒｒｙＪ，ＷｅｓｔＡ，ＷｉｌｌｉａｍｓＤＧ，ＹｅｐｅｚＥＡ．Ｍｅｃｈａｎｉｓｍｓｏｆ

ｐｌａｎｔｓｕｒｖｉｖａｌａｎｄｍｏｒｔａｌｉｔｙｄｕｒｉｎｇｄｒｏｕｇｈｔ：ｗｈｙｄｏｓｏｍｅｐｌａｎｔｓｓｕｒｖｉｖｅｗｈｉｌｅｏｔｈｅｒｓｓｕｃｃｕｍｂｔｏｄｒｏｕｇｈｔ？ＮｅｗＰｈｙｔｏｌｏｇｉｓｔ，２００８，１７８（４）：

７１９⁃７３９．

［４５］ 卫捷，陶诗言，张庆云．Ｐａｌｍｅｒ干旱指数在华北干旱分析中的应用．地理学报，２００３，５８（Ｓ１）：９１⁃９９．

［４６］ ＢｒｉｆｆａＫＲ，ＳｃｈｗｅｉｎｇｒｕｂｅｒＦＨ，ＪｏｎｅｓＰＤ，ＯｓｂｏｒｎＴＪ，ＳｈｉｙａｔｏｖＳＧ，ＶａｇａｎｏｖＥＡ．Ｒｅｄｕｃｅｄｓｅｎｓｉｔｉｖｉｔｙｏｆｒｅｃｅｎｔｔｒｅｅ⁃ｇｒｏｗｔｈｔｏｔｅｍｐｅｒａｔｕｒｅａｔ

ｈｉｇｈｎｏｒｔｈｅｒｎｌａｔｉｔｕｄｅｓ．Ｎａｔｕｒｅ，１９９８，３９１（６６６８）：６７８⁃６８２．

［４７］ ＢａｂｓｔＦ，ＢｏｕｒｉａｕｄＯ，ＰｏｕｌｔｅｒＢ，ＴｒｏｕｅｔＶ，ＧｉｒａｒｄｉｎＭＰ，ＦｒａｎｋＤＣ．Ｔｗｅｎｔｉｅｔｈｃｅｎｔｕｒｙｒｅｄｉｓｔｒｉｂｕｔｉｏｎｉｎｃｌｉｍａｔｉｃｄｒｉｖｅｒｓｏｆｇｌｏｂａｌｔｒｅｅｇｒｏｗｔｈ．

ＳｃｉｅｎｃｅＡｄｖａｎｃｅｓ，２０１９，５（１）：ｅａａｔ４３１３，ｄｏｉ：１０．１１２６／ｓｃｉａｄｖ．ａａｔ４３１３．

［４８］ 高露双，王晓明，赵秀海．长白山过渡带红松和鱼鳞云杉径向生长对气候因子的响应．植物生态学报，２０１１，３５（１）：２７⁃３４．

９１１９

２４期 刘亚玲 等：近４０年河北坝上地区杨树人工林径向生长对气候变化的响应差异