一种测算单株毛竹扩散能力及关键调节因子的方法技术

本发明专利技术公开了一种测算单株毛竹扩散能力及关键调节因子的方法，包括：选取有毛竹分布的林分样地；在样地内标记单株毛竹，标记满两年后设置样圆，计算单株毛竹扩散能力指标；计算样地林分结构指标：林分密度、平均胸径、平均树高、物种丰富度指数、玛格列夫丰富度指数、香农

【技术实现步骤摘要】
一种测算单株毛竹扩散能力及关键调节因子的方法


[0001]本专利技术涉及竹林经营领域，具体涉及一种测算单株毛竹扩散能力及关键调节因子的方法。

技术介绍

[0002]毛竹(Phyllostachys edulis)原产于我国亚热带地区，属于禾本科(Gramineae)竹亚科(Bambusoideae)刚竹属(Phyllostachys)，是典型的无性系植物，为单轴型散生竹类。竹养鞭、鞭生芽、芽育笋，笋成竹
……
周而复始，生生不息，逐年向四周扩散，对邻近森林群落结构组成、生物多样性、土壤理化性质等方面产生严重影响，受到中外学者高度关注并开展毛竹扩张方面的研究。然而缺乏对单株毛竹向周边扩散能力测算及关键抑制因子研究，而单株毛竹扩散是毛竹林向外扩张的源头，单株毛竹扩散能力及关键抑制因子的测算是摸清毛竹向外扩张的基础，为提出科学且有针对性的抑制技术提供依据。
[0003]关于毛竹扩散方面的研究，主要集中在毛竹扩张产生的影响，如对林地土壤肥力(吴家森，姜培坤，王祖良.天目山国家级自然保护区毛竹扩张对林地土壤肥力的影响[J].江西农业大学学报，2008，30(4)：689
‑
692)、凋落物水文功能(宋庆妮，杨清培，欧阳明，等.毛竹扩张的生态后效：凋落物水文功能评价[J].生态学杂志，2015，34(8)：2281
‑
2287)、周边森林群落菌根系统(潘璐，牟溥，白尚斌，等.毛竹林扩张对周边森林群落菌根系统的影响[J].植物生态学报，2015，39(4)：371
>‑
382)等的影响；同时，根据毛竹鞭根状态分析，提出关于抑制毛竹林扩散的研究，如蔡亮等(蔡亮，张瑞霖，李春福，等.基于竹鞭状态分析的抑制毛竹林扩散的方法.东北林业大学学报，2003，31(5)：68－70.)提出了挖沟、灌水的设想；另外，根据毛竹生物学特性和培育目的，国家林业局也发布了相关的技术规程《毛竹林丰产技术》(GB/T 20391
‑
2006)。综合上述毛竹相关研究，多是揭示毛竹扩张对其周边森林群落产生的影响，而单株毛竹扩散能力及关键调节因子的测算方法研究，未见报道。

技术实现思路

