生态学报  2022, Vol. 42 Issue (23): 9726-9735

文章信息

冯相艳, 赵文智, 蔺鹏飞, 王川
FENG Xiangyan, ZHAO Wenzhi, LIN Pengfei, WANG Chuan
祁连山北坡主要木本植物功能性状及其海拔分异
The main woody plant functional traits and altitude differences on the northern slope of Qilian Mountains
生态学报. 2022, 42(23): 9726-9735
Acta Ecologica Sinica. 2022, 42(23): 9726-9735
http://dx.doi.org/10.5846/stxb202110212974

文章历史

收稿日期: 2021-10-21
网络出版日期: 2022-07-27
祁连山北坡主要木本植物功能性状及其海拔分异
冯相艳1,2 , 赵文智1 , 蔺鹏飞1 , 王川1,2     
1. 中国科学院西北生态环境资源研究院, 中国生态系统研究网络临泽内陆河流域研究站, 中国科学院内陆河流域生态水文重点实验室, 兰州 730000;
2. 中国科学院大学, 北京 100049
摘要: 植物功能性状是植物对生境长期适应所形成的可测量特征, 受植物遗传特性和环境因子的共同影响。祁连山是我国干旱区的一条重要山脉, 北坡海拔介于2000—5000 m, 分布着青海云杉、祁连圆柏等10余种木本植物。深入了解祁连山区木本植物功能性状随海拔梯度分异特征对认识山地植物的适应性和植被垂直地带分布具有重要意义。选取祁连山北坡天然分布的11种木本植物, 其中灌木9种, 分别是金露梅(Potentilla fruticosa)、银露梅(Potentilla glabra)、水栒子(Cotoneaster multiflorus)、猫儿刺(Ilex pernyi)、吉拉柳(Salix gilashanica)、鬼箭锦鸡儿(Caragana jubata)、刚毛忍冬(Lonicera hispida)、高山绣线菊(Spiraea alpina)和鲜黄小檗(Berberis diaphana), 乔木2种, 分别是青海云杉(Picea crassifolia)和祁连圆柏(Juniperus przewalskii), 调查了木质密度、胡伯尔值、叶面积、叶干物质含量、叶含水率和比叶面积6个枝、叶功能性状。结果表明:(1) 灌木的木质密度、胡伯尔值、叶含水率和比叶面积显著高于乔木, 但叶面积和叶干物质含量显著低于乔木;(2)木质密度和叶干物质含量随海拔升高显著下降, 其他功能性状沿海拔分异不显著;(3)生长季气温与叶干物质含量正相关, 与胡伯尔值和叶含水率负相关, 生长季降水与各功能性状相关关系不显著。说明乔木和灌木植物采取不同的枝、叶资源利用策略来适应高山生境;随海拔升高, 降低结构性碳投资可能限制了高山植被的分布上限;生长季气温是引起性状沿海拔分异的重要环境因子, 但不能完全解释性状沿海拔分异规律。
关键词: 祁连山    植物功能性状    海拔    植被分布    木本植物    
The main woody plant functional traits and altitude differences on the northern slope of Qilian Mountains
FENG Xiangyan1,2 , ZHAO Wenzhi1 , LIN Pengfei1 , WANG Chuan1,2     
1. Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Linze Inland River Basin Research Station, Key Laboratory of Inland River Basin Ecohydrology, Lanzhou 730000, China;
2. University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: Plant functional traits are the measurable characteristics of plants which are formed by long-term adaptation to habitat environmental conditions. The characterization of plant functional traits is mainly influenced by both the plants specific genetic basis and the shifted environmental factors. Qilian Mountains, with an elevation ranging from 2000 to 5000 m, is one of the major mountains in the dry land in Northwest China, where naturally distributed more than 10 kinds of woody plants such as Picea crassifolia and Juniperus przewalskii on the northern slope. It is of great academic and practice significance to understand the acclimation strategy and the vertical distribution pattern of alpine woody plants by deeply investigating the woody plant functional traits variation and the behind driving mechanism along with the altitude gradients in the Qilian Mountains. In the present study, we investigated 11 naturally distributed woody plant species along with the altitude in the Qilian Mountains, which included 9 shrub species (Potentilla fruticosa, Potentilla glabra, Cotoneaster multiflorus, Ilex pernyi, Salix gilashanica, Caragana jubata, Lonicera hispida, Spiraea alpine, and Berberis diaphana, respectively) and 2 tree species (Picea crassifolia and Juniperus przewalskii, respectively). A total of six plant functional traits of shoot and leaf were investigated in this study, and which include woody density, huber value, leaf area, leaf dry matter content, leaf water content, and specific leaf area, respectively. The results showed that: (1) woody density, huber value, leaf water content, and specific leaf area of shrub species were significantly greater than those of tree species, and leaf area and leaf dry matter content of shrub species were significantly lower than those of tree species. (2) Woody density and leaf dry matter content declined significantly along with the altitude gradients increased, and the other functional traits changed insignificantly along with the altitude gradients increased. (3) Growing-season air temperature was positively correlated with leaf dry matter content, and was negatively correlated with huber value and leaf water content; however, non-significant relationships between growing-season precipitation and plant functional traits were found. The results showed that tree and shrub plants adapted to the alpine habitats by adopting different branch and leaf resources utilization strategies; the decreasing structure carbon investment with the altitude gradients increased possibly has limited the upshift of alpine plants; the growing-season air temperature was one of the key environmental factors tuning woody plant functional traits variation pattern along with the altitudes gradients but could not explain all the variance.
Key Words: Qilian Mountains    plant functional trait    altitude    plant distribution    woody plants    

