桉树人工林单木地上生物量空间效应分析

罗大鹏, 农明川, 李会朋, 王俊峰, 欧光龙

罗大鹏, 农明川, 李会朋, 等. 桉树人工林单木地上生物量空间效应分析[J]. 西南林业大学学报(自然科学), 2022, 42(2): 120–129 . DOI: 10.11929/j.swfu.202102012
引用本文: 罗大鹏, 农明川, 李会朋, 等. 桉树人工林单木地上生物量空间效应分析[J]. 西南林业大学学报(自然科学), 2022, 42(2): 120–129 . DOI: 10.11929/j.swfu.202102012
Luo Dapeng, Nong Mingchuan, Li Huipeng, Wang Junfeng, Ou Guanglong. Spatial Effect Analysis on Individual Aboveground Biomass of Eucalyptus Plantation[J]. Journal of Southwest Forestry University, 2022, 42(2): 120-129. DOI: 10.11929/j.swfu.202102012
Citation: Luo Dapeng, Nong Mingchuan, Li Huipeng, Wang Junfeng, Ou Guanglong. Spatial Effect Analysis on Individual Aboveground Biomass of Eucalyptus Plantation[J]. Journal of Southwest Forestry University, 2022, 42(2): 120-129. DOI: 10.11929/j.swfu.202102012

桉树人工林单木地上生物量空间效应分析

基金项目: 国家自然科学基金项目(31760206, 31660202)资助;云南省唐守正院士工作站(2018IC066)资助;云南省王广兴专家工作站(2018IC100)资助;云南省万人计划青年拔尖人才专项(YNWR–QNBJ–2018–184)
详细信息
    作者简介:

    罗大鹏(1996—),男,硕士研究生。研究方向:森林测计学。Email: 376421454@qq.com

    通讯作者:

    欧光龙(1983—),男,博士,教授,博士生导师。研究方向:森林测计学。Email: olg2007621@126.com

