高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

吉林北部山区长白落叶松林径向生长对气候干暖化的响应

陈彬杭 温晓示 张树斌 柴世品 孙晗 王襄平

陈彬杭, 温晓示, 张树斌, 柴世品, 孙晗, 王襄平. 吉林北部山区长白落叶松林径向生长对气候干暖化的响应[J]. 北京林业大学学报, 2018, 40(12): 18-26. doi: 10.13332/j.1000-1522.20180333
引用本文: 陈彬杭, 温晓示, 张树斌, 柴世品, 孙晗, 王襄平. 吉林北部山区长白落叶松林径向生长对气候干暖化的响应[J]. 北京林业大学学报, 2018, 40(12): 18-26. doi: 10.13332/j.1000-1522.20180333
Chen Binhang, Wen Xiaoshi, Zhang Shubin, Chai Shipin, Sun Han, Wang Xiangping. Growth responses of Larix olgensis forests to climatic drying-warming trend in the northern mountainous region of Jilin Province, northeastern China[J]. Journal of Beijing Forestry University, 2018, 40(12): 18-26. doi: 10.13332/j.1000-1522.20180333
Citation: Chen Binhang, Wen Xiaoshi, Zhang Shubin, Chai Shipin, Sun Han, Wang Xiangping. Growth responses of Larix olgensis forests to climatic drying-warming trend in the northern mountainous region of Jilin Province, northeastern China[J]. Journal of Beijing Forestry University, 2018, 40(12): 18-26. doi: 10.13332/j.1000-1522.20180333

吉林北部山区长白落叶松林径向生长对气候干暖化的响应

doi: 10.13332/j.1000-1522.20180333
基金项目: 

国家自然科学基金项目 31870430

科技部科技基础性工作专项 2015FY210200-8

国家重点研发计划课题 2017YFC0503901

详细信息
    作者简介:

    陈彬杭。主要研究方向:森林生态。Email: 229731286@qq.com  地址:100083  北京市海淀区清华东路35号北京林业大学林学院

    通讯作者:

    王襄平,教授,博士生导师。主要研究方向:森林生态。Email: wangxiangping@bjfu.edu.cn  地址:同上

  • 中图分类号: S718.51+2; S791.22