在山地生态系统中, 随海拔升高环境因子会形成明显的变化梯度, 如温度下降、大气压降低、光照增加以及降水波动加大等[14], 植物也会形成相应的功能性状特性来适应这种变化[57]。研究表明, 沿海拔升高比叶重增加[89]、叶氮含量降低[10]、光合速率下降[11]、木质部导管变细[12]。其中, 比叶重增加和木质部导管变细能够增强植物的低温抵抗能力[1314]。多数研究认为, 气温和降水是引起植物功能性状沿海拔分异的主要环境因子, 但具体的影响程度仍存在争议[1, 1517]。这是由于气温沿海拔呈线性下降趋势, 而降水沿海拔呈“驼峰型”变化趋势。因此, 无论是比叶重、叶氮含量还是导管直径沿海拔分异都与气温密切相关, 而与降水相关性较弱[10]

祁连山位于青藏高原、蒙新高原和黄土高原交汇带, 是我国季风和西风交汇的敏感区, 同时也具有典型的大陆性和高原气候特征[18]。山区海拔起伏大, 水热条件差异显著, 植被沿海拔梯度存在明显的地带性分布格局[1921]。沿海拔梯度, 断续分布着30余种草本植物和10种木本植物, 其中木本植物中灌木以金露梅为优势种, 乔木以青海云杉为优势种。由于适应能力差异, 不同物种分布于不同的海拔区间, 如青海云杉主要分布于2200—3300 m之间, 而金露梅主要分布于2200—3700 m之间。植被调查发现, 大多数植物适宜分布在海拔2700 m左右, 随着海拔继续上升, 物种数量与种群大小均显著减少[22]。已有研究从气候因素、生物竞争, 土壤养分和生理限制的角度分别探究了植被沿海拔分布格局的形成机制, 结果表明, 生长季土壤温度、大气温度、生物竞争、生理过程、土壤氮库与植被沿海拔分异密切相关[2325]。然而, 针对植物功能性状沿海拔分异及其对植被空间分布格局形成的影响却鲜见报道[2627]。前人分别在天山[26]和青海湖流域[27]开展了植物功能性状沿海拔分异的研究。上述研究均关注了草本植物功能性状沿海拔分异特征, 缺乏专门针对木本植物功能性状沿海拔分异的研究。因此, 在典型山地生态系统开展木本植物功能性状沿海拔分异将会加深对木本植物功能性状沿海拔分异的理解及其对植被垂直地带分布格局形成的潜在影响。

本研究假设在异质性如此明显的垂直地带分布格局下, 不同海拔分布的木本植物必然会形成特定的功能性状特性以适应生境, 选取了祁连山北坡沿海拔分布的11种具有代表性的木本植物作为研究对象, 测定了各物种在不同海拔的枝、叶的木质密度、胡伯尔值、叶面积、叶干物质含量、叶含水率和比叶面积6个功能性状。通过比较不同木本植物之间功能性状差异、各功能性状沿海拔分异规律以及功能性状与环境之间的相互关系, 以期从植物功能性状分异角度解释高山植被地带分布格局的形成机制。旨在较深入理解不同生长型木本植物适应高山生境的功能性状特性, 各植物功能性状沿海拔的分异规律, 引起功能性状沿海拔分异的主要环境因子等3个问题。

1 材料和方法 1.1 研究区概况

