江苏省天福国家湿地公园鸟类栖息地植物群落特征对鸻鹬类水鸟的影响研究

邵君学，胡昕欣，李欣，周婷婷，冯育青，刘宁，周敏军，谢冬,＊

(1.南京林业大学生物与环境学院，江苏 南京 210037；2.中华人民共和国国际湿地公约履约办公室 国家林业与草原局，北京 100714；3.太湖湿地生态系统国家定位观测站，江苏 苏州 215000；4.苏州市湿地保护管理站，江苏 苏州 215128；5.天福国家湿地公园保护管理中心，江苏 苏州 215000)

鸻鹬类(waders)水鸟是湿地生态系统鸟类多样性的重要组成部分，是“东亚-澳大利西亚”候鸟迁徙线上的重要水鸟类群。由于气候变化和人类活动的影响，大量湿地已经丧失和退化，这对水鸟多样性构成了严重威胁[1-5]。植物群落作为水鸟栖息地三大生境因子(植物、水、食物)之一，在鸻鹬类水鸟生境的恢复和保护工作中至关重要，因此植物群落特征对鸻鹬类水鸟组成和分布的影响机制日益受到关注。有关鸻鹬类水鸟的生境需求、保护和恢复的研究已有不少，如葛振鸣等[6]在长江口杭州湾对鸻形目鸟类群落生境选择的研究指出，海三棱藨草(Bolboschoenoplectus mariqueter(Tang &F.T.Wang) Tatanov)覆盖比例对鸻形目鸟类群落有正向影响，而芦苇(Phragmites communis(Cav.) Trin.ex Steud.)、互花米草(Spartina alternifloraLoisel.)密植区域则不利于鸻形目鸟类利用。胡越凯等[7]的研究指出，海三棱藨草生境破碎化趋势使得鸻鹬类鸟类群落有向其他生境转移的趋势。但有关这些研究，(1) 多是基于自然滨海湿地及河口滩涂湿地，而没有关注以池塘为典型生境类型的人工恢复湿地；(2)所讨论的植物群落特征大多局限于具体植物种类，而没有关注整体植物群落特征，以及在时间维度上伴随植物群落演替和人为干扰双重影响下植物群落的变化对鸻鹬类水鸟的影响。(3)仅关注鸻鹬类水鸟各种群特征，未从功能性状的水平上关注其物种集团的组成和分布。集团是指一群具有相似的方式探索和利用同类环境资源的物种集合和群落功能单位[8]。作为群落功能群的重要类型，通过鸟类群落集团的划分有助于从功能群的角度探究鸟类物种在群落中的功能定位[9]。这一概念早已成为鸟类群落功能群研究的常用方法[10-11]。但使用功能性状作为划分集团的变量还较少，本研究以鸻鹬类水鸟性状作为划分鸟类集团的变量，探究植物群落特征的季度变化对不同鸻鹬类鸟类集团的差异性影响，从而为鸻鹬类水鸟栖息地恢复和管理提供理论依据。

1 研究区域

天福国家湿地公园位于江苏省苏州市昆山市花桥经济开发区境内，地处太湖流域吴淞江水系中游，距离吴淞江干流5 km，地理坐标范围：121°05′22.32″～121°07′44.08″ E，31°19′18.28″～31°21′13.53″ N，湿地公园批复总面积779.54 hm2。年平均气温17.6℃；年平均降水量1 200.4 mm，年平均日照时间1 789.2 h。气候温和湿润，四季分明，光照充足，雨量充沛[12]。

鸟类栖息地恢复区为典型人工修复湿地，2017年改造恢复，遵循低影响管理原则，内部草甸和开阔水面为野生鸟类提供了栖息环境和食物资源，完整的生态系统服务功能让鸟类栖息地恢复区成为天福国家湿地公园的“一道独特的风景线”。

2 研究方法

2.1 植物群落样地调查

样地设置：根据研究区实际土地利用情况，划分45 块极少有人类活动用地进行调查。据方精云等[13]关于植物野外清查的研究方法，在每个样地内随机选取3个1 m×1 m 样方，共计135个样方，样方框内用绳索分割并固定成25个20 cm×20 cm的小方格以便于估算样方内植物盖度和数量。于2019年7月至2020年7月历经5个季度对研究区的植被群落进行调查。

