植物生态学报  2016 , 40 (7): 686-701 https://doi.org/10.17521/cjpe.2015.1043

Orginal Article

氮添加对山西太岳山天然油松林主要植物叶片性状的影响

肖迪1, 王晓洁1, 张凯1, 何念鹏2, 侯继华1*

1北京林业大学资源与生态系统过程北京市重点实验室, 北京 100083
2中国科学院地理科学与资源研究所生态系统网络观测与模拟重点实验室, 北京 100101

Effects of nitrogen addition on leaf traits of common species in natural Pinus tabuliformis forests in Taiyue Mountain, Shanxi Province, China

XIAO Di1, WANG Xiao-Jie1, ZHANG Kai1, HE Nian-Peng2, HOU Ji-Hua1*

1Key Laboratory for Forest Resources & Ecosystem Processes of Beijing Forestry University, Beijing 100083, China
and 2Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China

通讯作者:  *通信作者Author for correspondence (E-mail: houjihua@bjfu.edu.cn)

责任编辑:  XIAO DiWANG Xiao-JieZHANG KaiHE Nian-PengHOU Ji-Hua

收稿日期: 2015-01-29

接受日期:  2015-12-27

网络出版日期:  2016-07-10

版权声明:  2016 植物生态学报编辑部 本文是遵循CCAL协议的开放存取期刊,引用请务必标明出处。

基金资助:  中央高校基本科研业务费专项资金(TD- 2011-07)、国家自然科学基金青年基金(31000263)、国家自然科学基金面上项目(31470506)

展开

摘要

为探讨植物性状对大气氮沉降的响应与适应机制, 该文以中国特有的、在北方温性针叶林中广泛分布的天然油松(Pinus tabuliformis)林为研究对象, 在2009-2013年开展了氮添加对植物叶片性状影响的野外控制试验, 4个氮添加浓度分别为0 kg·hm-2·a-1 (CK)、50 kg·hm-2·a-1 (低氮)、100 kg·hm-2·a-1 (中氮)和150 kg·hm-2·a-1 (高氮)。试验过程中分别测定了油松、蒙古栎(Quercus mongolica)、茶条槭(Acer ginnala)、毛榛(Corylus mandshurica)、沙梾(Cornus bretschneideri)、绣线菊(Spiraea salicifolia)、金银忍冬(Lonicera maackii)、羊须草(Carex callitrichos)、龙常草(Diarrhena mandshurica)、大火草(Anemone tomentosa)和玉竹(Polygonatum odoratum)等11种主要植物的9种叶片性状, 包括叶厚度(LT)、比叶面积(SLA)、干物质含量(LDMC)、叶氮含量(LNC)、叶磷含量(LPC)等。结果表明: 1)在氮添加影响下, 玉竹等个别物种的LTSLA、绣线菊等部分物种的叶面积(LA)和LDMC差异显著, 上述所有物种的LNC与大多数物种的叶绿素含量(CC)、LPC显著增加, 油松等9种植物叶片N:P发生显著变化, 不同年龄、不同类型的植物叶片对氮添加的响应不同。2)叶性状之间普遍存在显著相关性, 如SLALNCLPC极显著正相关, LTLNCLPC极显著负相关, 且相关性随氮添加强度变化。3) 11种植物的叶片特征空间分布规律与叶经济谱的描述一致, 氮添加使植物在特征空间中的位置向叶片薄、生长快、叶寿命短的“快速投资-收益型”一端发生移动; 在垂直方向上, 阔叶乔木、灌木及草本的位置与针叶乔木的移动方向相反。当环境改变时, 植物会改变生存策略, 调整资源分配, 从而保证物种间相对位置和群落整体结构的稳定性。叶经济谱的形成不依赖于环境的变化, 而是植物一种固有的属性。

关键词: 大气氮沉降 ; 叶经济谱 ; 叶片功能性状 ; 天然油松林 ; 氮添加

Abstract

Aims Our objectives were to explore the effects of simulated nitrogen deposition on leaf traits and ecological stoichiometry characteristics of common species in natural Pinus tabuliformis forests.
Methods We conducted the experiment of nitrogen (N) addition from 2009 to 2013 in the natural Pinus tabuliformis forests in Taiyue Mountain, Shanxi, China. The levels of N addition were 0 (control), 50 (low-N), 100 (medium-N) and 150 (high-N) kg·hm-2·a-1, respectively. Eleven common plant species in 12 20 m × 20 m plots were selected, including Pinus tabuliformis, Quercus mongolica, Acer ginnala, Corylus mandshurica, Cornus bretschneideri, Spiraea salicifolia, Lonicera maackii, Carex callitrichos, Diarrhena mandshurica, Anemone tomentosa, and Polygonatum odoratum. Nine leaf traits were measured, including leaf thickness (LT), specific leaf area (SLA), leaf dry matter content (LDMC), leaf nitrogen content (LNC), leaf phosphorus content (LPC), and other four.
Important findings We found that: 1) LT and SLA of Polygonatum odoratum significantly differed among four levels of N addition. Leaf area (LA) and LDMC of several species, such as Spiraea salicifolia, had significant difference among the N addition concentration. LNC of all species, chlorophyll content (CC) and LPC of most species increased significantly with the addition of N. Leaf N:P of 9 species varied significantly, and leaves with different types and ages showed different responses to N addition. 2) Leaf traits were significantly correlated with each other. For instance, SLA was significantly positively correlated with LNC and LPC. In contrast LT was negatively connected with LNC and LPC. In addition, the degree of correlation changed with the level of N addition. 3) The pattern of species distribution in leaf trait space was consistent with the prediction from the theory of Leaf Economic Spectrum (LES). N addition drove species moving along axis 1 in the trait space, and propelled them towards different directions along axis 2, which indicated that these species tended to take the “fast investment-return” strategy. These results suggested that with the change of environmental conditions, plants changed their survival strategy and adjusted resource allocation to maintain the stability of communities. This is the inherent characteristic of plants, thus the formation of LES did not depend on the environment change.

Keywords: atmospheric nitrogen deposition ; leaf economic spectrum ; leaf functional traits ; natural Pinus tabuliformis forest ; nitrogen addition

0

PDF (1672KB) 元数据 多维度评价 相关文章 收藏文章

本文引用格式 导出 EndNote Ris Bibtex

肖迪, 王晓洁, 张凯, 何念鹏, 侯继华. 氮添加对山西太岳山天然油松林主要植物叶片性状的影响. 植物生态学报, 2016, 40(7): 686-701 https://doi.org/10.17521/cjpe.2015.1043

XIAO Di, WANG Xiao-Jie, ZHANG Kai, HE Nian-Peng, HOU Ji-Hua. Effects of nitrogen addition on leaf traits of common species in natural Pinus tabuliformis forests in Taiyue Mountain, Shanxi Province, China. Chinese Journal of Plant Ecology, 2016, 40(7): 686-701 https://doi.org/10.17521/cjpe.2015.1043

自工业革命以来, 大气氮沉降呈现急剧增加的趋势(Holland et al., 1999; Zheng et al., 2002; Adams et al., 2004; Suding et al., 2005), 这直接影响到人类和生态系统的健康、温室气体的平衡及生物的多样性, 引起了科学家们的广泛关注与担忧(Vitousek et al., 1997; Richter et al., 2005; Clark & Tilman, 2008; Liu et al., 2013)。研究表明, 氮素作为陆地生态系统中最主要的限制元素(贺金生和韩兴国, 2010), 对森林植物有重要影响(李德军等, 2005)。在其他养分成为限制因子之前, 氮的增加可以提高植物的净初级生产力(Aber et al., 1989; Tamm, 1990), 但长期、连续、过量的氮输入会降低植物的净光合速率和净初级生产力, 引起森林生产力的下降和森林生态系统功能的退化(Ulrich, 1995; Fenn et al., 1996; Magill et al., 2004; 方运霆等, 2005), 这被认为是森林衰退的主要原因(Nihlgard, 1985; Bobbink et al., 1998; 张维娜和廖周瑜, 2009)。因此, 研究大气氮沉降对森林植物的影响具有十分重要的意义。

植物通常会对自身资源进行配置、补偿和平衡, 以最大限度地减小环境变化对植物的不利影响, 植物功能性状(functional traits)就是这种调节机制的具体体现(Comelissen et al., 2003), 因此成为探索植物与环境之间关系的重要纽带, 也为预测在全球变化背景下的生态系统响应提供了方法和依据(Wright et al., 2005)。其中, 叶片性状对气候变化十分敏感, 能较为准确地反映植物对气候变化的响应与适应机制(Royer et al., 2010; Scoffoni et al., 2011)。研究表明, 氮沉降可导致植物比叶面积(SLA)、叶片氮含量(LNC)、叶片磷含量(LPC)及氮磷比(N:P)等发生显著改变(Pontes et al., 2007; 黄菊莹等, 2009; Ordoñez et al., 2009; Rose et al., 2013; 王常顺和汪诗平, 2015)。

Wright等(2004)以叶片的形态、生理、化学性状对叶片的资源权衡策略进行量化, 发现植物叶片性状沿着一条连续变化的功能性状组合谱(即叶经济谱, LES)有序排列, 其一端代表着比叶重小、含氮量高、光合速率大、呼吸速率高、叶寿命短的“快速投资-收益型”策略; 而另一端代表着寿命长、比叶重大、含氮量、光合速率和呼吸速率都偏低的“缓慢投资-收益型”策略。叶经济谱为解释碳、氮元素在植物体内的固定和利用、营养物质的流动, 以及植被的分布如何随环境条件变化而改变等问题提供了理论依据(Royer, 2009; 陈莹婷和许振柱, 2014)。之后, 关于植物的资源权衡关系及全球尺度上植物叶片性状随环境的变化规律的研究大量涌现(Reich et al., 2003; Freschet et al., 2010; Laughlin, 2011; Osnas et al., 2013)。

