长白落叶松解析木数据参数化3-PG模型

doi:10.11707/j.1001-7488.20220111

摘要/Abstract

摘要：

目的: 验证解析木数据参数化3-PG模型的可行性, 利用解析木数据标定的模型预测长白落叶松人工林生长变化, 为扩展参数化3-PG模型数据源提供依据, 为模型数据选择提供参考。方法: 以孟家岗林场长白落叶松人工林为研究对象, 基于2019年解析木数据模拟胸径连续观测数据, 结合相关公式计算叶、干、根生物量和蓄积量。根据参数敏感性分析结果, 采用直接计算、参考文献、迭代拟合和默认参数等方法对3-PG模型进行参数化, 利用密度试验林和固定样地数据对模型输出进行精度验证, 并对模型输出与地面观测值进行回归分析。结果: 验证结果表明, 模型拟合精度很高(P < 0.01, n=138), 可以较好反映解析木样地中林分生长变化。胸径、叶干生物量比、干生物量、总生物量和蓄积量的决定系数(R²)在0.95以上, 根生物量拟合精度略低(R²=0.88)。密度试验林数据(n=140)和固定样地数据(n=87)与模型输出之间具有较高相关性, 各计算量的R²在0.81~0.97之间(P < 0.01)。敏感性分析结果表明, 当胸径为20 cm时, 叶与干生物量分配比(pFS20)及初级生产力分配给根最大值(pRx)2个参数均对胸径、根、干生物量和蓄积量显示出较高敏感性。结论: 解析木数据参数化3-PG模型的拟合精度和预测精度均较高, 研究结果扩展了参数化3-PG模型的可用数据源, 可为利用3-PG模型模拟长白落叶松人工林生长变化提供新的依据。

关键词: 3-PG模型, 解析木数据, 参数化, 长白落叶松

Abstract:

Objective: In order to verify the usability of destructive sampling data in 3-PG(physiological principles in predicting growth) model parameterization, this study predicted the growth of Larix olgensis using the 3-PG model calibrated by destructive sampling data, which can help to expand the data source of 3-PG model parameterization and provide a reference for model data selection. Method: Taking L. olgensis plantation in Mengjiagang forest farm as the study object, the destructive sampling data collected in 2019 was calculated to simulate the continuous observation data. The biomass of foliage, stem, and root, and stand volume were calculated based on the relevant biomass formulas. According to the sensitivity analysis results, the 3-PG model was parameterized by direct calculation, reference, and iterative fitting. Then the modeling results were verified by continuous observation data and fixed plot data, and the model output was regressed with observed data. Result: The results showed that the fitting accuracy was very high, and the model output data could well reflect the growth of the stand (n=138). The determination coefficient (R²) values of DBH, the ratio of foliage to stem biomass, stem biomass, total biomass and volume were all above 0.95 (P < 0.01). In contrast, the root biomass had a relatively lower fitting accuracy (R²=0.88). The model's predicting accuracy was verified by continuous observation data and fixed plot data. High correlations were observed between the model output and continuous observation data (n=140) and fixed plot data (n=87). The R² of all values ranged from 0.81 to 0.97 (P < 0.01). According to the sensitivity analysis results, foliage∶stem partitioning ratio at DBH=20 cm (pFS20) and Maximum fraction of NPP to roots (pRx) showed high sensitivity to DBH, root biomass, stem biomass, and volume. Conclusion: 3-PG model calibrated by destructive sampling data had a relatively high fitting accuracy and prediction accuracy. The results of this study would expand the data source of 3-PG model parameterization and provide a new basis for the 3-PG model to simulate the growth of L. olgensis plantation.

Key words: 3-PG(physiological principles in predicting growth) model, destructive sampling data, parameterization, Larix olgensis

中图分类号:

白羽,庞勇,夏晓运,贾炜玮. 长白落叶松解析木数据参数化3-PG模型[J]. 林业科学, 2022, 58(1): 98-110.

Yu Bai,Yong Pang,Xiaoyun Xia,Weiwei Jia. 3-PG Model Parameterization Using Destructive Sampling Data of Larix olgensis[J]. Scientia Silvae Sinicae, 2022, 58(1): 98-110.

图/表 12

图1

表1

表2

表3

图2

表4

表5

图3

表6

图4

图5

图6

