全 文 :林业科学研究 2015,28(4):551 556
ForestResearch
文章编号:10011498(2015)04055106
基于单木过程模型的杉木林分生长模拟研究
马莉燕1,张怀清1,李永亮1,郭明春2,张晔?1
(1.中国林业科学研究院资源信息研究所,北京 100091;2.河北科技师范学院园艺科技学院,河北 秦皇岛 066004)
收稿日期:20141126
项目基金:“十二五”国家863计划课题(2012AA102002);国家自然科学基金(31170590)。
作者简介:马莉燕(1989—),女,在读硕士,主要研究方向:林分结构与生长可视化模拟。
通讯作者:研究员,硕士生导师,主要从事林业可视化模拟技术研究。
摘要:基于林分尺度的3PG模型原理,在考虑单木林分生理生态过程的基础上,改进了生物量积累和分配子模型,
建立适应于单木生长的生理生态过程模型,实现基于单木过程模型的林分生长预测。选择胸径和树高生长量2个
指标与实测值对比分析,结果表明:改进模型对胸径和树高生长量的模拟平均精度分别为95.57%和91.16%,均明
显高于3PG模型对应的89.4%和85.18%的平均模拟精度,改进模型的模拟精度显著提高,模拟效果更好。
关键词:3PG模型;单木生长模拟;杉木林分生长
中图分类号:S791.27 文献标识码:A
GrowthSimulationofChineseFirBasedonIndividualTreeProcessModel
MALiyan1,ZHANGHuaiqing1,LIYongliang1,GUOMingchun2,ZHANGYecheng1
(1.ResearchInstituteofForestResourceInformationTechniques,ChineseAcademyofForestry,Beijing 100091,China;
2.ColegeofHorticultureScienceandTechnology,HebeiNormalUniversityofScienceandTechnology,Qinhuangdao 066004,Hebei,China)
Abstract:Thebiomassproductionmodelanddistributionmodelwereimprovedbasedon3PGmodeltheory,soas
theindividualtreegrowthmodelcouldbeestablishedtoachievethegrowthforecastoftheChinesefir(Cunning
hamialanceolata)stand.Inthisstudy,thediameteratbreastheight(DBH)andheight(H)wereusedtoevaluate
andcomparetheresultsofthemodelsimulationandfielddata.ItshowedthattheaverageanalogprecisionsofDBH
andHoftheimprovedmodelwererespectively95.57% and91.16%,significantlyhigherthanthatofthe3PG
model,whoseprecisionswere89.4% and85.18%.Thesimulationaccuracyoftheimprovedsimulationmodelwas
significantlyhigherandthesimulationresultswerebeter.
Keywords:3PGmodel;individualtreegrowthmodel;Chinesefir;standgrowth
林分动态的三维图形模拟可将林分调查所得的
数据以立体图像的形式反映出来,既能真实的反映
实际林分的生长状况,又能满足人们的视觉要求[1]。
完善的模型能为合理地经营森林提供理论基础、操
作依据和量化信息,并对林分结构与生长的可视化
表达具有重要意义。3PG模型是由 Landsberg和
Waring在1997年开发的一个简单的林分水平上的
过程模型,充分考虑了气候、立地条件、肥力和林龄
等多种因素对生理过程的影响,以此来模拟林分的
动态生长过程[2]。该模型因参数录入量小且易于获
取,现已应用于多个国家和地区不同树种的研究
中[3]。Almeida等运用3PG模型作为管理位于巴西
的大面积速生桉树(Eucaluptusspp.)人工林的决策
支持工具[4]。Coops等开发了支持遥感数据的模拟
软件3PGS,利用遥感卫星图像作为输入信息来预
测位于澳大利亚及新西兰的样地特征,预测精度较
高[5]。Olga利用3PG模型结合 GIS来预测美国太
平洋西北地区杨树(Populusspp.)在不同管理和气
候条件下的生物质产量[6]。以上应用充分体现了3
PG模型预测大面积林分动态的能力。最早的3PG
林 业 科 学 研 究 第28卷
模型主要用于蓝桉(EucalyptusglobulesLabil.)和巨
桉(EucalyptusgrandisHilexMaiden)的生产力预
测[7],现在可以针对特定树种将模型参数化,使其应
用越来越广泛。但是,3PG模型在一些应用中仍存
在局限性:3PG模型作为一个林分尺度的生长预测
模型,并不能很好地反映每株林木的实际生长;该模
型一般在预测内部结构无明显差异的纯人工林生长
上表现突出,但针对比较复杂的原生态林分模拟效
果不佳。本研究采用3PG建模原理,建立一套充分
考虑林木生理和生态环境的过程模型,对林分单木
与环境交互作用下的生理生长过程进行分析,实现
林分单木的形态与生长过程模拟,通过统计林分单
木模拟结果就可推出整体林分的结构与生长过程
模拟。
1 试验材料
研究区位于湖南省株洲市攸县黄丰桥国有林
场,130°04′ 113°43′E,27.04′ 27.06′N,地貌以
中低山为主,海拔在115 1270m之间。该区属中
亚热带季风湿润气候,光热水资源丰富,年平均气温
17.8℃,1月平均气温4.2℃,7月平均气温26.5℃,
年平均日照时数 1612h,年平均降水量 1410.8
mm,降水主要集中在春季,占全年降水量的40%,年
平均无霜期265d,土壤主要为低山丘陵红壤、黄壤,
坡度一般在20° 25°之间[8]。林场的森林覆盖率为
86.24%,树种主要以杉木(Cunninghamialanceolata
(Lamb.)Hook.)、松类(Pinusspp.)为主,其中杉木
面积3197.6hm2,蓄积593738m3,全部为人工造林。
本文选用位于黄丰桥林场柏市分场的1999年杉木
人工林(30m×30m)固定样地,样地林木株数为
120株,以2008年林木调查数据作为数据源,林分
环境因子采集湖南攸县气象站2008年气象数据,其
中包括月平均温度、月平均最高温度、月平均最低温
度、太阳总辐射和月平均降水量等,模型验证数据取
自2013年该样地的单木实测数据。
2 研究方法———模型改进
3PG模型包括 4个子模型:生物量积累子模
型、生物量分配子模型、林木凋亡子模型和土壤水分
平衡子模型,本研究在考虑单木林分的生理生态过
程的基础上,充分考虑单木生理特征以及林分内部
结构对单木生长的综合影响,改进其中的生物量积
累子模型和生物量分配子模型,并通过数据反复模
拟效果进行模型优化,最终建立单木生理生态生长
模型。具体改进过程如下,改进模型与3PG模型对
比情况如表1所示。
表1 改进模型与3PG模型对比
子模型 3PG模型 改进模型 改进依据
光合有效辐射APAR: APAR=RAD(1-e-kLi/δ)δ APAR=RAD(1-e-kLi)e
cw
cwi
(f(o)-1)
考虑了空间结构对单木林分获
取光照养分的强度影响
冠层量子效益ac: ac=acxf(P)f(T)f(F)f(N) ac=acxf(P)f(T)f(F)f(N)f(C) 考虑了林木的竞争因子
根生物量分配比
例ηR:
ηR =
ηRxηRn
ηRn+(ηRx-ηRn)mφ
ηR =
ηRxηRn
ηRn+(ηRx-ηRn)mφ(1-n)
考虑了光照对根生物量分配的
影响
2.1 生物量积累子模型改进
生物量的积累利用林木生态机理模型,同时考
虑温度、降水、光照、土壤等多个环境因子来建立生
物量生产过程,允许不同时期有不同的模型限制因
子,对讨论林木自身和环境因素对林木生长的生理
生态过程影响进行定量研究。本研究在 3PG模
型[9]中经验生产子模型的基础上,充分考虑单木的
冠幅、开阔比数及林木竞争等具体影响因素,对光合
有效辐射模型和冠层量子效益模型进行了改进。
(1)光合有效辐射APAR:
APAR=RAD(1-e-kLi)e
cw
cwi
(f(o)-1) (1)
式中:RAD是太阳总辐射,Li是单木叶面积指
数,k是消光系数,f(o)是单木的开阔比数,cwi是第
i株单木的冠幅,cw是结构单元中的平均冠幅。
该子模型中底数方程 e
cw
cwi
(f(o)-1)是单株木获取
光照养分的强度分析,使得光合有效辐射计算更加
合理。单木的有效冠幅越大,获取有效光照就越
充分。
(2)冠层量子效益ac:
ac=acxf(P)f(T)f(F)f(N)f(C)
(2)
式中:acx是光能利用率的取值上限,f(P)是生
理修正因子,f(T)是温度修正因子,f(F)是冰冻修
正因子,f(N)是养分修正因子,f(C)是林木竞争强
255
第4期 马莉燕等:基于单木过程模型的杉木林分生长模拟研究
度修正因子。
该子模型在原有基础上增加林木竞争强度修正
因子,f(C)为单木竞争能力与邻体结构单元(包括
对象单木)中最大竞争能力的比值。实际林地中单
木的竞争能力有差异,竞争能力越大,获取土壤养分
就多,光合生产能力就大,这种差异在土壤贫瘠的林
分中表现更为明显[10]。
(3)净初级生产量NPP:
NPP=0.552 acAPARCpp (3)
式中:Cpp是净、总生产量之比。
2.2 生物量分配子模型改进
生物量分配主要指光合作用同化的碳,以碳水
化合物的形式在林木不同器官组织间分配,从而调
控林木的生长。3PG模型中碳量按一定比例分配
给树叶、树干和树根[11],高光照水平下,林木将加大
对根的生物量分配,低光照水平下增加对叶的生物
量分配[12]。因此,本研究具体考虑光照、土壤以及
林木竞争对生物量分配的影响,各器官生物量分配
模型公式如下:
(1)根生物量分配比例ηR:
ηR =
ηRxηRn
ηRn+(ηRx-ηRn)mφ(1-n)
(4)
其中,m=[m0+(1-m0)FR]q;n=n0[1
-e-k+e-kf(o)]
式中,ηR为根生物量分配比例,ηRx为根系生物
量分配比例取值上限,ηRn为根系生物量分配比例取
值下限,m、n分别为土壤肥力、光照对分配比例的影
响系数,φ为生理影响综合因子;FR为林分土壤肥
力等级,q为林分平均竞争与林木竞争之比,m0为
FR取零时m的取值,n0为每月实际 PAR与年平均
PAR的比值,f(o)为单木的开阔比数。
(2)叶生物量分配比例ηS:
ηS =
1-ηR
1+PFS
(5)
(3)干生物量分配比例ηF:
ηF =PFSηS (6)
式中,PFS为干、叶间生物量的比例。
3 模拟结果与分析
将2008年样地林木相关数据参数,应用改进模
型,得到2008—2013年单木生长量以及分配到各器
官的量,2008年样地每月林分各器官的增长生物量
如表2所示。由于实际测量中,单木各器官的生物
量无法直接测得,而胸径和树高是表征单木生长的
直接指标,可以反映出单木各器官的生物量分配与
增长的情况,并且容易实际测量,因此选择胸径和树
高两个指标进行分析。林木各组织所分配的生物量
进行形态结构生长的可视化过程通常采用生物量与
胸径、树高等形态参数之间的异速生长方程,本研究
中主要借鉴其他研究者得出的杉木异速生长方程,
如干生物量回归方程、叶生物量回归方程等[9,13]。
最后对比分析单木预测生长量与实际生长量的吻合
程度,评价预测模拟效果,2008—2013年样地单木
胸径及树高生长量模拟情况如表3和表4所示。
表2 2008年每月林分各器官的增长生物量
月份 干生物量/kg 叶生物量/kg 根生物量/kg
1 0.1987 0.0234 0.0418
2 0.1701 0.0183 0.0355
3 0.4189 0.0494 0.0886
4 0.5186 0.0560 0.1090
5 0.6523 0.0782 0.1388
6 0.7592 0.0986 0.1631
7 1.1296 0.1694 0.2474
8 0.7492 0.1048 0.1623
9 0.5813 0.0755 0.1246
10 0.2718 0.0326 0.0575
11 0.2039 0.0220 0.0395
12 0.2281 0.0246 0.0441
年总计量 5.8821 1.2527 0.7534
由表2可知,林分生长过程中生物量的增长在
单木各器官的分配比例不同,干 >根 >叶。各器官
的生物量增长量先逐月增加后逐月减少,最大相差
约7倍,7月份各器官的生物量增长最快。说明气
候条件对林分生物量增长影响显著,充足的光照和
降水以及较高的气温条件有利于林分快速生长。
由表 3可知,改进模型的模拟平均精度达到
9557%,3PG模型的模拟平均精度为89.4%。因
此,改进模型的模拟精度明显高于3PG模型,模拟
效果更好。
由表 4可知,改进模型的模拟平均精度达到
9116%,3PG模型的模拟平均精度为85.18%。因
此,改进模型的模拟精度明显高于3PG模型,模拟
效果更好。
将改进模型与3PG模型得出的单木胸径和树
高的生长量模拟值分别与对应的实测值进行回归分
析,结果如图1和图2所示。改进模型对单木胸径
生长量的模拟能够满足93.81%的相应实际生长量
355
林 业 科 学 研 究 第28卷
表3 2008—2013年单木胸径生长量模拟
实测值
2008年/cm 2013年/cm
改进模型
2013年模拟/cm 精度
3PG模型
2013年模拟/cm 精度
18.9 21.2 21.8 0.9717 22.3 0.9481
12.4 18.2 17.1 0.9396 15.8 0.8681
9.3 11.0 10.9 0.9909 12.7 0.8455
14.8 17.4 18.7 0.9253 18.2 0.9540
11.8 13.4 13.8 0.9701 15.2 0.8657
14.3 18.9 18.1 0.9577 17.7 0.9365
15.9 20.4 19.5 0.9559 19.3 0.9461
15.2 22.5 21.2 0.9422 18.6 0.8267
15.0 17.8 18.9 0.9382 18.4 0.9663
11.9 13.5 13.8 0.9778 15.3 0.8667
18.6 20.3 20.1 0.9901 22.0 0.9163
10.1 12.3 13.1 0.9350 13.5 0.9024
14.1 20.8 19.3 0.9279 17.5 0.8413
15.5 21.1 21.5 0.9810 18.9 0.8957
13.6 19.1 18.1 0.9476 17.0 0.8901
15.7 21.2 20.3 0.9575 19.1 0.9009
12.0 17.7 16.9 0.9548 15.4 0.8701
11.9 14.2 14.7 0.9648 15.3 0.9225
12.3 18.8 17.2 0.9149 15.7 0.8351
14.6 20.4 19.8 0.9706 18.0 0.8824
…
平均 0.9557 0.8940
表4 2008—2013年单木树高生长量模拟
实测值
2008年/m 2013年/m
改进模型
2013年模拟/m 精度
3PG模型
2013年模拟/m 精度
8.3 14.0 12.7 0.9071 12.6 0.9000
5.2 15.7 13.7 0.8726 11.5 0.7325
7.1 14.5 13.6 0.9379 11.4 0.7862
8.1 14.4 12.5 0.8681 12.4 0.8611
8.7 15.2 13.2 0.8684 13.0 0.8553
7.2 13.4 12.4 0.9254 11.5 0.8582
8.1 12.9 12.5 0.9690 11.4 0.8837
7.7 13.2 13.7 0.9621 12.0 0.9091
7.3 15.0 12.6 0.8400 11.6 0.7733
7.1 14.1 12.7 0.9007 11.4 0.8085
8.4 14.2 12.7 0.8944 12.7 0.8944
6.7 14.7 12.6 0.8571 11.0 0.7483
7.4 14.9 12.9 0.8658 11.7 0.7852
8.4 14.7 13.4 0.9116 12.7 0.8639
7.8 14.8 13.2 0.8919 12.6 0.8514
8.8 12.0 12.1 0.9917 13.1 0.9083
…
平均 0.9116 0.8518
变化,较之3PG模型对应的90.85%有显著提高;
改进模型对单木树高生长量的模拟能够满足
8447%的相应实际生长量变化,较之3PG模型对
应的81.34%也有显著提高。因此,基于单木生长
过程的改进模型基本能够适用于杉木林分的生长
模拟。
改进的生理生态模型对杉木林分生长的预测效
果可通过二维和三维可视化图像展示出来,生长后
的林分模拟图像如图3所示。二维图反映了林分单
木的胸径生长量及其在林地的分布情况,可以看出,
455
第4期 马莉燕等:基于单木过程模型的杉木林分生长模拟研究
图1 改进模型的模拟值与实测值回归分析
图2 3PG模型的模拟值与实测值回归分析
图3 生长后的林分模拟图
由于单木在林地的分布疏密不同,造成光照、水分、
土壤肥力等分配不均,单木胸径生长量也不同,相对
密集处的林分单木冠幅差异更大。三维图可以清楚
直观地反映林分的空间结构及其生长状况,说明改
进模型对杉木林分的生长模拟具有参考意义。
4 结论
运用改进模型对杉木生长过程进行模拟,选择
胸径和树高生长量2个指标与实测值对比分析,结
果表明:模型模拟值与观测值呈显著相关关系(P<
0.001),对应的模拟值能有效地解释胸径90%以上
的变化,对树高也能达到85%以上,说明基于单木
生长过程的改进模型适用于杉木林分生长预测。改
进模型对胸径、树高的平均模拟精度均明显高于3
PG模型,说明改进模型的模拟精度显著提高,模拟
效果更好。
模拟二维图反映了林分单木的胸径生长量及其
在林地的分布情况,三维图清楚直观地反映出林分
的空间结构及其生长状况,进一步验证了改进模型,
说明基于单木过程模型的杉木林分生长模拟对林分
生长预测以及可视化具有参考作用。
本研究中的改进模型,综合考虑了林木自身以
及空间结构带来的生长影响,对杉木林分中单木的
生长预测具有较高的精度,并对天然混交林中林木
的生长预测同样适用。但是,模型中参数仍然较多,
获取较复杂,常常参考其它研究中的参数,因此该模
型在林分生长模拟中的应用还有一定限制,进一步
优化模型具有重要意义。
参考文献:
[1]刘彦君,游先祥,叶 影.三维可视化技术在林业中的应用研究
进展[J].林业调查规划,2010,35(6):9-13.
555
林 业 科 学 研 究 第28卷
[2]LandsbergJJ,WaringRH.Ageneralisedmodelofforestproduetivi
tyusingsimPlifiedconceptsofradiationuseeficiency,carbonbal
anceandpartitioning[J].ForestEcologyandManagement,1997,
95:209-228.
[3]郭明春,刘建峰.杉木人工林林分生长模拟[J].福建林学院学
报,2011,31(1):65-68.
[4]AlmeidaAC,LandsbergJJ,SandsPJ.Parameterisationof3PG
forFastGrowingEucalyptusgrandisPlantations[J].ForestEcology
andManagement,2004,193(1-2):179-195.
[5]CoopsNC,WaringRH,LandsbergJJ.Estimationofpotentialfor
estproductivityacrosstheOregontransectusingsatelitedataand
monthlyweatherrecords[J].InternationalJournalofRemoteSens
ing,2001,22,3797-3812.
[6]OlgaP,QuinnH,JustinM.DesignofaGISBasedWebApplica
tionforSimulatingBiofuelFeedstockYields[J].ISPRSIntJGeo
Inf,2014,3(3),929-941.
[7]花利忠,江希钿.3PG模型在华南尾叶桉人工林的应用研究
[J].北京林业大学学报,2007,29(2):100-104.
[8]梁诗博.杉木林分生长与环境交互过程可视化模拟研究[D].
中国林业科学研究院,2012.
[9]赵梅芳.基于3PG机理模型的杉木林碳固定及蒸散量模拟研
究[D].中南林业科技大学,2008.
[10]方宗辉.福建省杉木人工林生物量及其分配研究[J].福建林
业科技,2005,32(3):82-85.
[11]朱智强,蒋菊生.3PG生长模型及其在橡胶树栽培领域的应用
[J].热带农业科学,2010,30(7):4-7.
[12]侯振宏,张小全,徐德应.杉木人工林生物量和生产力研究
[J].中国农学通报,2009,25(5):97-103.
[13]张 林,黄 永,罗天祥,等.林分各器官生物量随林龄的变化
规律———以杉木、马尾松人工林为例[J].中国科学院研究生
院学报,2005,22(2):171-178.
655