林龄和立地条件对杉木林乔木层生物量分配的影响

第４９卷第１期２０２２年３月福建林业科技
ＪｏｕｒｏｆＦｕｊｉａｎＦｏｒｅｓｔｒｙＳｃｉａｎｄＴｅｃｈ
Ｖｏｌ ４９　Ｎｏ １
Ｍａｒ ，２０２２
ｄｏｉ：１０．１３４２８／ｊ．ｃｎｋｉ．ｆｊｌｋ．２０２２．０１．００５
林龄和立地条件对杉木林乔木层生物量分配的影响
翁建宇
（福建省将乐国有林场，福建将乐３５３３００）
摘要：于２０１０—２０１９年，采用样地调查法对福建省三明市将乐县国有林场杉木人工林的地况因子、林况因子及乔木层的生物量进行调查；通过构建引入林龄和地位指数的相容性生物量模型，估算杉木林乔木层各组分生物量，研究各组分生物量分配规律与立地条件和林龄之间的关系。

结果表明：引入林龄和地位指数的相容性生物量模型，能够更好地反映林龄和立地条件对生物量异速生长规律和分配的影响；杉木各组分生物量的分配比例受林分林龄的影响显著，其中干材、树皮生物量占比随林龄的增大而增大，树叶、树根则相反，树枝生物量占比先增加后减小；除树枝外，立地条件不会改变各组分生物量占比随林龄的变化趋势，但会影响各组分生物量占比大小。

在立地条件较好的林分中，树干生物量占比较大，树叶、树根的分配则相反。

研究结果可为杉木人工林碳汇经营提供参考。

关键词：杉木人工林；生物量分配；立地条件；林龄
中图分类号：Ｓ７９１ ２７；Ｓ７１８ ５５＋６ 文献标识码：Ａ 文章编号：１００２－７３５１（２０２２）０１－００３０－１０
ＥｆｆｅｃｔｓｏｆＳｔａｎｄＡｇｅａｎｄＳｉｔｅＣｏｎｄｉｔｉｏｎｓｏｎＳｔａｎｄＢｉｏｍａｓｓＡｌｌｏｃａｔｉｏｎｓｏｆ
ＣｕｎｎｉｎｇｈａｍｉａｌａｎｃｅｏｌａｔａＰｌａｎｔａｔｉｏｎｓ
ＷＥＮＧＪｉａｎｙｕ
（Ｊｉａｎｇｌｅｓｔａｔｅ ｏｗｎｅｄｆｏｒｅｓｔｆａｒｍ，Ｊｉａｎｇｌｅ３５３３００，Ｆｕｊｉａｎ，Ｃｈｉｎａ）
Ａｂｓｔｒａｃｔ：Ｓｉｔｅｃｏｎｄｉｔｉｏｎｓ，ｆｏｒｅｓｔａｔｔｒｉｂｕｔｅｓａｎｄｂｉｏｍａｓｓｏｆｔｒｅｅｌａｙｅｒｗｅｒｅｍｅａｓｕｒｅｄｆｏｒ７７ｓａｍｐｌｅｐｌｏｔｓｏｆＣｕｎｎｉｎｇｈａｍｉａｌａｎｃｅｏｌａｔａｐｌａｎｔａｔｉｏｎｓｉｎＪｉａｎｇｌｅＳｔａｔｅＯｗｎｅｄＦｏｒｅｓｔＦａｒｍ，Ｓａｎｍｉｎｇ，ＦｕｊｉａｎＰｒｏｖｉｎｃｅ，ｉｎ２０１０—２０１９ ＴｈｅｂｉｏｍａｓｓｏｆｅａｃｈｃｏｍｐｏｎｅｎｔｉｎｓｔａｎｄｏｆＣｕｎｎｉｎｇｈａｍｉａｌａｎｃｅｏｌａｔａｐｌａｎｔａｔｉｏｎｓｗａｓｅｓｔｉｍａｔｅｄｂｙａｃｏｍｐａｔｉｂｌｅｂｉｏｍａｓｓｍｏｄｅｌｗｉｔｈａｇｅａｎｄｓｉｔｅｉｎｄｅｘ Ｔｈｅｒｅｌａｔｉｏｎｓｈｉｐｏｆｂｉ ｏｍａｓｓａｌｌｏｃａｔｉｏｎｓａｎｄｓｉｔｅｃｏｎｄｉｔｉｏｎｓａｎｄｓｔａｎｄａｇｅｗａｓａｎａｌｙｚｅｄ Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔ：Ｔｈｅｃｏｍｐａｔｉｂｌｅｂｉｏｍａｓｓｍｏｄｅｌｗｉｔｈａｇｅａｎｄｓｉｔｅｉｎｄｅｘｃｏｕｌｄｃｏｍｍｅｎｄａｂｌｙｄｅｓｃｒｉｂｅｔｈｅｅｆｆｅｃｔｓｏｆｓｔａｎｄａｇｅａｎｄｓｉｔｅｃｏｎｄｉｔｉｏｎｓｏｎｔｈｅａｌｌｏｍｅｔｒｙａｎｄａｌｌｏｃａｔｉｏｎｓｏｆｆｏｒｅｓｔｓｓｔａｎｄｂｉｏｍａｓｓ；Ｔｈｅｐｒｏｐｏｒｔｉｏｎｏｆｂｉｏｍａｓｓｏｆｅａｃｈｃｏｍｐｏｎｅｎｔｗａｓｓｉｇｎｉｆｉｃａｎｔｌｙａｆｆｅｃｔｅｄｂｙｓｔａｎｄａｇｅ Ｔｈｅｐｒｏｐｏｒｔｉｏｎｏｆｂｏｌｅａｎｄｂａｒｋｂｉｏｍａｓｓｉｎｃｒｅａｓｅｄｗｉｔｈｓｔａｎｄａｇｅ Ｔｈｅｐｒｏｐｏｒｔｉｏｎｏｆｌｅａｖｅｓａｎｄｒｏｏｔｓｂｉｏｍａｓｓｄｅｃｒｅａｓｅｄｗｉｔｈｓｔａｎｄａｇｅ Ｔｈｅｐｒｏｐｏｒｔｉｏｎｏｆｂｒａｎｃｈｂｉｏ ｍａｓｓｉｎｃｒｅａｓｅｄｆｉｒｓｔｌｙａｎｄｔｈｅｎｄｅｃｒｅａｓｅｄ；Ｅｘｃｅｐｔｆｏｒｂｒａｎｃｈｅｓ，ｓｉｔｅｃｏｎｄｉｔｉｏｎｓｄｉｄｎ’ｔｃｈａｎｇｅｔｈｅｔｒｅｎｄｏｆｐｒｏｐｏｒｔｉｏｎｏｆｂｉｏｍａｓｓｏｆｅａｃｈｃｏｍｐｏｎｅｎｔ Ｈｏｗｅｖｅｒ，ａｌａｒｇｅｒｐｒｏｐｏｒｔｉｏｎｏｆｓｔｅｍｂｉｏｍａｓｓａｎｄｓｍａｌｌｅｒｐｒｏｐｏｒｔｉｏｎｏｆｌｅａｖｅｓａｎｄｒｏｏｔｓｂｉｏｍａｓｓｗａｓｏｃｃｕｐｉｅｄｂｙｔｈｅｓｔａｎｄａｔｓｕｐｅｒｉｏｒｓｉｔｅｃｏｎｄｉｔｉｏｎｓ ＴｈｅｒｅｓｕｌｔｓｏｆｔｈｉｓｓｔｕｄｙｍｉｇｈｔｐｒｏｖｉｄｅｒｅｆｅｒｅｎｃｅｆｏｒｃａｒｂｏｎｓｅｑｕｅｓｔｒａｔｉｏｎｍａｎａｇｅｍｅｎｔｏｆＣｕｎｎｉｎｇ ｈａｍｉａｌａｎｃｅｏｌａｔａｐｌａｎｔａｔｉｏｎｓ
Ｋｅｙｗｏｒｄｓ：Ｃｕｎｎｉｎｇｈａｍｉａｌａｎｃｅｏｌａｔａｐｌａｎｔａｔｉｏｎｓ；Ｂｉｏｍａｓｓａｌｌｏｃａｔｉｏｎ；ｓｉｔｅｃｏｎｄｉｔｉｏｎ；ｓｔａｎｄａｇｅ
森林生态系统是陆地生态系统组分中，面积最大、最重要的生态系统之一［１］。

