刘旻霞

西北师范大学地理与环境科学学院,兰州 730070

植物功能性状是植物响应环境变化并对生态系统功能有一定影响的植物结构和生理特性,这些特性对生态系统结构特征具有一定的响应[1]。功能性状通过影响个体的生存、生长、生殖来间接地影响其适合度,可以作为预测有机体生长表现的重要指标[2]。物种间在功能性状间的差异既反映了其所处的生境特征,也体现了它们不同的资源利用方式。比如生长在水分不足生境中的植物一般具有较小的比叶面积(SLA)和较低的叶片含水量(LWC),体现了在控制水分散失和提高生长速率之间的权衡[3]。

叶片是植物与外界环境进行物质与能量交换的主要器官,是植物获取能量的最重要的组织结构,因此叶片水平功能性状的研究是植物生态学长期以来研究的重点之一[4- 5]。植物叶片的光合生理和形态特征(例如SLA、干物质量(LDMC)、净光合速率(Pn)、13C、相对叶绿素(SPAD)等)体现了植物的生长策略和资源利用方式,是植物与环境长期相互作用的结果,也是植物完成生活史周期和适应环境的重要生理生态特性[6- 7]。植物的光合生理和叶片形态特征除了受到诸多自然环境因子的影响外[8],在不同物种之间也有很大差异[7],比如杂草类植物一般具有较大的SLA和较低的水分利用效率(WUE),而禾草类植物则具有较高的WUE等。因此,研究同一地域不同植物物种光合和叶片形态特征之间的差异,对于理解和预测整个群落的资源利用及群落结构动态具有重要的意义[9]。

高寒草甸是青藏高原重要的草地类型,在这个特殊的地理单元上具有其独特的气候特征,被认为是气候变化的敏感区域。而坡向是青藏高原地区的主要地形因子之一,坡向影响了地面与风向的夹角和地面接受的太阳光辐射,这使得不同坡向之间的光照度(LI)、土壤温度(ST)、土壤含水量(SWC)、养分和植被分布均受到影响,是研究群落在微气候生境的发展、群落结构和功能变化,尤其是研究植物功能性状对生态系统功能响应的良好载体。研究认为,从北坡到南坡的生境梯度上,SWC在不断减少,ST和LI在不断增加,在此过程中,植物为适应减少的水分供给,其在功能性状方面必然也会表现出一定的响应[10]。在此情况下,我们提出如下科学假说：在甘南亚高寒草甸中,植物功能性状对不同坡向SWC、ST、LI及养分的响应具有一定的规律性,则可以推测：在北坡-南坡梯度上,随着环境因子的变化,植物的光合生理和叶片形态等也会发生相应变化。另外,植物需要从土壤中吸收水分和养分等用于光合作用和蒸腾作用,所以,叶片的SLA、LDMC、Pn、13C、SPAD值等与土壤因子存在直接关系。一般来说,随着坡向SWC的增加,ST和LI的减少,Pn、气孔导度等随坡向有降低的趋势[10]。因此,叶片的SLA、LDMC、Pn、13C、SPAD与这些土壤因子间理应存在显著的关系。为了验证以上科学假说,本研究选择在甘南亚高寒草甸的山地坡向间,通过研究植物功能性状和SWC、土壤养分、ST及LI的规律以及它们之间的相关关系,分析植物SLA、LDMC、13C值、Pn、SPAD及LWC等植物功能性状是否与其环境因子间具有显著关系。解决上述科学问题,将有助于揭示微气候生境植物适应环境因子变化的功能策略。

1 材料与方法

1.1 研究区概况

实验地位于甘肃省甘南藏族自治州合作市附近的兰州大学高寒草甸与湿地生态系统定位研究站,地理位置为N34°56′,E102°54′,海拔3010m。年均气温2.1℃,年均降水量558.9mm。最冷的12月至次年2月平均气温-8.8℃,最热的6月—8月平均气温为11.6℃,≥0℃的年积温为1732℃左右,南坡和西南坡主要以矮嵩草(Kobresiahumilis)、三刺草(Aristidatriseta)、及蒲公英(Taraxacumlugubre)等为主；西坡主要以黄花棘豆(Oxytropisochrocephala)、米口袋(Gueldenstaedtiaverna)等为优势种；西北坡和北坡主要以金露梅(Potentillafruticosa)灌丛及珠芽蓼(Polygonumvivipurum)、鹅绒萎陵菜(Potentillaanserine)等杂类草为主要优势种。土壤为高寒草甸土。

