基于野外调查和同质种植园实验的芦苇植物功能性状变异研究

ＡＣＴＡＥＣＯＬＯＧＩＣＡＳＩＮＩＣＡ

生态学报

Ｖｏｌ．４１，Ｎｏ．１０Ｍａｙ，２０２１

ＤＯＩ：１０．５８４６／ｓｔｘｂ２０２００９２２２４６８

马香艳，刘乐乐，尹美淇，宋慧佳，朱鹏程，于晓娜，杜宁，王仁卿，郭卫华．基于野外调查和同质种植园实验的芦苇植物功能性状变异研究．生态学报，２０２１，４１（１０）：３７５５⁃３７６４．

ＭａＸＹ，ＬｉｕＬＬ，ＹｉｎＭＱ，ＳｏｎｇＨＪ，ＺｈｕＰＣ，ＹｕＸＮ，ＤｕＮ，ＷａｎｇＲＱ，ＧｕｏＷＨ．ＶａｒｉａｔｉｏｎｉｎｐｌａｎｔｆｕｎｃｔｉｏｎａｌｔｒａｉｔｓｏｆＰｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓｂａｓｅｄｏｎｆｉｅｌｄｉｎｖｅｓｔｉｇａｔｉｏｎａｎｄｃｏｍｍｏｎｇａｒｄｅｎｅｘｐｅｒｉｍｅｎｔ．ＡｃｔａＥｃｏｌｏｇｉｃａＳｉｎｉｃａ，２０２１，４１（１０）：３７５５⁃３７６４．

基于野外调查和同质种植园实验的芦苇植物功能性状变异研究

马香艳１，２，３，刘乐乐１，２，３，尹美淇１，２，３，宋慧佳１，２，３，朱鹏程１，２，３，于晓娜１，２，３，杜　宁１，２，３，王仁卿１，２，３，郭卫华１，２，３，∗

１山东大学生命科学学院生态学与生物多样性研究所，青岛　２６６２３７２山东省植被生态工程技术研究中心，青岛　２６６２３７３山东青岛森林生态系统定位研究站，青岛　２６６２３７

摘要：研究表型可塑性和遗传变异在植物表型分化中的相对作用，有助于预测全球环境变化下的植物群落组成和生态系统功能的变化。芦苇（Ｐｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓ）是全球性广布的草本植物，种内变异丰富，在我国西北和东部均存在多个分化稳定的生态型，但中国芦苇在更大尺度上的表型研究还非常匮乏。将位于黄河上游的宁夏平原和黄河下游的黄河三角洲作为研究区域，通过野外调查和同质种植园实验对芦苇自然种群的植物功能性状变异进行观测。结果表明，无论在野外还是同质种植园，黄河三角洲芦苇的基径、叶长和叶宽均显著大于宁夏平原芦苇，说明两个地区的芦苇种群之间存在着受遗传决定的表型分化，这可能与两个地区间的降水等气候差异有关。在野外，宁夏芦苇的株高和叶厚显著大于黄河三角洲芦苇，但在同质园中差异消失或相反，说明株高、叶厚受环境影响较大，表型可塑性也是芦苇适应环境变化的重要机制。在同质种植园中，宁夏平原芦苇的叶片氮磷含量较低，但株数却显著多于黄河三角洲芦苇，反映了不同地区芦苇之间存在不同的适应策略，宁夏平原芦苇更偏向于高扩散率的杂草策略，而黄河三角洲芦苇更偏向于竞争策略。此外，宁夏平原芦苇的株高、叶长两个性状以及基径－比叶面积相关性在野外和同质园两个环境中存在一致性，表明了性状变异和权衡策略的遗传稳定性。综上，位于黄河上下游的芦苇种群间存在着适应性分化，这是表型可塑性和遗传变异共同作用的结果，不同来源芦苇对全球变化下的多重环境因子的响应还需要进一步研究。

关键词：芦苇；植物功能性状；表型可塑性；同质种植园

ＶａｒｉａｔｉｏｎｉｎｐｌａｎｔｆｕｎｃｔｉｏｎａｌｔｒａｉｔｓｏｆＰｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓｂａｓｅｄｏｎｆｉｅｌｄｉｎｖｅｓｔｉｇａｔｉｏｎａｎｄｃｏｍｍｏｎｇａｒｄｅｎｅｘｐｅｒｉｍｅｎｔ

ＭＡＸｉａｎｇｙａｎ１，２，３，ＬＩＵＬｅｌｅ１，２，３，ＹＩＮＭｅｉｑｉ１，２，３，ＳＯＮＧＨｕｉｊｉａ１，２，３，ＺＨＵＰｅｎｇｃｈｅｎｇ１，２，３，ＹＵＸｉａｏｎａ１，２，３，

１ＩｎｓｔｉｔｕｔｅｏｆＥｃｏｌｏｇｙａｎｄＢｉｏｄｉｖｅｒｓｉｔｙ，ＳｃｈｏｏｌｏｆＬｉｆｅＳｃｉｅｎｃｅｓ，ＳｈａｎｄｏｎｇＵｎｉｖｅｒｓｉｔｙ，Ｑｉｎｇｄａｏ２６６２３７，Ｃｈｉｎａ

ＤＵＮｉｎｇ１，２，３，ＷＡＮＧＲｅｎｑｉｎｇ１，２，３，ＧＵＯＷｅｉｈｕａ１，２，３，∗

２ＳｈａｎｄｏｎｇＰｒｏｖｉｎｃｉａｌＥｎｇｉｎｅｅｒｉｎｇａｎｄＴｅｃｈｎｏｌｏｇｙＲｅｓｅａｒｃｈＣｅｎｔｅｒｆｏｒＶｅｇｅｔａｔｉｏｎＥｃｏｌｏｇｙ，Ｑｉｎｇｄａｏ２６６２３７，Ｃｈｉｎａ

３ＱｉｎｇｄａｏＦｏｒｅｓｔＥｃｏｌｏｇｙＲｅｓｅａｒｃｈＳｔａｔｉｏｎｏｆＮａｔｉｏｎａｌＦｏｒｅｓｔｒｙａｎｄＧｒａｓｓｌａｎｄＡｄｍｉｎｉｓｔｒａｔｉｏｎ，Ｑｉｎｇｄａｏ２６６２３７，Ｃｈｉｎａ

Ａｂｓｔｒａｃｔ：Ｓｔｕｄｙｏｎｔｈｅｒｅｌａｔｉｖｅｃｏｎｔｒｉｂｕｔｉｏｎｓｏｆｐｈｅｎｏｔｙｐｉｃｐｌａｓｔｉｃｉｔｙａｎｄｇｅｎｅｔｉｃｖａｒｉａｔｉｏｎｔｏｐｌａｎｔｐｈｅｎｏｔｙｐｉｃｄｉｆｆｅｒｅｎｔｉａｔｉｏｎｃａｎｆａｃｉｌｉｔａｔｅｔｈｅｐｒｅｄｉｃｔｉｏｎｓｏｆｃｏｍｍｕｎｉｔｙｃｏｍｐｏｓｉｔｉｏｎａｎｄｅｃｏｓｙｓｔｅｍｆｕｎｃｔｉｏｎｉｎｇｕｎｄｅｒｇｌｏｂａｌｅｎｖｉｒｏｎｍｅｎｔａｌ

基金项目：国家重点研发计划课题（２０１７ＹＦＣ０５０５９０５）；国家自然科学基金项目（３１９７０３４７，３１７７０３６１）收稿日期：２０２０⁃０９⁃０２；　　修订日期：２０２１⁃０４⁃２０

∗通讯作者Ｃｏｒｒｅｓｐｏｎｄｉｎｇａｕｔｈｏｒ．Ｅ⁃ｍａｉｌ：ｗｈｇｕｏ＿ｓｄｕ＠１６３．ｃｏｍ

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

３７５６　生　态　学　报　　　４１卷　

ｃｈａｎｇｅｓ．Ａｓａｃｏｓｍｏｐｏｌｉｔａｎｇｒａｓｓ，ｃｏｍｍｏｎｒｅｅｄ（Ｐｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓ）ｈａｓａｈｉｇｈｉｎｔｒａｓｐｅｃｉｆｉｃｄｉｖｅｒｓｉｔｙａｎｄｈａｓｓｅｖｅｒａｌ

ｅｃｏｔｙｐｅｓｉｎｎｏｒｔｈｗｅｓｔａｎｄｅａｓｔＣｈｉｎａ．Ｈｏｗｅｖｅｒ，ｆｅｗｓｔｕｄｉｅｓｆｏｃｕｓｅｄｏｎｔｈｅｐｈｅｎｏｔｙｐｉｃｖａｒｉａｔｉｏｎｏｆｃｏｍｍｏｎｒｅｅｄａｔａｌａｒｇｅｒａｎｕｐｓｔｒｅａｍｗｅｔｌａｎｄ，ａｎｄｔｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａ，ａｄｏｗｎｓｔｒｅａｍｗｅｔｌａｎｄ，ｖｉａｆｉｅｌｄｉｎｖｅｓｔｉｇａｔｉｏｎａｎｄｃｏｍｍｏｎｇａｒｄｅｎ