参考文献 0

	董利虎, 李凤日, 贾炜玮, 等. 2019. 长白落叶松人工林立木生物量及含碳量计量方法: DB23/T 2333—2019. 哈尔滨: 黑龙江省市场监督管理局.
	Dong L H, Li F R, Jia W W, et al. 2019. Measurement method of standing biomass and carbon content of Larix olgensis plantation: DB23/T 2333—2019. Harbin: Administration for Market Regulation of Heilongjiang Province. [in Chinese]
	郭明春, 刘建锋. 杉木人工林林分生长模拟. 福建林学院学报, 2011, 31 (1): 65- 68. doi: 10.3969/j.issn.1001-389X.2011.01.013
	Guo M C , Liu J F . Modeling the stand growth of Chinese fir plantations. Journal of Fujian College of Forestry, 2011, 31 (1): 65- 68. doi: 10.3969/j.issn.1001-389X.2011.01.013
	凌向明, 陈昭. 浅析森林经营管理对森林碳汇的作用和影响. 低碳世界, 2019, 9 (3): 298- 300. doi: 10.3969/j.issn.2095-2066.2019.03.180
	Ling X M , Chen Z . The analysis of forest management function and influence on forest carbon sink. Low-Carbon World, 2019, 9 (3): 298- 300. doi: 10.3969/j.issn.2095-2066.2019.03.180
	马华文, 许翠清, 李海朝. 立地条件对长白落叶松光合特性的影响. 东北林业大学学报, 2008, 36 (8): 4- 7. doi: 10.3969/j.issn.1000-5382.2008.08.002
	Ma H W , Xu C Q , Li H C . Photosynthesis characteristics of Larix olgensis under different site conditions. Journal of Northeast Forestry University, 2008, 36 (8): 4- 7. doi: 10.3969/j.issn.1000-5382.2008.08.002
	孟宪宇. 测树学. 北京: 中国林业出版社, 2006.
	Meng X Y . Forest measurement. Beijing: China Forestry Publishing House, 2006.
	孙志虎. 2005. 长白落叶松人工用材林长期生产力维持的研究. 哈尔滨: 东北林业大学.
	Sun Z H. 2005. Study on long-term productivity maintenance of artificial timber stand of Larix olgensis. Harbin: Northeast Forestry University. [in Chinese]
	夏晓运, 庞勇, 黄庆丰, 等. 基于3-PG模型的长白落叶松生物量生长预测. 林业科学, 2021, 57 (3): 67- 78.
	Xia X Y , Pang Y , Huang Q F , et al. Prediction of biomass growth Larix olgensis based on 3-PG model. Scientia Silvae Sinicae, 2021, 57 (3): 67- 78.
	解雅麟, 王海燕, 雷相东. 基于过程模型的气候变化对长白落叶松人工林净初级生产力的影响. 植物生态学报, 2017, 41 (8): 826- 839.
	Xie Y L , Wang H Y , Lei X D . Effects of climate change on net primary productivity in Larix olgensis plantations based on process modeling. Chinese Journal of Plant Ecology, 2017, 41 (8): 826- 839.
	许晨璐, 孙晓梅, 张守攻. 日本落叶松与长白落叶松及其杂种光合特性比较. 北京林业大学学报, 2012, 34 (4): 62- 66.
	Xu C L , Sun X M , Zhang S G . Comparison in photosynthetic characteristics of Larix kaempferi, L. olgensis and their hybrids. Journal of Northeast Forestry University, 2012, 36 (8): 4- 7.
	赵梅芳, 项文化, 田大伦, 等. 基于3-PG模型的湖南会同杉木人工林蒸发散估算. 湖南农业科学, 2008, (3): 158- 162. doi: 10.3969/j.issn.1006-060X.2008.03.057
	Zhao M F , Xiang W H , Tian D L , et al. Simulating evapotranspiration of Chinese fir plantation in Huitong, Hunan using the 3-PG. Hunan Agricultural Sciences, 2008, (3): 158- 162. doi: 10.3969/j.issn.1006-060X.2008.03.057
	Coops N C , Waring R H , Landsberg J J . Assessing forest productivity in Australia and New Zealand using a physiologically-based model driven with averaged monthly weather data and satellite-derived estimates of canopy photosynthetic capacity. Forest Ecology and Management, 1998, 104 (1/3): 113- 127.
	Coops N C , Waring R H , Beverly E L . Assessing the past and future distribution and productivity of ponderosa pine in the Pacific northwest using a process model, 3-PG. Ecological Modelling, 2004, 183 (1): 107- 124.
	Dye P J , Jacobs S , Drew D . Verification of 3-PG growth and water-use predictions in twelve Eucalyptus plantation stands in Zululand, South Africa. Forest Ecology & Management, 2004, 193 (1/2): 197- 218.
	Esprey L J , Sands P J , Smith C W . Understanding 3-PG using a sensitivity analysis. Forest Ecology and Management, 2004, 193, 235- 250. doi: 10.1016/j.foreco.2004.01.032
	Forrester D I , Tang X . Analysing the spatial and temporal dynamics of species interactions in mixed-species forests and the effects of stand density using the 3-PG model. Ecological Modelling, 2015, 319, 233- 254.
	Gonzalez-Benecke C A , Jokela E J , Cropper W P , et al. Parameterization of the 3-PG model for Pinus elliottii stands using alternative methods to estimate fertility rating, biomass partitioning and canopy closure. Forest Ecology & Management, 2014, 327, 55- 75.
	Gupta R , Sharma L K . The process-based forest growth model 3-PG for use in forest management: a review. Ecological Modelling, 2019, 397, 55- 73. doi: 10.1016/j.ecolmodel.2019.01.007
	Landsberg J J , Waring R H . A generalised model of forest productivity using simplified concepts of radiation-use efficiency, carbon balance and partitoning. Forest Ecology and Management, 1997, 95 (3): 209- 228. doi: 10.1016/S0378-1127(97)00026-1
	Landsberg J J , Waring R H , Coops N C . Performance of the forest productivity model 3-PG applied to a wide range of forest types. Forest Ecology and Management, 2003, 172 (2): 199- 214.
	López-Serrano F R , Martínez-García E , Dadi T , et al. Biomass growth simulations in a natural mixed forest stand under different thinning intensities by 3-PG process-based model. European Journal of Forest Research, 2015, 134 (1): 167- 185. doi: 10.1007/s10342-014-0841-3
	Rodríguez R , Espinosa M , Real P , et al. Analysis of productivity of radiata pine plantations under different silvicultural regimes using the 3-PG process-based model. Australian Forestry, 2002, 65 (3): 165- 172. doi: 10.1080/00049158.2002.10674869
	Sands P J . Adaptation of 3-PG to novel species: guidelines for data collection ad parameter assignment. Cooperative Research Centre for Sustainable Production Forestry and CSIRO. Technical Report, 2004, 141, 34.
	Sands P J , Landsberg J J . Parameterisation of 3-PG for plantation grown Eucalyptus globulus. Forest Ecology and Management, 2002, 163, 273- 292. doi: 10.1016/S0378-1127(01)00586-2
	Santosh S , Thomas F , Randolph W . Determination of fertility rating (FR) in the 3-PG model for loblolly pine plantations in the southeastern United States based on site index. Forests, 2015, 33 (6): 256- 262.
	Song X , Bryan B A , Almeida A C , et al. Time-dependent sensitivity of a process-based ecological model. Ecological Modelling, 2013, 65, 114- 123.
	Wei L , John M , Zhang J W , et al. 3-PG simulations of young ponderosa pine plantations under varied management intensity: why do they grow so differently?. Forest Ecology and Management, 2014, 313, 69- 82. doi: 10.1016/j.foreco.2013.10.035
	Xie Y L , Wang H Y , Lei X D . Application of the 3-PG model to predict growth of Larix olgensis plantations in northeastern China. Forest Ecology and Management, 2017, 406, 208- 218. doi: 10.1016/j.foreco.2017.10.018

