植物生态学报  2016 , 40 (8): 788-797 https://doi.org/10.17521/cjpe.2016.0003

Orginal Article

张掖湿地芨芨草叶大小和叶脉密度的权衡关系

韩玲, 赵成章*, 徐婷, 冯威, 段贝贝, 郑慧玲

西北师范大学地理与环境科学学院, 甘肃省湿地资源保护与产业发展工程研究中心, 兰州 730070

Trade-off between leaf size and vein density of Achnatherum splendens in Zhangye wetland

HAN Ling, ZHAO Cheng-Zhang*, XU Ting, FENG Wei, DUAN Bei-Bei, ZHENG Hui-Ling

College of Geography and Environmental Science, Northwest Normal University, Research Center of Wetland Resources Protection and Industrial Development Engineering of Gansu Province, Lanzhou 730070, China

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

基金资助:  国家自然科学基金(41461013和91125014)

展开

摘要

叶大小-叶脉密度的权衡关系是植物叶经济谱理论的基础, 对理解资源竞争条件下植物叶片的物理构建与生理代谢的关系具有重要的意义。该文采用标准化主轴估计(standardized major axis estimation, SMA)的方法, 按芨芨草(Achnatherum splendens)株丛密度设置I (>12丛·m-2)、II (8-12丛·m-2)、III (4-8丛·m-2)和IV (<4丛·m-2) 4个密度梯度, 以叶面积和叶干质量分别表示叶大小, 对张掖洪泛平原湿地不同密度条件下芨芨草种群的叶大小和叶脉密度的关系进行研究。结果表明: 随着芨芨草株丛密度的降低, 湿地群落的土壤含水量逐渐减小、土壤电导率逐渐增加, 芨芨草的净光合速率(Pn)、蒸腾速率(Tr)和分枝数呈先增大后减小的趋势, 叶面积、叶干质量、比叶面积和株高呈逐渐减小趋势、光合有效辐射(PAR)和叶脉密度呈逐渐增加趋势; 芨芨草叶大小和叶脉密度在高密度(I)和低密度(IV)样地均呈极显著负相关关系(p < 0.01), 中密度(II、III)样地二者呈显著负相关关系(p < 0.05); 叶大小和叶脉密度回归方程的SMA斜率在不同密度样地均显著小于-1 (p < 0.05), 即芨芨草叶大小和叶脉密度呈“此消彼长”的权衡关系。在高密度湿地群落芨芨草倾向于大叶片低叶脉密度的叶片构建模式, 在低密度湿地群落选择小叶片高叶脉密度的异速生长模式, 体现了密度制约下湿地植物的生物量分配格局和资源利用对策。

关键词: 叶面积 ; 叶干质量 ; 叶脉密度 ; 权衡 ; 芨芨草 ; 张掖湿地

Abstract

Aims Trade-offs between leaf size and vein density are the basis of the theory of leaf economics spectrum, and are to understand the relationship between the physical build and physiological metabolism of plant leaves under different degrees of competition for resources. Our objective was to study the changes in the relationship between leaf size and vein density (leaf dry biomass and leaf area) in Achnatherum splendens populations with four plant bundle densities located in the flood plain wetland of Zhangye. Methods The study site was located at floodplain wetlands of Zhangye, Gansu Province, China. Survey and sampling were carried out in the communities that A. splendens dominated. According to the plant bundle density, the A. splendens communities were divided into four density gradients with “bundle” for the sampling units, high density (I, > 12 bundle·m-2), medium density (II, 8-12 bundle·m-2), medium density (III, 4-8 bundle·m-2) and Low density (IV, <4 bundle·m-2). According to the density of each combination, we chose seven (5 m × 5 m) A. splendens samples, resulting in a total of 28 samples (4 × 7). The soil physical and chemical properties of four density gradients were investigated and six samples of A. splendens were used to measure the leaf area, leaf dry biomass and vein density in laboratory, and biomass of different organs was measured after being dried at 85 °C in an oven. 28 plots were categorized into three groups: high, medium and low density, and the standardized major axis (SMA) estimation method was used to examine the allometric relationships between leaf area, leaf dry biomass and vein density. Important findings The results showed that with the population density changed from high, medium, to low, the soil moisture decreased, and soil electric conductivityincreased. The leaf area, leaf biomass and height of A. splendens decreased, and the vein density, specific leaf area and photosynthetically active radiation (PAR) increased gradually. In addition, leaf net photosynthetic rate (Pn), transpiration rate (Tr) and twig number firstly increased then decreased. There was a highly significantly negative correlation (p < 0.01) between the leaf size and vein density on the high- and low-level densities (I, IV), whereas less significant (p < 0.05) on the level of medium density (II, III). The SMA slope of regression equation in the scaling relationships between leaf size and vein density was significantly smaller than -1 (p < 0.05).

Keywords: leaf area ; leaf dry biomass ; vein density ; trade-off ; Achnatherum splendens ; Zhangye wetland

0

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

本文引用格式 导出 EndNote Ris Bibtex

韩玲, 赵成章, 徐婷, 冯威, 段贝贝, 郑慧玲. 张掖湿地芨芨草叶大小和叶脉密度的权衡关系. 植物生态学报, 2016, 40(8): 788-797 https://doi.org/10.17521/cjpe.2016.0003

HAN Ling, ZHAO Cheng-Zhang, XU Ting, FENG Wei, DUAN Bei-Bei, ZHENG Hui-Ling. Trade-off between leaf size and vein density of Achnatherum splendens in Zhangye wetland. Chinese Journal of Plant Ecology, 2016, 40(8): 788-797 https://doi.org/10.17521/cjpe.2016.0003

