城市中小河道沉水植物盖度和生物量计算方法——基于上海市的调查分析*

张 群，刘家霖，朱 义，张春松，宋 晴，崔心红

(1:城市困难立地生态园林国家林业和草原局重点实验室，国家林业和草原局城市困难立地绿化造林国家创新联盟，上海城市困难立地绿化工程技术研究中心，上海市园林科学规划研究院，上海 200232) (2:生物多样性与生态工程教育部重点实验室，长江河口湿地生态系统国家野外科学观测研究站，复旦大学生命科学学院生物多样性科学研究所，上海200438)

城市河流作为区域重要的生态空间，提供水文调控、气候调节、污染消纳、休闲游憩、生物多样性维持等多种生态系统功能与服务[1-4]，在维护城市与区域生态安全格局方面发挥着重要作用[5]. 因受到长期高强度人为活动干扰，城市河流比自然河流面临更高的渠道化、硬质化、水质恶化等威胁[6]，水域生物种类和数量明显降低，生态系统健康及其维持机制遭到严重破坏[7-8]，城市河流生态修复已成为该领域研究和管理实践的热点问题[3]. 城市中小河道数量众多，处在汇水区接纳、滞留和调蓄径流排水和点面源污染的前缘地带[9-10]，是现阶段水污染治理和生态修复的主要对象[11]. 近十余年来，随着中国黑臭水体治理、污染源整治及河湖长制等一系列水环境改善的重大政策实施，取得了水质恶化趋势明显遏制、水质指标初步改善的重要成效. 为进一步提升及修复中小河道水环境与水生态，恢复沉水植物成为重要举措之一.

沉水植物作为水生态系统的重要组分和初级生产者之一，在物质循环和能量流动中发挥重要的调控作用[12-13]，明确其盖度和生物量，是开展水生态调查、评价和修复的基础性内容[14-17]. 沉水植物覆盖度及生物量的调查方法主要包括样方法、遥感解译法和水下声呐法等，国内外学者主要在自然河湖水域开展相关工作.传统上现场调查选择代表性沉水植物分布区，采用方形或圆形沉水植物采样器，在自然湖泊设置多个样带和样方估测生物量[18-21]；利用高分遥感卫片、水下回声探测等技术手段，与现场调查监测数据比对，构建了多种沉水植物盖度和生物量的估测模型[22-26]. 由于城市中小河流的沉水植物主要以人工栽植恢复的单一优势种群为主，并且沿岸区段化用地状况差异造成的空间异质性较大，高楼、乔木等沿岸光影对狭长水域数据产品的影响，制约了遥感、声学等技术估测准确性，在城市区域中小河流开展采样布设方式、采样次数及单位面积生物量的相关关系等方面的研究较少.

本研究以上海中小河道为研究对象，基于2019年对全市45条中小河道、130个调查区段研究基础上，选择分布有沉水植物分布的4条典型城市中小河道，开展采样数量、采样方法、盖度计算等调查方法对沉水植物生物量估算准确性的影响研究，主要解决如下技术问题：(1)对于城市中小河道沉水植物的单位面积生物量估算，如何确定最小采样次数和样点布设方法？(2)探索利用多光谱现场调查技术，如何进行城市中小河道沉水植物盖度估测.以期为城市河道沉水植被调查和评价提供科学依据，为城市河道生态系统健康精准化管理提供技术方法.

1 材料与方法

1.1 研究区域概况

