云南省主要针叶树种树高曲线模型研建

doi:10.13466/j.cnki.lyzygl.2019.04.007

摘要/Abstract

摘要：

树高曲线在森林经营和收获调整中有着重要的作用,建立适用于云南省主要针叶树种的树高曲线模型,为其森林经营提供参考。基于云南省两期森林资源连续清查数据,以冷杉、思茅松、云南松、华山松和杉木作为研究对象,选用Richards等15种树高曲线模型作为备选模型,以决定系数(R ²)、均方根误差(RMSE)、平均绝对误差(MAE)、总体相对误差(TRE)和平均预估误差(MPE)作为模型拟合效果的评价指标。结果表明:Hossfeld方程能较好地描述冷杉和杉木的树高曲线,双曲线方程能较好地描述云南松和华山松的树高曲线,Wykoff方程能较好地描述思茅松的树高曲线;云南松最优树高曲线模型决定系数为0.676,其它针叶树种决定系数均大于0.710。独立样本数据检验表明,各树种最优树高曲线模型均有较好的适用性,对云南省主要针叶树种树高有较好的预测效果。

关键词: 云南省, 树高曲线, 针叶树种, 模型构建

Abstract:

Tree height curve plays an important role in the forest management and yield regulation,the tree height curve model for main conifer species in Yunnan Province was established to provide reference for its forest management.Based on two periods of national forest inventory data in Yunnan Province,this paper takes Abies fabri,Pinus kesiya var.langbianensis,Pinus yunnanensis,Pinus armandii,Cunninghamia lanceolata as research objects,Richards etc.fifteen tree height curve models were selected as alternative models,the coefficient of determination,root mean square error,mean absolute error,total relative error and mean predict error were used as the evaluation indicators of the pros and cons of the model.The results show that Hossfeld can better describe the tree height curve of Abies fabri and Cunninghamia lanceolata,hyperbola can better describe the tree height curve of Pinus yunnanensis and Pinus armandii,Wykoff can better describe the tree height curve of Pinus kesiya var.langbianensis.The coefficient of determination of the optimal tree height curve model of Pinus yunnanensis is 0.676,and the coefficient of determination of other conifer tree specieses are greater than 0.710,the independent sample data test shows that the optimal tree height curve model of each tree species has good applicability,and has a good prediction effect on the tree height of the main conifer tree species in Yunnan Province.

Key words: Yunnan Province, tree height curve, conifer species

中图分类号:

S758
S791

覃阳平, 李华, 李永亮, 张成程, 王蒙. 云南省主要针叶树种树高曲线模型研建[J]. 林业资源管理, 2019,(4): 46-51.

QIN Yangping, LI Hua, LI Yongliang, ZHANG Chengcheng, WANG Meng. Developing Models for Tree Height Curves of Main Conifer Species in Yunnan Province[J]. FOREST RESOURCES WANAGEMENT, 2019,(4): 46-51.

图/表 4

表1

表2

树高曲线模型

模型编号	方程名称	公式
1	幂函数	$H = 1.3 + a D b$
2	Wykoff等(1982)	$H = 1.3 + e (a + b D + 1)$
3	Bates and Watts(1980)	$H = 1.3 + aD b + D$
4	Loetsh等(1973)	$H = 1.3 + D 2 (a + bD) 2$
5	Schumacher(1939)	$H = 1.3 + a e - b D$
6	Mitscherlich(1919)	$H = 1.3 + a (1 - b e - cD)$
7	Ratkowsky(1990)	$H = 1.3 + a e - b D + c$
8	Logistic (1838)	$H = 1.3 + a 1 + b e - cD$
9	修正Weibull(Yang,1978)	$H = 1.3 + a (1 - e - b D c)$
10	Richards(1959)	$H = 1.3 + a (1 - e - cD) b$
11	Gompertz(1825)	$H = 1.3 + a e - b e - cD$
12	双曲线	$H = a - b D + c$
13	柯列尔(Rоляср,1878)	$H = 1.3 + a D b e - cD$
14	Hossfeld (1822)	$H = 1.3 + a (1 + b D - c)$
15	Korf (1939)	$H = 1.3 + a e - b D - c$