1.2 实验设计

实验于2015年7月上旬到8月上旬进行,在合作站附近选择一座南北坡分异明显的山地,测定坡向使用360°罗盘,在海拔3000m的山体中部(坡度见表1),采用顺时针方向布设5个研究样地,样地间距约25—35m,分别标记N(北坡)、NW(西北坡)、W(西坡)、SW(西南坡)及S(南坡)。在每个样地中部,沿山体垂直方向布设4个调查样方,大小均为50cm×50cm,间隔1m,调查每个样方的物种数、个体数及盖度。同时,用直径5cm规格的土钻在每个样方内采用对角线方法分别钻取5钻土壤(0—20cm)及(20—40cm),混合装入铝盒,用来做后续的水分分析。其样地布设见图1。

表1 不同坡向植被变化特征(平均值±标准误)

*为各坡向优势种

图1 研究区样地布置图Fig.1 Sample plot of the study area N:北坡North slope; NW: 西北坡North west slope; W: 西坡West slope; SW:西南坡West south slope; S:南坡South slope

1.2.1 环境因子测定

用烘干法重复3次测定土壤含水量,并用此公式 SWC=(massf-massd)/massf进行计算。测量土壤温度用EM50(Decagon,USA),并重复3次。测量各样地光照度用Sunfleck Ceptometer(型号ZDS- 10)进行测量,并与地面保持水平,测量高度距离地面60cm；设置测量时间从6:00am—20:00pm,每1个小时测量1次。在2015年的 6—9月份,在每个样地相隔20天左右测量土壤含水量、土温和光照度共4次。

1.2.2 土壤养分的测定

土壤全氮采用硫酸-K2SO4∶CuSO4∶Se催化法消煮,然后用SmartChem 200 化学分析仪器(WestCo Scientific Instruments, Brookfield, CT, USA)进行测试。土壤全磷用H2SO4-HClO4消煮后,采用钼锑抗法显色法测定。土壤有机碳采用重铬酸钾容量法。

1.2.3 物种多样性指数及重要值的计算公式

丰富度指数R=S

物种多样性指数的计算方法采用Shannon-Wiener指数计算法

S为i所在样方中物种的总数。

为了确定各个坡向的优势种,通过计算其重要值来判定。重要值计算方法如下：

重要值=(相对多度+相对盖度+相对频度)/3

1.2.4 植物叶片性状的测定

2015年7月上旬到8月上旬,选择晴朗天气(大气温度为25℃,光照强度为92K),使用便携式光合测定仪(LI- 6400,LI-COR,Lincoln,NE,USA),于上午9:00—11:30,测定植物叶片的气体交换参数。测定时选用开放气路,空气流速为500μmol/s,用红蓝光源提供1000μmol m-2s-1光强。随机选取长势良好的物种3—6株,取中端健康叶片进行测量,并做3—4次重复,数据稳定后读数。用便携式叶绿素仪(SPAD- 502,Minolta Camera Co,Osaka,Japan)测定植物叶片SPAD值,所测物种及重复数与测定光合特征值时的一致。SPAD- 502是在比较650nm和940nm波长的光线投射率的基础上测定的[11],因此所测叶绿素含量是一个相对值。

在各样地内摘取物种的健康叶片,标记好并装入自封袋带回实验室,扫描叶片保存图像并称其鲜重及饱和鲜重,之后烘干称重,计算叶面积用ImageJ软件。比叶面积(SLA)表示为叶面积/叶片干重；叶片干物质量(LDMC)用叶片干重/叶片饱和鲜重表示；叶片相对含水量表示为(叶片鲜重-叶片干重)/(叶片饱和鲜重-叶片干重)×100%。

1.2.5植物叶片13C值的测定

