黄河口湿地时空变化过程及其主控因素

李燕妮，李鹏，吴晓，毕乃双，王厚杰

中国海洋大学海洋地球科学学院，海底科学和探测技术教育部重点实验室，青岛 266100

湿地与海洋、森林并称地球的三大生态系统[1]，是珍贵的自然资源，具有不可替代的综合功能。湿地在蓄洪抗旱、促淤造陆、调节气候、降解污染等方面具有极为重要的作用，被称为“地球之肾”。湿地植物通过光合作用固定了大气中的二氧化碳，具有较强的存储二氧化碳和固碳能力，是重要的碳储库[2-3]，在全球碳循环中发挥了重要作用，其碳汇功能对于减缓全球变暖进程具有重要意义。

黄河以水少沙多、含沙量高而闻名于世，自1855 年改道渤海以来，巨量的河流泥沙塑造了陆地面积约为5000 余平方千米的现代黄河三角洲。黄河三角洲湿地地理位置优越，生态类型独特，是暖温带保存最完整、最典型、最年轻的湿地生态系统[4]。受到黄河频繁改道与入海水沙快速变化的影响，黄河三角洲河口湿地的分布格局、植被覆盖与群落结构等均随之快速变化。自1976 年黄河三角洲改道清水沟入海以来，现行黄河三角洲河口湿地呈现阶段性演化的特征，发育演变不断调整。1996 年黄河改道清8 汊以来，原清水沟废弃河道遭受快速侵蚀，由之前的向海淤进状态快速逆转为向陆蚀退状态，而现行清8 河口由于入海泥沙的堆积而不断淤积造陆。

因此，现行黄河口叶瓣兼具行水河口和废弃河口两类特征，为对比研究河口三角洲湿地变化提供了良好的试验场。研究黄河口湿地的时空变化过程，揭示其主要控制因素，阐明湿地植被时空分布特征，深化多因素影响下河口湿地发育演变的科学认识，可为黄河三角洲湿地管理与保护的相关政策制定提供有效的科学依据。

1 研究区概况

研究区位于现代黄河三角洲东部。根据《国际湿地公约》中的湿地分类系统[5]，结合本文研究需要，将潮间带滩涂湿地及草甸沼泽作为湿地面积变化的主要研究对象，盐田、库塘、养殖池塘等人工湿地主要受人为因素的影响，对河流输入变化和海洋动力环境变化不敏感，本文不予重点分析。

黄河自1976 年由刁口河道改道清水沟流路入海，清水沟流路行河以来先后经历了淤滩成槽、溯源冲刷和溯源淤积的阶段性演变过程。为有利于防洪安全和保障胜利油田的基础设施建设，1996 年5 月进行了一次小型人工改道，即清8 改汊入海工程，将河道向北改道至清8 航道。此后，清水沟河道被人工堤坝完全切断，与主干道断开，清水沟流路完全废弃，废弃河道长约20 km[6]。

本研究系统收集了1976—2020 年期间覆盖现行黄河三角洲叶瓣区域的Landsat 系列卫星遥感影像资料和黄河利津站水沙资料数据，对黄河口区域湿地演变的时空特征及影响因素进行了研究。为了对比研究清8 汊（行水）和清水沟（废弃）不同流路湿地演化的差异及主控因素，将研究区以黄河改道处为基准点O（图1），以清8 汊叶瓣与废弃清水沟叶瓣间“人”字区交点为固定点B，向南北两侧各呈60°进行分割，形成清8 汊叶瓣湿地区（Ⅰ区）及废弃清水沟叶瓣湿地区（Ⅱ区）。

图1 现行黄河三角洲遥感影像及湿地研究分区示意图2020 年10 月24 日 Landsat8 OLI 标准假彩色影像，UTM 投影。Fig.1 Image of the active Yellow River Delta and wetlandBase map: Landsat 8 false color image on October 24, 2020, UTM projection.

2 数据与方法

2.1 数据来源