ｓｃａｌｅｉｎＣｈｉｎａ．Ｈｅｒｅ，ｗｅｉｎｖｅｓｔｉｇａｔｅｄｔｈｅｉｎｔｒａｓｐｅｃｉｆｉｃｖａｒｉａｔｉｏｎｉｎｐｌａｎｔｆｕｎｃｔｉｏｎａｌｔｒａｉｔｓｏｆＰ．ａｕｓｔｒａｌｉｓｆｒｏｍＮｉｎｇｘｉａＰｌａｉｎ，ｅｘｐｅｒｉｍｅｎｔ．Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔｔｈｅｓｈｏｏｔｄｉａｍｅｔｅｒ，ｌｅａｆｌｅｎｇｔｈａｎｄｌｅａｆｗｉｄｔｈｏｆＰ．ａｕｓｔｒａｌｉｓｆｒｏｍｔｈｅＹｅｌｌｏｗＲｉｖｅｒｓｔａｂｌｅｐｈｅｎｏｔｙｐｉｃｄｉｆｆｅｒｅｎｔｉａｔｉｏｎｃｏｎｔｒｏｌｌｅｄｂｙｇｅｎｅｔｉｃｖａｒｉａｔｉｏｎｂｅｔｗｅｅｎｔｈｅｔｗｏｒｅｇｉｏｎｓ．Ｔｈｅｐｈｅｎｏｔｙｐｉｃｄｉｆｆｅｒｅｎｔｉａｔｉｏｎｃｏｕｌｄｂｅｅｘｐｌａｉｎｅｄｂｙｔｈｅａｄａｐｔａｔｉｏｎｔｏｌｏｃａｌｃｌｉｍａｔｉｃｃｏｎｄｉｔｉｏｎｓｓｕｃｈａｓｐｒｅｃｉｐｉｔａｔｉｏｎ．Ｉｎｔｈｅｆｉｅｌｄ，ｔｈｅｓｈｏｏｔｈｅｉｇｈｔａｎｄｔｈｅｌｅａｆ

ＤｅｌｔａｗｅｒｅｓｉｇｎｉｆｉｃａｎｔｌｙｌａｒｇｅｒｔｈａｎｔｈｏｓｅｆｒｏｍｔｈｅＮｉｎｇｘｉａＰｌａｉｎｉｎｂｏｔｈｔｈｅｆｉｅｌｄａｎｄｔｈｅｃｏｍｍｏｎｇａｒｄｅｎ，ｓｕｇｇｅｓｔｉｎｇａ

ｔｈｉｃｋｎｅｓｓｏｆＰ．ａｕｓｔｒａｌｉｓｆｒｏｍＮｉｎｇｘｉａＰｌａｉｎｗｅｒｅｓｉｇｎｉｆｉｃａｎｔｌｙｈｉｇｈｅｒｔｈａｎｔｈｏｓｅｆｒｏｍｔｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａ，ｂｕｔｔｈｅｄｉｆｆｅｒｅｎｃｅｄｉｓａｐｐｅａｒｅｄｏｒｒｅｖｅｒｓｅｄｉｎｔｈｅｃｏｍｍｏｎｇａｒｄｅｎ，ｗｈｉｃｈｓｕｇｇｅｓｔｓｔｈａｔｔｈｅｓｅｔｒａｉｔｓａｒｅｒｅｇｕｌａｔｅｄｄｉｒｅｃｔｌｙｂｙｔｈｅｃｕｒｒｅｎｔｅｎｖｉｒｏｎｍｅｎｔａｎｄｐｈｅｎｏｔｙｐｉｃｐｌａｓｔｉｃｉｔｙｉｓａｌｓｏａｎｉｍｐｏｒｔａｎｔｍｅｃｈａｎｉｓｍｆｏｒＰ．ａｕｓｔｒａｌｉｓｔｏａｄａｐｔｔｏｅｎｖｉｒｏｎｍｅｎｔａｌｌｏｗｅｒ，ｗｈｉｌｅｔｈｅｐｌａｎｔｄｅｎｓｉｔｙｗａｓｓｉｇｎｉｆｉｃａｎｔｌｙｌａｒｇｅｒｔｈａｎｔｈａｔｆｒｏｍｔｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａ．ＷｈｉｃｈｒｅｆｌｅｃｔｓｄｉｆｆｅｒｅｎｔａｄａｐｔｉｖｅｓｔｒａｔｅｇｉｅｓｔｈａｔＰ．ａｕｓｔｒａｌｉｓｆｒｏｍＮｉｎｇｘｉａＰｌａｉｎｐｒｅｆｅｒｓａｒｕｄｅｒａｌｓｔｒａｔｅｇｙｗｉｔｈａｈｉｇｈｒｅｐｒｏｄｕｃｔｉｏｎｒａｔｅ，ｗｈｉｌｅＰ．ｃｈａｎｇｅｓ．Ｉｎｔｈｅｃｏｍｍｏｎｇａｒｄｅｎ，ｔｈｅｃｏｎｔｅｎｔｏｆｎｉｔｒｏｇｅｎａｎｄｐｈｏｓｐｈｏｒｕｓｉｎｌｅａｖｅｓｏｆＰ．ａｕｓｔｒａｌｉｓｆｒｏｍＮｉｎｇｘｉａＰｌａｉｎｗａｓ

ａｕｓｔｒａｌｉｓｆｒｏｍｔｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａｆａｖｏｒｓａｃｏｍｐｅｔｉｔｉｖｅｓｔｒａｔｅｇｙ．Ｍｏｒｅｏｖｅｒ，ｔｈｅｒｅｗｅｒｅｃｏｒｒｅｌａｔｉｏｎｓｂｅｔｗｅｅｎｔｈｅｆｉｅｌｄａｎｄｔｈｅｃｏｍｍｏｎｇａｒｄｅｎｏｂｓｅｒｖａｔｉｏｎｓｉｎｓｈｏｏｔｈｅｉｇｈｔ，ｌｅａｆｌｅｎｇｔｈａｎｄｒｅｌａｔｉｏｎｓｈｉｐｂｅｔｗｅｅｎｓｈｏｏｔｄｉａｍｅｔｅｒａｎｄｓｐｅｃｉｆｉｃｌｅａｆａｒｅａ，ｗｈｉｃｈｉｎｄｉｃａｔｅｓｔｈｅｇｅｎｅｔｉｃｉｎｈｅｒｉｔａｎｃｅｏｆｔｒａｉｔｖａｒｉａｔｉｏｎａｎｄｔｒａｄｅ－ｏｆｆｓｔｒａｔｅｇｙ．Ｉｎｃｏｎｃｌｕｓｉｏｎ，ｔｈｅｒｅｉｓａｎａｄａｐｔｉｖｅｄｉｆｆｅｒｅｎｔｉａｔｉｏｎｂｅｔｗｅｅｎｔｈｅｕｐｓｔｒｅａｍａｎｄｄｏｗｎｓｔｒｅａｍｐｏｐｕｌａｔｉｏｎｓｏｆＰ．ａｕｓｔｒａｌｉｓｉｎｔｈｅＹｅｌｌｏｗＲｉｖｅｒｂａｓｉｎｄｕｅｔｏｐｈｅｎｏｔｙｐｉｃｐｌａｓｔｉｃｉｔｙａｎｄｇｅｎｅｔｉｃｖａｒｉａｔｉｏｎ．Ｆｕｒｔｈｅｒｓｔｕｄｉｅｓａｒｅｎｅｅｄｅｄｔｏｒｅｖｅａｌｔｈｅｒｅｓｐｏｎｓｅｓｏｆｐｌａｎｔｆｒｏｍｄｉｆｆｅｒｅｎｔｐｒｏｖｅｎａｎｃｅｓｔｏｆｕｔｕｒｅｃｈａｎｇｅｓｏｆｔｈｅｍｕｌｔｉｐｌｅｅｎｖｉｒｏｎｍｅｎｔａｌｆａｃｔｏｒｓ．

ＫｅｙＷｏｒｄｓ：Ｐｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓ；ｐｌａｎｔｆｕｎｃｔｉｏｎａｌｔｒａｉｔ；ｐｈｅｎｏｔｙｐｅｐｌａｓｔｉｃｉｔｙ；ｃｏｍｍｏｎｇａｒｄｅｎ