Growth responses of Larix olgensis forests to climatic drying-warming trend in the northern mountainous region of Jilin Province, northeastern China

  • 摘要: 目的中国东北地区的气候干暖化对长白落叶松林构成了严重威胁,有待进一步探究林分因子(如林龄、密度、蓄积量等)如何与气候条件共同作用、影响树木径向生长对气候变化响应的相关机制。方法本研究基于吉林长白山北部余脉至老爷岭间的长白落叶松林的树轮宽度标准年表,分析了气候条件及林分因子对长白落叶松林的年轮—气候关系中的相对作用大小,以研究长白落叶松在气候变化下的响应机制。结果近30年来,研究区气候干暖化趋势显著,长白落叶松的生长明显受到气候变化的影响,气候条件和林分因子共同作用于年轮—气候关系,但前者的影响要大于后者。Palmer干旱指数对主成分第1轴的解释力为50.27%,表明水分条件的影响更为重要;林分因子中的林龄有显著作用,林龄大的林分更易受到温度升高的限制作用,其余林分因子则在本研究中没有表现出显著影响。结论在未来气候干暖化情境下,林龄较小的长白落叶松的生长会得到一定促进,但若林龄过大则会受到温度上升导致水分亏缺的限制作用,所以有必要采取经营措施以减小升温对林分的不利影响,以保证林分的持续发展。
  • 图  1  研究地1986—2015年气候变化趋势

    Figure  1.  Climatic changing trend of research area from 1986 to 2015

    图  2  年轮—气候关系与林龄的关系

    Figure  2.  Relationship between growth-climate and stand age

    图  3  轮宽指数与季节气候指标相关系数的主成分分析

    Figure  3.  Principal components analysis of ring width index and seasonal climatic factors

    表  1  样地基本信息表

    Table  1.   Basic information of each sample plot

    样地号
    Sample plot
    纬度
    Latitude
    经度
    Longitude
    林龄
    Stand age/a
    立地条件
    Site
    condition
    林分密度/(株·hm-2)
    Stand density/
    (tree·ha-1)
    蓄积量/(m3·hm-2)
    Volume/
    (m3·ha-1)
    平均树高
    Mean tree
    height/m
    平均胸径
    Mean DBH/
    cm
    DPC142°51′N128°04′E31中Medium1 150238.5617.9716.93
    DPC242°51′N128°05′E41中Medium1 580379.5621.5516.19
    DPC342°50′N128°01′E24中Medium1 140204.7118.6915.35
    DPC442°50′N128°00′E22好Good1 960197.5914.6812.75
    DPC542°52′N128°01′E37中Medium820248.5419.3120.01
    JC0143°19′N130°34′E59差Bad340182.7120.2626.39
    JC0243°19′N130°34′E53差Bad430184.0219.5924.12
    JC0343°19′N130°33′E58差Bad1 820160.3012.1612.20
    JC0443°22′N130°35′E56中Medium1 480369.7116.4218.47
    LD0143°26′N130°50′E53好Good810427.9520.8726.04
    LD0243°26′N130°50′E57好Good1 600309.2815.5416.44
    LD0343°28′N130°39′E56好Good1 480172.9511.7114.26
    LD0443°28′N130°40′E72中Medium1 500231.6612.3215.63
    LD0543°31′N130°41′E29中Medium1 040285.3519.3319.27
    LD0643°31′N130°41′E28差Bad990223.4617.4818.06
    LD0743°31′N130°41′E28好Good1 090223.9716.9817.71
    LD0843°31′N130°41′E28差Bad1 020248.1618.3018.59
    注: DPC代表大蒲柴河; JC代表金苍; LD代表六道。下同。Notes: DPC is for Dapuchaihe; JC is for Jincang; LD is for Liudao. The same below.
    下载: 导出CSV

    表  2  样地标准年表统计参数

    Table  2.   Summary statistics for standard chronology of each plot

    样地
    Sample plot
    共同年区间
    Common
    interval time
    span of year
    平均敏
    感度
    Mean
    sensitivity
    标准差
    Standard
    deviation
    R1样本间平
    均相关系数
    Mean
    correlations
    among all
    radii
    R2不同树木
    之间平均相关系数
    Mean
    correlations
    between trees
    R3同一树木不同样
    本之间平均相关系数
    Mean correlations
    between trees and
    within trees
    信噪比
    Signal-to-noise
    ratio
    总体代表性
    Express
    population
    signal
    第1特征
    向量解释量
    PCA1/%
    DPC11990—20160.2030.2180.4470.6550.42221.1950.9550.408
    DPC21984—20160.2120.1880.3750.6390.36845.3420.9780.384
    DPC31997—20160.2150.1850.5120.5500.50621.6070.9560.449
    DPC42000—20160.2200.2180.3690.5720.3659.8840.9080.270
    DPC51992—20160.2480.2600.4490.6510.44627.5390.9650.388
    JC011969—20160.1780.2030.1800.5730.1727.1550.8770.511
    JC021974—20160.1810.2080.4440.5990.44015.9760.9410.487
    JC031969—20160.1830.1840.3670.6110.36326.0760.9630.548
    JC041979—20160.1900.1990.5050.7600.50112.7330.9270.288
    LD011979—20160.1030.1570.2170.5050.21327.8370.9650.395
    LD021966—20160.1350.2810.5010.7210.49413.5160.9310.429
    LD031968—20160.1360.1610.2980.4590.29516.9850.9440.419
    LD041961—20160.1500.1590.3400.6180.3376.4750.8660.527
    LD051994—20160.1550.1490.3260.5120.32243.1990.9770.491
    LD061993—20160.1560.2010.3510.5940.34352.4780.9810.520
    LD071993—20160.1620.3280.4430.5240.43555.6520.9820.531
    LD081993—20160.1690.1660.4930.6960.49056.1350.9820.528
    下载: 导出CSV

    表  3  PCA第1、2轴载荷

    Table  3.   Loading of seasonal climatic factors on PCA 1 and PCA 2

    季节气候指标
    Seasonal climatic factor
    PCA 1PCA 2
    P_Spr_T0.369
    P_Aut_T0.287-0.358
    P_Win_T-0.472
    C_Spr_T0.2720.216
    C_Sum_T0.275
    C_Aut _T-0.161-0.326
    P_Aut_T-0.375
    P_Spr_PDSI-0.165-0.331
    P_Sum_PDSI-0.364
    P_Aut_PDSI-0.356
    P_Win_PDSI-0.2290.241
    C_Spr_PDSI-0.382-0.106
    C_Sum_PDSI-0.381
    C_Aut_PDSI-0.3050.152
    注:P表示上年; C表示当年。Spr、Sum、Aut、Win分别表示春、夏、秋、冬季。下同。“—”表示系数绝对值过小。Notes: P means last year; C means current year. Spr, Sum, Aut, Win represent spring, summer, autumn, winter, respectively. The same below. “—” means absolute value of coefficients is too small.
    下载: 导出CSV

    表  4  样地气候条件及林分因子对PCA前两轴得分的解释力(R2)

    Table  4.   Explanatory power of climatic factors and stand conditions for the scores of PCA 1 and 2 axis (R2)

    变量VariablePCA 1PCA 2
    PDSI Palmer drought severity index-0.503**-0.006
    潜在蒸发量Potential evapotranspiration0.301*-0.024
    林龄Stand age-0.0650.163
    林分密度Stand density0.0080.004
    蓄积量Volume0.0140.164
    立地条件Stand factor0.1050.167
    注:*相关性较显著(P<0.05),**表示相关性显著(P<0.01),***表示相关性极显著(P<0.01)。下同。Notes: - means negative correlation. * represents significant correlation at P<0.05 level, ** represents significant correlation at P<0.01 level, *** represents significant correlation at P<0.001 level. The same below.
    下载: 导出CSV

    表  5  气候指标及林分条件对PCA前两轴得分的逐步回归分析方差分析

    Table  5.   Stepwise regression analysis of the scores of PCA 1 and PCA 2 with climatic factors and stand conditions

    项目ItemdfSS/%Sum SqMean SqFPr(>F)
    PCA 1
    PDSI Palmerdrought severity index150.2752.6352.6315.570.00***
    潜在蒸发量Potential evapotranspiration15.365.615.611.870.19
    林分密度Stand density17.187.517.512.510.14
    残差Residuals1337.1938.943.00
    PCA 2
    PDSI Palmerdrought severity index10.560.250.250.120.74
    林龄Stand age126.4812.0712.075.670.03*
    蓄积量Volume112.225.575.572.610.13
    残差Residual1360.7427.702.13
    下载: 导出CSV

    表  6  变量顺序不同的方差分析

    Table  6.   Appendix 1 variance analysis of variables in different order

    项目ItemdfSS/%Sum SqMean SqFPr(>F)
    PCA 1
    潜在蒸发量Potential evapotranspiration130.0531.4631.4610.500.01*
    PDSI Palmer drought severity index125.5826.7826.788.940.01*
    林分密度Stand density17.187.517.512.510.14
    残差Residual1337.1938.943.00
    PCA 2
    PDSI Palmer drought severity index10.560.250.250.120.74
    蓄积量Volume116.347.457.453.500.08
    林龄Stand density122.3610.1910.194.780.05
    残差Residual1360.7427.702.13
    下载: 导出CSV
  • [1] Hu F, Lee B, Kaufman D, et al. Response of tundra ecosystem in southwestern Alaska to Younger-Dryas climatic oscillation[J]. Global Change Biology, 2002, 8(11):1156 -1163. doi:  10.1046/j.1365-2486.2002.00550.x
    [2] Serreze M, Walsh J, Chapin Iii F, et al.Observational evidence of recent change in the northern high-latitude environment[J]. Climatic Change, 2000, 46(1-2): 159-207. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=9c0002f91192c93402544895933796d8
    [3] Kirschbaum M U F. Climate change impacts on forests[J]. Proceedings of the National Academy of Sciences of the United States of America, 2007, 104(50):19697-19702. doi:  10.1073/pnas.0701424104
    [4] 孙凤华, 任国玉, 赵春雨, 等.中国东北地区及不同典型下垫面的气温异常变化分析[J].地理科学, 2005, 25(2):167-171. doi:  10.3969/j.issn.1000-0690.2005.02.006

    Sun F H, Ren G Y, Zhao C Y, et al. An analysis of temperature abnormal change in northeast China and type underlying surface[J]. Scientia Geographica Sinica, 2005, 25(2):167-171. doi:  10.3969/j.issn.1000-0690.2005.02.006
    [5] 赵俊芳, 延晓冬, 贾根锁.东北森林净第一性生产力与碳收支对气候变化的响应[J].生态学报, 2008, 28(1):92-102. doi:  10.3321/j.issn:1000-0933.2008.01.011

    Zhao J F, Yan X D, Jia G S, et al. Simulating the responses of forest net primary productivity and carbon budget to climate change in Northeast China[J]. Acta Ecologica Sinica, 2008, 28(1):92-102. doi:  10.3321/j.issn:1000-0933.2008.01.011
    [6] 王宗明, 国志兴, 宋开山, 等.中国东北地区植被NDVI对气候变化的响应[J].生态学杂志, 2009, 28(6):1041-1048. http://d.old.wanfangdata.com.cn/Periodical/stxzz200906007

    Wang Z M, Guo Z X, Song K S, et al. Responses of vegetation NDVI in northeast China to climate change[J]. Chinese Journal of Ecology, 2009, 28(6):1041-1048. http://d.old.wanfangdata.com.cn/Periodical/stxzz200906007
    [7] 李腾, 何兴元, 陈振举.东北南部蒙古栎径向生长对气候变化的响应:以千山为例[J].应用生态学报, 2014, 25(7):1841-1848. http://d.old.wanfangdata.com.cn/Periodical/yystxb201407001

    Li T, He X Y, Chen Z J, et al. Tree-ring growth responses of Mongolian oak (Quercus mongolica) to climate change in southern Northeast: a case study in Qianshan Mountains[J]. Chinese Journal of Applied Ecology, 2014, 25(7):1841-1848. http://d.old.wanfangdata.com.cn/Periodical/yystxb201407001
    [8] 常永兴, 陈振举, 张先亮, 等.气候变暖下大兴安岭落叶松径向生长对温度的响应[J].植物生态学报, 2017, 41(3):279-289. http://d.old.wanfangdata.com.cn/Periodical/zwstxb201703001

    Chang Y X, Chen Z J, Zhang X L, et al. Responses of radial growth to temperature in Larix gmelinii of the Da Hinggan Ling under climate warming[J]. Acta Phytoecologica Sinica, 2017, 41(3):279-289. http://d.old.wanfangdata.com.cn/Periodical/zwstxb201703001
    [9] Andreu L, Gutierrez E, Macias M, et al. Climate increases regional tree- growth variability in Iberian pine forests[J]. Global Change Biology, 2007, 13:807-815. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=67752066cce37071f93ad3ce5c195f99
    [10] Sano M, Furuta F, Sweda T. Tree-ring-width chronology of Larix gmelinii, as an indicator of changes in early summer temperature in east-central Kamchatka[J]. Journal of Forest Research, 2009, 14(3):147-154. doi:  10.1007/s10310-009-0123-y
    [11] Wu X, Liu H, Guo D, et al. Growth decline linked to warming-induced water limitation in hemi-boreal forests[J/OL]. PLoS One, 2012, 7(8): e42619[2018-10-25]. https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0042619.
    [12] Bryukhanova M, Fonti P, Kirdyanov A, et al. Growth of trees on permafrost: habitat driven response to climate[C]//EGU General Assembly Conference. Vienna: EGU General Assembly Confer-ence Abstracts, 2013: 469-473.
    [13] 于大炮, 王顺忠, 唐立娜, 等.长白山北坡落叶松年轮年表及其与气候变化的关系[J].应用生态学报, 2005, 16(1):14-20. doi:  10.3321/j.issn:1001-9332.2005.01.003

    Yu D P, Wang S Z, Tang L N, et al. Relationship between tree-ring chronology of Larix olgensis in Changbai Mountains and the climate change[J]. Chinese Journal of Applied Ecology, 2005, 16(1):14-20. doi:  10.3321/j.issn:1001-9332.2005.01.003
    [14] 张先亮, 崔明星, 马艳军, 等.大兴安岭库都尔地区兴安落叶松年轮宽度年表及其与气候变化的关系[J].应用生态学报, 2010, 21(10):2501-2507. http://d.old.wanfangdata.com.cn/Periodical/yystxb201010008

    Zhang X L, Cui M X, Ma Y J, et al.Larix gmelinii tree-ring width chronology and its responses to climate change in Kuduer, Great Xing'an Mountains[J]. Chinese Journal of Applied Ecology, 2010, 21(10):2501-2507. http://d.old.wanfangdata.com.cn/Periodical/yystxb201010008
    [15] 孙毓, 王丽丽, 陈津.长白落叶松树轮生长对气候变化的响应以及对夏季气温的重建[J].地球环境学报, 2012, 3(3):889-899. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QKV20122015120800030018

    Sun Y, Wang L L, Chen J. Response of tree growth to climate change and reconstruction of summer temperature based on Korea larch[J]. Journal of Earth Environment, 2012, 3(3):889-899. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QKV20122015120800030018
    [16] 陈力, 尹云鹤, 赵东升, 等.长白山不同海拔树木生长对气候变化的响应差异[J].生态学报, 2014, 34(6):1568-1574. http://d.old.wanfangdata.com.cn/Periodical/stxb201406026

    Chen L, Yin Y H, Zhao D S, et al. Climate response of tree growth along an altitudinal gradient in the Changbai Mountains, Northeast China[J]. Acta Ecologica Sinica, 2014, 34(6):1568-1574. http://d.old.wanfangdata.com.cn/Periodical/stxb201406026
    [17] 于健, 徐倩倩, 刘文慧, 等.长白山东坡不同海拔长白落叶松径向生长对气候变化的响应[J].植物生态学报, 2016, 40(1):24-35. http://d.old.wanfangdata.com.cn/Periodical/zwstxb201601003

    Yu J, Xu Q Q, Liu W H, et al. Response of radial growth to climate change for Larix olgensis along an altitudinal gradient on the eastern slope of Changbai Mountain, Northeast China[J]. Acta Phytoecologica Sinica, 2016, 40(1):24-35. http://d.old.wanfangdata.com.cn/Periodical/zwstxb201601003
    [18] 吴征镒.中国植被[M].北京:科学出版社, 1995.

    Wu Z Y. Vegetation of China[M].Beijing:Science Press, 1995.
    [19] 梁鹏鸿, 王襄平, 吴玉莲, 等.黑龙江胜山保护区阔叶红松林不同演替阶段径向生长与气候变化的关系[J].植物生态学报, 2016, 40(5):425-435. http://d.old.wanfangdata.com.cn/Periodical/zwstxb201605001

    Liang P H, Wang X P, Wu Y L, et al. Growth responses of broad-leaf and Korean pine mixed forests at different successional stages to climate change in the Shengshan Nature Reserve of Heilongjiang Province, China[J]. Acta Phytoecologica Sinica, 2016, 40(5):425-435. http://d.old.wanfangdata.com.cn/Periodical/zwstxb201605001
    [20] 吴正方.东北阔叶红松林分布区生态气候适宜性及全球气候变化影响评价[J].应用生态学报, 2003, 14(5):771-775. doi:  10.3321/j.issn:1001-9332.2003.05.028

    Wu Z F. Assessment of eco-climatic suitability and climate change impacts of/on broad-leaved Korean pine forest in Northeast China[J]. Chinese Journal of Applied Ecology, 2003, 14(5):771-775. doi:  10.3321/j.issn:1001-9332.2003.05.028
    [21] Brooks J R, Flanagan L B, Ehleringer J R. Responses of boreal conifers to climate fluctuations: indications from tree-ring widths and carbon isotope analyses[J]. Revue Canadienne De Recherche Forestière, 1998, 28(4):524-533. doi:  10.1139/x98-018
    [22] 常锦峰, 王襄平, 张新平, 等.大兴安岭北部大白山高山林线动态与气候变化的关系[J].山地学报, 2009, 27(6):703-711. doi:  10.3969/j.issn.1008-2786.2009.06.010

    Chang J F, Wang X P, Zhang X P, et al. Alpine timber line dynamics in relation to climatic variability in the Northern Daxing'an Mountains[J]. Mountain Research, 2009, 27(6):703-711. doi:  10.3969/j.issn.1008-2786.2009.06.010
    [23] Wang X, Zhang Y, Mcrae D J. Spatial and age-dependent tree-ring growth responses of Larix gmelinii, to climate in northeastern China[J]. Trees Structure & Function, 2009, 23(4):875-885. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=e6a416bc8c67c4f70c8ff2ddc4e2513e
    [24] Hughes M K, Schweingruber F H, Cartwright D, et al. July-August temperature at Edinburgh between 1721 and 1975 from tree-ring density and width data[J]. Nature, 1984, 308:341-344. doi:  10.1038/308341a0
    [25] Fritts H C. Reconstructing large-scale climatic patterns from tree-ring data: t diagnostic analysis[J]. University of Arizona Press, 1991, 16(3):245-248.
    [26] 邵雪梅, 吴祥定.华山树木年轮年表的建立[J].地理学报, 1994, 49(2):174-181. doi:  10.3321/j.issn:0375-5444.1994.02.009

    Shao X M, Wu X D. Tree-ring chronologies for Pinus armandi Franch from Huashan, China[J]. Acta Geographica Sinica, 1994, 49(2):174-181. doi:  10.3321/j.issn:0375-5444.1994.02.009
    [27] 吴祥定.树木年轮与气候变化[M].北京:气象出版社, 1990.

    Wu X D. Tree ring and climatic change[M].Beijing:Meteorology Press, 1990.
    [28] 李建东.吉林植被[M].长春:吉林科学技术出版社, 2001.

    Li J D. Vegetation of Jilin[M].Changchun:Jilin Science and Technology Press, 2001.
    [29] Law B E, Sun O J, Campbell J, et al. Changes in carbon storage and fluxes in a chronosequence of ponderosa pine[J]. Global Change Biology, 2010, 9(4):510-524. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1046/j.1365-2486.2003.00624.x
    [30] Holmes R L. Computer-assisted quality control in tree-ring dating and measurement[J]. Tree-Ring Bull, 1983, 43(3):69-75.
    [31] Cook E, Holmes R. Users manual for program ARSTAN: laboratory of tree-ring research[M]. Tucson: University of Arizona, 1986.
    [32] 李江风, 袁玉江, 由希尧.树木年轮水文学研究与应用[M].北京:科学出版社, 2000.

    Li J F, Yuan Y J, You X Y. Dendrohydrology research and application[M].Beijing: Science Press, 2000.
    [33] Wang X P, Tang Z Y, Fang J Y. Climatic control on forests and tree species distribution in the forest region of Northeast China[J]. Journal of Integrative Plant Biology, 2006, 48(7):778-789. doi:  10.1111/jipb.2006.48.issue-7
    [34] O'Brien E. Biological relativity to water-energy dynamics[J]. Journal of Biogeography, 2010, 33(11):1868-1888. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=d214059adb461f71b4242dfd3a32fd94
    [35] Yamaguchi D, Cook E, Kairiukstis L. Methods of dendrochronology. applications in environmental science[J]. Arctic & Alpine Research, 1990, 23(1):120.
    [36] Yu D, Wang G G, Dai L, et al. Dendroclimatic analysis of Betula ermanii, forests at their upper limit of distribution in Changbai Mountain, Northeast China[J]. Forest Ecology & Management, 2007, 240(1-3):105-113. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=39e7ac05bd44d8ad307bdff1f7128897
    [37] 曾令兵, 王襄平, 常锦峰, 等.祁连山中段青海云杉高山林线交错区树轮宽度与气候变化的关系[J].北京林业大学学报, 2012, 34(5):50-56. http://j.bjfu.edu.cn/article/id/9799

    Zeng L B, Wang X P, Chang J F, et al. Alpine timberline ecotone tree growth in relation to climatic variability for Picea crassifolia forests in the middle Qilian Mountains, northwestern China[J]. Journal of Beijing Forestry University, 2012, 34(5):50-56. http://j.bjfu.edu.cn/article/id/9799
    [38] Allen C D, Breshears D D, Mcdowell N G. On underestimation of global vulnerability to tree mortality and forest die-off from hotter drought in the Anthropocene[J]. Ecosphere, 2016, 6(8):1-55. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1890/ES15-00203.1
    [39] Peterson D W, Peterson D L, Ettl G J. Growth responses of subalpine fir to climatic variability in the Pacific Northwest[J]. Canadian Journal of Forest Research, 2002, 32(9):1503-1517. doi:  10.1139/x02-072
    [40] Carrer M, Urbinati C. Age-dependent tree ring growth responses to climate of Larix decidua and Pinus cembra in the Italian Alps[J]. Ecology, 2004, 85(3):730-740. doi:  10.1890/02-0478
    [41] Martínez-Vilalta J, López B C, Loepfe L, et al. Stand- and tree-level determinants of the drought response of Scots pine radial growth[J]. Oecologia, 2012, 168(3):877-888. doi:  10.1007/s00442-011-2132-8
    [42] Rais A, Kuilen J W G V D, Pretzsch H. Growth reaction patterns of tree height, diameter, and volume of Douglas-fir (Pseudotsuga menziesii, [Mirb.] Franco) under acute drought stress in Southern Germany[J]. European Journal of Forest Research, 2014, 133(6):1043-1056. doi:  10.1007/s10342-014-0821-7
  • [1] 牛亦龙, 董利虎, 李凤日.  基于广义代数差分法的长白落叶松人工林地位指数模型 . 北京林业大学学报, 2020, 42(2): 9-18. doi: 10.12171/j.1000-1522.20190036
    [2] 贺梦莹, 董利虎, 李凤日.  长白落叶松−水曲柳混交林冠幅预测模型 . 北京林业大学学报, 2020, 42(7): 23-32. doi: 10.12171/j.1000-1522.20190250
    [3] 欧强新, 雷相东, 沈琛琛, 宋国涛.  基于随机森林算法的落叶松−云冷杉混交林单木胸径生长预测 . 北京林业大学学报, 2019, 41(9): 9-19. doi: 10.13332/j.1000-1522.20180266
    [4] 徐奇刚, 雷相东, 国红, 李海奎, 李玉堂.  基于多层感知机的长白落叶松人工林林分生物量模型 . 北京林业大学学报, 2019, 41(5): 97-107. doi: 10.13332/j.1000-1522.20190035
    [5] 白东雪, 刘强, 董利虎, 李凤日.  长白落叶松人工林有效冠高的确定及其影响因子 . 北京林业大学学报, 2019, 41(5): 76-87. doi: 10.13332/j.1000-1522.20190016
    [6] 罗梅, 陈绍志.  不同龄组长白落叶松种内及种间竞争研究 . 北京林业大学学报, 2018, 40(9): 33-44. doi: 10.13332/j.1000-1522.20180126
    [7] 沈剑波, 雷相东, 雷渊才, 李玉堂.  长白落叶松人工林地位指数及立地形的比较研究 . 北京林业大学学报, 2018, 40(6): 1-8. doi: 10.13332/j.1000-1522.20170400
    [8] 宋金凤, 李金博, 曹楷, 桑英, 崔晓阳.  草酸和柠檬酸提高长白落叶松对Pb胁迫的适应性 . 北京林业大学学报, 2017, 39(11): 18-27. doi: 10.13332/j.1000-1522.20170146
    [9] 姜礅, 孟昭军, 严善春.  用茉莉酸甲酯局部喷施长白落叶松苗对落叶松毛虫体内防御酶的影响 . 北京林业大学学报, 2017, 39(2): 58-63. doi: 10.13332/j.1000-1522.20160291
    [10] 姜礅, 王杰, 姜虹, 张文一, 孟昭军, 严善春.  茉莉酸甲酯局部诱导长白落叶松对舞毒蛾生长发育的影响 . 北京林业大学学报, 2016, 38(6): 67-71. doi: 10.13332/j.1000-1522.20150493
    [11] 孙志虎, 王秀琴, 陈祥伟.  不同抚育间伐强度对落叶松人工林生态系统碳储量影响 . 北京林业大学学报, 2016, 38(12): 1-13. doi: 10.13332/j.1000-1522.20160016
    [12] 张素芳, 张磊, 赵佳丽, 张莉, 张含国.  长白落叶松小RNA测序和其靶基因预测 . 北京林业大学学报, 2016, 38(12): 64-72. doi: 10.13332/j.1000-1522.20150404
    [13] 姚丹丹, 雷相东, 张则路.  基于贝叶斯法的长白落叶松林分优势高生长模型研究 . 北京林业大学学报, 2015, 37(3): 94-100. doi: 10.13332/j.1000-1522.20140221
    [14] 刘艳红, 马炜.  长白落叶松人工林可燃物碳储量分布及燃烧性 . 北京林业大学学报, 2013, 35(3): 32-38.
    [15]
    孙志虎, 毕永娟, 牟长城, 蔡体久
    基于FORECAST模型的长白落叶松人工林经营措施对长期生产力的影响 . 北京林业大学学报, 2012, 34(6): 1-6.
    [16] 许晨璐, 孙晓梅, 张守攻.  日本落叶松与长白落叶松及其杂种光合特性比较 . 北京林业大学学报, 2012, 34(4): 62-66.
    [17] 康瑶瑶, 刘勇, 马履一, 李国雷, 祝燕, 马跃.  施肥对长白落叶松苗木养分库氮磷吸收及利用的影响 . 北京林业大学学报, 2011, 33(2): 31-36.
    [18] 孙慧彦, 刘勇, 马履一, 贾忠奎, 康瑶瑶, 金虎范, 祝燕, 侯炳柱, 尹凤君.  长白落叶松苗木质量与造林效果关系的比较 . 北京林业大学学报, 2009, 31(6): 176-180.
    [19] 孙志虎, 牟长城, 张彦东.  用地统计学方法估算长白落叶松人工林凋落物现存量 . 北京林业大学学报, 2008, 30(4): 59-64.
    [20] 崔彬彬, 李贤军, 宗世祥, 赵俊卉, 肖化顺, 陈伟, 刘志军, 王志玲, 曹伟, 黄心渊, 张煜星, 周国模, 李国平, 江泽慧, 雷相东, 刘智, 施婷婷, 张展羽, 于寒颖, 周志强, 杜官本, 徐剑琦, 程金新, 雷霆, 程丽莉, 曹金珍, 关德新, 刘童燕, 张贵, 苏里坦, 吴家森, 骆有庆, 王正, 丁立建, 王正, 张则路, 张彩虹, 王海, 杨谦, 张璧光, 苏淑钗, 李云, 张璧光, 郭广猛, 郝雨, 黄群策, 雷洪, 李云, 张国华, 刘彤, 金晓洁], 吴家兵, 黄晓丽, 贺宏奎, 王勇, 张书香, 张慧东, 常亮, 秦岭, 方群, 秦广雍, 张佳蕊, 许志春, 张大红, 陈晓光, 宋南, 刘大鹏, 姜培坤, 李文军, 周晓燕, 李延军, 高黎, 刘海龙, 蔡学理, 陈燕, 姜静, 姜金仲, 张弥, 冯慧, 苏晓华, 于兴华, 张金桐, 刘建立, 王安志, 张冰玉, 尹伟伦, 陈绪和, 周梅, 王谦, 朱彩霞, 成小芳, 王德国, 陈建伟3, 聂立水, 亢新刚, 张连生, 张勤, 冯大领, 金昌杰, 梁树军, 崔国发, 韩士杰, 胡君艳, 姚国龙.  长白落叶松等几个树种冠幅预测模型的研究 . 北京林业大学学报, 2006, 28(6): 75-79.
  • 加载中