然而, 关于植物叶经济谱的形成动力仍存在较大的争议: 一方面, Maire等(2013)和Blonder等(2013)认为环境梯度引发植物资源权衡发挥作用, 通过自我调节响应环境变化, 因此, 叶经济谱受环境因子主导; 另一方面, Wright等(2004)、Leishman等(2007)、Wright和Sutton-Grier (2012)等则认为植物的生活史策略和叶经济谱是独立演变的, 受气候变化和群落类型的影响微之甚微。基于全球尺度的“叶经济谱”假说有一个重要前提假设, 是物种性状的单一平均值具有足够的代表性。然而, 这一重要假设的有效性并没有得到充分验证。众所周知, 达尔文生物进化论最基本的一点就是同种个体在许多性状上的差异是生物适应环境而进化的基础。利用单一性状值来描述不同地区的同一个物种掩盖了该物种沿环境梯度某一或多个性状的巨大变异(Albert et al., 2010a)。这种忽视种内性状差异的观点将严重限制我们对生物与环境关系的理解, 极大地阻碍现有模型的预测能力(McGill et al., 2006; Wright & Sutton-Grier, 2012)。因此, 在控制条件下研究局域种库内物种的性状如何响应环境变化可以尽量避免其他生态因子的干扰, 能更有效地揭示植物对环境变化的适应机制, 而目前这类研究还较为少见(Albert et al., 2010b, 2011)。

油松(Pinus tabuliformis)是我国的特有种, 在北方地区分布广泛。山西太岳山地区素有“油松之乡”的美称, 有关该地区大气氮沉降对森林植物功能性状影响的研究鲜有报道。本研究试图以该地天然油松林为研究对象, 通过人工氮添加控制试验对自然条件下的氮沉降进行模拟, 观察、预测氮素变化对局域种库内植物叶片功能性状的影响, 并为解释植物应对环境变化的响应机制提供理论依据。本研究的主要目的是: (1)检验氮添加是否对局域种库内植物的叶片性状产生显著影响, 并分析性状变异的来源; (2)比较不同种类植物叶片对氮沉降响应的差异; (3)检验叶片性状间的相关性、叶经济谱存在与否及其形成动力; (4)探究天然油松林内主要植物叶片性状在经济谱中的位置是否受到氮添加的影响。

1 材料和方法

1.1 研究区概况

本研究区位于山西省长治市沁源县, 太岳山森林经营局灵空山自然保护区内(36.00°-37.08° N, 111.75°-112.55° E, 海拔1150-2088 m), 是我国暖温带落叶阔叶林典型代表地区。该研究区属暖温带半干旱大陆性季风气候, 年平均气温6.2 ℃, 四季分明, 春季多风。年降水量662 mm, 降水集中在7-9月, 占全年降水量的60%以上。相对湿度60%-65%, 年日照时间2600 h, 年无霜期110-125天。土壤以褐土和棕壤为主。油松、蒙古栎(Quercus mongolica)、山杨(Populus davidiana)、白桦(Betula platyphylla)是本研究区最具代表性的树种, 灌木主要有毛榛(Corylus mandshurica)、沙棘(Hippophae rhamnoides)、黄刺玫(Rosa xanthina)、胡枝子(Lespedeza bicolor)等, 草类有薹草(Carex spp.)、莎草(Cyperus spp.)、蒿类等。

1.2 试验设计

以欧洲氮饱和实验NITREX项目(Wright & Rasmussen, 1998)和北美哈佛森林(Harvard Forest)等(Aber et al., 1998)研究设计为参照, 在未被破坏的天然油松林(112.02°-112.25° E, 36.52°-36.72° N, 平均海拔1500 m)中建立12个20 m × 20 m的固定试验样地, 设置对照(CK, 0 kg·hm-2·a-1)、低氮(LN, 50 kg·hm-2·a-1)、中氮(MN, 100 kg·hm-2·a-1)和高氮(HN, 150 kg·hm-2·a-1) 4个处理, 3次重复。样地林分特征及0-20 cm土壤理化性质详见表1。从2009年5月起, 用尿素(CO(NH2)2)作为氮源, 将年施用量平均分成12等份, 每月下旬分别将所需的尿素溶解于2 L水中, 用背式喷雾器在对应处理的20 m × 20 m样地内均匀喷洒(对照组喷施等量的水), 其他水源来自自然降水。

表1   天然油松林各处理样地的林分特征和土壤理化性质

Table 1   Stand characteristics and soil physical-chemical properties of the nitrogen-loaded plots in the natural forest of Pinus tabuliformis

处理
Treatment
林分特征
Stand characteristics
土壤理化性质
Physical-chemical properties of soil
林龄
Stand age
(a)
密度
Density
(plant·hm-2)
平均胸径
Mean diameter at breast height (cm)
平均株高
Mean plant height
(m)
坡度
Slope
(°)
土壤容重
Soil volume weight (g·cm-3)
土壤pH Soil pH土壤有机碳含量
Soil organic carbon content (mg·g-1)
土壤氮含量
Soil nitrogen content
(mg·g-1)
土壤磷含量
Soil phosphorus content
(mg·g-1)
对照 Control751 26723.917.7240.947.1245.301.700.43
低氮 Low-N751 56720.617.8210.977.1261.022.140.38
中氮 Medium-N751 20823.517.4251.077.1942.141.910.38
高氮 High-N751 22523.419231.067.2836.082.250.44

Control, 0 kg·hm-2·a-1 nitrogen; Low-N, 50 kg·hm-2·a-1 nitrogen; Medium-N, 100 kg·hm-2·a-1 nitrogen; High-N, 150 kg·hm-2·a-1 nitrogen. The data was partly derived from Li et al. (2014).Control, 氮添加浓度为0 kg·hm-2·a-1; Low-N, 氮添加浓度为50 kg·hm-2·a-1; Medium-N, 氮添加浓度为100 kg·hm-2·a-1; High-N, 氮添加浓度为150 kg·hm-2·a-1。本表部分数据引自李化山等(2014)。

新窗口打开

1.3 研究对象

根据群落调查结果计算样地内各植物种的重要值, 选取出现频度最高、重要值大于0.05的11个物种进行研究, 包括3种乔木、4种灌木和4种草本, 分别是针叶树种油松(Pt)、阔叶树种蒙古栎(Qm)和茶条槭(Acer ginnala, Ag); 毛榛(Cm)、沙梾(Cornus bretchneideri, Cb)、绣线菊(Spiraea salicifolia)和金银忍冬(Lonicera maackii, Lm); 羊须草(Carex callitrichos, Cc)、龙常草(Diarrhena mandshurica, Dm)、大火草(Anemone tomentosa, At)和玉竹(Polygonatum odoratum, Po)。

1.4 样品采集

样品采集在2013年7月完成。在每个样地中, 每个物种随机选择5株(部分草本一丛视为一株)长势良好的成熟个体。每株随机采集20片完全展开、没有病虫害的外层成熟叶片(油松当年生叶与多年生叶分别采样), 置于两片湿润滤纸之间放入自封袋, 用装有冰袋的保温箱尽快带回实验室。

1.5 性状测定

采集的样品当天进行处理。先将所有植物叶片表面洗净擦干, 剪去叶柄。采用SPAD 502-PLUS手持式叶绿素仪(Konica Minolta, Tokyo, Japan)测定各叶片的叶绿素含量(CC, 单位SPAD, 每片测量3个不同位置), 取平均值。每种叶片5个为一组沿主脉方向整齐叠放, 沿主脉约0.25 cm处均匀选择3个点, 用精度0.01 mm的数显游标卡尺测量厚度值, 取其平均值即为叶厚度(LT, μm)。阔叶树种、灌木、草本的叶片用HP Scanjet G3110扫描仪(Hewlett-Packard Development Company, Palo Alto, USA)扫描后, 应用Image J图像处理软件(National Institutes of Health, Bethesda, USA)计算每种植物叶片的平均叶面积(LA, cm2); 而针叶树种油松叶面积的测量, 首先用数显游标卡尺和直尺测量单松针的长度L (cm)、宽度W (mm)、厚度D (mm), 取平均值, 再根据公式LA (cm2) = 0.01× [π (W/2 + D) + 2W] × L计算得到单针叶面积(李轩然等, 2007)。用电子天平(精度为0.01 g)称取每物种20片叶片的饱和鲜质量(SM, g), 然后在85 ℃下烘干48 h至恒质量, 分别称量干质量(DM, g)。再据此计算其他指标: 比叶面积(SLA, m2·kg-1) = LA/DM; 叶干物质含量(LDMC, g·kg-1) = DM/SM。将烘干叶片磨粉过细筛后, 采用重铬酸钾外加热法测有机碳含量(LCC, mg·g-1); 经H2SO4- H2O2法消煮后, 用凯氏定氮法测定全氮含量(LNC, mg·g-1), 钼锑抗比色法测定全磷含量(LPC, mg·g-1)。

1.6 数据分析

叶性状数值由算术平均值和标准误差来表示。采用单因素方差分析检验氮添加对植物功能性状的影响, 采用Fisher最小显著差数(LSD)法检验各性状在不同处理间的显著性差异。为满足正态分布及方差齐性的假设, 对所有数据进行log10对数转换后, 利用Pearson相关性探讨叶性状间的相关关系。运用主成分分析(PCA)对不同植物多个叶性状进行综合分析, 并据此分析植物在叶经济谱上的分布状况以及氮添加对植物性状变异的解释力。数据分析与图表制作采用SPSS 18.0 (SPSS, Chicago, USA)和Sigma Plot 12.0 (Systat Software, San Jose, USA)软件完成。

2 实验结果

2.1 氮添加对叶片形态、叶绿素含量和叶干物质含量的影响