调查指标选取和测定：详细记录每个样方内草本植物的种类、数量、盖度以及植物群落高度，使用经纬仪记录样方所在地海拔、经纬度等指标。数量与盖度采用目测估算；高度：皮尺测量。

2.2 鸟类调查及鸟类性状测量

鸟类调查：采用“样带法”方式于2019年7月至8月，2020年7月至8月进行每月1 次的鸟类调查，每次持续2 天。鸟类观测选择晴朗少风的天气，并根据鸟类活动习性，在晨昏2个时间段（7:00—10:00；16:00—18:00）开展调查。根据《中国鸟类分类分布名录》（第3 版）[14]进行识别、记录，具体调查方法为：对样线或样点100 m 范围内的样带，通过双筒望远镜、单筒望远镜、专业单反相机对鸟类进行观察和辨别，并结合《江苏鸟类》等手册进行野外辨识，记录种类、数量，为控制调查系统误差，每次调查人员保持一致。

鸟类性状测量：鸟类性状数据采集自中国科学院动物研究所国家动物标本资源库，具体方法：取各鸻鹬类鸟种标本5 份，分别使用游标卡尺测量不同的性状（嘴裂长、翼长、跗跖长)。

2.3 数据处理与分析

2.3.1 植被群落特征统计及分析 分别统计植物群落高度（PCH）及调查区植被覆盖度（VC），并通过Pearson 相关性检验分析两者的相关关系，计算二者比值（PCH/VC），分析季度变化规律。植被群落高度：每个样方中所有植株高度的均值；植被覆盖度：样方中植株的投影面积占样方面积的百分数的均值。

2.3.2 鸟类数据分析及鸻鹬类水鸟集团划分 按季度统计鸟类个体总数（TNB）及鸻鹬类水鸟个体总数（TNW），并通过Pearson 相关性检验分析两者的相关关系。计算两者比值（TNW/TNB），结合植物群落特征PCH/VC，建立PCH/VC 和TNW/TNB的一元线性回归模型。

通过K-means 聚类方法分别将鸻鹬类水鸟翼长、嘴裂长、跗跖长3种性状数据划分成高-中-低3个水平，共获得9个鸟类集团变量，再对各鸟种就9个鸟类集团变量进行评估赋值，符合变量要求的赋值1，其余赋值0[15]。接着对各鸟种的得分情况进行主成分分析，截取特征值 ≥1 的主成分进行降维分析。在各主成分中选取载荷量绝对值 ≥0.6 的变量作为该主成分的显著变量[16]。根据显著变量的构成确定该成分主要反应的鸟类集团信息，将获得的若干主成分得分为变量两两为轴建立直角坐标系，绘制散点图，再对各主成分得分进行系统聚类分析，获得各鸟种主成分得分的平均联接（组间）树状图，散点图与树状图相结合划分出最终的鸟类集团[17]。最后分析不同鸟类集团个体数在鸻鹬类水鸟类个体总数(TNW)中的占比(即不同鸟类集团对TNW 的贡献值)，可得出植物群落特征PCH/VC对不同鸟类集团影响的差异性。

此外，依据上海气象局（http://soweather.com/）资料，定义3—5月为春季、6—8月为夏季、9—11月为秋季、12月至翌年2月为冬季[18]。本研究所获取数据均以IBM SPSS Statistics 20、Excel 2019软件进行数据处理和数据分析。

3 结果与分析

3.1 植物群落特征

在鸻鹬类水鸟频繁活动的水域周边共计12个样地、36个样方中，共记录到草本植物77种，分属31 科72 属；其中禾本科代表植物有：雀稗(Paspalum thunbergiiKunth ex Steud.)；菊科代表植物有：蒲公英(Taraxacum mongolicumHand.-Mazz.)。各季度植物群落高度（PCH）分别是32.32、48.81、25.44、24.84、32.39 cm；各季度调查区植被覆盖度（VC）分别是0.62、1.08、0.62、0.39、0.51；各季度PCH/VC分别是52.15、45.37、41.22、63.03、63.19（表1）。Pearson 相关性分析结果显示植物群落高度和植被覆盖度各季度的分布在0.05 水平（双侧）上显著相关（r= 0.916，p= 0.029），均在2019年秋季达到峰值，2020年春季达到最低（图1a、图1b）。PCH/VC比值在2019年夏季至冬季呈下降趋势，在2020年春季至夏季呈上升趋势(图1c)。