Spatial Effect Analysis on Individual Aboveground Biomass of Eucalyptus Plantation

  • 摘要: 以桉树人工林皆伐样地调查数据为研究对象,记录样地内林木的空间位置坐标,并实测每株林木的木材、树皮、树枝、树叶的生物量,采用Ripley's K函数描述林木地上部分各维度生物量的空间分布模式,通过全局和局部Moran's I描述其空间自相关特征,基于组内方差定量描述各维度生物量的空间异质性变化。结果表明:桉树人工林单木地上生物量变化具有尺度效应,即在不同的距离尺度下表现出不同的聚集或离散特征,除树叶生物量在18.5~21.5 m范围内呈现显著聚集分布外,其余均未表现出显著的聚集分布特征;桉树人工林地上各维度生物量的最佳距离带宽均为12.5 m,且各维度均存在一定的空间自相关,全局Moran's I范围在0.0160~0.0289;各维度局部空间自相关变化趋势相对一致,各维度均有近15%的林木有显著空间自相关;随着距离尺度的增加,桉树人工林各维度生物量空间异质性均逐渐增大,且在距离大于5 m后逐渐平稳,所有生物量维度中树叶生物量的空间变异程度最大。因此,桉树人工林各维度单木生物量存在一定程度的空间自相关和空间异质性,且各维度间差异不大。
    Abstract: Taking the clear-cutting plot survey data of Eucalyptus plantation as the research object, the location coordinates of each tree in the plot were recorded, and the biomass of wood, bark, branch and foliage for each tree was measured. Then, the spatial distribution patterns were described for the wood, bark, stem, branches, foliage, crown and aboveground biomass using Ripley's K function, the spatial auto-correlation characteristics were analyzed by both global and local Moran's I, and the spatial heterogeneity of each biomass component was quantitatively described using the intra-variance index. The results showed that the aboveground biomass of Eucalyptus plantations had scale effect, and it showed the different aggregation or dispersion characteristics at different distance scales. Except for the foliage biomass shows a significant aggregation distribution if the distance is in the range of 18.5–21.5 m, and spatial distribution characteristics of the others are insignificant aggregation. Moreover, the optimal distance bandwidth for all biomass components of Eucalyptus plantation are 12.5 m, and spatial auto-correlation could be found in each component at a certain content, and the global Moran's I from 0.0160 to 0.0289. Meanwhile, the local spatial auto-correlation change of all biomass components are similar relatively, and nearly 15% of the trees in each component have a significant spatial auto-correlation. The spatial heterogeneity of all biomass components increase gradually with the increases of the distance, and the spatial variation would be stable if the distance is greater than 5 m, the spatial variation of foliage biomass is the largest one in all biomass components. In sum, there is spatial auto-correlation and spatial heterogeneity in the individual aboveground biomass in the Eucalyptus plantation at a certain degree, and the differences among the different biomass components are not significant.
  • 森林生态系统作为地球上维系生态平衡和物质循环的主要生态系统,在全球大气循环研究中居重要地位[1]。森林生物量约占陆地生态系统植被总生物量的90%,在维持全球气候稳定、调节全球碳平衡、减缓大气温室气体浓度上升等方面具有不可替代的作用[2]。对森林生物量的研究关系到生态系统生产力、碳循环,是全球变化研究的基础参数,同时森林生物量也是表征植物活动的关键变量[3]。森林生物量普遍存在空间效应,忽视森林生物量的空间效应可能会带来生物量变化分析的有偏估计或者误差[4-5]

    空间效应分析可以用于检验森林生物量在空间上的相关性强弱以及在空间上的变异性[6],主要表现为空间异质性和空间自相关性两个方面[7-8]。空间异质性普遍存在于自然界中[9],它是指系统或系统属性在空间分布中的复杂性和变异性,系统属性可以是生物量,土壤含氮量等生态学变量[10];空间异质性在一定程度上反映森林生物量在空间上的差异,且空间异质性的评价方法随空间数据类型的不同而存在差异,空间异质性的研究理论常用的有地统计学方法,也有研究使用组内方差对空间异质性进行分析[11-12]。空间自相关是一种表现地理数据在空间上的依赖程度的统计方法[13-14],研究对象集聚或离散效果最显著时的距离阈值,分析变量间的空间自相关,可以大幅度提升空间格局分析的精度[15]。反映整个研究区域内总体的全局空间自相关指数,以及描述区域内研究对象与其相邻单元间的空间相关性的局域空间自相关指数(LISA)是分析描述空间自相关的两个重要方面[16-20]。关于森林生物量的空间效应分析,农明川等[5]对昆明市云南松(Pinus yunnanensis)生物量的空间分布规律进行研究,指出昆明市云南松地上生物量呈现显著的空间自相关关系;张维生[17]以空间自相关分析方法为基础,对黑龙江省的森林和地貌等因子进行了研究,指出该地区的森林、地貌等因子的存在空间正相关性;王维芳等[18]也基于相似的方法对帽儿山地区的森林生物量的空间自相关性进行了分析,同样得出了相似的结论。但是这些分析在一定程度上揭示在区域尺度上森林生物量的空间分布模式,以及森林生物量分布的空间自相关性;尤其是在林分尺度上,揭示林内单木生物量间的空间效应变化规律的研究较少,欧光龙等[20]采用空间回归模型构建思茅松(Pinus kesiya)不同维度的单木生物量模型,取得了较好的拟合效果,这在一定程度上说明了各维度单木生物量是存在空间效应的,但是没有对空间及邻近木关系进行研究;徐美玲等[21-22]对云南省思茅松林的混交度、大小比数和角尺度几个空间结构参数进行量化分析,结果表明思茅松天然林整体呈现均匀分布。这些研究表明林木单木各维度生物量存在空间效应,但是基于实测数据,准确描述林内林木各维度间的空间效应关系的研究仍然缺乏。

    桉树(Eucalyptus spp.)由于适应环境能力强、生长迅速,经营周期短,单位面积产量高等特点,成为了热带和亚热带地区最重要的人工林树种[23],也被认为是世界人工造林三大速生树种之一,是中国南部省区极其重要的速生丰产林造林树种。分析桉树人工林单木生物量各维度间的空间效应,揭示其生物量空间变化规律对于准确描述林木生物量生长变化规律至关重要。基于此,本研究以云南省楚雄州桉树人工林皆伐样地(100 m×30 m)的单木各维度生物量(木材、树皮、树干、树枝、枝叶、树冠以及地上生物量)的实测数据为基础,采用空间自相关和空间异质性分析方法,对楚雄州东华镇桉树人工林地上生物量的空间分布规律进行研究[24],以获知桉树人工林各维度生物量的空间分布差异。在此基础上深入分析各维度生物量空间分布差异的成因,并确定桉树人工林单木各维度生物量空间效应研究的最佳距离阈值,从而准确获取区域森林生物量的空间分布信息,以期为提高桉树人工林的经济效益和生态效益提供理论依据。

    东华镇位于云贵高原中部的云南省楚雄州西南部,地处东经101°21′29″~101°31′34″,北纬24°26′36″~25°01′16″,总面积448 km2,由于地处紫溪山南麓,整体地势呈东北低、西南高的特点。平均海拔1 895~1 990 m,属于亚热带季风气候,其年均温、年均降雨量、年均日照、平均降霜期分别为15.6 ℃、831 mm、2 422 h、96 d,楚雄州境内分布植被类型主要为中亚热带常绿阔叶林和针叶林[25]

    研究区桉树林为10年生人工纯林,林分保存密度为1667株/hm2。样地设置为矩形,样地面积为0.3 hm2,样地中心点海拔1783 m。桉树人工林样地,经过数据整理,共调查445株样木(图1)。样地调查记录径阶大于5 cm的林木的基本信息:位置、胸径、树高、冠幅、枝下高等信息。

    图  1  样地示意图
    圆圈大小与胸径成正比,横纵坐标均表示距离(单位:m)。
    Figure  1.  The map of study plots

    本研究所用数据为2019年12月于楚雄州实测调查,单木生物量的测定参考《林木生物量模型研究》[26]。树枝生物量和树叶生物量的测定采用全称重法,即先测出树枝和树叶质量,然后分别取样烘干后求出含水率进而推算出树枝和树叶的质量。对于胸径较小的树干采用全称重法测定生物量,而胸径较大的树干生物量则采用材积密度法测定,即通过对树干进行取样,测定样品的质量和体积(样品体积通过排水法测定)从而求算其密度,进而推算树干生物量[27],样木生物量信息见表1

    表  1  样木基本信息
    Table  1.  The basic characters of sample trees kg
    统计量木材
    生物量
    树皮
    生物量
    树干
    生物量
    树枝
    生物量
    树叶
    生物量
    树冠
    生物量
    地上
    生物量
    最小值1.980.122.100.090.120.222.64
    最大值487.1062.25549.3626.5720.9747.54596.90
    平均值71.688.7480.425.322.547.8688.27
    标准差72.449.2181.514.693.017.5588.89
    下载: 导出CSV 
    | 显示表格

    本研究采用的样地边缘校正方法为:Simulate Outer Boundary Values。此方法会在研究区域边界外,创建边界内所发现点的镜像点,以便校正边界附近的低估现象。将镜像与研究区域的边的最大距离范围相等的距离内的点,使用镜像的点会使边界点的相邻点估计更加精确[28]

    本研究选择Ripley's K函数来描述单木地上生物量的空间分布模式,K函数作为一种基于距离的概括要素空间分布累积特征的方法,常被用来描述点要素在空间范围内的相关性程度,反映要素随尺度变化而呈现出的空间分布模式[29],为了克服K函数的方差不稳定的问题,Besag[30]K函数进行开方线性转化为L函数,转化后的L函数公式如下:

    $ L\left(d\right)=\sqrt{\dfrac{{ A}\displaystyle {\sum }_{i=1}^{n}\displaystyle {\sum }_{j=1\text{,}j\ne i}^{n}k(i,j)}{\text{π}n(n-1)}} $

    (1)

    式中:A是研究区面积;d表示距离阈值;n为点要素数量;k(i,j)为权重,距离d的最大值取值应为最短边的一半,大于该距离可能导致分析结果不准确;通过比较真实点的L(d)和样地空间随机点(CSR)之间的d,若L(d)>0则表明在尺度为d时研究对象为聚集分布;L(d)<0则说明研究对象呈离散分布;而L(d)=0则表明研究对象呈随机分布。

    全局Moran's I是全局空间自相关分析的重要指标,它能反映出研究区域内的对象的总体空间聚类模式[31]

    $ I=\left(\dfrac{n}{\displaystyle \sum _{i=1}^{n}\displaystyle \sum _{j=1}^{n}{w}_{ij}}\right)\dfrac{\displaystyle \sum _{i=1}^{n}\displaystyle \sum _{j=1}^{n}{w}_{ij}\left({x}_{i}-\bar{x}\right)\left({x}_{j}-\bar{x}\right)}{\displaystyle \sum _{i=1}^{n}{\left({x}_{i}-\bar{x}\right)}^{2}} $

    (2)

    This page contains the following errors:

    error on line 1 at column 1: Start tag expected, '<' not found

    Below is a rendering of the page up to the first error.

    全局Moran's指数I的取值意义随Z(I)的阈值的不同而存在差异。

    $ Z\left(I\right)=\dfrac{I-E\left(I\right)}{\sqrt{\mathrm{V}\mathrm{a}\mathrm{r}\left(I\right)}} $

    (3)

    $ E\left(I\right)=\dfrac{-1}{n-1} $

    (4)

    式中:Z(I)是空间聚类模式强度的衡量指标;E(I)和Var(I)分别是指数值I的期望值和方差;在给定的置信水平下(假定α=0.1),若I<0,Z(I)<−1.96表示空间要素间存在显著的空间负相关性。若I<0,Z(I)<−2.58表示空间要素间存在极显著的空间负相关性,表现为相异聚集(也称异常值),即相邻要素的属性值间差异较大的聚集在一起;若I>0,Z(I)>1.96表示观测对象间存在显著的空间正相关关系。若I>0,Z(I)>2.58表示观测对象间存在极显著的空间正相关关系,呈现出相似聚集的现象,即相邻要素属性值间差异较小的相互聚集;若I=0或−1.96≤Z(I)≤1.96则表示观测对象间不存在空间相关性,随机分布于研究区域内[32]

    局部Moran's指数可用于反映局部区域内特定属性值间是否存在空间自相关关系。在本研究中,它是一种弥补全局Moran's指数仅能对总体的空间自相关情况进行分析而不能反映特定地点研究对象与其相邻要素间的空间相关性的缺陷的方法[33]。局部Moran's指数(Ii)的计算公式如下所示:

    $ {I}_{i}=\dfrac{{n}^{2}}{\displaystyle \sum _{i=1}^{n}\displaystyle \sum _{j=1}^{n}{w}_{ij}} \cdot \dfrac{\left({x}_{i}-\bar{x}\right)\displaystyle \sum _{j=1}^{n}{w}_{ij}\left({x}_{j}-\bar{x}\right)}{\displaystyle \sum _{i=1}^{n}{\left({x}_{i}-\bar{x}\right)}^{2}} $

    (5)

    This page contains the following errors:

    error on line 1 at column 1: Start tag expected, '<' not found

    Below is a rendering of the page up to the first error.

    空间要素的空间异质性可以通过组内方差进行描述,组内方差可定量地描述研究对象的局部空间异质性,是分组距离大小的函数,通常会随着分组距离的增加而增大[34],其表达式如下:

    $ {S}_{\mathrm{i}\mathrm{n}\mathrm{t}\mathrm{r}\mathrm{a}}=\dfrac{1}{B}\displaystyle \sum _{g=1}^{B}\dfrac{1}{{n}_{g}}\displaystyle \sum _{h=1}^{{n}_{g}}{({X}_{gh}-\overline{{X}_{g}})}^{2} $

    (6)

    This page contains the following errors:

    error on line 1 at column 1: Start tag expected, '<' not found

    Below is a rendering of the page up to the first error.

    使用Excel 2019软件对实测数据进行标准化处理和相关表格的制作,使用ArcGIS 10.6软件下的多距离空间聚类分析工具进行空间分布模式的分析,增量空间自相关工具进行全局空间自相关分析,聚类与异常值分析工具进行局部空间自相关分析,空间异质性的组内方差分析在SPSS 20.0中完成。

    从单木各维度生物量的Ripley's K函数经变换后的L函数的变化曲线可以看出,除树叶生物量(图2e)外,桉树人工林各维度生物量空间聚集模式呈现相似的规律,均以2 m为界,2 m以内均呈现离散分布,2~14 m左右呈现聚集分布,15 m以后均呈现离散分布(图2)。其中木材生物量空间分布格局在0~2 m范围内离散分布,2~14 m范围内表现为聚集分布,14~25 m范围内表现为离散分布(图2a);树皮生物量在0~2 m范围内离散分布,2~15 m范围内表现为聚集分布,15~25 m范围内表现为离散分布(图2b);树干生物量空间分布格局在0~2 m范围内离散分布,2~14.5 m范围内表现为聚集分布,14.5~25 m范围内表现为离散分布(图2c);树枝生物量空间分布格局在0~2 m范围内离散分布,2~14.5 m范围内表现为聚集分布,14.5~25 m范围内表现为离散分布(图2d);树叶生物量在0~1 m之间呈现聚集分布,1~3.5 m为离散分布状态,3.5~25 m呈现聚集分布趋势,且在18.5~21.5 m范围内呈现显著聚集分布(图2e);树冠生物量空间分布格局在0~2 m范围内离散分布,2~14.5 m范围内表现为轻度聚集分布,14.5~25 m范围内表现为离散分布(图2f);树枝生物量空间分布格局在0~2 m范围内离散分布,2~13.5 m范围内表现为轻度聚集分布,13.5~25 m范围内表现为离散分布(图2g)。总体来看,桉树人工林全林各维度生物量观测值均介于上下包迹线之间,且各维度生物量实际值更接近于上包迹线,聚集特性不是很明显。

    图  2  桉树人工林各器官维度生物量L函数变化图
    Figure  2.  L-function changes in different component biomass of Eucalyptus plantation

    桉树人工林全林的各维度生物量的增量空间自相关分析结果见图3。由图3可知,桉树人工林树枝生物量7 m、树叶生物量5.5 m和7 m、树冠生物量5.5 m和24~25.5 m时呈负的空间自相关,此时生物量呈现高值与低值相聚集的聚类模式。其他距离上均呈正空间自相关,生物量存在高值于高值或低值于低值聚集的聚类模式。木材、树皮、树干及地上生物量在5~30 m均呈正空间自相关,此时各维度生物量存在高值于高值或低值于低值聚集的聚类模式。

    观测到各维度生物量的Z值达到显著后的第一个聚集峰值所对应的距离均为12.5 m,也就是说当空间距离为12.5 m时空间相关性达到极限,因此,楚雄州桉树人工林的最佳分析距离尺度为12.5 m。随着距离尺度的增加,桉树人工林全林的木材、树皮、树干、树枝、树叶、树冠和地上生物量均表现出一定程度的空间自相关性且空间变化规律相似,但相关性不显著,Moran's I值范围0.0160~0.0289,单木各维度生物量间相关性较低。

    图  3  桉树人工林其他树种各器官维度生物量全局莫兰指数变化曲线
    Figure  3.  Curves of global Moran's index of different component biomass in Eucalyptus plantation

    图4可知,对于桉树全林各维度生物量而言,均表现出了不同程度的空间自相关关系,且各维度生物量的空间自相关变化特征相仿。在样地左侧基本呈现出明显的高值于低值聚集(HL)并伴有低值于低值聚集(LL)的情况;在样地右侧基本呈现出明显的低值于高值聚集(LH),在样地右侧还出现明显的高值于高值聚集(HH)并伴有低值于高值聚集(LH);由表2可知,从桉树人工林全林来看,尽管85%的生物量值没有呈现出显著的自相关性(NS),但是15%的生物量值显示出显著的空间自相关,其中4%的生物量值显示出极显著的空间自相关,说明各维度生物量在局部区域内存在显著聚集现象。

    图  4  桉树人工林各维度生物量局部Moran's 指数空间分布
    横纵坐标均表示距离(单位:m)。
    Figure  4.  Spatial distribution of local Moran's index of different component biomass in Eucalyptus plantation
    表  2  桉树人工林单木各维度生物量聚类模式统计
    Table  2.  Statistics of biomass clustering model of Eucalyptus plantation in each dimension
    显著性Z-score聚集模式木材生物量树皮生物量树干生物量树枝生物量树叶生物量树冠生物量地上生物量
    株数占比/%株数占比/%株数占比/%株数占比/%株数占比/%株数占比/%株数占比/%
    −1.96~1.96 NS 353 79.33 359 80.67 351 78.88 382 85.84 365 82.02 364 81.80 353 79.33
    HH 8 1.80 7 1.57 9 2.02 11 2.47 11 2.47 10 2.25 8 1.80
    HL 9 2.02 7 1.57 8 1.80 6 1.35 9 2.02 8 1.80 9 2.02
    LH 7 1.57 7 1.57 6 1.35 5 1.12 3 0.67 6 1.35 7 1.57
    LL 6 1.35 6 1.35 7 1.57 1 0.22 1 0.22 2 0.45 6 1.35
    −2.58~−1.96 HL 9 2.02 9 2.02 11 2.47 5 1.12 8 1.80 8 1.80 9 2.02
    LH 6 1.35 6 1.35 6 1.35 7 1.57 4 0.90 5 1.12 6 1.35
    1.96~2.58 HH 13 2.92 12 2.70 12 2.70 9 2.02 15 3.37 14 3.15 13 2.92
    LL 3 0.67 4 0.90 2 0.45 1 0.22 4 0.90 2 0.45 3 0.67
    ≤−2.58 HL 2 0.45 5 1.12 1 0.22 2 0.45 0 0.00 0 0.00 2 0.45
    LH 14 3.15 13 2.92 15 3.37 5 1.12 15 3.37 15 3.37 14 3.15
    ≥2.58 HH 14 3.15 9 2.02 16 3.60 6 1.35 9 2.02 8 1.80 14 3.15
    LL 1 0.22 1 0.22 1 0.22 5 1.12 1 0.22 3 0.67 1 0.22
    下载: 导出CSV 
    | 显示表格

    图5可知,以5 m为界限,在5 m以内随着距离的增加各维度生物量的空间变异性增加较为迅速,5 m以后增速较为平缓,但是一直保持缓慢增加的趋势。从维度来看,树叶生物量的空间变异程度明显大于其他维度的生物量。在0~30 m的小尺度范围内,桉树林人工林各维度生物量的空间变异较小,这说明桉树人工林各维度生物量随着距离尺度的增加,桉树人工林的空间变异逐渐增大,且在距离5 m后空间异质性变化趋于稳定。

    图  5  桉树人工林各维度生物量组内方差
    Figure  5.  Intra-block variance of different component biomass in Eucalyptus plantation

    已有学者就大尺度单一树种的森林生物量空间效应进行研究,但几乎没有研究基于实测数据描述单木水平维度生物量的变化特征[35-36]。本研究基于楚雄州0.3 hm2桉树人工林皆伐实测样地,分析了桉树人工林木材、树皮、树干、树枝、树叶、树冠和地上生物量7个维度生物量的空间分布规律。从空间分布模式看,桉树人工林全林各维度生物量聚集特性均不是很明显。这主要是由于桉树人工林林木分布相对均匀,样地距离尺度较小导致无大规模生物量聚集[37-40]

    从全局空间自相关分析结果来看,全林各维度生物量Morans' I值介于0.0160~0.0289,说明桉树各维度生物量均存在一定的空间自相关,桉树生长特性会影响其各维度生物量的空间变化,桉树人工林7~8 a即为成熟林,说明这一速生树种在成熟林阶段,林木已经出现了竞争和林木分化,这在一定程度上造成了林木生物量变化的空间自相关。Zhang等[41-42]采用地理加权回归(GWR)方法来探讨桉树人工林胸径和树高关系,指出桉树人工林的空间分布存在一定的空间自相关,表明桉树人工林单木地上部分各维度生物量和胸径变量在空间自相关变化特征是一致的。

    局部空间自相关结果表明,近15%样木在样地内呈现聚集分布,说明桉树人工林生物量存在局部聚集现象。生物量高值主要集中在样地右侧,生物量低值集中在样地左侧及中部,由于样地中部林木较为稀疏,固全林地上部分生物量值由高到低依次为样地右侧、样地左侧、样地中部,且各维度生物量空间分布差异不大。此外,桉树人工林单木生物量的空间分布还受到立地条件等因素的影响,调查样地右侧立地条件较好,林木长势也较好;样地中部土壤较为板结,且林木较为稀疏,长势相对较差。这也造成了立地条件好的右侧生物量高值聚集,立地条件较差的左侧和中部生物量低值聚集,这样的立地条件分布与单木地上生物量值分布一致。这在一定程度上也说明了要准确描述桉树人工林生物量变化规律,空间效应的分析是不可忽略的[43]

    桉树人工林各维度生物量空间变异性均随距离的增加而增大,均在距离5 m后逐渐平稳,且树叶生物量的空间变异程度最大。在给定的距离尺度下树叶生物量较其它维度生物量的空间变异性更大,这是由于树叶作为桉树光合作用固定能量的主要来源,相较于其他器官竞争更为明显,从而导致树叶生物量在空间上的变异程度更大,这也说明了林木各维度生物量在空间效应特征上存在差异。胥辉[27]同时用材积密度法和称重法对吉林省红石林业局的8个树种60株立木树干生物量进行测定,同时,考虑到树干不同位置的密度存在差异可能带来的生物量测定误差,本研究在使用材积密度法测定树干生物量时,以2 m为一个区分段分别对树干不同高度处的密度进行测定,从而提高树干生物量测定精度。

    本研究对桉树这一人工引进树种的单木地上各维度生物量空间效应的描述,可为今后其他森林类型的单木生物量空间效应变化分析提供理论支撑,还可将其应用于优化和调节森林结构以充分发挥森林的多种功能,为精准分析林业数据以及林业生产经营实践提供参考和科学依据。由于研究条件有限,本研究样地面积较小,树种单一,也没有对树根这一维度生物量进行研究,且对影响各维度生物量空间分布规律的成因和影响机制探讨深度有限,这些都是本研究的不足之处。今后可在条件允许的情况下扩大调查样地面积,获取更多的样本和数据,对林木生物量空间分布格局的形成机制开展更加深入的探讨,进一步研究更多树种的空间效应[37]

  • 图  1   样地示意图

    圆圈大小与胸径成正比,横纵坐标均表示距离(单位:m)。

    Figure  1.   The map of study plots

    图  2   桉树人工林各器官维度生物量L函数变化图

    Figure  2.   L-function changes in different component biomass of Eucalyptus plantation

    图  3   桉树人工林其他树种各器官维度生物量全局莫兰指数变化曲线

    Figure  3.   Curves of global Moran's index of different component biomass in Eucalyptus plantation

    图  4   桉树人工林各维度生物量局部Moran's 指数空间分布

    横纵坐标均表示距离(单位:m)。

    Figure  4.   Spatial distribution of local Moran's index of different component biomass in Eucalyptus plantation

    图  5   桉树人工林各维度生物量组内方差

    Figure  5.   Intra-block variance of different component biomass in Eucalyptus plantation

    表  1   样木基本信息

    Table  1   The basic characters of sample trees kg

    统计量木材
    生物量
    树皮
    生物量
    树干
    生物量
    树枝
    生物量
    树叶
    生物量
    树冠
    生物量
    地上
    生物量
    最小值1.980.122.100.090.120.222.64
    最大值487.1062.25549.3626.5720.9747.54596.90
    平均值71.688.7480.425.322.547.8688.27
    标准差72.449.2181.514.693.017.5588.89
    下载: 导出CSV

    表  2   桉树人工林单木各维度生物量聚类模式统计

    Table  2   Statistics of biomass clustering model of Eucalyptus plantation in each dimension

    显著性Z-score聚集模式木材生物量树皮生物量树干生物量树枝生物量树叶生物量树冠生物量地上生物量
    株数占比/%株数占比/%株数占比/%株数占比/%株数占比/%株数占比/%株数占比/%
    −1.96~1.96 NS 353 79.33 359 80.67 351 78.88 382 85.84 365 82.02 364 81.80 353 79.33
    HH 8 1.80 7 1.57 9 2.02 11 2.47 11 2.47 10 2.25 8 1.80
    HL 9 2.02 7 1.57 8 1.80 6 1.35 9 2.02 8 1.80 9 2.02
    LH 7 1.57 7 1.57 6 1.35 5 1.12 3 0.67 6 1.35 7 1.57
    LL 6 1.35 6 1.35 7 1.57 1 0.22 1 0.22 2 0.45 6 1.35
    −2.58~−1.96 HL 9 2.02 9 2.02 11 2.47 5 1.12 8 1.80 8 1.80 9 2.02
    LH 6 1.35 6 1.35 6 1.35 7 1.57 4 0.90 5 1.12 6 1.35
    1.96~2.58 HH 13 2.92 12 2.70 12 2.70 9 2.02 15 3.37 14 3.15 13 2.92
    LL 3 0.67 4 0.90 2 0.45 1 0.22 4 0.90 2 0.45 3 0.67
    ≤−2.58 HL 2 0.45 5 1.12 1 0.22 2 0.45 0 0.00 0 0.00 2 0.45
    LH 14 3.15 13 2.92 15 3.37 5 1.12 15 3.37 15 3.37 14 3.15
    ≥2.58 HH 14 3.15 9 2.02 16 3.60 6 1.35 9 2.02 8 1.80 14 3.15
    LL 1 0.22 1 0.22 1 0.22 5 1.12 1 0.22 3 0.67 1 0.22
    下载: 导出CSV
  • 毛学刚, 范文义, 李明泽, 等. 黑龙江长白山森林生物量的时空变化分析 [J]. 植物生态学报, 2011, 35(4): 371−379.
    王晓莉, 常禹, 陈宏伟, 等. 黑龙江省大兴安岭主要森林生态系统 生物量分配特征 [J]. 生态学杂志, 2014, 33(6): 1437−1444.
    李海奎, 雷渊才. 中国森林植被生物量和碳储量评估[M]. 北京: 中国林业出版社, 2010.

    Zhang L J, Ma Z H, Guo L. An evaluation of spatial autocorrelation and heterogeneity in the residuals of six regression models [J]. Forest Science, 2009, 55(6): 533−548.

    农明川, 付磊, 李晟, 等. 昆明市云南松林地上生物量空间效应分析 [J]. 西南林业大学学报(自然科学), 2020, 40(2): 117−124.
    罗君, 孙振亓, 张学斌. 基于ripley's K函数的绿洲景观格局演变分析: 以张掖市甘州区为例 [J]. 水土保持研究, 2019, 26(4): 224−231.
    杨志辉, 赵军, 朱国锋, 等. 含植被覆盖影响的石羊河流域土壤水分遥感估算及空间格局分析 [J]. 生态学报, 2020, 40(23): 8826−8837.

    Anselin L, Griffith D A. Do spatial effecfs really matter in regression analysis? [J]. Papers in Regional Science, 2005, 65(1): 11−34. DOI: 10.1111/j.1435-5597.1988.tb01155.x

    Wu J G, Jelinski D E, Luck M, et al. Multiscale analysis of landscape heterogeneity: scale variance and pattern metrics [J]. Geographic Information Sciences, 2000, 6(1): 6−19.

    Li H, Reynolds J F. On definition and quantification of heterogeneity [J]. Oikos, 1995, 73(2): 280. DOI: 10.2307/3545921

    刘畅. 黑龙江省森林碳储量空间分布研究[D]. 哈尔滨: 东北林业大学, 2014.
    李翔. 林分因子对凋落物性质空间异质性的影响[D]. 北京: 北京林业大学, 2019.
    应奎, 李旭东, 程东亚. 岩溶槽谷流域生态环境质量的遥感评定 [J]. 国土资源遥感, 2020, 32(3): 173−182.
    陈彦光. 基于moran统计量的空间自相关理论发展和方法改进 [J]. 地理研究, 2009, 28(6): 1449−1463.
    刘禹希, 王学锋, 吕少梁, 等. 南海北部海域大眼鲷空间自相关性 [J]. 水产学报, 2021, 45(8): 1361−1373.

    Anselin L. Local indicators of spatial association-LISA [J]. Geographical Analysis, 2010, 27(2): 93−115. DOI: 10.1111/j.1538-4632.1995.tb00338.x

    张维生. 黑龙江省森林空间自相关分析 [J]. 东北林业大学学报, 2008, 36(10): 16−18. DOI: 10.3969/j.issn.1000-5382.2008.10.005
    王维芳, 董薪明, 董小枫, 等. 森林生物量的空间自相关性研究 [J]. 森林工程, 2018, 34(2): 35−39. DOI: 10.3969/j.issn.1006-8023.2018.02.008
    刘磊. 基于多源数据的森林生物量与生产力估算研究 : 以江苏省森林资源调查数据为例[D]. 南京: 南京林业大学, 2010.
    欧光龙, 王俊峰, 肖义发, 等. 思茅松天然林单木生物量地理加权回归模型构建 [J]. 林业科学研究, 2014, 27(2): 213−218.
    徐美玲, 王俊峰, 胥辉, 等. 思茅松天然林空间结构与单木地上生物量分配关系 [J]. 云南大学学报(自然科学版), 2020, 42(2): 364−373.
    徐美玲, 农明川, 欧光龙. 思茅松天然林林分空间结构分析 [J]. 西南林业大学学报(自然科学), 2020, 40(3): 122−130.
    王豁然. 试论世界桉树栽培现状和我国桉树人工林发展策略 [J]. 世界林业研究, 1989, 2(3): 48−55.

    Zambrano-Monserrate M A, Ruano M A, Ormeño-Candelario V, et al. Global ecological footprint and spatial dependence between countries [J]. Journal of Environmental Management, 2020, 272: 111069. DOI: 10.1016/j.jenvman.2020.111069

    苏仕林, 周楠. 楚雄州彝族药用植物的民族植物学研究 [J]. 安徽农学通报, 2020, 26(16): 45−49. DOI: 10.3969/j.issn.1007-7731.2020.16.015
    胥辉, 张会儒. 林木生物量模型研究[M]. 昆明: 云南科学技术出版社, 2002.
    胥辉. 立木生物量模型构建及估计方法的研究[D]. 北京: 北京林业大学, 1998.
    刘帅, 张江, 李建军, 等. 森林空间结构分析中基于voronoi图的样地边缘校正 [J]. 林业科学, 2017, 53(1): 28−37. DOI: 10.11707/j.1001-7488.20170104

    Perry G L W. SpPack: spatial point pattern analysis in Excel using Visual Basic for Applications (VBA) [J]. Environmental Modelling & Software, 2004, 19(6): 559−569.

    Besag J. Spatial interaction and the statistical analysis of lattice systems [J]. Journal of the Royal Statistical Society: Series B (Methodological), 1974, 36(2): 192−225. DOI: 10.1111/j.2517-6161.1974.tb00999.x

    Moran P A P. Notes on continuous stochastic phenomena [J]. Biometrika, 1950, 37(1/2): 17−23. DOI: 10.2307/2332142

    孟斌, 王劲峰, 张文忠, 等. 基于空间分析方法的中国区域差异研究 [J]. 地理科学, 2005, 25(4): 393−400. DOI: 10.3969/j.issn.1000-0690.2005.04.002
    陈伟, 王女英, 柯文前. 顾及自然地表要素和陆地交通网络的空间可达性分析方法 [J]. 地理与地理信息科学, 2020, 36(5): 8−13. DOI: 10.3969/j.issn.1672-0504.2020.05.002

    Penttinen A. Statistics for Marked Point Patterns[M]// The Yearbook of the Finnish Statistical Society 2006, 2006: 70-91.

    张博, 欧光龙, 孙雪莲, 等. 空间效应及其回归模型在林业中的应用 [J]. 西南林业大学学报, 2016, 36(3): 144−152.
    何矣, 蒋瑞滨, 文敏. 基于landsat 8影像的汝城县森林地上生物量遥感估算研究 [J]. 湖南林业科技, 2020, 47(3): 59−64. DOI: 10.3969/j.issn.1003-5710.2020.03.010
    刘帅, 李建军, 李丹, 等. 林木空间分布格局分析方法比较及其适应性 [J]. 林业科学, 2019, 55(11): 73−84. DOI: 10.11707/j.1001-7488.20191109
    顾琪, 王平, 焦月潇, 等. 6个竹种的空间点格局分析 [J]. 安徽农业大学学报, 2019, 46(3): 432−439.
    欧光龙, 胥辉. 森林生物量模型研究综述 [J]. 西南林业大学学报(自然科学), 2020, 40(1): 1−11.
    王琰. 吕梁山森林生态系统碳密度及空间分布格局[D]. 太原: 山西大学, 2019.

    Zhang L J, Bi H Q, Cheng P F, et al. Modeling spatial variation in tree diameter-height relationships [J]. Forest Ecology and Management, 2004, 189(1/2/3): 317−329.

    Zhang L J, Gove J H, Heath L S. Spatial residual analysis of six modeling techniques [J]. Ecological Modelling, 2005, 186(2): 154−177. DOI: 10.1016/j.ecolmodel.2005.01.007

    刘畅, 胥辉, 欧光龙. 思茅松天然林单木含碳量空间异质性分析 [J]. 西南林业大学学报(自然科学), 2019, 39(4): 76−82.
图(5)  /  表(2)
计量
  • 文章访问数:  778
  • HTML全文浏览量:  312
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-02-08
  • 修回日期:  2021-08-27
  • 录用日期:  2021-11-08
  • 网络出版日期:  2021-11-25
  • 发布日期:  2022-03-19

目录

/

返回文章
返回