结果表明, 氮添加对不同物种的叶片形态、叶绿素含量和叶干物质含量有不同程度的影响(图1)。随着氮添加强度的增加, 龙常草和大火草的叶厚度(LT)显著减小, 玉竹的LT则显著增大(图1A); 氮添加对蒙古栎、毛榛、绣线菊、大火草和玉竹的叶面积(LA)有显著影响, 随氮添加它们的叶面积(LA)呈增大趋势(图1B), 但蒙古栎和玉竹的LA在高氮水平回降; 茶条枫比叶面积(SLA)先增后减, 玉竹SLA显著减小(图1C); 蒙古栎、茶条枫等8种植物的叶绿素含量(CC)随氮添加强度增大而显著升高, 其中茶条枫、金银忍冬和玉竹在高氮水平回降(图1D); 油松当年生叶和玉竹的叶干物质含量(LDMC)显著升高, 绣线菊和大火草的LDMC在中氮水平显著低于其他水平(图1E)。将油松当年生叶片(Pt-cy)与往年生叶片(Pt-fy)进行横向比较得出, 当年生叶具有较小的LTLACCLDMC和较大的SLA,LDMC受氮添加的影响更为显著(图1)。

图1   氮添加对叶厚度(LT, A)、叶面积(LA, B)、比叶面积(SLA, C)、叶绿素含量(CC, D)和叶干物质含量(LDMC, E)的影响(平均值±标准误差, n = 3)。Ag, 茶条槭; At, 大火草; Cb, 沙梾; Cc, 羊须草; Cm, 毛榛; Dm, 龙常草; Lm, 金银忍冬; Po, 玉竹; Pt-cy, 油松(当年生叶); Pt-fy, 油松(往年生叶); Qm, 蒙古栎; Ss, 绣线菊。CK, 氮添加浓度为0 kg·hm-2·a-1; LN, 氮添加浓度为50 kg·hm-2·a-1; MN, 氮添加浓度为100 kg·hm-2·a-1; HN, 氮添加浓度为150 kg·hm-2·a-1。不同小写字母表示氮添加处理间在p < 0.05水平上差异显著。

Fig. 1   Effects of nitrogen addition on leaf thickness (LT, A)、leaf area (LA, B)、specific leaf area (SLA, C)、chlorophyll content (CC, D) and leaf dry matter content (LDMC, E)(mean ± SE, n = 3). Ag, Acer ginnala; At, Anemone tomentosa; Cb, Cornus bretchneideri; Cc, Carex callitrichos; Cm, Corylus mandshurica; Dm, Diarrhena mandshurica; Lm, Lonicera maackii; Po, Polygonatum odoratum; Pt-cy, Pinus tabuliformis (current-year leaves); Pt-fy, P. tabuliformis (former-year leaves); Qm, Quercus mongolica; Ss, Spiraea salicifolia. CK, 0 kg·hm-2·a-1 nitrogen; LN, 50 kg·hm-2·a-1 nitrogen; MN, 100 kg·hm-2·a-1 nitrogen; HN, 150 kg·hm-2·a-1 nitrogen. Different lowercase letters indicate the significant difference at 5% level between nitrogen treatments.

2.2 氮添加对叶片碳、氮、磷元素含量的影响

分析表明, 氮添加对植物叶碳含量(LCC)没有显著影响(图2A)。而11种植物的叶氮含量(LNC)均随着氮添加强度增大而显著升高(图2B)。叶磷含量(LPC)对氮添加的响应十分显著, 油松当年生叶LPC在中氮水平显著升高, 在高氮水平又显著降低至对照水平以下; 油松往年生叶无显著变化; 其他植物的LPC均随氮添加强度增加而显著升高, 部分在高氮水平受到抑制(图2C)。蒙古栎和羊须草N:P对氮添加响应不显著, 油松、茶条枫和毛榛N:P受氮添加影响先增后减(氮添加组高于对照组), 大火草N:P显著增大, 其他植物叶片N:P均显著减小(图2D)。油松当年生叶片(Pt-cy)与往年生叶片(Pt-fy)相比, LCC和N:P较小, 而LNCLPC较大, 其中LPC对氮添加响应显著, 而往年生叶LPC无显著变化(图2)。

图2   氮添加对叶碳含量(LCC, A)、叶氮含量(LNC, B)、叶磷含量(LPC, C)及氮磷比(N:P, D)的影响(平均值±标准误差, n = 3)。Ag, 茶条槭; At, 大火草; Cb, 沙梾; Cc, 羊须草; Cm, 毛榛; Dm, 龙常草; Lm, 金银忍冬; Po, 玉竹; Pt-cy, 油松(当年生叶); Pt-fy, 油松(往年生叶); Qm, 蒙古栎; Ss, 绣线菊。CK, 氮添加浓度为0 kg·hm-2·a-1; LN, 氮添加浓度为50 kg·hm-2·a-1; MN, 氮添加浓度为100 kg·hm-2·a-1; HN, 氮添加浓度为150 kg·hm-2·a-1。不同小写字母表示氮添加处理间在p < 0.05水平上差异显著。

Fig. 2   Effects of nitrogen addition on leaf carbon content (LCC, A), leaf nitrogen content (LNC, B), leaf phosphorus content (LPC, C), and N-P ratio (N:P, D)(mean ± SE, n = 3). Ag, Acer ginnala; At, Anemone tomentosa; Cb, Cornus bretchneideri; Cc, Carex callitrichos; Cm, Corylus mandshurica; Dm, Diarrhena mandshurica; Lm, Lonicera maackii; Po, Polygonatum odoratum; Pt-cy, Pinus tabuliformis (current-year leaves); Pt-fy, P. tabuliformis (former-year leaves); Qm, Quercus mongolica; Ss, Spiraea salicifolia. CK, 0 kg·hm-2·a-1 nitrogen; LN, 50 kg·hm-2·a-1 nitrogen; MN, 100 kg·hm-2·a-1 nitrogen; HN, 150 kg·hm-2·a-1 nitrogen. Different lowercase letters indicate the significant difference at 5% level between nitrogen treatments.

2.3 氮添加对不同生活型植物叶片性状的影响比较

比较发现, 不同生活型植物叶性状对氮添加的响应不同(表2)。经氮添加处理, 各生活型植物的LTLCC未发生显著变化, 针叶乔木LALDMC变化显著, 阔叶乔木SLACC变化显著, 4种生活型植物的LNC均显著升高, 阔叶乔木、灌木及草本的LPC也极显著升高, 而灌木、草本植物N:P极显著下降。可以看出, 乔木树种对氮添加的响应包括形态结构指标(SLA、CC等)和N、P元素含量的变化两个方面, 而灌木及草本的响应则仅仅体现在N、P元素含量及配比的变化上。

表2   不同生活型植物叶片性状对氮添加的响应(平均值±标准误差)

Table 2   Variation of leaf traits in different life forms under different levels of nitrogen addition (mean ± SE)

叶性状
Leaf trait
处理
Treatment
生活型 Life form
针叶乔木 Coniferous tree阔叶乔木 Broadleaf tree灌木 Shrub草本 Herb
叶厚度
Leaf thickness (μm)
CK
LN
MN
HN
716.23 ± 15.16a
675.67 ± 27.71a
689.33 ± 27.24a
661.00 ± 17.85a
160.45 ± 10.44a
136.88 ± 10.45a
147.48 ± 5.23a
154.85 ± 9.76a
167.45 ± 18.16a
166.78 ± 13.89a
152.16 ± 7.41a
167.69 ± 16.35a
199.43 ± 28.20a
198.07 ± 23.30a
168.19 ± 20.76a
161.35 ± 17.11a
叶面积
Leaf area (cm2)
CK
LN
MN
HN
6.59 ± 0.14ab
7.13 ± 0.33a
7.18 ± 0.01a
6.45 ± 0.12b
22.91 ± 2.70a
23.92 ± 4.38a
29.33 ± 4.91a
26.49 ± 3.74a
25.59 ± 5.56a
23.79 ± 4.83a
25.97 ± 5.02a
30.02 ± 5.96a
39.75 ± 12.70a
39.36 ± 12.85a
44.70 ± 13.22a
51.44 ± 20.18a
比叶面积
Specific leaf area
(m2·kg-1)
CK
LN
MN
HN
12.42 ± 0.53a
12.62 ± 0.33a
11.54 ± 0.40a
12.67 ± 0.79a
24.32 ± 1.36ab
25.95 ± 0.91a
22.19 ± 1.43ab
21.73 ± 1.19b
36.06 ± 0.86a
35.15 ± 1.27a
32.66 ± 1.69a
32.00 ± 1.56a
34.22 ± 1.12a
33.30 ± 0.92a
32.98 ± 1.16a
31.60 ± 1.09a
叶绿素含量
Chlorophyll content
(SPAD)
CK
LN
MN
HN
45.73 ± 1.36a
40.37 ± 1.34a
44.03 ± 1.95a
42.23 ± 2.87a
43.78 ± 1.91b
46.87 ± 2.17ab
48.75 ± 1.82ab
49.47 ± 1.21a
35.20± 1.46a
37.03 ± 1.60a
39.65 ± 1.83a
39.71 ± 1.83a
30.78 ± 1.50a
31.81 ± 1.74a
34.86 ± 1.47a
33.62 ± 1.43a
叶干物质含量
Leaf dry matter content
(g·kg-1)
CK
LN
MN
HN
391.87 ± 5.43b
421.10 ± 7.77a
420.03 ± 3.91a
416.17 ± 7.31a
381.12 ± 15.37a
377.22 ± 3.14a
395.85 ± 11.58a
384.52 ± 5.94a
284.66 ± 12.41a
287.51 ± 11.45a
290.32 ± 15.77a
295.34 ± 11.52a
236.01 ± 16.81a
241.20 ± 16.62a
213.13 ± 23.40a
249.49 ± 11.31a
叶有机碳含量
Leaf carbon content
(mg·g-1)
CK
LN
MN
HN
587.59 ± 10.05a
541.67 ± 27.75a
553.05 ± 1.01a
548.83 ± 4.88a
473.87 ± 15.47a
486.62 ± 17.09a
481.64 ± 13.02a
487.60 ± 13.23a
484.68 ± 17.84a
450.92 ± 7.49a
454.37 ± 6.12a
468.28 ± 11.67a
452.51 ± 7.97a
441.45 ± 8.30a
449.41 ± 9.61a
450.11 ± 6.43a
叶氮含量
Leaf nitrogen content
(mg·g-1)
CK
LN
MN
HN
12.09 ± 0.39b
13.48 ± 0.47a
11.91 ± 0.50b
12.50 ± 0.38ab
25.25 ± 0.79Cc
26.60 ± 0.32BCb
27.56 ± 0.47Bb
29.78 ± 0.37Aa
22.47 ± 0.61Bc
25.63 ± 0.56Ab
25.81 ± 0.47Ab
26.93 ± 0.45Aa
22.96 ± 0.75Cc
28.04 ± 0.80Bb
28.20 ± 0.62Bb
30.47 ± 0.74Aa
叶磷含量
Leaf phosphorus content
(mg·g-1)
CK
LN
MN
HN
1.08 ± 0.12a
1.02 ± 0.10a
1.06 ± 0.11a
0.98 ± 0.07a
1.60 ± 0.11ABb
1.48 ± 0.05Bb
1.70 ± 0.07ABab
1.92 ± 0.17Aa
1.17 ± 0.05Cc
1.41 ± 0.03Bb
1.71 ± 0.07Aa
1.75 ± 0.05Aa
1.26 ± 0.04Bc
1.68 ± 0.08Ab
1.91 ± 0.08Aa
1.86 ± 0.08Aab
叶氮磷比
leaf N:P
CK
LN
MN
HN
12.54 ± 1.11a
14.75 ± 1.14a
12.44 ± 0.89a
14.40 ± 0.80a
16.39 ± 0.96a
18.16 ± 0.46a
16.44 ± 0.52a
16.70 ± 1.08a
19.86 ± 0.89Aa
18.47 ± 0.61ABa
15.80 ± 0.72Bb
15.77 ± 0.59Bb
18.50 ± 0.50Aa
17.26 ± 0.66ABab
15.16 ± 0.44Bc
16.95 ± 0.59ABb