同时，也是陆地生态系统中最大的碳库，具有维持生态系统碳平衡，减缓全球变暖、减少大气温室气体含量的功能［２］。

森林生物
量是计算碳储量和评估固碳能力的重要依据。

同时，森林生物量大小也直接关系到森林生产力的高低。

因此，准确测算森林生物量，对合理开展森林碳汇经营具有十分重要的意义。

森林生物量分配规律研究是森林生物量研究的重要组成部分，也是研究森林生态系统能量流动和物质转移的基础［３］。

由于森林生物量受立地条件的地形、水分、光照等自然因素以及胸径、树高等测树因子的影响，研究不同立地条件下森林生物量在各组分之间的分配规律十分重要［４］。

　收稿日期：２０２１－０９－２４；修回日期：２０２１－１１－１０
　基金项目：北京林业大学科技创新计划（立地条件对杉木人工林生物量的影响研究，ＢＬＸ２０１６１４）
　作者简介：翁建宇（１９８６—），男，福建省将乐国有林场林业工程师，从事森林经营管理、林业新业态产业工作。

Ｅ ｍａｉｌ：１０５７４１８４７＠ｑｑ ｃｏｍ。

第１期翁建宇：林龄和立地条件对杉木林乔木层生物量分配的影响
国内外学者已开展了很多森林生物量分配方面的研究。

Ｅｂｅｒｍｅｒｙｅｒ等［５
］早在１９世纪７０年代就开始
对几种森林的树枝落叶量和木材重量进行测定和研究，其后大量生物量模型相继建立，研究方法和准确性得到了发展和提升［６－９
］。

２０世纪８０年代起，学者们对森林生物量和林分调查因子两者之间的关系，以及森林生物量的动态变化开展了大量研究［１０－１５
］，而且森林生物量模型形式和模型参数的求解一直是相关领域的重点问题，与此同时，遥感技术的引入使得森林生物量模型在大尺度上的应用成为研究热点［１６－１９
］。

我国许多学者对不同地区、不同气候类型、不同树种的森林生物量和森林生产力进行大量研究，并取得了
丰硕成果［２０－３３
］。

２０世纪９０年代，
国内对森林生物量的主要研究方向是大尺度森林生物量模型的构建［２１，２８
］。

２１世纪以来，
我国许多学者也开始将遥感技术运用在森林生物量和生产力的研究上。

关于森林立地条件的研究，目前大多采用综合因子法，辅以样地调查和遥感技术，且大多侧重于对经济树种和极端立地条件的研究［３４－３８
］。

杉木人工林是我国主要用材林和生态防护林，准确估算其生物量是评价其碳汇功能和提高碳汇能力的基础。

然而，关于林龄和立地条件对其生物量积累速率和组分分配的影响规律尚不明确。

因此，本文通过构建引入林龄和地位指数的相容性生物量模型，估算杉木人工林乔木层各组分生物量，研究各组分生物量分配规律与立地条件和林龄之间的关系，以期为杉木人工林经营提供参考。

１　
研究区域与研究方法
１ １　
研究区概况
研究区位于福建省三明市将乐县国有林场，地处武夷山脉东南坡，海拔２００～６００ｍ，距离最近的海岸约２４０ｋｍ。

属于亚热带季风气候区，由大陆性气候主导，兼具海洋性气候的特点，四季分明，年均气温
１８ ６℃，
年均降水量１６７６ ３ｍｍ，全年降水集中在夏季，年均蒸发量１２０４ｍｍ，年均无霜期２９８ｄ。

土壤类型多为红壤和黄红壤，土壤厚度约４３～６２ｃｍ。

主要分布的乔木树种有杉木（Ｃｕｎｎｉｎｇｈａｍｉａｌａｎｃｅｏｌａｔａ）、马尾松（Ｐｉｎｕｓｍａｓｓｏｎｉａｎａ）、木荷（Ｓｃｈｉｍａｓｕｐｅｒｂａ）
、毛竹（Ｐｈｙｌｌｏｓｔａｃｈｙｓｅｄｕｌｉｓ）等；主要灌木种有盐肤木（Ｒｈｕｓｃｈｉｎｅｎｓｉｓ）、杨桐（Ａｄｉｎａｎｄｒａｍｉｌｌｅｔｔｉｉ）等；主要草本植物有黑莎草（Ｇａｈｎｉａｔｒｉｓｔｉｓ）、海金沙（Ｌｙｇｏｄｉｕｍｊａｐｏｎｉ ｃｕｍ）
等。

１ ２　数据采集与处理
本研究于２０１０—２０１９年设置并调查杉木人工林标准地７７块，其中２０１０年１０块、２０１１年１０块、２０１２年１２块、２０１３年７块、２０１４年３块、２０１５年１块、２０１７年１３块、２０１８年１１块，２０１９年１０块。

标准地面积
为０ ０６ｈｍ２（２０ｍ×３０ｍ）
，均匀地分布于不同年龄阶段和立地条件的林分中（表１、表２）。

调查标准地的地况因子（包括样地坐标、坡度、坡位、坡向、海拔等）和林况因子（包括样木胸径、树高等），按照福建杉木
地位指数表［１６
］进行立地条件评价。

在每个标准地内选取１株平均标准木，
伐倒测算生物量，按照平均标准木法计算样地生物量。