3.2 鸟类群落特征

调查期间共记录到鸟类119种，总个体数3 401 只；其中鸻鹬类水鸟共18种，个体数230只，代表种有金眶鸻(Charadrius dubiusScopoli.)、黑翅长脚鹬(Himantopus himantopusLinnaeus.)。各季度分别记录到鸟类个体总数(TNB) 及鸻鹬类水鸟个体总数(TNW) 分别为154、13；776、27；742、15；1 104、124；625、51，表明二者在各季度变化趋势一致(图1d、图1e)，但无显著相关性(r= 0.731，p= 0.160)。各季度二者的比值TNW/TNB 分别是0.08；0.03；0.02；0.11；0.08（表2），比值TNW/TNB 在2019年夏季至冬季呈下降趋势，在2019年冬季至2020年春季呈上升趋势(图1f)。

3.3 植物群落特征和鸻鹬类水鸟群落特征的关系

以植物群落特征(PCH/VC)为自变量，鸻鹬类水鸟群落特征(TNW/TNB)为因变量构建一元线性回归模型，发现模型汇总调整后R2= 0.764，p<0.05；自变量回归系数0.003，p<0.05；常量−0.116。最终构建回归方程：

式中：Y：鸻鹬类水鸟群落特征(TNW/TNB)；X：植物群落特征(PCH/VC)。

3.4 鸻鹬类水鸟集团划分

主成分分析显示，前3个主成分(PC1～PC3)的特征值大于1，累积贡献率达80.49%，因此选取PC1～PC3 进行下一步分析。PC1 的显著变量是高翼长、低翼长、高嘴裂长、低嘴裂长、高跗跖长、低跗跖长；PC2 的显著变量是中翼长、高跗跖长、中跗跖长；PC3 的显著变量是低跗跖长(表3)。对PC1～PC3 得分构建坐标系(图2a～c)及进一步的系统聚类结果(图2d)显示将组间联接平方距离定为7 时可合理地划分出4个鸟类集团。结合PC1～PC3 主成分的显著变量，将4个鸟类集团定义为：

图2 鸻鹬类水鸟PCA 分布及聚类分析Fig.2 Distribution in PCA axis and cluster analysis of waders species

表3 基于9个鸻鹬类水鸟集团变量的PCA 结果Table 3 PCA analysis results based on 9 waders guild variables

(1) 高翼长、高嘴裂长、高跗跖长集团，包括3种鸻鹬类水鸟，有黑翅长脚鹬(HimantopushimantopusLinnaeus.)、反嘴鹬(Recurvirostra avosettaLinnaeus.)、青脚鹬(Tringa nebulariaGunnerus.)。

(2) 低翼长、低嘴裂长、低跗跖长集团，包括6种鸻鹬类水鸟，有长趾滨鹬(Calidris subminutaMiddendorff.) 、矶鹬(Actitis hypoleucosLinnaeus.)、金眶鸻(Charadrius dubiusScopoli.)、青脚滨鹬(Calidris temminckiiLeisler.)、黑腹滨鹬(Calidris alpinaLinnaeus.)、红颈滨鹬(Calidris ruficollisPallas.)。

(3) 中翼长、高嘴裂长、低跗跖长集团，包括5种鸻鹬类水鸟，有针尾沙锥(Gallinago stenuraBonaparte.)、大沙锥(Gallinago megalaSwinhoe.)、扇尾沙锥(Gallinago gallinagoLinnaeus.)、白腰草鹬(Tringa ochropusLinnaeus.)、尖尾滨鹬(Calidris acuminataHorsfield.)。

(4) 中翼长、中嘴裂长、中跗跖长集团，包括4种鸻鹬类水鸟，有彩鹬(Rostratula benghalensisLinnaeus.)、泽鹬(Tringa stagnatilisBechstein.)、林鹬(Tringa glareolaLinnaeus.) 、鹤鹬(Tringa erythropusPallas.)。

3.5 植物群落特征对不同鸻鹬类水鸟集团的影响

通过植物群落特征和鸻鹬类水鸟之间的数学关系，结合4个鸻鹬类鸟类集团个体数在鸻鹬类水鸟个体总数(TNW) 中的占比可以看出：在研究区，植物群落特征的演替变化对4个鸻鹬类水鸟类集团具有差异性的影响。植物群落特征的季度变化对集团4 影响最大，达34%；对集团2 和集团3 的影响次之，达25%、27%；对集团1 的影响相对较小，仅14%。

