植物功能性状是指植物适应环境变化所形成的形态、生理和物候特征, 能反映植物对不同环境的响应和适应机制^[1-2]。植物功能性状不仅影响生态系统功能, 还能将环境、植物和生态系统过程与功能密切联系起来^[3]。叶片作为与外界环境进行物质交换的重要器官, 其功能性状易测量、数量多、可塑性较大, 特别是结构型性状(比叶面积、叶厚度、叶干物质含量等), 是植物生长需求对空间异质性或环境变化适应进化的结果^[4-6]。近年来, 国内外学者开展了典型植物群落叶结构型功能性状与光照、气候、温度、降水、地形地貌、土壤养分等差异化环境因子间的相关性研究^[7-9]。如, 前人研究表明, 植物在林下充足光照环境下倾向于资源获取策略, 在弱光环境中生存策略由资源获取策略转变为资源保守策略^[7]。坡向、坡位、海拔等地形因子的变化会间接影响温度、降水、光照等因子产生再分配, 导致气候环境的空间异质性, 影响植物的生长发育并对植物功能性状的变化产生间接作用^[9]。因此, 致力于揭示植物叶功能性状与环境因子的关系, 有助于更好理解植物应对环境变化的适应策略及分配利用资源能力。

海岛由于地理位置隔离, 成为独立的生态系统, 相比内陆生态系统, 其更具敏感性和脆弱性。地理隔离显著影响植物的表型性状^[10], 如, 黄柳菁等调查了大陆性海岛平潭岛的种子植物, 并分析了其植物功能性状特征, 总结出海岛植物结构功能性状与内陆的区别, 强调海岛植物为了适应强风、干旱、贫瘠的环境, 呈现出以草本植物为主、小型叶植物占比大、叶被毛植物居多等结构功能性状特征^[11]。国内一些典型海岛(海南岛和平潭岛)不同演替阶段植物功能性状与环境因子的变化规律已有研究^[12-13]。如, 张增可等研究了平潭岛4个不同演替阶段群落植物功能性状对环境因子的响应, 结果表明, 演替早期植物主要分布在土壤pH值、容重高的贫瘠环境, 拥有较高比叶面积(SLA)、叶磷含量(LPC)、茎氮含量、茎磷含量的性状组合；演替后期植物主要分布在土壤养分和水分高的肥沃环境, 拥有较高的茎组织密度、叶干物质含量(LDMC)、叶碳含量(LCC)、叶氮含量(LNC)的性状组合。土壤有机质和全氮是影响海岛植物功能性状的主要土壤因子, 坡位和坡度是影响海岛植物功能性状的主要地形因子^[14]。海岛典型植物群落优势种叶功能性状与环境因子的相关性也有研究, 如, 琅岐岛4种优势植物叶功能性状主要受坡度、海拔、土壤电导率、土壤含水量等环境因子的影响^[15]。综上所述, 关于海岛不同典型植物群落优势种叶功能性状与环境因子的相关已有研究, 但典型植物群落同一优势种适应不同异化生境叶功能性状的变动还没有研究。

滨柃(Eurya emarginata), 山茶科(Theaceae)柃属(Eurya)常绿灌木。在中国呈地带性间断分布于浙江、福建和台湾等东南沿海及海岛, 且能形成以滨柃为优势种的滨海群落。由于滨柃群落分布特性, 关于滨柃群落结构特征已有研究。聂森研究发现, 滨柃群落物种丰富度和多样性与海拔呈负相关^[16]。李加庆对群落优势种滨柃的平均株高、平均地径与环境因子的相关性研究表明, 滨柃的平均株高与pH值、速效钾、全盐量、海拔呈负相关, 滨柃的平均地径与全钾呈负相关^[17]。滨柃多生于滨海山坡灌丛中及基岩海岛、海岸边岩石缝中, 在海岛恶劣生境中表现出较强抗逆性和适应性, 因此, 关于滨柃的抗逆性等方面国内外展开了较多研究^[18-19], 如, 滨柃在150 mmol/L的盐胁迫环境下能正常生长, 可作为海岛环境修复的耐盐优势树种^[18]。滨柃的存活率和生物量与干旱胁迫程度呈负相关^[19]。滨柃生境特殊, 在构建海岛特色景观, 改善海岛生态环境上具重要作用。因此, 关于滨柃的组培应用也有研究^[20], 但有关滨柃叶功能性状对环境因子的响应、权衡关系及适应机制的研究鲜见报道。因此, 本研究以福建省福州市平潭大练岛、海坛岛和大怀屿三地分布的滨柃为研究对象, 对滨柃叶功能性状指标与环境因子的相关性进行分析, 探究滨柃对资源的分配利用能力、揭示其对异质环境的适应机制, 为滨柃资源的保护及开发利用提供理论参考。

