摘要
滨海湿地是全球气候变化与人类活动的敏感区和脆弱区。为探究滨海湿地浮游植物群落时空差异及影响因子,本研究在上海鹦鹉洲湿地与外围城市河道共设置14个采样点分4个季节进行浮游植物样品采集,共鉴定浮游植物7门97种,硅藻门和绿藻门种类最为丰富。鹦鹉洲湿地与外围河道共有种有61种(占总种类的62.9%),湿地特有种19种、外围河道特有种17种。种类数的季节差异显著,表现为夏季(89种)>冬季(68种)>秋季(49种)>春季(42种),四季共有种23种,夏季特有种12种,其他季节特有种仅2~3种。从生物量来看,鹦鹉洲湿地浮游植物群落类型比外围河道表现出更多样的季节演替:湿地从春季到冬季表现为硅藻-蓝|裸|甲藻-隐藻-裸|绿藻的演替规律,而外围河道则表现为隐|硅藻-蓝藻-隐|硅藻-隐|硅藻的季节演替规律。鹦鹉洲湿地浮游植物主要优势种为细小隐球藻(Aphanocapsa elachista)、尖尾蓝隐藻(Chroomonas acuta);外围河道主要优势种为啮蚀隐藻(Cryptomonas erosa)、尖尾蓝隐藻、细小隐球藻、歪头颤藻(Oscillatoria curviceps)以及颗粒直链藻(Melosira granulata),优势种季节更替明显。水温、溶解氧、透明度、盐度及营养盐是影响上海鹦鹉洲湿地及外围河道浮游植物群落分布最主要的环境因子,其中盐度是区分湿地与河道浮游植物群落的关键因子,水温与氨氮是区分四季浮游植物群落的关键因子。变差分解显示,环境因子对功能离散度(FDiv)的解释率(18.8%)显著高于对Shannon-Wiener多样性的解释率(10%),表明环境对物种生态位的筛选强于对物种个体的筛选。受环境影响不显著的功能均匀度(FEve)却受时间因子的显著影响,可能与季节更替过程中气象条件、水体生态系统营养级结构的变化、浮游生物群落季节性演替过程中种间关系等因素有关。功能丰富度(FRic)表明秋季是四季中浮游植物群落抗干扰能力最弱的季节,且湿地的抗干扰能力显著强于河道。环境因子对FRic的解释量(41.2%)显著高于对物种丰富度的解释量(16%),表明FRic的环境敏感性比物种丰富度高。本研究将对滨海区域生物多样性保护、生态系统功能恢复与管理提供科学依据。
Abstract
Coastal wetlands are sensitive and vulnerable areas to global climate change and anthropogenic pressures. To study the spatio-temporal differences and driving factors of phytoplankton communities in coastal wetlands, we collected phytoplankton samples from 14 sites for four seasons in the Yingwuzhou Wetland (YWW) and its surrounding urban canals (SUC). A total of 97 phytoplankton species belonging to 7 phyla were identified during the study period, most of which are Bacillariophyta and Chlorophyta. A total of 61 phytoplankton species (62.9% of total species) were identified in both sites. Besides, 19 species were only identified in the YWW, while 17 species were specific to SUC. The richness of phytoplankton species exhibited significantly seasonal changes, with species numbers ranked in summer (89 species)>winter (68 species)>fall(49 species)>spring (42 species). There are 23 shared species in all seasons, 12 species only identified in summer, and 2-3 typical species identified in other seasons. More diverse seasonal succession was found in YWW than those in SUC in terms of biomass in phytoplankton community types. The communities in YWW showed a succession from Bacillariophyta type(spring) to Cyanophyta-Euglenophyta-Chlorophyta type (summer), to Cryptophyta type (fall), and to Euglenophyta-Chlorophyta type (winter). While in SUC communities showed the other pattern from Cryptophyta-Bacillariophyta type (spring), to Cyanophyta-Euglenophyta-Pyrrophyta type (summer), to Cryptophyta-Bacillariophyta type (fall), and to Cryptophyta-Bacillariophyta type (winter). Seasonal dominant species varied significantly in both YWW and SUC. The dominant species of phytoplankton community in YWW were Aphanocapsa elachista and Chroomonas acuta, while the main dominant species of phytoplankton community in SUC were Cryptomonas erosa, Chroomonas acuta, Aphanocapsa elachista, Oscillatoria curviceps, and Melosira granulata. Statistical analyses revealed water temperature(WT), dissolved oxygen, secchi depth, salinity, and nutrients were the main physicochemical factors regulating the distribution of phytoplankton communities in this area. Salinity was the key factor in distinguish the phytoplankton community in the YWW from SUC. WT and NH3-N were the key factors in differentiating phytoplankton communities in different seasons. Variance partitioning analysis showed that the explained variance by environmental factors to the functional divergence (18.8%) was significantly higher than to the Shannon-Wiener diversity index (10%), which indicated that environmental filtering on niches was stronger than filtering on individuals. The functional evenness was not significantly influenced by the environmental factors but was significantly influenced by the temporal factors, which were possibly related to climate change, trophic structure shift of aquatic ecosystems, or biological interaction alter of seasonal plankton communities. The functional richness of both sites indicated that fall was most vulnerable to the disturbance, and the resilience of the wetland was significantly higher than that of its surrounding canals. The explained variance by environmental factors to the functional richness (41.2%) was significantly higher than to the species richness (16%), which indicated that functional richness was more sensitive to environmental changes than species richness. This study may provide scientific references for biodiversity conservation, ecosystem function restoration, and management in coastal areas.
滨海湿地作为陆地、河流和海洋生态系统的生态交错区,具有固碳能力强、生物多样性高、对环境变化极为敏感的特点,是滨海河口区域生态系统结构与功能维持不可或缺的组成部分。然而,近半个世纪以来,气候变化以及沿海地区人类活动导致的生境破坏、环境污染等因素使得滨海湿地面积不断缩小、碳汇锐减、生物多样性降低,陆地与海洋生态系统服务功能均受到严重威胁[1]。作为水域生态系统的一部分,湿地与周边城市化水系的物质、能量、信息交流频繁,联系紧密。入海河流是连接陆地与海洋生态系统的重要廊道,不断为滨海湿地及近海海域输送邻近陆域人类活动产生的营养物质、微塑料等污染物,导致近岸海域有害藻类水华暴发,滨海湿地生物多样性及生态系统功能受到威胁。
浮游植物作为水生态系统中的初级生产者,在滨海河口区域生态系统中居重要地位。浮游植物个体小、世代周期短、对环境变化敏感,能够对湿地水环境状况做出迅速的响应。营养盐、透明度、浊度、盐度等环境因子[2-4]的变化会导致湿地浮游植物群落结构发生改变,进而影响湿地生态系统的结构与功能。
鹦鹉洲生态湿地是正在建设的上海首个国家海洋公园的重要组成部分,是上海市“十三五”期间开展海洋生态保护修复取得的重要成果。该湿地是杭州湾北岸以提升生态系统服务为目标的人工恢复盐沼湿地,通过综合的海岸带生态工程恢复原有的受损滨海生态系统,已成为兼具生态修复和科普游憩等功能的城市滨海湿地,在加强区域海洋环境保护、改善海洋生态环境、维持生物多样性等方面发挥着积极作用。该人工恢复滨海湿地的环境独特且复杂,弄清此区域生物群落动态的驱动因子可以为盐沼湿地的生态修复与恢复、生物多样性维持等提供重要依据。目前关于鹦鹉洲湿地的研究主要集中于湿地生态修复[5-6]、景观格局变化[7]、鸟类多样性研究[8]等方面,有关浮游植物群落的相关研究相对较少,缺乏对湿地与周边城市水系及入海河流浮游植物群落动态变化及驱动因子的研究。本研究将上海鹦鹉洲湿地及外围河道中的浮游植物群落作为研究对象,从时间、空间和环境3个维度探究滨海湿地与周边区域城市河道水体中浮游植物群落差异的主要影响因子,以及物种多样性、功能多样性对环境因子和时间因子的响应差异,为滨海河口区域生态系统功能维持以及水生生物多样性保护提供科学依据。
1 材料与方法
1.1 研究区域概况
上海鹦鹉洲湿地及外围河道(30°41′~30°46′N,121°17′~121°23′E)位于上海市金山区西南部杭州湾北岸,人口密集,工、农业活动频繁。研究区域属于黄浦江上游支流水系,年平均降雨量为1159.4 mm,为亚热带季风气候。湿地面积为23.2 hm2,平均水深为1.8 m,由湿地净化展示区、盐沼湿地恢复区、自然湿地引鸟区三大部分组成。上海鹦鹉洲湿地在固碳增汇方面优于周边自然湿地,同时也是金山区发现鸟类种数最多的区域,生态系统服务价值极高[9]。鹦鹉洲湿地自东侧金山城市沙滩海域进水,湿地进水口涵洞设置水阀,湿地内水位不受邻近海域潮汐影响,海水经湿地净化后重新流入金山城市沙滩海域。外围河道取邻近鹦鹉洲湿地北部的部分金山区城市河道,金山区是典型平原河网地区,区辖内河道总长度为2500.26 km,面积为42.68 km2,占区域面积的7.28%,河水出海汇入杭州湾金山城市沙滩海域后与鹦鹉洲湿地连通。
1.2 采样点布设
本研究在上海金山鹦鹉洲生态湿地以及外围河道共设置14个采样点,于2022年8月(夏季)、2022年11月(秋季)、2023年2月(冬季)、2023年5月(春季)进行浮游植物定性、定量样品和水样采集。其中上海金山鹦鹉洲湿地设置7个采样点(W1~W7),依次位于湿地进水口到出水口之间的不同湿地处理单元;外围河道设置7个采样点(R1~R7),分别位于鹦鹉州湿地上游的金卫城河(R1、R7)、张泾河(R2)、红旗港(R3)、龙泉港(R4)、龙泉港出海闸(R5)、老龙泉港(R6)城市河道及其交汇处(图1)。

