基于beta回归的迎春5号杨树树干密度混合效应模型.pdf
《基于beta回归的迎春5号杨树树干密度混合效应模型.pdf》由会员分享,可在线阅读,更多相关《基于beta回归的迎春5号杨树树干密度混合效应模型.pdf(12页珍藏版)》请在咨信网上搜索。
1、DOI:10.12171/j.10001522.20220450基于 beta 回归的迎春 5 号杨树树干密度混合效应模型吴新华苗铮郝元朔董利虎(东北林业大学林学院,黑龙江哈尔滨150040)摘要:【目的】探究迎春 5 号杨树在树干纵向上的木材密度影响因子和变异规律,构建迎春 5 号杨树边材、心材、树皮和树干密度混合效应 beta 回归模型,为树干生物量预测和木材材性研究提供参考。【方法】以黑龙江省尚志市 90 株迎春 5 号杨树解析木数据为基础,构建迎春 5 号杨树边材、心材、树皮和树干密度的混合效应 beta 回归模型。采用相关性分析和最优子集法筛选 beta 回归基础模型的变量;利用负二
2、倍的对数似然值、赤池信息准则、贝叶斯信息准则、调整确定系数(Ra2)、似然比检验对收敛模型进行拟合优度的评价,利用留一交叉验证法对模型进行检验,指标为平均绝对误差(MAE)和平均绝对百分比误差;结合两种抽样方式(方案:不限定相对高;方案:限定相对高在 0.1 以下)对模型进行校正。【结果】边材、心材、树皮和树干密度不仅受到相对高的影响,还分别与胸径平均生长量、年龄、胸径密切相关,基于林木因子建立的混合效应 beta 回归模型的 Ra2分别为 0.53、0.52、0.52、0.63,MAE0.05g/cm3,与基础模型相比均提高了预测精度。边材和心材密度从树干基部往上先减小后增大,在相对高 0.
3、2 处有拐点;树皮密度从树干基部到树梢先增大后减小,在相对高 0.6 处有拐点;树干密度沿着树干向上逐渐增大。固定相对高时,边材、心材密度都与胸径平均生长量呈负相关,树皮、树干密度分别与年龄、胸径呈负相关。在不限定相对高的情况下,沿着树干随机抽取 4 个圆盘的密度测量值来校准模型得到稳定的预测精度;限定取样高度在相对高 0.1(2.0m)以下时,对边材、心材、树皮和树干分别抽取一个圆盘(对应高度为 1.0、1.3、2.0、1.0m)的密度测量值,得到与最优抽样组合相似的预测精度。相对高、胸径平均生长量、年龄和胸径是迎春 5 号杨树木材密度的显著影响因子。【结论】beta 回归模型可对(0,1)
4、区间的迎春 5 号杨树树干密度直接模拟,引入随机效应可提高模型的预测精度。边材、心材、树皮和树干密度在树干纵向上的变化规律不同,构建的混合效应 beta 回归模型可为迎春 5 号杨树树干生物量估算和木材性质研究奠定基础。关键词:迎春 5 号杨树;木材密度;beta 回归;广义线性混合模型中图分类号:S781.31;S792.11;S758.1文献标志码:A文章编号:10001522(2023)05006712引文格式:吴新华,苗铮,郝元朔,等.基于 beta 回归的迎春 5 号杨树树干密度混合效应模型 J.北京林业大学学报,2023,45(5):6778.WuXinhua,MiaoZheng,
5、HaoYuanshuo,etal.MixedeffectmodelofstemdensityofPopulus nigra P.simoniibasedonbetaregressionJ.JournalofBeijingForestryUniversity,2023,45(5):6778.Mixed effect model of stem density of Populus nigra P.simoniibased on beta regressionWuXinhuaMiaoZhengHaoYuanshuoDongLihu(SchoolofForestry,NortheastForestr
6、yUniversity,Harbin150040,Heilongjiang,China)Abstract:Objective This paper aims to explore the influencing factors and variation rules of wooddensityinthelongitudinalstemofPopulus nigraP.simonii,sobetaregressionmodelswithmixedeffectofsapwood,heartwood,barkandstemdensityofthepoplarwereconstructed,whic
7、hwasusedasareferenceforstembiomasspredictionandwoodtimberproperties.MethodMixedeffectbetaregressionmodelsforsapwood,heartwood,barkandstemdensityofP.nigra P.simonii wereestablished,whichbasedontheanalyticaldataof90treesofP.nigrax P.simoniiplantationinShangzhiCity,HeilongjiangProvinceof收稿日期:20221107修回
8、日期:20230315基金项目:中央高校基本科研业务费专项(2572020DR03),黑龙江头雁创新团队计划项目(森林资源高效培育技术研发团队)。第一作者:吴新华。主要研究方向:木材密度研究。Email:地址:150040 黑龙江省哈尔滨市香坊区和兴路 26 号东北林业大学林学院。责任作者:董利虎,博士,教授。主要研究方向:林分生长与收获模型、生物量、碳储量。Email:地址:同上。本刊网址:http:/;http:/第45卷第5期北京林业大学学报Vol.45,No.52023年5月JOURNALOFBEIJINGFORESTRYUNIVERSITYMay,2023northeasternCh
9、ina.Usingcorrelationanalysisandoptimalsubsetmethodstoscreenthevariablesofthebetaregressionbasemodel,andthegoodnessoffitoftheconvergencemodelwasevaluatedby2log-likehoodvalue,akaikeinformationcriterion,bayesianinformationcriterion,adjustedcertaintycoefficient(Ra2)andlikelihoodratiotest.Theleave-one-ou
10、t-cross-validationwasusedtotestthemodel,theindexesweremeanabsoluteerror(MAE)andmeanabsoluteerrorpercentage.Twosamplingmethodswerecombined(scheme:norelativeheight;scheme:limitrelativeheightbelow0.1)tocorrectthemodel.ResultThedensitiesofsapwood,heartwood,barkandstemwerenotonlyaffectedbyrelativeheight,
11、butalsocloselyrelatedtotheaveragegrowthofDBH,ageandDBH,respectively.Ra2ofthemixed-effectbetaregressionmodelbasedontreefactorswas0.53,0.52,0.52,0.63,respectively,andtheMAE0.05g/cm3.Sapwooddensityandheartwooddensitydecreasedfirstandthenincreasedfromthebasetothetopofthestem,withaninflectionpointatarela
12、tiveheightof0.2.Barkdensityfirstincreasedandthendecreasedfromthebaseofthestemtothetopofthetree,andtherewasaninflectionpointattherelativeheightof0.6.Thestemdensityincreasedgraduallyalongthestem.Whenfixedrelativeheight,thedensitiesofsapwoodandheartwoodwerebothnegativelycorrelatedwiththeaveragegrowthof
13、DBH.ThedensitiesofbarkandstemwerenegativelycorrelatedwithageandDBH,respectively.Withoutlimitingtherelativeheight,thewooddensityvaluecorrespondingtotheheightof4discsrandomlysampledalongthestemwascalibratedtoobtainstablepredictionaccuracy.Whenthesamplingheightwaslimitedto0.1(2.0m)orless,therewaslittle
14、differenceinthepredictionaccuracybetweentheoptimalsamplingcombinationandthedensityvalues(1.0,1.3,2.0,1.0m,respectively)ofsapwood,heartwood,barkandstematadischeight.Relativeheight,averagegrowthof DBH,age and DBH were significant influencing factors of wood density of P.nigra P.simonii.ConclusionThebe
15、taregressionmodelcandirectlysimulatethestemdensityofP.nigra P.simoniiinthe(0,1)interval,andtherandomeffectcanimprovethepredictionaccuracyofthemodel.Thelongitudinalvariationsofsapwood,heartwood,barkandstemdensityaredifferent.Theconstructedmixed-effectbetaregressionmodelcanlayafoundationforbiomassesti
16、mationandwoodpropertystudyofP.nigra P.simonii.Key words:Populus nigra P.simonii;wooddensity;betaregression;generalizedlinearmixedmodel木材密度又称为木材基本密度,通常用绝干材质量除以生材体积来表示。木材密度与碳储量、木材机械强度、水分运输、耐旱性等密切相关12,在许多研究利用树干木材密度乘以树干体积来计算树干生物量34,进而估算森林碳储量5。碳储量与生态环境变化息息相关,所以木材密度对研究森林生态系统固碳功能具有重要作用6。木材密度也是林木生长率和死亡率的重要预
17、测因子,可以反映林木的生长和性状表现7,通常低木材密度代表高生长量810。准确估算木材密度不仅有助于阐明地上生物量的空间分配格局,还能反映树干内部的机械稳定性和生长速率的变化1112。作为森林碳汇计量中的重要生物量估算因子13,木材密度主要受树种、相对高、树龄、生长量、胸径和环境条件等因素的影响1417。有研究表明杨树(Populusspp.)的木材密度沿树干自下往上逐渐增大18,而落叶松(Larix gmelinii)却呈现出相反的规律。樟子松(Pinus sylvestris)、日本花柏(Chamaecyparispisifera)等树种木材密度沿树干自下往上先减小后增大1920,这表明不
18、同树种间的木材密度变化规律具有较大差异。Iida 等21的研究表明树高一定时,较小的胸径具有较高的木材密度。木材密度还与年轮宽度,尤其与晚材比例密切相关,在环孔材中更加明显22。研究表明干旱条件会抑制树木生长进而增加晚材比例,获得较大的木材密度,进一步说明低含水率导致较高的木材密度23。目前对松属(Pinus)、杉科(Taxodiaceae)以及硬阔树种木材密度的研究较多19,2425,暂未见迎春 5 号杨树(P.nigra P.simonii)木材密度的相关报道。各研究显示不同树种之间木材密度有所差异,同一树种内部木材密度也不尽相同,例如边材、心材和树皮木材密度;径向和纵向木材密度;幼树和成
19、熟树木材密度26。准确把握树干各部分的木材密度变化规律,合理预测木材密度,有利于生物量的精确预估。早期的研究者常利用线性和非线性回归来构建木材密度模型,种类有多项式、双曲线、指数方程等27。但木材密度数据位于(0,1)区间内,在数据结构上属于连续有界数据,利用上述模型对该类数据68北京林业大学学报第45卷进行拟合时,容易产生超过其上下限的值,降低预测精度28。此外,(0,1)区间的数据存在分布不对称的问题,无法基于正态性假设来进行推断预测。广义beta 回归模型非常灵活,能够直接应用于连续有界数据,并且可以适应“J”形、“L”形、“单峰”、“U”形和矩形的各种类型的分布数据29,解决了数据不对
20、称分布问题。另外,木材密度在不同树木、同一树木的不同位置间具有较大的变异性,且数据一般来源于固定样地或样木,数据之间存在一定的相关性,并不满足传统回归模型中独立同分布的假设条件,故利用混合效应模型对木材密度进行建模,通过随机效应方差协方差结构来表现个体间的差异和个体内的自相关性,提高模型预测精度3032。迎春 5 号杨树是小叶杨(P.simonii)和黑杨(P.nigra)杂交而成的优良无性系。该品种自 1982 年从林科院引入牡丹江林区,是东北地区主要造林树种之一,同时也是培育速生丰产林及短轮伐期工业原料林的优良品种。基于经济效益,前人对于迎春5 号杨树的研究主要集中在造林方法、种植密度和经
21、营技术方面,暂未见其木材密度以及生物量的研究。本研究从生态效益出发,构建迎春 5 号杨树木材密度模型,进而估算其树干生物量,为杨树人工林碳储量和木材材性研究奠定基础,具有生态学意义。1研究区概况与研究方法 1.1 研究区概况研究区位于黑龙江省尚志市的尚志林场(450008452328N,12739321293509E)、元宝林场(450826451826N,12810171283307E)、老街基林场(440906450401N,12744091280506E)和苇河林场(440109445201N,12810521284704E)。各林场隶属于尚志国有林场管理局,地处长白山脉张广才岭西坡,地
22、形以山地丘陵为主,地势西低东高,土壤以暗棕壤为主。环境气候条件相似,均为中温带大陆性季风气候,冬季寒冷干燥,夏季温暖湿润,降水量充沛。森林资源丰富,尚志国有林场森林总面积为152737hm2,其中人工林面积 49168hm2,占比32.19%。植被属长白山植物区系,人工林主要树种为长白落叶松(Larix olgensis)、樟子松、红松(Pinuskoraiensis)和杨树。2020 年 9 月在尚志市上述林场选择不同年龄、不同林分密度、不同海拔的迎春 5 号杨树人工林设置标准地共 18 块,其中尚志林场 11 块,元宝林场3 块,老街基林场 1 块,苇河林场 3 块。标准地的形状为矩形。由
23、于林分疏密程度不同,所设标准地面积 0.060.10hm2不等。各标准地所在林场的相关立地因子如表 1 所示。在每块标准地中进行每木检尺,按照等断面积标准木法将每块标准地内的林木划分为 5 级,径阶从大到小依次为优势木、一级木、二级木、三级木和劣势木。以每级木的平均胸径为标准,在每块标准地附近选取 5 棵胸径大小误差为5%左右的林木作为解析木,对其进行树干解析,共获得 90 棵解析样木。表1各林场相关的立地因子Tab.1Sitefactorsrelatedtoeachforestfarm林场Forestfarm地形Terrain海拔Altitude/m坡度Slope坡向Aspect坡位Slop
24、eposition尚志林场ShangzhiForestFarm山坡Hillside20131905、515阳、阴Sunny,shady中下Lower-middle元宝林场YuanbaoForestFarm平地Flat20624505平Flat老街基林场LaojiejiForestFarm平地Flat30005平Flat苇河林场WeiheForestFarm山坡Hillside26032405、510阳、阴Sunny,shady中下Lower-middle在进行树干解析前,要在树干上标记北向和1.3m 处,并测量胸径。解析木伐倒后,测量其树高和年龄等因子。从解析样木根颈(0m)处每隔 1.0m截
25、取 2cm 厚的圆盘样品,截取过程中避免节子,胸径(1.3m)处单独取盘,直至距离树梢不足 1.0m 为止。按高度对圆盘进行编号,在非工作面标注编号和北向,将不同解析样木的圆盘分别装袋带回实验室。各圆盘以髓心为定点,取适当大小的楔形或全部圆盘作为该圆盘的样品。由于迎春 5 号杨树的边材、心材和树皮颜色差异明显,因此将样品分为边材、心材和树皮 3 部分。若样品没有心材,则将其分为边材和树皮两部分。利用排水法分别测量边材、心材和树皮样品入水前后的质量,利用水的密度分别计算排出水的体积,即为边材、心材、树皮样品的生材体积。最后将样品置于 85 烘箱 24h,测量样品的绝干质量,木材密度表示为每单位生
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 beta 回归 迎春 杨树 树干 密度 混合 效应 模型
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。