研究区位于祁连山中段大野口流域, 地理坐标为100°17′E—100°18′E, 38°31′N—38°33′ N, 海拔范围为2300—3800 m。研究区年均气温为0.4 ℃, 最低温-30℃, 最高温36 ℃;年降水量415.8 mm, 季节性差异现明显, 约80%的降水发生在5—10月, 每年的11月至次年4月为结冰期;年蒸发量1011.3 mm。土壤类型主要为森林灰褐土, 土壤pH值介于7—8, 土壤厚度介于0.3—2 m, 成土母质主要是钙质岩、泥炭岩和砾岩等。植被以乔木、灌木和草本植物为主, 乔木以青海云杉(Picea crassifolia)为主, 混杂有少量祁连圆柏(Juniperus przewalskii);灌木物种相对丰富, 主要包括鬼箭锦鸡儿(Caragana jubata)、金露梅(Potentilla fruticosa)、银露梅(Potentilla glabra)、高山绣线菊(Spiraea alpina)、吉拉柳(Salix gilashanica)、鲜黄小檗(Berberis diaphana)、水栒子(Cotoneaster multiflorus)、刚毛忍冬(Lonicera hispida)和猫儿刺(Ilex pernyi);草本主要为珠芽蓼(Polygonum viviparum)、针茅(Stipa capillata)、马蔺(Iris lactea)、山羽藓(Abietinella abietina)和匐灯藓(Plagiomnium cuspidatum)等苔藓植物组成。

1.2 植物材料

以祁连山北坡大野口流域分布的11种木本植物作为研究材料。包括2种乔木(青海云杉和祁连圆柏)和9种灌木(鬼箭锦鸡儿、金露梅、银露梅、高山绣线菊、吉拉柳、鲜黄小檗、水栒子、刚毛忍冬和猫儿刺)。

1.3 试验方法

共选取7个采样点(20 m ×20 m), 样点海拔范围为2300—3650 m, 各采样点沿海拔间隔为100—400 m(图 1)。在2020—2021年生长季收集了各样点的环境因子并测定了11种主要木本植物的功能性状。在设定样点采集可获取的成年阳生叶与枝条分别测定6个枝、叶功能性状, 包括木质密度(WD)、胡伯尔值(HV)、叶面积(LA)、叶干物质含量(LDMC)、叶含水率(LWC)和比叶面积(SLA)。在每个样点内, 每个物种至少从3株成年个体采样, 且每个指标不少于6个重复。

图 1 祁连山北坡主要木本植物采样点分布图 Fig. 1 Distribution of sampling sites of main woody plants on the northern slope of Qilian Mountains
1.3.1 环境因子获取

环境因子的数据来源于中国生态系统研究网络临泽内流河流域研究站大野口观测点, 选取海拔(ALT)、生长季气温(Tgs)、生长季降水(Pgs)作为环境因子。使用同年5—10月各样点的日均温作为Tgs, 总降水量作为Pgs。其中, ALT使用商用GPS测定;Tgs和Pgs使用气温计和雨量计测定。其中, 气温计设在距地表 1 m处;雨量计设在林窗处, 消除了冠层截留的干扰。

1.3.2 枝功能性状测定

在每个样点, 每个物种选取3个成年植株, 在每个植株上随机选取2个枝条用来测定木质密度和胡伯尔值。木质密度测量方法为, 用枝剪截取中间5—8 cm茎段, 剥掉树皮, 放置在蒸馏水中浸泡12 h使其达到饱和状态。随后将茎断面用少量凡士林密封。利用排水法测定茎段体积。随后将茎段放置在80℃烘箱中烘72 h, 用精度为0.0001 g的电子天平(ME204E, 梅特勒托利多, 苏黎世, 瑞士)称重得到茎段干重。用得到的茎段干重除以茎段体积计算得到木质密度。

