近30 年红河三角洲岸线演变

杨宏达，乔璐璐，刘雪，Le DucCuong，张蕾，李珏，印萍

1. 中国海洋大学海洋地球科学学院，青岛 266100

2. 海底科学与探测技术教育部重点实验室，青岛 266100

3. 越南科学院海洋环境与资源研究所，海防 180000

4. 中国地质调查局青岛海洋地质研究所，青岛 266100

红河三角洲位于越南北部（图1），面积约为17 000 km2[1]，受波浪、潮流、河流共同控制[2]。红河一共有8 条支流，年均径流量达120 km3[3]，其中大型的河口主要有白藤（Bach Dang）河口、太平（Thai Binh）河口、巴拉特（Balat）河口、宁科（Ninh Co）河口、天（Day）河口，其中巴拉特河口的径流量约占总量的65%[4]。红河年均入海泥沙量约8.2×107m3[5]，其中90%在雨季输送[6]，主要沉积在北部湾西南近海区域[7]。

红河三角洲所位于的北部湾海域以全日潮为主，平均潮差2.5～3.5 m；波浪在旱季为东—东北方向，在雨季为东—东南方向，平均和最大波高分别为0.7～1.3 m 和3.5～4.5 m[8]。根据前人研究，北部湾1990—2010 年海岸线的总长度在持续减小，20 年间减小了118.50 km，其中后10 年减小幅度明显大于前10 年[9]，且岸线趋向于平直化[10]。北部湾岸线的变迁主要和人类活动有关[11-12]，海岸线被大量开发为港口、码头、养殖区等[13]。砂质、淤泥和基岩岸线的占比逐渐减小[14]，到1998 年人工堤坝已占北部湾岸线总长度的39.46%[15]。除人类活动影响外，河口岸线变化也较大[16-17]。

关于红河三角洲岸线，前人主要分析了岸线变化特征以及河口沙坝、上游大坝等对其的影响：红河三角洲岸线在1999 年之后的10 年总体上向陆迁移[5]，岸线变化频繁剧烈，受多种因素共同影响。关于河口沙坝的影响，Do Minh Duc 等发现，河流输送的泥沙主要聚集在河口，使海岸线快速淤积（可达100 m/a），而其他岸段因缺少泥沙供应，导致侵蚀[4,18-19]。关 于 上 游 大 坝 的 影 响，Tran Duc Thanh 等发现，和平大坝建立之后，大坝下游河水中泥沙含量减少，宁科（Ninh Co）河口附近岸线的年平均后退速率从1965—1990 年的0.41 m/a 提高到了1990—1998 年的14.5 m/a[20]，但对主要河口巴拉特的岸线变化速率在大坝建立前后的变化缺少研究。另外，入海泥沙量变化导致的河口演变[18]、波浪斜向入射[21]、潮流和波浪流[22]都会对海岸侵蚀产生重要影响。近年来岸线趋于平直，各种人工建筑物也会影响波浪和潮汐的动力作用，从而影响岸线侵蚀[13,23]。虽然红树林能很好地应对海岸侵蚀[24]，但由于人类的砍伐，红树林也受到了很大程度的破坏。加剧的岸线侵蚀将导致建筑物不稳定等各种灾害发生，三角洲生态服务功能也会逐渐减弱[25-27]。

前人对红河三角洲岸线演变进行了定性描述，但缺少定量化分析，而且研究年份仅截止于2001年，对人类活动频繁的近20 年的海岸变化相关报道较少。近年来，由于全球变暖，海平面上升，台风发生频率增加[28]，海洋动力发生了重要变化，对岸线变化产生了重要影响。在岸线变化机制上，前人的研究分析了海洋动力和河流入海泥沙对岸线变化的影响，但是对于北部湾海域每年夏季高频次的台风作用分析较少。前人已在山东威海海滩观测发现，台风梅花（1109 号）过境后，天鹅湖沙滩后退了3.6 m[29]，热带气旋对岸线变化产生了剧烈影响。因此，本文基于1987—2015 年Landsat 遥感数据，收集1935—1985 年历史数据、红河多年径流量与输沙量，以及该海域多年热带气旋数量等数据，对红河三角洲变化最频繁剧烈的巴拉特河口的岸线演变开展研究，从而为岸线防护以及海岸带经济发展提供科学支撑。