图(3) / 表 (6)
计量
  • 文章访问数:  783
  • HTML全文浏览量:  1006
  • PDF下载量:  20
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-10-25
  • 修回日期:  2018-11-27
  • 刊出日期:  2018-12-01

吉林北部山区长白落叶松林径向生长对气候干暖化的响应

doi: 10.13332/j.1000-1522.20180333
    基金项目:

    国家自然科学基金项目 31870430

    科技部科技基础性工作专项 2015FY210200-8

    国家重点研发计划课题 2017YFC0503901

    作者简介:

    陈彬杭。主要研究方向:森林生态。Email: 229731286@qq.com  地址:100083  北京市海淀区清华东路35号北京林业大学林学院

    通讯作者: 王襄平,教授,博士生导师。主要研究方向:森林生态。Email: wangxiangping@bjfu.edu.cn  地址:同上
  • 中图分类号: S718.51+2; S791.22

摘要: 目的中国东北地区的气候干暖化对长白落叶松林构成了严重威胁,有待进一步探究林分因子(如林龄、密度、蓄积量等)如何与气候条件共同作用、影响树木径向生长对气候变化响应的相关机制。方法本研究基于吉林长白山北部余脉至老爷岭间的长白落叶松林的树轮宽度标准年表,分析了气候条件及林分因子对长白落叶松林的年轮—气候关系中的相对作用大小,以研究长白落叶松在气候变化下的响应机制。结果近30年来,研究区气候干暖化趋势显著,长白落叶松的生长明显受到气候变化的影响,气候条件和林分因子共同作用于年轮—气候关系,但前者的影响要大于后者。Palmer干旱指数对主成分第1轴的解释力为50.27%,表明水分条件的影响更为重要;林分因子中的林龄有显著作用,林龄大的林分更易受到温度升高的限制作用,其余林分因子则在本研究中没有表现出显著影响。结论在未来气候干暖化情境下,林龄较小的长白落叶松的生长会得到一定促进,但若林龄过大则会受到温度上升导致水分亏缺的限制作用,所以有必要采取经营措施以减小升温对林分的不利影响,以保证林分的持续发展。

