高级检索

留言板

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

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

基于Voronoi图的阔叶红松林空间结构特征分析

宋语涵 张晨 蔡体久 琚存勇

宋语涵, 张晨, 蔡体久, 琚存勇. 基于Voronoi图的阔叶红松林空间结构特征分析[J]. 北京林业大学学报, 2021, 43(1): 20-26. doi: 10.12171/j.1000-1522.20200056
引用本文: 宋语涵, 张晨, 蔡体久, 琚存勇. 基于Voronoi图的阔叶红松林空间结构特征分析[J]. 北京林业大学学报, 2021, 43(1): 20-26. doi: 10.12171/j.1000-1522.20200056
Song Yuhan, Zhang Chen, Cai Tijiu, Ju Cunyong. Quantitative analysis of spatial structural characteristics of broadleaved Korean pine forest based on Voronoi diagram[J]. Journal of Beijing Forestry University, 2021, 43(1): 20-26. doi: 10.12171/j.1000-1522.20200056
Citation: Song Yuhan, Zhang Chen, Cai Tijiu, Ju Cunyong. Quantitative analysis of spatial structural characteristics of broadleaved Korean pine forest based on Voronoi diagram[J]. Journal of Beijing Forestry University, 2021, 43(1): 20-26. doi: 10.12171/j.1000-1522.20200056

基于Voronoi图的阔叶红松林空间结构特征分析

doi: 10.12171/j.1000-1522.20200056
基金项目: “十三五”国家重点研发计划(2018YFC0507302),东北林业大学大学生创新训练项目(201810225122)
详细信息
    作者简介:

    宋语涵。主要研究方向:森林生态。Email:1074271550@qq.com 地址:150040 黑龙江省哈尔滨市和兴路26号东北林业大学林学院

    通讯作者:

    琚存勇,副教授。主要研究方向:生态学与3S应用。Email:qucy09@nefu.edu.cn 地址:同上

  • 中图分类号: S758.5+3