本文所采用的数据包括卫星遥感数据及黄河入海水沙资料两部分。其中，卫星遥感数据为Landsat 系列卫星影像，数据来源于地理空间数据云（http://www.gscloud.cn/），选取了1976—2020 年每年影像获取时间相近、研究区云量少、植被状态较好的37 幅遥感影像（表1）。首先，通过ENVI 软件对遥感影像进行影像配准、辐射定标、大气校正等预处理；然后，在此基础上提取湿地面积、植被覆盖度、植被覆盖类型等信息。

表1 Landsat 遥感影像数据信息（1976—2020）Table 1 Information of Landsat imagery (1976—2020)

黄河入海水沙数据来源于黄河水利委员会公开发布的黄河利津站水沙观测资料，包括1976—2019 年逐年径流量、输沙量等数据。最近70 年来黄河入海水沙的变化过程及控制机制、河口沉积效应等，可参见Wang[7-8]和Wu[6,9]等研究。

2.2 湿地面积信息提取

湿地面积提取采用分层次分类的方法[10]，建立解译标志（表2），利用监督分类提取滩涂湿地面积，利用归一化植被指数（normal difference vegetation index, NDVI）提取草甸沼泽面积，利用改进的归一化 水 体 指 数（modified normal difference water index,MNDWI）提取人工水体面积。

2.2.1 滩涂及沼泽湿地信息提取

本研究区内的滩涂主要为潮间带滩涂，滩涂湿地面积较大且与除水体外其他地物在颜色方面有明显差别，可在ENVI 中采用最大似然法进行监督分类提取，在ArcGIS 中结合目视解译对异常区域剔除。

研究区内的沼泽湿地覆有植被，以草本沼泽和沼泽化草甸为主，可采用NDVI 进行提取，沼泽湿地的NDVI 值为大于0.2 的正值。

NDVI 值可指示影像像元中绿色植被的数量，较高的NDVI 值对应较多的绿色植被，可用来监测区域的植被覆盖变化。通过比值处理，NDVI 可以部分消除与太阳高度角、卫星观测角、地形、云及阴影和大气有关的条件变化影响，增强了对植被的响应能力[11]。

?

计算公式如下：

其中， ρ(NIR)表示遥感影像中的近红外波段的反射值， ρ(RED)表示红光波段的反射值，在Landsat 系列影像中各波段对应如表3 所示。

表3 不同传感器中NDVI 各波段对应值Table 3 Band of NDVI in different Landsat sensors

2.2.2 人工湿地信息提取

研究区域内存在大量的人工湿地，如盐田、库塘、养殖池等，其含水量较多，可通过MNDWI 进行提取。MNDWI 表达式为：

其中， ρ(MIR)表示遥感影像中的中红外波段的反射值， ρ(GREEN)表示绿光波段的反射值，在Landsat 系列影像中各波段对应如表4 所示。

表4 不同传感器中MNDWI 各波段对应值Table 4 Band of MNDWI in different Landsat sensors

2.2.3 非湿地区域信息剔除

研究区内的非湿地区主要为北部的孤东油田。在1987 年孤东海堤建成后，该区域不再被海水周期性覆盖，因此可归为非湿地区域。由于孤东油田形状比较规则，可直接通过目视解译将其进行裁剪剔除。

2.3 植被覆盖度信息提取

像元二分模型是一种简单实用、应用广泛的遥感估算模型[12]。该模型假设单个影像像元内的地表由无植被覆盖与有植被覆盖组成，即纯土壤与纯植被。通过遥感传感器观测到的光谱信息也由这两个组分因子线性加权合成，各组分因子的权重是各自面积在像元中所占的比例，其中植被的权重可表征该像元的植被覆盖度[13]。

根据像元二分模型，在一个像元中NDVI 值可以表达为由无植被覆盖部分所贡献的信息NDVIsoil和绿色植被所贡献的信息 NDVIveg两部分组成，表达式为：

其中， NDVIveg和 NDVIsoil分别为纯植被像素和纯裸地像素的NDVI 值。纯裸土区域的NDVI 值理论上应该为0 或接近于0，纯植被区域的NDVI 值应该接近于1，但在实际环境中由于受各种因素影响，NDVIsoil和 NDVIveg的值会因地理位置和植被类型不同而改变，因此一般情况下无固定值，需结合本区域NDVI 概率分布值而确定。结合前人在黄河三角洲地区的研究，本文将NDVI 统计数据中频率为5%的NDVI 值作为 NDVIsoil，把植被覆盖区的内取频率95%的NDVI 值作为 NDVIveg[11]。

