中国针叶林优势树种叶片氮磷钾生态化学计量特征及内稳态分析

邓博文，许瑶瑶，陈逸飞，张 硕，郑德祥＊，廖晓丽，林仁忠，靳少非,4

(1. 福建农林大学林学院，福建 福州 350002；2. 闽江学院地理科学系，福建 福州 350108；3. 顺昌县国有林场，福建 顺昌 353200；4. 华东师范大学崇明生态研究院，上海 202162)

生态化学计量学是研究生态系统中养分平衡的科学，为研究植物-环境相互作用提供了良好的指标和研究方向[1-3]。探究植物叶片N、P、K与环境因子的相关性，有利于挖掘植物关键生物元素循环规律，为研究植物的元素平衡吸收利用提供理论依据[4-5]。在气候变化、植物生理生化作用、土壤地球化学分布特征、植物群落组成等多重因素影响下[6-8]，森林生态化学计量存在着纬向变化特征，如对全球452个点、1 280种植物、5 087个观测值的研究发现，随着纬度的增加，平均温度降低，植物叶片N、P含量增加，N/P比下降[7]。对中国区域及中国东部南北样带植物N、P含量的研究结果也显示，叶片N、P含量呈现随纬度降低而减少的趋势，但N/P无显著变化[9-11]。目前，对于另一关键化学计量元素K与其他营养元素的生态化学计量分析未得到广泛研究[12]。

针叶林在我国各气候带内均有分布，共6科30属约200种，其广阔的跨气候带的分布特征为研究针叶林生态化学计量空间大尺度变化特征提供了研究素材[13]，对探究针叶树种叶片化学计量特征纬向变化具有重要意义。植物叶片生态化学计量的纬向变化特征是多环境影响因素的整体效应，温度、水分、土壤养分等环境因子也对植物生态化学计量特征产生一定的影响，N、P、K是影响植物生长发育的重要指标，反映出土壤生产力的强弱[14-15]。随着温度升高，植物叶片中N/P呈增加趋势[7,11]，闫晓俊等[16]在对杉木幼苗增温实验中发现，大气和土壤单独增温均可增加叶片N含量，大气和土壤同时增温显著增加了叶片P含量，但对N/P影响不显著；在干旱区域，叶片的化学计量与降水呈显著相关，而在热带多雨地区，并未发现显著关系[17]。自20世纪以来，人工合成氮肥技术的发展导致氮肥广泛应用、化石燃料大规模开采利用，增加了大气N沉降量[18]。刘文飞等[19]在对杉木进行大气氮沉降模拟实验发现，增加N输入可提高植物叶片N含量，降低植物P、K含量，进而增加了N/P、N/K；过量的N输入改变了植物生态化学计量特征，表现为叶片N含量增加，增加N/P比，导致植物生长更易受到P限制作用[20-22]。土壤生态化学计量特征可以直接影响叶片生态化学计量特征，并通过内稳态调节机制应对土壤养分的变化[23-25]。目前，对于探究多环境因子对叶片化学计量变化的研究较少，因此，探究不同环境因子对植物叶片化学计量特征影响的贡献及植物内稳态状态对于综合了解叶片生态化学计量特征具有重要的意义。

本研究选取中国区域针叶林优势树种为对象，分析中国区域优势针叶林土壤－植物N、P、K含量以及化学计量特征，探究不同环境因子对叶片化学计量特征的影响，分析中国针叶树种叶片化学计量稳态性特征，该研究结果可为综合研究针叶林生物地球化学循环特征与全球环境变化响应提供数据基础和理论支持。

1 材料和方法

1.1 数据来源与处理

本研究中所使用的数据均来自于已有文献数据[26-33]，经过匹配，共81个取样点的土壤养分数据被纳入，计262个土壤剖面，每个剖面设置3个土壤深度：0～20、20～40、40～60 cm。

本研究中收集了包括西伯利亚落叶松(Larix sibiricaLedeb.)、黄花落叶松(L. olgensisHenry.)、鱼鳞云杉(Picea jezoensisCheng.)、臭冷杉(AbiesnephrolepisMaxim.)、石松(Lycopodium clavammThunb.)、扁枝石松(L. complanatumHolub.)、玉柏石松(L. obscurumLinn.)、杉蔓石松(L. annotinumL.)、赤松(Pinus densifloraSieb.)、西伯利亚松(P.sylvestrisLitv.)、青海云杉(P. crassifoliaKom.)、油松(P. tabuliformisCarriere.)、白皮松(P. bungeanaZucc.)、马尾松(P. massonianaLamb.)、云南松(P.yunnanensisFranch.)、 华 山 松 (P. armandiiFranch.)、杉木(Cunninghamia lanceolata(Lamb.)Hook.)、侧柏(Platycladus orientalisFranco.)、柏木(Cupressus funebrisEndl.)、干香柏(C. duclouxianaHichel.)20种具代表性的针叶林树种。