4 讨论

4.1 植物群落特征季度变化和鸻鹬类水鸟群落动态

鸻鹬生境的优劣与植被覆盖度和植物群落高度相关[19]。据唐承佳等[20]的研究，湿地植被覆盖度大于60%时，鸻鹬基本不出现，当植被覆盖度为10%～20%时，数量最多。而不同的植物群落高度则会吸引着不同性状的鸻鹬类水鸟，根据《野生生物湿地管理》[21]，低矮丛生植被更加吸引大型鸻鹬类，湿润的低矮草地则更加吸引中小型鸻鹬类前来觅食和栖息。研究结果表明，植物群落特征(PCH/VC)变化趋势和鸻鹬类水鸟群落特征(TNW/TNB)的变化趋势保持一致(图1c，1f)，且都是冬季最低，春季最高。这一结果在和前人的研究结果相互印证的基础上，可定量看出PCH 和VC 是如何影响鸻鹬类以及整体鸟类季度分布。鸟类对栖息地的选择的主要影响因素有水、食物和隐蔽性[22]。由于集团4 鸻鹬类水鸟对植物群落特征的季度变化的敏感性相对其他鸟类集团最大，因此推测，研究区目前更具有提供并调控中等翼长嘴裂长和跗跖长的鸻鹬类水鸟生境能力。集团2 和集团3 的敏感性处于中等，推测研究区拥有适合这两个集团鸻鹬类水鸟栖息的生境，但需加强这种生境的维护和管理。集团1 的突出性状特征在于高嘴裂长、高跗跖，这种性状特征更有利于获取较深水域里的食物，推测集团1 鸻鹬类水鸟在研究区对水深的敏感性可能大于对植物群落季度变化的敏感性。

本研究仅从植物群落特征的季度变化角度分析和探讨了植物群落特征对鸻鹬类的影响机制，还不能完全解释在食物资源相对缺乏的春季，鸻鹬类个体总数达到峰值的原因。由此也可进一步推测研究区的植物群落特征并非通过提供食物进而影响鸻鹬类季度动态变化。

4.2 栖息地管理措施对植物群落和鸟类群落季度变化的影响

草地植物群落特征季度变化一方面遵循植物群落演替规律，另一方面，也受到人类活动的干扰。通过走访调查，自2020年1月开始，天福国家湿地公园开始刈割研究区高茎草本植物地上部分。显然人为干扰显著影响植物群落高度，同时大量植物叶片被人为地从湿地生态系统中移出，这对植被覆盖度也产生了重要影响。两方面原因都解释了冬季和春季PCH 和VC 都处于变化曲线的波谷位置(图1a，1b)。同时结果显示，冬季之后，春季的植物群落高度和植被覆盖度的比值显著增加，说明管理方目前采取的管理措施对植物群落高度的抑制作用小于对植被覆盖度的抑制作用。一定程度上增加了鸻鹬类在整体鸟类中的比重。

根据式（1）可推测，草地植物群落正常演替不受人为活动干扰，则VC 会显著增加，同时PCH在短期内不会增长(会随着季相更替规律变化)，由此导致植物群落高度和植被覆盖度的比值减小，意味着鸻鹬类在鸟类中的占比有极大可能性会降低，因此管理工作须对植物群落的自然演替有更加充分的认识和警觉。

5 结论

通过实地调查和数据分析，将两个植物群落特征(植物群落高度、植被覆盖度)转换成比值，结合鸻鹬类水鸟个体总数在整体鸟类个体总数上的占比，构建回归模型，两者显著正相关；并进一步划分鸻鹬类鸟类集团，得出不同性状的鸻鹬类对植物群落特征的季度变化响应程度有差异，其中集团4 即中翼长、中嘴裂长、中跗跖长的鸻鹬类水鸟受植物群落特征的季度变化影响大于其他鸻鹬类鸟类集团；目前无区别化的栖息地管理措施须结合植被群落演替规律及不同性状鸻鹬类水鸟的生境需求而改进。

致谢：承蒙苏州市湿地保护管理站及天福国家湿地公园保护管理中心大力支持；在鸟类功能性状数据研究方面，感谢西交利物浦大学李黎及中国科学院动物研究所国家动物标本资源库贺鹏的悉心指导和协助。再次致以谢忱！