汪晓红　谯媛媛　亓新海　张乐乐　解小锋　韩敬　付德刚

摘 要 为提升赤松人工林质量，发挥其生态效益，以山东省济南市莱芜区国有华山林场的赤松人工林为研究对象，以胸径分布、树高、枝下高、冠幅、生物量及其分配为指标，探究林分密度对赤松人工林生长特征的影响。结果表明，1）当林分密度为600株·hm-2、900株·hm-2时，树木株数随着树木胸径的增大先增加后减少；当林分密度超过1 200株·hm-2时，小径级树木株数随着林分密度的增大而增加。2）随着林分密度的增大，赤松人工林树高表现为先升高后降低的变化趋势，当林分密度为1 200～1 500株·hm-2时，赤松人工林平均树高最高，林木之间的竞争效应最强；赤松人工林第一活枝高随林分密度的增大而升高，第一死枝高变化趋势不明显。3）林分密度越大，赤松单株冠幅越小，光合生产力越低，树木衰退现象越明显。4）赤松人工林生物量的分配受林分密度的制约，其中赤松单株总生物量随着林分密度的增大而显著减少，根生物量随林分密度的增大而增加。

关键词 赤松；林分密度；生长特征

中图分类号：S718.4 文献标志码：A DOI：10.19415/j.cnki.1673-890x.2023.17.061

林分密度是制约林木获取资源数量的重要因子，不仅影响树木的胸径、树高、冠幅等生长指标，还影响树木单株生物量的分配策略[1-3]。合理的林分密度不仅能改善林木生长环境、提高树木生产力及物种多样性，还有助于提高森林的碳汇能力，对提升森林生态系统功能有重要意义[4-5]。

林木的生长指标主要包括胸径、树高、生物量、蓄积量等。有学者研究发现，森林植被密度对森林乔木层的生长指标具有重要影响[6-7]。Woeste等研究表明，低林分密度下美国黑核桃（Juglans nigra）的树高和地径明显高于中高林分密度[8]；杜振宇等发现，在林分密度较低时，侧柏林木的胸径和树高均有大幅提高[9]；李洁等研究表明，银合欢（Leucaena leucocephala）的林分密度越大，其平均胸径越大[10]；张雷发现，林分密度可对树木的胸径、生物量产生影响[11]。因此，研究林木生长特征对林分密度的响应规律，对合理确定不同林木的林分密度，以及指导森林经营管理具有重要意义。

赤松（Pinus densiflora）为常绿乔木，树皮呈橘红色，高可达30 m，其因树干可做木材、种子能食用、叶可提取芳香油等特点，在我国多地备受喜爱[12]。赤松具有抗风力强、耐贫瘠等特点，是山东省森林生态系统中重要的优势碳汇树种之一[13-14]。但由于大部分赤松人工林经营管理粗放、抚育不及时，出现林分结构不合理、生长不良等问题，严重影响赤松人工林多种效应的发挥，急需对其进行改造提升[15]。为此，本文以赤松人工林为研究对象，探究赤松响应林分密度变化的生长特征，为赤松人工林林分结构改造和可持续发展与经营提供科学依据。

1  材料与方法

1.1  研究区概况

研究区位于山东省济南市莱芜区国有华山林场（北纬36°22′～36°28′、东经117°22′～117°33′），地处鲁中山区，泰山东麓，属温带大陆性季风气候，年平均降水量800 mm，有乔木50科200多种。

选取的赤松人工林，林龄62年，林分为公益林，禁伐。该林区土壤为棕壤，土层厚度在10～30 cm，表层土多沙砾。随着时间的推移，受自然稀疏、人为间伐和旅游开发等因素干扰，林分密度逐渐产生差异。

1.2  样地设置及数据调查

2022年5月，通过典型样地法选取现存林分密度分别为（600±60）株·hm-2、（900±90）株·hm-2、（1 200±90）株·hm-2、（1 500±90）株·hm-2、（1 800±90）株·hm-2的5片赤松人工林，分别命名为D600、D900、D1200、D1500、D1800，每种林分密度各设置4个25.8 m

×25.8 m的标准样地，并对样地内的乔木进行每木检尺，记录每株树木的胸径、冠幅、树高等指标。样地基本情况如表1所示。