植物是一个功能和结构平衡的生命体, 在漫长的进化过程中通过逐步调整结构, 使其功能和结构达到权衡状态, 这对植物提高光合收益和优化生存策略具有重要的意义(Brodribb & Jordan, 2011; Reich & Cornelissen, 2014)。叶片形态及叶脉网络功能性状能够反映植物适应特定生境的基本行为和光合生理功能, 体现了植物提高叶片光截取、碳获取以及水分输导效率的生态策略(金鹰和王传宽, 2015)。叶脉为广泛分布于整个叶片的网络系统结构(Brodribb et al., 2010), 叶脉密度反映了叶脉与叶肉的接触程度和叶片的碳投资模式, 其大小影响光合蒸腾过程中水分、养分及光合产物运输的效率与安全性(Sack & Frole, 2006; Sack & Scoffoni, 2013); 叶大小的变化与某一给定的叶生物量分配和总叶面积密切相关(Niklas et al., 2007), 决定了植物对光的拦截效率和碳获取能力(王常顺和汪诗平, 2015), 直接影响着植物光合生产力和生存适合度(Castro-Díez et al., 2000); 比叶面积是一定干物质投资所展开的捕光表面积数量, 对植物的相对生长速率具有重要影响(Wright et al., 2004)。叶片结构性状和光合生理活性间相互联系、相互依存的内在关系, 反映了植物的结构特性和生物学过程对资源利用水平的决定作用(金鹰和王传宽, 2015)。随着叶片的展开, 植物的净光合速率(Pn)和叶片蒸腾速率(Tr)会随之增加(Sack et al., 2012), 较高的叶脉密度能够满足蒸腾带来的大量水分和能量消耗, 有助于植物调节叶片温度和维持体内水分平衡(Nardini et al., 2010), 有限的光合产物不可能同时兼顾叶肉、叶脉等组织间的资源投入, 导致在某种情况下叶脉密度又会限制叶面积和光合组织密度最大化, 反映了叶片蒸腾成本与光合收益之间的相互匹配与经济权衡关系(Sack & Scoffoni, 2013), 对植株构建最佳的生理和功能状态发挥着重要的调控作用。因此, 叶大小与叶脉密度之间相互影响、相互制约的生态学机制是理解植物生理活动中叶片功能多样化和光合碳固定的重要依据(杨冬梅等, 2012; 李乐等, 2013)。密度是自然界植物种群普遍存在的一种重要选择压力 (Japhet et al., 2009), 随着种群密度的增加, 植物生存所需空间、光照和土壤水分等资源可能受到邻株的物理或几何上的拥挤效应, 从而引起不同程度的隐蔽度和竞争强度的干扰(张明娟等, 2012), 改变了种群内每个植株可获得性资源的数量, 植物的光合作用将被抑制(史元春等, 2015), 促使植物权衡叶大小和叶脉密度的生物量分配、调整叶片构型构建, 实现叶片生理特性和环境空间异质性的精细协调, 体现了最大化地吸取异质性分布资源的觅食行为和可塑性反应, 也是植物最经济的生境适应性和风险规避对策(Silvertown & Charlesworth, 2001)。因此, 研究不同密度间植物叶大小和叶脉密度的关系及成因, 对于理解干旱区湿地群落密度制约下植物叶经济谱的形成机制具有重要意义。

芨芨草(Achnatherum splendens)是我国西北地区常见的多年生密丛性草本植物, 草秆坚硬、直立丛生、根系发达、须根较多, 叶片拥有较强的机械组织和输导组织, 具有较强的耐旱、耐盐碱的表型可塑性, 通常发育于轻度盐化的河湖边缘以及低洼地等隐域性生境, 是一种优良的饲用植物、纤维植物和水土保持植物(张翼飞等, 2012)。近年来, 学术界对植物叶大小的形成机制和环境适应性(Mcdonald et al., 2003; 刘明虎等, 2013; 党晶晶等, 2014)、叶脉网络功能性状与叶片光合(Brodribb et al., 2010; Sack et al., 2012)、水分利用能力(Zhang et al., 2012; Sack & Scoffoni, 2013)、碳构建(Reich & Cornelissen, 2014)的关系以及叶脉密度与叶大小的比例关系(sack et al., 2012)等已有大量研究, 而叶片各功能性状间关联性的形成可能依赖于不同研究尺度和环境因子的差异。目前对芨芨草的研究集中在不同密度芨芨草空间格局对环境胁迫的响应(张明娟等, 2012)、生物学特征和抗盐性(张雅琼等, 2010)及水分利用策略(Wang et al., 2004)等方面, 但芨芨草叶片形态和叶脉性状的生长规律研究不深入, 尤其湿地群落不同密度生境叶大小和叶脉密度的生长关系对植物生理功能和异质性资源配置的响应机理尚不清晰。鉴于此, 本文用叶面积和叶干质量分别表示叶大小, 研究张掖湿地群落不同密度间芨芨草叶大小-叶脉密度的关系及差异, 旨在理解密度制约下植物叶片的物理构建方式, 有助于揭示湿地异质生境中植物功能结构耦合与权衡的生理生态学机制。

1 材料和方法

1.1 研究区域和样地概况

研究区位于甘肃省张掖市甘州区西城驿黑河干流河床边缘缓坡地带的洪泛平原湿地, 地理位置为38.56° N, 100.26° E, 海拔1482.7 m, 属温带大陆性气候, 雨热同期, 年平均气温7.8 ℃, ≥0 ℃积温为2734 ℃, 年降水量132.6 mm, 主要集中在6-9月, 年蒸发量1986.5 mm, 年日照时间3077 h, 土壤以砾石土、灰棕荒漠土和草甸土为主, 湿地内土壤类型的多样化决定了旱生、盐生、湿生等植物群落的形成, 植被分布具有明显的水平分异特征, 以湿生植物、盐生植物与陆生乔灌林木为主。主要植物有 芨芨草、芦苇(Phragmites australis)、冰草(Agropyron cristatum)、节节草(Equisetum ramosissimum)、宽苞水柏枝(Myricaria bracteata)、多枝柽柳(Tamarix ramosissima)和沙枣(Elaeagnus angustifolia)等。该研究区属于张掖黑河国家湿地保护区核心区, 不存在家畜放牧和人类活动等干扰因素, 在丰水期季节性和间歇性洪水能够满足植物生长的水分需求, 并局部改善湿地的水盐状况; 在枯水期除一些洼地存有水塘外, 大部分区域没有地表水, 植被主要依靠地下潜水生存。在降水稀少、蒸发强烈的大环境影响下, 湿地土壤的盐渍化程度高, 土壤水分空间分异明显, 湿地植被群落的结构组成单一、物种多样性低、寡种优势现象较为普遍, 芨芨草凭借较强的生态可塑性, 形成了不同于河岸林、河流湿地和荒漠的洪泛湿地单优势种群植被景观(占玉芳等, 2012)。

1.2 实验材料