1 研究区概况

平潭岛地处福建省东部, 隶属于平潭综合试验区, 地理位置位于25°15′-25°45′N, 119°32′-120°10′E, 由以海坛岛为主的126个岛屿和近千个岩礁组成, 海岸线长408 km, 陆地面积为324.12 km², 为我国第五大海岛。研究区属于典型的南亚热带海洋性季风气候。年均气温为19.6 ℃, 年均降水量900-2100 mm, 年均日照1700-1980 h, 年均风速9.0 m/s, 一年中125 d以上是大风天气(7级以上), 夏秋季节常受台风影响。土壤类型多为风沙土、盐土、砖红壤, 土层瘠薄养分含量较低^[21]。植被调查中发现, 滨柃灌丛在平潭大练岛、海坛岛、大怀屿等海岛相对恶劣环境中均有一定面积和数量的分布, 因此, 本文以大练岛、海坛岛、大怀屿3个海岛的滨柃进行试验, 其中大练岛和海坛岛为有居民海岛, 大怀屿为季节性有居民海岛, 各海岛概况见表 1。

2 材料与方法 2.1 样地设置

2020年7月-9月采用标准样方法对大练岛(围营山)、海坛岛(君山)、大怀屿植物群落中的滨柃叶片进行采集, 共设置30个样方。记录每个样方的经纬度、海拔(Elevation, E)、坡度、坡向及坡位。采用便携式测光仪测定郁闭度(Crown density, CD)。各岛的环境因子概况见表 2。

2.2 滨柃样品采集与叶功能性状测定

于每个样方中四角及中心位置选取5株生长发育成熟的滨柃, 每株从东南西北四个方向采样。每一株滨柃取10片成熟、无病虫害、完全展开、边缘完整(均不带叶柄)装入标记好的自封袋, 带回实验室处理。本次共收集到30份滨柃叶片样本, 每份样本含50片, 其中大练岛有样本12份, 海坛岛9份, 大怀屿9份。采用精度为0.001 g电子天平称量叶鲜重, 将叶片于5 ℃的黑暗环境中浸泡12 h后, 用滤纸吸干叶表面水分, 称量叶饱和鲜重^[5]。采用精度为0.01 mm电子游标卡尺测量叶厚度(Leaf thickness, LT), 取叶片的叶脉旁上、中、下三个点的平均值。通过CanoScan LiDE 300扫描仪扫描叶片, 运用Image J软件计算叶长(Leaf length, LL)、叶宽(Leaf width, LW)、叶形指数(Leaf index, LI)、叶周长(Leaf perimeter, LP)、叶面积(Leaf area, LA)5个指标。将叶片于85 ℃烘箱烘干至恒重, 48 h后称量叶干重。叶碳含量和叶氮含量采用碳氮元素分析仪测定。叶磷含量采用钼锑抗比色法测定。基于相关文献资料^{[1, 5-6]}, 整理各性状指标的计算公式及生态学意义见表 3。

2.3 土壤采样与化学指标测定