自然界中的植物个体表型会受到遗传因素与环境因素的共同影响，从发育和遗传机制来看，种内性状可以分为遗传变异主导的性状和表型可塑性主导的性状两类［１⁃２］，前者产生的变异是由遗传基因决定且可稳定遗传的，后者产生的变异是由环境直接塑造的，即同一基因型的植物个体在不同的生长环境中表现出不同的表型［３］，但是短时间内产生的这种表型变异通常不能被后代继承。虽然这两种机制存在较大的差别，但它们协同共存，都是生物适应异质环境的重要手段［４］。在全球气候变化背景下，植物对多重环境变化的响应及适应机制成为生态学研究的热点［５⁃６］，植物表型变异及其内部机制的研究不仅有助于解析植物对环境的适应策略，还能对物种未来分布和进化的预测模型提供理论依据。

植物功能性状是植物应对环境变化形成的一系列形态、生理、物候等特征属性，这些属性通常是在植物中比较常见的表型，可以直接或者间接地反映生存环境的变化［７⁃８］。比如，随着气温升高，植物叶片的氮磷含量下降，氮磷比呈现上升趋势［９］；与湿润生境相比，处于干旱生境中的植物叶片单位面积叶干重较大，叶片气孔４８个植物物种的出叶时间、叶片大小、叶片单位面积质量等叶片性状间的相关性进行研究时，发现了其性状受到了物种系统发育差异的强烈影响［１１］，遗传因素在植物生长的整个过程中也发挥着重要的作用。此外，当环境发生变化时，植物性状的改变往往不是单一的，而是多个植物功能性状协同变化的，性状间存在复杂的协同和权衡关系，共同维持和调节植物的生命活动［４，１２］，植物的叶片氮含量、比叶面积等与资源获取相关的叶性状之间联系紧密，叶片氮含量升高，往往伴随着比叶面积的升高［１３］。植物功能性状既相互联系又存在较大的种内与种间差异，普定喀斯特地区的木本植物相比于同纬度其他生境植物来说，有着较小的叶面积和较大的叶组织密度等性状特征，这些有利于减少蒸腾和增加养分储存的性状组合可能是其适应干旱贫瘠的喀斯特环境而采取的生态策略［１４］。植物功能性状变异性及性状相关性的研究对于理解植物应对全球变化、环境监

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

导度较低，有利于植物提高水分利用效率［１０］。同时，植物功能性状也受谱系影响，研究者对中国东部森林中

　１０期　　　马香艳　等：基于野外调查和同质种植园实验的芦苇植物功能性状变异研究　３７５７

测与评价、生态保护和恢复等理论研究和实践活动具有重要意义［１５］。

广泛分布的物种所处的环境多样，往往具有大量的表型变异，由于局部适应还会形成稳定的种内分

化［１６］。芦苇（Ｐｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓ）作为一种世界广布的湿地植物，主要通过风媒异交的种子进行有性生殖，也可以通过根状茎进行营养繁殖，具有非常高的种内遗传多样性和表型可塑性，而且其生态适应幅度很宽，可以生长在内陆湖泊、河口和滨海滩涂等浅水湿地，也可以生长在盐碱、荒漠等恶劣的环境中［１７⁃１９］。芦苇生境的多样性导致了其丰富的生态型，在中国河西走廊地区根据生境不同将芦苇划分为沙丘芦苇、沼泽芦苇、盐化草３个形态类型［２１］。在全球尺度上，芦苇的株高和单株最大生物量随着纬度的增加而减少，而这两个性状的表型可塑性随着纬度的增加而增加，这种性状和表型可塑性变异都是对环境适应的结果［２２］。然而，现在对中国芦苇的研究只局限于小区域，缺少对中国芦苇在较大尺度上的对比研究，这限制了对中国芦苇表型分化与生态适应的深入认识。黄河上游的宁夏平原和黄河下游的黄河三角洲都处于黄河流域，具有连通性，同时跨越中国西部和东部不同的气候区域，又具有气候差异性，这为我们在更大尺度上的研究提供了良好的地理条件。同时，芦苇作为一种世界广布的植物，其高度的表型变异和丰富的生态型为我们研究植物种内表型差异提供了良好的实验材料。目前已有较多研究探讨了芦苇适应异质生境的策略［２３］，但是从表型可塑性和遗传变异两方面阐明其适应机制的报道却相对较少。

本实验以芦苇为研究对象，将野外调查测量实验和同质园移栽实验结合起来，比较分析黄河上游宁夏平原和黄河下游黄河三角洲两个区域芦苇的植物功能性状差异，揭示芦苇在不同气候区域的表型分化和不同性状之间的潜在联系，加深对芦苇的适应性进化机制的认识。１　研究方法１．１　野外调查和采样

宁夏平原位于黄河上游（３５．８３°Ｎ—３９．３８°Ｎ，１０４．２８°Ｅ—１０７．２０°Ｅ），处于中国西部内陆地区，为典型的温带大陆性气候，昼夜温差和气温年较差大，７月气温最高，平均为２０．０—２３．１℃，１月气温最低，平均为－８．２—Ｎ—３８．１２°Ｎ，１１８．０７°Ｅ—１１９．１８°Ｅ），处于中国东部临海区域，具有中国暖温带最广阔、最完整的新生河口湿地，为温带大陆季风性气候，年平均气温约为１２．１℃，年平均降水量约为５５１．６ｍｍ，四季分明，光照充足，雨热同期［２６］，而且水盐生境异质性强，具有许多的芦苇单优群落［２１］。

２０１８年７月下旬和８月上旬，分别在黄河三角洲和宁夏平原两个区域各调查３０个芦苇种群，种群之间距－６．６℃，平均年降水量为２６２．３—３６３．２ｍｍ，干旱少雨，生态环境脆弱［２４⁃２５］。黄河三角洲位于黄河下游（３６．５５°甸芦苇等［２０］，在黄河下游根据芦苇种群形态的变异规律，将该地区芦苇划分为盐生芦苇、淡水芦苇、巨型芦苇

离大于１０ｋｍ，每个种群内选取５棵健康植株，个体间相距约２０ｍ，即在宁夏平原和黄河三角洲各选取１５０株植株个体进行植物功能性状的测定。所测的功能性状和测量方法如下，植株地面位置到植株最大高度的距离为株高；靠近地面或水面的芦苇节间的直径为基径；每个区域中抽穗的芦苇种群数量除以总的种群数量为抽穗占比。在所选植株上随机选取３片完全展开的、生长健康的叶片，叶片最后端到叶尖间的距离记为叶长；叶片最宽处的距离为叶宽；测完叶长、叶宽后将３片叶子叠在一起，对折两次后测量叶片平均厚度为叶厚。通过扫描、图像处理程序ＩｍａｇｅＪｖ１．８．０确定叶面积，然后将叶片放入８０℃烘箱中４８ｈ烘至恒重得到叶干重，计算得出芦苇叶片的比叶面积（Ｓｐｅｃｉａｌｌｅａｆａｒｅａ，ＳＬＡ）。最后沿着地面或水面收获所测植株的地上部分，烘至恒重得到植株干重，再加上之前测得的３片叶子的叶干重即为芦苇整株地上生物量。根据野外测量地点的情况，收集部分种群中所测芦苇的根状茎，将其带回实验室用于后续的同质种植园试验。１．２　同质种植园实验

芦苇同质园位于山东省青岛市山东大学青岛校区（３６．３６°Ｎ，１２０．６９°Ｅ），属于典型的温带季风性气候，年平均气温为１３．２℃，年平均降水量为６６６．７ｍｍ，夏季高温多雨，冬季寒冷干燥［２７］。我们共采集了宁夏平原１８

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

个种群的根状茎（每个种群选取１个个体），黄河三角洲１６个种群的根状茎（每个种群选取１个个体），于

３７５８　生　态　学　报　　　４１卷　

２０１８年８月，将这些芦苇根状茎分别种植在内径约２５ｃｍ、高约３３ｃｍ、容量约１６Ｌ的塑料桶中并放置于温室。于２０１９年５月，将所有发芽芦苇转移到露天同质种植园中，桶间距为１ｍ，确保不同种群芦苇间彼此不受影响。在芦苇生长期间，定期对芦苇进行除虫、施肥、人工除草，定期浇水，保持水分等资源充足。