表2

表3

云南省主要针叶树种树高曲线模型的误差分析

树种	评价指标		模型编号
树种	评价指标		1		2		3		4		5		6		7		8		9		10		11		12		13		14	15
冷杉	R²	0.694		0.752		0.724		0.744		0.753		0.751		0.753		0.746		0.751		0.753		0.752		0.749		0.744		0.755		0.753
	RMSE/m	4.67		4.20		4.44		4.27		4.20		4.23		4.21		4.27		4.23		4.21		4.22		4.24		4.28		4.20		4.21
	MAE/m	3.82		3.22		3.57		3.32		3.24		3.26		3.24		3.30		3.25		3.24		3.24		3.27		3.34		3.23		3.25
	TRE/%	-0.88		0.01		-1.13		-0.55		0.13		0.00		0.11		-0.21		-0.16		0.01		-0.05		0.00		-0.30		0.07		0.20
	MPE/%	3.70		3.33		3.51		3.38		3.33		3.35		3.34		3.38		3.35		3.34		3.34		3.36		3.39		3.32		3.34
思茅松	R²	0.693		0.712		0.696		0.708		0.711		0.711		0.712		0.709		0.712		0.712		0.711		0.711		0.712		0.712		0.711
	RMSE/m	2.22		2.15		2.21		2.16		2.15		2.15		2.15		2.16		2.15		2.15		2.15		2.15		2.15		2.15		2.15
	MAE/m	1.81		1.71		1.80		1.75		1.71		1.72		1.72		1.74		1.72		1.72		1.72		1.72		1.72		1.72		1.72
	TRE/%	-0.26		0.02		-0.43		-0.17		0.07		0.00		0.01		-0.07		-0.02		-0.01		-0.02		0.00		-0.02		-0.01		0.01
	MPE/%	1.50		1.45		1.49		1.46		1.45		1.45		1.45		1.46		1.45		1.45		1.45		1.45		1.45		1.45		1.45
云南松	R²	0.665		0.674		0.665		0.676		0.671		0.676		0.676		0.673		0.676		0.676		0.676		0.676		0.676		0.676		0.676
	RMSE/m	2.16		2.13		2.16		2.12		2.14		2.12		2.12		2.14		2.12		2.12		2.13		2.12		2.12		2.12		2.12
	MAE/m	1.69		1.65		1.69		1.65		1.66		1.65		1.65		1.66		1.65		1.65		1.65		1.65		1.65		1.65		1.65
	TRE/%	-0.34		0.21		-0.61		-0.09		0.31		0.00		0.00		-0.11		-0.03		0.00		-0.03		0.00		-0.02		0.00		0.02
	MPE/%	0.86		0.85		0.86		0.85		0.86		0.85		0.85		0.85		0.85		0.85		0.85		0.85		0.85		0.85		0.85
华山松	R²	0.748		0.735		0.750		0.746		0.730		0.751		0.747		0.728		0.750		0.750		0.739		0.752		0.750		0.750		0.752
	RMSE/m	1.81		1.85		1.80		1.81		1.87		1.80		1.81		1.88		1.80		1.80		1.84		1.80		1.80		1.80		1.80
	MAE/m	1.43		1.45		1.43		1.42		1.46		1.41		1.43		1.50		1.43		1.42		1.46		1.41		1.43		1.42		1.42
	TRE/%	-0.32		0.31		-0.35		0.19		0.40		0.00		-0.09		-0.36		-0.21		-0.19		-0.20		0.00		-0.20		-0.19		-0.10
	MPE/%	2.30		2.36		2.29		2.31		2.38		2.29		2.30		2.39		2.29		2.29		2.34		2.28		2.29		2.29		2.28