CK, 0 kg·hm-2·a-1 nitrogen; LN, 50 kg·hm-2·a-1 nitrogen; MN, 100 kg·hm-2·a-1 nitrogen; HN, 150 kg·hm-2·a-1 nitrogen. Different capital or lowercase letters indicate the significant difference at 1% level or 5% level among treatments. Data with significant difference were in bold.CK, 氮添加浓度为0 kg·hm-2·a-1; LN, 氮添加浓度为50 kg·hm-2·a-1; MN, 氮添加浓度为100 kg·hm-2·a-1; HN, 氮添加浓度为150 kg·hm-2·a-1。不同小写字母(a, b, c)和大写字母(A, B, C)分别表示不同氮添加水平下叶片性状在p < 0.05水平和p < 0.01水平上差异显著。差异显著的数据用加粗字体显示。

新窗口打开

2.4 氮添加对植物叶性状及其相互关系的影响

Pearson相关性分析表明, 叶性状之间普遍存在着显著或极显著相关关系(表3)。SLALTCCLDMCLCC极显著负相关, 与LALNCLPC和N:P极显著正相关, 相关系数值LT > LNC > LCC > LDMC > LPC > LA > CC > N:P。N:P与LTLPC呈极显著负相关, 与LASLALNC极显著正相关, 与CC显著正相关, 相关系数值表现为LPC > LT > SLALNC > LA > CC。氮添加处理强度(N-Tr.)与SLA和N:P显著负相关, 与CCLNCLPC显著正相关, 相关系数数值LPC > LNC > CC > N:P > SLA; 所有叶性状与物种(Sp.)均显著相关, 且氮添加处理和物种对所有叶性状的交互作用(N-Tr. × Sp.)显著。

表3   叶性状间Pearson相关系数及氮添加处理和物种的影响效应分析

Table 3   Pearson correlation coefficients among leaf traits and the effects of N-treatment and species

LTLASLACCLDMCLCCLNCLPCN:PN-Tr.Sp.Sig. (N-Tr.× Sp.)
LT1-0.06-0.50**0.00
LA-0.17**10.060.37**0.00
SLA-0.73**0.42**1-0.12*0.71**0.01
CC0.100.25**-0.34**10.18**-0.30**0.00
LDMC0.11-0.43**-0.56**0.35**10.04-0.72**0.00
LCC0.44**-0.55**-0.64**0.22**0.52**1-0.04-0.58**0.00
LNC-0.73**0.56**0.71**0.13*-0.37**-0.60**10.25**0.62**0.00
LPC-0.46**0.42**0.50**0.03-0.40**-0.54**0.71**10.34**0.39**0.00
N:P-0.36**0.17**0.29**0.14*0.08-0.060.29**-0.41**1-0.17**0.26**0.00

CC, chlorophyll content; LA, leaf area; LCC, leaf carbon content; LDMC, leaf dry matter content; LNC, leaf nitrogen content; LPC, leaf phosphorus content; LT, leaf thickness; N:P, N-P ratio; SLA, specific leaf area. N-Tr., N addition treatment; Sp., Species; Sig. (N-Tr.×Sp.), the significance of N-Tr. by Sp. interactions. Data were log10-transformed before the analysis. *, p < 0.05; **, p < 0.01.CC, 叶绿素含量; LA, 叶面积; LCC, 叶有机碳含量; LDMC, 叶干物质含量; LNC, 叶氮含量; LPC, 叶磷含量; LT, 叶厚度; N:P, 氮磷比; SLA, 比叶面积。N-Tr., 氮添加处理; Sp., 物种; Sig. (N-Tr. × Sp.), 氮添加处理与物种交互作用的显著性。分析前对数据进行了log10转化。*, p < 0.05; **, p < 0.01。

新窗口打开

对于与构成叶经济谱相关的关联性状LNC vs. SLALPC vs. SLALNC vs. LTLPC vs. LT进行分析发现, 氮添加强度增大时, LNCLPC vs. SLALT的拟合线斜率均发生了明显变化, 且LPC vs. SLALT的斜率变化幅度大于LNC vs. SLALT的变化幅度(图3)。此外, 随着氮添加浓度的增加, LNC:SLALPC:SLALNC:LTLPC:LT的值均显著增大(表4), 说明这4组相关关系受氮添加影响显著。

图3   氮添加对叶氮含量(LNC)与比叶面积(SLA) (A)、LNC与叶厚度(LT) (B)、叶磷含量(LPC)与SLA (C)和LPCLT (D)关系的影响。CK, 氮添加浓度为0 kg·hm-2·a-1; LN, 氮添加浓度为50 kg·hm-2·a-1; MN, 氮添加浓度为100 kg·hm-2·a-1; HN, 氮添加浓度为150 kg·hm-2·a-1

Fig. 3   Effects of nitrogen addition on leaf nitrogen content (LNC) vs. specific leaf area (SLA) (A), LNC vs. leaf thickness (LT) (B), leaf phosphorus content (LPC) vs. SLA (C) and LPC vs. LT (D). CK, 0 kg·hm-2·a-1 nitrogen; LN, 50 kg·hm-2·a-1 nitrogen; MN, 100 kg·hm-2·a-1 nitrogen; HN, 150 kg·hm-2·a-1 nitrogen.

表4   氮添加处理对叶氮含量、叶磷含量与比叶面积、叶厚度比值影响的单因素方差分析(平均值±标准误差, n = 36)

Table 4   Results of one-way ANOVA for the effect of nitrogen addition on the ratio of leaf nitrogen content (LNC) over specific leaf area (SLA) and leaf thickness (LT), and leaf phosphorus content (LPC) over specific leaf area (SLA) and leaf thickness (LT) (mean ± SE, n = 36)

处理
Treatment
叶氮含量/比叶面积
LNC:SLA
叶磷含量/比叶面积
LPC:SLA
叶氮含量/叶厚度
LNC:LT
叶磷含量/叶厚度
LPC:LT
CK0.74 ± 0.022d0.05 ± 0.002b0.12 ± 0.007c0.01 ± 0.000b
LN0.85 ± 0.018c0.05 ± 0.002b0.15 ± 0.009bc0.01 ± 0.001b
MN0.92 ± 0.027b0.06 ± 0.002a0.16 ± 0.009ab0.01 ± 0.001a
HN1.05 ± 0.033a0.07 ± 0.003a0.18 ± 0.012a0.01 ± 0.001a

CK, 0 kg·hm-2·a-1 nitrogen; LN, 50 kg·hm-2·a-1 nitrogen; MN, 100 kg·hm-2·a-1 nitrogen; HN, 150 kg·hm-2·a-1 nitrogen. Different lowercase letters indicate the significant difference at 5% level among treatments.CK, 氮添加浓度为0 kg·hm-2·a-1; LN, 氮添加浓度为50 kg·hm-2·a-1; MN, 氮添加浓度为100 kg·hm-2·a-1; HN, 氮添加浓度为150 kg·hm-2·a-1。不同字母表示不同氮添加强度下叶片性状差异显著(p < 0.05)。

新窗口打开

2.5 氮添加对植物在叶经济谱中位置的影响