Quantitative analysis of spatial structural characteristics of broadleaved Korean pine forest based on Voronoi diagram

  • 摘要:   目的  森林空间结构是森林生物物理过程和动态演替的结果,又决定着森林的生长过程和生态服务功能的发挥。小兴安岭地区作为多样性保护与水源涵养的重要生态功能区,森林退化已经非常严重,有必要对现存的地带性顶极群落的结构特征进行分析,为退化森林的结构优化与调整提供理论支撑,促进其发挥应有的森林生态服务功能。  方法  本文以细鳞河国家级自然保护区内阔叶红松林为研究对象,根据野外调查数据,利用ArcGIS软件建立Voronoi多边形,确定空间结构单元,计算聚集指数、大小比数、竞争指数、多样性混交度等空间结构参数,分析阔叶红松林的空间结构特征。  结果  阔叶红松林林分水平格局为聚集分布,聚集指数与Voronoi多边形边数标准差的判定结果一致;针叶树比阔叶树有更大的竞争优势,阔叶树的种内竞争更激烈,但Hegyi竞争指数与大小比数二者结果并不完全一致;树种多样性混交度表明林分处于中度混交的水平,但每一树种的多样性混交度并不相同,针叶树的树种隔离程度比阔叶树的要高,并且这种隔离程度的强弱与树种在林分中的丰度没有直接关系。  结论  小兴安岭地区阔叶红松林以聚集分布、混交程度中等、竞争程度中等为其特征。红松与冷杉所在的结构单元混交度高、二者在各自结构单元内有明显的竞争优势;而枫桦、青楷槭、暴马丁香所在的结构单元混交度低,他们受到的竞争压力大,更多表现为种内竞争,尤以暴马丁香最具代表性。这些特征与文献报道的长白山区阔叶红松林的不尽相同,这既有选用参数不同的原因,也可能是地形气候差异的自然选择结果,还需要进一步调查研究。
  • 图  1  基于树木位置生成的Voronoi图

         圆点为树木位置。D表示两棵树之间的距离。Dots represent tree location. D represents the distance between two trees.

    Figure  1.  Voronoi diagram created based on tree location

    图  2  样地内乔木树种径级分布

    Figure  2.  Distribution patterns of DBH class of tree species in sample plots

    图  3  Voronoi多边形边数的频数分布

    Figure  3.  Frequency distribution of Voronoi polygon edge number

    图  4  优势树种大小比数的频率分布

    Figure  4.  Frequency distribution of size ratio of dominant tree species

    图  5  优势树种的平均大小比数(U)与竞争指数(CI/10)

    Figure  5.  Size ratio (U) and Hegyi’s competitionindex of dominant tree species (CI/10)

    图  6  林分与优势树种的多样性混交度分布频数

    Figure  6.  Frequency distribution of diversity mingling of dominant tree species and forest stand

    表  1  样地内优势树种的胸径信息

    Table  1.   Information about DBH of dominant tree species in sample plots

    主要树种
    Dominant tree
    species
    株数
    Plant
    number
    胸径 DBH/cm
    平均值
    Mean
    最小值
    Min. value
    最大值
    Max. value
    红松 Pinus koraiensis 51 38.1 4.7 84.6
    冷杉 Abies fabri 76 21.9 4.6 70.5
    青楷槭 Acer tegmentosum 82 12.3 4.7 45.0
    枫桦 Betula costata 94 11.7 5.3 30.7
    暴马丁香
    Syringa reticulata var.
    amurensis
    108 8.1 4.5 24.8
    下载: 导出CSV
  • [1] 汤孟平. 森林空间结构研究现状与发展趋势[J]. 林业科学, 2010, 46(1):117−122. doi:  10.11707/j.1001-7488.20100119.

    Tang M P. Advances in study of forest spatial structure[J]. Scientia Silavae Sinicae, 2010, 46(1): 117−122. doi:  10.11707/j.1001-7488.20100119.
    [2] 惠刚盈, von Gadow K, 赵中华, 等. 结构化森林经营原理[M]. 北京: 中国林业出版社, 2016.

    Hui G Y, von Gadow K, Zhao Z H, et al. Principles of structure-based forest management[M]. Beijing: China Forestry Publishing House, 2016.
    [3] Pommerening A. Evaluating structural indices by reversing forest structural analysis[J]. Forest Ecology and Management, 2006, 224(3): 266−277.
    [4] 卿东升, 李建军. 天然林空间结构多目标优化评价模型探究[J]. 林业科技通讯, 2019(2):3−6.

    Qing D S, Li J J. Research on multi-objective optimization evaluation model of natural forest spatial structure[J]. Forest Science and Technology, 2019(2): 3−6.
    [5] 夏富才, 姚大地, 赵秀海, 等. 长白山北坡阔叶红松林空间结构[J]. 东北林业大学学报, 2009, 37(10):5−7. doi:  10.3969/j.issn.1000-5382.2009.10.002.

    Xia F C, Yao D D, Zhao X H, et al. Horizontal spatial structure of Korean pine broad-leaved forest in northern slope of Changbai Mountains[J]. Journal of Northeast Forestry University, 2009, 37(10): 5−7. doi:  10.3969/j.issn.1000-5382.2009.10.002.
    [6] Ghalandarayeshi S, Nord-Larsen T, Johannsen V K, et al. Spatial patterns of tree species in Suserup Skov: a semi-natural forest in Denmark[J]. Forest Ecology and Management, 2017, 406: 391−401.
    [7] Gonçalves F M P, Revermann R, Cachissapa M J, et al. Species diversity, population structure and regeneration of woody species in fallows and mature stands of tropical woodlands of southeast Angola[J]. Journal of Forestry Research, 2018, 29: 1569−1579.
    [8] 李际平, 封尧, 赵春燕, 等. 基于Voronoi图的杉木生态公益林空间结构量化分析[J]. 北京林业大学学报, 2014, 36(4):1−7.

    Li J P, Feng Y, Zhao C Y, et al. Quantitative analysis of stand spatial structure of Cunninghamia lanceolata non-commercial forest based on Voronoi diagram[J]. Journal of Beijing Forestry University, 2014, 36(4): 1−7.
    [9] 胡艳波, 惠刚盈, 戚继忠, 等. 吉林蛟河天然红松阔叶林的空间结构分析[J]. 林业科学研究, 2003, 16(5):523−530. doi:  10.3321/j.issn:1001-1498.2003.05.002.

    Hu Y B, Hui G Y, Qi J Z, et al. Analysis of the spatial structure of natural Korean pine broadleaved forest[J]. Forest Research, 2003, 16(5): 523−530. doi:  10.3321/j.issn:1001-1498.2003.05.002.
    [10] Yang M, Cai T J, Ju C Y, et al. Evaluating spatial structure of a mixed broad-leaved/Korean pine forest based on neighborhood relationships in Mudanfeng National Nature Reserve, China[J]. Journal of Forestry Research, 2019, 30(4): 1375−1381.
    [11] 孙清芳, 刘滨凡, 马燕娥. 山河屯林业局红松阔叶混交林林分空间结构特征[J]. 森林工程, 2019, 35(6):1−5. doi:  10.3969/j.issn.1006-8023.2019.06.001.

    Sun Q F, Liu B F, Ma Y E. Spatial structure spatial structure characteristics of Pinus koraiensis mixed broad-leaved forest in Shanhetun Forestry Bureau[J]. Forest Engineering, 2019, 35(6): 1−5. doi:  10.3969/j.issn.1006-8023.2019.06.001.
    [12] Ali A. Forest stand structure and functioning: current knowledge and future challenges[J]. Ecological Indicators, 2019, 98: 665−677.
    [13] 郝月兰, 张会儒, 唐守正. 基于空间结构优化的采伐木确定方法研究[J]. 西北林学院学报, 2012, 27(5):163−168. doi:  10.3969/j.issn.1001-7461.2012.05.31.

    Hao Y L, Zhang H R, Tang S Z. Determination method of cutting tree based on forest stand structure optimization[J]. Journal of Northeast Forestry University, 2012, 27(5): 163−168. doi:  10.3969/j.issn.1001-7461.2012.05.31.
    [14] 侯红亚, 王立海. 小兴安岭红松针阔混交林林分空间结构分析[J]. 森林工程, 2012, 28(1):1−5. doi:  10.3969/j.issn.1001-005X.2012.01.001.

    Hou H Y, Wang L H. Spatial structure analysis of coniferous-broadleaved Korean pine mixed forest in Xiaoxing’anling Mountains[J]. Forest Engineering, 2012, 28(1): 1−5. doi:  10.3969/j.issn.1001-005X.2012.01.001.
    [15] 朱颖, 高路, 倪红伟. 黑龙江省国有重点林区森林资源动态变化及原因分析[J]. 森林工程, 2018, 34(3):32−38. doi:  10.3969/j.issn.1006-8023.2018.03.006

    Zhu Y, Gao L, Ni H W. Analysis on dynamic change and the reason of forest resources in the state-owned forest area of Heilongjiang Province[J]. Forest Engineering, 2018, 34(3): 32−38. doi:  10.3969/j.issn.1006-8023.2018.03.006
    [16] 汤孟平, 周国模, 陈永刚, 等. 基于 Voronoi 图的天目山常绿阔叶林混交度[J]. 林业科学, 2009, 45(6):1−5. doi:  10.3321/j.issn:1001-7488.2009.06.001

    Tang M P, Zhou G M, Chen Y G, et al. Mingling of evergreen broad-leaved forests in Tianmu Mountain based on Voronoi diagram[J]. Scientia Silavae Sinicae, 2009, 45(6): 1−5. doi:  10.3321/j.issn:1001-7488.2009.06.001
    [17] 汤孟平, 陈永刚, 施拥军, 等. 基于Voronoi图的群落优势树种种内种间竞争[J]. 生态学报, 2007, 27(11):4707−4716. doi:  10.3321/j.issn:1000-0933.2007.11.039.

    Tang M P, Chen Y G, Shi Y J, et al. Intraspecific and interspecific competition analysis of community dominant plant populations based on Voronoi diagram[J]. Acta Ecologica Sinica, 2007, 27(11): 4707−4716. doi:  10.3321/j.issn:1000-0933.2007.11.039.
    [18] 刘帅, 吴舒辞, 王红, 等. 基于 Voronoi 图的林分空间模型及分布格局研究[J]. 生态学报, 2014, 34(6):1436−1443.

    Liu S, Wu S C, Wang H, et al. The stand spatial model and pattern based on voronoi diagram[J]. Acta Ecologica Sinica, 2014, 34(6): 1436−1443.
    [19] 惠刚盈, von Gadow K, Albert M. 一个新的林分空间结构参数: 大小比数[J]. 林业科学研究, 1999, 12(1):1−6. doi:  10.3321/j.issn:1001-1498.1999.01.001.

    Hui G Y, von Gadow K, Albert M. A new parameter for stand spatial structure: neighbourhood comparison[J]. Forest Research, 1999, 12(1): 1−6. doi:  10.3321/j.issn:1001-1498.1999.01.001.
    [20] 汤孟平, 唐守正, 雷相东, 等. 两种混交度的比较分析[J]. 林业资源管理, 2004, 4(8):25−27.

    Tang M P, Tang S Z, Lei X D, et al. Comparison analysis on two minglings[J]. Forest Resource Management, 2004, 4(8): 25−27.
    [21] Clark P J, Evans F C. Distance to nearest neighbor as a measure of spatial relationships in populations[J]. Ecology, 1954, 35(4): 445−453.
    [22] 赵中华, 惠刚盈, 胡艳波, 等. 角尺度判断林木水平分布格局的新方法[J]. 林业科学, 2016, 52(2):10−16.

    Zhao Z H, Hui G Y, Hu Y B, et al. The new method judged horizontal distribution pattern by uniform angle index[J]. Scientia Silvae Sinicae, 2016, 52(2): 10−16.
    [23] Contreras M A, Affleck D, Chung W. Evaluating tree competition indices as predictors of basal area increment in western Montana forests[J]. Forest Ecology and Management, 2011, 262(11): 1939−1949.
    [24] 徐丽娜, 金光泽. 小兴安岭凉水典型阔叶红松林动态监测样地: 物种组成与群落结构[J]. 生物多样性, 2012, 20(4):470−481.

    Xu L N, Jin G Z. Species composition and community structure of a typical mixed broad-leaved-Korean pine (Pinus koraiensis) forest plot in Liangshui Nature Reserve, Northeast China[J]. Biodiversity Science, 2012, 20(4): 470−481.
    [25] 王琴香, 沈海龙, 和春庭, 等. 红松人工林和相邻次生林林下红松天然更新种群格局分析[J]. 森林工程, 2018, 34(2):16−20. doi:  10.3969/j.issn.1006-8023.2018.02.004.

    Wang Q X, Shen H L, He C T, et al. Analysis of natural regeneration pattern of Pinus koraiensis seedlings regenerated in Pinus koraiensis plantation and its adjacent Quercus mongolica natural forest[J]. Forest Engineering, 2018, 34(2): 16−20. doi:  10.3969/j.issn.1006-8023.2018.02.004.
    [26] 赵春燕, 李际平, 李建军. 基于 Voronoi 图和 Delaunay 三角网的林分空间结构量化分析[J]. 林业科学, 2010, 46(6):78−84. doi:  10.11707/j.1001-7488.20100612.

    Zhao C Y, Li J P, Li J J. Quantitative analysis of forest stand spatial structure based on Voronoi diagram & Delaunay triangulated network[J]. Scientia Silvae Sinicae, 2010, 46(6): 78−84. doi:  10.11707/j.1001-7488.20100612.
    [27] 张弓乔, 惠刚盈. Voronoi 多边形的边数分布规律及其在林木格局分析中的应用[J]. 北京林业大学学报, 2015, 37(4):1−7.

    Zhang G Q, Hui G Y. Analysis and application of polygon side distribution of Voronoi diagram in tree patterns[J]. Journal of Beijing Forestry University, 2015, 37(4): 1−7.
    [28] Canham C D, Lepage P T, Coates K D. A neighborhood analysis of canopy tree competition: effects of shading versus crowding[J]. Canadian Journal of Forest Research, 2004, 34(4): 778−787.
    [29] 惠刚盈, 胡艳波, 赵中华. 结构化森林经营研究进展[J]. 林业科学研究, 2018, 31(1):85−93.

    Hui G Y, Hu Y B, Zhao Z H. Research progress of structure-based forest management[J]. Forest Research, 2018, 31(1): 85−93.
    [30] Maleki K, Kiviste A, Korjus H. Analysis of individual tree competition on diameter growth of silver birch in Estonia[J/OL]. Forest Systems, 2015, 24(2): e023 [2019−12−25]. https://revistas.inia.es/index.php/fs/article/view/5742.
  • [1] 张岗岗, 刘瑞红, 惠刚盈, 张弓乔, 赵中华, 胡艳波.  林分空间结构参数N元分布及其诠释 . 北京林业大学学报, 2019, 41(4): 21-31. doi: 10.13332/j.1000-1522.20180228
    [2] 马得利, 孙永康, 杨建英, 史常青, 张艳, 张春禹.  基于无人机遥感技术的废弃采石场立地条件类型划分 . 北京林业大学学报, 2018, 40(9): 90-97. doi: 10.13332/j.1000-1522.20180051
    [3] 邱邦桂, 李永霞, 杨小林, 马和平.  拉萨半干旱河谷宜林地土壤水分时空分布格局 . 北京林业大学学报, 2016, 38(11): 9-15. doi: 10.13332/j.1000-1522.20150515
    [4] 郑俊强, 郭瑞红, 李东升, 李东, 李金功, 朱保坤, 韩士杰.  氮沉降和干旱对阔叶红松林凋落物分解的影响 . 北京林业大学学报, 2016, 38(4): 21-28. doi: 10.13332/j.1000-1522.20150464
    [5] 刘帅, 肖翠, 王均伟, 侯嫚嫚, 廖嘉星, 范秀华.  长白山阔叶红松林乔木幼苗年际动态及影响因素 . 北京林业大学学报, 2016, 38(11): 57-66. doi: 10.13332/j.1000-1522.20160012
    [6] 李一琳, 丁长青.  基于GIS和MaxEnt技术对濒危物种褐马鸡的保护空缺分析 . 北京林业大学学报, 2016, 38(11): 34-41. doi: 10.13332/j.1000-1522.20160134
    [7] 彭密, 国庆喜.  胜山阔叶红松林群落空间结构的最小面积研究 . 北京林业大学学报, 2016, 38(12): 21-27. doi: 10.13332/j.1000-1522.20150519
    [8] 吕寻, 杨双宝, 刘文桢, 郭小龙, 李安民, 袁一超, .  小陇山锐齿栎原始群落林木径阶空间结构特征 . 北京林业大学学报, 2015, 37(5): 11-18. doi: 10.13332/j.1000-1522.20140238
    [9] 陈亚南, 杨华, 马士友, 任玫玫.  长白山2种针阔混交林空间结构多样性研究 . 北京林业大学学报, 2015, 37(12): 48-58. doi: 10.13332/j.1000-1522.20150171
    [10] 汪平, 贾黎明, 魏松坡, 王奇峰.  基于Voronoi图的侧柏游憩林空间结构分析 . 北京林业大学学报, 2013, 35(2): 39-44.
    [11] 张震, 刘萍, 丁易, 刘黎明.  天山云杉林不同发育阶段种群分布格局研究 . 北京林业大学学报, 2010, 32(3): 75-79.
    [12] 韩路, 王海珍, 彭杰, 梁继业, 马春晖.  不同生境胡杨种群径级结构与格局动态研究 . 北京林业大学学报, 2010, 32(1): 7-12.
    [13] 王原, 吴泽民, 张浩, 赵霞.  基于RS和GIS的马鞍山市分区城市森林景观格局综合评价 . 北京林业大学学报, 2008, 30(4): 46-52.
    [14] 刘盛, 李国伟.  林分碳贮量测算方法的研究 . 北京林业大学学报, 2007, 29(4): 166-169.
    [15] 王芳, 吴彩燕, 黄荣凤, 李瑞, 胡万良, 殷亚方, 何亚平, 张璧光, 秦爱光, 颜绍馗, 魏潇潇, 杨平, 周永学, 白岗栓, 高黎, 张莉俊, 郑小贤, 刘杏娥, 邓小文, 毛俊娟, 王费新, 胡胜华, 袁怀文, 张洪江, 王正, 杜社妮, 王晓欢, 王兆印, 樊军锋, 张克斌, 李猛, 费世民, NagaoHirofumi, 汪思龙, 张岩, 乔建平, 赵天忠, 崔赛华, 常旭, 戴思兰, 谭学仁, 孙向阳, 王小青, 刘燕, 王胜华, 罗晓芳, 李昀, 刘云芳, KatoHideo, 陈放, 张双保, 龚月桦, 范冰, 徐嘉, 江泽慧, 韩士杰, 孔祥文, 江玉林, 高荣孚, 张旭, 李华, 张占雄, 王海燕, , 陈宗伟, 李媛良, 刘秀英, 任海青, IdoHirofumi, 侯喜录, 郭树花, 李晓峰, 陈秀明, , 丁磊, 常亮, 杨培华, , 张代贵, 李考学, 陈学平, , 徐庆祥, 蒋俊明, 薛岩, 高建社, 张桂兰, 费本华, 金鑫, 李雪峰, 涂代伦, 王晓东, , 刘永红, 续九如, 张红丽, 丁国权, , .  基于GIS与信息量模型的地层因素对三峡库区滑坡发育的影响关系 . 北京林业大学学报, 2007, 29(6): 138-142.
    [16] 李雪华, 王兰珍, 李国雷, 李义良, 方升佐, 李生宇, 汪杭军1, 张冬梅, 赵铁珍, 薛康, 刘鹏举, 周传艳, 韦艳葵, 王立海, 王旭, 耿玉清, 任强, HUALi_zhong, 吴丽娟, 段文霞, 刘剑锋, 朱小龙, 党文杰, 雷加强, 宋永明, 崔同林, 韩士杰, 周亮, 余新晓, 李建章, 杨慧敏, 朱波, 苏晓华, 高岚, 黎明, 刘勇, JIANGXi_dian, 方陆明, 刘勇, 尹光彩, 李振基, 何茜, 阎秀峰, 周宇飞, 周国逸, 杨娅, 鹿振友, 喻理飞, 程云清, 宗文君, 徐新文, 王清文, 张冰玉, 唐小明, 周国逸, 沈熙环, 孙向阳, 王春林, 王新杰, 徐扬, 刘锐, 玲, 赖志华, 李吉跃, HEXiu_bin, 虞木奎, 柯水发, 李晓兰, 王伟宏, 茹广欣, 张志毅, 陈实, 郭蓓, 周晓梅, 孙阁, 温亚利, , 宋爱琴, 李俊清, 3, 齐涛, 国庆, 陈培金, 李丙文, 陈峻崎, 张可栋, 王旭, 蒋德明, 唐旭利, 王晓静, 王建林, 姚永刚, 长山, 刘志明, 周玉平, 赵双荣, 关少华, 宋湛谦, 陈放, 王春林, 杨伟伟, 闫俊华, 郑凌峰.  基于GIS的喀斯特森林生态仿真系统的构建 . 北京林业大学学报, 2007, 29(2): 152-155.
    [17] 胡胜华, 袁怀文, 颜绍馗, 胡万良, 高黎, 杨平, 黄荣凤, 李瑞, 周永学, 白岗栓, 张洪江, 张莉俊, 王费新, 邓小文, 毛俊娟, 王芳, 魏潇潇, 吴彩燕, 张璧光, 郑小贤, 殷亚方, 刘杏娥, 何亚平, 秦爱光, 赵天忠, 谭学仁, 王正, 张岩, 孙向阳, 费世民, 王兆印, 王晓欢, 乔建平, 杜社妮, 樊军锋, 刘燕, 常旭, NagaoHirofumi, 戴思兰, 王胜华, 罗晓芳, 王小青, 张克斌, 崔赛华, 汪思龙, 李猛, 王海燕, 陈放, 江泽慧, 徐嘉, 张占雄, 范冰, 龚月桦, , 张旭, KatoHideo, 刘云芳, 韩士杰, 李昀, 张双保, 江玉林, 高荣孚, 李华, 孔祥文, , 李媛良, 常亮, 刘秀英, 杨培华, 丁磊, 侯喜录, 陈秀明, IdoHirofumi, 陈宗伟, 郭树花, 任海青, 李晓峰, 张代贵, 费本华, , , 蒋俊明, 陈学平, 张桂兰, 徐庆祥, 李考学, 高建社, 薛岩, 王晓东, 李雪峰, 金鑫, 涂代伦, , 刘永红, 续九如, , 张红丽, , 丁国权.  干扰对人工诱导的阔叶红松林群落结构及高等植物多样性的影响 . 北京林业大学学报, 2007, 29(6): 72-78.
    [18] 王勇, 张春雨, 何恒斌, 赵广亮, 于占源, 李秀芬, 于格, 何利娟, 马履一, 李笑吟, 赖巧玲, 杨永福, 张亚利, 吕兆林, 王华, 王献溥, 张力平, 李长洪, 钟健, 姜春宁, 曹金珍, 许月卿, 贾彩凤, 史军辉, 林峰, 郭小平, 李悦, 郝玉光, 朱清科, 王希群, 尚晓倩, 孙长霞, 王骏, 毕华兴, 习宝田, 李鸿琦, 赵博光, 邵晓梅, D.PascalKamdem, 黄忠良, 郭惠红, 贾桂霞, 于顺利, 鲁春霞, 赵秀海, 胥辉, 杨培岭, 杨明嘉, 王继兴, 曾德慧, 郑彩霞, 朱教君, 周金池, 尚宇, 包仁艳, 丁琼, 崔小鹏, 刘燕, 张榕, 陈宏伟, 姜凤岐, 李黎, 任树梅, 甘敬, 郑景明, 杨为民, 费孛, 王庆礼, 张志2, 欧阳学军, 谢高地, 丁琼, 王秀珍, 朱金兆, 蔡宝军, , 沈应柏, 贾桂霞, 张中南, 刘足根, 贾昆锋, 何晓青, 范志平, 沈应柏, 刘鑫, 张池, 刘艳, , 纳磊, 李凤兰, 陈伏生, 鹿振友, 周金池, 李林, 毛志宏, 张方秋, 唐小明, , 周小勇, 赵琼, 马玲, 申世杰, .  基于GIS和RUSLE的土壤侵蚀量计算——以贵州省猫跳河流域为例 . 北京林业大学学报, 2006, 28(4): 67-71.
    [19] 王华, 许月卿, 吕兆林, 郭小平, 林峰, 姜春宁, 赖巧玲, 张力平, 何恒斌, 王献溥, 马履一, 于格, 张春雨, 李笑吟, 王勇, 何利娟, 赵广亮, 李秀芬, 曹金珍, 史军辉, 李长洪, 张亚利, 贾彩凤, 钟健, 杨永福, 于占源, 杨培岭, 于顺利, 李悦, D.PascalKamdem, 朱教君, 郝玉光, 贾桂霞, 王希群, 曾德慧, 赵博光, 郑彩霞, 尚晓倩, 毕华兴, 孙长霞, 黄忠良, 王继兴, 朱清科, 邵晓梅, 郭惠红, 王骏, 赵秀海, 习宝田, 胥辉, 李鸿琦, 杨明嘉, 鲁春霞, 杨为民, 尚宇, 崔小鹏, 谢高地, 丁琼, 张榕, 王庆礼, 刘燕, 王秀珍, 欧阳学军, 张志2, 丁琼, 包仁艳, 甘敬, 姜凤岐, 郑景明, 任树梅, 陈宏伟, 费孛, 周金池, 朱金兆, 李黎, 刘足根, 范志平, 何晓青, 沈应柏, , 贾桂霞, 贾昆锋, 沈应柏, 张池, 张中南, 刘鑫, 蔡宝军, 刘艳, , 李凤兰, 毛志宏, 张方秋, 鹿振友, 陈伏生, 唐小明, 李林, 周金池, 纳磊, , 马玲, 周小勇, 赵琼, 申世杰, .  长白山阔叶红松林林隙大小结构研究 . 北京林业大学学报, 2006, 28(4): 34-38.
    [20] 黄华国, 邵海荣, 李黎, 赵晓松, 田晶会, 曹世雄, 张德荣, 贺庆棠, 王小丹, 董运斋, 王瑞刚, 李雪玲, 戴松香, 华丽, 马宇飞, 黄荣凤, 高岩, 郭明辉, 习宝田, 金幼菊, 李俊清, 贺庆棠, 于志明, 贺康宁, 王四清, 张晓丽, 关德新, 陈少良, 陈少良, 古川郁夫, 李文彬, 陈斌如, 闫丽, 冷平生, 阎海平, 刘力源, 李海英, 任云卯, 王百田, 邹祥旺, 王蕾, 鲍甫成, 高双林, 陈莉, 吴家兵, 高攀, 阎海平, 杨永福, 李建章, 侯智, 王金满, 金昌杰, 陈华君, 侯智, 张卫强, 金小娟, 郝志勇, 程根伟, 陈源泉, 赵有科, 韩士杰, 尹婧, 翁海娇, 杜建军, 高旺盛, 李涛, 李鹤, 赵琼, 杨爽, 段杉.  长白山阔叶红松林通量观测的footprint及源区分布 . 北京林业大学学报, 2005, 27(3): 17-23.
  • 加载中