1 数据与方法

本文采用的基础数据来自Landsat 系列卫星的影像数据，从美国地质勘探局官网（网址：https://earthexplorer.usgs.gov/）下载，影像分辨率为30 m。为保证影像精度及多年研究中的可对比性，选取云量小于10%、高潮时的遥感影像进行对比分析。基于ENVI 软件，进行辐射校正、拼接、大气校正、研究区裁剪、几何精校正等预处理[30]。预处理之后，根据“我国近海海洋综合调查与评价项目（908 专项）”海岸带遥感调查和前人岸线解译经验[31-32]，在MapInfo 软件中提取大潮高潮线作为岸线，共提取了1987、1990、1995、2000、2005、2010、2015近30 年共7 张红河三角洲的岸线，结果如图1 所示。

图1 1987—2015 年红河三角洲海岸线Fig.1 Red River Delta coastline changes from 1987 to 2015

2 多年岸线变化特征

2.1 1930—1985 年巴拉特地区岸线变化

根据收集的资料[8]，前人已对1930—1985 年巴拉特地区岸线变化进行了研究（图2），可以看出在巴拉特河口地区前进与后退剧烈，在靠近河口部位发生明显前进，而远离河口部位发生明显后退。

图2 1930—1985 年巴拉特河口岸线变化[8]Fig.2 Changes of the Balat estuary shoreline in 1930-1985[8]

2.2 1987—1990 年红河三角洲岸线变化

1987—1990 年，红河三角洲南部岸线有多处明显前进，而北部出现多处后退，中部变化不大（图1），主要是巴拉特地区变化明显。通过计算得到如图3 所示的巴拉特地区主岸线沿岸平均年变化率，可以看出巴拉特地区沿岸大部分岸线以后退为 主。河口位置发生明显后退，从遥感图像上对比发现，该岸线后退可能与入海泥沙量减少和沙坝侵蚀有关。定量对比发现，1987—1990 年巴拉特地区年平均前进速率为39.19 m/a，年平均后退速率为30.68 m/a。

图3 1987—1990 巴拉特河口岸线变化Fig.3 Changes of the Balat estuary from 1987 to 1990

2.3 1990—1995 年红河三角洲岸线变化

1990—1995 年，红河三角洲除河口岸线以前进为主外，河口的南部和北部岸线均出现了多处明显后退（图1）。巴拉特地区变化明显，通过计算得到如图4 所示的巴拉特地区主岸线沿岸平均年变化率，可以看出巴拉特地区沿岸大部分岸线以前进为主，而河口位置发生明显后退，从遥感图像上对比发现，在这一时期，由于波浪等导致沙坝侵蚀或者由于海平面上升，主岸线与沙坝之间的人工养殖池长时间被水淹没，主岸线发生明显后退。1990—1995 年巴拉特地区年平均前进速率为28.26 m/a，年平均后退速率为87.25 m/a。

2.4 1995—2000 年红河三角洲岸线变化

1995—2000 年，红河三角洲岸线变化不大，南部出现了几处明显前进，北部也出现几处明显前进，中部变化不大（图1）。对于巴拉特地区，岸线的后退与前进交替出现，在河口位置发生了一处明显前进（图5），从遥感图像上对比发现，该岸线前进与人工养殖池有关。1995—2000 年巴拉特地区年平均前进速率为23.63 m/a，年平均后退速率为11.57 m/a。

2.5 2000—2005 年红河三角洲岸线变化

2000—2005 年，红河三角洲南部岸线发生一处明显前进，而北部出现几处后退，中部变化不大（图1）。巴拉特地区岸线的后退与前进交替出现，变化不剧烈（图6）。2000—2005 年巴拉特地区年平均前进速率为5.93 m/a，年平均后退速率为5.06 m/a。

图4 1990—1995 巴拉特河口岸线变化Fig.4 Changes of the Balat estuary from 1990 to 1995

图5 1995—2000 巴拉特河口岸线变化Fig.5 Changes of the Balat estuary from 1995 to 2000