主成分分析结果表明, 第一主成分和第二主成分的方差解释率分别为52.04%和15.62%, 合计67.66% (表5)。其中主成分1主要与LNC、SLA、LCC、LPCLA等因子相关, 相当于LES理论中的“投资-收益”策略轴, 而主成分2与N:P、LPC、LT等相关度较高(图4A)。其中, 氮添加处理(N-Tr.)对第一主成分的解释率约为-19.9%, 与其他因子相比解释力偏低; 而对于第二主成分的解释率则高达-56.8%, 解释力仅次于N:P。这说明N-Tr.促使植物叶性状在第一、第二主成分上的因子得分均有所下降, 但主要作用于第二主成分。此外, 进一步计算得到氮添加在针叶乔木、阔叶乔木、灌木、草本中第一、第二主成分载荷向量之和的绝对值分别是0.77、0.73、0.78和0.83, 说明它对4种生活型植物叶片性状的影响力大小为草本>灌木>针叶乔木>阔叶乔木。

表5   叶性状在主成分分析中的载荷

Table 5   Loadings of leaf traits in principal components analyses

载荷
Loading
叶厚度
Leaf thickness
叶面积
Leaf area
比叶面积
Specific leaf
area
叶绿素含量
Chlorophyll
content
叶干物质含量
Leaf dry matter
content
叶碳含量
Leaf carbon
content
叶氮含量
Leaf nitrogen
content
叶磷含量
Leaf phosphorus
content
氮磷比
N-P ratio
氮添加处理
N addition treatment
主成分1
Principal component 1
0.70-0.64-0.830.090.580.80-0.91-0.80-0.11-0.20
主成分2
Principal component 2
-0.400.080.160.320.330.100.24-0.430.94-0.57

The total explanation of two principal component is 67.66%, the explanation of principal component 1 and component 2 are 52.04% and 15.62%, respectively.两个主成分总的解释力为67.66%, 其中主成分1和主成分2解释力分别为52.04%和15.62%。

新窗口打开

图4   天然油松林主要植物的叶性状主成分分析结果: 天然油松林主要植物的叶性状载荷(A)、因子得分(B)及对氮添加的响应(C)。分析前对数据进行了log10转化。CC, 叶绿素含量; LA, 叶面积; LCC, 叶有机碳含量; LDMC, 叶干物质含量; LNC, 叶氮含量; LPC, 叶磷含量; LT, 叶厚度; N:P, 氮磷比; SLA, 比叶面积。Ag, 茶条槭; At, 大火草; Cb, 沙梾; Cc, 羊须草; Cm, 毛榛; Dm, 龙常草; Lm, 金银忍冬; Po, 玉竹; Pt-cy, 油松(当年生叶); Pt-fy, 油松(往年生叶); Qm, 蒙古栎; Ss, 绣线菊。C, 针叶乔木; B, 阔叶乔木; H, 草本; S, 灌木。CK, 氮添加浓度为0 kg·hm-2·a-1; LN, 氮添加浓度为50 kg·hm-2·a-1; MN, 氮添加浓度为100 kg·hm-2·a-1; HN, 氮添加浓度为150 kg·hm-2·a-1

Fig. 4   Results of principle component analyses (PCA) of leaf traits across common species in the natural Pinus tabuliformis forest. Loadings of leaf traits across primary species in the natural forest of P. tabuliformis (A), factor scores of different species (B) and responses to nitrogen addition (C). Data were log10-transformed before analysis. CC, chlorophyll content; LA, leaf area; LCC, leaf carbon content; LDMC, leaf dry matter content; LNC, leaf nitrogen content; LPC, leaf phosphorus content; LT, leaf thickness; N:P, N-P ratio; SLA, specific leaf area. Ag, Acer ginnala; At, Anemone tomentosa; Cb, Cornus bretchneideri; Cc, Carex callitrichos; Cm, Corylus mandshurica; Dm, Diarrhena mandshurica; Lm, Lonicera maackii; Po, Polygonatum odoratum; Pt-cy, Pinus tabuliformis (current-year leaves); Pt-fy, P. tabuliformis (former-year leaves); Qm, Quercus mongolica; Ss, Spiraea salicifolia. C, coniferous tree; B, broadleaf tree; H, herb; S, shrub. CK, 0 kg·hm-2·a-1 nitrogen; LN, 50 kg·hm-2·a-1 nitrogen; MN, 100 kg·hm-2·a-1 nitrogen; HN, 150 kg·hm-2·a-1 nitrogen.

分析11种植物的因子得分情况(图4B), 油松的第一因子得分最高, 其次是金银忍冬和羊须草, 得分最低的是玉竹和大火草, 即油松、金银忍冬和羊须草具有较高的LTLDMCLCC, 而玉竹和大火草等具有更大的LASLALNCLPC。第二因子得分最高的是金银忍冬、绣线菊和毛榛, 与得分最低的沙梾和油松当年生叶相比具有更高的CC和N:P。对各生活型植物叶性状在不同氮添加强度下的因子得分进行比较, 经过氮添加处理的阔叶乔木(B)、灌木(S)及草本(H)的位置相比于对照向左下方发生偏移(第一因子、第二因子得分均降低); 而针叶乔木(C)的施氮组与对照组相比则向左上方偏移(第一因子得分下降, 而第二因子得分升高) (图4C)。也就是说, 模拟氮沉降对植物叶片的经济投资策略产生了一定影响。

3 讨论

3.1 氮添加对植物叶片性状的影响

本试验表明, 氮添加对天然油松林11种植物的叶片性状产生了显著影响: 龙常草、大火草、玉竹的LTSLA, 蒙古栎、毛榛、绣线菊、大火草和玉竹LA, 以及油松当年生叶、绣线菊、大火草和玉竹的LDMC都产生了显著性差异, 9种植物叶片的N:P随氮添加浓度发生显著变化, 全部植物的LNC与大多数植物的CCLPC显著升高。

一般认为, 氮添加会导致SLA的增大(Rose et al., 2013), 因为氮添加促进了植物的生长而提高了植物的光合效率(Ordoñez et al., 2009)。万宏伟等(2008)和黄菊莹等(2009)研究发现, 羊草(Leymus chinensis)等几种草本植物的SLA因氮添加而显著增大; Bassin等(2009)的氮添加实验也表明, 虽然大多数植物叶片的比叶重(LMA, SLA的倒数)无显著变化, 但两种植物叶片的LMA显著降低(即SLA显著增大)。然而, 并非所有研究结果都是如此, 赵新风等(2014)发现, 钠猪毛菜(Salsola nitraria)受氮添加影响, SLA呈下降趋势, LDMC呈增大趋势, 与上述结果相反。研究表明, 较小的SLA和较高的LDMC使得植物具有更好的资源获取能力、碳氮固持效率及水分利用效率, 在逆境中更具优势(Wilson et al., 1999; Aerts & Chapin, 2000; Zhang et al., 2007)。因此, 这两项指标常被用来指示植物对环境变化的响应策略, SLA的减小和LDMC的增大往往意味着环境条件的恶化, 因为植物开始趋向于保守的生长策略(Rose et al., 2013)。本研究结果表明, 大部分植物的SLALDMC并没有表现出显著差异, 只有茶条枫和玉竹的SLA呈显著减小趋势, 油松当年生叶和玉竹叶片的LDMC显著增大, 与赵新风等(2014)的结果相一致。这说明不同植物对氮添加的响应不同。并且, 氮添加对油松、茶条枫和玉竹这3种植物的生长是不利的。

另外, 本研究中11种植物叶片的CCLNC增大、LCC保持不变的趋势与前人研究结果一致, 但多数植物LPC显著上升、部分植物叶片N:P显著下降的变化趋势却与其他研究中LPC降低而N:P上升的结果相反(李德军等, 2005; Rose et al., 2013; 刘洋等, 2013)。有研究表明, 在N供应受限制的环境中, 氮添加会提高土壤中活性氮的含量, 缓解N的限制作用, 促进植物的生长, 并使N:P显著增大、P逐渐短缺(陈凌云, 2010; 刘洋等, 2013)。遆超普和颜晓元(2010)提供的数据显示, 山西的铵态氮、硝态氮沉降量与排放值的比值分别为1.58和0.26, 氮湿沉降总量由1980年的0.20 Tg·a-1增长到2007年的0.47 Tg·a-1, 超过全国平均水平, 并且山西是全国氮沉降增速最快的地区。以本试验中对照组数值做参考, 11种植物LNC的平均值为21.37 mg·g-1, 略高于全球和全国陆生植物的平均水平(分别为20.1和20.2 mg·g-1); LPC的平均值为1.26 mg·g-1, 低于全球、全国水平(分别为1.8和1.46 mg·g-1); N:P的均值为17.6,明显高于全球陆生植物的平均水平13.8及全国平均水平16.3 (Reich & Oleksyn, 2004; Han et al., 2005)。以叶片N:P作为判断植物生长限制因子的指标, 其中N:P < 14表明氮限制、N:P > 16代表磷限制、处于二者之间则为氮磷共同限制或者均不受限(Aerts & Chapin, 2000; Güsewell, 2004), 相对于全球和全国而言, 山西的植物表现为较明显的磷元素缺乏, 而非受到氮素限制。因此, 施加过量的氮素可能不会像其他在氮素缺乏地区进行的实验一样起到促进植物生长的作用, 还可能不利于植物的生长(这与SLALDMC预测的结果一样)。目前, 氮添加导致植物叶片LPC显著上升和N:P显著下降的类似研究结果还未见报道。据推测, 这种现象的出现可能与本试验地区土壤氮素供应充足而磷供应缺乏的自然条件、较高的施氮水平, 以及研究对象本身的物种特性不同于其他实验有关。而关于叶片N:P显著下降的罕见趋势, 可以从N:P与LNCLPC相关关系的研究中找到答案。大量研究表明, 植物LPC与N:P存在显著的负相关关系, 而LNC与N:P没有显著的相关性或相关系数很小(Güsewell, 2004; He et al., 2006, 2008)。本研究中, N:P虽然与LNCLPC都存在极显著的相关性, 但前者的相关系数远小于后者。由此可见, 植物叶片N:P主要受P含量的影响。因此, 虽然绝大部分植物的LNCLPC受到氮添加的影响均显著上升, 但其中5种植物叶片的N:P因主要受到LPC升高的影响而显著下降。