利用ENVI 软件的快速统计功能可提取各时相影像数据的NDVI 统计值，获得频率在5% 和95%时对应的 NDVImin及 NDVImax，对各时相影像进行波段运算并对异常值进行处理。根据现行黄河三角洲叶瓣植被覆盖度的变化，将植被大致划分为5 个等级（表5）。

表5 植被覆盖度等级划分标准Table 5 Rank of fractional vegetation cover

?

2.4 植被覆盖类型信息提取

本文以黄河三角洲湿地较为典型的四种植物：芦苇、互花米草、碱蓬和柽柳为研究对象，选择了1980、1995、2000、2004、2011、2016、2020 年共7 期影像进行分类研究。首先结合谷歌地球历史图像及图新地球软件的高分辨图像进行先验学习，建立湿地植被的解译标志（表6）；在此基础上，使用最大似然法和支持向量机两种监督分类方法进行分类，利用地面真实感兴趣区或真实分类图像进行混淆矩阵分析生成精度评价表，选择总体精度和Kappa系数都较高的支持向量机分类结果。

3 结果与讨论

3.1 现行黄河三角洲叶瓣湿地时空变化特征及影响因素

根据卫星遥感影像提取现行黄河三角洲湿地面积信息，建立了自1976 年黄河改道清水沟以来的河口不同类型湿地面积的时间序列数据（图2），其时序变化总体分为以下4 个阶段（图3）。

图2 现行黄河三角洲叶瓣湿地面积时序变化Fig.2 Temporal variation of the wetland area in the active Yellow River Delta lobe

（1）快速淤进阶段（1976—1981 年）：湿地面积整体呈现快速增长趋势，平均增长速率约为77.42 km2/a。湿地面积与黄河累积入海水沙量均呈现显著正相

关，决定系数R2分别为0.92 和0.86（图3）。在黄河改道初期，入海流路尚未成槽，频繁分汊摆动，河流漫流沉积明显[14]，沉积物向外海的输送能力较弱，湿地呈扇形分布向海淤进（图4a）。陆地面积快速增长为植被生长提供了栖息地，且低流速的淡水补给为植被提供了良好的生长环境，植被覆盖面积随之增加，植被类型以碱蓬等抗盐性植被为主，在河道周边多覆盖喜淡水植物芦苇（图5a）。

（2）慢速淤进阶段（1983—1995 年）：湿地面积增长相对缓慢，平均增长速率约为31.61 km2/a。在此期间，滩涂湿地与草甸沼泽呈现出同步变化的趋势。湿地面积与黄河入海累积输沙量和径流量亦呈现显著正相关，其线性关系的决定系数R2分别为0.68 和0.72（图3）。1983—1986 年期间，黄河清水沟河道逐渐归一成槽，呈单河道入海，径流泥沙向海输运效率明显提高，湿地随河道延伸而不断向东淤进，河流并汊后淡水的补给主要集中于河道两侧。在远离河道处，淡水补给缺少不利于植被发育，植被生长速度较慢，植被覆盖度分布呈现出由河道两侧向南北两岸递减，南北两岸出现裸地或无植被覆盖区（图4b）。1987 年湿地北侧由于孤东海堤的阻隔而不再被海水周期性覆盖，湿地面积下降，湿地呈现“楔”状（图4c），河道北岸植被覆盖度较低（图4d），淡水补给主要在河道附近，芦苇沿河道两岸分布，柽柳、碱蓬等抗盐类植物生长在芦苇外侧（图5b）。

图3 黄河口湿地面积随黄河累积入海输沙量（a）及累积径流量（b）的阶段性变化Qsa 为累积输沙量，Qa 为累积径流量。Fig.3 Stepwise variations of wetland area in the active Yellow River delta with cumulative sediment load （a） and cumulative water discharge （b） from the Yellow RiverQsa is cumulative sediment load, Qa is cumulative water discharge.

图4 不同时期河口湿地植被覆盖度空间分布Fig.4 Spatial distribution of wetland fractional vegetation cover at different stages