2.6 2005—2010 年红河三角洲岸线变化

2005—2010 年，红河三角洲南部岸线有多处明显前进，而北部出现多处明显后退，中部出现多处前进（图1）。巴拉特地区沿岸岸线以后退为主（图7），从遥感图像上（图8）对比发现，该岸线后退可能与2009 年 宣 光（Tuyen Quang）大 坝 和2010 年 桑 拉（Sonla）大坝建成导致入海泥沙量减少有关，因此沙坝侵蚀，人工养殖池被水淹没，主岸线发生明显后退。2005—2010 年巴拉特地区年平均前进速率为20.87 m/a，年平均后退速率为56.19 m/a。

2.7 2010—2015 年红河三角洲岸线变化

2010—2015 年，红河三角洲南部岸线以明显前进为主，北部也以明显前进为主，但出现几处后退，中部以前进为主（图1）。巴拉特地区沿岸岸线以前进为主，后退仅限于河口东北部（图9），从遥感图像上对比发现，该岸线后退可能与入海泥沙量减少，沙坝侵蚀有关。2010—2015 年巴拉特地区年平均前进速率为22.13 m/a，年平均后退速率为89.72 m/a。

图6 2000—2005 巴拉特河口岸线变化Fig.6 Changes of the Balat estuary from 2000 to 2005

图7 2005—2010 巴拉特河口岸线变化Fig.7 Changes of the Balat estuary from 2005 to 2010

图8 2005 年和2010 年巴拉特河口遥感图像Fig.8 Remote sensing images of the Balat estuary in 2005 and 2010

2.8 1987—2015 年红河三角洲岸线总体变化

图9 2010—2015 巴拉特河口岸线变化Fig.9 Changes of the Balat estuary from 2010 to 2015

观察多年岸线变化图（图1），可以看出沿岸地区变化不一致，部分岸段前进与后退交替变化特征明显，在巴拉特河口地区淤进与蚀退频繁，变化剧烈，河口北部岸线后退明显，主要受沙坝侵蚀影响。计算得到1987—2015 年巴拉特地区主岸线年平均前进速率为3.94 m/a，年平均后退速率为29.43 m/a。前人计算得到长江口河口岸线1974 年到2015 年的平均变化速率为34 m/a[20]，与其对比，发现巴拉特地区主岸线年平均变化速率大小与其类似（图10）。

图10 1987—2015 年巴拉特地区主岸线沿岸平均年变化率Fig.10 Average annual change rate of the main shoreline in Balatregion from 1987 to 2015

3 岸线多年变化原因分析

通过调查资料，发现红河三角洲海岸线演变受热带气旋以及河流上游水坝等多种因素共同影响。

3.1 河口沙坝的影响

海岸线的向海推进与沙坝的形成和扩大密切相关。大量入海沉积物在河口前形成沙坝，沙坝的快速增长降低了波浪的破坏作用，为海岸线和沙坝之间的河道中沉积物沉降提供了良好的环境[9]。河口沙坝在红河三角洲的巴拉特地区作用非常明显。当沙坝被侵蚀时，沙坝与主岸线之间的人工养殖池也被水淹没，因此，主岸线明显后退，如1990—1995年和2005—2010 年。

3.2 河流上游建坝的影响

图11 1960—2010 年4 个站位径流量和输沙量变化[33]Fig.11 Variation in runoff and sediment budget at four stations from 1960 to 2010[33]

红河入海径流量和输沙量受堤坝建设影响明显（图11），进而会对岸线前进与后退速率产生影响[20]。前人在1972 年第一个大坝塔婆建立之后，在山西站（Son Tay，图11e）观测到的年输沙量下降[33]。根据前人对1965—1990 年宁科（Ninh Co）地区岸线变化的研究结果，可以发现在和平大坝建立之后，宁科河口岸线的年平均后退速率从1965—1990 年的0.41 m/a 提高到了1990—1998 年的14.5 m/a[21]。而前人对1930—1985 年巴拉特地区岸线变化的研究结果表明，在塔婆大坝建立之后巴拉特河口前进速率相应发生了降低[3]。本文对比发现，在2009 年宣光（Tuyen Quang）大坝和2010 年桑拉（Sonla）大坝建成之后，沉积物供应减少，2005—2015 年巴拉特河口岸线后退速率有所提高。