3.2 植物叶片应对环境变化的资源权衡策略

黄海侠等(2013)发现, 同种植物叶片的物质分配策略在异龄叶间存在显著差异: 当年生叶的SLA显著大于往年生叶、LDMC比往年生叶小, 且当年生叶性状的变异系数往往高于往年生叶。本研究也从油松当年生叶和往年生叶的性状比较中, 发现了类似的结果。这是因为叶片在生长初期往往需要较多的资源来优化同化器官、最大程度地获取光资源(Sultan, 2000); 但在老化过程中, 内部养分元素(如氮和磷)发生迁移且光合速率降低(Bleecker, 1998; Ono et al., 2001), 结构充分形成并固化, 对环境变化的响应就变得相对迟缓(Oguchi et al., 2003); 老叶为忍受多年环境变化, 在机械和防御构造上的额外投资较多(Chabot & David, 1982)。此外, 新叶对环境变化敏感, 轻微的环境扰动就能诱导其性状发生改变, 表现出较强的可塑性(Wyka et al., 2000; Frak et al., 2001)。因此, 叶片的物质投资和分配格局是随叶龄变化的, 这是植物适应环境的一个重要途径(Silvertown et al., 1997)。

通常, 具有较高SLALNCLPC的叶片, 往往具有较高的光合速率和生长速率, 并且更薄、更可口、更易啃食, 寿命也更短, 具有更大的LTLDMC的植物则相反(Crawley, 1983; Kikuzawa, 1991)。大量研究表明, 代表光能利用效率的SLA (Westoby et al., 2002)与代表光合速率的LNC (Grime, 2001)及间接反映生长速率的N:P呈极显著正相关关系, 而与代表物理支撑结构的LTLDMCLCC呈显著负相关关系(Wilson et al., 1999)。这是因为植物在特定的环境条件下可利用的资源总量是有限的, 从养分配置角度讲, 若一方面的资源投入大, 则其他方面的投入就要减少, 从而使植物性状之间保持普遍的协同变化关系(Shipley et al., 2006)。本研究结果与上述结论相一致, SLALTCCLDMCLCC极显著负相关, 与LALNCLPC和N:P极显著正相关, 作为一个核心性状将所有性状联系起来, 体现出叶片在功能效率和结构建成两个方面“此消彼长”的资源权衡投资策略(Wright et al., 2005; Ordoñez et al., 2009)。从各种植物在主成分分析(PCA)第一主成分轴上的分布情况来看, 得分最高的油松和得分最低的玉竹、大火草, 分别具有叶厚、光合能力弱、生长缓慢、寿命长和叶薄、光合能力强、生长快速、寿命短的特点, 恰好反映出两种完全不同的生存策略, 验证了LES理论(陈莹婷和许振柱, 2014)。此外, LNCLPCSLALT之间相关拟合曲线的斜率在氮添加影响下均显著增大。类似地, Fonseca等(2000)发现SLALA的关系在生境内部和不同生境之间都不相同; Wright等(2001)发现SLALNC的相关关系在不同水位处理和氮素处理下都有显著差异; Pensa等(2010)也发现, 同种植物LNC与叶寿命的关系随着纬度梯度不断增长。这说明, 植物叶片性状之间的相关关系并非一成不变, 而是随生境变化而变化的。因此, 我们在试图利用叶性状普遍存在的相关关系, 由部分易获取的功能性状来预测其他相对难以测量的性状值及变化趋势时, 需要综合考虑局域气候-地理特点或者综合的植物功能型属性。

PCA结果表明, 植物在多元特征空间中的位置在氮添加的作用下向左侧发生了较为明显的横向迁移(图4C), 即趋向于叶片更大、更薄、生长更快和寿命更短的“快速投资-收益型”策略。而在垂直方向上, 阔叶乔木、灌木和草本植物均朝下移动, 针叶乔木则向上移动, 说明不同植物对于同一环境因素的变化有不同的响应策略, 即不同植物叶片特征对环境变化的响应是异质的(Donovan et al., 2011; Wright & Sutton-Grier, 2012)。Wright等(2004)指出, 全球尺度上的温度和降水因子对植物的影响总体上是十分微弱的, 植物通过不断调节自身的功能性状及其权衡关系, 稳固地占据着一定的生态位空间, 从而部分抵消了温度和降水变化所带来的影响。在本研究中, 虽然11种植物在特征空间中都发生了明显的位移, 但物种之间的相对位置保持不变, 在群落水平上维持了叶片经济的内稳态。并且, 在综合各项性状的大幅度变异之后, 植物在第一主轴(即LES)上的位移仅为-0.20, 说明在局域范围内, 氮添加对植物生存策略的影响也非常有限。也就是说, 尽管局域种库内植物对环境的响应是多样化的, 但多种多样的性状变异并不足以改变植物在第一主轴上的相对位置, 也没有大幅度改变整个群落在第一主轴上的位置, 植物依旧维持着相当稳定的生态位。因此, 无论是在全球尺度还是局域尺度上, 当环境发生变化时, 植物将自发改变生活策略, 发挥权衡作用对资源进行重新调配, 在很大程度上维持种间相对位置和群落整体结构的稳定性。这种调节和适应机制是植物本身所固有的, 并不依赖于环境气候变化而产生、发展, 也不受环境因素所主导, 故第一主轴的形成是基于植物根本性的生态学权衡策略独立演变而来的(Wright et al., 2004, 2005; Shipley et al., 2006; Wright & Sutton-Grier, 2012)。

综上所述, 氮添加对山西太岳山天然油松林叶片性状的影响显著, 不同类型、不同年龄的植物叶片对氮添加的响应各不相同, 叶性状间的相关性随氮添加强度而改变, 植物叶经济谱及其涉及的叶性状权衡关系被证实存在, 植物叶片在氮添加的影响下更趋向于“快速投资-收益型”策略, 第一主轴的形成基于植物的生态学权衡策略而不受环境因素主导。本研究对于揭示植物性状对大气氮沉降等气候变化的响应与适应机制具有一定意义。

致谢 感谢中国科学院地理研究所“可桢杰出青年学者”项目(2013RC102)资助。感谢生物多样性与生态工程教育部重点实验室开放课题资助。感谢山西太岳山灵空山林场对研究工作的支持。感谢北京林业大学资源与生态系统过程北京市重点实验室赵秀海老师对实验方案的指导, 以及秦书琪、马冰倩、刘宇三位同学在野外工作中给予的帮助。

The authors have declared that no competing interests exist.

作者声明没有竞争性利益冲突.


参考文献

[1] Aber JD, McDowell W, Nadelhoffer K, Magill A, Berntsen G, Kamakea M, McNulty S, Currie W, Rustad L, Fernandez I (1998).

Nitrogen saturation in temperate forest ecosystems: Hypothesis revisited

.BioScience, 48, 921-934.

[本文引用: 1]     

[2] Aber JD, Nadelhoffer KJ, Steudler P, Melillo JM (1989).

Nitrogen saturation in northern forest ecosystems

. BioScience, 39, 378-386.

[本文引用: 1]     

[3] Adams MA, Ineson P, Binkley D, Cadisch G, Tokuchi N, Scholes M, Hicks K (2004).

Soil functional responses to excess nitrogen inputs at global scales

. AMBIO, 33, 530-536.

[本文引用: 1]     

[4] Aerts R, Chapin SF (2000).

The mineral nutrition of wild plants revisited: A re-evaluation of processes and patterns

.Advances in Ecological Research, 30, 1-67.

[本文引用: 2]     

[5] Albert CH, Grasseina F, Schurrd FM, Vieilledente G, Violle C (2011).

When and how should intraspecific variability be considered in trait-based plant ecology?

Perspectives in Plant Ecology, Evolution, and Systematics, 13, 217-225.

[本文引用: 1]     

[6] Albert CH, Thuiller W, Yoccoz NG, Douzet R, Aubert S, Lavorel S (2010a).

A multi-trait approach reveals the structure and the relative importance of intra- vs. interspecific variability in plant traits

.Functional Ecology, 24, 1192-1201.

[本文引用: 1]     

[7] Albert CH, Thuiller W, Yoccoz NG, Soudant A, Boucher F, Saccone P, Lavorel S (2010b).

Intraspecific functional variability: Extent, structure and sources of variation

.Journal of Ecology, 98, 603-613.

[本文引用: 1]     

[8] Bassin S, Werner RA, Sörgel K, Volk M, Buchmann N, Fuhrer J (2009).

Effects of combined ozone and nitrogen deposition on the in situ properties of eleven key plant species of a subalpine pasture

.Oecologia, 158, 747-756.

[9] Bleecker AB (1998).

The evolutionary basis of leaf senescence: Method to the madness?

Current Opinion in Plant Biology, 1, 73-78.

[本文引用: 1]     

[10] Blonder B, Violle C, Enquist BJ, Cornelissen H (2013).

Assessing the causes and scales of the leaf economics spectrum using venation networks in Populus tremuloides

. Journal of Ecology, 101, 981-989.

[11] Bobbink R, Hornung M, Roelofs JGM (1998).

The effects of air-borne nitrogen pollutants on species diversity in natural and semi-natural European vegetation

.Journal of Ecology, 86, 717-738.

[本文引用: 1]     

[12] Chabot BF, David JH (1982).

The ecology of leaf life spans

.Annual Review of Ecology and Systematics, 13, 229-259.

[本文引用: 1]     

[13] Chen LY (2010).