上海市2019年全市河道46392条，总长度29862.14 km，平均长度为0.644 km，河网密度达4.71 km/km2，其中市、区、镇(乡)管河道3242条，其余93%以上的河道属于村管及其他河道[9]. 本研究的中小河道是指除了通航、水资源调度、排洪等主导功能以外的城市河道，一般通过设置涵闸进行片区内水利控制. 基于Google Earth影像目视判读、上海市水利分片划分及现场踏查，并根据2019年在全市选择45条中小河道，每条河道3～4个区段，每个调查区段为沿河长度50 m，共计130个调查区段，其中记录到沉水植物调查区段共计27个，仅占总调查区段的20.8%. 并测定沉水植物总盖度(多盖度等级法，采用目测法估计河道内所有沉水植物单位面积内植物茎叶对土地面积的覆盖程度，通常划分为5个等级和2个辅助等级[27])、电导率、pH、蓝绿藻、浊度等指标. 分布有沉水植物区段河道的浊度介于0～85.57 NTU之间，且沉水植物总盖度为COP1及以上(盖度介于5%～25%)的河道浊度小于30 NTU，浊度大于30 NTU的区段，其沉水植物仅零星分布，总盖度均小于5%. 于2020年9-10月及2021年5月(沉水植物生物量较多及较少的两季[18-19,28])，选取进木港(350 m)、三友河(640 m，两季)、北潮港区段(1000 m)、梅陇港区段(1100 m)4条以苦草(Vallisnerianatans)为优势种，竹节水松(Cabombacaroliniana)、黑藻(Hydrillaverticillata)、穗状狐尾藻(Myriophyllumspicatum)为伴生种的中小河道(图1)开展调查，不同河道内各沉水植物比例不同. 上述河道除雨季排涝间歇性人工降低水位外，水流大多处于静止状态. 2020年6月，4条河道5个区段水质调查指标(电导率、pH、浊度、叶绿素a、高锰酸盐指数、TN、TP)见表1.

图1 调查河道现场图(北潮港为南北方向区段，梅陇港为东西方向区段)Fig.1 Field map of river course (Beichao Port is a north-south section, Meilong Port is an east-west section)

表1 4条调查河道基本水化学指标Tab.1 Basic hydrochemical indexes of four rivers

1.2 最少采样数量计算

沉水植物在生长过程中，密度及生物量随时间的变化而变化，一般在每年夏季达到最大值[19]，且由于植株长度、密度等差异，导致单位面积的沉水植物生物量之间也存在差异[18]，为准确计算沉水植物单位面积生物量所需的最少采样数量，结合理论和实际操作两方面，利用功效分析法(power analysis)[29]及数据合并法(data binning)[30]进行计算.

图2 沉水植物采样器(1驱动柄、 2固定杆、3铰接座、4弹簧、 5驱动杆、6连杆、7铰接杆、 8伸缩杆、9采样钳、11手柄)Fig.2 Submerged plant sampling apparatus (1 driving handle, 2 fixed rod, 3 hinged seat, 4 spring, 5 driving rod, 6 connecting rod, 7 hinged rod, 8 telescopic rod, 9 sampling clamp, 11 handle)

本研究使用Cohen提出的适用于t检验的较大效应值基准[31]来计算最少采样次数，然后据此设计采样方案以估算沉水植物单位面积的平均生物量. 同时，在主要优势种为苦草或竹节水松的3条(或区段)中小河道内利用自制沉水植物采样器(图2)采集沉水植物. 采样器的其固定杆和采样钳通过伸缩杆铰接连接，采样面积固定，拉动驱动柄和手柄后采样钳闭合，切断沉水植物地上部分. 采集57～65个0.2 m× 0.2 m的样方，采样次数是估算中型湖泊沉水植物生物量样方的2～3倍[19-20]，分析采样次数与单位面积沉水植物生物量之间的关系，估算最少采样次数. 沉水植物离开水体后，失水较快，鲜重随时间变化差异较大[32]. 为减少差异，将夹起的沉水植物按种类分开后，适当沥干表面水分，迅速称重.

1.3 断面设置及采样方法

河道沉水植物采样点的设置，一般采用系统布点规则，即在多条河流交汇的河段随机布点[33]；也有按照河岸地形地貌在上中下游分别设置采样点[34]. 但是针对单条河道的沉水植物分布，由于河道长度及植物分布的差异，不同采样方法可能对生物量的估算产生影响. 因此，本研究比较等距样线法(根据河道长度等分，在中间的等分点设置样方)和多断面样点法(根据沉水植物长势将河道划分为10个断面，每个断面设置4～5个样方)估算同一条河道单位面积生物量的准确性. 在北潮港内先以等距离分隔采集43个样方，记为T1组；再根据河道内沉水植物的目视估算盖度，选定10个断面，每个断面采集3～5个样方，共计42个，记为T2组.

1.4 沉水植物盖度计算