图1上海鹦鹉洲湿地与外围河道采样点分布
Fig.1Distribution of sampling sites in the Yingwuzhou Wetland and its surrounding canals in Shanghai
1.3 样品采集与分析
1.3.1 水体理化样品采集与测定
利用哈希便携式多参数水质分析仪(HQ30d)现场测定水温(WT)、酸碱度(pH)、溶解氧(DO);采用塞氏盘测定水体透明度(SD);采用盐度计(MASTER-S/Millα)测定盐度(Sal);采用浊度计(WGZ-500B)测定浊度(Turb)。使用有机玻璃采水器在水面以下0.5 m处采集水样500 mL于采样瓶中,水样冷藏保存立即送往实验室测定叶绿素a(Chl.a)、硝态氮(NO-3-N)、亚硝态氮(NO-2-N)、氨氮(NH3-N)、溶解态无机氮(DIN)、溶解态无机磷(DIP)、化学需氧量(COD)等指标,测定方法参考《海洋调查规范第4部分:海水化学要素调查》(GB/T12763.4—2007)。 NO-3-N浓度采用锌镉还原法、NO-2-N浓度采用N-(1-萘基)-乙二胺分光光度法、NH3-N浓度采用次溴酸盐氧化法、溶解态无机磷采用钼酸铵分光光度法、化学需氧量采用碘化钾碱性高锰酸钾分光光度法进行测定。DIN浓度为NO-3-N、NO-2-N、NH3-N浓度之和。
1.3.2 浮游植物样品采集
浮游植物定性样品采集:采用25#浮游生物网(孔径64 μm)在水体表层至水面以下0.5 m处拖动数次,待水滤去后打开浮游生物网底部阀门,将过滤后得到的浮游植物样品置于50 mL标本瓶中,在瓶身贴上标签注明采样时间、地点等基本信息,加入4%甲醛溶液保存待检。浮游植物定量样品采集:采用有机玻璃采水器在水面以下0.5 m处采集水样置于1 L采样瓶中,加入1.5%鲁哥试剂,静置48 h后利用虹吸法缓缓吸去上清液,将浮游植物样品浓缩至50 mL,并加入4%甲醛溶液保存后待镜检计数。
1.3.3 浮游植物镜检鉴定
参照《中国淡水藻志》《中国淡水藻类——系统、分类及生态》等文献[10-12],在显微镜(Motic BA200)下对上述浮游植物定性样品进行种类鉴定。将浓缩后的定量样品充分混合后吸取0.1 mL于计数框中,在10×40倍显微镜下对浮游植物细胞进行计数,每个样品重复2次。取平均值计算每个样品中的浮游植物细胞数量生物量,并换算为重量生物量。
1.4 数据处理与统计分析
使用ArcGIS绘制采样点图;优势度指数Y≥0.02确定为优势种;将优势种相对丰度大于10%的功能类群定义为优势功能类群,功能类群的划分依据Reynolds等[13-15]提出的FG功能群分类方法确定。采用物种丰富度(SRic)、Shannon-Wiener多样性指数(H′)以及基于Simpson多样性指数计算的Simpson均匀度指数(Jsi)作为浮游植物群落的物种多样性特征。
基于浮游植物丰度以及细胞大小、几何形态、是否群体、行为特征、生殖方式、营养模式等形态、生理、行为三类功能性状特征分别计算功能离散度(FDiv)、功能均匀度(FEve)以及功能丰富度(FRic)3个功能多样性指数。其中功能离散度数值越高说明浮游植物生态位差异越大,种间资源竞争越弱;功能均匀度数值越高说明功能空间资源利用越充分,反之数值越低表明群落中物种功能性状越相似,资源利用率越低;功能丰富度值越高则反映群落中浮游植物功能性状占据的功能空间越大,群落对干扰的承受能力越强。
(1)
(2)
(3)
(4)
(5)
式中,S为物种数,Pi为物种i个体数在总个体数中的比例,d为多度权重离散度,为所有出现物种与重心距离的平均值,l为分支长,PEWl为局部加权均匀度,xi为物种i的功能性状值。
为探究鹦鹉洲湿地及外围河道浮游植物群落差异的主要影响因子,开展了非度量多维尺度分析(non-metric multidimensional scaling,NMDS)、分类回归树分析(classification and regression tree analysis,CART)、冗余分析(redundancy analysis,RDA),数据统计分析、可视化、群落指标计算及指示种分析均在R(version 4.3.1)语言中完成。选择在该季节及区域优势度指数≥0.02的物种进行排序分析,排序分析前先进行去趋势对应分析(detrended correspondence analysis,DCA),根据梯度长度值(前两轴均<4)采用RDA排序,并将所有测得的环境因子经前项选择并选出显著的环境因子分析其与浮游植物优势种的关系。此外,将挑选出的所有显著环境因子作为一组环境变量,取样时间(季节)作为一组时间变量,分别对浮游植物物种多样性和功能多样性进行变差分解。定量分析物种与功能多样性对环境变量和时间变量的响应。
2 结果与分析
2.1 鹦鹉洲湿地与外围河道水体理化指标特征
上海鹦鹉洲湿地与外围河道水体理化因子存在显著的季节与空间差异(P<0.01,表1)。鹦鹉洲湿地及外围河道水体pH和水温季节变化一致,水体均呈弱碱性,四季水温变化范围为9.10~32.14℃,河道水温略高于鹦鹉洲湿地。两个水体的浊度、盐度和透明度季节差异均较大,外围河道的浊度均高于鹦鹉洲湿地,鹦鹉洲湿地的盐度和透明度则均高于外围河道。春季鹦鹉洲湿地的叶绿素a浓度最高,夏季最低;河道的叶绿素a浓度夏季最高,秋季最低。湿地和河道的溶解氧浓度均表现为秋、冬季较高,春、夏季较低。夏季外围河道的叶绿素a以及溶解氧浓度均高于湿地,其他季节则低于湿地。四季中,鹦鹉洲湿地的NO-3-N、NO-2-N、NH3-N、DIN以及DIP浓度均低于河道,COD则高于河道。冬季鹦鹉洲湿地及外围河道DIN浓度较高,春季鹦鹉洲湿地及外围河道DIP浓度较高。
表1上海鹦鹉洲湿地及外围河道水体理化指标
Tab.1 Physical and chemical indicators of water bodies of Yingwuzhou Wetland and its surrounding canals in Shanghai