Effects of N, P Addition on N:P Stoichiometry of Different Functional Groups in Potentilla fruticosa Community in a Sub-alpine Meadow. Master degree dissertation, Lanzhou University, Lanzhou

. (in Chinese with English abstract)
[陈凌云 (2010).

添加氮磷对亚高寒草甸金露梅群落各功能群化学计量学特征的影响

. 硕士学位论文, 兰州大学, 兰州.]

[本文引用: 1]     

[14] Chen YT, Xu ZZ (2014).

Review on research of leaf economics spectrum

. Chinese Journal of Plant Ecology, 38, 1135-1153. (in Chinese with English abstract)
[陈莹婷, 许振柱 (2014).

植物叶经济谱的研究进展

. 植物生态学报, 38, 1135-1153.]

[本文引用: 2]     

[15] Clark CM, Tilman D (2008).

Loss of plant species after chronic low-level nitrogen deposition to prairie grasslands

.Nature, 451, 712-715.

[本文引用: 1]     

[16] Comelissen JHC, Cerabolini B, Castro-Diez P, Villar-Salvador P, Montserrat-Marti G, Puyravaud JP, Maestro M, Werger MJA, Aerts R (2003).

Functional traits of woody plants: Correspondence of species rankings between field adults and laboratory-grown seedlings?

Journal of Vegetation Science, 14, 311-322.

[本文引用: 1]     

[17] Crawley MJ (1983).

Herbivory: The Dynamics of Animal-Plant Interactions

. Blackwell Scientific Publications, Oxford, UK.

[本文引用: 1]     

[18] Donovan LA, Maherali H, Caruso CM, Huber H, de Kroon H (2011).

The evolution of the worldwide leaf economics spectrum

. Trends in Ecology Evolution, 26, 88-95.

[本文引用: 1]     

[19] Fang YT, Mo JM, Zhou GY, Xue JH (2005).

Response of diameter at breast height increase to N addition in forests of Dinghushan Biosphere Reserve

.Journal of Tropical and Subtropical Botany, 13(3), 198-204. (in Chinese with English abstract)
[方运霆, 莫江明, 周国逸, 薛璟花 (2005).

鼎湖山主要森林类型植物胸径生长对氮沉降增加的初期响应

. 热带亚热带植物学报, 13(3), 198-204.]

[本文引用: 1]     

[20] Fenn ME, Poth MA, Johnson DW (1996).

Evidence for nitrogen saturation in the San Bernardino Mountains in Southern California

.Forest Ecology and Management, 82, 211-230.

[本文引用: 1]     

[21] Fonseca CR, Overton JM, Collins B, Westoby M (2000).

Shifts in trait combinations along rainfall and phosphorus gradients

.Journal of Ecology, 88, 964-977.

[22] Frak E, Roux XL, Millard P, Dreyer E, Jaouen G, Saint-Joanis B, Wendler R (2001).

Changes in total leaf nitrogen and partitioning of leaf nitrogen drive photosynthetic acclimation to light in fully developed walnut leaves

.Plant, Cell & Environment, 24, 1279-1288.

[本文引用: 1]     

[23] Freschet GT, Cornelissen JHC, van Logtestijn RSP, Aerts R (2010).

Evidence of the ‘plant economics spectrum’ in a subarctic flora

.Journal of Ecology, 98, 362-373.

[本文引用: 1]     

[24] Grime JP (2001).

Plant Strategies, Vegetation Processes and Ecosystem Properties

. John Wiley & Sons, Chichester, UK.

[本文引用: 1]     

[25] Güsewell S (2004).

N:P ratios in terrestrial plants: Variation and functional significance

.New Phytologist, 164, 243-266.

[本文引用: 2]     

[26] Han WX, Fang YT, Guo DL, Zhang Y (2005).

Leaf nitrogen and phosphorus stoichiometry across 753 terrestrial plant species in China

. New Phytologist, 168, 377-385.

[本文引用: 1]     

[27] He JS, Fang JY, Wang ZH, Guo DL, Flynn DFB, Geng Z (2006).

Stoichiometry and large-scale patterns of leaf carbon and nitrogen in the grassland biomes of China

.Oecologia, 149, 115-122.

[本文引用: 1]     

[28] He JS, Han XG (2010).

Ecological stoichiometry: Searching for unifying principles from individuals to ecosystems

. Chinese Journal of Plant Ecology, 34, 2-6. (in Chinese with English abstract)
[贺金生, 韩兴国 (2010).

生态化学计量学: 探索从个体到生态系统的统一化理论

. 植物生态学报, 34, 2-6.]

[本文引用: 1]     

[29] He JS, Wang L, Flynn DFB, Wang XP, Ma WH, Fang JY (2008).

Leaf nitrogen:phosphorus stoichiometry across Chinese grassland biomes

.Oecologia, 155, 301-310.

[本文引用: 1]     

[30] Holland EA, Dentener FJ, Braswell BH, Sulzman JM (1999).

Contemporary and pre-industrial global reactive nitrogen budgets

.Biogeochemistry, 46, 7-43.

[本文引用: 1]     

[31] Huang HX, Yang XD, Sun BW, Zhang ZH, Yan ER (2013).

Variability and association of leaf traits between current-year and former-year leaves in evergreen trees in Tiantong, Zhejiang, China

.Chinese Journal of Plant Ecology, 37, 912-921. (in Chinese with English abstract)
[黄海侠, 杨晓东, 孙宝伟, 张志浩, 阎恩荣 (2013).

浙江天童常绿植物当年生与往年生叶片性状的变异与关联

. 植物生态学报, 37, 912-921.]

[32] Huang JY, Yuan ZY, Li LH (2009).

Changes in N, P and specific leaf area of green leaves of Leymus chinensis along nitrogen, phosphorus and water gradients

.Chinese Journal of Plant Ecology, 33, 442-448. (in Chinese with English abstract)
[黄菊莹, 袁志友, 李凌浩 (2009).

羊草绿叶氮、磷浓度和比叶面积沿氮、磷和水分梯度的变化

. 植物生态学报, 33, 442-448.]

[本文引用: 1]     

[33] Kikuzawa K (1991).

A cost-benefit analysis of leaf habit and leaf longevity of trees and their geographical pattern

. The American Naturalist, 138, 1250-1263.

[本文引用: 1]     

[34] Laughlin DC (2011).

Nitrification is linked to dominant leaf traits rather than functional diversity

.Journal of Ecology, 99, 1091-1099.

[本文引用: 1]     

[35] Leishman MR, Haslehurst T, Ares A, Baruch Z (2007).

Leaf trait relationships of native and invasive plants: Community- and global-scale comparisons

.New Phytologist, 176, 635-643.

[36] Li DJ, Mo JM, Fang YT, Li ZA (2005).

Effects of simulated nitrogen deposition on biomass production and allocation in Schima superba and Cryptocarya concinna seedlings in subtropical China

.Acta Phytoecologica Sinica, 29, 543-549. (in Chinese with English abstract)
[李德军, 莫江明, 方运霆, 李志安 (2005).

模拟氮沉降对南亚热带两种乔木幼苗生物量及其分配的影响

. 植物生态学报, 29, 543-549.]

[本文引用: 2]     

[37] Li HS, Wang JS, Liu X, Jiang SS, Zhang CY, Zhao XH (2014).

Effects and its sustained effect of simulated nitrogen deposition on soil respiration in Pinus tabulaeformis forests in the Taiyue Mountain, China

.Acta Scientiae Circumstantiae, 34, 238-249. (in Chinese with English abstract)
[李化山, 汪金松, 刘星, 蒋思思, 张春雨, 赵秀海 (2014).

模拟氮沉降对太岳山油松林土壤呼吸的影响及其持续效应

. 环境科学学报, 34, 238-249.]

[38] Li XR, Liu QJ, Cai Z, Ma ZQ (2007).

Specific leaf area and leaf area index of conifer plantations in Qianyanzhou Station of subtropical China. Journal of Plant Ecology (

Chinese Version), 31, 93-101. (in Chinese with English abstract)
[李轩然, 刘琪璟, 蔡哲, 马泽清 (2007).

千烟洲针叶林的比叶面积及叶面积指数

. 植物生态学报, 31, 93-101.]

[本文引用: 1]     

[39] Liu XJ, Zhang Y, Han WX, Tang AH, Shen JL, Cui ZL, Vitousek P, Erisman JW, Goulding K, Christie P, Fangmeier A, Zhang FS (2013).

Enhanced nitrogen deposition over China

.Nature, 494, 459-462.

[本文引用: 1]     

[40] Liu Y, Zhang J, Chen YM, Chen L, Liu Q (2013).

Effect of nitrogen and phosphorus fertilization on biomass allocation and C:N:P stoichiometric characteristics of Eucalyptus grandis seedlings

.Chinese Journal of Plant Ecology, 37, 933-941. (in Chinese with English abstract)
[刘洋, 张健, 陈亚梅, 陈磊, 刘强 (2013).

氮磷添加对巨桉幼苗生物量分配和C:N:P化学计量特征的影响

. 植物生态学报, 37, 933-941.]

[本文引用: 2]     

[41] Magill AH, Aber JD, Currie WS, Nadelhoffer KJ, Martin ME, McDowell WH, Melillo JM, Steudler P (2004).

Ecosystem response to 15 years of chronic nitrogen additions at the Harvard Forest LTER, Massachusetts, USA

.Forest Ecology and Management, 196, 7-28.

[本文引用: 1]     

[42] Maire V, Gross N, Hill D, Martin R, Wirth C, Wright IJ, Soussana JF (2013).

Disentangling coordination among functional traits using an individual-centred model: Impact on plant performance at intra- and inter-specific levels

. PLOS ONE, 8, e77372. doi: 10.1371/journal.pone.00773-72.eCollection 2013.

[43] McGill BJ, Enquist BJ, Weiher E, Westoby M (2006).

Rebuilding community ecology from functional traits