样地号 Plot ID	林班号 Compartment ID	小班号 Sub-compartmentID	树龄 Age/a	面积 Area/hm²	立地指数 Site index	株数 Stem number/hm^-2	平均胸径 Mean DBH/cm	平均树高 Meanheight/m	蓄积量 Volume/(m³·hm^-2)
LYS720	64	11	52	0.06	21.48	642	24.82	26.67	337.62
LYS03	68	13	37	0.06	23.06	1 017	19.81	23.92	370.47
LYS04	81	2	36	0.06	26.22	1 100	18.59	22.35	324.16
LYS13	70	5	22	0.06	21.94	1 167	13.21	15.65	191.64
LYS14	70	5	35	0.06	21.85	1 467	13.82	16.02	198.77
LYS15	70	5	21	0.06	22.32	1 167	14.37	16.20	153.40

类型Type	参数Parameters	参数值Value	来源Source
异速生长关系与分配 Allometric relationships and partitioning	胸径2 cm时叶与干生物量分配比 Foliage∶stem partitioning ratio at DBH = 2 cm	0.8	拟合 Fitting
	胸径20 cm时叶与干生物量分配比Foliage∶stem partitioning ratio at DBH = 20 cm	0.6	拟合 Fitting
	干生物量与胸径关系中常数 Constant in the stem biomass v. diam. relationship	0.032 52	董利虎等(2019)
	干生物量与胸径关系中指数 Power in the stem biomass v. diam. relationship	2.742 6	董利虎等(2019)
	初级生产力分配给根最大值 Maximum fraction of NPP to roots	0.5	拟合 Fitting
	初级生产力分配给根最小值 Minimum fraction of NPP to roots	0.15	拟合 Fitting