1.2 数据处理

1.2.1 内稳态指数 内稳态指数（H）采用下式计算：

式中：x为土壤生态化学计量特征(如N、N/P等)，y为叶片生态化学计量，c为拟合常数。1/H表示log10(x)和log10(y)回归斜率，其绝对值范围为0～1，生物体内稳态越稳定越高，其H值越大。

对内稳态按照稳定性进行分类级别划分[34]，回归分析显著性检验采用α= 0.1进行单侧检验，若回归关系不显著(P＞ 0.1)，被定义为“绝对稳态”；若回归关系显著，当1/H= 1，则植物体不具有内稳态性质；0 ＜ 1/H＜ 1的数据集分别定义为：0 ＜1/H＜ 0.25为“稳态”，0.25 ＜ 1/H＜ 0.5为“弱稳态”，0.5 ＜ 1/H＜ 0.75为“弱敏感态”，1/H＞ 0.75为“敏感态”。

1.2.2 环境因子相关性 为探究不同环境因子对叶片生态化学计量的影响，本研究利用结构方程模型分析纬度、经度、年降水量、年均气温、最暖月均气温、最冷月均气温、年均气温差、水热指数、年蒸发量、水分亏欠指数、土壤不同层次生态化学计量等参数对中国优势针叶树种叶片化学计量的影响。在分析过程中，设置4个潜变量，分别为温度相关变量、降水相关变量、地理区域变量、土壤养分供应；1个显变量为水分亏欠指数，分析不同变量对叶片化学计量影响。具体的模型见图1。

在结构方程模型分析的基础上，利用R软件relaimpo包计算各个变量对叶片化学计量的相对贡献；采用One-way ANOVA分析不同土壤层次中N、P、K、N/P、N/K、K/P等化学计量学特征；利用一元线性回归分析针叶林叶片中N、P、K、N/P、N/K、K/P与土层次中的相关关系，并分析其稳态性。研究中所采用的统计软件为R (R version 3.5.3)。

图 1 结构方程模型结构Fig. 1 The frame of structure equation model

2 结果与分析

2.1 中国针叶林优势树种叶片及土壤生态化学计量分析

中国区域针叶林叶片N、P、K平均含量分别为13.16、1.02、6.82 g·kg-1，N/P、N/K、K/P平均比分别为15.48、2.17、7.77。土壤中N、N/P、N/K在不同土层内存在差异（图2），其中，土层0～20 cm N、P、K含量及N/P、N/K最高；0～20、20～40、40～60 cm土层的N平均含量分别为1.73、0.99、0.72 g·kg-1；0～20、20～40、40～60 cm土层的N/P分别为3.10、2.12、1.74，N/K分别为0.10、0.06、0.05。

2.2 不同环境要素对中国针叶林优势树种叶片生态化学计量贡献分析

结构方程模型表明：多种环境因子与我国优势针叶林树种叶片化学计量特征存在显著关系（图1）。由表1可知：纬度对叶片N含量贡献率最高，为19.18%；水热指数对叶片P含量贡献率最高，为12.75%；纬度对叶片K含量贡献率最高，为25.08%；对N/P、N/K、P/K生态化学计量贡献率最高的分别为年降雨量（贡献率15.65%）、土壤40～60 cm深度N/K值（贡献率16.76%）及年降雨量（贡献率18.21%）。此外，土壤因子中，N、P、K、N/P、N/K、K/P对叶片化学计量贡献率分别为17.01%、11.20%、8.35%、16.27%、30.79%、27.52%。

2.3 中国针叶林优势树种叶片化学计量稳态性特征

对我国优势针叶树叶片中化学计量学特征随土壤环境的内稳态特征进行分析，结果（表2）表明：叶片P、K、N/P、N/K、P/K利用稳态性模型模拟结果不显著，属绝对稳态型；而叶片N内稳态指数为0.261，属于弱稳态型。

3 讨论

3.1 中国针叶林优势树种叶片生态化学计量分析

本研究利用文献资料探究了中国针叶林叶片N、P、K含量及其生态化学计量学特征，结果表明，中国针叶林优势树种叶片N ＞ K ＞ P，这与之前对兴安落叶松（Larix gmelinii Rupr.）及东部针叶树、黄土高原地区叶片养分含量的研究结果相似（表3）[13，35-36]。由于数据来源具有差异性，本研究系统对比了其他学者研究结果发现，叶片K含量及与之相对应的研究过少，因此，对N、P与之相对应的化学计量重点研究。

植物叶片中N和P含量是研究关键生境要素相对限制性的重要指标。依据已有的文献指标，当N/P ＜ 14时，植物生长受到N限制；当N/P ＞16时，则主要受到P限制；当N/P介于14～16时，同时受到N、P限制[37]，这也是我国总体样带植物的叶片化学计量特征[9]。植物营养元素含量的多寡受到土壤中营养条件的影响，并表现为对土壤营养元素环境的适应性[1-2,38]。