表１　标准地各调查因子按照年龄的分布情况
年龄／ａ标准地数平均胸径／ｃｍ平均高／ｍ密度／
（株·ｈｍ－２
）生物量／（ｔ·ｈｍ－２
）总计干材树皮树枝树叶
树根
５～１０１１９ ９±２ ６６ ５±２ ２２０６０±６３４３５ １±１９ ９１３ ９±９ １３ １±１ ９４ ０±２ ６４ ５±２ ４９ ６±４ ６１１～２０２１１４ ７±３ ５１２ ３±３ １１７４４±４９９９９ ７±３９ ７５２ ７±２４ ７８ ８±３ ４１２ ２±４ ９６ ４±１ ３１９ ６±５ ８２１～３０２５１７ ４±５ ０１５ ５±３ ９１８４８±８０９１６０ ６±６５ ０９１ ９±４３ ８１３ ８±５ ３１９ ４±７ ５７ ８±１ ３２７ ７±７ ７３１～４０１８２１ ４±５ ３１７ ９±３ ８１４３６±７２９２０２ ５±７４ ６１１９ ６±５１ ０１７ ３±６ ０２４ ５±８ ５８ ５±１ ４３２ ５±８ ３４１～４９
２２１ ９±１ ８１７ ３±４ ０１３６３±１６１２１０ ５±５３ ７１２１ ２±４１ ４１８ ４±３ ７２６ １±５ ５９ ８±０ ７３５ １±３ ８总计　
７７
１６ ６±５ ７
１３ ９±５ ０
１７４１±６９９
１３７ １±７９ ５７７ ３±５１ ７
１１ ８±６ ６
１６ ６±９ ５
７ ２±２ ０
２４ ２±１０ ３
标准木的生物量采用全称重法测算。

标准木伐倒后，从树基部开始按照１ｍ为单位进行区分，梢头不足１ｍ的部分作为梢头段。

树干：分段称取鲜重，在区分段的中央部位截取厚度为２～５ｃｍ的圆盘，分别称取圆盘的干材部分和树皮部分鲜重，并按照二者的比值计算该段干材和树皮鲜重，每段分别对干材和树
·
１３·
福建林业科技第４９卷
皮取样约５０ｇ；树枝和树叶：分段称取带叶枝（树枝＋树叶）鲜重，每段分别选取１个生长良好且大小适中的标准枝，分别称取标准枝的裸枝和叶鲜重，并按照二者的比值计算该段裸枝和树叶鲜重，每段分别对树
枝和树叶取样约５０ｇ；树根：采用全根挖掘法，测定根桩、大根（≥２ｃｍ）、粗根（０ ５～２ｃｍ）、中根（
０ ２～０ ５ｃｍ）、小根（≤０ ２ｃｍ）的鲜重，并分别取样５０～１００ｇ。

将树干、树皮、树枝、树叶和根所取样品在８５℃恒温下烘干至恒重，此即为样品干重。

按照样品的干重与鲜重比值（干湿比）分别计算各段各组分干重，合计得干材、树皮、树枝、树叶、树根和全树总生物量。

表２　标准地各调查因子按照地位指数级的分布情况
地位指数级标准地数平均胸径／ｃｍ平均高／ｍ密度／
（株·ｈｍ－２
）生物量／（ｔ·ｈｍ－２
）总计干材树皮树枝树叶
树根
１２８１１ ６±３ ６９ ４±３ ３２４１８±８５５７２ ８±３８ ０３５ ４±２１ ９６ ４±３ ４８ ８±４ ９６ ０±１ ８１６ ２±６ ３１４１８１４ ９±３ ４１２ ４±２ ４１９７２±７６５１１３ ４±３９ ７５９ ７±２３ ４１０ ０±３ ５１４ ０±５ ０７ ３±１ ７２２ ３±６ ５１６１４１６ ４±６ １１３ ２±５ ８１６２７±５９４１２０ ６±７５ ８６７ ３±４７ ９１０ ５±６ ５１４ ７±９ ４６ ６±２ ０２１ ６±１０ ４１８１４１９ １±５ ６１６ ５±４ ６１５２５±４４０１８２ １±７７ １１０７ １±４９ ６１５ ５±６ ４２２ ０±９ ３７ ９±２ １２９ ７±１０ １２０１４１９ ０±６ ５１６ ５±５ １１５７１±５４４１７８ ７±９３ ６１０５ １±６２ ４１５ ２±７ ６２１ ５±１０ ７７ ８±２ １２９ １±１１ ５２２
９１７ ７±５ ９１４ ３±５ ９１４５６±７４９１３２ ８±９８ ３７５ ７±６５ ２１１ ４±８ １１５ ９±１１ ５６ ８±２ ５２２ ９±１２ １总计
７７
１６ ６±５ ７
１３ ９±５ ０
１７４１±６９９
１３７ １±７９ ５７７ ３±５１ ７１１ ８±６ ６
１６ ６±９ ５
７ ２±２ ０
２４ ２±１０ ３
１ ３　
引入林龄和立地条件的相容性生物量模型以生物量总量的独立模型为基础，采用分级控制联合估计法［２０
］建立林分尺度的相容性生物量模型。

本研究采用３级控制，首先，以生物量总量作控制，对地上部分生物量与根生物量进行联合估计。

然后，以地上部分生物量作控制，对树干生物量和树冠生物量进行联合估计。

最后，以树干生物量作控制，对干材生物量和树皮生物量进行联合估计；以树冠生物量作控制，对树枝生物量和树叶生物量进行联合估计，得到非线性联立方程组（式（１））。

这样，就满足了生物量总量与各组分生物量估计值之间的相容性，即Ｗ∧１
＝Ｗ∧
２
＋Ｗ∧
３
，Ｗ∧
２
＝Ｗ∧
４
＋Ｗ∧
５
，Ｗ∧
４
＝Ｗ∧
６
＋Ｗ∧
７
，Ｗ∧
５
＝Ｗ∧
８
＋Ｗ∧
９。

Ｗ１
＝ｆ１
（ｘ）Ｗ２
＝１
１＋ｇ１（ｘ）·Ｗ∧
１
Ｗ３
＝ｇ１
（ｘ）１＋ｇ１
（ｘ）·Ｗ∧
１
Ｗ４
＝１
１＋ｇ２（ｘ）·Ｗ∧
２
Ｗ５
＝ｇ２
（ｘ）１＋ｇ２
（ｘ）·Ｗ
∧
２
Ｗ６
＝１１＋ｇ３（ｘ）·Ｗ
∧
４
Ｗ７
＝ｇ３
（ｘ）１＋ｇ３
（ｘ）·Ｗ
∧
４
Ｗ８
＝
１
１＋ｇ４（ｘ）·Ｗ∧
５
Ｗ９
＝ｇ４
（ｘ）１＋ｇ４
（ｘ）·Ｗ∧
５
（１）
式中：Ｗ１
为生物量总量（ｔ·ｈｍ－２
）；ｆ１
（ｘ）为生物量总量的独立估计模型；Ｗ２
、Ｗ３
分别为地上部分生物量（ｔ·ｈｍ－２
）、根生物量（ｔ·ｈｍ－２
）；Ｗ∧
１
为生物量总量估计值（ｔ·ｈｍ－２
）；ｇ１
（ｘ）为根生物量与地上部分生物
量的比值函数；Ｗ４
、Ｗ５
分别为树干生物量（ｔ·ｈｍ－２
）、树冠生物量（ｔ·ｈｍ－２
）；
Ｗ∧
２
为地上部分生物量估计·２３·
第１期翁建宇：林龄和立地条件对杉木林乔木层生物量分配的影响
值（ｔ·ｈｍ－２）；ｇ２（ｘ）为树冠与树干生物量的比值函数；Ｗ６、Ｗ７、Ｗ８、Ｗ９分别为干材生物量（ｔ·ｈｍ－２）、树皮生物量（ｔ·ｈｍ－２）、树枝生物量（ｔ·ｈｍ－２）、树叶生物量（ｔ·ｈｍ－２）；Ｗ∧４、Ｗ∧５分别为树干生物量、树冠生物量估计值（ｔ·ｈｍ－２）；ｇ３（ｘ）为树冠与树干生物量的比值函数；ｇ４（ｘ）为树叶与树枝生物量的比值函数。