“/”表示缺失值。
2.2 鹦鹉洲湿地和外围河道浮游植物种类组成
上海鹦鹉洲湿地和外围河道的四季调查共鉴定出浮游植物97种,隶属于7门28科62属。其中硅藻门(Bacillariophyta)的种类最为丰富,为25属40种(占总物种数的41.24%);绿藻门(Chlorophyta)次之,为21属33种(占比34.02%);裸藻门(Euglenophyta)5属11种(占比11.34%);蓝藻门(Cyanophyta)6属7种(占比7.22%);其他门种类较少(图2)。
浮游植物门水平上的组成在鹦鹉洲湿地和外围河道无明显差异,均以硅藻门和绿藻门为主。鹦鹉洲湿地共检出浮游植物80种,主要包括硅藻门34种、绿藻门27种、裸藻门8种。外围河道共检出浮游植物78种,主要包括硅藻门37种、绿藻门23种、裸藻门8种。但浮游植物物种组成则在湿地和河道存在一定差异性。两类水体共有种有61种(图3a),鹦鹉洲湿地特有种有19种,主要包括绿藻门10种,有胶球藻(Coccomyxa sp.)、柠檬形顶棘藻(Lagerheimia citriformis)、细小单针藻(Monoraphidium minutum)等;硅藻门3种,有优美旭氏藻(Schroederella delicatula)、近缘桥弯藻(Cymbella affinis)、尖头桥弯藻(Cymbella cuspidata);裸藻门3种,有梨形扁裸藻(Phacus pyrum)、三脊扁裸藻(Phacus trimarginatus)、细尾鳞孔藻(Lepocinclis gracilicauda);甲藻门1种,为裸甲藻(Gymnodinium aeruginosum);蓝藻门2种,有居式腔球藻(Coelosphaerium kuetzingianum)、微小色球藻(Chroococcus minutus)。外围河道特有种有17种,主要为绿藻门6种,有空球藻(Eudorina elegans)、汉斯集星藻(Actinastrum hantzschii)等;硅藻门6种,有颗粒直链藻(Melosira granulata)、肘状肘形藻(Ulnaria ulna)等;裸藻门3种,有旋转囊裸藻(Trachelomonas volvocina)、曲尾扁裸藻(Phacus lismorensis)、陀螺藻(Strombomonas sp.);蓝藻门2种,有惠氏微囊藻(Microcystis wesenbergii)、螺旋浮丝藻(Planktothrix spiroides)。鹦鹉洲湿地中浮游植物种类数存在空间差异,湿地进水口(21种)种类数低于出水口(33种);种类组成上,湿地进水口隐藻门、蓝藻门、甲藻门浮游植物比例相对较高。