采用五点交叉取样法进行土壤采样, 钻取5个25 cm深的土壤样品, 将其混合装入编号的自封袋后带回实验室进行处理, 以供土壤化学指标测定。根据鲁如坤^[22]测定土壤化学指标的方法, 对土壤的pH值(pH value, pH)、有机质(Soil organic matter, SOM)、全氮含量(Total nitrogen, T-N)、全磷含量(Total phosphorus, T-P)、全钾含量(Total potassium, T-K)、碱解氮(Alkaline hydrolyzable nitrogen, A-N)、有效磷(Available phosphorus, A-P)、速效钾(Available potassium, A-K)、水溶性盐总量(Total content of water-soluble salt, S)9项指标均进行3次重复测定。

2.4 数据处理与分析

采用Excel 2019、SPSS 25.0和Canoco 5.0软件进行数据处理和统计。其中采用变异系数(Coefficient of variation, CV)和单因素方差法(One-way analysis of variance, ANOVA)对3个海岛间滨柃各叶性状和土壤因子进行差异显著性分析, 并采用最小显著差异法(Least significant difference, LSD)比较不同数据组间的差异^[23]。利用冗余分析法(Redundancy analysis, RDA)分析环境因子对滨柃叶功能性状的影响。

3 结果与分析 3.1 不同海岛土壤因子的差异性

根据变异系数法分析结果发现(表 4), 土壤因子各指标CV介于21.61%-54.16%。单因素方差分析表明(表 4), 不同海岛的土壤因子指标间, pH、S、SOM、T-N、T-K、A-N、A-P、A-K存在极显著的差异性(P < 0.01), T-P存在显著的差异性(P < 0.05), 其中S、SOM、T-N、T-P、A-N均以大练岛最高, 以海坛岛最低。T-K、A-P、A-K以大怀屿最高, 海坛岛最低。pH以大怀屿最高, 以大练岛最低。

3.2 不同海岛滨柃叶功能性状的差异性

根据变异系数法分析结果发现(表 5), 滨柃的各叶性状CV介于2.26%-25.78%, 以叶体积(Leaf volume, LV)变异程度最大, 变异程度最小的是叶片含水量(Leaf water content, LWC)。单因素方差分析表明(表 5), 不同海岛滨柃的叶功能性状指标间差异较大, 比叶重(Specific leaf weight, SLW)、相对水分亏缺(Relative water deficit, RWD)、LL、LW、LI、LP、LA、SLA、LV和LPC存在极显著的差异性(P < 0.01), 叶组织密度(Leaf tissue density, LD)、LNC存在显著的差异性(P < 0.05), 说明这些滨柃叶功能性状在群体间多样性丰富, 存在着广泛变异, 具有较好的育种潜力, 其中LL、LW、LT、LP、LA、LV均以大怀屿最高, 以海坛岛最低。SLW以大怀屿为最高, 大练岛最低。RWD以海坛岛最高, 以大怀屿最低。SLA、LPC、LNC以大练岛最高, 大怀屿最低。LD以海坛岛最高, 大练岛最低。不同海岛滨柃的LI、LDMC、LWC和LCC差异不显著。

3.3 不同海岛样方的土壤因子聚类分析

采用系统聚类法对3个海岛的30份土壤样品(表 2中的9个土壤因子)进行聚类分析。根据聚类结果绘制聚类树(图 1), 不同海岛样方的9个土壤因子聚为一类。

3.4 滨柃叶功能性状聚类分析

采用系统聚类法对3个海岛的30份滨柃叶功能性状(表 6)进行聚类分析。根据聚类结果绘制聚类树(图 2), 在欧氏距离为22时, 3个海岛的30份滨柃样品被分为两大类, 即大练岛、海坛岛的滨柃样品聚为一类(即样方W1-W12, J1-J9), 大怀屿的滨柃样品聚为一类(即样方D1-D9), 其中大练岛、海坛岛为有居民海岛, 在地理位置相互毗邻, 遗传距离相近, 其叶片的功能性状最为接近, 聚为一类。

3.5 滨柃叶功能性状与环境因子间的冗余分析