生物量总量的独立模型和各比值模型均以ＣＡＲ函数（式（２））为基本形式，以林分公顷断面积和林分平均高为自变量。

ｙ＝ａ·Ｇｂ·Ｈｃ（２）式中：ｙ为生物量总量或各生物量比值；Ｇ为林分公顷断面积（ｍ２·ｈｍ－２）；Ｈ为林分平均高（ｍ）；ａ、ｂ、ｃ为
参数。

为了分析地位指数（ＳＩ）对异速生长关系的影响，在参数中引入ＳＩ因子，公式（３）为以参数ａ为例的改进方程。

ｙ＝（ａ０＋ａ１·ＳＩ）·Ｇｂ·Ｈｃ（３）为了分析杉木林乔木层各组分生物量分配随年龄的动态变化，以Ｃｈａｐｍａｎ Ｒｉｃｈａｒｄｓ方程为基本形式
（式（４）），建立林分每公顷胸高断面积（Ｇ）、林分平均高（Ｈ）与林龄（Ａ）之间的关系模型，并在模型的尺度参数ａ中引入地位指数ＳＩ（式（５）），与相容性生物量方程组（式（１））构成引入林龄和立地条件的相容性生物量模型系统，从而进一步分析立地条件对杉木林乔木层各组分生物量分配动态的影响。

模型拟合和计算在Ｒ软件中实现。

ｙ＝ａ·（１－ｅ－ｂ·Ａ）ｃ（４）
ｙ＝（ａ０＋ａ１·ＳＩ）·（１－ｅ－ｂ·Ａ）ｃ（５）式中：ｙ为林分每公顷胸高断面积（ｍ２·ｈｍ－２）或林分平均高（ｍ）；Ａ为林龄（ａ）；ａ、ａ０、ａ１、ｂ、ｃ为参数。

１ ４　模型评价
采用决定系数（Ｒ２）、赤池信息准则（ＡＩＣ）、对数似然值（ｌｏｇＬｉｋ）评价模型的拟合优度，Ｒ２越大、ＡＩＣ越小、ｌｏｇＬｉｋ越大则表明模型优度越高，并通过留一交叉验证方法计算各指标值。

采用最大似然比检验法检验改进模型与基础模型之间是否有明显改进。

２　结果与分析
２ １　模型构建与评价
２ １ １　生物量模型构建及评价　关于林分生物量与胸高断面积、树高之间的异速生长关系是否受林龄的影响，不同研究中的结论不尽相同。

大多数认为，对于同一棵树而言，各组分生物量与胸径、树高之间的异
速生长是一定的，也就是说，不受林龄的影响。

故本研究不考虑林龄对此的影响。

考虑地位指数的影响，分别在基础模型的参数ａ、ｂ、ｃ中引入地位指数ＳＩ，构建各组分生物量与林分公顷断面积和平均高异速生长方程；采用对数似然比检验，若改进后的模型相对于基础模型有显著的改进，则说明ＳＩ对异速生长关系有显著影响，反之亦然。

结果表明，基础模型和改进模型均表现良好，但改进模型相对于基础模型Ｒ２变化不大，大部分引入ＳＩ的改进模型ＡＩＣ和低于基础模型，而ｌｏｇＬｉｋ高于基础模型；但从对数似然比检验看，除在ａ参数中引入ＳＩ的总生物量方程显著优于基础模型外，其它改进模型均无显著优势（表３）。

这就说明，立地条件对各组分生物量与林分公顷断面积和平均高之间的异速生长关系无显著影响，因此，在相容性生物量模型系统中仍采用基础模型作为总生物量和各组分生物量模型的模型形式。

表３　组分生物量独立模型拟合统计量
模型名称Ｒ２ＡＩＣｌｏｇＬｉｋ最大似然比检验
χ２Ｐ
Ｍｏｄ＿Ｗ
１
０ ９４４６８３ ４５４－３３７ ７２７－－
Ｍｏｄ＿Ｗ
１
ａ０ ９４７６８１ ５５１－３３５ ７７６３ ９０２０ ０４８
Ｍｏｄ＿Ｗ
１
ｂ０ ９４７６８１ ６０１－３３５ ８００３ ８５３０ ０５０
Ｍｏｄ＿Ｗ
１ｃ０ ９４７６８１ ７７１－３３５ ８８５３ ６８３０ ０５５
·３３·
福建林业科技第４９卷·４３·
表３（续）
模型名称Ｒ２ＡＩＣｌｏｇＬｉｋ最大似然比检验
χ２Ｐ
Ｍｏｄ＿Ｗ
０ ９４８６５１ ６４６－３２１ ８２３－－２
ａ０ ９４８６５３ ０５０－３２１ ５２５０ ５９７０ ４４０
Ｍｏｄ＿Ｗ
２
Ｍｏｄ＿Ｗ
ｂ０ ９４８６５３ ２２４－３２１ ６１２０ ４２２０ ５１６
２
ｃ０ ９４８６５３ ３４４－３２１ ６７２０ ３０２０ ５８２
Ｍｏｄ＿Ｗ
２
Ｍｏｄ＿Ｗ
０ ９１９３８３ １５５－１８７ ５７７－－３
ａ０ ９２０３８４ ９８１－１８７ ４９００ １７４０ ６７７
Ｍｏｄ＿Ｗ
３
ｂ０ ９２０３８４ ９４３－１８７ ４７２０ ２１１０ ６４６
Ｍｏｄ＿Ｗ
３
ｃ０ ９２０３８４ ９０１－１８７ ４５１０ ２５４０ ６１５
Ｍｏｄ＿Ｗ
３
Ｍｏｄ＿Ｗ
０ ９６８６０２ ０９７－２９７ ０４９－－４
ａ０ ９６８６０２ ０６１－２９６ ０３０２ ０３７０ １５４
Ｍｏｄ＿Ｗ
４
ｂ０ ９６８６０２ ０６８－２９６ ０３４２ ０２９０ １５４
Ｍｏｄ＿Ｗ
４
ｃ０ ９６８６０１ ９３９－２９５ ９６９２ １５９０ １４２
Ｍｏｄ＿Ｗ
４
Ｍｏｄ＿Ｗ
０ ９３７３９０ ０９２－１９１ ０４６－－５
ａ０ ９３８３９１ ３２８－１９０ ６６４０ ７６４０ ３８２
Ｍｏｄ＿Ｗ
５
ｂ０ ９３８３９０ ９３５－１９０ ４６７１ １５７０ ２８２
Ｍｏｄ＿Ｗ
５
ｃ０ ９３８３９１ ００３－１９０ ５０２１ ０８９０ ２９７
Ｍｏｄ＿Ｗ
５
Ｍｏｄ＿Ｗ
０ ９５１６０４ ９６３－２９８ ４８２－－６
ａ０ ９５２６０４ ９８０－２９７ ４９０１ ９８３０ １５９
Ｍｏｄ＿Ｗ
６
Ｍｏｄ＿Ｗ
ｂ０ ９５２６０４ ４３５－２９７ ２１８２ ５２８０ １１２
６
ｃ０ ９５２６０４ ６４０－２９７ ３２０２ ３２４０ １２７
Ｍｏｄ＿Ｗ
６
Ｍｏｄ＿Ｗ
０ ９４３３００ ２１４－１４６ １０７－－７
ａ０ ９４５３００ ４２８－１４５ ２１４１ ７８６０ １８１
Ｍｏｄ＿Ｗ
７
ｂ０ ９４５３００ １９１－１４５ ０９５２ ０２４０ １５５
Ｍｏｄ＿Ｗ
７
Ｍｏｄ＿Ｗ
ｃ０ ９４５３００ ２７７－１４５ １３９１ ９３７０ １６４
７
Ｍｏｄ＿Ｗ
０ ９３５３５３ ２１６－１７２ ６０８－－８
Ｍｏｄ＿Ｗ
ａ０ ９３６３５４ １５５－１７２ ０７８１ ０６１０ ３０３
８
ｂ０ ９３６３５３ ９５０－１７１ ９７５１ ２６６０ ２６０
Ｍｏｄ＿Ｗ
８
ｃ０ ９３６３５４ ００４－１７２ ００２１ ２１２０ ２７１
Ｍｏｄ＿Ｗ
８
Ｍｏｄ＿Ｗ
０ ８３１２０８ ３７７－１００ １８８－－９
ａ０ ８３３２０９ ３０１－９９ ６５１１ ０７６０ ３００
Ｍｏｄ＿Ｗ
９
Ｍｏｄ＿Ｗ
ｂ０ ８３４２０８ ８１４－９９ ４０７１ ５６２０ ２１１
９
ｃ０ ８３４２０８ ５９４－９９ ２９７１ ７８３０ １８２
Ｍｏｄ＿Ｗ
９
：Ｍｏｄ＿Ｗ１为总生物量的基础模型，Ｍｏｄ＿Ｗ１ａ、Ｍｏｄ＿Ｗ１ｂ、Ｍｏｄ＿Ｗ１ｃ分别为在总生物量基础模型的参数ａ、ｂ、ｃ中引入ＳＩ后的改进模型；其它组分生物量模型命名与总生物量相同； 为Ｐ＜０ ０５，即改进模型性能显著优于基础模型；下同。