图5 不同时期河口湿地植被覆盖类型空间分布Fig.5 Spatial distribution of wetland vegetation types at different stages

（3）快速蚀退阶段（1996—2001 年）：1996 年5 月黄河入海流路改道清8 汊河，形成现行河口和废弃河口并存的格局。区域湿地面积出现逐年减小趋势，平均蚀退速率约为27.4 km2/a。该阶段湿地面积与黄河入海累积输沙量和径流量之间呈现显著的负相关关系，而同时段的清8 汊叶瓣湿地面积却与入海累积水沙无显著相关关系（图6a、6b）。上述反差表明，在该阶段河口湿地面积变化受控于南部废弃清水沟河口湿地的强烈蚀退，其蚀退速率远高于北部清8 河口新生湿地的发育速率。究其原因，在该时段黄河入海水沙量快速减少，下游多次出现断流情况（1997 年黄河断流时间高达226 天），导致清8 河口湿地发育缓慢，难以弥补海洋动力的侵蚀破坏；同时，废弃清水沟河口水沙供应断绝，湿地面积减小（图6c、6d），河口湿地受到强烈的海洋动力侵蚀，成为该阶段河口湿地时空变化的主导因素。

（4）稳定发育阶段（2004—2020 年）：研究区湿地面积处于相对稳定增加的状态，但湿地面积与黄河累积入海水沙量无显著相关关系，R2仅为0.04 和0.02（图3），而同期清8 河口湿地面积变化与黄河累积入海水沙量则呈现显著正相关，这表明南侧废弃清水沟河口区域的侵蚀速率与现行河口新生湿地发育的速率基本相当。在此阶段，伴随着南部废弃清水沟河口的持续蚀退，北部清8 河口向北淤进且堆积速率加快。其根本原因是2002 年以来黄河实施调水调沙，入海水沙以脉冲方式集中在6—7 月期间入海，期间黄河入海径流量增加，下游河道冲刷导致大量粗粒级泥沙被输送入海并迅速在河口处沉降，促进了清8 河口湿地面积增加[9]。2007 年河口向北出汊（图4g），在北部生成新的沙嘴，而原东向沙嘴发生侵蚀向陆后退，湿地冲淤共存，导致清8 河口湿地面积呈现波动性变化。此外，在此阶段黄河口南部区域潮滩湿地及小部分草甸沼泽湿地被人为开发，盐田、养殖池、库塘等人工湿地明显增加，使自然湿地面积减小，也可能对湿地发育造成一定影响。

3.2 清8 汊湿地时空变化及影响因素

现行清8 河口区自1996 年黄河改道后持续向海淤进，并伴随着河口湿地的快速扩展。在清8 河口区，无论是湿地面积还是湿地植被面积均与黄河累积入海输沙量及径流量呈现显著正相关性（图6—7），尤其是2002 年黄河调水调沙以来，入海粗颗粒泥沙量增加，河口造陆速率加快，为植被生长提供了栖息地，而入海径流增加则为植被生长提供了充足的淡水。

图6 湿地面积与黄河累积入海输沙量及径流量相关性分析Fig.6 Stepwise variation of wetland area along with the cumulative sediment load and cumulative water discharge from the Yellow River

图7 清8 汊湿地植被面积与黄河入海输沙量（a）及径流量（b）相关性分析Fig.7 Stepwise variation of vegetation area in the active river mouth （Q8） along with the cumulative sediment load （a） and cumulative water discharge （b） from the Yellow River

清8 河口湿地植被覆盖面积总体上呈现出阶段性变化特征（图8）：第一阶段（1996—2001 年），清8 河口发育初期，湿地状态不稳定，加之该阶段黄河断流频繁，淡水供应不足，植被覆盖面积相对较小；第二阶段（2004—2010 年），黄河实施调水调沙，入海径流增加同时入海泥沙颗粒粗化，促进了河口湿地的快速发育；第三阶段（2011—2020 年），湿地植被覆盖面积明显增加，高植被覆盖区和中植被覆盖区面积增加尤为突出。清水沟流路自2008 年起在每年6—7 月的调水调沙期进行河口三角洲湿地生态补水调度，通过洪水满溢和侧渗增加湿地水面面积，2010年后每年累积补水量增加至2000万m3以上[15]，淡水的补充为河口湿地植被提供了良好的生长环境。