图2上海鹦鹉洲湿地及其外围河道浮游植物种类组成
Fig.2The composition of phytoplankton species in the Yingwuzhou Wetland and its surrounding canals in Shanghai
本研究区域内浮游植物种类组成季节差异显著(P < 0.01,图3b)。4个季节浮游植物物种丰富度从高到低依次为:夏季(89种)、冬季(68种)、秋季(49种)、春季(42种)。在种类最为丰富的夏季,鹦鹉洲湿地及外围河道均检出浮游植物63种,以硅藻门(湿地27种,河道28种)及绿藻门(湿地22种,河道20种)为主。研究区域内四季共有种23种,主要为绿藻门衣藻属(Chlamydomonas)等4属;硅藻门脆杆藻属(Fragilaria)等10属;裸藻门裸藻属(Euglena);隐藻门蓝隐藻属(Chroomonas)、隐藻属(Cryptomonas)的种类。春、秋季浮游植物共有种最少(29种),夏、冬季共有种最多(66种)。夏季特有种12种,均为蓝藻门和绿藻门的种类;春季特有种3种,秋、冬季特有种均为2种。

图3上海鹦鹉洲湿地及其外围河道浮游植物物种组成的空间(a)与季节(b)差异
Fig.3Spatial (a) and seasonal (b) variation of phytoplankton species composition in the Yingwuzhou Wetland and its surrounding canals in Shanghai
2.3 鹦鹉洲湿地及其外围河道浮游植物生物量
研究区域内浮游植物年平均生物量为0.61 mg/L。鹦鹉洲湿地及其外围河道浮游植物年平均生物量相当,分别为0.58和0.62 mg/L。鹦鹉洲湿地内进、出水口的浮游植物年平均生物量差异较大,分别为0.62、0.25 mg/L。鹦鹉洲湿地及外围河道浮游植物生物量季节差异显著(P < 0.01,图4),湿地较河道的季节变幅小。鹦鹉洲湿地与外围河道浮游植物生物量季节变化均表现为夏季最高(分别为0.87、1.89 mg/L),秋季(0.77、0.12 mg/L)和冬季(0.63、0.46 mg/L)次之,春季最低(0.03、0.008 mg/L)。夏季,鹦鹉洲湿地浮游植物生物量以蓝藻门(0.27 mg/L,占总生物量的31%)、裸藻门(0.19 mg/L,占21.3%)、甲藻门(0.15 mg/L,占16.9%)为主,而外围河道仅以蓝藻门为主(1.59 mg/L,占84.1%)。秋季,湿地和河道浮游植物生物量均转变为以隐藻门(分别为0.64、0.05 mg/L,占总生物量的82.7%、43.73%)为主。冬季,湿地浮游植物生物量转变为以裸藻门(0.20 mg/L,占32.36%)和绿藻门(0.12 mg/L,占19.2%)为主,而河道则以隐藻门(0.25 mg/L,占54.21%)和硅藻门(0.14 mg/L,占30.8%)为主(图4)。