[0004]本专利技术提供了一种测算单株毛竹扩散能力及关键调节因子的方法，该方法较现有的人工开挖地下竹鞭计算单株毛竹扩散能力的方法更省力省时，效率明显提高，而且不会产生因开挖竹鞭产生对土壤结构、微生物的变化及影响。
[0005]一种测算单株毛竹扩散能力及关键调节因子的方法，包括步骤：
[0006]S1.选取有毛竹分布的林分样地；
[0007]S2.在样地内标记各单株毛竹，标记满两年以标记的单株毛竹为圆心、10米为半径设置样圆，作为圆心的单株毛竹记作中心竹株，测定样圆内新增毛竹株数及各新增毛竹分别距中心竹株的距离，计算单株毛竹扩散能力指标；
[0008]所述单株毛竹扩散能力指标的计算公式如下：
[0009][0010][0011][0012]式中，R
i
为第i株中心竹株的扩散速率，D
ij
为第i株中心竹株与第j株新增毛竹的距离，a为样地调查间隔期2年，k
i
为第i株中心竹株10米半径样圆内的新增毛竹株数，j为第i株中心竹株10米半径样圆内新增的第j株毛竹，R1为样地内n株中心竹株的平均扩散速率即毛竹平均扩散速率，R2为单株中心竹株平均每年新增的毛竹株数；
[0013]S3.确定样地林分内毛竹和/或木本植物的物种名称，测量得到样地林分内毛竹和/或木本植物的胸径和树高，并采用相邻网格法记录毛竹和/或木本植物的位置坐标(x、y)；
[0014]S4.根据所述物种名称、胸径、树高和位置坐标(x、y)测量数据，计算样地林分结构指标：林分密度、平均胸径、平均树高、物种丰富度指数、玛格列夫丰富度指数、香农
‑
维纳多样性指数、皮卢均匀度指数、树种简单混交度、多样性混交度、竞争指数和聚集指数；
[0015]S5.将R1、R2对应记作Y1、Y2归为一组变量组Y，将林分密度、平均胸径、平均树高、物种丰富度指数、玛格列夫丰富度指数、香农
‑
维纳多样性指数、皮卢均匀度指数、树种简单混交度、多样性混交度、竞争指数和聚集指数对应记作X1、X2、X3、X4、X5、X6、X7、X8、X9、X
10
、X
11
归为一组变量组X，经标准化后进行典型相关分析，得到二组典型变量V1‑
U1和V2‑
U2；其中：
[0016]V1＝
‑
3.387Y1+0.649Y2，
[0017]U1＝0.487X1+0.251X2‑
0.037X3+0.374X4‑
0.155X5+3.537X6+3.385X7‑
0.644X8+6.920X9+0.043X
10
‑
0.537X
11
；
[0018]V2＝
‑
1.181Y1+3.331Y2，
[0019]U2＝1.166X1+1.632X2‑
1.733X3‑
3.153X4‑
5.206X5‑
0.206X6‑
0.377X7‑
5.821X8+0.308X9+0.012X
10
+0.107X
11
；
[0020]S6.根据二组典型变量V1‑
U1和V2‑
U2得出单株毛竹扩散能力与关键调节因子的关系：
[0021]毛竹平均扩散速率与香农
‑
维纳多样性指数、皮卢均匀度指数、多样性混交度呈负相关；
[0022]单株中心竹株平均每年新增的毛竹株数与物种丰富度指数、玛格列夫丰富度指数、树种简单混交度呈负相关；
[0023]所述物种丰富度指数、玛格列夫丰富度指数、香农
‑
维纳多样性指数、皮卢均匀度指数、树种简单混交度、多样性混交度作为调节单株毛竹扩散能力的关键调节因子。
[0024]为了达到更好的专利技术效果，进行以下优选：
[0025]S2步骤中，新增毛竹计算原则包括：
[0026](a)当两个以上样圆交集时，交集内的新增毛竹按照就近原则计算；
[0027](b)当两个以上样圆交集时，样圆相交点上的新增毛竹，以算术平均法计算。以提高测算准确度。
[0028]进一步优选地，S2步骤中，距样地四周边界向样地内设置10米范围作为缓冲区，以
进一步提高测算准确度。
[0029]S2步骤中，新增毛竹计算原则还包括：(c)缓冲区内只统计新增竹株。
[0030]所述两组典型变量V1‑
U1和V2‑
U2在表达林分结构与单株毛竹扩散能力的相关关系的贡献率和显著性水平为V1‑
U1大于V2‑
U2。
[0031]本专利技术测算单株毛竹扩散能力及关键本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种测算单株毛竹扩散能力及关键调节因子的方法，其特征在于，包括步骤：S1.选取有毛竹分布的林分样地；S2.在样地内标记各单株毛竹，标记满两年以标记的单株毛竹为圆心、10米为半径设置样圆，作为圆心的单株毛竹记作中心竹株，测定样圆内新增毛竹株数及各新增毛竹分别距中心竹株的距离，计算单株毛竹扩散能力指标；所述单株毛竹扩散能力指标的计算公式如下：所述单株毛竹扩散能力指标的计算公式如下：所述单株毛竹扩散能力指标的计算公式如下：式中，R
i
为第i株中心竹株的扩散速率，D
ij
为第i株中心竹株与第j株新增毛竹的距离，a为样地调查间隔期两年，k
i
为第i株中心竹株10米半径样圆内的新增毛竹株数，j为第i株中心竹株10米半径样圆内新增的第j株毛竹，R1为样地内n株中心竹株的平均扩散速率即毛竹平均扩散速率，R2为单株中心竹株平均每年新增的毛竹株数；S3.确定样地林分内毛竹和/或木本植物的物种名称，测量得到样地林分内毛竹和/或木本植物的胸径和树高，并采用相邻网格法记录毛竹和/或木本植物的位置坐标(x、y)；S4.根据所述物种名称、胸径、树高和位置坐标(x、y)测量数据，计算样地林分结构指标：林分密度、平均胸径、平均树高、物种丰富度指数、玛格列夫丰富度指数、香农
‑
维纳多样性指数、皮卢均匀度指数、树种简单混交度、多样性混交度、竞争指数和聚集指数；S5.将R1、R2对应记作Y1、Y2归为一组变量组Y，将林分密度、平均胸径、平均树高、物种丰富度指数、玛格列夫丰富度指数、香农
‑
维纳多样性指数、皮卢均匀度指数、树种简单混交度、多样性混交度、竞争指数和聚集指数对应记作X1、X2、X3、X4、X5、X6、X7、X8、X9、X
10
、X
11
归为一组变量组X，经标准化后进行典型相关分析，得到二组典型变量V1‑
U1和V2‑
U2；其中：V1＝
‑
3.387Y1+0.649Y2，U1＝0.487X1+0.251X2‑
0.037X3+0...

【专利技术属性】
技术研发人员：王秀云，陈卓梅，宋绪忠，杨华，
申请(专利权)人：浙江省林业科学研究院，
类型：发明
国别省市：