图8 清8 汊叶瓣区各等级植被覆盖区面积阶段变化Fig.8 Temporal distribution of Q8 wetland fractional vegetation cover at different stages

1996—2020 年间7 期影像的解译数据反映了河口湿地植被覆盖的空间变化（图9），总体来看，低植被覆盖度区与高植被覆盖度区变化差异最为明显。1996 年低植被覆盖区域较为分散，分布在靠改道点处以及潮滩外侧，中高、高植被覆盖度区集中在改道点处；2001—2008 年期间低植被覆盖区域主要沿河道向海延伸至河口处，河道两侧有少部分中、中高植被度覆盖区；至2011 年，低河道两侧中高、高植被覆盖度区域增加明显，同时在潮滩外侧区域发育；2016—2020 年中高、高植被覆盖度区“包围”式覆盖于低植被覆盖度区外侧，且在沙嘴处有高植被覆盖度区域显著发育。

图9 不同时期清8 汊湿地植被覆盖度空间分布Fig.9 Spatial distribution of Q8 wetland Fractional Vegetation Cover during different periods

从植被覆盖类型来看（图10），2000 年低植被覆盖区域主要覆盖有碱蓬及少量柽柳，在改道处有芦苇分布；2004 年河口处碱蓬覆盖面积明显增大，北侧出现少量互花米草；2011 年河道两岸有芦苇分布，且叶瓣东岸覆盖有大量互花米草；2020 年河道两岸芦苇发育较好，柽柳与碱蓬混生，被互花米草“包围”其中。

总体而言，现行河口湿地的植被覆盖度和植被覆盖类型的空间分布均呈现明显的阶段性变化，其主要原因有两个方面：① 黄河淡水补给增加为植被生长提供了良好的繁殖环境，尤其是2008 年起进行的清水沟流路生态补水扩大了淡水的覆盖面积，河道两侧湿地为喜淡水的植被（例如芦苇）提供了良好的栖息环境；②互花米草的大量繁殖增加了植被覆盖面积（图10），尤其是2016 年以来，中高、高植被覆盖度区域面积占比明显增加，其中互花米草扩展速度最快。互花米草是一种光合效率高、种群密度大的盐生植物，生长繁殖迅速[16]，且由于是外来物种，在河口沙嘴及潮间带处的爆发增加了植被的覆盖度和覆盖面积。

图10 不同时期清8 汊湿地植被覆盖类型空间分布Fig.10 Spatial distribution of Q8 wetland vegetation types during different periods

3.3 废弃清水沟湿地时空变化及影响因素

废弃清水沟叶瓣位于研究区东南侧，1997 年南部清水沟河道被人工堤坝完全切断，与主干道断开。本部分研究中选取的时间段为1989—2020 年，跨越清水沟河道废弃前、中、后三个阶段。

该区域湿地面积经历了先上升后下降再趋于稳定的阶段性变化：第1 阶段为20 世纪80 年代到1996 年河道废弃前，湿地面积随着黄河累积入海泥沙量的增加而不断向海淤进。第2 阶段为1996 年清水沟河道废弃后，黄河入海水沙供应断绝，湿地持续萎缩，海洋动力导致的侵蚀成为影响湿地的主控因素。第3 阶段为2006 年后，废弃清水沟叶瓣湿地面积趋于稳定，平均面积为108.58 km2，湿地逐步趋于相对平衡状态。

从植被覆盖面积构成来看（图11），各年份湿地植被覆盖面积构成差别较大，但总体来看，低植被覆盖度区域在湿地中占比最大，中高、高植被覆盖度区域面积占比较小。1989 年，中高、高植被覆盖度区在黄河河道两侧零散分布；1996 年，随着河道继续向东南方向延伸，中高、高植被覆盖度区域向河道两侧靠拢，集中于河道两侧呈长条状分布；2001 年，在河道顶端植被茂盛，中高、高植被覆盖度区域在此处聚集，在河道两侧集中覆盖。2020 年，中高、高植被覆盖区覆盖范围趋于稳定，主要分布于废弃河道顶端及两侧（图12）。