图4上海鹦鹉洲湿地(a)及其外围河道(b)浮游植物生物量(mg/L)组成及季节差异
Fig.4Composition and seasonal variation of phytoplankton biomass (mg/L) in the Yingwuzhou Wetland (a) and its surrounding canals (b) in Shanghai
2.4 鹦鹉洲湿地及其外围河道浮游植物优势种
鹦鹉洲湿地共确定浮游植物优势种6种,分别为蓝藻门的细小隐球藻(Aphanocapsa elachista),硅藻门的小环藻(Cyclotella sp.)、劳氏角毛藻(Chaetoceros lorenzianus),甲藻门的裸甲藻,隐藻门的尖尾蓝隐藻(Chroomonas acuta)和啮蚀隐藻(Cryptomonas erosa)。外围河道共确定浮游植物优势种9种,分别为蓝藻门的细小隐球藻、惠氏微囊藻、歪头颤藻(Oscillatoria curviceps),硅藻门的小环藻、脆杆藻(Fragilaria sp.)、劳氏角毛藻、颗粒直链藻,隐藻门的尖尾蓝隐藻、啮蚀隐藻、卵形隐藻(Cryptomonas ovata)。其中,细小隐球藻、尖尾蓝隐藻、啮蚀隐藻是湿地与河道的共有优势种(表2)。
细小隐球藻是鹦鹉洲湿地夏、冬季以及外围河道夏季的绝对优势种(Y分别为0.40、0.48、0.17);尖尾蓝隐藻是鹦鹉洲湿地春、秋季以及外围河道春、冬季的绝对优势种(Y分别为0.46、0.64、0.32、0.29);啮蚀隐藻是外围河道春、秋、冬季的绝对优势种(Y分别为0.38、0.20、0.29)。
2.5 鹦鹉洲湿地及其外围河道浮游植物物种多样性及功能多样性
鹦鹉洲湿地及外围河道浮游植物物种多样性与功能多样性存在显著的时空差异(图5)。季节上,外围河道浮游植物物种多样性与功能多样性的季节差异均比湿地的季节差异大,即湿地的季节波动相对稳定。湿地与外围河道浮游植物Shannon-Wiener多样性指数与Simpson均匀度指数的季节变化表现为春季>夏季、冬季>秋季,外围河道的功能离散度与功能均匀度也表现出相似的季节变化规律,但湿地的功能离散度与功能均匀度却呈现相反的季节变化规律,均在秋季最高,春、夏较低。此外,功能丰富度则与物种丰富度表现出相似的季节变化。空间上,外围河道浮游植物Shannon-Wiener多样性指数和Simpson均匀度指数各季节均略高于湿地,而物种丰富度则表现为湿地略高于河道。
表2上海鹦鹉州湿地及外围河道浮游植物优势种和优势度指数*
Tab.2 Dominant species and dominance index of phytoplankton in the Yingwuzhou Wetland and its surrounding canals in Shanghai

*“-”表示在该季节没有成为优势种。

图5上海鹦鹉洲湿地及其外围河道浮游植物物种多样性和功能多样性
Fig.5Species diversity and functional diversity indices of phytoplankton in Yingwuzhou Wetland and its surrounding canals in Shanghai
2.6 鹦鹉洲湿地及其外围河道浮游植物与环境因子的关系
非度量多维尺度分析结果显示,鹦鹉洲湿地和外围河道4个季节浮游植物群落差异性较大(图6),均具有明显的群落特征(表3)。鹦鹉洲湿地春季浮游植物群落位于排序轴的左上部分,群落优势功能类群为X2、B、Y,硅藻门是浮游植物生物量的主要组成部分。外围河道春季浮游植物群落位于排序轴的左下部分,群落优势功能群为Y、X2,浮游植物生物量以隐藻门为主。夏季鹦鹉洲湿地及其外围河道浮游植物群落位于排序轴的右下部分,夏季鹦鹉洲湿地浮游植物优势功能类群为K,生物量以蓝藻门为主,指示种分析显示,此群落指示种为裸甲藻。夏季外围河道浮游植物优势功能类群为K、M、MP,生物量以蓝藻门为主,小环藻、歪头颤藻、螺旋浮丝藻是该群落的指示种。秋季鹦鹉洲湿地及外围河道浮游植物群落位于排序轴y轴正轴两侧。秋季湿地群落中浮游植物优势功能群为X2、Y,以隐藻门生物量为主,指示种为硅藻门二头脆杆藻(Fragilaria biceps)和钝脆杆藻(Fragilaria capucina)。秋季外围河道浮游植物群落优势功能类群为Y、P、X2,以隐藻门生物量为主,卵形隐藻是该群落的指示种。冬季鹦鹉洲湿地及其外围河道浮游植物群落位于排序轴的右上部分,其中湿地中浮游植物优势功能类群为K,生物量以裸藻门为主。绿藻门的球衣藻(Chlamydomonas globosa)、胶球藻,硅藻门的优美旭氏藻、劳氏角毛藻以及蓝藻门的细小隐球藻是该群落的指示种;而外围河道浮游植物优势功能类群为X2、P、Y,生物量以隐藻门为主,指示种有硅藻门的颗粒直链藻、隐藻门的啮蚀隐藻以及金藻门的黄群藻(Synura sp.)。
表3上海鹦鹉州湿地及外围河道浮游植物功能类群、群落类型及指示种
Tab.3 Functional groups, community types and indicator species of phytoplankton in the Yingwuzhou Wetland and its surrounding canals in Shanghai