芨芨草为多年生草本植物, 与典型盐生植物相比, 芨芨草并非喜盐植物, 仅是生态幅较宽的耐盐旱生植物, 在海拔较低、地下水位较高的碱性湿地平原以至高达5000 m的青藏高原, 从干草原带一直到荒漠区均有分布(张雅琼等, 2010)。叶丛茂盛, 叶鞘无毛, 较疏松, 具膜质边缘; 叶舌三角形或尖披针形, 长5-10 mm; 叶片扁平或边缘内卷, 质坚韧, 长30-60 cm, 宽7-12 mm, 叶片上表皮可见6-7条沿大小脉呈深波状脉纹突起, 微粗糙有表皮毛, 叶脉大小维管束相间排列, 无主脉, 叶片的中脉与侧脉平行排列, 维管束鞘两层, 下表皮光滑无毛(吴征镒和郭本兆, 1989)。

1.3 实验设置

1.3.1 植物采样

2015年8月初, 在实地考察的基础上, 选择一块芨芨草为优势种的湿地群落, 以“丛”为取样单位进行实验, 按芨芨草株丛(即由种子产生, 由无性繁殖而生成的无性系)的密度状况(张明娟等, 2012), 将芨芨草株丛的密度水平分为4个梯度, 分别为: 高密度(I, >12丛·m-2)、中密度(II, 8-12丛·m-2)和(III, 4-8丛·m-2)、低密度(IV, <4丛·m-2)。按已划分的芨芨草种群的密度梯度各选取7个5 m × 5 m样方, 总计28个样方。首先在每一个样方各选择大中小2丛芨芨草, 共6丛, 用卷尺测量从基部到最高分枝的高度, 定位株高, 并记录芨芨草株丛分枝数, 然后选取每个株丛最外层4个方位充分伸展且健康完整的2-3个叶片, 做好标记后进行以下步骤:

(1)植物光合参数测定。于8月15-18日(代表植物生长旺盛期), 选择晴朗天气的9:00-12:00进行气体交换参数的测定。光合测定使用GFS-3000便携式光合测量系统(Heinz Walz GmbH, Bavaria, Germany), 测量过程中使用人工红蓝光源, 光合有效辐射(PAR)为1200 μmol·m-2·s-1, CO2浓度约为340 μmol·mol-1, 相对湿度(RH)保持在40%-50%, 流速设定为750 μmol·s-1, 叶室面积为8 cm2, 对做好标记的每个叶片记录5组数据用于统计分析, 分别测定叶片的PnTr等参数(任青吉等, 2015), 对于不能充满叶室的叶片, 则保存于湿润的封口袋中, 带回实验室, 用便携式激光叶面积仪(CI-202, Walz, Camas, USA)测定叶面积, 进而计算PnTr

(2)选择晴朗无云的天气用手持光量子计(3415F, 3415F, Walz, Plainfield, USA)于9:00-11:00对4个密度梯度内所有芨芨草进行PAR测定, 测定位置分别在植株上方(距冠层15 cm)、中层(株高1/2处)和地表(距地表15cm处)。

(3)叶脉密度和叶片形态特征的测定。基于每个芨芨草被标记的2-3个叶片, 采集后保存于湿润的封口袋中, 带回实验室用福尔马林-乙酸溶液(37%甲醛溶液、50%乙醇和13%冰醋酸溶液)固定。每片叶片先测量叶脉密度数据, 用含5% NaOH的乙醇进行化学清理, 再用番红-固绿染色(Berlyn & Miksche, 1976), 用水包埋呈透明薄膜状在体视显微镜(SMZ168-BL, Motic, Hong Kong, China)下放大10倍后拍照, 每个叶片拍摄10个视野; 再用Motic Images Plus 2.0软件统计分析叶脉的总长度, 叶脉密度用单位叶面积(单位: mm2)的叶脉总长度(单位: mm)表示(Sack et al., 2012); 最后将叶片放在80 ℃的烘箱中烘干48 h至恒质量, 用电子天平(精度为0.0001 g)称取质量, 记录叶干质量。

1.3.2 土壤取样与理化性质测量

(1)土壤含水量采样。土壤取样于2015年8月下旬完成, 取样时间内基本无明显降水, 土壤水盐状况相对稳定。在每个小样方随机选取3个样点挖掘1 m × 1 m × 0.5 m土壤剖面, 用环刀分5层间隔10 cm采取土样, 现场编号、称鲜质量后带回实验室在105 ℃的烘箱内烘12 h, 取出称质量, 计算出各样地0-50 cm土层土壤质量含水量。

(2)土壤盐分采样。用电导法测定土壤含盐量, 在室温下称取过2 mm筛的风干土样10 g, 加50 mL去CO2蒸馏水(水土比为5:1), 取浸出液, 置振荡机上振荡5 min。将布氏漏斗与抽气系统相连后把悬浊的土浆缓缓倒入漏斗, 直至抽滤完毕, 滤液倒入三角瓶备用(张雅琼等, 2010)。用上海雷磁仪器厂生产的DDS-11C便携式电导仪测浸出液的电导率, 3次重复, 取平均值。

1.4 数据分析

对芨芨草种群样地划分出的4个密度梯度进行各性状统计分析。根据各密度梯度已测定的叶面积和叶干质量的平均值, 计算出芨芨草的比叶面积, 计算公式为: 比叶面积(单位: cm2·g-1) =叶面积/叶干质量。芨芨草叶大小分别指单叶面积和单叶干质量, 对芨芨草叶大小和叶脉密度的实验数据先进行以10为底的对数转换, 使之符合正态分布后再进行分析。对任意两个功能关系的研究, 我们采用y = axb, 线性转换成log(y) = blog(x) + log(a), 式中xy表示两个特征参数, b为斜率, 即异速生长参数或相对生长的指数, 当|b| = 1时, 表示两者是等速变化关系; 当|b|显著偏离1时, 两者间为异速变化关系(Harvey & Pagel, 1991)。密度间异速生长方程的参数估计, 主要采用标准化主轴估计(SMA)方法(Warton et al., 2006), 由软件(S)MATR Version 2.0 (Falster et al., 2012)计算完成, 每一个回归斜率的置信区间根据Pitman (1939)方法计算, 并采用Warton和Weber (2002)的方法对每个密度的斜率进行异质性测试。不同密度间植物功能性状平均值的差异比较采用单因素方差分析(one-way ANOVA, α = 0.05), 实验所用数据采用SPSS 16.0软件进行处理分析, 用SigmaPlot 10.0和Excel软件绘图。