图(6) / 表 (1)
计量
  • 文章访问数:  119
  • HTML全文浏览量:  28
  • PDF下载量:  28
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-03-20
  • 修回日期:  2020-04-08
  • 网络出版日期:  2021-01-11
  • 刊出日期:  2021-02-05

基于Voronoi图的阔叶红松林空间结构特征分析

doi: 10.12171/j.1000-1522.20200056
    基金项目:  “十三五”国家重点研发计划(2018YFC0507302),东北林业大学大学生创新训练项目(201810225122)
    作者简介:

    宋语涵。主要研究方向:森林生态。Email:1074271550@qq.com 地址:150040 黑龙江省哈尔滨市和兴路26号东北林业大学林学院

    通讯作者: 琚存勇,副教授。主要研究方向:生态学与3S应用。Email:qucy09@nefu.edu.cn 地址:同上
  • 中图分类号: S758.5+3

摘要:   目的  森林空间结构是森林生物物理过程和动态演替的结果,又决定着森林的生长过程和生态服务功能的发挥。小兴安岭地区作为多样性保护与水源涵养的重要生态功能区,森林退化已经非常严重,有必要对现存的地带性顶极群落的结构特征进行分析,为退化森林的结构优化与调整提供理论支撑,促进其发挥应有的森林生态服务功能。  方法  本文以细鳞河国家级自然保护区内阔叶红松林为研究对象,根据野外调查数据,利用ArcGIS软件建立Voronoi多边形,确定空间结构单元,计算聚集指数、大小比数、竞争指数、多样性混交度等空间结构参数,分析阔叶红松林的空间结构特征。  结果  阔叶红松林林分水平格局为聚集分布,聚集指数与Voronoi多边形边数标准差的判定结果一致;针叶树比阔叶树有更大的竞争优势,阔叶树的种内竞争更激烈,但Hegyi竞争指数与大小比数二者结果并不完全一致;树种多样性混交度表明林分处于中度混交的水平,但每一树种的多样性混交度并不相同,针叶树的树种隔离程度比阔叶树的要高,并且这种隔离程度的强弱与树种在林分中的丰度没有直接关系。  结论  小兴安岭地区阔叶红松林以聚集分布、混交程度中等、竞争程度中等为其特征。红松与冷杉所在的结构单元混交度高、二者在各自结构单元内有明显的竞争优势;而枫桦、青楷槭、暴马丁香所在的结构单元混交度低,他们受到的竞争压力大,更多表现为种内竞争,尤以暴马丁香最具代表性。这些特征与文献报道的长白山区阔叶红松林的不尽相同,这既有选用参数不同的原因,也可能是地形气候差异的自然选择结果,还需要进一步调查研究。