城市河道沉水植物盖度是城市河道生态修复技术中的重要指标，但以往的调查方法与技术主要依靠样方取样后的测定与目测估算. 随着技术的发展，多光谱技术在大尺度已经广泛实现了地表不同生态系统植被盖度的高精度测算. 因此，受Song等研究的启发[35]，在调查河道上方2 m处，采用RedEgde多光谱相机拍摄照片，利用植被指数等空间计算叠加方法提取沉水植物分布，进而计算其单位面积的盖度.

由于植物均进行光合作用，所以各类绿色植物具有很相似的反射波谱特性，其特征是在可见光波段0.55 μm(绿光)附近有反射率为10%～20%的一个波峰. 充分利用该特性，以及绿光对水的穿透性较强，我们选择绿光照片进行沉水植物盖度识别与提取. 沉水植物盖度测算步骤如下：(1)利用Rededge-M多光谱相机(图3)拍摄5波段影像；(2)提取出其中绿光波段的数据，同时转换遥感影像的数据格式和地图投影，并利用ArcGIS提取波段数据信息，进行质量控制；(3)利用Arcmap工具，对绿光原始波段影像进行重分类，提取沉水植物的色彩阈值范围；(4)目视选择影像中沉水植物分布高、中、低的3个区域，生成3个对应矢量图层；(5)以矢量图层为边界提取重分类后的绿波段影像，并进行重采样(resample)；(6)重采样时仅保留沉水植物敏感的像素阈值范围并赋值为1，其余像素数值均赋值为0；(7)针对每个小样方，计算像素值为1的占整个样方像素点综合的比例，即为该样方沉水植物盖度.

如图3b案例所示，高盖度=16605/(16605+34695)=32.4%，中盖度=10026/(10026+31878)=23.9%，低盖度=5198/(5198+37324)=12.2%. 进而，计算得出该影像拍摄区域的沉水植物盖度为22.8%±8%.

图3 Rededge-M多光谱相机(a)；绿光波段照片中沉水植物分布高、中、低3个区域盖度计算案例示意图(b)Fig.3 Rededge-M multispectral camera (a); Schematic of coverage calculation for submerged plants in the green band images: high, medium, and low coverages (b)

图4 三友河(a)和梅陇港、进木港(b)取样次数与沉水植物生物量的关系Fig.4 The relationship between the sampling frequency and the biomass of submerged plants in Sanyou River (a)， Meilong and Jinmu Ports (b)

1.5 数据处理及分析

1.5.1 最小采样数 对于功效分析法，一般的功效阈值为0.8[36]，本研究设置功效值(置信度水平)为90%，采用0.8作为t检验的较大效应值基准[31]，在具有95%的准确性下，计算最少采样次数以估算河道内沉水植物单位面积的平均生物量.

数据合并法，首先计算采样次数对应的单位沉水植物生物量，即第10次为前10次(以采样顺序编号的次数)生物量总和的均值，第20次为前20次生物量总和的均值，依此类推；其次，将所得均值以4为单位分段获得14～16组数据，计算每组数据的均值及标准差；再次，利用Sigmaplot作图软件，以单位面积生物量作为Y轴，取样次数作为X轴画出散点图(图4)；最后，添加所有数据的均值作为参考线，探索取样次数与单位面积生物量之间的关系，计算最少取样数.

1.5.2 不同采样方法的差异比较 对获取于北潮港的两组数据合并后，进行正态性及方差齐性检验，发现数据不符合正态分布，但方差满足齐次性. 因此，为检验T1与T2的组间差异，采用威尔科克森(wilcox)秩和检验进行分析.

所有的分析工作在R 4.0.2软件中进行[37].

2 结果与分析

2.1 中小河道采样数及采样方法的确定

采用功效分析法，在90%置信区间下，以效应值0.8为基准值，n≈34，即样方取34能较为准确地代表整个河道的沉水植物生物量. 以徐汇区三友河为例，该河道主要分布人工种植的苦草，并伴生少量的自然生长的黑藻及金鱼藻(Ceratophyllumdemersum). 通过在0.2 m × 0.2 m地块上连续采集沉水植物地上生物量，建立采样次数均值生物量与所有样本生物量之间的关系，以表征小样方的空间异质性[38]. 结果表明，当取样次数达到28～32次后的预测生物量(图4a箭头所示)与该河道沉水植物单位面积生物量均值较为接近，该取样次数可较为准确地估算河道沉水植物生物量. 以相同方法，对梅陇港和进木港的最小采样数进行计算，结果如表2、图4b所示. 最少取样数介于32～36次之间. 结合能效分析方法的结果(34次)，因此建议对上海市单条中小河道沉水植物生物量估算的最少取样次数为30～35次. 同时，研究发现等距样线法和多断面样点法对估算河道单位面积生物量的准确性差异不显著(P=0.993)(图5).