.Trends in Ecology Evolution, 21, 178-185

[本文引用: 1]     

[44] Nihlgard B (1985).

The ammonium hypothesis—An additional explanation to the forest dieback in Europe

.AMBIO, 14, 2-8

[本文引用: 1]     

[45] Oguchi R, Hikosaka K, Hirose T (2003).

Does the photosynthetic light-acclimation need change in leaf anatomy?

Plant, Cell & Environment, 26, 505-512.

[本文引用: 1]     

[46] Ono K, Nishi Y, Watanable A, Terashima I (2001).

Possible mechanisms of adaptive leaf senescence

.Plant Biology, 3, 234-243.

[本文引用: 1]     

[47] Ordoñez JC, van Bodegom PM, Witte JPM, Wright IJ, Reich PB, Aerts R (2009).

A global study of relationships between leaf traits, climate and soil measures of nutrient fertility

. Global Ecology and Biogeography, 18, 137-149.

[本文引用: 3]     

[48] Osnas JLD, Lichstein JW, Reich PB, Pacala SW (2013).

Global leaf trait relationships: Mass, area, and the leaf economics spectrum

.Science, 340, 741-744.

[本文引用: 1]     

[49] Pensa M, Karu H, Luud A, Kund K (2010).

Within-species correlations in leaf traits of three boreal plant species along a latitudinal gradient

.Plant Ecology, 208, 155-166.

[50] Pontes LDS, Soussana JF, Louault F, Andueza D, Carrère P (2007).

Leaf traits affect the above-ground productivity and quality of pasture grasses

.Functional Ecology, 21, 844-853.

[本文引用: 1]     

[51] Reich PB, Oleksyn J (2004).

Global patterns of plant leaf N and P in relation to temperature and latitude

.Proceedings of the National Academy of Sciences of the United States of America, 101, 11001-11006.

[本文引用: 1]     

[52] Reich PB, Wright IJ, Cavender-Bares J, Craine JM, Oleksyn J, Westoby M, Walters MB (2003).

The evolution of plant functional variation: Traits, spectra, and strategies

.International Journal of Plant Sciences, 164(S3), S143-S164.

[本文引用: 1]     

[53] Richter A, Burrows JP, Nüss H, Granier C, Niemeier U (2005).

Increase in tropospheric nitrogen dioxide over China observed from space

.Nature, 437, 129-132.

[本文引用: 1]     

[54] Rose L, Rubarth MC, Hertel D, Leuschner C (2013).

Management alters interspecific leaf trait relationships and trait-based species rankings in permanent meadows

.Journal of Vegetation Science, 24, 239-250.

[本文引用: 4]     

[55] Royer DL (2009).

Nutrient turnover rates in ancient terrestrial ecosystems

.Palaios, 23, 421-423.

[本文引用: 1]     

[56] Royer DL, Miller IM, Peppe DJ, Hickey LJ (2010).

Leaf economic traits from fossils support a weedy habit for early angiosperms

.American Journal of Botany, 97, 438-445.

[本文引用: 1]     

[57] Scoffoni C, Rawls M, McKown A, Cochard H, Sack L (2011).

Decline of leaf hydraulic conductance with dehydration: Relationship to leaf size and venation architecture

.Plant Physiology, 156, 832-843.

[本文引用: 1]     

[58] Shipley B, Lechowicz MJ, Wright I, Reich PB (2006).

Fundamental trade-offs generating the worldwide leaf economics spectrum

. Ecology, 87, 535-541.

[本文引用: 2]     

[59] Silvertown J, Franco FB, Harper JL (1997). Plant Life Histories: Ecology, Phylogeny and Evolution. Cambridge University Press, Cambridge, UK.

[本文引用: 1]     

[60] Suding KN, Collins SL, Gough L, Clark C, Cleland EE, Gross KL, Milchunas DG, Pennings S (2005).

Functional- and abundance-based mechanisms explain diversity loss due to N fertilization

.Proceedings of the National Academy of Sciences of the United States of America, 102, 4387-4392.

[本文引用: 1]     

[61] Sultan SE (2000).

Phenotypic plasticity for plant development, function and life history

.Trends in Plant Science, 5, 537-542.

[本文引用: 1]     

[62] Tamm CO (1990).

Nitrogen in Terrestrial Ecosystems: Questions of Productivity, Vegetational Change, and Ecological Stability

. Springer, Berlin. 116.

[本文引用: 1]     

[63] Ti CP, Yan XY (2010).

Estimation of atmospheric nitrogen wet deposition in China mainland from based on N emission data

.Journal of Agro-Environment Science, 29, 1606-1611. (in Chinese with English abstract)
[遆超普, 颜晓元 (2010).

基于氮排放数据的中国大陆大气氮素湿沉降量估算

. 农业环境科学学报, 29, 1606-1611.]

[64] Ulrich B (1995).

The history and possible causes of forest decline in central Europe, with particular attention to the German situation

.Environmental Reviews, 3, 262-276.

[本文引用: 1]     

[65] Vitousek PM, Aber JD, Howarth RW, Likens GE, Matson PA, Schindler DW, Schlesinger WH, Tilman DG (1997).

Human alteration of the global nitrogen cycle: Sources and consequences

.Ecological Applications, 7, 737-750.

[本文引用: 1]     

[66] Wan HW, Yang Y, Bai SQ, Xu YH, Bai YF (2008).

Variations in leaf functional traits of six species along a nitrogen addition gradient in Leymus chinensis steppe in Inner Mongolia. Journal of Plant Ecology (Chinese Version)

, 32, 611-621. (in Chinese with English abstract)
[万宏伟, 杨阳, 白世勤, 徐云虎, 白永飞 (2008).

羊草草原群落6种植物叶片功能特性对氮素添加的响应

. 植物生态学报, 32, 611-621.]

[67] Wang CS, Wang SP (2015).

A review of research on responses of leaf traits to climate change

.Chinese Journal of Plant Ecology, 39, 206-216. (in Chinese with English abstract)
[王常顺, 汪诗平 (2015).

植物叶片性状对气候变化的响应研究进展

. 植物生态学报, 39, 206-216.]

[本文引用: 1]     

[68] Westoby M, Falster DS, Moles AT, Vesk PA, Wright IJ (2002).

Plant ecological strategies: Some leading dimensions of variation between species

.Annual Review of Ecology and Systematics, 33, 125-159.

[本文引用: 1]     

[69] Wilson PJ, Thompson K, Hodgson JG (1999).

Specific leaf area and leaf dry matter content as alternative predictors of plant strategies

.New Phytologist, 143, 155-162.

[本文引用: 2]     

[70] Wright IJ, Reich PB, Cornelissen JHC, Falster DS, Groom PK, Hikosaka K, Lee W, Lusk CH, Niinemet Ülo, Oleksyn J, Osada N, Poorter H, Warton DI, Westoby M (2005).

Modulation of leaf economic traits and trait relationships by climate

.Global Ecology and Biogeography, 14, 411-421.

[本文引用: 3]     

[71] Wright IJ, Reich PB, Westoby M (2001).

Strategy shifts in leaf physiology, structure and nutrient content between species of high- and low-rainfall and high- and low-nutrient habitats

.Functional Ecology, 15, 423-434.

[72] Wright IJ, Reich PB, Westoby M, Ackerly DD, Baruch Z, Bongers F, Cavender-Bares J, Chapin T, Cornelissen JHC, Diemer M, Flexas J, Garnier E, Groom PK, Gulias J, Hikosaka K, Lamont BB, Lee T, Lee W, Lusk C, Midgley JJ, Navas ML, Niinemets U, Oleksyn J, Osada N, Poorter H, Poot P, Prior L, Pyankov VI, Roumet C, Thomas SC, Tjoelker MG, Veneklaas EJ, Villar R (2004).

The worldwide leaf economics spectrum

.Nature, 428, 821-827.

[本文引用: 1]     

[73] Wright JP, Sutton-Grier A (2012).

Does the leaf economic spectrum hold within local species pools across varying environmental conditions?

Functional Ecology, 26, 1390-1398.

[本文引用: 3]     

[74] Wright RF, Rasmussen L (1998).

Introduction to the NITREX and EXMAN projects

.Forest Ecology and Management, 101, 1-7.

[本文引用: 1]     

[75] Wyka T, Robakowski P, Zytkowiak R (2000).

Leaf age as a factor in anatomical and physiological acclimative responses of Taxus baccata L. needles to contrasting irradiance environments

. Photosynthesis Research, 95, 87-99.

[本文引用: 1]     

[76] Zhang JL, Zhu JJ, Cao KF (2007).

Seasonal variation in photosynthesis in six woody species with different leaf phenology in a valley savanna in southwestern China

. Trees, 21, 631-643.

[本文引用: 1]     

[77] Zhang WN, Liao ZY (2009).

Research advances in the effect of nitrogen deposition on forest plants

. Environmental Science Survey, 28(3), 21-24. (in Chinese)
[张维娜, 廖周瑜 (2009).

氮沉降增加对森林植物影响的研究进展

. 环境科学导刊, 28(3), 21-24.]

[本文引用: 1]     

[78] Zhao XF, Xu HL, Zhang P, Zhang QQ (2014).

Influence of nutrient and water additions on functional traits of Salsola nitraria in desert grassland

.Chinese Journal of Plant Ecology, 38, 134-146. (in Chinese with English abstract)
[赵新风徐海量, 张鹏, 张青青 (2014).

养分与水分添加对荒漠草地植物钠猪毛菜功能性状的影响

. 植物生态学报, 38, 134-146.]

[79] Zheng XH, Fu CB, Xu XK, Yan XD, Huang Y, Han SH, Hu F, Chen GX (2002).

The Asian nitrogen cycle case study

.AMBIO, 31(2), 79-87.

[本文引用: 1]     

/