2 结果分析

2.1 不同密度条件下湿地芨芨草光合生理和生物学特征变化

不同密度条件下芨芨草的光合生理和生物学特征的变化如表1所示。随着芨芨草株丛密度的逐渐降低, 芨芨草的株高呈逐渐减小趋势, 从密度(I)到密度(IV)减少了18.84%, 中密度(II、III)无显著差异(p > 0.05); 芨芨草分枝数、PnTr呈先增大后减小的趋势, 从高密度(I)到低密度(IV), 分枝数增加了45.89%, PnTr从高密度(I)到低密度(IV)分别减少12.65%和7.44%, 中密度(II、III)均出现较大值且无显著差异(p > 0.05); 芨芨草的PAR呈逐渐增大趋势, 由高密度(I)的636.30 μmol·m-2·s-1增加到低密度(IV)的1105.10 μmol·m-2·s-1, 增加了0.74倍, 不同密度条件下芨芨草PAR存在显著差异(p < 0.05)。

2.2 不同密度条件下湿地芨芨草叶性状和土壤理化性质的变化

不同密度条件下芨芨草种群叶性状和土壤理化性质的变化如图1所示。由单因素方差分析得出, 芨芨草叶面积、叶干质量和土壤电导率在不同密度存在显著差异(p < 0.05), 叶脉密度、比叶面积和土壤含水量在中密度(II、III)无显著差异(p > 0.05)。随着芨芨草株丛密度的逐渐降低, 芨芨草叶面积和叶干质量呈逐渐减小趋势, 叶脉密度呈逐渐增大趋势, 从高密度(I)到低密度(IV), 叶面积和叶干质量分别减小了47.18%和40%、叶脉密度增加了13.33%; 芨芨草比叶面积随株丛密度的降低呈逐渐减小趋势, 从最大值到最小值减小了12.5%; 土壤电导率随密度的降低呈逐渐增加趋势, 从高密度(I)的3.14 ms·cm-1增加到低密度(IV)的10.90 ms·cm-1, 增加了2.47倍; 土壤含水量呈逐渐减小趋势、从高密度(I)到低密度(IV)减小了71.2%。无论是用叶面积, 还是叶干质量表示叶大小, 不同密度生境叶大小和叶脉密度均呈“此消彼长”的权衡关系。

表1   不同密度条件下芨芨草种群光合生理和生物学特征(平均值±标准误差)

Table 1   Photosynthetic physiological and biological characteristics of Achnatherum splendens under different densities (mean ± SE)

密度
Density
株高
Plant high (cm)
分枝数
Twig number
Pn
(μmol CO2·m-2·s-1)
Tr
(mmol H2O·m-2·s-1)
PAR
(μmol·m-2·s-1)
> 12 bundle·m-2 (I)161.00 ± 10.24a113.33 ± 5.29d13.20 ± 0.12c6.45 ± 0.02c636.30 ± 14.18d
8-12 bundle·m-2 (II)149.67 ± 8.62b203.67 ± 12.34b13.83 ± 0.17a6.58 ± 0.09a839.20 ± 27.95c
4-8 bundle·m-2 (III)140.34 ± 7.50b273.33 ± 15.83a14.12 ± 0.18a6.67 ± 0.10a918.80 ± 36.94b
<4 bundle·m-2 (IV)130.67 ± 5.56c165.33 ± 11.90c13.77 ± 0.13b6.53 ± 0.03b1 105.10 ± 40.62a

PAR, photosynthetically active radiation; Pn, net photosynthetic rate; Tr, transpiration rate. Different lowercase letters in the same column indicate significant differences of different density treatments (p < 0.05).PAR, 光合有效辐射; Pn, 净光合速率; Tr, 蒸腾速率。同列不同小写字母表示不同密度间差异显著(p < 0.05)。

新窗口打开

图1   不同密度条件下芨芨草叶性状和土壤理化性质的变化(平均值±标准误差)。不同小写字母表示相同构件在不同密度处理下差异显著(p < 0.05)。I, 高密度(> 12 bundle·m-2); II, 中密度(8-12 bundle·m-2); III, 中密度(4-8 bundle·m-2); IV, 小密度(<4 bundle·m-2)。

Fig. 1   Change of leaf morphological traits and soil physical and chemical properties of Achnatherum splendens on different densities (mean ± SE). Different lowercase letters indicate significant differences of different density treatments (p < 0.05). I, high density (>12 bundle·m-2); II, medium density (8-12 bundle·m-2); III, medium density (4-8 bundle·m-2); IV, low density (< 4 bundle·m-2).

2.3 不同密度条件下湿地芨芨草叶面积、叶干质量与叶脉密度的关系

2.3.1 不同密度条件下湿地芨芨草叶面积与叶脉密度的关系

图2所示, 芨芨草叶面积与叶脉密度在高密度(I)和低密度(IV)呈极显著负相关关系(yI = -2.37x + 1.61, R2 = 0.88, p < 0.01; yIV = -1.92x + 1.47, R2 = 0.85, p < 0.01), 在中密度(II、III)呈显著负相关关系(yII = -2.25x + 1.47, R2 = 0.78, p < 0.05; yIII = -2.17x + 1.54, R2 = 0.75, p < 0.05); 随着株丛密度逐渐降低, 芨芨草叶脉密度和叶面积回归方程的SMA斜率分别为-2.37 (95%的置信区间CI = (-2.55, -2.27))、-2.25 (95%的置信区间CI = (-2.36, -2.18))、-2.17 (95%的置信区间CI = (-2.27, -1.94))和-1.92 (95%的置信区间CI = (-2.19, -1.70)), 不同密度间二者SMA斜率均与-1存在极显著差异(p < 0.01), 且呈逐渐减小趋势。表明随着株丛密度的降低, 叶面积越小, 叶脉密度越大, 叶面积的生长速率显著小于叶脉密度的增长速度, 且在不同密度生境呈现不同的权衡关系。