表2 3条河道最少取样次数Tab.2 Minimum sampling numbers of three rivers

图5 北潮港2种取样方法的分组差异Fig.5 Group differences between the two sampling methods in Beichao Port

2.2 中小河道沉水植物的盖度

与沉水植物取样同期同步，对三友河沉水植物盖度进行光谱学测算. 如图6a所示，在3 km长的河道上，设置了9个拍摄点，有效照片12张，拍摄高度约2 m，按分辨率折算单张照片覆盖面积约7.5～9.5 m2，总体覆盖面积接近100 m2. 该面积已超出人为取样可覆盖面积和可目视测算的有效面积，具有较高的空间代表性. 此外，在计算过程中，由于水面反射或水体有杂质，该研究在单张照片内设置多个有代表性且便于计算的小样方，提取9张受拍摄角度和光影反射相对较小且相对清晰反映沉水植物的多光谱照片，对9份绿光波段照片进行解析，人工监视下目视识别并提取27个小样框，计算小样框内沉水植物盖度的相对比例(去量纲化)，并进行综合统计. 如图6b所示，三友河沉水植物盖度在44%水平上，标准差约25%. 与现场人工评估的沉水植物盖度结果相比(约50%)，光谱影像计算与现场目视评估结果相近.

图6 9个拍摄断面的位置(a)和进木港沉水植物盖度影像学提取结果(b)Fig.6 Nine locations for imaging sections (a) and image extraction results of submerged plant coverage in Jinmu Port (b)

3 讨论

3.1 河道长宽对单位面积生物量估算的影响

城市河道生物量的计算应首先对河道长度、驳岸类型、河道内沉水植物的种类及分布进行全面的调查和记录，再采用随机或等距方式设置断面及样点. 上海市2019年全市河道46392条，总长度29862.14 km，平均长度为0.644 km，其中城市中小河道仅包含在市、区、镇(乡)管河道的3242条之内，其平均长度均在1 km之内，且不具有通航功能的中小河道河宽一般小于8 m(调查数据，未发表)，河宽差异较小. 因此对于上海绝大多数的中小河道，长度和宽度对取样次数的影响较小，单条河道30～35次的取样次数适用性较高. 但是，对于特殊的长度超过1 km的河道中小河道，该取样次数的适用性还需要进一步研究.

3.2 自然及人工水域沉水植物生物量研究方法的差异

自然水域例如湖泊的底部较为平坦，沉水植物分布连续且均匀，通过随机的10次采样，即能较为准确地估算单种沉水植物生物量[18]. 而对于多种沉水植物生物量的动态研究，则需设置多个断面，根据断面的长度等距设置2～8个样方，共计30个样方，亦能准确估算[19]. 但上海城市中小河道属平原河网河流结构，受自然条件和城市化的双重影响，其中的井网结构是人为构筑的河流结构[11]，河道底部除底泥外，还掺杂有少量的石块、水泥、建筑垃圾等硬质材料. 对上海4条河道调查时发现，硬质材料一般分布于河道两侧，呈不连续性和随机性，可能是河道施工、河岸改造过程中遗留的少量建筑废弃物，影响了沉水植物的定植与存活. 因此，利用人工样方法调查城市中小河道沉水植物的生物量，增加采样次数可降低异质性带来的误差. 但是，随着采样次数、河道长度、河网密度的增加，导致河道调查的人力、时间等投入提高，为有效获取及监测反映河道沉水植物的生物量指标数据，需要权衡投入与产出、成本与效益之间的关系.在明确调查目标的基础上，结合研究区域的河网遥感图像及水利片区划分，可选取不同用地类型、交汇数量多、具有明显差异的能代表河网特征的河道，再根据现场状况选取适宜的断面及样点进行调查和采样.

3.3 沉水植物优势种覆盖度的多光谱监测方法