在实验室内把摘回的植物叶片用超纯水洗净并自然晾干,70℃烘干48h,使样品完全干燥,粉碎后,过100目筛,使样品均一化。取植物粉末2mg置入TOC/TNb分析仪中将样品燃烧成CO2,并收集。再将收集到的CO2气体注入碳同位素分析仪(CCIA- 36d)进行稳定碳同位素值分析[12]。分析精度± 0.1‰。其中,稳定碳同位素(13C)的计算公式可表示为：

1.3 数据分析

采用Excel2007统计计算数据,用SPSS18.0软件进行one-way ANOVA差异性检验和相关性分析,并用sigmplot10.0作图,RDA分析应用国际标准通用软件CANOCO4.5.1进行。

2 结果与分析

2.1 SWC、ST及LI随坡向的变化

从北坡到南坡,环境因子均发生了变化,其中SWC在3个月的4次测量中差异最明显,从测量的平均值来分析,不同坡向之间均有显著差异(P<0.05),表现是从南坡到北坡的递增趋势；ST和LI则表现出了相反的变化趋势,从南坡到北坡呈递减趋势,7月份的土温高于其他两个月,各测量月平均值在不同坡向之间有显著差异(P<0.05)；LI北坡只有580lx,而南坡最高为730lx,从不同测量月份来看,8、9月份的LI在不同坡向均低于7月份的。各月测量平均值南坡与北坡之间有显著差异(P<0.05),其他坡向差异不大(图2)。

图2 土壤水分(SWC)、温度(ST)及光照度(LI)随坡向的变化(Mean±SE)Fig.2 Change in SWC, ST and LI on different slope aspects (Mean±SE) N: 北坡；NW: 西北坡；W: 西坡；SW: 西南；S: 南坡

2.2 土壤养分在不同坡向的变化

0—20cm的土层中,土壤全氮(STN)和土壤有机碳含量(SOC)都是北坡>南坡,西坡相对较低,STN北坡、西北坡、西南坡、南坡与西坡之间有显著差异(P<0.05)。SOC是北坡、西北坡和西坡、西南坡、南坡之间均有显著差异(P<0.05)。而土壤全磷含量(STP)的变化则不同,自南坡、西南坡、西坡、西北坡到北坡依次增大,对不同坡向STP进行差异显著性检验,0—20cm STP在不同坡向均表现出显著差异(P<0.05)。20—40cm的STN、SOC、STP的变化趋势与0—20cm一致,但总体都低于0—20cm土壤,表现出表层土壤养分富集现象。20—40cm STN的变化与0—20cm不一致,趋势为北坡>西北坡>南坡>西坡>西南坡,但其值在不同坡向低于表层土壤(图3)。

图3 不同坡向土壤养分的变化特点(Mean±SE)Fig.3 Variation of the soil nutrients in different slope aspects (Mean±SE)N: 北坡；NW: 西北坡；W: 西坡；SW: 西南；S: 南坡; 同一部分具有不同字母的表示差异显著(P<0.05)

2.3 不同坡向植物群落结构的变化

从表1可以看出,优势种及群落结构在坡向梯度上发生了显著的变化：南坡、西南坡的物种分布较少,主要物种有莎草科的矮嵩草(Kobresiahumilis)与菊科的蒲公英(Taraxacumlugubre)等；西坡的豆科类物种占据比较大的优势；而北坡、西北坡则主要以金露梅(Potentillafruticosa)灌丛为主,并伴生一些杂类草(即非豆科杂草)。物种丰富度、多样性指数及盖度自北坡到南坡呈递减趋势,而且南坡的坡度(32.08°)明显大于北坡(20.02°)。

2.4 光合生理和叶片形态特征的生境差异