２ １ ２　生长模型的构建与评价　在参数ａ中引入ＳＩ的林分公顷胸高断面积（Ｇ）和林分平均高（Ｈ）与林龄（Ａ）之间的关系方程与基础方程相比，Ｒ２明显升高，ＡＩＣ显著降低，ｌｏｇＬｉｋ也显著升高；对数似然比检验
结果也表明，改进模型相对于基础方程差异显著（表４）。

也就是说，立地条件对Ｇ和Ｈ的生长过程有显著影响，因此，在相容性生物量模型系统中采用引入ＳＩ的改进模型作为Ｇ和Ｈ生长方程的模型形式。

２ ２　立地条件对林分胸高断面积和平均高生长过程的影响
不同立地条件林分每公顷胸高断面积（Ｇ）与林龄（Ａ）之间的关系见图１。

结果表明，Ｇ随Ａ的增大而增大，且地位指数越高，增长趋势越显著；不同立地条件林分的平均高（Ｈ）与林龄（Ａ）之间的关系也表现出相同规律，Ｈ随Ａ的增大而增大，且地位指数越高，增长趋势越显著。

立地条件对Ｇ－Ａ和Ｈ－Ａ关系的
第１期翁建宇：林龄和立地条件对杉木林乔木层生物量分配的影响
影响将间接地影响生物量分配动态变化趋势。

表４　林分公顷胸高断面积和林分平均高生长方程参数拟合结果
模型名称
Ｒ
２
ＡＩＣｌｏｇＬｉｋ最大似然比检验χ２
ＰＭｏｄ＿Ｇ０ ５４９５６７ ４６６－２７９ ７３３－－
Ｍｏｄ＿Ｇａ０ ８２１４９８ ４０１－２４４ ２００７１ ０６５０ ０００
Ｍｏｄ＿Ｈ０ ５５２４１１ ３２３－２０１ ６６１－
－Ｍｏｄ＿Ｈａ０ ９３４
２６６ ３２２
－１２８ １６１
１４７ ００１
０ ０００
：Ｍｏｄ＿Ｇ
为Ｇ与Ａ的基础模型，
Ｍｏｄ＿Ｇａ为在Ｇ与Ａ基础模型的参数ａ中引入ＳＩ后的改进模型；Ｈ与Ａ的模型命名方式同上。

图１　不同地位指数林分每公顷胸高断面积（左）、林分平均高（右）与林龄的关系曲线
２ ３　
林龄和立地条件对生物量及各组分生物量分配的影响
不同立地条件的杉木人工林生物量随着林龄的变化趋势相同。

以ＳＩ＝２０为例，林龄对生物量变化的影响见图２左；随着林龄的增长，林分总生物量和各组分生物量均随之增长；且生物量随着林龄的增长呈
Ｓ型曲线，即在幼龄林时期，林分生物量增长较缓慢，而中龄林时期增速逐渐升高，到生长后期则有减缓的趋势，总体呈现出快—慢—快的变化规律。

在相同林龄时，不同立地条件林分总生物量随着地位指数级的增加而增加，也就是说，地位指数大的林分拥有更多的生物量，且不同地位指数级林分之间的差异随林龄增加而逐渐增大（图２右）。

对比地上部分生物量（
Ｗ２
）、根生物量（Ｗ３
）占总生物量（Ｗ１
）的比值（图３）可知，在各个年龄阶段，地上生物量占比均大于根生物占比，且随着林龄的增大，地上生物量占比不断提高，根生物量则相反。

不同地位指数级的林分之间，二者占比有所差异，地位指数低的林地有着更高比例的根生物量，这可能是因为立地条件较差时，林内树木需要更加发达的根系来吸收养分和水分以供自身生长所需。

然而，随着林龄增大，不同地位指数级根生物量占比差异逐渐缩小，这就表明，幼龄林和中龄林的生物量在地上部分和根之间的分配比例更易受到立地条件的影响，而随着年龄增大，这种影响呈减缓趋势。

对比树干生物量（Ｗ４
）
、树冠生物量（Ｗ５
）占总生物量的比值（图４）可知，各地位指数级的林分中，树干生物量占比在林龄较小时增长迅速，随着林龄增长其增速逐渐减缓，为全树生物量比例最大的部分；树冠生物量占比则相反。