质园芦苇的株高、基径、叶长、叶宽、叶厚、比叶面积等性状，并观察数出每桶中芦苇的株数。同时，分别记录桶中芦苇是否抽穗，用每个区域中抽穗的芦苇种群数量除以每个区域的总种群数量得到抽穗占比。除了以上形态性状外，在同质园中还测量了叶片气体交换、叶绿素荧光、叶片含氮量、叶片含磷量等生理指标。在晴朗无云的天气条件下，在早上９点到１１点，使用Ｌｉ⁃６８００便携式光合仪（ＬＩ⁃ＣＯＲ，美国）进行测量，设置其光强为１８００μｍｏｌｍ－２ｓ－１，二氧化碳浓度为４００μｍｏｌ／ｍｏｌ，温度为２５℃，湿度为５０％，空气流速为４００μｍｏｌ／ｓ，测得芦苇的净光合速率、蒸腾速率、气孔导度和胞间二氧化碳浓度。同时，使用ＰＡＭ⁃２５００便携式调制叶绿素荧光仪（ＷＡＬＺ，德国）对植物叶片测定瞬时荧光值（Ｆｔ）和最大荧光值（Ｆｍ’），然后通过计算得到ＰＳＩＩ的实际量子效率、非光化学淬灭、电子传递效率、ＰＳＩＩ的最大量子效率等荧光参数。光合和荧光指标测量完毕后，收获同质园中芦苇的地上部分，烘干测定每桶芦苇地上生物量。新鲜叶片的全氮含量使用凯氏定氮法获得。植物样品经过浓硫酸消解后，用钼锑抗比色法测定总磷含量。１．３　统计分析

测量的数据使用ＭｉｃｒｏｓｏｆｔＥｘｃｅｌ２０１６进行整理，使用ＩＢＭＳＰＳＳＳｔａｔｉｓｔｉｃｓＶ２１．０数据分析软件对所测的性状数据进行ｔ检验和Ｄｕｎｃａｎ检验，以比较在野外和同质园环境下的宁夏平原和黄河三角洲芦苇性状的差异。为检验性状的可遗传性与可塑性，对同质园中宁夏平原芦苇的株高、基径、叶长、叶宽、叶厚、比叶面积等性状分别与相对应的野外种群性状进行Ｐｅａｒｓｏｎ相关性和线性回归分析。为检验性状之间的协同与权衡关系，对所测的主要性状进行Ｓｐｅａｒｍａｎ相关性分析。绘图工具采用ＧｒａｐｈＰａｄＰｒｉｓｍ８．４．２和Ｒ软件。２　结果

２．１　不同种源间芦苇植物功能性状比较

在野外环境中，来自宁夏平原和黄河三角洲的芦苇在株高、基径、叶长、叶宽、叶厚性状上存在显著差异２）。来自宁夏平原的芦苇无论是在野外还是同质园环境下，其基径、叶长、叶宽都显著低于黄河三角洲芦苇（图１），但宁夏平原芦苇的抽穗占比要高于黄河三角洲芦苇（表１；表２）。在野外环境中，宁夏平原芦苇的株高、叶厚显著高于黄河三角洲芦苇（图１；表１）。在同质园环境中，两个种源间的芦苇在光合、呼吸等生理性状上均无差异（表３），但在叶片化学元素氮、磷含量和植株密度方面存在显著差异，来自宁夏平原的芦苇叶片含氮量和含磷量较低（表２），在植株数量（密度）上显著多于黄河三角洲芦苇（表２）。２．２　芦苇植物功能性状间的关联分析

芦苇的各个植物功能性状之间具有不同的关联，研究结果表明，无论是野外还是同质园环境中，比叶面积和基径、叶厚性状间都存在显著负相关关系（图２）。大部分芦苇都表现出叶长和叶宽两个性状的显著正相关，然而大多数芦苇的株高与基径、叶长、叶宽等性状间无显著相关性（图２）。２．３　野外和同质园中芦苇植物功能性状相关性分析

宁夏平原芦苇在野外和同质园两种环境下表现出了不同的表型，对于株高、叶长两个性状来说，两个环境下生长的芦苇之间存在显著正相关关系（表４），而在基径、叶宽、叶厚、比叶面积等性状上不存在显著的相关关系（表４）。３　讨论与分析

来源于不同环境的植物在株高、基径、叶片性状等植物功能性状上会产生差异，这是植物为适应不同环境

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

（图１；表１）；在同质园环境下，来自两个地区的芦苇在基径、叶长、叶宽、叶厚性状上也存在显著差异（图１；表

２０１９年９月，将每桶芦苇作为统计单位，从每桶芦苇中选取５棵健康植株，用与野外相同的方法测得同

　１０期　　　马香艳　等：基于野外调查和同质种植园实验的芦苇植物功能性状变异研究　３７５９

图１　宁夏平原和黄河三角洲芦苇在野外和同质园环境中植物功能性状的比较（数值为平均值±标准误差）

Ｆｉｇ．１　ＣｏｍｐａｒｉｓｏｎｏｆｐｌａｎｔｆｕｎｃｔｉｏｎａｌｔｒａｉｔｓｏｆＰｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓｉｎＮｉｎｇｘｉａＰｌａｉｎａｎｄｔｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａｉｎｔｈｅｆｉｅｌｄａｎｄｉｎｔｈｅｃｏｍｍｏｎｇａｒｄｅｎ（ｍｅａｎ±ｓｔａｎｄａｒｄｅｒｒｏｒ）

不同小写字母代表各组之间的显著性差异（Ｄｕｎｃａｎ检验，显著性水平Ｐ＜０．０５）

表１　野外宁夏平原和黄河三角洲芦苇的植物功能性状比较

Ｔａｂｌｅ１　ＣｏｍｐａｒｉｓｏｎｏｆｐｌａｎｔｆｕｎｃｔｉｏｎａｌｔｒａｉｔｓｏｆＰｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓｉｎＮｉｎｇｘｉａＰｌａｉｎａｎｄｔｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａｉｎｔｈｅｆｉｅｌｄ

性状Ｔｒａｉｔｓ

单株地上生物量Ｂｉｏｍａｓｓｐｅｒｐｌａｎｔ／ｇ株高Ｐｌａｎｔｈｅｉｇｈｔ／ｃｍ基径Ｂａｓａｌｄｉａｍｅｔｅｒ／ｃｍ叶长Ｌｅａｆｌｅｎｇｔｈ／ｃｍ叶宽Ｌｅａｆｗｉｄｔｈ／ｃｍ叶厚Ｌｅａｆｔｈｉｃｋｎｅｓｓ／ｃｍ比叶面积ＳＬＡ／（ｃｍ２／ｇ）抽穗占比Ｐｅｒｃｅｎｔａｇｅｏｆｓｐｉｋｅｓ

宁夏平原ＮｉｎｇｘｉａＰｌａｉｎ１１．５０±６．００１４０．００±４５０．４１±０．１５２２．１１±７．５６１．５１±０．５３０．０２±０．０５１３２．２９±１７．０１９６．６７％

黄河三角洲

ＴｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａ

１２．１０±５．１５１１８．００±４２０．４７±０．１３２４．７２±７．２６１．８４±０．５３０．０２±０．０６１３２．４３±１７．６３３．３３％

ｔ－０．９２１４．３００－４．０００－３．０００－５．３３０２．８４１－０．４２３－

Ｐ０．３５８＜０．００１＜０．００１０．００３＜０．００１０．００５０．６７３－

　　表格第２、３列中的数值为性状平均值±标准差，表格第４、５列为ｔ检验的结果

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

３７６０　生　态　学　报　　　

表２　同质园中宁夏平原和黄河三角洲芦苇的形态及化学元素含量性状比较

４１卷　

Ｔａｂｌｅ２　ＣｏｍｐａｒｉｓｏｎｏｆｍｏｒｐｈｏｌｏｇｉｃａｌａｎｄｃｈｅｍｉｃａｌｃｏｎｔｅｎｔｔｒａｉｔｓｏｆＰｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓｉｎＮｉｎｇｘｉａＰｌａｉｎａｎｄｔｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａｉｎｔｈｅｃｏｍｍｏｎｇａｒｄｅｎ性状Ｔｒａｉｔｓ总地上生物量Ｔｏｔａｌｂｉｏｍａｓｓ／ｇ株高Ｐｌａｎｔｈｅｉｇｈｔ／ｃｍ基径Ｂａｓａｌｄｉａｍｅｔｅｒ／ｃｍ叶长Ｌｅａｆｌｅｎｇｔｈ／ｃｍ叶宽Ｌｅａｆｗｉｄｔｈ／ｃｍ叶厚Ｌｅａｆｔｈｉｃｋｎｅｓｓ／ｃｍ比叶面积ＳＬＡ／（ｃｍ２／ｇ）株数Ｓｈｏｏｔｎｕｍｂｅｒ

叶含氮量Ｌｅａｆｎｉｔｒｏｇｅｎｃｏｎｔｅｎｔ

叶含磷量Ｌｅａｆｐｈｏｓｐｈｏｒｕｓｃｏｎｔｅｎｔ／（μｇ／ｍＬ）抽穗占比Ｐｅｒｃｅｎｔａｇｅｏｆｓｐｉｋｅｓ

宁夏平原ＮｉｎｇｘｉａＰｌａｉｎ３６．３１±８．１０８８．００±１９０．２２±０．６１１３．２２±３．９９０．６６±０．２１０．０１±０．００３１４７．２９±１７．１０１０３．００±４３０．３７±０．０６０．６３±０．１２８８．８９％

黄河三角洲

ＴｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａ

６１．４３±５４．３３９０．００±２８０．３４±０．６６１９．４１±４．６４１．０７±０．２８０．０２±０．００３１５１．１３±１５．０２６１．００±４１０．６７±０．３６０．７６±０．１７６２．５０％

ｔ１．８３００．５５２１２．５２１９．３９３１０．６００２．４００－２．８９０３．５３０２．５５０－０．７７８

Ｐ０．０８６０．５８６＜０．００１＜０．００１＜０．００１０．０２３０．００７０．００３０．０１６－０．４４２

　　表格第２、３列中的数值为性状平均值±标准差，表格第４、５列为ｔ检验结果

表３　同质园中宁夏平原和黄河三角洲芦苇的生理性状比较

Ｔａｂｌｅ３　ＣｏｍｐａｒｉｓｏｎｏｆｐｈｙｓｉｏｌｏｇｉｃａｌｔｒａｉｔｓｏｆＰｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓｉｎＮｉｎｇｘｉａＰｌａｉｎａｎｄｔｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａｉｎｔｈｅｃｏｍｍｏｎｇａｒｄｅｎ性状Ｔｒａｉｔｓ

净光合速率

Ｎｅｔｐｈｏｔｏｓｙｎｔｈｅｔｉｃｒａｔｅ／（μｍｏｌｍ－２ｓ－１）蒸腾速率

Ｔｒａｎｓｐｉｒａｔｉｏｎｒａｔｅ／（ｍｏｌｍ－２ｓ－１）

宁夏平原ＮｉｎｇｘｉａＰｌａｉｎ１５．５２±５．８００．０１±０．０２２６６．００±１０２１．１５±１．９４０．４９±０．０５０．７５±０．３７２９．０１±２．７１０．７７±０．０２

黄河三角洲

ＴｈｅＹｅｌｌｏｗＲｉｖｅｒＤｅｌｔａ

１５．５２±５．８００．０１±０．０１２９７．００±８８０．６８±０．６３０．５０±０．０６０．７０±０．２８２９．３２±３．６００．７５±０．０３

ｔ－０．０１７－０．６９９０．９４７－０．９２３０．３４４－０．５３００．３３５－１．５９０

Ｐ０．９８７０．４８９０．３５１０．３６３０．７３３０．６０００．７４００．１２３

胞间二氧化碳浓度

Ｉｎｔｅｒｃｅｌｌｕｌａｒｃａｒｂｏｎｄｉｏｘｉｄｅ／（μｍｏｌ／ｍｏｌ）气孔导度

Ｓｔｏｍａｔａｌｃｏｎｄｕｃｔａｎｃｅ／（ｍｏｌｍ－２ｓ－１）ＰＳＩＩ的实际量子效率

ＴｈｅａｃｔｕａｌｑｕａｎｔｕｍｅｆｆｉｃｉｅｎｃｙｏｆＰＳＩＩ非光化学淬灭

Ｎｏｎ－ＰｈｏｔｏｃｈｅｍｉｃａｌＱｕｅｎｃｈｉｎｇ

电子传递效率

Ｅｌｅｃｔｒｏｎｔｒａｎｓｆｅｒｅｆｆｉｃｉｅｎｃｙ／（μｍｏｌｍ－２ｓ－１）ＰＳＩＩ的最大量子效率

ＭａｘｉｍｕｍｑｕａｎｔｕｍｅｆｆｉｃｉｅｎｃｙｏｆＰＳＩＩ

　　表格第２、３列为性状平均值±标准差，表格第４、５列为ｔ检验结果

表４　宁夏平原芦苇在同质园和野外的植物功能性状Ｐｅａｒｓｏｎ相关性和线性回归分析

Ｔａｂｌｅ４　ＰｅａｒｓｏｎｃｏｒｒｅｌａｔｉｏｎａｎｄｌｉｎｅａｒｒｅｇｒｅｓｓｉｏｎａｎａｌｙｓｉｓｏｆｐｌａｎｔｆｕｎｃｔｉｏｎａｌｔｒａｉｔｓｏｆＰｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓｉｎｔｈｅｃｏｍｍｏｎｇａｒｄｅｎａｎｄｆｉｅｌｄｏｆＮｉｎｇｘｉａＰｌａｉｎ性状Ｔｒａｉｔｓ

株高Ｐｌａｎｔｈｅｉｇｈｔ／ｃｍ基径Ｂａｓａｌｄｉａｍｅｔｅｒ／ｃｍ叶长Ｌｅａｆｌｅｎｇｔｈ／ｃｍ叶宽Ｌｅａｆｗｉｄｔｈ／ｃｍ叶厚Ｌｅａｆｔｈｉｃｋｎｅｓｓ／ｃｍ　　∗Ｐ＜０．０５，∗∗Ｐ＜０．０１比叶面积ＳＬＡ／（ｃｍ２／ｇ）

Ｐｅａｒｓｏｎ相关性系数（ｒ）

Ｐｅａｒｓｏｎｃｏｒｒｅｌａｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔ（ｒ）

∗

０．６３８∗

线性回归系数（Ｒ２）

Ｌｉｎｅａｒｒｅｇｒｅｓｓｉｏｎｃｏｅｆｆｉｃｉｅｎｔ（Ｒ２）

∗

０．４０７∗

０．４０９０．５６７∗０．３２６０．４４６－０．１１２

０．１９３０．１０６０．１９８０．１２４

０．３２２∗

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　１０期　　　马香艳　等：基于野外调查和同质种植园实验的芦苇植物功能性状变异研究　３７６１

图２　芦苇植物功能性状间Ｓｐｅａｒｍａｎ相关性分析

Ｆｉｇ．２　ＳｐｅａｒｍａｎｃｏｒｒｅｌａｔｉｏｎａｎａｌｙｓｉｓｏｆｆｕｎｃｔｉｏｎａｌｔｒａｉｔｓｏｆＰｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓ

∗代表Ｐ＜０．０５，∗∗代表Ｐ＜０．０１，∗∗∗代表Ｐ＜０．００１

而采取的生存策略［４］。来自宁夏平原和黄河三角洲的芦苇在野外和同质园中都在株高、基径、叶长、叶宽、叶厚等形态性状方面存在显著性差异，在叶片含氮量和含磷量等化学成分含量等性状上也有显著不同，但是在叶片气体交换、叶绿素荧光等生理性状和比叶面积性状方面并无显著差异。芦苇性状间也存在协变现象，通过性状间相互权衡、相互影响来使植物达到资源利用最优化［２８］，维持其正常的生长。当生存环境发生很大的变化时，宁夏平原芦苇的株高和叶长还是主要受到遗传因素的影响，基径、叶宽、叶厚、比叶面积等性状受到环境的影响较大，遗传和环境共同影响着植物的功能性状。在未来全球变化的大背景下，植物可能会不断通过遗传变异和表型可塑性两种机制在新环境中生存繁殖和定居［３］。

在不同的气候等环境条件下，芦苇种群会发生本地适应，产生稳定的表型分化。中国东西部不同来源植

物间的稳定表型分化的形成可能和降水有关，对内蒙古草原与荒漠区中４种一年生植物进行模拟研究时就发现，随着有效降水梯度的增加，植物的生长期、植物高度、叶长、生物量等均随之增加［２９］。在本研究中，黄河三角洲靠近海岸，年平均降水量为５５１．６ｍｍ，降水丰沛，地质复杂且异质性强［２６］，而宁夏平原处于内陆地区，年平均降水量为２６２．３—３６３．２ｍｍ，降水不足，气候干燥［２５］，在正常年份，黄河三角洲年平均降水量约为宁夏平原的两倍，水分可能成为限制宁夏平原芦苇叶片性状的主要因素之一。对于草本植物来说，在抽穗、开花期，降水量越多，物候发生时间越提前［３０］。但是有的研究却发现，土壤水分的减少可能导致植物快速完成营养生长而提前进入生殖生长阶段，长期干旱造成了植物物候的提前［２９，３１］。在同质园中，我们发现宁夏平原芦苇抽穗时间早于黄河三角洲芦苇，宁夏平原芦苇抽穗比例也大于黄河三角洲芦苇（表２），这与在野外观察到的宁夏

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

３７６２　生　态　学　报　　　４１卷　

平原芦苇抽穗植株占比远远大于黄河三角洲芦苇相一致（表１），长期对不同水分条件的适应形成了两个地区芦苇由遗传控制的稳定表型分化。在本研究中，野外宁夏平原芦苇的株高显著高于黄河三角洲芦苇（图１），可能是由于黄河三角洲盐分含量过高，抑制了芦苇的生长，这与在加盐处理的控制实验中得到的结果相一致，加盐处理下，河口型和内陆型两种生态型芦苇的株高都呈现出下降趋势，都遭受到了盐胁迫，株高下降也是芦苇针对高盐生境的一种表型适应策略［３２］。两个地区芦苇的株高在野外条件下也存在显著差异，但是在同质园环境下，这种株高差异消失（图１），表明了环境对植物有着很大的影响，表型可塑性也是芦苇产生表型分化以适应环境变化的重要机制。因此，芦苇在不同气候环境下的植物功能性状变异受遗传变异和表型可塑性共同影响。