凋落物与根更新 Litterfall and root turnover	最大凋落物速率 Maximum litterfall rate	0.06	拟合 Fitting
凋落物与根更新 Litterfall and root turnover	根月平均凋落速率 Mean monthly root turnover rate	0.005	拟合 Fitting
生长调节 Growth modifier	生长最低温度 Minimum temperature for growth	-2	许晨露等(2012)
	生长最适温度 Optimum temperature for growth	17	孙志虎(2005)
	生长最高温度 Maximum temperature for growth	38	许晨露等(2012)
	植物生长最大林龄 Maximum stand age used in age modifier	100	拟合 Fitting
树干死亡与自疏 Stem mortality and self-thinning	最大林龄的死亡率 Mortality rate for large tree age	0.6	计算 Calculation
	在1 000株·hm^-2每株树最大树干质量 Max. stem mass per tree @ 1 000 trees·hm^-2	240	计算 Calculation
	林分郁闭年龄 Age at full canopy cover	15	拟合 Fitting
冠层电导 Conductance	冠层量子效率 Canopy quantum efficiency	0.05	拟合 Fitting

数据类型 Data type	数据来源 Data source	数据获取时间 Data acquisition time	数据特点 Data characteristics
解析木数据 Destructive sampling data	6块解析木样地 Six destructive sampling plots	2019年 Observed in 2019	模拟连年观测数据，用于标定模型 Used to simulate continuous observation data and calibrate model
密度林数据 Continuous observation data	5块密度林样地 5 continuous observation plots	1975—2020年 Observed from 1975 to 2020	覆盖林龄较广，检验大时间跨度预测精度 Covered a wide range of stand age, used to verify the prediction accuracy of large time span
固定样地数据 Fixed plot data	24块固定样地 24 fixed plots	2012、2014、2016和2020年 Observed in 2012，2014，2016 and 2020	覆盖地理范围广，检验大区域范围预测精度 Covered a wide spatial range, used to verify the prediction accuracy of large area

统计量 Statistic value	R²	斜率 Slope	截距 Intercept	ME	MAE	RMSE	MRE(%)
胸径DBH/cm	0.96	0.98	0.19	-0.12	0.57	0.27	-0.81
叶干生物量比Foliage∶stem bionass ratio	0.91	1.20	-0.01	0.00	0.01	0.05	-7.08
干生物量Stem biomass/(t·hm^-2)	0.93	1.10	1.29	10.63	11.78	2.51	11.01
根生物量Root biomass/(t·hm^-2)	0.83	0.75	10.36	3.38	4.21	1.41	12.29
总生物量Total biomass/(t·hm^-2)	0.93	0.90	15.48	1.49	9.43	0.94	1.06
蓄积量Volume/(m³·hm^-2)	0.93	1.12	-38.27	-11.48	19.20	2.61	-5.15

统计量 Statistic value	R²	斜率break/>Slope	截距break/>Intercept	ME	MAE	RMSE	MRE(%)
胸径DBH/cm	0.97	0.97	0.16	-0.29	0.43	0.53	-1.93
叶干生物量比Foliage∶stem bionass ratio	0.81	0.98	-0.02	-0.02	0.02	0.14	-36.56
干生物量Stem biomass/(t·hm^-2)	0.95	1.22	-7.73	13.91	14.28	3.64	14.00
根生物量Root biomass/(t·hm^-2)	0.85	0.71	8.43	0.38	3.28	0.60	1.36
总生物量Total biomass/(t·hm^-2)	0.94	0.98	0.71	-2.30	8.92	1.48	-1.57
蓄积量Volume/(m³·hm^-2)	0.94	1.20	-63.79	-15.67	23.99	3.86	-6.47