由图4可知,北坡-南坡梯度上,4种代表物种的6个性状均发生了不同程度的变化,其中矮嵩草在西南坡具有最大的LDMC(0.38),北坡、西北坡与西坡、西南坡及南坡之间有显著差异(P<0.05)。蒲公英在西北坡具有最小的LDMC(0.17),北坡、西坡、西南坡与西北坡及南坡之间有显著差异(P<0.05)。米口袋的LDMC北坡、西北坡、西坡与西南坡、南坡之间有显著差异(P<0.05)。金露梅北坡与西北坡及西坡之间有显著差异(P<0.05)。Pn的最大值出现在南坡的矮嵩草(18.69μmolCO2m-2s-1),而最小值是西坡的金露梅(8.01μmolCO2m-2s-1),总体来看,物种在北坡-南坡梯度上Pn不断增大。SLA的最大值是北坡的蒲公英(252.73cm2/g),最小值是南坡的矮嵩草(78.26cm2/g)。在北坡-南坡梯度上,SLA呈减小趋势。4个共有种的叶片13C值在坡向梯度上均有显著差异(P<0.05),其中,矮嵩草的13C值在南坡最大为(-23.02‰),金露梅的13C值在北坡最小为(-31.01‰)。LWC从北坡-南坡呈逐渐减少趋势,其中蒲公英的LWC在北坡最大为(7.05%),最小值是西坡的矮嵩草(2.21%)。北坡-南坡梯度上,4个物种的SPAD值不断增大,其中米口袋在北坡具有最大的SPAD值(43.41),最小值是西北坡的金露梅(23.04)。对同一物种在不同坡向的性状特征值进行比较发现,LDMC、Pn、13C及SPAD值在北坡-南坡梯度上呈递增趋势,而SLA、LWC呈递减趋势。

图4 不同坡向物种在叶干物质量(LDMC)、净光合速率(Pn)、比叶面积(SLA)、叶片同位素(13C)、叶片相对含水量(LWC)及相对叶绿素(SPAD) 的差异(平均值±标准误差)Fig.4 Comparison of species for LDMC, net photosynthesis rate (Pn), specific leaf area (SLA), 13C, relative leaf water content (LWC) and SPAD in difference slope aspects (P<0.05). 不同小写字母表示同一物种在不同坡向间差异显著(P< 0.05);AS：矮嵩草Kobresia humilis；MKD：米口袋Gueldenstaedtia verna；PGY：蒲公英Taraxacum lugubre；JLM：金露梅Potentilla fruticosa

2.5 植物叶性状与环境因子的关系

本研究采用冗余分析(RDA)方法,按5个坡向统一处理数据进行分析,将植物叶片的6个主要性状作为物种(species),将6项土壤因子作为环境因素(environments),并将2个变量组进行RDA分析(图5)。图中实线箭头所处的象限表示环境因子与排序轴之间的正负相关性,箭头连线的长度代表着某个环境因子与研究对象相关程度的大小,连线越长,影响越大。箭头与箭头之间的夹角越小,相关性越大,反之越小。实线箭头与排序轴的夹角代表着某个环境因子与排序轴的相关性大小。图5显示,矮嵩草、米口袋、蒲公英及金露梅的SLA、LWC总体上与SWC、STN、STP及SOC呈正相关,与ST、LI呈负相关,性状指标SPAD、Pn、13C、LDMC则与SWC、STN、STP及SOC呈负相关关系,与ST、LI呈正相关相关系,且相关性大小表现为：SWC>STP >ST>SOC >LI> STN。表2为各物种RDA第一、二排序轴与环境及土壤因子的相关系数,由表2可以看出,环境因子对植物叶片性状的前2轴的累积解释量占特征值的总和均达 99% 以上,表明排序轴可有效反映甘南亚高寒草甸微气候生境上土壤及其他环境因子的梯度变化特征。

表2 各物种RDA排序轴与环境因子的相关系数

3 结论与讨论

3.1 南北坡梯度上环境因子变化特征