其中，比叶面积是最重要的性状之一，能反映植物对不同生境的适应特征，使其成为植物功能性状研究中的重要指标［３４］。许多研究认为，在资源丰富的环境中生长的物种往往有较高的比叶面积，而较低的比叶面积和大境中的芦苇（图１），因为在同质园中，芦苇生长环境相对优越，与野外生境相比，水热光资源较为充裕，较大的比叶面积有利于其充分利用光照，快速获取资源。在对木本植物的研究中也证明了这点，研究表明，来自湿润环境的栎树具有高的比叶面积和低的叶片厚度，同时具有较高的光合速率，但生长速率因时期而异，是一种提高竞争能力的资源获取策略［３６］。黄河三角洲芦苇的叶长、叶宽等叶片性状都显著大于宁夏平原，叶片中的氮磷含量也较高，但是黄河三角洲芦苇的密度却远远小于宁夏平原芦苇（表２），这反映了植物适应环境的不同策略，黄河三角洲芦苇偏向于对资源快速获取的方式来适应环境，而宁夏平原芦苇则通过增加株数、提高繁殖力来维持生长，都是植物长期适应环境形成的结果［３７］。

植物对环境的适应并不是通过单一性状的改变而实现的，植物的各个性状之间是相互关联、相互影响植物的叶片性状与植物利用资源的能力密切相关，而且也能反映出植物应对不同环境的生存策略［３３］。

的叶片厚度有利于植物适应资源匮乏的环境［３５］。在本研究中，同质园中的芦苇比叶面积要显著高于野外环

的［２８］。现在研究的对象也从植物叶片拓展到了根、茎、花、种子等其他器官及植物整体性状［１２］，植物的不同器官（如叶、茎和花）对资源环境变化的响应往往是相互作用的，以便其进行有效的生长和繁殖，例如在武夷山５９种木本植物的树皮、茎干、叶片功能性状特征及关联研究中发现，树皮厚度与茎干含水率、茎干氮含量呈极显著负相关，树皮碳含量与茎干碳含量呈显著正相关［３８］，功能性状间的关联性有利于植物快速调整生存策即植物比叶面积越大，基径越小，芦苇通过比叶面积和基径两个性状间的协变来共同适应环境，体现性状间的相互权衡。这种相关性不仅体现在植物不同器官之间，同一个器官内也具有很强的协变性［３９］，例如叶片形态大小和叶片表面温度间存在某种关联，在高温干旱区域，小的叶片利于降低叶片表面温度，这是植物对特殊生境的一种适应策略［４０］。这种性状协变性的大小还可能会限制表型可塑性，有研究表明不同的环境条件会改变或导致性状协变性的变化，从而使表型可塑性得到改变，但是在研究中这种协变性相对较为稳定，所以对表型可塑性影响较小［４１］。

同种植物的植物功能性状变异不仅受到表型可塑性的影响，也受到遗传因素的控制，二者可以同时作为略以适应环境。在本研究中，比叶面积和基径两个植物性状是协同变化的，两者间呈显著负相关关系（图２），

某些植物耐受性强、生态型丰富的重要原因。当植物的生存条件突然改变，植物的表型可塑性导致其产生与原来不同的功能性状来适应新环境，例如本研究中芦苇在同质园环境中比叶面积大大增加。这种机制能使植物快速的改变生存策略，在最短的时间内获得最高的资源利用效率。表型可塑性也可能是物种避免灭绝或迁徙的一种策略，通过增强个体的适应能力来扩大其对多种环境更广泛的耐受性［２２］。植物如果长期生存在某一特定环境中，遗传基因可能发生改变，导致植物表型产生稳定的变异，这是植物长期适应其生存环境的结果，就像本研究中宁夏平原芦苇的株高和叶长受到遗传因素的影响，当环境改变时仍然保持相同的生长策略，在野外和同质园环境下都具有显著相关性。当然，在全球气候变化的大背景下，更多不同来源的芦苇对多重环境因子的响应机制还需要得到持续的关注和探究［４２］。

本论文探究了黄河上下游芦苇在野外和同质种植园中的植物功能性状变异规律，揭示了遗传分化和表型

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　１０期　　　马香艳　等：基于野外调查和同质种植园实验的芦苇植物功能性状变异研究　３７６３

可塑性在芦苇环境适应中的重要作用，对进一步研究广布植物物种对全球环境变化的响应和适应机制具有重要意义，也为芦苇湿地的保护、恢复和管理提供了科学依据。

参考文献（Ｒｅｆｅｒｅｎｃｅｓ）：

［１］　唐青青，黄永涛，丁易，臧润国．亚热带常绿落叶阔叶混交林植物功能性状的种间和种内变异．生物多样性，２０１６，２４（３）：２６２⁃２７０．

２０２０，４２（１）：１⁃６．

［２］　ＸｕＷＭ，ＴｏｍｌｉｎｓｏｎＫＷ，ＬｉＪ．Ｓｔｒｏｎｇｉｎｔｒａｓｐｅｃｉｆｉｃｔｒａｉｔｖａｒｉａｔｉｏｎｉｎａｔｒｏｐｉｃａｌｄｏｍｉｎａｎｔｔｒｅｅｓｐｅｃｉｅｓａｌｏｎｇａｎｅｌｅｖａｔｉｏｎａｌｇｒａｄｉｅｎｔ．ＰｌａｎｔＤｉｖｅｒｓｉｔｙ，［３］　杨雪，申俊芳，赵念席，高玉葆．不同基因型羊草数量性状的可塑性及遗传分化．植物生态学报，２０１７，４１（３）：３５９⁃３６８．

［４］　ＡｇｒａｗａｌＡＡ．Ａｓｃａｌｅ⁃ｄｅｐｅｎｄｅｎｔｆｒａｍｅｗｏｒｋｆｏｒｔｒａｄｅ⁃ｏｆｆｓ，ｓｙｎｄｒｏｍｅｓ，ａｎｄｓｐｅｃｉａｌｉｚａｔｉｏｎｉｎｏｒｇａｎｉｓｍａｌｂｉｏｌｏｇｙ．Ｅｃｏｌｏｇｙ，２０２０，１０１（２）：ｅ０２９２４．

８６３⁃８７６．

［５］　杨继鸿，李亚楠，卜海燕，张世挺，齐威．青藏高原东缘常见阔叶木本植物叶片性状对环境因子的响应．植物生态学报，２０１９，４３（１０）：［６］　ＬｉａｏＨＸ，ＬｉＤＪ，ＺｈｏｕＴ，ＨｕａｎｇＢ，ＺｈａｎｇＨＪ，ＣｈｅｎＢＭ，ＰｅｎｇＳＬ．Ｔｈｅｒｏｌｅｏｆｆｕｎｃｔｉｏｎａｌｓｔｒａｔｅｇｉｅｓｉｎｇｌｏｂａｌｐｌａｎｔｄｉｓｔｒｉｂｕｔｉｏｎ．Ｅｃｏｇｒａｐｈｙ，［７］　张增可，郑心炫，林华贞，林欣，黄柳菁．海岛植物不同演替阶段植物功能性状与环境因子的变化规律．生态学报，２０１９，３９（１０）：［８］　ＭａｅｓＳＬ，ＰｅｒｒｉｎｇＭＰ，ＤｅｐａｕｗＬ，Ｂｅｒｎｈａｒｄｔ⁃ＲöｍｅｒｍａｎｎＭ，ＢｌｏｎｄｅｅｌＨ，ＢｒūｍｅｌｉｓＧ，ＢｒｕｎｅｔＪ，ＤｅｃｏｃｑＧ，ｄｅｎＯｕｄｅｎＪ，ＧｏｖａｅｒｔＳ，Ｈäｒｄｔｌｅ

′Ｍ，ＬａｎｄｕｙｔＤ，ＭáｌｉšＦ，ＶａｎｎｅｓｔｅＴ，ＷｕｌｆＭ，ＶｅｒｈｅｙｅｎＫ．Ｗ，ＨéｄｌＲ，ＨｅｉｎｋｅｎＴ，ＨｅｉｎｒｉｃｈｓＳ，ＨｅｒｔｚｏｇＬ，ＪａｒｏｓｚｅｗｉｃｚＢ，ＫｉｒｂｙＫ，Ｋｏｐｅｃｋｙ

ＰｌａｎｔｆｕｎｃｔｉｏｎａｌｔｒａｉｔｒｅｓｐｏｎｓｅｔｏｅｎｖｉｒｏｎｍｅｎｔａｌｄｒｉｖｅｒｓａｃｒｏｓｓＥｕｒｏｐｅａｎｔｅｍｐｅｒａｔｅｆｏｒｅｓｔｕｎｄｅｒｓｔｏｒｅｙｃｏｍｍｕｎｉｔｉｅｓ．ＰｌａｎｔＢｉｏｌｏｇｙ，２０２０，２２（３）：３７４９⁃３７５８．

２０２１，４４（４）：４９３⁃５０３．

［９］　ＲｅｉｃｈＰＢ，ＯｌｅｋｓｙｎＪ．ＧｌｏｂａｌｐａｔｔｅｒｎｓｏｆｐｌａｎｔｌｅａｆＮａｎｄＰｉｎｒｅｌａｔｉｏｎｔｏｔｅｍｐｅｒａｔｕｒｅａｎｄｌａｔｉｔｕｄｅ．ＰｒｏｃｅｅｄｉｎｇｓｏｆｔｈｅＮａｔｉｏｎａｌＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ［１０］　董廷发，冯玉龙，类延宝，张丽坤．干旱和湿润生境中主要优势树种叶片功能性状的比较．生态学杂志，２０１２，３１（５）：１０４３⁃１０４９．

Ｏｅｃｏｌｏｇｉｃａ，２００６，３０（２）：２１２⁃２２２．

ｏｆｔｈｅＵｎｉｔｅｄＳｔａｔｅｓｏｆＡｍｅｒｉｃａ，２００４，１０１（３０）：１１００１⁃１１００６．

４１０⁃４２４．

［１１］　ＳｕｎＳＣ，ＪｉｎＤＭ，ＬｉＲＪ．ＬｅａｆｅｍｅｒｇｅｎｃｅｉｎｒｅｌａｔｉｏｎｔｏｌｅａｆｔｒａｉｔｓｉｎｔｅｍｐｅｒａｔｅｗｏｏｄｙｓｐｅｃｉｅｓｉｎＥａｓｔ⁃ＣｈｉｎｅｓｅＱｕｅｒｃｕｓｆａｂｒｉｆｏｒｅｓｔｓ．Ａｃｔａ［１２］　何芸雨，郭水良，王喆．植物功能性状权衡关系的研究进展．植物生态学报，２０１９，４３（１２）：１０２１⁃１０３５．

［１３］　ＤｅＢａｔｔｉｓｔｉＤ，ＦｏｗｌｅｒＭＳ，ＪｅｎｋｉｎｓＳＲ，ＳｋｏｖＭＷ，ＢｏｕｍａＴＪ，ＮｅｙｌａｎｄＰＪ，ＧｒｉｆｆｉｎＪＮ．Ｍｕｌｔｉｐｌｅｔｒａｉｔｄｉｍｅｎｓｉｏｎｓｍｅｄｉａｔｅｓｔｒｅｓｓｇｒａｄｉｅｎｔｅｆｆｅｃｔｓ［１４］　钟巧连，刘立斌，许鑫，杨勇，郭银明，许海洋，蔡先立，倪健．黔中喀斯特木本植物功能性状变异及其适应策略．植物生态学报，２０１８，［１５］　ＹａｎｇＹ，ＬｉｕＭ，ＰａｎＹＦ，ＨｕａｎｇＨＹ，ＰａｎＸＹ，ＳｏｓａＡ，ＨｏｕＹＰ，ＺｈｕＺＣ，ＬｉＢ．Ｒａｐｉｄｅｖｏｌｕｔｉｏｎｏｆｌａｔｉｔｕｄｉｎａｌｃｌｉｎｅｓｉｎｇｒｏｗｔｈａｎｄｄｅｆｅｎｃｅｏｆ［１６］　ＰｉｔｈｅｒＪ．Ｃｌｉｍａｔｅｔｏｌｅｒａｎｃｅａｎｄｉｎｔｅｒｓｐｅｃｉｆｉｃｖａｒｉａｔｉｏｎｉｎｇｅｏｇｒａｐｈｉｃｒａｎｇｅｓｉｚｅ．ＰｒｏｃｅｅｄｉｎｇｓｏｆｔｈｅＲｏｙａｌＳｏｃｉｅｔｙＢ：ＢｉｏｌｏｇｉｃａｌＳｃｉｅｎｃｅｓ，２００３，２７０［１７］　ＰａｃｋｅｒＪＧ，ＭｅｙｅｒｓｏｎＬＡ，ＳｋáｌｏｖáＨ，ＰｙšｅｋＰ，ＫｕｅｆｆｅｒＣ．ＢｉｏｌｏｇｉｃａｌＦｌｏｒａｏｆｔｈｅＢｒｉｔｉｓｈＩｓｌｅｓ：Ｐｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓ．ＪｏｕｒｎａｌｏｆＥｃｏｌｏｇｙ，２０１７，［１８］　ＣｈａｍｂｅｒｓＲＭ，ＭｅｙｅｒｓｏｎＬＡ，ＳａｌｔｏｎｓｔａｌｌＫ．ＥｘｐａｎｓｉｏｎｏｆＰｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓｉｎｔｏｔｉｄａｌｗｅｔｌａｎｄｓｏｆＮｏｒｔｈＡｍｅｒｉｃａ．ＡｑｕａｔｉｃＢｏｔａｎｙ，１９９９，６４［１９］　庄瑶，孙一香，王中生，杨琳璐，邓自发，姚志刚，安树青．芦苇生态型研究进展．生态学报，２０１０，３０（８）：２１７３⁃２１８１．

（３／４）：２６１⁃２７３．１０５（４）：１１２３⁃１１６２．（１５１４）：４７５⁃４８１．

ａｎｉｎｖａｓｉｖｅｗｅｅｄ．ＮｅｗＰｈｙｔｏｌｏｇｉｓｔ，２０２１，２３０（２）：８４５⁃８５６．４２（５）：５６２⁃５７２．

ｏｎｐｌａｎｔｂｉｏｍａｓｓａｌｌｏｃａｔｉｏｎ，ｗｉｔｈｉｍｐｌｉｃａｔｉｏｎｓｆｏｒｃｏａｓｔａｌｅｃｏｓｙｓｔｅｍｓｅｒｖｉｃｅｓ．ＪｏｕｒｎａｌｏｆＥｃｏｌｏｇｙ，２０２０，１０８（４）：１２２７⁃１２４０．

［２０］　朱学艺，王锁民，张承烈．河西走廊不同生态型芦苇对干旱和盐渍胁迫的响应调节．植物生理学通讯，２００３，３９（４）：３７１⁃３７６．［２１］　张淑萍，王仁卿，张治国，郭卫华，刘建，宋百敏．黄河下游湿地芦苇形态变异研究．植物生态学报，２００３，２７（１）：７８⁃８５．

［２２］　ＲｅｎＬＪ，ＧｕｏＸ，ＬｉｕＳＮ，ＹｕＴ，ＧｕｏＷＨ，ＷａｎｇＲＱ，ＹｅＳＹ，ＬａｍｂｅｒｔｉｎｉＣ，ＢｒｉｘＨ，ＥｌｌｅｒＦ．ＩｎｔｒａｓｐｅｃｉｆｉｃｖａｒｉａｔｉｏｎｉｎＰｈｒａｇｍｉｔｅｓａｕｓｔｒａｌｉｓ：［２３］　陈琳，张俪文，刘子亭，路峰，冯光海，颜坤，韩广轩．黄河三角洲河滩与潮滩芦苇对盐胁迫的生理生态响应．生态学报，２０２０，４０（６）：［２４］　李帅，马文超，顾艳文，魏虹，彭月，李昌晓．宁夏黄河流域景观破碎化时空变化特征．生态学报，２０１６，３６（１１）：３３１２⁃３３２０．［２６］　宗秀影，刘高焕，乔玉良，林松．黄河三角洲湿地景观格局动态变化分析．地球信息科学学报，２００９，１１（１）：９１⁃９７．［２７］　夏子瑶，万夫敬．１９８７⁃２０１６年青岛市气候变化特征分析．现代农业科技，２０１８，（４）：１７１⁃１７１，１７３⁃１７３．

２０９０⁃２０９８．

Ｃｌｉｎａｌａｄａｐｔｉｏｎｏｆｆｕｎｃｔｉｏｎａｌｔｒａｉｔｓａｎｄｐｈｅｎｏｔｙｐｉｃｐｌａｓｔｉｃｉｔｙｖａｒｙｗｉｔｈｌａｔｉｔｕｄｅｏｆｏｒｉｇｉｎ．ＪｏｕｒｎａｌｏｆＥｃｏｌｏｇｙ，２０２０，１０８（６）：２５３１⁃２５４３．

［２５］　王豫，赵小艳，李艳春，王咏薇，朱婷婷．宁夏平原湿地面积动态演变对局地气候效应的影响．生态环境学报，２０１８，２７（７）：１２５１⁃１２５９．

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