地位指数越大的林分，树干生物量比例越大，树冠生物量比例越小，二者占比在林分生长过程中随林龄变化较为稳定。

表明立地条件较好时，树木或林分可以利用更少的树冠支撑更多的树干生长，更加倾向于对树干进行生物量分配。

树干生物量由干材生物量（Ｗ６
）和树皮生物量（Ｗ７
）组成，其中干材占绝大部分。

由干材生物量占总
生物量的比值可知，其随林龄的变化规律与树干生物量占比相似。

然而，树皮生物量占总生物量的比例则
·
５３·
福建林业科技第４９卷
不同，总体上随林龄的增大逐渐升高，且在林龄较小时增长迅速，最终趋于不变。

地位指数级越高的林分其树皮生物量占比也越大，但随着林龄的增大，不同地位指数级林分的树皮生物量占比趋于一致，介于８ ８％～９ ２％之间见图５。

树冠生物量由树枝生物量（Ｗ８
）和树叶生物量（Ｗ９
）组成。

树枝生物量占比随林龄的变化趋势受立地条件的影响显著，在林龄较小时，地位指数高的林分其树枝生物量占比也较高，随着林龄增加先升高后下降，且地位指数越高，下降越快；最终在成熟阶段，地位指数高的林分反而树枝生物量占比更小。

树叶生物量占比随林龄增长逐渐下降，在林龄较小时，各地位指数级林分的树叶生物量占比差异较大，地位指数越高，占比越低，且随地位指数的下降逐渐上升，这种差异随着林龄的增大而逐渐缩小，但在成熟之后仍有明显差异见图６。

图２　生物量与林龄（左）、地位指数（右）的关系
图３　地上部分（Ｗ２
）和根（Ｗ３
）生物量占总生物量（Ｗ１
）的比例随林龄的变化趋势
·６３·
第１期翁建宇：林龄和立地条件对杉木林乔木层生物量分配的影响
·７３·图４　树干（Ｗ４）和树冠（Ｗ５）生物量占总生物量（Ｗ１）的比例随林龄的变化趋势
图５　干材（Ｗ６）和树皮（Ｗ７）生物量占总生物量（Ｗ１）的比例随林龄的变化趋势
图６　树枝（Ｗ８）和树叶（Ｗ９）生物量占总生物量（Ｗ１）的比例随林龄的变化趋势
３　讨论
本研究表明，杉木人工林乔木层生物量及其各组分生物量之间符合异速生长关系，基于分级控制的相容性异速生长方程组适用于生物量各组分生物量估算［６－７，９－１１］。

林龄和立地条件对树木各组分生物量的
福建林业科技第４９卷·８３·
异速生长关系影响显著，引入地位指数的相容性生物量异速生长模型系统的精度明显改善。

很多关于白松（Ｐｉｎｕｓａｒｍａｎｄｉｉ）、红松（Ｐｉｎｕｓｋｏｒａｉｅｎｓｉｓ）等林分各组分生物量估算的研究也表明，加入林龄对树枝、树叶的生物量的模型估计更加精确，与本研究的结论一致［１２－１３］。

林龄对杉木人工林乔木层各组分生物量占比的影响显著。

随着林龄的增加，地上部分生物量比例逐渐增加，而根生物量则相反，到了成熟龄阶段二者的比例则趋于稳定，这是因为处于幼龄阶段的林分需要旺盛的根系，为树木生长提供充足的水分和养分［３１］；一些对日本落叶松林的研究也表明，林分幼龄阶段树根比例较大，随着林龄的增大，根比例逐渐下降，最后稳定在２０％左右［１４－１７］。

在本研究中，树叶生物量占比随林龄的增加而逐渐降低，树枝的生物量占比随林龄增大先上升后缓慢下降。

黄志宏等［４１］在湖南会同第１代杉木人工林生物量分配动态的研究中，也发现了叶的生物量比例与林龄呈负相关，与本研究结论一致。

树叶是树木有机物质合成的器官，叶生物量占比随着林龄的下降使得有机物质的积累速率不断下降，这也验证经典理论中对树木生长率随林龄（或个体大小）的增大而逐渐下降的假设［１８－１９］。

立地质量对杉木人工林乔木层各组分生物量占比的影响显著。

地位指数越高，树干生物量所占比例越高，树叶和根生物量所占比例越低，秦建华等［４２］在研究杉木林的生物量分配变化规律中也得到相同的结果。

这就说明，立地条件好的林分，吸收大气ＣＯ２并合成有机物质的树叶和吸收养分、水分的树根能够发挥更大的潜力，以更小的占比支撑更大机体的生长。

因此，在培育杉木人工林时，特别是大径材用材林，应选择立地条件较好的林地，以获得更高的木材产量和质量。

４　结论
引入林龄和地位指数的相容性生物量模型，与基础模型相比，模型性能更高，具有更明确的生态学意义，能够更好地探究林龄和立地条件对生物量异速生长规律和分配的影响。

杉木各组分生物量的分配规律受所处的立地条件和林龄的影响显著。

当立地条件相同时，干材、树皮、生物量占比随林龄的增大而增大，树叶、树根则相反，树枝生物量占比先增加后减小。

除树枝外，立地条件不会改变各组分生物量占比随林龄的变化趋势，但会影响各组分生物量占比随林龄变化的高低。

在林龄相同时，立地条件较好的林分其树干生物量占比较大，树叶、树根的分配则相反。

因此，在估算杉木人工林乔木层生物量及其在各组分中的分配规律时，要充分考虑林龄和立地条件的影响。

