高山松林生态系统研究进展

杜燕, 包维楷

1. 中国科学院山地生态恢复与生物资源利用重点实验室, 生态恢复与生物多样性保育四川省重点实验室, 中国科学院成都生物研究所, 四川 成都 610041；2. 中国科学院大学, 北京 100049

高山松（Pinus densataMast.）林是我国西南山区特有的森林类型，分布范围大致为北纬28°～33°，东经93°～104°，东起四川岷江流域，西迄西藏朗县，北起四川道孚，南至云南永胜，在川西地区分布海拔为2 000～3 800 m，在滇西北为3 000～3 400 m，藏东南为2 600～3 500 m[1-3]。根据第六次全国森林资源清查（1999～2003年）数据，高山松林面积为180.47 ×104 hm2[4]，其蓄积量、生物量密度和碳密度在中国森林类型中均处于较高水平[5-6]，是西南地区重要的碳库，在固碳释氧、固土保肥、水源涵养等方面发挥着重要作用[7-9]，具有较高的生态保育价值。

系统认识高山松林的群落类型与结构、生态系统服务功能等是非常有必要的，对于高山松林的科学保护与管理以及资源的合理开发利用具有重要的理论价值。系统查阅了1980～2020年发表的高山松林相关文献资料，包括文献58篇，专著5本。综述如下方面的研究现状：（1）高山松林群落类型与特征；（2）群落生物量与生产力；（3）生态系统水源涵养能力；并进一步梳理出当前需要聚焦的方向，为高山松林的深入研究提供基础。

1 群落类型与特征

1.1 群落类型

以《中国植被志》研编内容与规范[10]为群落分类标准，对已有资料进行梳理，将高山松林划分为6个群系9个群丛组6个群丛（见表1）。除上述群落类型外，高山松还可与丽江云杉（Picea likiangensis（Franch.） E. Pritzel）、川西云杉（Picea likiangensisvar. rubescensRehder & E. H. Wilson）、白桦（Betula platyphyllaSuk.）、帽斗栎（Quercus guajavifoliaH.Léveillé）等形成混交林，灌木层优势种还包括云南杜鹃（Rhododendron yunnanenseFranch.）、大白杜鹃（R. decorumFranch.）、腋花杜鹃（R. racemosumFranch.）、圆锥山蚂蟥（Desmodium elegansCandolle）等，草本层优势种还包括秦岭槲蕨（Drynaria baroniiDiels）、金茅（Eulalia speciosa（Debeaux） Kuntze）等[1, 14-16, 19]。

1.2 群落结构特征

少量研究对群落优势种高山松的胸径、高度和冠幅结构进行了分析。卢杰等[20]对林芝八一镇高山松种群调查发现，其径级结构、高度结构和冠幅结构均呈反“J”形，为增长型种群。在林芝地区高山松天然次生林中，幼龄林的胸径分布呈偏左的近似正态分布，随着林龄的增加，偏度和峰度均变小，中龄林为中间高两边低的近似正态分布[21]。对西藏林芝地区高山松种群格局的分析结果表明，其在幼苗-幼树-立木的发育过程中，空间格局由集群分布转变为随机或均匀分布[22-23]。左政等[24]对香格里拉的高山松林直径结构分析发现，其林分直径结构以迈耶负指数分布函数的拟合效果最好。

表 1 高山松林群落分类Tab. 1 The classification of Pinus densata forest

2 群落生物量与生产力

2.1 群落生物量与碳储量

在林分水平上，少量研究以蓄积量或林龄为自变量，建立了高山松林乔木层生物量回归模型（见表2）；还有多个研究以香格里拉地区的高山松林为研究对象，构建了高山松林地上生物量遥感估测模型[28-33]，但精度较低，均不超过80%。对高山松林群落生物量已有调查数据进行收集（见表3），结果表明，高山松成熟林群落生物量为81.24～318.79 t·hm-2，其中乔木层为79.39～311.53 t·hm-2，灌木层为0.105～1.49 t·hm-2，草本层为0.035～0.3318 t·hm-2，凋落物层为0.515～90.84 t·hm-2。群落细根生物量的研究仅1个，何永涛等[42]测定了林芝八一镇高山松林0～50 cm土层中细根生物量，为431.2 g·m-2，其中活细根生物量为326.1 g·m-2，死细根为105.1 g·m-2；在垂直分布上，69.6%的活细根生物量集中在0～10 cm，而40～50 cm土层内无细根分布。

对高山松林碳储量的研究包括树种含碳率和群落碳密度两个方面。王金亮等[43]测定了香格里拉高山松林不同林龄不同器官的含碳率，树种平均含碳率为51.31%，与张坤[44]基于1978—1994年期间关于中国生物量的文献和全国第三次森林资源清查（1984—1988年）数据的计算结果接近，为50.09%。杨阳[37]测定了林芝地区不同林龄高山松林的草本和凋落物含碳量，幼龄林、中龄林和成熟林草本地上含碳量分别为441.44 g·kg-1，412.84 g·kg-1和433.30 g·kg-1；草本地下含碳量分别为389.12 g·kg-1，370.54 g·kg-1和393.12 g·kg-1；凋落物含碳率分别为500.90 g·kg-1，507.35 g·kg-1和480.37 g·kg-1。对高山松林群落碳储量已有调查数据进行收集，结果表明，高山松林乔木层碳密度为49.543～103.24 t·hm-2，灌木层碳密度为10.964 t·hm-2，草本层碳密度为0.868 t·hm-2，凋落物碳密度为1.43～18.68 t·hm-2，土壤层（0～100 cm）碳密度为216.274 t·hm-2[39,45-50]。

表 2 高山松林生物量密度估算模型Tab. 2 Biomass density estimation model of Pinus densata forest