图2   不同密度条件下芨芨草叶面积与叶脉密度的关系。I, 高密度(>12 bundle·m-2); II, 中密度(8-12 bundle·m-2); III, 中密度(4-8 bundle·m-2); IV, 低密度(< 4 bundle·m-2)。

Fig. 2   Relationship between leaf area and vein density of Achnatherum splendens under different levels of densities. I, high density (>12 bundle·m-2); II, medium density (8-12 bundle·m-2); III, medium density (4-8 bundle·m-2); IV, low density (<4 bundle·m-2).

图3   不同密度条件下芨芨草叶干质量与叶脉密度的关系。I, 高密度(>12 bundle·m-2); II, 中密度(8-12 bundle·m-2); III, 中密度(4-8 bundle·m-2); IV, 低密度(<4 bundle·m-2)。

Fig. 3   Relationship between leaf dry mass and vein density of Achnatherum splendens under different levels of densities. I, high density (>12 bundle·m-2); II, medium density (8-12 bundle·m-2); III, medium density (4-8 bundle·m-2); IV, low density (<4 bundle·m-2).

2.3.2 不同密度条件下湿地芨芨草叶干质量与叶脉密度的关系

图3所示, 芨芨草叶脉密度与叶干质量在高密度(I)和低密度(IV)呈极显著异速负相关(yI = -1.53x - 0.36, R2 = 0.89, p < 0.01; yIV = -1.13x - 0.59, R2 = 0.85, p < 0.01), 在中密度(II、III)呈显著负相关(yII = -1.41x - 0.47, R2 = 0.75, p < 0.05; yIII = -1.38x - 0.36, R2 = 0.77, p < 0.05); 随着株丛密度的逐渐降低, 芨芨草叶干质量和叶脉密度回归方程的SMA斜率分别为-1.53 (95%的置信区间CI = (-1.88, -1.36))、-1.41 (95%的置信区间CI = (-1.48, -1.26))、-1.38 (95%的置信区间CI = (-1.49, -1.25))和-1.13 (95%的置信区间CI = (-1.23, -1.01)), 叶干质量和叶脉密度在不同密度的SMA斜率均与-1存在显著差异(p < 0.05), 且呈逐渐减小趋势。表明随着密度的降低, 叶干质量越小、叶脉密度越大, 叶干质量的增加速率小于叶脉密度的生长速率, 密度生境中资源的限制加剧了芨芨草叶脉密度和叶干质量的资源权衡机制。

3 讨论

植物的生物量分配格局和生理活性对植株密度的变化具有敏锐的响应性(Cipollini & Bergelson, 2001; 黎磊和周道玮, 2011), 植物个体往往通过调整光合构件的生长速率等长期的表型可塑性来提高植物异质性资源的利用效率(武高林等, 2008)。叶片是植物碳水耦合权衡的重要器官, 叶大小和叶脉密度的权衡机制体现了叶性状对不同密度制约、邻体干扰等过程中的异质性资源的权衡和组合, 从而响应和适应植物的光合蒸腾生理需求, 进而达到生存和生长之间的经济平衡(Funk & Vitousek, 2007)。本研究发现, 高密度(I)和低密度样地(IV)芨芨草的叶大小与叶脉密度之间存在相反的异速生长关系, 且均达到极显著性水平(p < 0.01), 中密度样地(II、III)是芨芨草叶大小逐步减小、叶脉密度逐步增大的过渡区域, 叶大小和叶脉密度仅呈显著负相关关系(p < 0.05)(图1-3), 与Sack等(2012)对485个物种的叶脉密度和叶大小关系的研究结果基本一致。不同密度条件下芨芨草叶大小和叶脉密度生长关系的差异性, 既是密度制约对叶片形态和功能的选择塑造, 也是湿地植物根据自身的生存需求, 在碳投入-收益的权衡下选择的最佳资源配置模式。

植物叶片构型构建与植物生理功能所需的光照、碳投入和水分运输的投资收益相契合, 是植物提高光合收益的一个重要的生活史策略(Westoby et al., 2002; 孙素静等, 2015)。种群密度的制约性调节是决定植物生物量分配格局的关键因素(黎磊和周道玮, 2011), 能够促使植物权衡有限光合产物在叶肉、叶脉等组织间的合理分配(Uhl & Mosbrugger, 1999), 以期实现光合碳获取最大化的目标。研究区属于河西走廊温带干旱区的天然湿地生态系统, 土壤水和地下水是芨芨草可能的主要利用水源(吴华武等, 2015)。处于湿地群落高密度(I)样地的芨芨草种群占有的地上环境空间相对窄, 为了保证有限空间下光拦截能力和光合碳获取能力最大化, 芨芨草以牺牲叶脉密度的资源分配为代价, 选择了增大叶面积和叶干质量的权衡策略和叶片构型构建机制(图1-3)。主要原因有: (1)来自邻株植物的遮阴改变了植物所处光环境的红外/远红外光比率(R/FR), 限制了可获得性PAR (表1), 芨芨草做出增加株高、减少分枝数的遮阴躲避反应(表1), 将更多的生物量投资于单位分枝上的叶片, 形成最大的叶面积和叶肉组织密度, 从而减弱了相邻植株的光抑制作用, 保证了最大的光合作用面积和同化速率, 能够将更多的光能转化为用于自身生长的化学能(黎磊和周道玮, 2011); (2)较高的株丛密度避免了土壤水分的潜在蒸散(吴华武等, 2015), 土壤含水量较高且盐分胁迫弱(图1), 芨芨草进行光合作用的水分输送阻力较小, 较小的叶脉密度也能满足植株较低的光合和蒸腾速率的水分需求(表1), 将更多的生物量用于光转换和同化结构的构建能够极大地缓冲资源有限性约束; (3)该生境芨芨草比叶面积最大(图1), 说明低光胁迫下芨芨草种群为了捕捉更多的光能会倾向用有限的营养物质创造出更大的受光面积以利吸收更多光能, 最大限度地获取光照和控制水分散失之间的平衡(Wright et al., 2004)。因此, 高密度(I)条件下芨芨草叶脉密度和叶大小形成了极显著负相关关系(p < 0.01)。与West等(1997)提出的“自然选择的压力必然使其向最优化的方向进化, 使得与外界环境的交换面积和新陈代谢效率最大化, 并使营养物质运输距离和时间最小化”一致。