3.3 气候变化影响

厄尔尼诺现象可造成世界气候的变化，使局部地区降雨量过多，极大地影响岸线的变化速率。调查资料发现，红河上游的大东勇站、中游的元江站、下游的蛮耗站的年径流序列在1993 年均由枯水期突变为丰水期[34]，这可能与1991—1994 年发生的厄尔尼诺带来的强降雨有关[4]，沙坝与主岸线之间的人工养殖池被水淹没，1990—1995 年巴拉特河口岸线的后退速率明显增长。同时，2006—2007 年、2009—2010 年均发生了厄尔尼诺现象，对应的2005—2010 年巴拉特河口岸线的后退速率也有明显增长。

海平面上升对海岸侵蚀速率的影响也很大。根据前人研究结果，1965—1995 年期间，海平面上升对红河三角洲海岸侵蚀速率增加的作用约占34%，而1995—2005 年为12%。

3.4 热带气旋的影响

热带气旋是越南北部地区典型的天气事件，台风数量与岸线年平均后退速率有很好相关性（图12），台风数量越多，平均后退速率越大。1990—1995年经过红河三角洲台风数量增多，年平均后退速率有明显提高；1995—2005 年经过红河三角洲台风数量减少，年平均后退速率对应有明显降低；2005—2015 年经过红河三角洲台风数量增多，年平均后退速率对应有明显提高。海后（Hai Hau）潮汐站曾测得当波高为4.25 m，持续时间2.4 h 的情况下，岸线后退了7.1 m[8]。

3.5 海洋动力的影响

海流和波浪是影响海岸侵蚀的重要因素[5]。前者在红河三角洲海域呈现出明显的季节性变化：冬季流向南、夏季流向北[7]，并可能通过诱发增水等影响岸线变化。后者在近岸破碎并诱生沿岸流，进而引起沿岸输沙作用，造成海岸侵蚀[8]。

无论旱季或者雨季，巴拉特河口南部海域平均波高始终高于河口北部（图13a，图13b），河口南部的海洋动力始终强于北部。从沉积物输运方向（图14）也可以看出：在巴拉特、拉赫（Lach）和天河口，5～10 m 等深线海域，沉积物主要向东南输运，5 m 等深线以浅的区域，沉积物沿岸向西南输运；在特雷（Tra Ly）河口，沉积物向东北输运；在太平河口，沉积物向东输运[35]。因此，巴拉特河口北部和南部岸线呈现不一样的变化特征；北部岸线后退与前进交替出现，总体以前进为主，可能与复杂的岸线和动力条件及河流入海物质变化有关；南部岸线持续后退，与波浪导致的沿岸南向输沙密切相关。

图12 1987—2015 年红河三角洲年均后退速率与过境台风数量之间的关系Fig.12 The relationship between the annual averaged retreat rate of the Red River Delta and the numbers of typhoons passed through the delta from 1987 to 2015

图13 冬(旱)季(a)与夏(雨)季(b)平均波高等值线[36]Fig.13 Average wave height in winter/summer（dry/rain）seasons[36]

图14 红河三角洲近岸沉积物输运路径[35]Fig.14 Sediment transport pathways along the coast of the Red River Delta[35]

4 结论

（1）在时间变化上，红河三角洲岸线表现为前进与后退交替变化特征。南部岸线变化趋势较为稳定，除1990—1995 年、2010—2015 年以前进为主外，1935—2015 年其他年份均以后退为主，与波浪导致的沿岸南向输沙密切相关。红河三角洲北部岸线由于岸线复杂且为河流入海物质控制区，表现为后退与前进交替出现，总体以前进为主。

（2）红河主要支流之一巴拉特河口地区淤进与蚀退变化剧烈，1987—2015 年巴拉特河口主岸线年平均前进速率为3.94 m/a，年平均后退速率为29.43 m/a。该河口岸线的频繁变化主要与上游建坝造成入海沙量变化、台风造成的侵蚀有关。