光照强度容易受地形因素的影响,而土壤温度、水分等又会受到光照强度的影响。本研究中,从北坡-南坡的坡向梯度上,环境因子发生了很大变化,土壤含水量在4次测量中都呈现出了明显的下降趋势,南坡明显低于北坡,而土壤温度和光照度则呈明显的上升趋势(图2)。在北半球,太阳辐射较多的是南坡,日照时数长,蒸发量大,因此热而干旱；而北坡受到了较少的太阳辐射,日照时数短,蒸发量小,因此比较潮湿、阴冷[13]。这些因素使得南坡矿化作用较强,养分含量低,比较贫瘠；而北坡有机质降解较低、累积较多、土壤腐殖化较强。本研究显示：在北坡-南坡梯度上,土壤全氮、全磷及有机碳含量总体是北坡较高,南坡较低。有研究表明,湿度越大,土壤中的氮素越容易积累；温度越高,土壤有机质分解加速,氮素越难积累[14]。以前的研究也表明水分是这个梯度上的主要限制因子[15]。西坡的土壤全氮及有机碳含量比南坡更低,这可能是由于西坡处于放牧通道,动物的踩踏行为导致部分土壤养分含量减少。总的来说,一个地区土壤环境因子不但受水热因子影响,同时也受到土壤特性、放牧干扰和群落类型等的综合作用。

3.2 南北坡梯度上植物群落结构变化特征

由于不同坡向的光照、热量、水分及养分等一些自然因素的不同,植被的空间分布格局也不同,同时有一定的规律性[16]。另外,坡向梯度上的土温、土壤含水量、物种分布等和经纬度及海拔等大尺度上的变化一致[15]。在北坡-南坡梯度上,物种的组成类型发生了明显的梯度变化(表1)。本研究结果显示禾草类(包括矮嵩草和三刺草)是南坡的优势种,这说明禾草类具有比较高的营养利用率[17],能够在南坡贫瘠的生境中,凭借着对地下养分资源较强的竞争能力而占据优势。而灌木(金露梅)和杂类草则是北坡的优势种,这表明在北坡充足的资源条件下,杂类草和灌木快速地获取资源而生长迅速,在竞争中占据优势。与禾草和杂类草及灌木相比,豆科物种占有相对较小的数量,在西坡分布相对较多(表1)。沿着北坡-南坡的微生境格局,物种的Shannon-Wiener指数、Richness指数和盖度都显著的降低了(表1)。有研究表明：在干旱到湿润的生境梯度上,物种多样性指数和丰富度指数不断增加[18]。这与本研究结果一致,这是因为物种多样性和丰富度一般被认为是由资源供给所决定的[19]。这表明在高寒草甸地区,坡向梯度上物种的多样性变化是由可利用资源所决定的,北坡丰富的资源能够使得更多的物种存活下来。

图5 植物叶性状与环境因子的RDA二维排序图Fig.5 RDA two-dimensional ordination diagram of the plant leaf traits and the environmental factorsa:矮嵩草 b:米口袋 c:蒲公英 d:金露梅; N:北坡, NW:西北坡, West:西坡, SW:西南, South:南坡;SWC:土壤含水量; ST:土壤温度; LI:光照度; STP:土壤全磷; SOC:土壤有机碳; STN:土壤全氮; LWC:叶片相对含水量; SLA:比叶面积; LDMC:叶片干物质量; Pn:净光合速率; SPAD:相对叶绿素; C13:13C值

3.3 坡向梯度上物种性状变化特征

许多研究表明,植物在叶片形态特征(如SLA、LWC等)和光合生理特征(如Pn、13C等)上存在着显著的种间和种内差异[20],即使它们生长在非常相似的环境中[21],被认为是长期自然选择的结果。我们对甘南亚高寒草甸4个主要物种的叶片形态和光合生理特征的分析表明,不同物种之间在LWC、SLA、LDMC、Pn、SPAD和13C之间有非常明显的差异,说明这些物种虽然生长在非常相似的环境中,但是在叶片形态和光合策略上有着显著的不同[20],以不同的方式适应相似的生境条件。

叶片形态特征的中心性状SLA能充分反映温度、水分、光照等气候特点。本研究中的4种物种的SLA分布范围是在78.26—252.73cm2/g之间。由此可见,相对于北坡的物种,亚高寒草甸南坡的物种叶片结构趋向于单位叶面积分配更多的生物量,表现出低水分和营养元素地带植物的生长策略特点,其具有阳生叶特点(低SLA,高LDMC)的叶片更有利于对光照的利用和自身的生存。LDMC可反映叶片对干旱气候适应程度,LDMC越高则叶密度越大,细胞壁越厚,具有更强的抵抗胁迫的能力[22]。本研究中的4种物种从北坡-南坡的梯度上,LDMC都有不同程度的增大,其中南坡的矮嵩草LDMC最大,表明它抵抗干旱胁迫的能力更强。LWC随北坡南坡变化而减少,说明随土壤含水量减少、土壤温度和光照度增加,植物叶片含水量也随之减少(图4)。一般而言,土壤水分越充足,植物叶片含水量也就越高。说明植物叶片含水量随着土壤水分条件的改善而增加。