最优化分配理论预测, 在无邻体植物竞争时, 植物应该增大获取限制资源的器官分配, 而减少获取非限制资源的器官分配(Weiner, 2004), 即在特定的环境下, 植物为了满足生存需求会形成不同的形态结构和特定的水分供应策略, 使植物更好地调整和适应异质性资源(龚容和高琼, 2015)。湿地芨芨草种群处于较低密度(IV)样地时, 生存空间扩大、光合有效辐射最大(表1)、伴随着较高的光合速率和蒸腾速率(表1), 为了进一步拓展生存生态位, 芨芨草优先倾向于增加分枝数(表1), 从而选择了减少单位分枝上叶面积和叶干质量的资源投资, 将更多的光合产物用于增加叶脉密度以保证高效率的水分输导系统的生存模式(图1-3)。主要原因有: (1)芨芨草为了避免较强的蒸腾作用引起自身失水过多而大量减小叶面积(图1), 这不仅达到了小叶易进行热量和物质交换的目的, 而且减小了蒸腾作用面积, 有利于保持生命活动的水分平衡(党晶晶等, 2014), 形成了最小的叶面积和叶干质量(图1); (2)该生境中土壤含水量最低、土壤电导率大(图1), 盐分胁迫造成的“生理干旱”现象阻碍了植物的吸水作用(周洪华和李卫红, 2015), 芨芨草将更多的叶生物量投资于机械支持与微管结构, 使得该生境的芨芨草叶脉密度较高密度(I)增加了13.33% (图1), 为水分转移至叶肉提供更大的接触面积和便捷的路径, 因此增加叶片中水分运输的量与速率也许是在该生境维持生存繁衍的关键性策略(Nardini et al., 2010); (3)该生境芨芨草比叶面积较高密度(I)减小了12.5% (图1), 即当土壤含水量不足、叶片蒸腾速率较大时, 芨芨草选择较小的叶面积有利于植物获取异质性资源、提高适应贫瘠环境的能力和强光下的自我保护构造(Wright et al., 2004), 这是湿地芨芨草对环境胁迫所产生的一种安全性经济权衡策略。因此, 低密度(IV)样地芨芨草的叶脉密度和叶大小呈极显著的负相关关系(p < 0.01), 这与Sack等(2012)关于“小叶中高的叶脉密度提供冗余的‘高速公路’来运输水分, 有利于提高植物的耐旱能力, 增强叶片抵御虫食干扰能力”的研究结果一致。

植物个体间的相互作用实际上是可获得性资源的有限性导致的构件水平的竞争(黎磊和周道玮, 2011), 在不同密度的资源竞争强度下, 植物的表型可以通过改变不同器官间生物量分配的可塑性反应进行调整(Cheplick, 2006)。芨芨草种群处于中密度(II、III)样地时, 地表土壤理化过程和强度均发生了变化(图1), 种内对空间和光的竞争有所减弱, 株高居中, 芨芨草分枝数、净光合和蒸腾速率最大(表1), 为了增强资源竞争能力和生境适合度, 芨芨草选择了折中光合面积和同化结构的稳步投资, 兼顾了植株的光截取和水分输送等功能需求, 选择了叶大小和叶脉密度均衡生长模式, 可以达到更大限度占有地上部分资源和拓展空间的目的, 同时也能减少植株内部的资源竞争, 这对于种群的稳定与扩张都是有利的。因此, 生长在中密度(II、III)样地的芨芨草叶大小和叶脉密度的相关性仅达到显著水平(p < 0.05, 图2, 图3), 是叶大小由大到小、叶脉密度由小到大转变的区域, 这是由于光合产物在不同构件或功能间权衡分配的结果。

植物叶大小和叶脉密度关系不仅是植物个体生理过程的表现, 也体现了植物叶形态的微生境依赖规律及表型可塑性。湿地群落中高密度(I)样地的芨芨草种群为了适应强烈的光竞争, 权衡叶片组织和水分运输系统的碳投入, 选择减少叶脉密度而增加叶大小资源投资比例的异速生长格局, 以保证稳定的光合收益; 低密度(IV)样地的芨芨草趋向于以牺牲叶大小为代价构建较大的叶脉密度, 保证畅通高效的水分运输路径应对水分胁迫和增加自身的竞争力; 不同种群密度所选择的生存模式均体现了重要的碳-水投资器官如何有效调节资源分配, 即水分运输能力和光合能力的匹配问题。这种随着环境变化的感应机制和“趋利避害”的表型可塑性是湿地植物与环境长期相互协同进化与适应的结果, 可以提高密度制约下植物的生存、生长和繁殖适合度, 也是顺利完成生活史周期的最佳生理和功能状态。本文仅从密度角度分析了湿地生态系统中植物叶大小-叶脉密度的变化特征, 它们的关系可能还受到地下生物量、光合色素含量和物种遗传特性等其他因子的影响, 还需要进一步的实验验证。

The authors have declared that no competing interests exist.

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


参考文献

[1] Berlyn GP, Miksche JP (1976). Botanical Microtechnique and Cytochemistry. Iowa State University Press, Ames, USA.

[本文引用: 1]     

[2] Brodribb TJ, Feild TS, Sack L (2010).

Viewing leaf structure and evolution from a hydraulic perspective.

Functional Plant Biology, 37, 488-498.

[本文引用: 2]     

[3] Brodribb TJ, Jordan GJ (2011).

Water supply and demand remain balanced during leaf acclimation of Nothofagus cunninghamii trees.

New Phytologist, 192, 437-448.

[本文引用: 1]     

[4] Castro-Díez P, Puyravaud JP, Cornelissen JHC (2000).

Leaf structure and anatomy as related to leaf mass per area variation in seedlings of a wide range of woody plant species and types.

Oecologia, 124, 476-486.

[本文引用: 1]     

[5] Cheplick GP (2006).

A modular approach to biomass allocation in an invasive annual (Microstegium vimineum; Poaceae).

American Journal of Botany, 93, 539-545.

[本文引用: 1]     

[6] Cipollini DF, Bergelson J (2001).

Plant density and nutrient availability constrain constitutive and wound-induced expression of trypsin inhibitors in Brassica napus.

Journal of Chemical Ecology, 27, 593-610.

[本文引用: 1]     

[7] Dang JJ, Zhao CZ, Li Y, Hou ZJ, Dong XG (2014).