1）树高测定。用树高仪（CGQ-1）测定树木地面根茎到树梢之间的高度。2）胸径测定。用胸径尺測定距离地面1.3 m处的树干直径。3）冠幅测定。用米尺分别测定树冠南北方向和东西方向的宽度，并取其平均值。

1.3  乔木层生物量测定

采用平均标准木法，在每个标准样地内选择生长正常，能够代表平均生长状况的赤松作为测定乔木生物量的样木。将样木伐倒、挖根，烘干至恒重，测定树干、枝、叶、根的生物量。

1.4  数据统计分析

采用SPSS 18.0统计分析软件对数据进行单因素方差分析（One-Way ANOVA），组间比较采用最小显著差异法（LSD），在5%水平上检验不同样地测定结果的差异性。

2  结果与分析

2.1  不同林分密度下赤松人工林胸径分布

如图1所示，不同林分密度下，赤松人工林胸径分布趋势不尽相同。其中，林分密度为600株·hm-2、900株·hm-2时，树木株数随着胸径的增大呈先增加后减少的分布趋势；D600组中，树木的胸径（单位为cm）主要分布在[16，18），占组内树木总株数的22.9%；D900组中，胸径处于[10，12）组的赤松株数最多，[12，14）次之，分别占组内树木总株数的18.4%、15.8%。当林分密度超过1 200株·hm-2时，小径级（胸径小于10 cm）树木株数随着林分密度的增大而增加，其中胸径小于4 cm的赤松株数分别占其组内树木总株数的17.3%（D1200）、18.5%（D1500）、25.0%（D1800）。

总体来看，林分密度越大，赤松人工林中小径级的树木株数越多；树木株数随着胸径大小的分布趋势与天然杂木林的分布曲线存在一定的差异[16-17]。以上数据说明，该赤松人工林一定程度上存在人为干扰。

2.2  林分密度对赤松人工林树高和枝下高的影响

树高是指树木从地面上根茎到树梢之间的高度，是衡量林分生长最重要因子之一，同一种树木在不同的生长环境下树高不同[18-20]。如图2所示，随着林分密度的增大，赤松人工林树高表现出先升高后降低的变化趋势，其中林分密度为1 200株·hm-2时，赤松人工林平均树高最高，为8.9 m，与D1500之间差异不显著（p＞0.05），但显著高于D600、D900、D1800（p＜0.05）。当林分密度为600株·hm-2时，赤松平均树高最低，分别为D900、D1200、D1500、D1800的90.5%、83.6%、86.6%、91.4%，显著低于其他林分密度（p＜0.05）。由此可见，林分密度是影响赤松人工林树高的重要因素。

枝下高在一定程度上可反映树冠的衰退情况，当树木开始衰退时，底部的枝条会逐渐死亡，第一活枝高逐渐升高[21]。如图2所示，赤松人工林第一活枝高随林分密度的增大而升高，其中D1200、D1500、D1800平均第一活枝高分布在3.7～4.0 m，各林分密度之间差异不显著（p＞0.05），但显著高于D600（p＜0.05）。由此可见，只有林分密度变化较大时，才会引起赤松人工林第一活枝高发生显著变化。

第一死枝高随林分密度的增大无明显变化规律，其中D1200组赤松人工林平均第一死枝高最高，但与D900、D1500、D1800差异不显著（p＞0.05）。D600第一死枝高最低，分别为D900、D1200、D1500、D1800的93.2%、85.4%、88.6%、85.8%。

总体而言，林分密度增大，树木拥挤程度加剧，林木为了获取更多的光照资源而趋向有利于自身生长的方向生长。

2.3  林分密度对赤松人工林冠幅的影响

如图3所示，随着林分密度的增大，赤松人工林冠幅呈现减小的变化趋势。其中，林分密度为600株·hm-2的赤松人工林平均冠幅最大，为5.1 m；林分密度為1 800株·hm-2的冠幅最小，仅为D600的71.0%。通过单因素方差分析可以看出，D900和D600、D1200之间差异不显著（p＞0.05），但显著大于D1500和D1800（p＜0.05）。由此可见，林分密度对赤松人工林冠幅产生较大影响，林分密度越大，林木冠幅越小，光合生产力越低，树木衰退现象越明显。

2.4  林分密度对赤松单株生物量分配的影响

植物生物量分配不仅与植物的遗传发育有关，还受到生长环境的限制[22]，其中生物量在不同器官间的分配是植物对环境变化响应的重要生理过程[23]。如表2所示，赤松单株总生物量随着林分密度的增大而减少，各林分密度之间差异显著（p＜0.05）。由此可见，林分密度是影响赤松人工林生物量的重要环境因素之一。

根系是植物吸收水分和养分的重要器官，当植物营养受限时，往往会增大根的生物量[24-25]。如表2所示，根生物量占比（RMR）随林分密度的增大而升高，其中D1800的RMR最大，分别为D600、D900、D1200、D1500的1.9倍、1.7倍、1.6倍、1.1倍。观测数据在一定程度上说明林分密度越大，林木之间的资源竞争越激烈，赤松将更多的物质输送到根部以促进根的生长。

总体而言，赤松生物量的分配受林分密度的制约，林分密度越大，赤松分配到地下的生物量越多，地上的生物量越少，这也是赤松适应林分密度变化的策略之一。

3  结论与讨论

胸径和树高不仅是森林资源调查和监测中的关键内容，还是反映森林生长状况、生产力和生态功能等方面的重要参数指标[18-19]。树高和胸径的关系决定了树高曲线在计算材积、生物量、碳储量、立地指数、林分生长和收获等方面具有不可替代的作用[26]。本研究中，赤松人工林胸径随林分密度的变化趋势不尽相同。当林分密度为600株·hm-2、900株·hm-2时，树木株数随着树木胸径的增大呈现先增加后减少的分布趋势；当林分密度超过1 200株·hm-2时，小径级树木株数随着林分密度的增大而增加。树木株数随着胸径大小的分布趋势与天然杂木林的分布存在差异，分析可能与赤松人工林受到人工干扰有关。

林分密度越大，树木的生长空间越拥挤，林木为了获取更多的资源而趋向有利于自身生长的方向生长[27]。本研究发现，赤松人工林树高随着林分密度的增大表现为近似抛物线状的分布趋势。由此可见，随着林分密度的增大，树木之间的竞争作用增强，当林分密度为1 200～1 500株·hm-2时，赤松人工林树高之间的竞争效应最强；而当林分密度为1 800株·hm-2时，赤松人工林树高降低，分析原因可能是个体营养分配趋向根部而造成的。

生物量是绿色植物通过光合作用积累的物质和能量，是反映植物生产力的重要指标[28]。本研究中，赤松单株生物量随着林分密度的增大而显著减少，可见林分密度是影响赤松人工林生物量的重要环境因素之一。有学者研究发现，在弱光照条件下，植物将更多的资源分配给茎、叶的生长，而在低营养条件下，植物往往将资源用以提高根的产量[24]。本研究发现，根生物量占比随林分密度的增大而降低，说明赤松人工林分密度越大，植株可获取的营养越少，赤松将更多的有机物分配到地下以促进根的生长。

参考文献：

[1] 黎磊，周道玮，盛连喜.密度制约决定的植物生物量分配格局[J].生态学杂志，2011，30（8）：1579-1589.

[2] ANNING A K， MCCARTHY B C. Competition， size and age affect tree growth response to fuel reduction treatments in mixed-oak forests of Ohio[J].Forest Ecology and Management，2013，307：74-83.

[3] 郑超凡，周驰宇，贺斌，等.云南松生物量特征与苗木质量指数间相对生长关系的分析[J].山西农业大学学报（自然科学版），2023，43（2）：121-128.

[4] 方精云，王襄平，沈泽昊，等.植物群落清查的主要内容、方法和技术规范[J].生物多样性，2009，17（6）：533-548.

[5] 张阳锋，尹光天，杨锦昌，等.造林密度对米老排人工林初期生长的影响[J].林业科学研究，2018，31（4）：83-89.

[6] 沈国舫.森林培育学[M].北京：中国林业出版社，2001.

[7] 潘文婷，夏莘，夏良放，等.造林密度对近熟期鹅掌楸生长和材质的影响[J].南京林业大学学报（自然科学版），2018，42（5）：46-52.

[8] WOESTE K E， JACOBS D F， MCKENNA J R. Half-sib seed source and nursery sowing density affect black walnut （Juglans nigra） growth after 5 years[J].New Forests，2011，41（2）：235-245.

[9] 杜振宇，葛忠强，王清华，等.鲁中山地侧柏人工林土壤质量评价及林分密度的影响作用[J].水土保持研究，2020，27（2）：128-135.

[10] 李洁，列志旸，许松葵，等.不同密度的银合欢林生长分析[J].中南林业科技大学学报，2016，36（6）：70-74.

[11] 张雷.祁连山青海云杉林结构和树木生长随海拔的变化[D].北京：中国林业科学研究院，2015.

[12] 杜华.昆嵛山赤松林不同林型结构特征与生产力的研究[J].北京林业大学学报，2012，34（1）：19-24.

[13] 马琳，梁军，周国英，等.昆嵛山天然赤松次生林空间结构特征的研究[J].中南林业科技大学学报，2013，33（9）：50-54.

[14] 孟鹏，安宇宁，白雪峰.沙地赤松光合及叶绿素a快相荧光动力学特性[J].生态学报，2016，36（11）：3469-3478.

[15] 孫志强，张星耀，林琳，等.赤松纯林林分特征对昆嵛山鳃扁叶蜂发生量的影响[J].生态学报，2010，30（4）：857-866.

[16] 朱文泉，何兴元，陈玮，等.城市森林结构的量化研究：以沈阳树木园森林群落为例[J].应用生态学报，2003（12）：2090-2094.

[17] 李辉，谢会成，赵春仙，等.济南市城市森林结构特征分析[J].西北林学院学报，2013，28（2）：213-217.

[18] KRUMLAND B E， WENSEL L C. A generalized height-diameter equation for coastal California species[J].Western Journal of Applied Forestry，1988，3（4）：113-115.

[19] 曾翀，雷相东，刘宪钊，等.落叶松云冷杉林单木树高曲线的研究[J].林业科学研究，2009，22（2）：182-189.

[20] 何春霞.树高极限的水分限制机理研究[D].北京：北京林业大学，2009.

[21] 吕乐，董利虎，李凤日.黑龙江省东部地区天然椴树单木冠幅预测模型[J].东北林业大学学报，2019，47（7）：37-42.

[22] 丁俊祥，范连连，李彦，等.古尔班通古特沙漠6种荒漠草本植物的生物量分配与相关生长关系[J].中国沙漠，2016，36（5）：1323-1330.

[23] YUE K， FORNARA D A， LI W， et al. Nitrogen addition affects plant biomass allocation but not allometric relationships among different organs across the globe[J].Journal of Plant Ecology， 2021，14（3）：361-371.

[24] BLOOM R G， MALLIK A U. Indirect effects of black spruce （Picea mariana） cover on community structure and function in sheep laurel （Kalmia angustifolia） dominated heath of eastern Canada[J].Plant and Soil， 2004，265：279-293.

[25] TRUMBORE S E， GAUDINSKI J B. The secret lives of roots[J].Science，2003，302（5649）：1344-1345.

[26] 樊伟，许崇华，崔珺，等.基于混合效应的大别山地区杉木树高−胸径模型比较[J].应用生态学报，2017，28（9）：2831-2839.

[27] RIOFIO J， MIREN R， MAGUIRE A D， et al. Species mixing effects on height-diameter and basal area increment models for Scots pine and Maritime pine[J].Forests，2019，10（3）：249.

[28] WEINER J. Allocation， plasticity and allometry in plants[J].Perspectives in Plant Ecology， Evolution and Systematics，2004，6（4）：207-215.

（责任编辑：张春雨  丁志祥）

收稿日期：2023-07-09

基金项目：山东省财政资金项目“国有林场森林质量精准提升技术研究和示范工程（国有林场森林抚育技术研究）”（SDGP370000000202302001097）。

作者简介：汪晓红（1977—），女，山东莒县人，本科，高级工程师，研究方向为林木生理生态。

*为通信作者，E-mail：hanjing@lyu.edu.cn。