通过对滨柃叶片功能性状与环境因子进行RDA排序, 筛选出CD、T-K、A-P、E、A-K、A-N共6个显著相关的环境因子(P < 0.05), 共同解释量为80.70%, 其中CD、T-K、A-P、E与滨柃叶功能性状呈极显著相关(P < 0.01), 说明CD、T-K、A-P、E是影响滨柃叶功能性状的主要环境因子。由图 3发现: CD、E、A-N分别与LL、LW、LP、LA、SLA、RWD、LWC、LNC、LPC呈正相关, 与其他叶性状参数呈负相关；T-K与LV、LI、LT、SLW、LD呈正相关, 与其他叶性状参数呈负相关；A-P与LL、LW、LI、LP、LA、SLA、LV、LWC、LNC、LPC呈正相关, 与其他叶性状参数呈负相关；A-K与RWD、LDMC、LCC、LD、SLW、LT呈正相关, 与其他叶性状参数呈负相关。

4 讨论 4.1 滨柃采集样方土壤理化因子的变动

海岛生境空间异质化明显, 已有的研究表明, 不同的海拔高度其水热条件不同, 土壤中的水分可能随着海拔升高逐渐降低^[14], 海岛典型常绿阔叶林和灌草丛因枯落物产量及分解速度不同, 土壤理化性质也不同^[13], 不同坡度的植物群落由于受雨水冲刷程度不同, 土壤理化性质可能差异较大^[15]。本研究, 同一海岛不同样方土壤理化因子的变异系数中, 大练岛围营山在6.85%-34.95%之间, 海坛岛君山在4.67%-53.20%之间, 大怀屿在26.07%-58.37%之间, 3个不同海岛全部样方土壤理化因子的变异系数在21.61%-54.16%之间, 同一海岛不同样方间, 3个海岛全部样方间都是pH变异系数最小, S变异系数最大。其他土壤有机质, 土壤养分含量(A-N、A-K、A-P等), 同一海岛不同样方间和3个不同海岛全部样方间也在不同程度变动着, 聚类分析结果表明, 不同海岛样方的9个土壤理化因子可能聚为一类(图 1), 因此, 样方土壤因子的数理统计分析结果表明, 同一海岛的不同样方可能因位置不同而区划为异质土壤生境。冗余分析结果表明, 3个海岛30个全部样方的土壤养分指标(A-N、A-K、A-P)与海拔高度和样方郁闭度具有相关性(图 3)。本研究中, 郁闭度大的群落为常绿阔叶林, 乔木层以台湾相思为建群种, 灌木层以滨柃为优势种, 是海岛植被演替后期的群落类型, 而郁闭度小的群落主要为常绿灌丛群落, 以滨柃为建群种的群落, 是海岛植被演替中后期的群落类型, 有研究表明, 随着演替的进行, 土壤养分和水分逐渐增加, 土壤pH逐渐下降^{[13, 24]}。可以推测, 本研究中相对于常绿灌丛群落, 常绿阔叶林群落中的台湾相思的根瘤菌具有较好的固氮作用, 能增加土壤的氮元素, 台湾相思冠幅大, 降低强风影响, 枯落物产量高且分解快, 有效维持土壤肥力及改善土壤理化性质。因此, A-N和A-P与CD呈正相关关系。而A-N和A-P与海拔呈正相关可能是因为海拔高处的滨柃采样样方刚好是台湾相思林。A-K与CD和海拔呈负相关可能是因为, 海水中富含钾离子, 能以盐雾的形式经风吹散落在土壤表层, 而海拔高可能离海岸线距离更远, 盐雾不易吹落入土壤中, 郁闭度高的常绿阔叶林乔木层的树冠能有效阻断盐雾散落在土壤中。

4.2 滨柃叶功能性状的变动