Variations with slope in stem and leaf traits of Melica przewalskyi in alpine grassland.

Chinese Journal of Plant Ecology, 38, 1307-1314. (in English with Chinese abstract)
[党晶晶, 赵成章, 李钰, 侯兆疆, 董小刚 (2014).

高寒草地甘肃臭草茎-叶性状的坡度差异性

. 植物生态学报, 38, 1307-1314]

[8] Falster DS, Warton DI, Wright IJ (.

URL      [本文引用: 3]     

[9] Funk JL, Vitousek PM (2007).

Resource-use efficiency and plant invasion in low-resource systems.

Nature, 446, 1079-1081.

[本文引用: 1]     

[10] Gong R, Gao Q (2015).

Research progress in the effects of leaf hydraulic characteristics on plant physiological functions.

Chinese Journal of Plant Ecology, 39, 300-308. (in English with Chinese abstract)
[龚容, 高琼 (2015).

叶片结构的水力学特性对植物生理功能影响的研究进展

. 植物生态学报, 39, 300-308.]

[11] Harvey PH, Pagel MD (1991). The Comparative Method in Evolutionary Biology. Oxford University Press, Oxford, UK.

[本文引用: 2]     

[12] Japhet W, Zhou DW, Zhang HX, Zhang HX, Yu T (2009).

Evidence of phenotypic plasticity in the response of Fagopyrum esculentum to population density and sowing date.

Journal of Plant Biology, 52, 303-311.

[13] Jin Y, Wang CK (2015).

Trade-offs between plant leaf hydraulic and economic traits.

Chinese Journal of Plant Ecology, 39, 1021-1032. (in Chinese with English abstract)
[金鹰, 王传宽 (2015).

植物叶片水力与经济性状权衡关系的研究进展

. 植物生态学报, 39, 1021-1032.]

[本文引用: 1]     

[14] Li L, Zeng H, Guo DL (2013).

Leaf venation functional traits and their ecological significance.

Chinese Journal of Plant Ecology, 37, 691-698. (in English with Chinese abstract)
[李乐, 曾辉, 郭大立 (2013).

叶脉网络功能性状及其生态学意义

. 植物生态学报, 37, 691-698.]

[15] Li L, Zhou DW (2011).

Density-dependent regulation of above- and below-ground modules in Allium cepa var. proliferum populations.

Chinese Journal of Plant Ecology, 35, 284-293. (in English with Chinese abstract)
[黎磊, 周道玮 (2011).

红葱种群地上和地下构件的密度制约调节

. 植物生态学报, 35, 284-293.]

[本文引用: 2]     

[16] Liu MH, Xin ZM, Xu J, Sun F, Dou LJ, Li YH (2013).

Influence of leaf size of plant on leaf transpiration and temperature in arid regions.

Chinese Journal of Plant Ecology, 37, 436-442. (in English with Chinese abstract)
[刘明虎, 辛智鸣, 徐军, 孙非, 窦立军, 李永华 (2013).

干旱区植物叶片大小对叶表面蒸腾及叶温的影响

. 植物生态学报, 37, 436-442.]

[17] Mcdonald PG, Fonseca CR, Overton JM, Westoby M (2003).

Leaf-size divergence along rainfall and soil-nutrient gradients: Is the method of size reduction common among clades?

Functional Ecology, 17, 50-57.

[本文引用: 1]     

[18] Nardini A, Raimondo F, Lo Gullo MA (2010).

Leafminers help us understand leaf hydraulic design. Plant,

Cell & Environment, 33, 1091-1100.

[19] Niklas KJ, Cobb ED, Niinemets Ü, Reich PB, Sellin A, Shipley B, Wright IJ (2007).

“Diminishing returns” in the scaling of functional leaf traits across and within species groups.

Proceedings of the National Academy of Sciences of the United States of America, 104, 8891-8896.

[本文引用: 4]     

[20] Pitman EJG (1939).

A note on normal correlation.

Biometrika, 31, 9-12.

[21] Reich PB, Cornelissen H (2014).

The world-wide ‘fast-slow’ plant economics spectrum: A traits manifesto.

Journal of Ecology, 102, 275-301.

[本文引用: 1]     

[22] Ren QJ, Li HL, Bu HY (2015).

Comparison of physiological and leaf morphological traits for photosynthesis of the 51 plant species in the Maqu alpine swamp meadow.

Chinese Journal of Plant Ecology, 39, 593-603.
[任青吉, 李宏林, 卜海燕 (2015).

玛曲高寒沼泽化草甸51种植物光合生理和叶片形态特征的比较

. 植物生态学报, 39, 593-603.]

[本文引用: 1]     

[23] Sack L, Frole K (2006).

Leaf structural diversity is related to hydraulic capacity in tropical rain forest trees.

Ecology, 87, 483-491.

[本文引用: 2]     

[24] Sack L, Scoffoni C (2013).

Leaf venation: Structure, function, development, evolution, ecology and applications in the past, present and future.

New Phytologist, 198, 983-1000.

[本文引用: 1]     

[25] Sack L, Scoffoni C, McKown AD (2012).

Developmentally based scaling of leaf venation architecture explains global ecological patterns.

Nature Communications, 3, 837.

[本文引用: 1]     

[26] Shi YC, Zhao CZ, Song QH, Du J, Chen NJ, Wang JW (2015).

Slope-related variations in twig and leaf traits of Robinia pseudoacacia in the northern mountains of Lanzhou.

Chinese Journal of Plant Ecology, 39, 362-370. (in English with Chinese abstract)
[史元春, 赵成章, 宋清华, 杜晶, 陈静, 王继伟 (2015).

兰州北山刺槐枝叶性状的坡向差异性

. 植物生态学报, 39, 362-370.]

[本文引用: 2]     

[27] Silvertown J, Charlesworth D (2001).

Introduction to Plant Population Biology. Blackwell

, London.

[本文引用: 1]     

[28] Sun SJ, Li FL, Bao WK (2015).

Advances on construction of leaf venation system and its significance of phylogeny.

Journal of Tropical and Subtropical Botany, 23, 353-360. (in English with Chinese abstract)
[孙素静, 李芳兰, 包维楷 (2015).

叶脉网络系统的构建和系统学意义研究进展

. 热带亚热带植物学报, 23, 353-360.]

[本文引用: 1]     

[29] Uhl D, Mosbrugger V (1999).

Leaf venation density as a climate and environmental proxy: A critical review and new data.

Palaeoecology, 149, 15-26.

[本文引用: 3]     

[30] 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.]