胡伯尔值由枝条边材面积与总叶面积的比值计算得到。将所测枝条过夜复水后(水势≥0.3 Mpa)将所有叶片取下, 用扫描仪(CanoScan LiDE300, 佳能, 东京, 日本)和Image J软件(http://imagej.nih.gov/ij/)扫描计算得到总叶面积。枝条边材面积使用染色法测定, 将稀释红色墨水从枝条的基部流入茎段, 随后从中间截断茎段, 用扫描仪(CanoScan LiDE300, 佳能, 东京, 日本)扫描断面, 染为红色的部分为边材, 中间未染色的部分为心材。

1.3.3 叶功能性状测定

使用同一批采集的样本同时测定与计算得到叶面积、叶含水率、叶干物质含量和比叶面积。首先, 随机选取6—8片成年叶(水势≥0.3 Mpa), 使用精度为万分之一克的分析天平称重, 得到叶鲜重值。然后用扫描仪(CanoScan LiDE300, 佳能, 东京, 日本)将叶面积扫描, 使用Image J软件(http://imagej.nih.gov/ij/)计算得到叶面积。将叶按编号分别装入信封, 置于80℃烘箱中烘72 h, 烘至质量不再减少时用电子天平(ME204E, 梅特勒托利多, 苏黎世, 瑞士)称量得到叶干重。用叶鲜重与叶干重之差除以叶干重计算得到叶含水率;用叶干重除以叶鲜重计算得到叶干物质含量;用叶面积除以叶干重计算得到比叶面积。

1.4 数据分析

利用单因素方差分析(ANOVA)检验各功能性状在不同生长型植物间的差异, 检验水平设定为α=0.05; 用线性回归分析了6个植物功能性状沿海拔的分异;用Pearson相关分析检验各功能性状之间以及各功能性状与环境因子的相关关系, 相关性的显著性使用Student′s t检验法进行检验;最后, 用主成分分析(PCA)分析各性状之间的关系以及不同生长型植物对不同主成分的贡献。所有的统计分析与绘图均在R 4.0.2软件[28]中完成。

2 结果与分析 2.1 11种木本植物的功能性状特征

通过对11种主要木本植物功能性状数据描述性统计分析发现(表 1), LA偏度与峰度均最大, 分别为4.15和18.63。6个功能性状的变异系数范围在16.81%—149.21%。其中, WD变异系数最小, LA变异系数最大(表 1)。比较乔、灌各功能性状均值发现, 两种生长型木本植物的各功能性状之间存在显著差异。乔木的LA、LDMC显著高于灌木, WD、HV、LWC和SLA显著低于灌木(图 2)。

表 1 祁连山北坡主要木本植物功能性状的描述性统计 Table 1 Descriptive statistics of functional traits of main woody plants on the northern slope of Qilian Mountains
功能性状
Functional trait
样本量
Sample size
平均值±
标准偏差
Mean±SD
最小值
Minimum
最大值
Maximum
中值
Median
方差
Variance
偏度
Skewness
峰度
Kurtosis
变异系数
Coefficient of variation/%
木质密度
Woody density/(g/cm3)
154 0.55±0.12 0.74 0.24 0.6 0.02 -0.66 -0.74 0.23
胡伯尔系数
Huber value/(cm2/cm2)
171 4.6±3.29 0.48 17.93 3.98 10.79 1.31 1.95 71.35
叶含水率
Leaf water content/(g/g)
154 1.85±0.54 0.93 3.54 1.78 0.29 1.1 1.04 29.05
叶面积  Leaf area/cm2) 154 8.68±12.95 1.42 88.3 4.72 167.81 4.15 18.63 149.21
比叶面积
Specific leaf area/(cm2/g)
154 158.10±28.1 28.10 335.57 156.60 5315.85 0.11 -0.58 46.12
叶干物质含量
Leaf dry mass content/(g/g)
154 0.36±0.06 0.22 0.52 0.36 0 -0.23 -0.42 16.81

图 2 祁连山北坡主要木本植物功能性状在不同生长型物种之间差异 Fig. 2 Differences in plant functional traits of main woody plants on the northern slope of Qilian Mountains among different life forms *** P < 0.001, ** P < 0.01, * P < 0.05
2.2 植物功能性状沿海拔分异规律

各植物功能性状沿海拔分异规律存在差异(图 3)。线性回归分析表明, WD与LDMC沿海拔呈显著下降趋势, 其他功能性状沿海拔分异趋势不显著。各功能性状沿海拔分异存在种间差异, 种间比较发现, 尽管其他物种的WD与LDMC沿海拔均呈显著下降趋势, 但青海云杉的WD与LDMC沿海拔变化则不明显。鬼箭锦鸡儿、吉拉柳和银露梅的LWC沿海拔呈明显上升趋势, 其他物种变化不明显。不同物种的HV、LA和SLA沿海拔分异趋势不一致。

图 3 祁连祁连山北坡主要木本植物功能性状沿海拔分异特征 Fig. 3 Differential characteristics along altitude of main woody plant functional traits on the northern slope of Qilian Mountains 黑色实线表示线性拟合显著(P < 0.05), 黑色虚线表示线性拟合不显著(P>0.05)
2.3 植物功能性状间及与环境因子的关系

Pearson相关性检验结果表明各功能性状与环境因子存在不同程度的相关关系(图 4)。其中, ALT与LDMC显著负相关, 与HV和LWC显著正相关;Tgs与HV和LWC负相关, 与LDMC正相关;Pgs与各性状均不相关。性状间相关性分析表明, 枝功能性状HV与LDMC负相关, 与LWC正相关;叶功能性状SLA、LWC和LDMC之间均存在显著相关关系;LA与SLA显著正相关, 与其他功能性状不相关。

图 4 祁连山北坡主要木本植物功能性状与环境因子的Pearson相关性系数 Fig. 4 Pearson coefficients of main woody plant functional traits and environmental factors on the northern slope of Qilian Mountains WD:木质密度, Woody density;HV:胡伯尔值, huber value;LA:叶面积, leaf area;LDMC:叶干物质含量, leaf dry mass content;LWC:叶含水率, leaf water content;SLA:比叶面积, specific leaf area;ALT:海拔, altitude;Tgs:生长季气温, growing-season temperature;Pgs:生长季降水, growing-season precipitation. *** P < 0.001, ** P < 0.01, * P < 0.05, P < 0.1

主成分分析揭示了不同功能性状和生长型木本植物对各主成分的贡献(图 5)。其中, 第一主成分解释了功能性状变异的52.5%, 第二主成分解释了功能性状变异的26.5%, 前两个主成分累积解释接近80%, 代表了功能性状变异的绝大部分。LWC、HV、LDMC和SLA是第一主成分的主要贡献变量, WD和LA是第二主成分的主要贡献变量。主成分分析可以将灌木与乔木明显区分开, 分别位于第一主轴的右侧和左侧。乔木主要贡献了LA与LDMC的变异, 灌木则主要贡献了SLA、LWC和HV的变异。

图 5 祁连山北坡木本植物功能性状的主成分分析(PCA) 结果 Fig. 5 Principal component analysis (PCA) results of woody plant functional traits on the northern slope of Qilian Mountains 红色三角和绿色圆圈分别代表乔木和灌木
3 讨论 3.1 祁连山北坡不同生长型木本植物功能性状特征

植物功能性状是植物在特定生境下所形成的适应性特征, 分布于不同生境的植物类群往往会形成特异的单个或多个性状特征来适应当地的环境选择压力[1011, 2930]。在山地生态系统中, 灌木相比于乔木可以在更大的海拔范围分布, 因此灌木需要进化出相应的高海拔性状适应特性[23, 31]。对比不同生长型木本植物枝性状发现, 相比于乔木, 灌木具有较大的木质密度和胡伯尔值, 这反映了灌木枝的结构性碳投资较大, 有利于其抵御高海拔低温生境, 并且, 由于具有较大的胡伯尔值, 灌木单位叶面积得到的水力支持更大, 使其在生长期获得充足的水分供给, 从而实现在高海拔短暂的生长季内快速生长发育[32]。对比不同生长型木本植物叶性状发现, 相比乔木, 灌木的叶面积、叶干物质含量较小, 叶含水率、比叶面积较大。灌木具有较低的叶干物质含量表明其用于叶片结构建设的碳投资较低, 这与其枝较大的结构性碳投资相反。这表明乔、灌木的枝、叶采取不同的资源利用策略来适应高山生境。相较之下, 灌木枝和乔木叶采取较为“保守型”的资源利用策略, 而灌木叶和乔木枝则采取“获取型”的资源利用策略。这可能是由于祁连山北坡分布的两种优势乔木青海云杉和祁连圆柏均为常绿植物, 因此需要将更多的碳分配到叶结构建设以抵抗冬季低温;而灌木均为落叶植物, 通过落叶来“躲避”低温胁迫, 因而无需将过多的碳分配到叶结构建设, 而是用于枝结构建设。

3.2 植物功能性状沿海拔分异规律

探究植物功能性状沿海拔分异有利于理解植被如何适应高山生境, 从而为高山植被边界形成提供理论依据[7, 33]。分析祁连山北坡木本植物功能性状沿海拔分异发现, 木质密度和叶干物质含量沿海拔升高显著下降, 其他功能性状沿海拔分异不显著, 表明不同海拔分布的木本植物形成了不同的功能性状特性以适应不同海拔生境。木质密度和叶干物质含量沿海拔下降反映了沿海拔升高植物结构性碳投资减少[3435], 表明低海拔分布的木本植物形成了较高的结构性碳投资的功能性状特性, 随着海拔升高, 木本植物的功能性状的结构性碳投资逐渐降低。较低结构性碳投资使高海拔分布的木本植物在相对较短的生长季内将更多的资源用于生长发育。然而, 较低的结构性碳投资使高海拔分布的木本植物很难抵御低温胁迫, 因此可能限制了木本植物向更高的海拔“攀爬”[3638]。进一步分析表明, 植物功能性状沿海拔分异在乔、灌间也存在明显差异。本研究发现灌木的木质密度和叶干物质含量沿海拔降低非常明显, 而优势乔木种青海云杉的木质密度和叶干物质含量变化则不明显。类似的规律也在全球其他地区观察到[39]。这表明灌木性状对海拔变化更加敏感。因此, 在全球变暖情景下, 高山灌木的可能表现出更强的环境适应性, 反之, 乔木将面临环境变暖所带来的存活风险[4041]。植物功能性状沿海拔分异的种间差异同样体现在叶面积、胡伯尔值和比叶面积上, 这表明除环境因子外, 植物自身的遗传特性也会影响植物功能性状沿海拔分异, 从而使不同物种形成了特异的高海拔性状适应机制。

3.3 植物功能性状与环境因子的关系

在山地生态系统中, 海拔是引起水热条件变化主要原因[4243]。本研究发现海拔与生长季气温和降水均存在显著相关关系, 其中, 生长季气温与海拔显著负相关且相关性性系数接近-1(相关性系数为-0.96), 这表明沿海拔上升, 生长季气温呈显著下降趋势。相反, 生长季降水与海拔呈正相关关系, 但相关性系数远低于生长季气温(相关性系数为0.55), 表明生长季降水沿海拔虽然有增加的趋势, 但其对海拔依赖性远低于生长季气温。这可能与生长季降水沿海拔呈“驼峰型”变化趋势有关[1]。各环境因素与植物功能性状相关性分析发现, 海拔与胡伯尔值、叶干物质含量和叶含水率显著相关;生长季气温同样与胡伯尔值、叶干物质含量和叶含水率相关;而生长季降水与各功能性状均不相关。这表明在山地生态系统, 生长季气温对植物功能性状沿海拔分异的驱动作用更大, 而生长季降水作用较小[9, 15, 44]。此外, 本研究发现尽管生长季气温与各性状具有相关性, 但相关性程度远低于海拔(图 4)。这表明仅考虑生长季气温并不能完全解释性状沿海拔分异, 还应考虑其他环境因子以及物种自身的遗传特性对植物功能性状沿海拔分异的多重影响。

4 结论

祁连山北坡乔、灌木的枝、叶性状差异显著。灌木枝和乔木叶采取“保守型”的资源利用策略, 而灌木叶和乔木枝则采取“获取型”的资源利用策略。木本植物在不同海拔形成了不同的功能性状特性以适应生境。木质密度和叶干物质含量沿海拔下降明显, 表明相较于低海拔分布的木本植物, 高海拔分布的木本植物降低了功能性状的结构性碳投资。这种资源利用策略有助于高海拔植物在较短的生长季快速生长发育, 但低结构性碳投资同时也限制了木本植物在更高的海拔分布。相较于生长季降水, 生长季气温与植物功能性状沿海拔分异相关性更强, 但生长季气温不是唯一环境因子。未来研究应系统地从不同角度揭示植物功能性状沿海拔分异机制以及其对植被分布的直接或间接影响, 这将对理解气候变暖情景下高山植被分布格局动态变化具有重要意义。

参考文献
[1]
Körner C. The use of 'altitude' in ecological research. Trends in Ecology & Evolution, 2007, 22(11): 569-574.
[2]
Egorov P I, Agadzhan N, Korotaev M M, Baranova V P, Benevole T, Boikova O I, Vasilev V K, Zavadovs A, Kuzmin M P, Reutova M B, Tsyganov N, Yakovlev I. Effect of high-altitude acclimatization on resistance to environmental factors. Environmental Space Sciences-Ussr, 1969, 3(1): 69.
[3]
Li K, Peng W, Zhou Y L, Ren Y, Zhao J H, Fu X S, Nie Y Q. Host genetic and environmental factors shape the composition and function of gut microbiota in populations living at high altitude. BioMed Research International, 2020, 2020: 1482109.
[4]
Wakui A, Sueyoshi M, Shimokawabe A, Kudo G, Morimoto J, Nakamura F. Environmental factors determining the distribution of highland plants at low-altitude algific talus sites. Ecological Research, 2017, 32(2): 183-191. DOI:10.1007/s11284-016-1429-9
[5]
Ackerly D D, Reich P B. Convergence and correlations among leaf size and function in seed plants: a comparative test using independent contrasts. American Journal of Botany, 1999, 86(9): 1272-1281. DOI:10.2307/2656775
[6]
Swenson N G, Enquist B J. Ecological and evolutionary determinants of a key plant functional trait: wood density and its community-wide variation across latitude and elevation. American Journal of Botany, 2007, 94(3): 451-459. DOI:10.3732/ajb.94.3.451
[7]
Nicotra A B, Atkin O K, Bonser S P, Davidson A M, Finnegan E J, Mathesius U, Poot P, Purugganan M D, Richards C L, Valladares F, van Kleunen M. Plant phenotypic plasticity in a changing climate. Trends in Plant Science, 2010, 15(12): 684-692. DOI:10.1016/j.tplants.2010.09.008
[8]
Körner C, Neumayer M, Menendez-Riedl S P, Smeets-Scheel A. Functional morphology of mountain plants). Flora, 1989, 182(5/6): 353-383.
[9]
Körner C, Paulsen J. A world-wide study of high altitude treeline temperatures. Journal of Biogeography, 2004, 31(5): 713-732. DOI:10.1111/j.1365-2699.2003.01043.x
[10]
Read Q D, Moorhead L C, Swenson N G, Bailey J K, Sanders N J. Convergent effects of elevation on functional leaf traits within and among species. Functional Ecology, 2014, 28(1): 37-45. DOI:10.1111/1365-2435.12162
[11]
Wright I J, Reich P B, Westoby M, Ackerly D D, Baruch Z, Bongers F, Cavender-Bares J, Chapin T, Cornelissen J H C, Diemer M, Flexas J, Garnier E, Groom P K, Gulias J, Hikosaka K, Lamont B B, Lee T, Lee W, Lusk C, Midgley J J, Navas M L, Niinemets V, Oleksyn J, Osada N, Poorter H, Poot P, Prior L, Pyankov V I, Roumet C, Thomas S C, Tjoelker M G, Veneklaas E J, Villar R. The worldwide leaf economics spectrum. Nature, 2004, 428(6985): 821-827. DOI:10.1038/nature02403
[12]
Yang D, Wang A Y, Zhang J L, Bradshaw C J A, Hao G Y. Variation in stem xylem traits is related to differentiation of upper limits of tree species along an elevational gradient. Forests, 2020, 11(3): 349. DOI:10.3390/f11030349
[13]
Tyree M T, Davis S D, Cochard H. Biophysical Perspectives of Xylem Evolution: is there a Tradeoff of Hydraulic Efficiency for Vulnerability to Dysfunction?. IAWA Journal, 1994, 15(4): 335-360. DOI:10.1163/22941932-90001369
[14]
李志民, 王传宽. 木本植物木质部的冻融栓塞应对研究进展. 植物生态学报, 2019, 43(8): 635-647.
[15]
Mellert K H, Fensterer V, Küchenhoff H, Reger B, Kölling C, Klemmt H J, Ewald J. Hypothesis-driven species distribution models for tree species in the Bavarian Alps. Journal of Vegetation Science, 2011, 22(4): 635-646. DOI:10.1111/j.1654-1103.2011.01274.x
[16]
Payette S. Contrasted dynamics of northern Labrador tree lines caused by climate change and migrational lag. Ecology, 2007, 88(3): 770-780. DOI:10.1890/06-0265
[17]
Lu X M, Camarero J J, Liang E Y. Threshold responses of juniper tree growth and regeneration to climate warming and drought stress at alpine treeline. Trees, 2021, 35(4): 1081-1083. DOI:10.1007/s00468-021-02135-6
[18]
马蓉蓉, 黄雨晗, 周伟, 周际, 白中科, 官炎俊, 郑连福, 詹培元, 杨正, 张艳. 祁连山山水林田湖草生态保护与修复的探索与实践. 生态学报, 2019, 39(23): 8990-8997.
[19]
王国宏. 祁连山北坡中段植物群落多样性的垂直分布格局. 生物多样性, 2002, 10(1): 7-14. DOI:10.3321/j.issn:1005-0094.2002.01.002
[20]
Grytnes J A, Kapfer J, Jurasinski G, Birks H H, Henriksen H, Klanderud K, Odland A, Ohlson M, Wipf S, Birks H J B. Identifying the driving factors behind observed elevational range shifts on E uropean mountains. Global Ecology and Biogeography, 2014, 23(8): 876-884. DOI:10.1111/geb.12170
[21]
Kearney M, Porter W. Mechanistic niche modelling: combining physiological and spatial data to predict species' ranges. Ecology Letters, 2009, 12(4): 334-350. DOI:10.1111/j.1461-0248.2008.01277.x
[22]
李坤鹏, 胡辉, 罗超, 强平, 邓媛. 祁连山乌鞘岭地区植被组成与分布特征研究. 环境生态学, 2020, 2(1): 79-84, 88.
[23]
Liang E Y, Wang Y F, Piao S L, Lu X M, Camarero J J, Zhu H F, Zhu L P, Ellison A M, Ciais P, Peñuelas J. Species interactions slow warming-induced upward shifts of treelines on the Tibetan Plateau. PNAS, 2016, 113(16): 4380-4385. DOI:10.1073/pnas.1520582113
[24]
Liu J G, Gou X H, Gunina A, Long X E, Zhang F, Zhang J Z. Soil nitrogen pool drives plant tissue traits in alpine treeline ecotones. Forest Ecology and Management, 2020, 477: 118490. DOI:10.1016/j.foreco.2020.118490
[25]
Feng X Y, Lin P F, Zhao W Z. The physiological constraints of alpine treeline in Qilian Mountains. Forest Ecology and Management, 2022, 503: 119761. DOI:10.1016/j.foreco.2021.119761
[26]
孙雪娇, 常顺利, 张毓涛, 李吉玫. 天山森林植物功能性状与碳库沿海拔梯度的变化. 生态学报, 2018, 38(14): 4994-5005.
[27]
向响, 黄永梅, 杨崇曜, 李泽卿, 陈慧颖, 潘莹萍, 霍佳璇, 任梁. 海拔对青海湖流域群落水平植物功能性状的影响. 植物生态学报, 2021, 45(5): 456-466.
[28]
R: A language and environment for statistical computing. Vienna, Austria, R Foundation for Statistical Computing, 2020.
[29]
Wright I J, Dong N, Maire V, Prentice I C, Westoby M, Díaz S, Gallagher R V, Jacobs B F, Kooyman R, Law E A, Leishman M R, Niinemets V, Reich P B, Sack L, Villar R, Wang H, Wilf P. Global climatic drivers of leaf size. Science, 2017, 357(6354): 917-921. DOI:10.1126/science.aal4760
[30]
Osnas J L D, Lichstein J W, Reich P B, Pacala S W. Global leaf trait relationships: mass, area, and the leaf economics spectrum. Science, 2013, 340(6133): 741-744. DOI:10.1126/science.1231574
[31]
Guo Z W, Lin H, Chen S L, Yang Q P. Altitudinal patterns of leaf traits and leaf allometry in bamboo Pleioblastus amarus. Frontiers in Plant Science, 2018, 9: 1110. DOI:10.3389/fpls.2018.01110
[32]
Mencuccini M, Rosas T, Rowland L, Choat B, Cornelissen H, Jansen S, Kramer K, Lapenis A, Manzoni S, Niinemets V, Reich P, Schrodt F, Soudzilovskaia N, Wright I J, Martínez-Vilalta J. Leaf economics and plant hydraulics drive leaf: wood area ratios. The New Phytologist, 2019, 224(4): 1544-1556. DOI:10.1111/nph.15998
[33]
Anderson J T, Gezon Z J. Plasticity in functional traits in the context of climate change: a case study of the subalpine forbBoechera stricta(Brassicaceae). Global Change Biology, 2015, 21(4): 1689-1703. DOI:10.1111/gcb.12770
[34]
Momo S T, Ploton P, Martin-Ducup O, Lehnebach R, Fortunel C, Sagang L B T, Boyemba F, Couteron P, Fayolle A, Libalah M, Loumeto J, Medjibe V, Ngomanda A, Obiang D, Pélissier R, Rossi V, Yongo O, Sonké B, Barbier N. Leveraging signatures of plant functional strategies in wood density profiles of African trees to correct mass estimations from terrestrial laser data. Scientific Reports, 2020, 10: 2001. DOI:10.1038/s41598-020-58733-w
[35]
Preston K A, Cornwell W K, Denoyer J L. Wood density and vessel traits as distinct correlates of ecological strategy in 51 California coast range angiosperms. The New Phytologist, 2006, 170(4): 807-818. DOI:10.1111/j.1469-8137.2006.01712.x
[36]
Kiaei M, Moosavi V, Ebadi S E. Effects of altitude on density and biometric properties of hornbeam wood (Carpinus betulus). Forest Systems, 2019, 28(2): e011.
[37]
Giagli K, Vavrcik H, Fajstavr M, Cerny J, Novosadova K, Martinik A. Stand factors affecting the wood density of naturally regenerated young silver birch growing at the lower altitude of the Czech Republic region. Wood Research, 2019, 64(6): 1011-1022.
[38]
Ramananantoandro T, Ramanakoto M F, Rajoelison G L, Randriamboavonjy J C, Rafidimanantsoa H P. Influence of tree species, tree diameter and soil types on wood density and its radial variation in a mid-altitude rainforest in Madagascar. Annals of Forest Science, 2016, 73(4): 1113-1124.
[39]
Fajardo A. Insights into intraspecific wood density variation and its relationship to growth, height and elevation in a treeline species. Plant Biology: Stuttgart, Germany, 2018, 20(3): 456-464.
[40]
Pandey S, Cherubini P, Saurer M, Carrer M, Petit G. Effects of climate change on treeline trees in Sagarmatha (Mt. Everest, Central Himalaya). Journal of Vegetation Science, 2020, 31(6): 1144-1153.
[41]
McDowell N G, Allen C D. Darcys law predicts widespread forest mortality under climate warming. Nature Climate Change, 2015, 5(7): 669-672.
[42]
Sartor J D, Cannon T W. Effect of ice nucleation altitude on production of precipitation in a convective circulation. Transactions-American Geophysical Union, 1975, 56(12): 993-993.
[43]
Quiroz R S, Miller A J. Height-lag correlations of density with pressure and temperature at rocket altitudes of the stratosphere. Journal of the Atmospheric Sciences, 1968, 25(1): 104-112.
[44]
Normand S, Treier U A, Randin C, Vittoz P, Guisan A, Svenning J C. Importance of abiotic stress as a range-limit determinant for European plants: insights from species responses to climatic gradients. Global Ecology and Biogeography, 2009, 18(4): 437-449.