杉木	R²	0.690		0.750		0.703		0.734		0.754		0.757		0.758		0.755		0.756		0.758		0.758		0.756		0.749		0.759		0.759
	RMSE/m	1.90		1.71		1.86		1.76		1.70		1.69		1.69		1.70		1.69		1.68		1.69		1.69		1.72		1.68		1.68
	MAE/m	1.50		1.30		1.48		1.35		1.29		1.28		1.28		1.28		1.28		1.28		1.28		1.29		1.29		1.28		1.28
	TRE/%	-0.71		-0.30		-1.13		-0.68		-0.20		0.00		0.03		-0.09		-0.09		0.00		-0.01		0.00		-0.20		0.01		0.05
	MPE/%	1.97		1.77		1.93		1.83		1.76		1.75		1.75		1.76		1.75		1.75		1.75		1.75		1.78		1.74		1.74

表3

表4

参考文献 28

[1]	Temesgen H, Hann D,J, Monleon D,J . Regional Height-Diameter Equations for Major Tree Species of Southwest Oregon[J]. 2008,22:213-219.
[2]	张鹏, 王新杰, 许昊 . 将乐地区马尾松标准树高曲线的研究[J]. 中南林业科技大学学报, 2015,35(3):69-73.
[3]	赵俊卉, 亢新刚, 刘燕 . 长白山主要针叶树种最优树高曲线研究[J]. 北京林业大学学报, 2009,31(4):13-18.
[4]	孟宪宇 . 测树学[M]. 3版. 北京: 中国林业出版社, 2006.
[5]	Larjavaara H C, Muller-Landau M . Measuring tree height:A quantitative comparison of two common field methods in a moist tropical forest[J]. Methods in Ecology and Evolution, 2013,21(5):9.
[6]	Adame P, Del RíO M, Cañellas I . A mixed nonlinear height-diameter model for pyrenean oak(Quercus pyrenaica Willd.)[J]. Forest Ecology and Management, 2008,256(1):88-98.
[7]	代忠迪, 姜立春 . 大兴安岭不同生态区域兴安落叶松树高曲线的研究[J]. 植物研究, 2015,35(4):583-589.
[8]	段光爽, 李学东, 冯岩 , 等. 华北落叶松天然次生林树高曲线的混合效应模型[J]. 南京林业大学学报:自然科学版, 2018,42(2):163-169
[9]	刘鑫, 王海燕, 雷相东 , 等. 基于BP神经网络的天然云冷杉针阔混交林标准树高-胸径模型[J]. 林业科学研究, 2017,30(3):368-375.
[10]	赵俊卉, 亢新刚, 张慧东 , 等. 长白山3个主要针叶树种的标准树高曲线[J]. 林业科学, 2010,46(10):191-194.
[11]	陈飞, 王健敏, 孙宝刚 , 等. 云南松的地理分布与气候关系[J]. 林业科学研究, 2012,25(2):163-168.
[12]	孙雪莲, 舒清态, 欧光龙 , 等. 基于随机森林回归模型的思茅松人工林生物量遥感估测[J]. 林业资源管理, 2015(1):71-76.
[13]	汪家社 . 杉木生态系统生物量与固碳能力的分析与评价[J]. 福建林业科技, 2008(2):1-4.
[14]	李元玖, 陈奇伯, 熊好琴 , 等. 滇中高原华山松人工林碳储量及固碳释氧效益[J]. 中南林业科技大学学报, 2015,35(2):79-84.
[15]	李艳丽, 杨华, 亢新刚 , 等. 长白山云冷杉种群结构和动态分析[J]. 北京林业大学学报, 2014,36(3):18-25.
[16]	王明亮, 唐守正 . 标准树高曲线的研制[J]. 林业科学研究, 1997(3):36-41.
[17]	胥辉, 全宏波, 王斌 . 思茅松标准树高曲线的研究[J]. 西南林学院学报, 2000(2):74-77.
[18]	王冬至, 张冬燕, 张志东 , 等. 基于非线性混合模型的针阔混交林树高与胸径关系[J]. 林业科学, 2016,52(1):30-36.
[19]	樊伟, 许崇华, 崔珺 , 等. 基于混合效应的大别山地区杉木树高-胸径模型比较[J]. 应用生态学报, 2017,28(9):2831-2839.
[20]	张懋功 . 云南年鉴[M]. 昆明: 云南年鉴社, 2017.
[21]	Temesgen H, Gadow K V . Generalized height-diameter models—an application for major tree species in complex stands of interior British Columbia[J]. European Journal of Forest Research, 2004,123(1):45-51.
[22]	马武, 雷相东, 徐光 , 等. 蒙古栎天然林单木生长模型的研究——Ⅱ.树高-胸径模型[J]. 西北农林科技大学学报:自然科学版, 2015,43(3):83-90.
[23]	曾伟生, 唐守正 . 立木生物量方程的优度评价和精度分析[J]. 林业科学, 2011,47(11):106-113.
[24]	董云飞, 孙玉军, 王轶夫 , 等. 基于BP神经网络的杉木标准树高曲线[J]. 东北林业大学学报, 2014,42(7):154-156.
[25]	董云飞, 孙玉军, 许昊 , 等. 基于非线性混合模型的杉木标准树高曲线[J]. 东北林业大学学报, 2014,42(11):72-76.
[26]	杜志, 甘世书 . 基于BP神经网络的杉木和马尾松树高曲线模型研究[J]. 中南林业调查规划, 2017,36(4):36-39.
[27]	高东启, 邓华锋, 程志楚 , 等. 基于度量误差模型方法建立的林分相容性树高曲线方程组[J]. 西北农林科技大学学报:自然科学版, 2015,43(5):65-70.
[28]	丁贵杰 . 马尾松人工林标准树高曲线模型的研究[J]. 浙江林学院学报, 1997(3):15-20.