“*”表示该季节无指示种。
前项选择结果表明,WT、Turb、pH、Chl.a、DIP、Sal、DO、SD是造成以上群落差异性的显著环境因子(图6)。其中影响群落差异最主要的环境因子是WT、DO、SD,其次为Sal、Turb、pH、DIP、Chl.a。WT、Turb、DIP 与春、夏季鹦鹉洲湿地和外围河道浮游植物群落呈正相关关系,pH与夏季鹦鹉洲湿地和外围河道浮游植物群落呈正相关,Chl.a与春季群落呈正相关关系;Sal、DO、SD与秋、冬季鹦鹉洲湿地和外围河道浮游植物群落呈正相关关系。
分类回归树结果表明上海鹦鹉洲湿地及外围河道浮游植物群落时空分布主要受Sal、WT、NH3-N的影响(图7)。仅考虑空间维度时,Sal是造成上海鹦鹉洲湿地及外围河道浮游植物群落差异的关键因子(图7a),Sal值<3.2‰时划分为外围河道浮游植物群落;Sal值≥3.2‰时划分为鹦鹉洲湿地浮游植物群落。仅考虑时间维度时,WT和NH3-N是造成上海鹦鹉洲湿地及外围河道浮游植物群落季节差异的关键因子(图7b),当19℃≤WT值<28℃时划分为春季浮游植物群落;WT值≥28℃时划分为夏季浮游植物群落;WT值<19℃且NH3-N值<0.003 mg/L时划分为秋季浮游植物群落;WT值<19℃且NH3-N值≥0.003 mg/L时划分为冬季浮游植物群落。
冗余分析结果表明,WT、SD、Sal、DO、DIP、NH3-N、NO-2-N是鹦鹉洲湿地及外围河道浮游植物优势种分布具有显著解释性的影响因子(P<0.01),前两轴的累计解释量为56.08%。第一排序轴与WT呈正相关,与Sal呈负相关;第二排序轴与SD、DO呈正相关,与DIP、NO-2-N呈负相关。优势种中,歪头颤藻与NH3-N、WT呈显著正相关,与Sal呈显著负相关;颗粒直链藻、裸甲藻与NH3-N呈显著正相关;脆杆藻与SD、DO、NH3-N呈正相关;啮蚀隐藻、尖尾蓝隐藻、细小隐球藻、劳氏角毛藻与SD、DO、Sal呈显著正相关;小环藻、卵形隐藻、惠氏微囊藻与WT、DIP、NO-2-N呈显著正相关(图8)。

图6上海鹦鹉洲湿地及其外围河道浮游植物群落NMDS分析
Fig.6Non-metric multidimensional scale analysis (NMDS) of phytoplankton communities in the Yingwuzhou Wetland and its surrounding canals in Shanghai

图7上海鹦鹉洲湿地及其外围河道浮游植物群落空间维度(a)和时间维度(b)分类回归树分析
Fig.7Classification and regression tree analysis of phytoplankton communities divided by space (a) and season (b) in the Yingwuzhou Wetland and its surrounding canals in Shanghai
2.7 鹦鹉洲湿地及外围河道浮游植物物种多样性、功能多样性与环境和时间因子的关系
环境和时间(季节)两组因子对鹦鹉洲湿地及外围河道浮游植物多样性不同指数的解释率差异较大。首先,两组因子对物种丰富度和功能丰富度的总解释率最高(分别解释了29.3%和44.6%),环境因子单独、极显著地解释了16%(P<0.01)和41.2%(P<0.001)的物种和功能丰富度变异,时间因子则主要与环境因子交互解释了物种丰富度14.1%和功能丰富度2.3%的变异量。其次,两组因子对Simpson均匀度指数和功能离散度的总解释率较高(分别解释了25.4%和27.9%),环境因子显著解释了22.1%和18.8%的Simpson均匀度和功能离散度变异,时间因子主要与环境因子交互解释了3%的Simpson均匀度和9.1%的功能离散度变异量。此外,两组因子对Shannon-Wiener多样性和功能均匀度的总解释率最低(分别解释了21.5%和10.1%)。其中,环境因子仅显著解释了10%的Shannon-Wiener多样性变异量,同时解释了3.8%的功能均匀度,时间因子单独、显著解释了6.3%的功能均匀度,并与环境因子交互解释了11.5%的Shannon-Wiener多样性变异量(图9)。

图8上海鹦鹉洲湿地及其外围河道浮游植物优势种与环境因子的RDA分析
Fig.8Redundancy analysis (RDA) of dominant species of phytoplankton and environmental factors in the Yingwuzhou Wetland and its surrounding canals in Shanghai