参考文献：
［１］ＦＡＯａｎｄＵＮＥＰ Ｔｈｅｓｔａｔｅｏｆｔｈｅｗｏｒｌｄ’ｓｆｏｒｅｓｔｓ：ｆｏｒｅｓｔｓ，ｂｉｏｄｉｖｅｒｓｉｔｙａｎｄｐｅｏｐｌｅ［Ｒ］ Ｒｏｍｅ，２０２０
［２］ＨｏｕｇｈｔｏｎＲＡ Ａｂｏｖｅｇｒｏｕｎｄｆｏｒｅｓｔｂｉｏｍａｓｓａｎｄｔｈｅｇｌｏｂａｌｃａｒｂｏｎｂａｌａｎｃｅ［Ｊ］ ＧｌｏｂａｌＣｈａｎｇｅＢｉｏｌｏｇｙ，２００５，１１（６）：９４５－９５８ ［３］ＮｉｋｌａｓＫＪ，ＥｎｑｕｉｓｔＢＪ Ｃａｎｏｎｉｃａｌｒｕｌｅｓｆｏｒｐｌａｎｔｏｒｇａｎｂｉｏｍａｓｓｐａｒｔｉｔｉｏｎｉｎｇａｎｄａｎｎｕａｌａｌｌｏｃａｔｉｏｎ［Ｊ］ ＡｍｅｒｉｃａｎＪｏｕｒｎａｌｏｆＢｏｔａ
ｎｙ，２００２，８９（５）：８１２－８１９
［４］李鑫，欧阳勋志，刘琪瞡 江西省２００１—２００５年森林植被碳储量及区域分布特征［Ｊ］ 自然资源学报，２０１１，２６（４）：
６５５－６６５
［５］ＥｂｅｒｍｅｒｙｅｒＥ ＤｉｅｇｅｓａｍｔｅＬｅｈｒｅｄｅｒＷａｌｄｓｔｒｅｕｍｉｔＲｕｃｋｓｉｃｈｔａｕｆｄｉｅｃｈｅｍｉｓｃｈｅｓｔａｔｉｋｄｅｓＷａｌｄｂａｕｅｓ［Ｍ］ Ｂｅｌｉｎ：Ｊ Ｓｐｒｉｎｇｅｒ，
１８７６：１１６
［６］ＷｅｉｎｅｒＪ Ａｌｌｏｃａｔｉｏｎ，ｐｌａｓｔｉｃｉｔｙａｎｄａｌｌｏｍｅｔｒｙｉｎｐｌａｎｔｓ［Ｊ］ ＰｅｒｓｐｅｃｔｉｖｅｓｉｎＰｌａｎｔＥｃｏｌｏｇｙ，ＥｖｏｌｕｔｉｏｎａｎｄＳｙｓｔｅｍａｔｉｃｓ，２００４，６（４）：２０７－２１５
［７］ＺｉａｎｉｓＤ，ＭｅｎｃｕｃｃｉｎｉＭ Ｏｎｓｉｍｐｌｉｆｙｉｎｇａｌｌｏｍｅｔｒｉｃａｎａｌｙｓｅｓｏｆｆｏｒｅｓｔｂｉｏｍａｓｓ［Ｊ］ ＦｏｒｅｓｔＥｃｏｌｏｇｙａｎｄＭａｎａｇｅｍｅｎｔ，２００４，１８７（２－３）：３１１－３３２
［８］ＰｉｌｌｉＲ，ＡｎｆｏｄｉｌｌｏＴ，ＣａｒｒｅｒＭ Ｔｏｗａｒｄｓａｆｕｎｃｔｉｏｎａｌａｎｄｓｉｍｐｌｉｆｉｅｄａｌｌｏｍｅｔｒｙｆｏｒｅｓｔｉｍａｔｉｎｇｆｏｒｅｓｔｂｉｏｍａｓｓ［Ｊ］ ＦｏｒｅｓｔＥｃｏｌｏｇｙａｎｄ
Ｍａｎａｇｅｍｅｎｔ，２００６，２３７（１－３）：５８３－５９３
［９］ＧａｒｇａｇｌｉｏｎｅＶ，ＰｅｒｉＰＬ，ＲｕｂｉｏＧ ＡｌｌｏｍｅｔｒｉｃｒｅｌａｔｉｏｎｓｆｏｒｂｉｏｍａｓｓｐａｒｔｉｔｉｏｎｉｎｇｏｆＮｏｔｈｏｆａｇｕｓａｎｔａｒｃｔｉｃａｔｒｅｅｓｏｆｄｉｆｆｅｒｅｎｔｃｒｏｗｎ
ｃｌａｓｓｅｓｏｖｅｒａｓｉｔｅｑｕａｌｉｔｙｇｒａｄｉｅｎｔ［Ｊ］ ＦｏｒｅｓｔＥｃｏｌｏｇｙａｎｄＭａｎａｇｅｍｅｎｔ，２０１０，２５９（６）：１１１８－１１２６
［１０］ＧａｒｃíａＭｏｒｏｔｅＦＡ，ＬóｐｅｚＳｅｒｒａｎｏＦＲ，ＡｎｄｒéｓＭ，ｅｔａｌ Ａｌｌｏｍｅｔｒｉｅｓｂｉｏｍａｓｓｓｔｏｃｋｓａｎｄｂｉｏｍａｓｓａｌｌｏｃａｔｉｏｎｉｎｔｈｅｔｈｅｒｍｏｐｈｉｌｉｃ
第１期翁建宇：林龄和立地条件对杉木林乔木层生物量分配的影响
·９３·
ＳｐａｎｉｓｈｊｕｎｉｐｅｒｗｏｏｄｌａｎｄｓｏｆＳｏｕｔｈｅｒｎＳｐａｉｎ［Ｊ］ ＦｏｒｅｓｔＥｃｏｌｏｇｙａｎｄＭａｎａｇｅｍｅｎｔ，２０１２（２７０）：８５－９３
［１１］ＳｉｎｇｎａｒＰ，ＤａｓＭＣ，ＳｉｌｅｓｈｉＧＷ，ｅｔａｌ Ａｌｌｏｍｅｔｒｉｃｓｃａｌｉｎｇ，ｂｉｏｍａｓｓａｃｃｕｍｕｌａｔｉｏｎａｎｄｃａｒｂｏｎｓｔｏｃｋｓｉｎｄｉｆｆｅｒｅｎｔａｇｅｄｓｔａｎｄｓｏｆｔｈｉｎ ｗａｌｌｅｄｂａｍｂｏｏｓＳｃｈｉｚｏｓｔａｃｈｙｕｍｄｕｌｌｏｏａ，ＰｓｅｕｄｏｓｔａｃｈｙｕｍｐｏｌｙｍｏｒｐｈｕｍａｎｄＭｅｌｏｃａｎｎａｂａｃｃｉｆｅｒａ［Ｊ］ ＦｏｒｅｓｔＥｃｏｌｏｇｙａｎｄ
Ｍａｎａｇｅｍｅｎｔ，２０１７（３９５）：８１－９１
［１２］ＰｅｉｃｈｌＭ，ＡｒａｉｎＭＡ Ａｌｌｏｍｅｔｒｙａｎｄｐａｒｔｉｔｉｏｎｉｎｇｏｆａｂｏｖｅ ａｎｄｂｅｌｏｗｇｒｏｕｎｄｔｒｅｅｂｉｏｍａｓｓｉｎａｎａｇｅ ｓｅｑｕｅｎｃｅｏｆｗｈｉｔｅｐｉｎｅｆｏｒ
ｅｓｔｓ［Ｊ］ ＦｏｒｅｓｔＥｃｏｌｏｇｙａｎｄＭａｎａｇｅｍｅｎｔ，２００７，２５３（１－３）：６８－８０
［１３］ＫｉｎｇＪＳ，ＧｉａｒｄｉｎａＣＰ，ＰｒｅｇｉｔｚｅｒＫＳ，ｅｔａｌ Ｂｉｏｍａｓｓｐａｒｔｉｔｉｏｎｉｎｇｉｎｒｅｄｐｉｎｅ（Ｐｉｎｕｓｒｅｓｉｎｏｓａ）ａｌｏｎｇａｃｈｒｏｎｏｓｅｑｕｅｎｃｅｉｎｔｈｅ
ＵｐｐｅｒＰｅｎｉｎｓｕｌａｏｆＭｉｃｈｉｇａｎ［Ｊ］ ＣａｎａｄｉａｎＪｏｕｒｎａｌｏｆＦｏｒｅｓｔＲｅｓｅａｒｃｈ，２００７，３７（１）：９３－１０２
［１４］ＶａｎｎｉｎｅｎＰ，ＹｌｉｔａｌｏＨ，ＳｉｅｖｎｅｎＲ，ｅｔａｌ ＥｆｆｅｃｔｓｏｆａｇｅａｎｄｓｉｔｅｑｕａｌｉｔｙｏｎｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆｂｉｏｍａｓｓｉｎＳｃｏｔｓｐｉｎｅ（Ｐｉｎｕｓｓｙｌ
ｖｅｓｔｒｉｓＬ ）［Ｊ］ Ｔｒｅｅｓ，１９９６，１０（４）：２３１－２３８
［１５］ＢｒｏｗｎＳ Ｍｅａｓｕｒｉｎｇｃａｒｂｏｎｉｎｆｏｒｅｓｔｓ：ｃｕｒｒｅｎｔｓｔａｔｕｓａｎｄｆｕｔｕｒｅｃｈａｌｌｅｎｇｅｓ［Ｊ］ ＥｎｖｉｒｏｎｍｅｎｔａｌＰｏｌｌｕｔｉｏｎ，２００２，１１６（３）：３６３－
３７２
［１６］ＤｉｃｋｓｏｎＲＥ Ｃａｒｂｏｎａｎｄｎｉｔｒｏｇｅｎａｌｌｏｃａｔｉｏｎｉｎｔｒｅｅｓ［Ｊ］ ＡｎｎａｌｓｏｆＦｏｒｅｓｔＳｃｉｅｎｃｅ，１９８９，４６（Ｓ１）：６３１－６４７
［１７］ＴｏｂｉｎＢ，ＮｉｅｕｗｅｎｈｕｉｓＭ ＢｉｏｍａｓｓｅｘｐａｎｓｉｏｎｆａｃｔｏｒｓｆｏｒＳｉｔｋａｓｐｒｕｃｅ（Ｐｉｃｅａｓｉｔｃｈｅｎｓｉｓ（Ｂｏｎｇ ）Ｃａｒｒ ）ｉｎＩｒｅｌａｎｄ［Ｊ］ Ｅｕｒｏｐｅ
ａｎＪｏｕｒｎａｌｏｆＦｏｒｅｓｔＲｅｓｅａｒｃｈ，２００７，１２６（２）：１８９－１９６
［１８］ＷｅｉｎｅｒＪ Ａｌｌｏｃａｔｉｏｎ，ｐｌａｓｔｉｃｉｔｙａｎｄａｌｌｏｍｅｔｒｙｉｎｐｌａｎｔｓ［Ｊ］ ＰｅｒｓｐｅｃｔｉｖｅｓｉｎＰｌａｎｔＥｃｏｌｏｇｙ，ＥｖｏｌｕｔｉｏｎａｎｄＳｙｓｔｅｍａｔｉｃｓ，２００４，６（４）：２０７－２１５
［１９］ＩｌｏｍｋｉＳ，ＮｉｋｉｎｍａａＥ，ＭｋｅｌＡ Ｃｒｏｗｎｒｉｓｅｄｕｅｔｏｃｏｍｐｅｔｉｔｉｏｎｄｒｉｖｅｓｂｉｏｍａｓｓａｌｌｏｃａｔｉｏｎｉｎｓｉｌｖｅｒｂｉｒｃｈ［Ｊ］ ＣａｎａｄｉａｎＪｏｕｒｎａｌｏｆ
ＦｏｒｅｓｔＲｅｓｅａｒｃｈ，２００３，３３（１２）：２３９５－２４０４
［２０］唐守正，张会儒，胥辉 相容性生物量模型的建立及其估计方法研究［Ｊ］ 林业科学，２０００，３６（Ｓ１）：１９－２７
［２１］肖兴威 中国森林生物量与生产力的研究［Ｄ］ 哈尔滨：东北林业大学，２００５
［２２］胥辉，刘伟平 相容性生物量模型研究［Ｊ］ 森林与环境学报，２００１，２１（１）：１８－２３
［２３］李磊 闽北杉木生长分析及地上部分生物量研究［Ｄ］ 保定：河北农业大学，２０１４
［２４］王天博 福建将乐马尾松人工林单木生物量模型研究［Ｄ］ 北京：北京林业大学，２０１２
［２５］蔡兆伟 福建杉木人工林生物量模型研究［Ｄ］ 北京：北京林业大学，２０１４
［２６］曾伟生，唐守正 利用度量误差模型方法建立相容性立木生物量方程系统［Ｊ］ 林业科学研究，２０１０（６）：７９７－８０３ ［２７］曾伟生，唐守正 立木生物量方程的优度评价和精度分析［Ｊ］ 林业科学，２０１１（１１）：１０９－１１６
［２８］曾伟生 全国立木生物量方程建模方法研究［Ｄ］ 北京：中国林业科学研究院，２０１１
［２９］高志雄，王新杰，王廷蓉，等 福建杉木单木生物量分配及相容性模型的应用［Ｊ］ 南京林业大学学报（自然科学版），
２０１５，３９（４）：１５７－１６２
［３０］符利勇，雷渊才，曾伟生 几种相容性生物量模型及估计方法的比较［Ｊ］ 林业科学，２０１４，５０（６）：４２－５４
［３１］董利虎，李凤日，贾炜玮 东北林区天然白桦相容性生物量模型［Ｊ］ 林业科学，２０１３，４９（７）：７５－８５
［３２］孟盛旺 大兴安岭主要树种地上生物量研究［Ｄ］ 北京：北京林业大学，２０１８
［３３］李忠国，孙晓梅，陈东升，等 基于哑变量的日本落叶松生长模型研究［Ｊ］ 西北农林科技大学学报：自然科学版，
２０１１，３９（８）：６９－７４
［３４］张复礼 宿松县森林立地类型的划分和评价［Ｊ］ 安徽林业科技，２００４（２）：１０－１１
［３５］邹吉东 辽西半干旱地区困难立地条件造林技术探讨［Ｊ］ 内蒙古林业调查设计，２０１３，３６（３）：３４－３４
［３６］王光恩，林杰，吴敬，等 福建省杉木实生林数量化地位指数表编制方法的探讨［Ｊ］ 林业资源管理，１９８０（１）：２６－３７ ［３７］朱光玉，康立，何海梅，等 基于树高—年龄分级的杉木人工林多形立地指数曲线模型研究［Ｊ］ 中南林业科技大学学报，２０１７，３７（７）：１８－２９
［３８］余其芬，唐德瑞，董有福 基于遥感与地理信息系统的森林立地分类研究［Ｊ］ 西北林学院学报，２００３，１８（２）：８７－９０ ［３９］杨传桃 不同年龄杉木林分结构与生长过程［Ｄ］ 长沙：中南林业科技大学，２０１９
［４０］温远光，梁宏温，蒋海平 广西杉木人工林生物量及分配规律的研究［Ｊ］ 广西农业大学学报，１９９５（１）：５５－６４ ［４１］黄志宏，田大伦，康文星，等 湖南会同第１代杉木人工林生物量分配动态［Ｊ］ 中南林业科技大学学报，２０１１，３１（５）：
３７－４３
［４２］秦建华，姜志林 杉木林生物量及其分配的变化规律［Ｊ］ 生态学杂志，１９９６（１）：１－７。