３７６４　生　态　学　报　　　４１卷　

［２８］　ＰｉｇｌｉｕｃｃｉＭ．Ｐｈｅｎｏｔｙｐｉｃｉｎｔｅｇｒａｔｉｏｎ：ｓｔｕｄｙｉｎｇｔｈｅｅｃｏｌｏｇｙａｎｄｅｖｏｌｕｔｉｏｎｏｆｃｏｍｐｌｅｘｐｈｅｎｏｔｙｐｅｓ．ＥｃｏｌｏｇｙＬｅｔｔｅｒｓ，２００３，６（３）：２６５⁃２７２．［３０］　符瑜，潘学标．草本植物物候及其物候模拟模型的研究进展．中国农业气象，２０１１，３２（３）：３１９⁃３２５．

ＳｏｃｉｅｔｙｏｆＪａｐａｎ．Ｓｅｒ．Ⅱ，２００７，８５（３）：３５９⁃３６７．性．生态学报，２０１９，３９（２）：５４２⁃５４９．

［２９］　闫建成，梁存柱，付晓玥，王炜，王立新，贾成朕．草原与荒漠一年生植物性状对降水变化的响应．草业学报，２０１３，２２（１）：６８⁃７６．［３１］　ＳｈｉｎｏｄａＭ，ＩｔｏＳ，ＮａｃｈｉｎｓｈｏｎｈｏｒＧＵ，ＥｒｄｅｎｅｔｓｅｔｓｅｇＤ．ＰｈｅｎｏｌｏｇｙｏｆＭｏｎｇｏｌｉａｎｇｒａｓｓｌａｎｄｓａｎｄｍｏｉｓｔｕｒｅｃｏｎｄｉｔｉｏｎｓ．ＪｏｕｒｎａｌｏｆｔｈｅＭｅｔｅｏｒｏｌｏｇｉｃａｌ［３２］　许振伟，宋慧佳，李明燕，张廷靖，郭霄，ＦｒａｎｚｉｓｋａＥ，ＨａｎｓＢ，杜宁，侯文轩，郭卫华．不同生态型芦苇种群对盐胁迫的生长和光合特［３３］　ＰｉｅｒｃｅＳ，ＢｒｕｓａＧ，ＶａｇｇｅＩ，ＣｅｒａｂｏｌｉｎｉＢＥＬ．ＡｌｌｏｃａｔｉｎｇＣＳＲｐｌａｎｔｆｕｎｃｔｉｏｎａｌｔｙｐｅｓ：ｔｈｅｕｓｅｏｆｌｅａｆｅｃｏｎｏｍｉｃｓａｎｄｓｉｚｅｔｒａｉｔｓｔｏｃｌａｓｓｉｆｙｗｏｏｄｙａｎｄ［３４］　梁士楚，刘润红，荣春艳，常斌，姜勇．漓江河岸带植物功能性状变异与关联．植物生态学报，２０１９，４３（１）：１６⁃２６．

ｈｅｒｂａｃｅｏｕｓｖａｓｃｕｌａｒｐｌａｎｔｓ．ＦｕｎｃｔｉｏｎａｌＥｃｏｌｏｇｙ，２０１３，２７（４）：１００２⁃１０１０．

［３５］　ＷｒｉｇｈｔＩＪ，ＲｅｉｃｈＰＢ，ＷｅｓｔｏｂｙＭ，ＡｃｋｅｒｌｙＤＤ，ＢａｒｕｃｈＺ，ＢｏｎｇｅｒｓＦ，Ｃａｖｅｎｄｅｒ⁃ＢａｒｅｓＪ，ＣｈａｐｉｎＴ，ＣｏｒｎｅｌｉｓｓｅｎＪＨＣ，ＤｉｅｍｅｒＭ，ＦｌｅｘａｓＪ，

ＧａｒｎｉｅｒＥ，ＧｒｏｏｍＰＫ，ＧｕｌｉａｓＪ，ＨｉｋｏｓａｋａＫ，ＬａｍｏｎｔＢＢ，ＬｅｅＴ，ＬｅｅＷ，ＬｕｓｋＣ，ＭｉｄｇｌｅｙＪＪ，ＮａｖａｓＭ⁃Ｌ，ＮｉｉｎｅｍｅｔｓÜ，ＯｌｅｋｓｙｎＪ，Ｏｓａｄａｓｐｅｃｔｒｕｍ．Ｎａｔｕｒｅ，２００４，４２８（６９８５）：８２１⁃８２７．

Ｎ，ＰｏｏｒｔｅｒＨ，ＰｏｏｔＰ，ＰｒｉｏｒＬ，ＰｙａｎｋｏｖＶＩ，ＲｏｕｍｅｔＣ，ＴｈｏｍａｓＳＣ，ＴｊｏｅｌｋｅｒＭＧ，ＶｅｎｅｋｌａａｓＥＪ，ＶｉｌｌａｒＲ．Ｔｈｅｗｏｒｌｄｗｉｄｅｌｅａｆｅｃｏｎｏｍｉｃｓ

［３６］　Ｒａｍíｒｅｚ⁃ＶａｌｉｅｎｔｅＪＡ，ＬóｐｅｚＲ，ＨｉｐｐＡＬ，ＡｒａｎｄａＩ．Ｃｏｒｒｅｌａｔｅｄｅｖｏｌｕｔｉｏｎｏｆｍｏｒｐｈｏｌｏｇｙ，ｇａｓｅｘｃｈａｎｇｅ，ｇｒｏｗｔｈｒａｔｅｓａｎｄｈｙｄｒａｕｌｉｃｓａｓａｒｅｓｐｏｎｓｅ［３７］　ＤａｙｒｅｌｌＲＬＣ，ＡｒｒｕｄａＡＪ，ＰｉｅｒｃｅＳ，ＮｅｇｒｅｉｒｏｓＤ，ＭｅｙｅｒＰＢ，ＬａｍｂｅｒｓＨ，ＳｉｌｖｅｉｒａＦＡＯ．Ｏｎｔｏｇｅｎｅｔｉｃｓｈｉｆｔｓｉｎｐｌａｎｔｅｃｏｌｏｇｉｃａｌｓｔｒａｔｅｇｉｅｓ．［３８］　范瑞瑞．武夷山５９种木本植物树皮、茎干、叶片功能性状特征及其关联研究［Ｄ］．福州：福建师范大学，２０１８．

ＦｕｎｃｔｉｏｎａｌＥｃｏｌｏｇｙ，２０１８，３２（１２）：２７３０⁃２７４１．

ｔｏｐｒｅｃｉｐｉｔａｔｉｏｎａｎｄｔｅｍｐｅｒａｔｕｒｅｒｅｇｉｍｅｓｉｎｏａｋｓ（Ｑｕｅｒｃｕｓ）．ＮｅｗＰｈｙｔｏｌｏｇｉｓｔ，２０２０，２２７（３）：７９４⁃８０９．

［３９］　靳莎，闫淑君，黄柳菁，陈莹，马雯雯，王云霄，王喆．植物叶功能性状间的权衡研究进展．四川林业科技，２０１９，４０（５）：９６⁃１０３．［４０］　李永华，李臻，辛智鸣，刘明虎，李艳丽，郝玉光．形态变化对叶片表面温度的影响．植物生态学报，２０１８，４２（２）：２０２⁃２０８．［４１］　ＧｉａｎｏｌｉＥ，Ｐａｌａｃｉｏ⁃ＬóｐｅｚＫ．Ｐｈｅｎｏｔｙｐｉｃｉｎｔｅｇｒａｔｉｏｎｍａｙｃｏｎｓｔｒａｉｎｐｈｅｎｏｔｙｐｉｃｐｌａｓｔｉｃｉｔｙｉｎｐｌａｎｔｓ．Ｏｉｋｏｓ，２００９，１１８（１２）：１９２４⁃１９２８．

ｄｏｍｉｎａｎｃｅｏｆａｎｅｘｏｔｉｃｐｌａｎｔｗｉｔｈａｎｏｖｅｌｔｒａｄｅ⁃ｏｆｆｐａｔｔｅｒｎ．ＪｏｕｒｎａｌｏｆＥｃｏｌｏｇｙ，２０２０，１０８（５）：１９８９⁃１９９９．

［４２］　ＱｉｕＳＹ，ＬｉｕＳＳ，ＷｅｉＳＪ，ＣｕｉＸＨ，ＮｉｅＭ，ＨｕａｎｇＪＸ，ＨｅＱ，ＪｕＲＴ，ＬｉＢ．Ｃｈａｎｇｅｓｉｎｍｕｌｔｉｐｌｅｅｎｖｉｒｏｎｍｅｎｔａｌｆａｃｔｏｒｓａｄｄｉｔｉｖｅｌｙｅｎｈａｎｃｅｔｈｅ

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