图9环境因子和时间因子对上海鹦鹉洲湿地及其外围河道浮游植物物种多样性和功能多样性的变差分解(*P<0.05,**P<0.01,***P<0.001)
Fig.9Variance partitioning analysis of temporal and environmental factors for species diversity indices and functional diversity indices of phytoplankton in the Yingwuzhou Wetland and its surrounding canals in Shanghai (*P<0.05, **P<0.01, ***P<0.001)
3 讨论
3.1 鹦鹉洲湿地与外围河道浮游植物群落特征及时空差异
本研究在位于杭州湾北岸的上海鹦鹉洲湿地与外围河道的四季调查中共鉴定浮游植物97种,物种丰富度高于Zhang等[18]2004—2010年夏季在杭州湾报道的84种、刘炜等[19]2018年秋季在杭州湾北岸人工潟湖报道的54种;与杨娜等[20]2018年夏季在杭州湾北岸人工潟湖鉴定到的96种接近;略低于刘巧等[21]2017年7—12月在杭州湾生物修复区与非修复区记录到的103种。不同时间和空间尺度上物种丰富度差异在河口区域研究中较为常见,杭州湾河口区域水动力条件决定了水环境条件极不稳定,水体环境因子显著的时空变化使得浮游植物时空差异较大。上海鹦鹉洲湿地与外围河道浮游植物种类组成与已有河口区域浮游植物研究[22-23]类似,均以硅藻门为主,绿藻门次之。绿藻门、裸藻门浮游植物在总种类数中的比例高于同区域其他研究,蓝藻门、甲藻门的组成比例则相对较低[18-20]。鹦鹉洲湿地和外围河道共有种有61种,占总种类数的62.9%,可能是因为鹦鹉洲湿地与外围河道具有间接连通性,外围河道中包括了入海河流,它是连接陆地与海洋生态系统的重要生态廊道,在不断为滨海湿地及近海海域输送物质的同时,也加强了水生生物的输入。
鹦鹉洲湿地与外围河道不同的水动力条件可能是主导两种水体浮游植物生物量组成空间差异的因素之一。外围河道中水体流速较快并且高度混合,具有较硬硅质外壳的硅藻等重量系数较大的藻类更能抵抗水流冲刷等物理环境[24-25]而大量生长,使得浮游植物生物量较高。而湿地水体流速相对较缓、水体更加稳定,偏好弱紊动水体的小型绿藻、裸藻等具有较小重量系数的类群占比较高,从而导致鹦鹉洲湿地浮游植物生物量低于外围河道。从生物量季节演替的角度看,鹦鹉洲湿地从春季到冬季表现出硅藻-蓝|裸|甲藻-隐藻-裸|绿藻的演替规律,而外围河道则表现出隐|硅藻-蓝藻-隐|硅藻-隐|硅藻的季节演替规律。鹦鹉洲湿地浮游植物群落类型的季节更替比外围河道更明显、更多样,可能与湿地本身作为季节性积水地带/水域会随季节性气候更替形成不同、但季节内相对稳定的水环境有关。与鹦鹉洲湿地不同,除夏季外外围河道均是以隐藻生物量为主的群落类型,可能是隐藻喜生于氮丰富的水体中[26]且外围河道的氮浓度远高于湿地所致。湿地春季和外围河道冬季硅藻生物量也较高,与春、冬季水温较低以及光照不充分限制了很多藻类生长、但硅藻仍然能获得生长优势有关[27-29]。此外,冬季水量少、流速缓、水力停留时间长、颗粒物沉降等因素使得水体透明度高于其他季节,促进浮游植物光合作用,造成冬季两种水体的浮游植物生物量高于春、秋季,与Wang等在杭州湾海域的研究结果一致[30],同时也与冬季水体扰动减弱、为浮游植物提供稳定的生长环境有关。
鹦鹉洲湿地与外围河道浮游植物均为多优势种型群落,但优势种组成存在差异。外围河道是狭长的城市生态廊道,与周围陆地生态系统互动频繁,受周边人为活动影响较大,生境异质性较高,可以支持更多类型浮游植物大量生长繁殖,因而外围河道较鹦鹉洲湿地的浮游植物优势种组成更为多样,各优势种的优势度指数也相对较低。尖尾蓝隐藻和啮蚀隐藻是春、秋、冬三季在鹦鹉洲湿地与外围河道均占有优势地位的浮游植物种类。隐藻作为混合营养型生物以及具有特殊光合作用的浮游植物,可以通过切换不同的营养机制在低光照强度下维持种群数量[15,26,31]。夏季,蓝藻门的细小隐球藻则成为鹦鹉洲湿地与外围河道的共有优势种,这可能与蓝藻门在高温、强光条件下可以充分占据可利用空间、高效利用水体中的可获取资源而大量生长[32-33]有关。劳氏角毛藻是秋季鹦鹉洲湿地与外围河道均出现的优势种,作为一种我国近岸海域常见的广布暖水种,劳氏角毛藻常在夏、秋季海域中形成优势种[34]。
3.2 影响鹦鹉洲湿地及其外围河道浮游植物群落的主要环境因子
本研究中NMDS结合envfit分析结果显示,WT、DO、SD、Sal、Turb、pH、DIP、Chl.a是对群落时空差异具有显著解释性的环境因子;分类回归树分析显示,当只考虑季节或空间单个因子时,Sal是区分鹦鹉洲湿地和外围河道浮游植物群落的关键因子,WT和NH3-N是区分研究区域浮游植物群落季节差异的关键因子。RDA分析结果显示,WT、SD、Sal、DO、DIP、NH3-N、NO-2-N是对优势种分布具有显著解释性的影响因子。综上得出,WT、DO、SD、营养盐及Sal是上海鹦鹉洲湿地及外围河道浮游植物群落分布最主要的影响因子。
水温是具有明显季节特征的影响因子,水温上升是浮游植物代谢和大量繁殖过程的主要驱动力[35-37],较高的水温有助于蓝藻达到最大生长速率形成优势从而影响群落结构[38-40]。夏季鹦鹉洲湿地及其外围河道浮游植物群落优势种惠氏微囊藻、歪头颤藻等蓝藻受水温的影响显著(图8),这几种优势种所属浮游植物群落分布也与水温呈显著正相关(图6),同时分类回归树分析显示,水温是区分夏季与其他季节浮游植物群落的关键因子(图7)。本研究中SD、DO与秋、冬季鹦鹉洲湿地及冬季外围河道浮游植物群落分布呈正相关关系,与夏季鹦鹉洲湿地及外围河道浮游植物群落分布呈负相关。透明度高的水体环境有利于浮游植物获得光照大量繁殖[41-43];而透明度低的水体一方面由于水下光照不足从而限制浮游植物生长,另一方面由于悬浮于水中的大颗粒物的扰动或物理作用对浮游植物的生长不利。DO与浮游植物的正相关关系也在相关研究中得以证实[44-45],本研究中夏、冬季鹦鹉洲湿地优势种细小隐球藻、尖尾蓝隐藻以及秋、冬季外围河道优势种啮蚀隐藻、尖尾蓝隐藻、劳氏角毛藻均与SD和DO呈正相关。
营养盐是浮游植物生长与繁殖的重要物质基础,营养盐浓度变化会影响浮游植物的群落分布[46-48]。磷是浮游植物光合作用、能量储存以及细胞分裂等生长过程的重要元素[49-50],NH3-N是浮游植物最容易利用的氮源[51],一定浓度范围内的NO-2-N也可为浮游植物生长所用。本研究中DIP、NO-2-N与春季鹦鹉洲湿地和外围河道以及秋季外围河道浮游植物群落呈显著正相关,这些群落的优势类群均为隐藻和硅藻,两者均对低温及高营养盐的环境条件有较好的适应性而发展为优势类群[52-54]。分类回归树还表明,NH3-N是区分秋、冬季节浮游植物群落的关键因子,这与冬季枯水期水量减少导致营养盐浓度增大有关。盐度是影响河口区域浮游生物群落结构与分布的决定性因素之一[55-56],特殊的水动力条件使得河口区域淡水与海水充分混合,成为较独特的海水-淡水生态交汇区。特殊的盐度条件筛选出具有能适应一定盐度种类的浮游植物群落,这与城市内河淡水条件下的浮游植物群落差异较大,因此本研究中Sal是区分鹦鹉洲湿地与外围河道浮游植物群落差异的关键因子。除以上环境因子外,本研究未涵盖的其他环境因子也可能影响鹦鹉洲湿地以及外围河道的浮游植物群落格局,如电导率、总溶解固体、重金属含量等[57-58]。不同的环境组成可为浮游植物群落提供不同的生态位空间,使其群落组成发生改变。未来研究可进一步探究环境因子与浮游植物种群之间的具象关系,明晰种群对环境因子的具体响应机制。
除了WT、Sal、营养盐等环境因素,滨海湿地作为季节性积水地带/水域和城市入海河流的独特水文动态特征、水体生物因子及周边人类活动等因素也可成为浮游植物群落的重要影响因子。此外,扩散限制等空间过程也可能在滨海湿地的浮游植物群落构建中发挥作用,对应浮游植物群落组成之间的相似度随着样点间地理距离发生改变。
3.3 浮游植物物种多样性与功能多样性对环境和时间因子的响应差异
外围河道浮游植物物种多样性与功能多样性的季节差异均比湿地的季节差异大,即湿地生物多样性的季节波动相对河道稳定(图5)。这可能与湿地的调节功能有关[59],其对季节性水文条件的调节及对季节变化所带来扰动的调节,对生物多样性起到一定的维稳作用。
物种多样性和功能多样性对环境与时间因子(季节)的响应存在较大差异(图9)。功能离散度阐释了浮游植物生态位差异及与种间资源竞争的强弱[60]。本研究中,环境因子显著解释了Shannon-Wiener多样性和功能离散度的变异,且对功能离散度的解释率(18.8%)显著高于对Shannon-Wiener多样性的解释率(10%),这表明环境对物种生态位的筛选要强于对物种个体的筛选。一般来说,如优势种分布在功能空间的边缘,则功能离散度高,表明生态位分化明显,种间资源竞争较弱;反之,如优势种接近功能空间中心,则功能离散度低,表明种间资源竞争激烈[17,61]。本研究结果中,外围河道春、冬季与湿地秋季的功能离散度明显高于其他季节,表明群落种间资源竞争较弱,可能原因在于这3个群落均以啮蚀隐藻和尖尾蓝隐藻占绝对优势,两个种的生态幅均较宽[62-63],两者个体大小差异较大、且为混合营养型生物以及具有特殊光合作用,可以通过切换不同的营养机制实现生态位分化。功能均匀度值越高说明功能空间资源利用越充分,反之数值越低表明群落中物种功能性状越相似,资源利用率越低[63]。本研究中,功能均匀度集中分布在0.6~0.8之间,表明功能空间的资源利用较为充分,资源利用率较高。本研究结果表明环境因子与时间因子对功能均匀度和Simpson均匀度的影响不同。环境因子仅对Simpson均匀度的影响显著,对功能均匀度的影响不显著,仅解释了少量的功能均匀度变异,这与王松波等[64]的研究结果一致。而功能均匀度受时间因子的影响显著,可能与季节更替过程中气象条件、水体生态系统营养级结构的变化、浮游生物群落季节性演替过程中种间关系等因素有关。功能丰富度越高则反映群落中浮游植物功能性状占据的功能空间越大,群落对干扰的承受能力越强[60]。本研究中鹦鹉洲湿地和河道功能丰富度均表现为秋季远低于其他季节,且秋季湿地的功能丰富度显著高于外围河道,表明秋季是该研究区域浮游植物群落四季中抗干扰能力最弱的季节,且湿地的抗干扰能力显著高于河道。虽然功能丰富度的大小受物种丰富度的直接影响[64-65],但两者对环境的响应程度差异较大。变差分解结果表明,环境因子均显著解释了物种丰富度和功能丰富度的变异,但对功能丰富度的解释量(41.2%)显著高于对物种丰富度的解释量(16%)。这一结果表明功能丰富度对环境的响应显著大于物种丰富度,也就是说功能丰富度的环境敏感性比物种丰富度高,在一定程度上佐证了“环境筛选的是物种的功能,而不是物种本身”这一观点。
4 结论
2022年8月-2023年5月分四季对上海鹦鹉洲湿地与外围河道进行浮游植物群落采集,共确定浮游植物97种。整个调查期内,鹦鹉洲湿地与外围河道水体中浮游植物种类组成较为相似、平均生物量基本相当,但优势种结构和群落类型的演替规律存在较大差异。浮游植物群落显著的季节与空间差异主要由WT、DO、SD、NH3-N、DIP、Sal环境因子驱动,其中,Sal是区分湿地与河道浮游植物群落的关键因子,WT和NH3-N是区分四季浮游植物群落的关键因子。环境对物种生态位的筛选要强于对物种个体的筛选,研究区域内功能空间的资源利用较为充分,秋季是浮游植物群落抗干扰能力最弱的季节,湿地的抗干扰能力显著高于河道,功能丰富度的环境敏感性比物种丰富度高。浮游植物群落结构以及功能多样性等特征能够较为有效地反映生态交错带水环境的时空变化,未来可用于滨海区域生物多样性保护以及生态系统功能恢复效果的评估与预测。
致谢:感谢华东师范大学陈雪初教授在野外样品采集工作中的帮助与支持。