图 2 中国优势针叶树种不同土壤发生层中N、K、P含量以及化学计量分析Fig. 2 Box plots for the contents of N, P, and K, and their stoichiometry in different soil layers in China’s coniferous forests

表 1 不同环境因子对我国优势针叶树种叶片化学计量特征贡献率Table 1 Contribution of different environmental factors to leaf stoichiometry characteristics of dominant coniferous species in China %

表 2 中国优势针叶林叶片生态化学计量内稳态指数分析Table 2 Homeostasis analysis for the stoichiometry in China’s coniferous forests

表 3 中国针叶林N、P含量及N/P研究对比Table 3 Comparison of the N, P, and N/P among different studies

为研究中国区域针叶林土壤养分对植物养分的影响，选取3个土层的N、P、K数据用于研究植物养分供应。除K、P外，土壤N含量以及对应化学计量均随土壤深度增加而降低，其中，N/P、N/K随土壤深度的增加显著降低，而P、K含量呈相对稳定的分布特征。由此推断，随着中国N沉降的加剧，且伴随着降雨导致土壤淋溶增加，中国针叶林土壤N/P、N/K有增大趋势，进而使中国针叶林土壤P限制潜势继续增加。

3.2 不同环境要素对中国优势针叶林优势树种叶片生态化学计量贡献分析

在本研究中，中国针叶林优势树种叶片生态化学计量特征与环境因子之间具有显著的关系。不同于前人的研究以纬度、温度及降水之间的关系为主[7, 9, 11]，本研究中增加了更多的环境因子，如土壤养分化学计量特征、年平均温度、年降雨量、水热指数等。本研究发现，随着纬度增加，植物叶片中P含量呈显著增加，这表明我国针叶树种叶片P含量在空间上呈显著的纬向分布特征；此外，水热指数和土壤能够解释12.75%和11.20%的P贡献率，表明叶片中土壤P含量以及降水、温度三者对叶片P含量有显著影响。任书杰等[35]研究东北高纬度针叶林N、P含量随纬度变化规律时发现，在一个相对小的高纬度地区，随着纬度增加，叶片N含量会随着纬度增加反而呈显著减少的趋势，而P含量未达到统计学显著变化，这与中国区域其他功能群植物叶片变化特征不同[9-11]，本研究未发现叶片中N含量随着纬度显著增加，推测产生这种现象的原因，一方面是由于本研究所选取的天然林中的针叶树相对于阔叶树生长缓慢，且土壤因子对叶片N的贡献合计达17.01%，表明中国优势针叶林叶片N含量主要受到了土壤N含量的影响；另一方面是由于针对兴安落叶松的研究尺度较小，其他影响因素对叶片化学计量的影响超过了纬度的贡献。本研究发现，针叶林叶片中K含量随着纬度增加呈显著增加，推测这与中国区域土壤K含量分布的规律一致，推测原因为在南方高温多雨的地区，风化造成K容易流失[39]，而北方温度低且降雨少，土壤中K易于积累，不易成为限制元素。因此，中国南北气候带差异也造成了植物叶片内K含量南低北高的空间分布规律。

3.3 中国针叶林优势树种叶片化学计量稳态性特征

植物养分含量和化学计量特征与土壤养分供应之间存在着显著关系[1,38]，在土壤养分环境发生变化时，植物可通过内稳态调节机制来保持体内化学计量特征稳定性。本研究通过分析中国针叶林优势树种内稳态调节机制发现，我国针叶林优势树种除叶片N含量易受到土壤N含量变化影响外，其他的叶片生态化学计量特征均呈绝对稳态，这表明我国优势针叶树种具有良好的土壤环境适应特征，但目前人类工业活动释放的过量N以及全球温度增加对土壤中氮素含量变动产生了一定影响[16]。已有研究证明[19]，在模拟大气氮沉降增加情景下，针叶叶片中N含量增加，能够显著影响与之相关的其他化学计量特征，因此，需要多加关注土壤N含量变动对叶片生态化学计量特征产生的影响。

4 结论

中国区域针叶林叶片N、P、K平均含量分别为13.16、1.02、6.82 g·kg-1，N/P、N/K、K/P平均比分别为15.48、2.17、7.77。纬度、土壤、降水和温度因子共同对我国针叶林优势树种叶片化学计量特征产生影响，纬度对叶片N和K贡献率最高，水热指数对叶片P含量贡献率最高。对N/P、N/K、P/K影响最大的分别为年降雨量、土壤40～60 cm深度的N/K值以及年降雨量。除N含量外，中国针叶林优势树种叶片生态化学计量特征均呈绝对稳态，这表明我国优势针叶树种具有良好的环境适应特征。