English Abstract

陈彬杭, 温晓示, 张树斌, 柴世品, 孙晗, 王襄平. 吉林北部山区长白落叶松林径向生长对气候干暖化的响应[J]. 北京林业大学学报, 2018, 40(12): 18-26. doi: 10.13332/j.1000-1522.20180333
引用本文: 陈彬杭, 温晓示, 张树斌, 柴世品, 孙晗, 王襄平. 吉林北部山区长白落叶松林径向生长对气候干暖化的响应[J]. 北京林业大学学报, 2018, 40(12): 18-26. doi: 10.13332/j.1000-1522.20180333
Chen Binhang, Wen Xiaoshi, Zhang Shubin, Chai Shipin, Sun Han, Wang Xiangping. Growth responses of Larix olgensis forests to climatic drying-warming trend in the northern mountainous region of Jilin Province, northeastern China[J]. Journal of Beijing Forestry University, 2018, 40(12): 18-26. doi: 10.13332/j.1000-1522.20180333
Citation: Chen Binhang, Wen Xiaoshi, Zhang Shubin, Chai Shipin, Sun Han, Wang Xiangping. Growth responses of Larix olgensis forests to climatic drying-warming trend in the northern mountainous region of Jilin Province, northeastern China[J]. Journal of Beijing Forestry University, 2018, 40(12): 18-26. doi: 10.13332/j.1000-1522.20180333
  • 气候变化已经成为21世纪人类共同面临的重大环境问题,大气中CO2等温室气体浓度的增加,致使全球气候系统正在发生着有史以来从未有过的急剧变化,主要表现在全球变暖、降水格局变化、海平面上升等方面[1]。全球变暖的影响在北半球中高纬度的大陆内部表现的尤为明显[2]。森林作为陆地生态系统的重要组成部分,一直以来都是气候变化研究中的重点对象[3]。东北地区是我国纬度最高的区域,该区域内的森林对气候变化的响应受到全球气候变暖的影响,还具有明显的地域性特点[4]。众多研究表明,温度、降水对该地区的森林都有一定的影响[5-7],其中高温造成的负面影响尤为显著[8]

    长白落叶松(Larix olgensis)是北半球高纬度地区分布较为广泛的针叶树种之一,在海拔500~1 950 m均有分布并常形成大面积纯林。国外已有不少针对落叶松相关种属的生长与气候变化关系的研究在与我国临近的西伯利亚地区开展[9-12]。我国已有较多关于东北地区长白落叶松对气候变化响应的研究[13-15],但大多集中于长白山一带且多是针对高海拔或林线处林分[13, 16-17],在其他区域的研究还相对较少。吉林北部山区接近长白落叶松在我国的分布北界[18],在该区域开展研究对理解长白落叶松林对气候变化的响应有着重要的意义[19]。同时,该区域气候干暖化趋势明显,对树木生长已经产生了显著影响[20],急需开展相关研究为该地区森林的经营和保护提供理论依据。

    树木生长对气候变化的响应受到树种、气候及林分因子等多方面的影响。如Brooks等[21]提出气候变化对不同树种的影响有差异,温度上升可能严重地影响某一树种的生长,却对另一树种影响不大。针对长白落叶松、兴安落叶松(L. gmelinii)进行的研究发现,不同海拔(温度、降水梯度)的落叶松对气候变化的响应存在明显差异[13, 22]。除树种、气候之外,还存在其他因素的影响。例如不同树龄的兴安落叶松的生长对气候变化的响应存在明显差异,老龄林更易受到水分胁迫的影响[23]。因此,气候条件、林分因子对森林与气候变化之间的响应关系都存在影响,但鲜有研究同时探究二者共同的影响及其各自的相对作用大小。

    目前有关森林对气候变化响应的研究,多采用植被覆盖指数(Normalized difference vegetation index, NDVI)、碳收支、生产力等指标作为反映[5-6]。而用树木年轮(简称树轮)来研究生态系统对气候变化的响应,是全球气候变化研究的热点[24-26],相对其他方法,具有定年准确、连续性强、分辨率高等优点[27]。研究树木径向生长与气候因子的关系有助于理解一个地区森林生产力的生态限制因素,以及气候变化条件下森林生态系统的变化机理[8-9]

    本文以吉林北部山区不同林龄的长白落叶松天然林为研究对象,通过研究树木径向生长与历史气候变化的关系,探讨了气候条件、林分因子(林龄、林分密度、蓄积量、立地条件等)对年轮—气候关系的影响,以明确影响长白落叶松林对气候变化响应方式的主要影响因素。这对于在未来全球气候变暖的大背景下为该地区的长白落叶松林提供合理的保护抚育措施具有参考价值。

    • 研究区为吉林省北部山区,位于长白山北部余脉至老爷岭之间(42°49′48″~43°31′12″ N、128°00′36″~130°41′24″ E),沿纬度梯度自南向北包括敦化市的大蒲柴河林场(DPC)、汪清县的金苍林场(JC)和六道林场(LD)等3个研究地点,海拔550~800 m。研究区内为典型的温带季风气候,具有冬季冷且长、夏季热而短且四季分明的特点,年均气温1.0~3.0 ℃,最冷月均温-18 ℃,年降水量690~750 mm。该地区地带性植被为阔叶红松(Pinus koraiensis)林(分布于海拔500~1 100 m),也分布着各种次生森林类型,以及长白落叶松等人工林[28]。研究区内的长白落叶松林混生有少量白桦(Betula platyphylla)、春榆(Ulmus davidiana)、蒙古栎(Quercus mongolica)等阔叶树种,同时还伴生有少量欧亚大陆针叶林中的寒温性树种,如鱼鳞云杉(Picea jezonensis var. microsperma)、红皮云杉(Picea koraiensis)等。土壤以棕黑色的壤土或砂壤土为主。

    • 野外调查于2016年8—9月份进行,在上述3个研究地点分别设置5、4和8个20 m×50 m的长白落叶松样地(表 1)。对样地中所有胸径大于3 cm的乔木个体挂牌并进行每木检尺,测量树高、冠幅等。各样地的林龄采取国内外相关研究中常用的方法进行估测[29],使用生长锥钻取样地中胸径最大的5株乔木的完整年轮样品(至髓心),以这5个个体的平均年龄作为该样地的林龄,这样可以防止出现由于个别个体树龄较大而导致对林分年龄的估计偏高。另外,蓄积量采用相关蓄积量方程计算得到,立地条件则参考土壤表层含水率并结合实地调查确定,划分为好、中、差3个等级。

      表 1  样地基本信息表

      Table 1.  Basic information of each sample plot

      样地号
      Sample plot
      纬度
      Latitude
      经度
      Longitude
      林龄
      Stand age/a
      立地条件
      Site
      condition
      林分密度/(株·hm-2)
      Stand density/
      (tree·ha-1)
      蓄积量/(m3·hm-2)
      Volume/
      (m3·ha-1)
      平均树高
      Mean tree
      height/m
      平均胸径
      Mean DBH/
      cm
      DPC142°51′N128°04′E31中Medium1 150238.5617.9716.93
      DPC242°51′N128°05′E41中Medium1 580379.5621.5516.19
      DPC342°50′N128°01′E24中Medium1 140204.7118.6915.35
      DPC442°50′N128°00′E22好Good1 960197.5914.6812.75
      DPC542°52′N128°01′E37中Medium820248.5419.3120.01
      JC0143°19′N130°34′E59差Bad340182.7120.2626.39
      JC0243°19′N130°34′E53差Bad430184.0219.5924.12
      JC0343°19′N130°33′E58差Bad1 820160.3012.1612.20
      JC0443°22′N130°35′E56中Medium1 480369.7116.4218.47
      LD0143°26′N130°50′E53好Good810427.9520.8726.04
      LD0243°26′N130°50′E57好Good1 600309.2815.5416.44
      LD0343°28′N130°39′E56好Good1 480172.9511.7114.26
      LD0443°28′N130°40′E72中Medium1 500231.6612.3215.63
      LD0543°31′N130°41′E29中Medium1 040285.3519.3319.27
      LD0643°31′N130°41′E28差Bad990223.4617.4818.06
      LD0743°31′N130°41′E28好Good1 090223.9716.9817.71
      LD0843°31′N130°41′E28差Bad1 020248.1618.3018.59
      注: DPC代表大蒲柴河; JC代表金苍; LD代表六道。下同。Notes: DPC is for Dapuchaihe; JC is for Jincang; LD is for Liudao. The same below.

      年轮采样选取了位于林冠层且健康无病害的长白落叶松30~40株,使用生长锥在胸径处(离地面1.3 m)沿东西、南北两个方向各取1个树轮样品。所采样品装入塑料吸管内保存。在实验室内按照标准处理程序进行固定、干燥并用砂纸打磨至年轮界限在显微镜下清晰可见。然后使用LINTAB 6树轮测量系统(Rinntech, Germany)以0.01 mm分辨率精确测量逐年轮宽序列。之后使用COFECHA程序[30]对同一样地中的所有年轮宽度序列进行相关性检验,剔除与主序列相关性差的序列,最后使用ARSTAN程序为剩余序列建立标准年表(表 2)[31]。基于此,得到了所有样地的轮宽指数(Ring width index, RWI)序列用于表征林分的逐年径向生长情况。由于各样地林龄不同,为了分析结果的可比较性,截取了所有样地的RWI序列共有的区间(1986—2015年)用于后续的数据分析。

      表 2  样地标准年表统计参数

      Table 2.  Summary statistics for standard chronology of each plot

      样地
      Sample plot
      共同年区间
      Common
      interval time
      span of year
      平均敏
      感度
      Mean
      sensitivity
      标准差
      Standard
      deviation
      R1样本间平
      均相关系数
      Mean
      correlations
      among all
      radii
      R2不同树木
      之间平均相关系数
      Mean
      correlations
      between trees
      R3同一树木不同样
      本之间平均相关系数
      Mean correlations
      between trees and
      within trees
      信噪比
      Signal-to-noise
      ratio
      总体代表性
      Express
      population
      signal
      第1特征
      向量解释量
      PCA1/%
      DPC11990—20160.2030.2180.4470.6550.42221.1950.9550.408
      DPC21984—20160.2120.1880.3750.6390.36845.3420.9780.384
      DPC31997—20160.2150.1850.5120.5500.50621.6070.9560.449
      DPC42000—20160.2200.2180.3690.5720.3659.8840.9080.270
      DPC51992—20160.2480.2600.4490.6510.44627.5390.9650.388
      JC011969—20160.1780.2030.1800.5730.1727.1550.8770.511
      JC021974—20160.1810.2080.4440.5990.44015.9760.9410.487
      JC031969—20160.1830.1840.3670.6110.36326.0760.9630.548
      JC041979—20160.1900.1990.5050.7600.50112.7330.9270.288
      LD011979—20160.1030.1570.2170.5050.21327.8370.9650.395
      LD021966—20160.1350.2810.5010.7210.49413.5160.9310.429
      LD031968—20160.1360.1610.2980.4590.29516.9850.9440.419
      LD041961—20160.1500.1590.3400.6180.3376.4750.8660.527
      LD051994—20160.1550.1490.3260.5120.32243.1990.9770.491
      LD061993—20160.1560.2010.3510.5940.34352.4780.9810.520
      LD071993—20160.1620.3280.4430.5240.43555.6520.9820.531
      LD081993—20160.1690.1660.4930.6960.49056.1350.9820.528
    • 本研究使用了2个来源的气候数据。用于分析树木径向生长与历史气候之间关系所使用的数据来自于东安格利亚大学的CRU TS 3.22高分辨率气候数据集,我们从中提取了1986—2015年间的月均温及帕尔默干旱指数(Palmer drought severity index, PDSI)。其中PDSI是一个基于降水量和温度的反映干燥程度的指标,能够很好反映水分的盈缺。为使分析结果更清晰,按照季节对一年的12个月进行了划分,其中上年12月至当年2月为冬季,3—5月为春季,6—8月为夏季,9—11月为秋季,并将对应月份的温度、PDSI取平均作为季节气候数据。CRU数据的空间分辨率为0.5°(研究区所在纬度,栅格边长约为41 km),由于金苍和六道两个研究地点的距离较近,位于同一栅格内,因此在分析时两地使用的气候数据相同。考虑到树木当年的生长与上年的气候之间存在“滞后效应”[32],因此在数据分析时除了使用当年春季到秋季的气候数据,还使用了上年四季的气候数据。

      但上述CRU的气候数据无法反映出不同样地间的气候条件差异。采用Wang等[33]利用国家气象台站数据建立的模型,基于样地经、纬度和海拔估算了各样地的年均温(Annual mean temperature, AMT)、年均降水量(Annual precipitation, AP)。然后使用估算的年均温、降水计算了各样地的潜在蒸散量(Potential evapotranspiration, PET)[33]。研究表明,能量和水分的有效性是影响生产力的关键气候要素[34],这里采用PET和PDSI分别反映能量和水分的有效性。由于无法获取计算PDSI需要的一些数据,因此采用从CRU数据集中提取的PDSI数据(取1986—2015年间PDSI的均值)。

    • 首先分析了研究区近30年的年均温、PDSI的变化趋势,然后对各样地的RWI序列与对应的各季节(上年春季到当年秋季)的平均温度、PDSI序列进行相关性分析。由于季节气候指标较多,因此RWI与季节气候指标的相关系数只能用以初步判断影响年轮—气候关系的因素。为能明确起主要作用的因素,对上述得到的相关系数矩阵进行主成分分析。为研究气候条件和林分因子的相对影响大小,计算了各样地的气候及林分因子对PCA前2轴得分的解释力(R2),同时建立多元线性模型并使用逐步回归筛选得到对年轮—气候关系影响最显著的变量。主成分分析和逐步回归分析均用R软件完成。

    • 本研究中不同地点在空间上有一定距离,也表现在气候上的差异(图 1)。大蒲柴河(DPC)温度较金苍(JC)和六道(LD)林场低,PDSI则较金苍和六道低。即相对金苍和六道,大蒲柴河的气候条件更为干冷,但近30年(1986—2015年)的年气候指标变化仍表现出相似的趋势。两地均以0.014 ℃/a的速率变暖,PDSI也以近似速率下降。结合上述分析的温度和PDSI数据变化趋势,不难发现研究区在近30年呈现明显的干暖化趋势。

      图  1  研究地1986—2015年气候变化趋势

      Figure 1.  Climatic changing trend of research area from 1986 to 2015

    • 建立标准年表,其统计参数如表 2所示,年表的平均敏感度介于0.103~0.248之间,标准差则介于0.149~0.281之间,样本总体代表性高达0.866~0.982,均超过了样本总体代表性(EPS)的阈值0.85[35]。这完全符合吴祥定[27]对于高质量年表的要求,适合进一步用以年轮与气候变化关系的分析。

      对标准年表轮宽指数与季节温度、PDSI进行的相关分析表明,落叶松径向生长明显受到气候变化的影响,与上年和当年多个季节的温度和PDSI都有显著的相关性。更为重要的是,不同林龄的林分对气候变化的响应存在明显差异(图 2),温度升高对林龄小的林分有促进生长的作用,而对林龄过大的林分则会产生抑制作用,且林龄越大抑制作用越明显。

      图  2  年轮—气候关系与林龄的关系

      Figure 2.  Relationship between growth-climate and stand age

      在对年轮—气候关系进行主成分分析后,发现PCA 1轴和2轴共解释了63.1%的变异,表明PCA前2轴能较好地反映年轮—气候关系中的主要信息,因此可用于进一步探讨不同样地的气候条件、林分因子等因素的影响。由PCA前2轴载荷系数(表 3)中可以看出,PCA 1的季节PDSI的载荷系数均为负数,且绝对值明显高于温度指标,这说明PCA 1主要反映径向生长对水分历史变化的响应,径向生长受到水分的限制作用,PCA 1得分为负时值越小,受到水分限制作用越显著。而PCA 2的季节温度的载荷系数也多为负数且绝对值明显高于PDSI,主要反映轮宽与温度指标的历史变化的相关性,PCA 2得分的负值越小表示受温度的限制作用越强。

      表 3  PCA第1、2轴载荷

      Table 3.  Loading of seasonal climatic factors on PCA 1 and PCA 2

      季节气候指标
      Seasonal climatic factor
      PCA 1PCA 2
      P_Spr_T0.369
      P_Aut_T0.287-0.358
      P_Win_T-0.472
      C_Spr_T0.2720.216
      C_Sum_T0.275
      C_Aut _T-0.161-0.326
      P_Aut_T-0.375
      P_Spr_PDSI-0.165-0.331
      P_Sum_PDSI-0.364
      P_Aut_PDSI-0.356
      P_Win_PDSI-0.2290.241
      C_Spr_PDSI-0.382-0.106
      C_Sum_PDSI-0.381
      C_Aut_PDSI-0.3050.152
      注:P表示上年; C表示当年。Spr、Sum、Aut、Win分别表示春、夏、秋、冬季。下同。“—”表示系数绝对值过小。Notes: P means last year; C means current year. Spr, Sum, Aut, Win represent spring, summer, autumn, winter, respectively. The same below. “—” means absolute value of coefficients is too small.
    • 从PCA分析结果(图 3)可以看到,由于气候条件的差异,能够区分出各个地点所对应的样地,DPC的样地居于图右侧,而LD的样地则位于左侧,JC的样地多数处于二者之间。结合各个样地的PET数据,PCA 1轴上越往右,对应样地的PET越大,PDSI的箭头也均朝向左边。这表明PCA1轴反映了样地气候条件差异对年轮—气候关系的影响(表 3)。通过计算样地气候条件、林分因子等对PCA轴的单因素解释力(表 4),发现PET和PDSI都与PCA 1得分显著相关,且两者作用方向相反,而林龄、林分密度等林分因子与PCA 1、PCA 2得分的相关性均不显著。

      图  3  轮宽指数与季节气候指标相关系数的主成分分析

      Figure 3.  Principal components analysis of ring width index and seasonal climatic factors

      表 4  样地气候条件及林分因子对PCA前两轴得分的解释力(R2)

      Table 4.  Explanatory power of climatic factors and stand conditions for the scores of PCA 1 and 2 axis (R2)

      变量VariablePCA 1PCA 2
      PDSI Palmer drought severity index-0.503**-0.006
      潜在蒸发量Potential evapotranspiration0.301*-0.024
      林龄Stand age-0.0650.163
      林分密度Stand density0.0080.004
      蓄积量Volume0.0140.164
      立地条件Stand factor0.1050.167
      注:*相关性较显著(P<0.05),**表示相关性显著(P<0.01),***表示相关性极显著(P<0.01)。下同。Notes: - means negative correlation. * represents significant correlation at P<0.05 level, ** represents significant correlation at P<0.01 level, *** represents significant correlation at P<0.001 level. The same below.

      通过逐步回归筛选相关气候因子和林分因子后的结果(表 5),可以看到PDSI对PCA 1的解释力为50.27%,说明样地的水分条件对年轮—气候关系有重要作用; 对于PCA 2,林龄也显著解释了25.58%的变异,说明林龄也显著影响了径向生长对历史气候变化的响应,进一步支持了图 2的结果。综上结果表明,年轮—气候关系主要受到样地气候条件的影响,其次受到林龄的影响,而蓄积量、林分密度等其余林分因子的影响则相对较弱。

      表 5  气候指标及林分条件对PCA前两轴得分的逐步回归分析方差分析

      Table 5.  Stepwise regression analysis of the scores of PCA 1 and PCA 2 with climatic factors and stand conditions

      项目ItemdfSS/%Sum SqMean SqFPr(>F)
      PCA 1
      PDSI Palmerdrought severity index150.2752.6352.6315.570.00***
      潜在蒸发量Potential evapotranspiration15.365.615.611.870.19
      林分密度Stand density17.187.517.512.510.14
      残差Residuals1337.1938.943.00
      PCA 2
      PDSI Palmerdrought severity index10.560.250.250.120.74
      林龄Stand age126.4812.0712.075.670.03*
      蓄积量Volume112.225.575.572.610.13
      残差Residual1360.7427.702.13
    • 本研究对处于不同气候和林分条件下的长白落叶松林的年轮—气候关系进行了分析,得到的结果与以往的相关研究存在一些异同。

    • 此前已有长白落叶松的相关研究表明,其径向生长对气候变化的响应与气候密切相关。如于大炮等[13]在长白山北坡的研究发现长白落叶松的年表—气候关系随海拔不同而有显著差异,并与林线位置有关。一般认为,长白落叶松的生长在高海拔主要受气温影响,在低海拔则主要受降水量影响[16, 36]。对兴安落叶松等其他树种的研究,也广泛发现了这种现象,这与低海拔温度较高、蒸散较强,导致的水分有效性较低有关[22, 37]。在本研究中,大 蒲柴河、六道和金苍经纬度差异并不是很大,但PDSI仍然解释了PCA 1高达50.27%的变异(表 5),这充分反映了气候条件确实对年轮—气候关系存在重要影响。本研究中所有样地海拔均不太高(550~800 m),结果表明,PDSI的作用远高于PET,即中低海拔长白山落叶松的年表—气候关系确实主要受水分条件的影响,与上述研究一致。不过,这并不能说明热量不起作用。在表 5中,如果将PET放在PDSI之前进入模型,则PET可显著解释PCA 1的30.05%的变异(表 6)。事实上,所谓气候干暖化常常是由于气候变暖导致的蒸散增加、有效水分减少引起的,未必一定是由于降水量的下降导致[38]。这一点在PCA排序图中也很明显,图 3中各个样地的温度朝PCA 1的右方增大,而PDSI箭头方向则刚好相反,说明PCA 1右侧样地的PDSI较低是由于其较高的温度(蒸散)所导致。综上,本研究中不同地点长白落叶松对历史气候变化响应的差异主要受水分条件的影响,但水分亏缺则随着热量的上升而加剧。这也意味着,随着东北地区今后进一步的气候变暖,长白落叶松(尤其是老龄林)的生长将会受到严重威胁。

      表 6  变量顺序不同的方差分析

      Table 6.  Appendix 1 variance analysis of variables in different order

      项目ItemdfSS/%Sum SqMean SqFPr(>F)
      PCA 1
      潜在蒸发量Potential evapotranspiration130.0531.4631.4610.500.01*
      PDSI Palmer drought severity index125.5826.7826.788.940.01*
      林分密度Stand density17.187.517.512.510.14
      残差Residual1337.1938.943.00
      PCA 2
      PDSI Palmer drought severity index10.560.250.250.120.74
      蓄积量Volume116.347.457.453.500.08
      林龄Stand density122.3610.1910.194.780.05
      残差Residual1360.7427.702.13
    • 研究表明,过去30年的气候变暖对不同林龄落叶松的影响有着很大差异(图 2):林龄较小的落叶松,RWI与各季节历史温度变化多呈正相关,说明气候变暖能够促进其生长; 而林龄较大的落叶松则多呈负相关关系,说明气候变暖对其径向生长有明显的负效应。

      年轮与历史气候变化的关系除了受气候条件影响之外,还受到林分条件的影响。具体来说,树木径向生长也受竞争以及林分株数密度、树种、林龄等因素的影响[23, 39]。而在本研究中,单从表 3的单因素解释力的结果来看,这些林分因子似乎没有影响。而多元回归模型的结果却表明,在同时考虑其他气候、林分因子时,林龄对年轮—气候关系存在显著影响(表 5)。

      一些研究表明,水分对径向生长的限制作用会随林龄增加而加剧,一般认为大树对气候的变化更加敏感[40]。但是,关于林龄的影响机制还存在争议。一方面,随着林龄的增加林分生产力下降是常见的现象,可能说明随着林龄增大、树木活力下降,因此对环境胁迫更加敏感。但随着林龄的增长,还伴随着树木大小(蓄积量)的增加,大树由于蒸腾作用强、对水分和养分的需求大,也可导致其对气候变化更敏感[41]。本研究结果表明,林龄对PCA2轴的解释力最大,解释了26.48%的变异,这和上述研究一致,说明林龄是影响年轮—气候关系的重要因素。表 5中,蓄积量的解释力不显著。但即使是将蓄积量放在林龄之前进入模型(表 6),其解释力在0.05水平上也不显著。因此,本研究结果更多地支持上述前一种解释,即林龄自身对长白落叶松的气候变化响应有着重要作用,而树木大小的作用相对是较弱的。不过,本文的结果与Martínez-Vilalta[41]的结果相反,因此,关于林龄和树木大小的影响,还有待于今后进一步研究才能确定。

      还有研究表明,林分密度是影响年轮—气候关系的一个重要因素[42]。而在本研究中,林分密度的解释力较弱且不显著(表 5)。本文中林分密度在340~1 960株/hm2之间,样地间差异足够大,因此可以排除是由于林分密度的梯度不够大导致其解释力不显著。即林分密度对年轮—气候关系的影响相对较弱,这一结果也与以往一些研究不同。因此,林分密度对径向生长的气候变化响应究竟有多大影响,也有待进一步探讨。

    • 图 1表明吉林北部山区的确存在明显的干暖化趋势,与以往东北地区研究一致。气候干暖化对东北树木生长的不利影响也得到众多研究证实。比如,对大兴安岭库都尔地区兴安落叶松的研究表明,区域性的干暖化趋势已经影响到当地主要树木的生长[14]。对阔叶红松林的研究也表明,虽然近几十年的气候干暖化对其不同演替阶段的径向生长影响复杂,正负效应都存在,但已经导致红松林生长的明显下降[19]。随着今后干暖化趋势的加剧,可以预料到蒸发、蒸腾作用的加剧从而导致土壤水分有效性的降低。这一威胁对于林龄较大的林分尤为严重,可能导致其严重的生长衰退甚至死亡[38],而成熟林却是天然林保护的重要对象。因此,需要深入研究如何采取适当的措施以减少干暖化的不利影响,以有效保护东北地区的森林生态系统。

参考文献 (42)

目录

    /

    返回文章
    返回