海岛植物叶片作为暴露于大气环境中的植物器官对生境因子的变动比较敏感, 主要表现为叶功能性状的变动, 以适应空间异质性^[25]。本研究中, 同一海岛不同样方的滨柃叶功能性状16项指标的变异系数中, 大练岛围营山在3.25%-27.04%之间, 海坛岛君山在1.02%-16.98%之间, 大怀屿在1.61%-10.24%之间, 3个不同海岛全部30个样方滨柃叶功能性状16项指标的变异系数在2.66%-22.99%之间, 同一海岛不同样方间, 3个海岛全部样方间都是LWC变异系数最小, 大练岛围营山12个样方采集的滨柃, 叶功能性状指标变异系数最大的是SLA, 海坛岛君山9个样方采集的滨柃, 叶功能性状指标变异系数最大的是LV, 大怀屿9个样方采集的滨柃, 叶功能性状指标变异系数最大的是LPC, 3个海岛全部30个样方采集的滨柃, 叶功能性状指标变异系数最大的是SLA。相比于30个滨柃采样样方土壤理化因子的变动系数, 30个样方滨柃样品叶功能性状的变异系数更小, 聚类分析结果表明, 大练岛和海坛岛的滨柃样品叶功能性状聚为一类, 大怀屿滨柃样品叶功能性状聚为独立的一类(图 2), 这与不同海岛样方的9个土壤理化因子可能聚为一类不一致。出现这现象的原因在于: (1)相关研究表明功能性状的种内变异系数不超过30%^{[5, 26]}, 说明植物在环境中有一定的生态幅, 其性状变异幅度是有限的^[27]；(2)不同生境因子对叶功能性状的影响不同, 可能一个因子是正向的, 另一个因子却是负向的, 如, 本研究中, 冗余分析表明, A-N和A-P与LL、LW、LP、LA、SLA、LWC、LNC、LPC呈正相关, 而A-K与这些叶功能性状指标呈负相关。

4.3 三个海岛30个样方叶功能性状与土壤、地形因子和群落郁闭度的相关性

植物根据环境差异化进行不同功能性状间的调控以形成最佳的性状组合, 来实现最大程度利用资源, 避免逆境影响, 以提高生境适合度和扩宽生态幅^[28]。本研究表明, CD、E和A-N表现出强烈的协同性, CD、E和A-N分别与滨柃的LL、LW、LP、LA、SLA、RWD、LWC、LNC、LPC呈正相关, 与滨柃其他叶性状参数(LI、LT、LDMC、SLW、LV、LD和LCC)呈负相关(图 3)。土壤养分(如本研究的A-N)越高, 滨柃叶功能性状倾向于有更大的LA、更高的SLA, 更低的LT、LDMC和LD, 因此, 在养分更充足的情况下, 滨柃倾向于更快速生长的策略。这与许多研究结果类似^{[5, 15, 29]}。值得一提的是, 土壤养分(A-N)的提高, 可能与植物群落结构有关, 本研究中, 部分样方为台湾相思群落, 其中滨柃为灌木层优势种, 上有台湾相思遮荫, 光照较弱, CD高, 而台湾相思因具有固氮菌, A-N可能也相对高, 另外一部分样方为滨柃灌丛, 这些灌丛的滨柃暴露于强光照下, CD低, 群落结构中因无台湾相思, 因而土壤A-N相对较低。同时, 相对于滨柃灌丛群落, 滨柃为灌木层优势种的台湾相思群落中, 台湾相思冠幅大, 枯落物产量高且分解快, 有效提高了养分含量。海拔高度的变化能改变生境的水热条件, 进而影响植物的功能性状, 一般大陆性植物功能性状的研究认为, 海拔升高, 植物倾向于提高养分利用率和水分保存率的保守功能性状^[29-30], 但有很多另外^[31-32], 因此, 海拔的变动影响植物功能性状首先需要具体分析, 海拔变动与水热条件变化的关系。关于海岛植物功能性状的研究, 张增可等^[14]和张秀芳等^[15]的研究发现, 叶片性状表现出随海拔上升的资源保守策略。本研究中, 海拔的升高, 对植物功能性状的影响与土壤A-N和CD对植物功能性状的影响高度一致, 随海拔升高, 植物功能性状倾向有利于捕获资源的开放策略, 而不是提高养分利用率和水分保存率的保守策略, 原因可能是, (1)低海拔滨柃的土壤生境更易受到海水盐度的影响, 如海水盐分以盐雾形式通过风作用沉降于土壤^[33], 本研究中, 海拔和土壤总水溶性盐呈负相关(偏相关分析, P < 0.05)；(2)有研究表明降雨量随海拔升高而增大^[34]。因此, 与高海拔滨柃生境相比, 低海拔滨柃生长于风大、干旱、盐分高和贫瘠的恶劣生境, 受人为干扰程度高, 滨柃叶功能性状倾向于具有更小的LA, 更低的SLA和LWC, 更高的LDMC和LD, 这些功能性状的组合有利于防止水分流失和资源消耗。综上所述, 30个研究样方中, 光照不同, 植物群落结构特征可能影响土壤养分, 海拔变动可能导致土壤盐分和水热条件变化, 是CD、E和A-N影响滨柃功能性状指标具有协同性的主要原因。