图11 废弃清水沟叶瓣区湿地植被覆盖度时序变化Fig.11 Temporal distribution of Qingshuigou wetland Fractional Vegetation Cover

图12 不同时期清水沟湿地植被覆盖度空间分布Fig.12 Spatial distribution of Qingshuigou wetland Fractional Vegetation Cover during different periods

从植被覆盖类型来看（图13），1995 年南部河道尚未废弃，植被自改道点O 向海呈扇形延伸，由河道向两侧依次分布芦苇-柽柳-碱蓬，河道两侧芦苇分布茂盛，2000 年以后芦苇分布区呈长条状集中在废弃河道两侧，而柽柳和碱蓬在靠近海的外侧。相较于北侧清8 河口湿地植被，南侧废弃河口湿地仅在2011 年以后零星分布少量互花米草，以芦苇为代表的喜淡水、高植被覆盖度的湿生植被面积变化较小，而以柽柳、碱蓬为代表的植被覆盖度较低的盐生植被面积变化较大，受到海洋动力环境影响，面积持续下降。

图13 不同时期废弃清水沟湿地植被覆盖类型空间分布Fig.13 Spatial distribution of abandoned Qingshuigou wetland vegetation types during different periods

综上所述，在1996 年之前，废弃清水沟湿地时空变化主要受控于黄河入海水沙以及河口的向海淤积延伸；而在1996 年之后，由于入海水沙供应断绝，海洋动力引起的侵蚀成为湿地时空变化的主控因素。同时，由于潮汐不对称而引起的废弃河道再充填以及海水入侵加强，改变了废弃河道区域湿地植被的生境，由废弃前的淡水环境逐渐转变为废弃后的海水环境，喜淡水的芦苇植被受到逆境胁迫而导致面积减小。

4 结论

（1）基于1976—2020 年的长时间序列卫星遥感影像研究表明，黄河口湿地发育演变呈现显著的阶段性和空间差异性变化，经历了快速增加、稳定增长、快速蚀退、相对稳定四个阶段：① 1976—1981 年，河水分汊漫流，湿地以“扇形”扩张，面积快速增加；② 1983—1995 年，河道逐渐归一成槽，沙嘴先向东后向东南延伸，湿地形态沿河道方向呈现典型的西北宽东南窄的三角“楔形”，不断向东南方向稳定淤进；③ 1996—2001 年，黄河改道至清8 汊入海，清8 汊湿地开始发育，废弃清水沟湿地出现蚀退，北部增长速度不及南部蚀退速度，湿地面积不断减小；④ 2004—2020 年，北部清8 汊叶瓣湿地向海淤进，南部废弃清水沟叶瓣原突出的沙嘴逐渐被侵蚀，面积变化较小，空间上呈现清8 汊叶瓣湿地与废弃清水沟叶瓣湿地“双扇形”的形态。

（2）河口湿地植被的空间分布表现出集群特征：以中高、高植被覆盖度类型为主的植被主要分布在河道两侧淡水供应充足的区域，典型植被为喜淡水的湿生植被芦苇；低植被覆盖度类型植被主要分布于潮滩湿地处，距离河道越远植被覆盖度越低，受海洋动力改造作用较大，面积具有明显的年际变化，典型植被为盐生植被碱蓬。

（3）黄河口湿地时空变化与黄河来水来沙、海洋侵蚀、人类活动等因素密切相关，但影响北部清8 河口湿地与废弃清水沟湿地变化的主控因素截然不同：① 清8 河口湿地时空变化主要受控于黄河入海径流量和输沙量，湿地面积和植被覆盖与累积入海水沙量呈显著正相关，植被覆盖度与人类活动高度相关，如黄河人为改道、调水调沙、生态补水等；② 废弃清水沟叶瓣湿地的时空变化主要受控于海洋动力的侵蚀改造作用，同时潮汐不对称作用下废弃河道再充填以及海水入侵在一定程度上改变了河道两侧湿地的植被生境，致使植被面积逐渐减小，植被群落结构发生演替。