树种	数据类型		样本数量		胸径/cm									树高/m
树种	数据类型		样本数量		平均值		最大值		最小值			标准差		平均值		最大值	最小值		标准差
冷杉	建模数据	172		32.7		101.8		6.5		17.4	19.0		36.0		4.0			8.4
冷杉	检验数据	53		30.7		81.3		7.0		15.6	18.2		32.5		5.9			7.8
思茅松	建模数据	563		15.9		33.8		5.6		4.9	12.3		23.5		3.3			4.0
思茅松	检验数据	146		16.1		34.5		5.1		5.5	12.4		23.0		3.8			4.2
云南松	建模数据	2772		13.6		36.3		5.0		4.9	9.3		26.6		2.2			3.7
云南松	检验数据	640		13.8		34.3		5.0		4.9	9.4		24.6		2.3			3.8
华山松	建模数据	329		13.0		41.1		5.4		5.9	8.5		22.5		2.8			3.6
华山松	检验数据	97		13.1		41.1		5.6		6.5	8.3		22.0		2.8			3.6
杉木	建模数据	474		11.7		35.2		5.1		4.9	8.7		19.6		1.9			3.4
杉木	检验数据	125		11.1		29.6		5.3		4.3	8.4		20.5		2.4			3.5

树种	模型编号	模型参数			评价指标
树种	模型编号	a	b	c	R²	RMSE/m	MAE/m	TRE/%	MPE/%
冷杉	14	31.22	606.11	2.01	0.754	3.95	3.01	-1.49	6.00
思茅松	2	3.47	-17.53		0.707	2.27	1.83	1.19	2.99
云南松	12	50.69	2522.02	47.72	0.671	2.21	1.70	-0.46	1.82
华山松	12	61.73	4922.05	79.88	0.745	1.83	1.40	-2.77	4.44
杉木	14	15.07	403.39	2.51	0.758	1.72	1.31	0.06	3.62