光合作用决定着植物物质积累能力的高低,一定程度上也决定了植物生长的快慢。北坡-南坡梯度上,4个物种的Pn和SPAD逐渐增大(图4),与前人的研究一致[21]。一般而言,阳生植物的净光合速率高于阴生植物,较高光合速率的物种也被认为在能量获取和竞争力上有优势[23]。本研究中,4个物种在北坡的光合能力较弱(图4),可能是养分等限制的原因[24- 26]。在一个土壤水分、养分逐渐升高北坡(图4),养分和水分对植物的限制作用逐渐减小,4个物种的Pn和SPAD值的降低,反映了不同物种通过对营养物质的不同分配策略来保护匮乏资源。植物叶片13C值是反应植物水分利用效率最直接的指标,其含量的多少不仅与植物本身有关,还与植物生长的环境有密切关系,其中最主要的就是土壤含水量[18]。本文研究结果显示,南坡的优势物种矮嵩草叶13C含量明显高于北坡代表物种金露梅(图 4),表明矮嵩草比金露梅具有较高的水分利用效率,更能适应南坡的贫瘠生境。

3.4 坡向梯度上物种性状与环境因子的相关关系

植物叶性状值能清楚的反映环境变化对植物的影响。冗余分析进一步研究了坡向梯度上的环境因子和各物种性状之间的相关关系,探索了环境因子对物种组成及性状的影响大小。6个环境因子均对群落物种分布及性状变化起了作用(图5)。分析结果表明,环境因子对物种组成及性状影响顺序为：SWC> STP> ST> SOC> LI> STN。表明土壤含水量的影响力最大,是坡向梯度上群落结构变化、物种组成及性状值变化的主导因子,这一结果与其他研究结果一致,即SWC是坡向梯度上物种组成及多样性的主要限制性因素[18]。坡向梯度上物种的SLA及LWC与第一轴正相关,与土壤温度及光照度呈负相关,这进一步说明了物种的SLA与LWC与可利用资源有关,北坡丰富的水分资源和养分资源有利于灌木金露梅和杂草类生存。南坡由于光照及土壤温度导致的蒸发量较大,且坡度较大(均值32.08°),土壤含水量相对北坡较少,因此南坡物种的叶片13C、Pn、SPAD及LDMC相对较高,反映了植物为适应南坡贫瘠的生境,其功能性状具有明显的响应特征。

[1] Cavender-Bares J, Kitajima K, Bazzaz F A. Multiple trait associations in relation to habitat differentiation among 17 Floridian oak species. Ecological Monographs, 2004, 74(4): 635- 662.

[2] Violle C, Navas ML, Vile D, Kazakou E, Fortunel C, HummelI, Garnier E. Let the concept of trait be functional! Oikos, 2007, 116(5): 882- 892.

[3] Nicotra A B, Hermes J P, Jones C S, Schlichting C D. Geographic variation and plasticity to water and nutrients inPelargoniumaustrale. New Phytologist, 2007, 176(1): 136- 149.

[4] Smith EC, Griffiths H, Wood L, Gillon J. Intra-specific variation in the photosynthetic responses of cyanobiont lichens from contrasting habitats. New Phytologist, 1998, 138(2): 213- 224.

[5] 王海珍, 韩路, 徐雅丽, 牛建龙, 于军. 干旱胁迫下胡杨光合光响应过程模拟与模型比较. 生态学报, 2017, doi: 10.5846/stxb201511242373. (in press)