English Abstract

宋语涵, 张晨, 蔡体久, 琚存勇. 基于Voronoi图的阔叶红松林空间结构特征分析[J]. 北京林业大学学报, 2021, 43(1): 20-26. doi: 10.12171/j.1000-1522.20200056
引用本文: 宋语涵, 张晨, 蔡体久, 琚存勇. 基于Voronoi图的阔叶红松林空间结构特征分析[J]. 北京林业大学学报, 2021, 43(1): 20-26. doi: 10.12171/j.1000-1522.20200056
Song Yuhan, Zhang Chen, Cai Tijiu, Ju Cunyong. Quantitative analysis of spatial structural characteristics of broadleaved Korean pine forest based on Voronoi diagram[J]. Journal of Beijing Forestry University, 2021, 43(1): 20-26. doi: 10.12171/j.1000-1522.20200056
Citation: Song Yuhan, Zhang Chen, Cai Tijiu, Ju Cunyong. Quantitative analysis of spatial structural characteristics of broadleaved Korean pine forest based on Voronoi diagram[J]. Journal of Beijing Forestry University, 2021, 43(1): 20-26. doi: 10.12171/j.1000-1522.20200056
  • 森林空间结构体现了森林中树木及其属性在空间上的分布格局[1],即决定着森林的生长过程和生态服务功能的发挥,又是森林生物物理过程和动态演替的结果[2-3],它的优劣很大程度上决定着森林生态系统是否健康和稳定[4]。森林空间结构虽然暂无统一定义,但一般都强调与树木的空间位置有关,通常从树木的水平分布格局、树种的混交关系和树木的大小关系或竞争关系等3个方面进行描述[1]。自21世纪以来,不同森林的空间结构特征在国内外得到了广泛研究[4-11],结果表明不同地区不同森林类型的结构特征有明显差异,甚至同一地理分区的同一森林类型的结构特征也不相同,如夏富才等发现长白山北坡阔叶红松林呈现团状分布[5],而胡艳波等研究了同属长白山系位于吉林蛟河的6块红松阔叶林样地,林分分布格局有2块样地是均匀分布,4块样地为随机分布[9]。这些研究对正确评价森林的生态服务功能[12]、采伐木的科学选择[13]、次生林的优化和栖息地的重建[3]提供了可能。

    《林业发展“十三五”规划》在分区施策的基础上提出:在东北生态保育区要“强化森林经营,调整和优化森林结构”,原生林要“维护良好的森林结构和功能”,次生林要“调整和优化林分树种结构、龄组结构、径级结构和密度结构”。原生林的森林结构是当地气候地理条件下长期自然进化的结果,是最健康、最稳定的结构[5-7],要达到次生林结构优化的目标必须师法自然,学习参考原生林的结构格局[2-3,8]。阔叶红松林作为小兴安岭地区的地带性顶极植被,具有怎样的空间结构特征还不是很清楚,本文采用Voronoi图确定空间结构单元,计算聚集指数、大小比数、多样性混交度等结构参数和竞争指数,分析阔叶红松林群落空间结构现状特征与树木间竞争关系,以期为当地次生林的结构优化和调整提供参考,促进森林生态系统结构的稳定和生态服务功能的发挥。

    • 黑龙江细鳞河国家级自然保护区地处小兴安岭东南坡,为小兴安岭向三江平原的过渡地带,地理坐标为47°26′58″ ~ 47°37′37″N,130°03′26″ ~ 130°24′05″E,属于温带湿润大陆性气候,年平均温度为2.5 ℃,年降水量600 mm左右,无霜期147 d左右,土壤有暗棕壤、白浆土、草甸土等。该地区的地带性植被类型为阔叶红松林,主要树种有红松(Pinus koraiensis)、冷杉(Abies fabri)、水曲柳(Fraxinus mandshurica)、紫椴(Tilia amurensis)、枫桦(Betula costata)、春榆(Ulmus davidiana var. japonica)、青楷槭(Acer tegmentosum)、暴马丁香(Syringa reticulata var. amurensis)等[14-15]

    • 通过查询林相图与实地调查,在细鳞河国家级自然保护区核心区内选择阔叶红松林保存完好的区域设立100 m × 100 m的样地,并将样地划分成10 m × 10 m的相邻网格,以每个网格为调查单元,调查其中所有胸径大于4 cm乔木树种的胸径、树高及相对于固定角点的坐标等信息。

    • 林分空间结构特征是在分析林分每个空间结构单元的基础上进行的,空间结构单元通常由1株目标树与n株近邻树组成,至于n取多大一直存在着争议[8,13,16]。事实上,基于距离确定的n株树,并不一定都是目标树的近邻树,也可能包含了不与目标树直接竞争的个体[17]。而根据每株树的坐标生成的Voronoi图是按照对象最近属性将空间划分成许多单元区域,图中任一凸多边形所包围区域的点到该凸多边形的对象点的距离小于该点到其他任何对象点的距离[13]。因此利用Voroni图确定相邻木,就可以避免因近邻树数目选取过小而导致无法兼顾周围所有近邻树、或者选取过大而将实际上并非近邻的树纳入考量范围。在Arcmap分析工具箱(analysis tools)的邻域分析模块(proximity)中可创建Voronoi图(creat thiessen polygons)(图1),而相邻多边形功能(polygon neighbors)可以输出目标树及对应的近邻树。

      图  1  基于树木位置生成的Voronoi图

      Figure 1.  Voronoi diagram created based on tree location

    • 处在样地边缘的目标树,其近邻树可能位于样地之外,为消除这种边缘效应,将所在Voronoi多边形的一边为样地边缘的树木仅作为其他目标树的近邻树,而不作为计算结构参数的目标树[18],如图1中左下角最外侧的树。

    • 大小比数是衡量目标树在空间结构单元中优势程度的一个指标[2,19],它是结构单元中胸径大于目标树的近邻树株数占所考察的全部近邻树的比例,用公式表示为:

      $$U_{i} = \frac{1}{n}\sum\limits_{j = 1}^n {K_{i \! j}} $$

      式中:n为结构单元i中相邻树株数;若相邻树j比目标树i的胸径小,则Kij = 0,否则Kij = 1。则目标树的大小比数越小,表示该树种在结构单元中越有优势。

      目前用于描述树木空间隔离程度的常用指标是Gadow提出的混交度[4-6,8,13,16],此指标仅关注近邻树与目标树的树种异同,而忽略了近邻树的种间差异,因此本文采用汤孟平提出的多样性混交度描述空间结构单元中的树种隔离程度[20]

      $${{M_i = }}\frac{{{{{n}}_{{i}}}}}{{{{{n}}^{{2}}}}}\sum\limits_{{{j}} = 1}^n {{v_{i \! j}}} $$

      式中:Mi表示目标树i的树种多样性混交度,ni为目标树in株最近邻树中树种的个数,n为目标树的最近邻树株数,当目标树i与近邻树 j属不同树种时vij = 1,同种则为0。多样性混交度越大,表明树种之间的隔离程度越大、各树种株数越接近。从林分的角度看,多样性混交度为0,表明该林分为纯林;多样性混交度为1,则证明该林分各树种株数相等且完全隔离。当然,完全隔离的情况在天然情况下几乎不存在[20]

      水平分布格局采用聚集指数判定,其计算方法为:

      $$ {{R}} = {\bar r_{\rm{A}}}/{\bar r_{\rm{E}}}$$

      式中:$ {\bar {r}}_{{\rm{A}}}=\dfrac{1}{N}\displaystyle\sum _{i=1}^{N}{r}_{i} $${\bar r_{\rm{E}}} = 0.5\sqrt {F/N} $R为聚集指数,ri为第i株目标树到它最近邻树的距离,N为样地内树木株数,F为样地面积。再根据正态分布检验的原则,$\; \mu =\dfrac{{\bar {r}}_{{\rm{A}}}-{\bar {r}}_{{\rm{E}}}}{\sigma } $$ \mathrm{\sigma }=\dfrac{0.261 \; 36}{\sqrt{{N}^{2}/F}} $,可判断样地分布格局,即按0.05的显著水平,若$ \left|\mu \right| $ < 1.96,则可判断为随机分布,若$ \left|\mu \right| $ > 1.96,当R < 1则判断为聚集分布,当R > 1则为均匀分布[21-22]

    • 目标树胸径的大小虽然一定程度上能够说明在结构单元内受到的竞争强弱,但是距离的远近也是竞争是否存在的重要决定因素[23]。本文以Hegyi竞争指数描述树木竞争状况,其计算公式如下[1,23]

      $${\rm{CI}}_{i} = \sum\limits_{i = 1}^n {\left( {\frac{1}{{{L_{i \! j}}}} \times \frac{{{D_{{j}}}}}{{{D_i}}}} \right)} $$

      式中:CIi表示目标树i的竞争指数,Lij表示目标树i和近邻树j之间的距离,DiDj分别表示目标树i和近邻树j的胸径。CI越大,表明目标树受到的竞争越激烈,反之,受到的竞争越弱。

    • 1 hm2样地内共调查到胸径 > 4 cm的乔木(含亚乔木)721株,分属24个种,平均胸径16.0 cm,平均树高12.1 m,胸高断面积为25.32 m2/hm2。所有乔木的径级分布如图2所示,胸径与其株数近似呈倒J型分布,随着胸径的增大,其株数越来越小。这一特征与小兴安岭其他地区原始红松林相一致[24-25]。去除边缘效应后,还有659株乔木可以参与空间结构参数的分析,其中数量排前5位树种的统计特征如表1所示。

      表 1  样地内优势树种的胸径信息

      Table 1.  Information about DBH of dominant tree species in sample plots

      主要树种
      Dominant tree
      species
      株数
      Plant
      number
      胸径 DBH/cm
      平均值
      Mean
      最小值
      Min. value
      最大值
      Max. value
      红松 Pinus koraiensis 51 38.1 4.7 84.6
      冷杉 Abies fabri 76 21.9 4.6 70.5
      青楷槭 Acer tegmentosum 82 12.3 4.7 45.0
      枫桦 Betula costata 94 11.7 5.3 30.7
      暴马丁香
      Syringa reticulata var.
      amurensis
      108 8.1 4.5 24.8

      图  2  样地内乔木树种径级分布

      Figure 2.  Distribution patterns of DBH class of tree species in sample plots

    • 根据Voronoi多边形所确定的结构单元中所含林木最多为12株,即一株目标树的周围有11株近邻树,但只有一个此类型结构单元;所含林木最少的结构单元有4株树构成,以含6株树木的结构单元为最多(图3)。所有结构单元中近邻树株数的平均值为5.85,约为6株,这与其他研究相吻合[2,18,26],说明不同林型间用Voronoi图确定的结构单元的大小规律几乎没有差异。

      图  3  Voronoi多边形边数的频数分布

      Figure 3.  Frequency distribution of Voronoi polygon edge number

      经过计算,阔叶红松林样地的聚集指数为0.736,$ \left|\mu \right| $为12.963,说明样地内的林木整体表现为聚集分布。张弓桥与惠刚盈指出[27],虽然Voronoi多边形边数不因分布格局的不同而有明显差异,总是近似呈正态分布,但不同分布格局的多边形边数的标准差之间存在一定差异,利用这一规律也可判断林分的水平分布格局,并且随机分布的标准差分布范围为[1.264,1.402]。此阔叶红松林的Voronoi多边形边数标准差为1.481,说明此林分的水平分布格局为聚集分布,这与聚集指数的判断结果相一致。

    • 不同树种的大小比数分布趋势明显不同,针叶树种特别是红松处于优势的情形更多(图4),说明它们在演替过程中获得了更多的水热资源,它们的近邻树多数都是胸径小于它们的个体,阔叶树尤其是暴马丁香生长处于明显劣势,其大小比数的平均值达到了0.645,针阔叶树种大小比数频数分布的这种规律与其他地区阔叶红松林相似[5,9-10]

      图  4  优势树种大小比数的频率分布

      Figure 4.  Frequency distribution of size ratio of dominant tree species

      针叶树在林分中的竞争优势从竞争指数上也能反映出来(图5,图中竞争指数已经缩小10倍),竞争指数越小,说明其处在一个更有利于生长的环境。但竞争指数与大小比数反应出的规律并不完全相同,从大小比数看,青楷槭比暴马丁香受到的竞争威胁要小,但竞争指数却反应了相反的趋势。由于竞争指数不仅考虑了结构单元中树木的胸径大小而且还把树木间的距离作为考量因子,因此更能准确地衡量目标树与近邻树的竞争关系,而经常被作为树木生长模型的重要参数[23]

      图  5  优势树种的平均大小比数(U)与竞争指数(CI/10)

      Figure 5.  Size ratio (U) and Hegyi’s competitionindex of dominant tree species (CI/10)

    • 整个林分的平均多样性混交度为0.497,处于中等混交的水平。5个丰度最高的树种中,红松与冷杉的混交度分布趋势基本一致,混交程度都是中等偏上的居多,说明它们的近邻树至少有一半以上与它们不同;而阔叶树的混交程度则趋势相反,其近邻树多属同一树种,树种隔离程度较低,尤其以暴马丁香最为突出(图6),说明它们的生存压力更多的是来自种内竞争。进一步分析样地内仅有的37株水曲柳的混交度,分布在5个混交度区间中的前两个区间相对频数为40.5%,大于后两个区间相对频数24.3%(其中最后一个区间分布频率为0%),说明阔叶树种内竞争较多的这种趋势与树种丰度并没有直接关系。虽然全林分与阔叶树种混交度程度偏低的现象与胡艳波等[9]、Yang等[10]和孙清芳等[11]报告的不同,我们认为主要是由于混交度计算方式的不同引起的。当前普遍采用的Gadow提出的混交度往往会夸大混交程度[20]。如在一个含4棵近邻树的结构单元中,只要与目标树不是一个物种,近邻树3株同种或不同种得到的Gadow混交度都是0.75,而用本文的计算方法,3株同种的情形多样性混交度只有0.25。

      图  6  林分与优势树种的多样性混交度分布频数

      Figure 6.  Frequency distribution of diversity mingling of dominant tree species and forest stand

    • 森林结构量化分析的基础是对象树与近邻树组成的结构单元,结构单元的大小目前还没有达成统一。惠刚盈等研究了1 000块完全随机分布林分的角尺度均值与近邻树株数的关系,认为近邻树株数为4时最经济又不会误判水平格局类型。尽管林分中树种与个体大小为非随机分布时,结构参数中混交度和大小比数都随着近邻树株数而变化,但考虑到结构分析应该建立在同一个结构单元上,建议选4株近邻树进行空间结构分析[2]。这一原则在大多数的森林结构分析中得到了应用[5-6,9-11]。这种固定数量的近邻树选取方法会把非直接相邻的树木作为近邻树(如图1中的左上部3号树会被选为1号的近邻树),违背了近邻树的选取初衷,有的学者尝试用能够确定目标树唯一利用空间的Voronoi图来选择近邻树[8,16-18,26],并且已有研究发现除了角尺度差异较大,其他指标两种方法计算结果具有一致性[13]。然而,基于Voronoi图确定近邻树也会把距离较远的树木选为近邻树或竞争树,如图1中编号4和5的两株,它们的间距超过了实际可能发生对称性竞争关系的距离[28],这样就会夸大近邻树对目标树的影响。在林分密度低、聚集分布趋势明显时,用Voronoi图确定近邻树应慎重。

      Voronoi图确定的结构单元包含的树木株数不同,这为角尺度的应用带来不便,本文利用聚集指数与Voronoi多边形边数标准差判断林分分布格局为聚集分布(团状分布),符合根据图1的直观判断。黑龙江细鳞河国家级自然保护区阔叶红松林的这种分布格局与夏富才等报告的长白山北坡的相同[5],但与小兴安岭凉水地区[14]、长白山系张广才岭地区[9]和老爷岭地区[10]的不一样。这说明师法自然的次生林结构调整过程中需要注意因地制宜,不能搞一刀切。

      大小比数作为衡量目标树在结构单元中竞争优势的一个参数,正确反映了黑龙江细鳞河国家级自然保护区阔叶红松林中红松与冷杉在林分中的优势生态位。在森林的结构优化中,可以用它调整树木间竞争关系[29],郝月兰等把大小比数超过林分平均值的目标树确定为择伐对象,使保留木能更好生长,并指出择伐对象的选择除了考虑大小比数,也应结合临界距离进行确定[13]。换句话说,大小比数结合竞争指数在确定择伐木时可能更合适。国外的不少研究表明以距离为参数的竞争指数能更好地描述复杂林分结构多变的近邻竞争关系,它比非距离相关的竞争指数能更准确地预测林木生长的竞争压力[23,30]

      本文发现,同样是阔叶红松林,特别是细鳞河自然保护区与凉水自然保护区又同属小兴安岭地区,但林分水平分布格局却不相同,是微地形小气候的原因还是林分演替阶段不同的原因造成的还不清楚,这需要进一步扩大样本量进行研究。在不确定顶极植被群落空间结构特征尺度效应的情况下,要调整优化当地次生林甚至人工林的森林结构、提高其生态服务功能还需要做更多细致的前期工作。

参考文献 (30)

目录

    /

    返回文章
    返回