N、P和K是植物必须的养分元素, 本研究30个样方中, 土壤A-N与LNC, 土壤A-P与LPC呈正相关, 而土壤A-N和A-P含量对植物功能性状的影响也呈现较为一致的协同性, A-N和A-P分别与滨柃的LL、LW、LP、LA、SLA、LWC、LNC、LPC呈正相关, 而与滨柃其他叶性状参数(LT、LDMC、SLW、LD和LCC)呈负相关(图 3), 仅土壤A-N与RWD呈正相关, 而土壤A-P与RWD呈负相关, 土壤A-N与LI和LV呈负相关, 而土壤A-P与LI和LV呈正相关。本研究土壤A-N与RWD呈正相关, 而土壤A-P与RWD呈负相关的原因可能是, 水分与氮磷在叶各项生理活动中的互作方式不同, 有研究表明, 遭受水分胁迫时, 叶氮含量的增加有利于气孔关闭, 增进叶细胞保水^[35], 而叶磷含量增加能增大气孔导度, 主要促进光合作用和呼吸作用^[36], 具体的机理机制需进一步研究。K, 与N和P一样, 是植物必需的养分元素, 本研究中, 随土壤A-K的增加, 叶保水功能性状增强, 如RWD、LT、LDMC、LD和LCC增大, 而与生长相关的功能性状降低, 如LA和SLA等。土壤中的K作为可为植物吸收的必需元素, 没有与土壤N和P一样, 对叶功能性状的影响表现出一致的协同性, 可能原因在于, 土壤A-K与土壤总水溶性盐具正相关(偏相关分析P < 0.05), 土壤总水溶性盐升高, 植物易遭受水分胁迫, 叶倾向于具有保水分和养分的功能性状^[37]。

5 结论

在本研究中, 滨柃叶性状具有较大的环境可塑性, 其中以LV变异程度最大, 变异程度最小的是LWC。不同海岛间滨柃大部分叶性状指标差异显著, 各岛的环境因子均存在显著差异。滨柃叶片功能性状通过一定的性状变异和性状组合适应海岛生境的变化, 其中显著影响滨柃各叶性状变异的主要环境因子为CD、T-K、A-P、E、A-K、A-N, 不同海岛的滨柃叶性状对不同异化生境响应和适应程度各不相同, 表明各海岛的滨柃对生境的适应策略存在较大差异, 使得滨柃成为海岛植被中广泛分布的优势群落格局。通过研究滨柃叶功能性状特征及其与环境因子的关系, 有利于更好理解海岛植物适应环境所采取的生态策略, 可为滨柃资源的引种驯化、选育工作及开展其遗传改良等提供理论指导。但由于本文未对人为干扰因子、风因子、土壤含水量、容重等数据收集可能会造成研究结果出现偏差。