[本文引用: 3]     

[31] Wang J, Cheng JM, Wan HE, Fang F (2004).

Study on soil moisture characteristics and water use efficiency of Achnatherum splendens grassland in Loess Plateau.

Arid Meteorology, 22, 51-55.

[32] Warton DI, Weber NC (2002).

Common slope tests for bivariate errors-in-variables models.

Biometrical Journal, 44, 161-174.

[本文引用: 1]     

[33] Warton DI, Wright IJ, Falster DS, Westoby M (2006).

Bivariate line-fitting methods for allometry.

Biological Reviews, 81, 259-291.

[本文引用: 1]     

[34] Weiner J (2004).

Allocation, plasticity and allometry in plants. Perspectives in Plant Ecology,

Evolution and Systematics, 6, 207-215.

[35] West GB, Brown JH, Enquist BJ (1997).

A general model for the origin of allometric scaling laws in biology.

Science, 276, 122-126.

[本文引用: 1]     

[36] Westoby M, Falster DS, Moles AT (2002).

Plant ecological strategies: Some leading dimensions of variation between species.

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

[本文引用: 1]     

[37] Wright IJ, Reich PB, Westoby M, Ackerly DD, Baruch Z, Bongers F, Cavender-Bares J, Chapin T, Cornelissen JH, 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 world- wide leaf economics spectrum.

Nature, 428, 821-827.

[38] Wu GL, Chen M, Du GZ (2008).

Effects of nutrient and lighton seedlings morphological plasticity of four Saussurea species with different ecological breadth.

Chinese Journal of Applied Ecology, 19, 1708-1713. (in Chinese with English abstract)
[武高林, 陈敏, 杜国祯 (2008).

营养和光照对不同生态幅风毛菊属植物幼苗形态可塑性的影响

. 应用生态学报, 19, 1708-1713.]

[本文引用: 1]     

[39] Wu HW, Li XY, Jing ZY, Li J, Zheng XR, Zhao DZ (2015).

Variations in water use for Achnatherum splendens in Lake Qinghai watershed based on δD and δ18O.

Acta Ecologica Sinica, 35, 8174-8183. (in English with Chinese abstract)
[吴华武, 李小雁, 蒋志云, 李静, 郑肖然, 赵殿智 (2015).

基于δD、δ18O的青海湖流域芨芨草水分利用来源变化研究

. 生态学报, 35, 8174-8183.]

[本文引用: 1]     

[40] Wu ZS, Guo BZ (1989). Flora Repubulicae Popularis Sinicae. (Tomus 9). Science Press, Beijing. 1087. (in Chinese) [吴征镒, 郭本兆 (1989). 中国植物志(第九卷). 科学出版社, 北京. 1087.]

[本文引用: 1]     

[41] Yang DM, Zhang JJ, Zhou D, Qian MJ, Zheng Y, Jin LM (2012).

Leaf and twig functional traits of woody plants and their relationships with environmental change.

Chinese Journal of Ecology, 31, 702-713. (in English with Chinese abstract)
[杨冬梅, 章佳佳, 周丹, 钱敏杰, 郑瑶, 金灵妙 (2012).

木本植物茎叶功能性状及其关系随环境变化的研究进展

. 生态学杂志, 31, 702-713.]

[本文引用: 1]     

[42] Zhan YF, Ma L, Li XY, Teng YF, Lu YF, Tian XP (2012).

Niches of dominant populations of wetland plants in Zhangye section of Heihe River.

Journal of Northeast Forestry University, 40(10), 61-66. (in English with Chinese abstract)
[占玉芳, 马力, 李小燕, 滕玉风, 鲁延芳, 田晓萍 (2012).

黑河流域(张掖段)湿地植物群落优势种群生态位

. 东北林业大学学报, 40(10), 61-66.]

[本文引用: 1]     

[43] Zhang MJ, Liu MS, Xu C, Chi T, Hong C (2012).

Spatial pattern responses of Achnatherum splendens to environmental stress in different density levels.

Acta Ecologica Sinica, 32, 595-604. (in English with Chinese abstract)
[张明娟, 刘茂松, 徐驰, 池婷, 洪超 (2012).

不同密度条件下芨芨草空间格局对环境胁迫的响应

. 生态学报, 32, 595-604.]

[44] Zhang SB, Guan ZJ, Sun M, Zhang JJ, Cao KF, Hu H (2012).

Evolutionary association of stomatal traits with leaf vein- density in Paphiopedilum, Orchidaceae.

PLOS ONE, 7, e40080. doi: 10.1371/journal.pone.0040080.

[本文引用: 1]     

[45] Zhang YF, Wang W, Liang CZ, Wang LX, Pei H, Wang CY, Wang WF (2012).

Suitable habitat for the Achnatherum splendens community in typical steppe region of Inner Mongolia.

Acta Ecologica Sinica, 32, 1193-1201. (in English with Chinese abstract)
[张翼飞, 王炜, 梁存柱, 王立新, 裴浩, 王成燕, 王伟峰 (2012).

内蒙古典型草原区芨芨草群落适生生境

. 生态学报, 32, 1193-1201.]

[本文引用: 3]     

[46] Zhang YQ, Liang CZ, Wang W, Wang LX, Peng JT, Yan JC, Jia JC (2010).

Soil salinity and Achnatherum splendens distribution.

Chinese Journal of Ecology, 29, 2438-2443. (in English with Chinese abstract)
[张雅琼, 梁存柱, 王炜, 王立新, 彭江涛, 闫建成, 贾成朕 (2010).

芨芨草群落土壤盐分特征

. 生态学杂志, 29, 2438-2443.]

[47] Zhou HH, Li WH (2105).

Responses and adaptation of xylem hydraulic conductivity to salt stress in Populus euphratica.

Chinese Journal of Plant Ecology, 39, 81-91. (in English with Chinese abstract)
[周洪华, 李卫红 (2015).

胡杨木质部水分传导对盐胁迫的响应与适应

. 植物生态学报, 39, 81-91.]

[本文引用: 1]     

/