[6] Cunningham S A, Summerhayes B, Westoby M. Evolutionary divergences in leaf structure and chemistry, comparing rainfall and soil nutrient gradients. Ecological Monographs, 1999, 69(4): 569- 588.

[7] 牛书丽, 蒋高明, 高雷明, 李永庚, 刘美珍. 内蒙古浑善达克沙地97种植物的光合生理特征. 植物生态学报, 2003, 27(3): 318- 324.

[8] Lawlor D W. Photosynthesis, productivity and environment. Journal of Experimental Botany, 1995, 46: 1449- 1461.

[9] Wang R Z. Photosynthesis, transpiration, and water use efficiency of vegetative and reproductive shoots of grassland species from north-eastern China. Photosynthetica, 2001, 39(4): 569- 573.

[10] 刘旻霞, 马建祖. 甘南高寒草甸植物功能性状和土壤因子对坡向的响应. 应用生态学报, 2012, 23(12): 3295- 3300.

[11] Manetas Y, Grammatikopoulos G, Kyparissis A. The use of the portable, non-destructive, SPAD- 502 (Minolta) chlorophyll meter with leaves of varying trichome density and anthocyanin content. Journal of Plant Physiology, 1998, 153(3/4): 513- 516.

[12] Swap R J, Aranibar J N, Dowty P R, Gilhooly IIIWP, MackoSA. Natural abundance of13C and15N in C3and C4vegetation of southern Africa: patterns and implications. Global Change Biology, 2004, 10(3): 350- 358.

[13] Sternberg M, Shoshany M. Influence of slope aspect on Mediterranean woody formations: comparison of a semiarid and an arid site in Israel. Ecological Research, 2001, 16(2): 335- 345.

[14] 林丽, 李以康, 张法伟, 郭小伟, 韩道瑞, 李婧, 曹广民. 青藏高原高寒矮嵩草草甸退化演替主成分分析. 中国草地学报, 2012, 34(1): 24- 30.

[15] Badano E I, Cavieres L A, Molina-Montenegro M A, Quiroz C L. Slope aspect influences plant association patterns in the Mediterranean matorral of central Chile. Journal of Arid Environments, 2005, 62(1): 93- 108.

[16] 方精云, 沈泽昊, 崔海亭. 试论山地的生态特征及山地生态学的研究内容. 生物多样性, 2004, 12(1): 10- 19.

[17] Kull O, Aan A. The relative share of graminoid and forb life-forms in a natural gradient of herb layer productivity. Ecography, 1997, 20(2): 146- 154.

[18] De Bello F, Lepš J, Sebastià M T. Variations in species and functional plant diversity along climatic and grazing gradients. Ecography, 2006, 29(6): 801- 810.

[19] Morin P J. Biodiversity’s ups and downs. Nature, 2000, 406(6795): 463- 464.

[20] Yu Q, Chen Q, Elser J J. Linking stoichiometric homoeostasis with ecosystem structure, functioning and stability. Ecology letters, 2010, 13(11): 1390- 1399.

[21] Bassow S L, Bazzaz F A. Intra-and inter-specific variation in canopy photosynthesis in a mixed deciduous forest. Oecologia, 1997, 109(4): 507- 515.

[22] 毛伟, 李玉霖, 赵学勇, 黄迎新, 罗亚勇, 赵玮. 3种藜科植物叶特性因子对土壤养分、水分及种群密度的响应. 中国沙漠, 2009, 29(3): 468- 473.

[23] Givnish T J. Adaptation to sun and shade: a whole-plant perspective. Australian Journal of Plant Physiology, 1988, 15(2): 63- 92.

[24] 安慧, 上官周平. 黄土高原植被不同演替阶段优势种的光合生理特性. 应用生态学报, 2007, 18(6): 1175- 1180.

[25] 代存芳, 易映彤, 刘妍, 巴秋爽, 范亚文. 扎龙湿地硅藻植物群落季节变化及其对环境的响应. 生态学报, 2017, doi: 10.5846/stxb201603100424. (in press)

[26] 孟令军, 张利敏, 张丽梅, 冯仲科. 五大连池药泉山大型土壤动物对旅游踩踏的响应. 生态学报, 2016, 36(20): 6607- 6617.