表 3 高山松林群落生物量（t·hm-2）Tab. 3 Community biomass of Pinus densata forest

2.2 群落生产力

高山松林生产力的相关研究很少，仅3篇。吴兆录等[51]建立了高山松林乔木器官净第一性生产力的优化回归模型，对香格里拉吉迪林场林龄40年和100年的高山松林生产力进行了估算，结果分别为12.192 t·hm-2·a-1和10.013 t·hm-2·a-1，其中乔木层分别为12.160 t·hm-2·a-1和9.980 t·hm-2·a-1，主要分配在 叶 和 树 干 中；灌 木 层 分 别 为0.018 t·hm-2·a-1和0.021 t·hm-2·a-1，草 本 层 分 别 为0.014 t·hm-2·a-1和0.012 t·hm-2·a-1。唐佳[52]估算出工布自然保护区高山松林年生产力为9.85 t·hm-2·a-1。基于四川、云南和西藏17个样地数据，高山松成熟林全林和乔木层生产力分别为5.48～18.07 t·hm-2·a-1和4.29～14.23 t·hm-2·a-1；在地理分布格局上，包括高山松林在内的温性松林生物生产力随纬度的增加而递减，随经度的增加而递增；而在水热分布格局上，呈现一种双曲面的递增函数，随降水量的增加呈自然对数递增，随温度的增加呈线性递增，随温暖指数和潜在增散量的增加也呈递增趋势[34]。

3 水源涵养能力

3.1 林冠截留量

高山松林林冠截留量仅在色季拉山和纳帕海开展过研究，结果表明高山松林林冠截留率为24.32%～28.87%，且林冠截留量与降雨量间以幂函数拟合效果最好[53-55]。

3.2 苔藓层持水特性

高山松林林下苔藓层持水特性研究仅1例，陈甲瑞和王小兰[56]对色季拉山东坡高山松林的调查发现，林下苔藓层平均厚度为1.13 cm，生物量为1.02 t·hm-2，自然含水量为1.83 t·hm-2，最大持水量为8.69 t·hm-2；浸水实验表明，苔藓在前0.25 h内吸水速率最大，可达654.279 g·kg-1·h-1，随后迅速下降，2 h后下降速度减缓，8 h后吸水速率基本稳定并趋向于零，吸水速率与浸泡时间呈幂函数关系。

3.3 凋落物持水特性

已有研究对色季拉山成熟林以及纳帕海幼龄林和中龄林凋落物层持水能力和吸水特性进行了测定（见表4和表5），结果表明高山松成熟林凋落物层最大持水量为117.27 t·hm-2；浸水实验发现，在第5 min时吸水速率最大，可达6 210 g·kg-1·h-1，随后迅速下降，8～10 h后吸水量基本达到最大值，吸水速率与浸泡时间呈幂函数关系[59]。

3.4 土壤持水特性

土壤层持水特性的研究集中在林芝地区和纳帕海（见表6），结果表明高山松成熟林0～30 cm土层的土壤最大持水量为380.98 t·hm-2[61]。

表 4 高山松林凋落物层持水量、持水率和有效拦蓄量Tab. 4 Water-holding capacity, water-holding rate, and effective interception of litter layer in Pinus densata forest

表 5 高山松林凋落物层持水量、吸水速率和持水率的拟合方程Tab. 5 Fitting equation of water-holding capacity, water absorption rate and water-holding rate of litter layer in Pinus densata forest

表 6 高山松林土壤层最大持水量、持水率和拦蓄量Tab. 6 Water-holding capacity, water-holding rate, and water storage capacity of soil layer in Pinus densata forest

4 结论与展望

4.1 群落类型与结构

样方调查数据是认识群落类型和结构的基础资料，现有调查资料多形成于上世纪50—70年代，近期的调查集中在林芝市、香格里拉县等高山松林集中分布区；调查对象主要针对高山松纯林，对混交林的关注较少；研究空白点较多，基础数据严重缺乏。在群落类型的划分上，当前研究多基于群落外貌对群落类型进行划分，仅罗建[13]采用了数量分类方法，分类和命名标准不统一，无法形成完整的高山松林群落分类系统。在物种组成方面，已有研究仅列出各层片主要组成物种，缺乏对物种科属统计、生活型和地理区系成分特征的分析。在群落结构上，仅在胸径结构方面有少量研究，垂直结构[20]仅见1个报道。因此，当前对高山松林群落类型、生境条件、物种组成、群落结构等的认识是不全面的，后续需要对高山松林开展野外调查，明确其分布范围和地理分布格局，形成完整的群落分类系统，补全群落特征的相关认识。

4.2 群落生物量与生产力

森林生物量、碳储量和生产力是研究森林生态系统结构和功能的重要基础数据，目前针对高山松林生物量、碳储量和生产力的研究集中在乔木层，且调查数据集中在林芝地区和香格里拉市，目前仅有1个研究[19]对高山松林各层片的生物量数据均进行了测算。针对高山松林这一群落类型的生物量和生产力空间分布格局和影响因素未见报道。因此，后续应对生物量数据缺失的区域开展补充调查，重点关注林下植被和地下部分的生物量测定；阐明高山松林生物量和生产力在各地理梯度上的变化规律及其驱动因子。

4.3 水源涵养能力

涵养水源是森林生态系统的重要功能之一，当前高山松林水源涵养的研究集中在林芝地区和香格里拉市；研究对象以幼龄林和中龄林居多，成熟林水源涵养能力的数据严重缺乏；层次上以林冠层、凋落物层和土壤层的调查为主，苔藓层相关研究仅1例[56]，已有研究对土壤层的调查深度标准不一。因此，后续应对成熟林水源涵养能力开展补充调查，重点关注土壤层和苔藓层的持水能力特性。