利用卫星数据对水体水生植被的监测主要集中于大型的浅水湖泊[39-41]，且部分方法基于全年遥感数据的统计分析，不适用于实时监测.而对于中小型湖泊及水库水生植被类群的遥感监测，可利用高光谱分辨率的Sentinel-2数据进行统计分析[24]，但沉水植物的冠层深度及水环境因素会导致遥感分类出现偏差. 水声学技术对浅水湖泊沉水植物的调查也取得较多的进展及成果[42-44]，还可通过回声探测技术研究浅水湖泊沉水植物的分布[26]，但是这些方法需要大量的专业设备，并对水体深度有一定的要求.

城市中小河道宽度较窄，分辨率较低的遥感影像难于识别和区分. 以上海为例，城市河道水体的平均深度介于1.5～2.0 m间(未发表数据)，因此基于卫星或无人机的深空和高空遥感技术以及目前在自然河湖地带采用的水声学技术等监测方法较难应用于城市中小河道沉水植物的监测及评估[45]. 本研究提出基于河面上方的人工控制的主动式多光谱相机拍摄，可以充分地发挥多光谱技术对沉水植物优势种覆盖度的监测能力，并有效地避免天气变化、云层覆盖度、水深异质性等陆面本底状况对深空和高空遥测技术的阻碍[46]. 但是，本研究提出的人工控制的主动式多光谱相机监测技术尚处于开发状态，还需要辅助以代表性断面的人工取样方法予以校正，后续相关研究需要针对以下两点进行调研与探索，以更好地服务于城市中小河道沉水植物调查技术的完善与发展：(1)如何从多光谱照片中有针对性地对优势种沉水植物覆盖度和非优势沉水植物覆盖度进行区分？(2)如何调整R/G/B/NIR波段计算公式以更好地表征高浊度(>30 NTU)下城市河道优势沉水植物覆盖度？

3.4 中小河道沉水植物的保护和管理建议

水生植物是水体生态系统的初级生产者，其在生长期可吸收水体中大量的营养盐，对水体有较好的净化作用[47]. 但是沉水植物过度繁殖可导致植物衰亡分解产生大量的氮、磷等营养物，进而引起水质的恶化及二次污染[48-49]，因此当河道内沉水植物大量繁殖后，需对沉水植物的生物量进行控制. 适宜的控制措施能去除植物体中的部分营养物质，促进沉水植物的恢复和生长，维持水体的正常水质[43]. 但是不同沉水植物的繁殖和扩散能力不同，建议在河道管理过程中，主要控制和打捞黑藻、穗状狐尾藻(Myriophyllumspicatum)、竹节水松(Cabombacaroliniana)、菹草(Potamogetoncrispus)等依赖茎段散布营养繁殖的种类[18]. 此外，不同时期控制沉水植物生物量，对水体的影响存在差异，在生长期收割沉水植物可降低水体的营养盐水平[50]. 因此，对中小河道沉水植物的收割管理，应权衡种类、强度及时间阶段的相互影响，制定科学合理的措施及方案.

随着城镇黑臭水体治理、河湖长制深入推进，准确快速地掌握沉水植物特征对水生态修复具有重要指导价值.由于目前缺乏在城镇河道人工恢复的沉水植物生长响应与机制深入研究，后续应明确其主要影响要素和典型类型，通过不断累积城镇典型片区的河道沉水植物调查数据，提高单位面积生物量估算和近地无人机多光谱测定盖度的准确性，构建城镇典型片区河道沉水植物分布与生物量估算模型，从而建立快速、准确、高效的城镇河道沉水植物调查方法和技术.

4 结论

1)通过对上海中小河道沉水植物生物量和盖度的调查研究，建议在单位长度为1000 m的城镇河道，开展沉水植物单位面积生物量估算中，最少需取样30～35次，可随机或等距设置采样点.

2)利用多光谱相机进行河道沉水植物总盖度估算，通过5～8次的5波段影像拍摄，平均后可获得整条河道的沉水植物盖度. 利用单位面积生物量及总盖度，可准确估算河道沉水植物生物量.

致谢：感谢复旦大学邱士允博士、王新新博士对本文写作的指导，感谢上海应用技术大学高涛、万一砾同学以及杨震林在河道调查数据获取过程中提供的帮助.