闽江河口潮滩季节性冲淤变化格局及其控制机制

王爱军，李海琪，叶翔，梁灏深，张望泽,5，吴水兰,6，冉畅，陶舒琴，刘子同,7，于谦

1.自然资源部第三海洋研究所海洋与海岸地质研究室，厦门 361005

2.福建省海洋物理与地质过程重点实验室，厦门 361005

3.自然资源部海峡西岸海岛海岸带生态系统野外科学观测研究站，厦门 361005

4.上海海洋大学海洋科学学院，上海 200090

5.福州大学先进制造学院，泉州 362251

6.南京信息工程大学海洋科学学院，南京 210044

7.汕头大学海洋科学研究院，汕头 515063

8.南京大学海岸与海岛开发教育部重点实验室，南京 210023

潮滩是在潮汐作用显著、细颗粒沉积物供应丰富条件下形成的一种地貌类型，位于陆地、海洋、大气及人类活动相互作用的敏感地带。潮滩广泛分布在除基岩海岸和高纬度（70°～73°以上）海岸外的沿海地区，并且随着潮差、波浪、物质供应、植被类型等不同而表现出显著的空间差异[1-4]。目前全球 潮 滩 面 积 约 为127 921 km2，主 要 分 布 在 亚 太 地区，如印尼、中国、印度等国家[5]。我国的潮滩规模大且分布广泛，其中长江口-杭州湾及以北沿海地区主要为平原型潮滩，而长江口-杭州湾以南的浙、闽、粤、桂等沿海地区主要为港湾型潮滩[6]。从20 世纪初开始，各国学者对潮滩沉积和地貌特征、潮滩沉积物输运和堆积的过程、盐沼植被对潮滩沉积过程的影响、极端事件对潮滩地貌演化影响等方面进行了系统研究和总结[7-16]。进入21 世纪，随着现场观测技术和数值模拟的不断改进，我国科研人员开始关注潮滩沉积过程及其控制机制、植被-泥沙-地貌耦合的生物地貌过程、极端事件影响下的潮滩沉积物输运与地貌演化过程、潮滩沉积层序成因机制及潮滩演化趋势模拟等方面，取得了丰硕的成果[17-33]。

潮滩作为重要的海岸带生态系统，不仅在蓝色碳汇方面发挥着重要作用[34-35]，而且在海岸防护方面起到至关重要的作用[36]。然而，随着海平面不断上升、地面沉降、流域及海岸带人类活动的不断加剧，潮滩面积逐渐减少[5]，其中人类活动是主要驱动因素[37]。潮滩面积的减少不仅导致蓝色碳汇能力减弱[5,35]，而且引起诸多地区海岸遭受不同程度侵蚀[38]。高强度的人类活动显著影响了潮滩局部沉积物输运过程及地貌演化过程[20,24]，而流域输沙量的锐减则引起了潮滩冲淤格局的改变[39]。随着流域及河口地区人类活动强度的不断增大，河口系统状态正在发生转换[40]，直接影响着河口潮滩未来演化趋势及其功能发挥。

本文选择我国东南沿海典型的山溪性中小型河流河口潮滩作为研究区域，通过沉积物采样、定点冲淤观测和沉积动力学观测等手段，初步分析了河口潮滩的季节性冲淤变化过程及其控制机制，以期为科学评估潮滩的海岸防护和潮滩生态系统保护与修复能力提供科学依据。

1 研究区概况

闽江发源于福建与江西省交界的武夷山东麓及仙霞岭南麓，总体呈NW-SE 向横贯福建中北部，经福州汇入东海南部，全长2 872 km，其中干流全长近577 km，流域面积为60 992 km2，多年平均径流量为605.5×108m3，多年平均输沙量为750×104t，并且表现出明显的洪、枯季变化，即洪季（3—8 月）径流量约占全年的76%，输沙量约占全年的92%[41]。近年来随着流域水库建设的不断增加，平均输沙量已经下降为原来的三分之一[42]。

受地转偏向力及台湾海峡地形效应作用，闽江河口区潮汐作用强，实测最大潮差达到7.04 m，平均潮差4.46 m（梅花站），属正规半日潮，河口内潮流基本呈往复流，其流向与河槽线一致，河口外潮流具有一定的旋转流特征[41]。闽江口外海面开阔，夏季以西南风为主，冬季以东北风为主，波浪作用强，波浪的主要形式为风浪及涌浪同时存在的混合浪，多年平均波高为1.1 m，历史最大波高为6.5 m[41]。此外，闽江口海域是我国东南沿海地区遭受台风影响最为严重的区域之一，平均每年影响的台风次数为5.7 个，最大台风增水值可达2.52 m[43]。

闽江河口地区的潮滩主要分布在川石岛南部、琅岐岛东部以及梅花水道南侧海岸（图1a），其中琅岐岛潮滩宽度为600～1 100 m（图1b），滩面平缓。现场调查结果显示（图1c），琅岐岛潮滩由岸向海主要分布有互花米草、海三棱藨草、光滩等，其中互花米草盐沼滩面平均坡度约为0.36°，海三棱藨草盐沼滩面平均坡度约为0.14°，潮间带中下部的光滩平均坡度约为0.08°，河道边缘地区的滩面坡度约为1.84°。

图1 研究区位置（a）、沉积物采样站位（b）、剖面形态及观测站位图（c）Fig.1 Sketch map of study area (a), sediment sampling sites (b), morphological character of profile and observation sites of tidal flat of Langqi Island (c)

2 研究方法

2.1 野外采样与观测

在琅岐岛东部潮滩选择一条剖面利用RTK 对剖面进行了滩面高程测量，并在剖面上由岸向海设置了10 个采样站位（图1b），分别于2022 年8 月和2023 年2 月采集了表层1 cm 厚度的沉积物。为了掌握潮间带不同空间位置的冲淤变化特征，在互花米草盐沼边缘（L03 站）、海三棱藨草盐沼中部（L05 站）及光滩外缘（L09 站）设置了3 个站位，安装了冲淤观测架，在观测架上设置12 个点，定期测量每个点距离滩面的高度，然后对各站位每次的测量数据进行平均，获得各站位的冲淤变化信息（图1c）。

为了解研究区水动力状况，在海三棱藨草内部靠近互花米草边缘设置一个沉积动力学观测站位（图1c，L04 站），利用Nortek 公司生产的声学多普勒海流计（ADV，VECTOR 6MHz）分别于2022 年1 月25 日至2 月1 日、7 月21—31 日开展近底部潮流观测，仪器探头距离滩面30 cm，观测点位于滩面以上15 cm 处；观测设置为脉冲模式进行采样，采样间隔为10 min，其中冬季采样频率为4 Hz，每次采集1 024 组数据，夏季采样频率为16 Hz，每次采集4 096组数据；冬季利用北京海洲赛维科技有限公司生产的波潮仪（TWaves）进行淹没水深观测，以脉冲模式进行采样，以10 min 为采样间隔，采样频率设置为8 Hz，每次采集时长为5 min，取平均值。

2.2 样品分析与数据处理

2.2.1 沉积物粒度分析

在实验室内将采集回来的常温保存样品充分混合，取约2 g 重的沉积物加入浓度为5%的H2O2溶液以去除有机质，再加入浓度为1 mol/L 的HCl溶液以去除碳酸盐，用蒸馏水清洗至溶液为中性，加入浓度为5‰的六偏磷酸钠溶液（(NaPO3)6）浸泡24 h，并利用超声波震荡1 min，使样品充分混和、分散后，用英国马尔文公司生产的 Mastersizer 2000 型激光粒度仪进行粒度分析，并用仪器自带软件导出中值粒径和各粒级组分数据，采用谢帕德分类体系对沉积物进行分类[44]。

2.2.2 悬沙浓度校对

利用在L04 站采集的表层沉积物在室内配置不同悬沙浓度（SSC）的水体与ADV 采集的声学信号进行相关分析，结果发现两者之间具有很好的相关性（图2）。因此，根据该回归方程可以将ADV 观测获得的信噪比（Signal-Noise Ratio，SNR）转化成时间序列的SSC。

图2 悬沙浓度与ADV 观测获得的信噪比的相关关系Fig.2 Relationship between SSC and SNR recorded by ADV

2.2.3 沉积动力学参数计算

利用仪器自带软件将观测数据导出，然后再利用该软件的数据转换功能，将原始数据转换成数据文件以便计算水平方向上的平均流速和波浪参数。

潮间带地区作为极浅水环境，整个水层都属于边界层范围内，平均水平流速在垂向上的分布遵循对数分布[45]：

其中，U(z)是平均水平流速度，u*是摩阻速度；κ是卡门常数（κ=0.4)；z是距底床高度，z0是底部粗糙长度，τb是潮流引起的底部切应力，ρw是海水的密度。一般在边界层范围内，可以利用对数剖面方法计算潮流引起的底部切应力。

有效波高（Hs）采用如下公式计算[46]：

式中，Sη是水位的功率谱。将ADV 测得的压力转换到水位时，需要考虑衰减，根据线性波理论，衰减系数可以表示为[47]：

其中，k是波数，h是平均水深，z是压力传感器的深度（负值）。波浪引起的底部切应力（τw）可以用如下公式计算[45]：

其中，fw为波浪摩擦因子，Uw为波浪轨迹速度，可以表达为波高（Hs）、波周期（T）和水深（h）的函数[45]：

其中k为波数，可以表达为：k= 2π/L，L为波长，可以用图解法计算得到[48]。波浪摩擦因子可以表达为：

式中A为半轨迹距离，可以表达为：A=UwT/2π。

由流速和悬沙浓度数据计算近底层悬沙瞬时水平输运率f(t)及潮周期内悬沙净水平输运通量F(t)[45]：

其中U(z,t)和C(z,t) 分别表示距离海底高度z(m)处t时刻的瞬时水平流速（ua，m/s）和悬沙浓度(SSC, mg/L)，Δt表示实测时刻代表时间长度。

研究区表层沉积物冬季以粗颗粒沉积物为主，夏季以细颗粒沉积物为主，因此，在计算沉积物临界侵蚀切应力时分别选择适合非粘性沉积物和粘性沉积物的计算公式。非粘性沉积物临界侵蚀切应力时选择如下计算公式[45]：

其 中，θ是Shields 参 数，θcr是 沉 积 物 临 界 起 动Shields 参 数，s=ρs/ρw，ρs是 沉 积 物 密 度（取 值 为2 650 kg/m3），ρw是海水密度（取值为1 025 kg/m3），d是 沉 积 物 粒 径， g 是 重 力 加 速 度（ 取 值 为9.81 m/s2）,D*是一个与沉积物粒径有关的无量纲参数，ν是水体的运动黏滞系数，取值为1.36×10-6m2/s。对于黏性沉积物，其临界侵蚀切应力计算用如下公式[49]：

其中，β是与颗粒粒径有关的系数，根据试验资料[49]，在中值粒径小且表层沉积物含水量较大时，取β=0.3，W为含水量，定义为样品中水的质量与沉积物质量之比，根据采集的现场表层底质样品的分析获得。

3 研究结果

3.1 沉积物粒度的季节分布特征

粒度分析结果显示（图3、图4），琅岐岛潮滩表层沉积物粒度组分以粉砂和砂为主，黏土含量相对较少，呈现出显著的时空差异。从季节方面看，冬季表层沉积物砂含量显著增大，黏土含量显著减小，沉积物中值粒径Φ 值较小；夏季则相反，砂含量显著减小，黏土含量和粉砂明显增大，沉积物中值粒径Φ 值较大（表1）。从空间上看，潮间带上部的互花米草盐沼内表层沉积物组分以粉砂为主，平均含量为60.6%，砂含量平均为21.0%，黏土含量平均为18.4%，中值粒径平均为5.60 Φ；在潮间带中部的海三棱藨草盐沼，沉积物组分以砂为主，平均含量为57.1%，粉砂含量平均为32.4%，黏土含量平均为10.5%，中值粒径平均为4.31 Φ；在潮间带下部的光滩地区，表层沉积物组分以粉砂为主，平均含量为58.7%，砂含量平均为23.0%，黏土含量平均为18.3%，中值粒径平均为5.56 Φ（图3、图4）。

表1 琅岐岛潮间带不同季节表层沉积物粒度组分及中值粒径Table 1 Seasonal variations of composition and median grain-size of surficial sediment

图3 琅岐岛潮滩夏季（左）和冬季（右）表层沉积物粒度组成Fig.3 Surficial sediment composition along the profile of tidal flat in summer (left) and winter (right)

图4 琅岐岛潮滩不同季节表层沉积物中值粒径Fig.4 Seasonal variations of median grain-size of surficial sediment along the profile of tidal flat

根据Shepard 三角分类图，琅岐岛潮间带表层沉积物类型以黏土质粉砂（YT）为主，并且呈现出显著的季节变化，其中夏季表层沉积物黏土质粉砂占80%，仅在潮间带中部出现砂质粉砂（ST）和砂-粉砂-黏土（STY）；冬季以砂（S）和砂质粉砂为主，其次为粉砂质砂（TS）（图5）。

图5 琅岐岛潮滩不同季节表层沉积物类型Fig.5 Seasonal variations of surficial sediment type along the profile of tidal flat based on Shepard classification

3.2 潮滩侵蚀-淤积的时空变化特征

潮间带不同空间位置的冲淤观测结果显示，琅岐岛潮间带滩面冲淤表现出显著的时空变化特征。互花米草盐沼边缘滩面表现出总体持续淤积的变化特征，仅在冬末至初春季节出现轻微侵蚀（0.55 cm），观测期间年淤积厚度为1.22 cm（图6a）。海三棱藨草中部滩面在初夏至初秋季节表现为持续淤积，淤积厚度为14.13 cm；初秋至冬末季节表现为持续侵蚀，侵蚀深度为20.13 cm；冬末至春末期间呈现出缓慢淤积，淤积厚度为2.54 cm；整个初夏至翌年春末表现为净侵蚀，侵蚀深度为3.42 cm（图6b）。光滩外缘滩面在初夏至秋末季节表现为持续侵蚀，侵蚀深度为5.90 cm；秋末至翌年春末季节表现为持续淤积，侵蚀深度为5.38 cm，整个初夏至翌年春末表现为净侵蚀，侵蚀深度为0.52 cm（图6c）。

图6 潮间带不同站位冲淤变化过程正值为淤积，负值为冲刷。Fig.6 Variations of erosion and accretion along the section in different stations Positive value means accretion, negative value means erosion.

3.3 潮间带水动力特征

琅岐岛潮间带中上部L04 站位淹没水深观测结果显示（图7a、8a），冬季和夏季观测期间大潮最大淹没水深分别为2.20 和1.86 m，小潮最大淹没水深分别为1.37 和0.90 m。潮流观测结果显示（图7b、8b），冬季小潮观测期间近底部潮流流速为0.008 ～0.127 m/s，大潮期间为0.004 ～0.187 m/s；夏季小潮观测期间近底部潮流流速为0.001 ～0.046 m/s，大潮期间为0.001 ～0.060 m/s；无论冬季还是夏季，潮周期内潮流流速最大值均出现在落潮后期。波浪计算结果显示（图7c、8c），冬季和夏季观测期间有效波高最大值分别为0.63 和0.20 m，平均值分别为0.26 和0.07 m，冬季波浪明显较夏季强。底部切应力计算结果显示（图7d、8d），冬季和夏季观测期间潮流引起的底部切应力最大值分别为0.551 和0.057 N/m2，平均值分别为0.067 和0.003 N/m2；波浪引起的底部切应力最大值分别为1.885 和0.293 N/m2，平均值分别为0.652 和0.097 N/m2，冬季潮流和波浪引起的底部切应力明显大于夏季。

图7 潮滩中上部光滩冬季近底部水动力时间序列变化过程a：淹没水深，b：近底部流速，c：有效波高，d：底部切应力，e：近底部悬沙浓度。Fig.7 Time-series of near-bottom hydrodynamic processes in winter on the mid-upper tidal flat a: Inundation height, b: near-bottom mean current speed, c: significant wave height, d: bottom shear stress induced by tidal current and wave, e: near-bottom suspended sediment concentration.

3.4 潮间带近底部悬沙浓度变化

悬沙浓度计算结果显示，琅岐岛潮间带中上部L04 站位近底部悬沙浓度呈现出显著的时间变化特征（图7e、8e）。在季节尺度上，冬季观测期间近底部悬沙浓度为17.4 ～844.5 mg/L，平均值为306.8 mg/L，而夏季观测期间近底部悬沙浓度为3.4 ～607.8 mg/L，平均值为138.5 mg/L，冬季悬沙浓度显著高于夏季。在大-小潮周期尺度上，小潮期间近底部悬沙浓度相对较低，大潮期间相对较高，但潮周期平均值相差不大，如冬季大潮和小潮观测期间近底部悬沙浓度平均值分别为296.5 和268.6 mg/L，而夏季大潮和小潮观测期间近底部悬沙浓度平均值分别为158.1 和123.7 mg/L。在潮周期尺度上，涨潮初期和落潮后期相对较高，高平潮期间悬沙浓度相对较低。

3.5 潮间带近底部悬沙输运

近底部悬沙东向水平输运率计算结果呈现出显著的时间变化特征（图9），冬季观测期间近底部悬沙在东向上水平输运率为-0.118 ～0.051 kg·m-2·s-1，平均值为0.011 kg·m-2·s-1；夏季观测期间近底部悬沙水平输运率为-0.006 ～0.009 kg·m-2·s-1，平均值为0.001 kg·m-2·s-1。

每个潮周期内的东-西方向上的净悬沙输运量计算结果显示（图9），琅岐岛潮间带中部潮周期内近底部悬沙水平净输运通量呈现出明显的时间变化特征。在大-小潮时间尺度上，小潮期间潮周期内悬沙净水平输运通量显著小于大潮，而净输运方向则呈现出季节性差异：小潮期间冬季呈现出向岸净输运和向海净输运交替出现，而夏季则整体呈现向海净输运的变化特征；中潮期间冬季总体表现为向海净输运，而夏季则整体呈现向陆净输运的变化特征；大潮期间冬季总体呈现向岸净输运，而夏季则呈现向岸净输运和向海净输运交替出现的变化特征。在季节尺度上，冬季潮周期内悬沙净水平输运通量为-124.30～243.24 kg/m，观测期间净输运通量为289.49 kg/m；夏季周期内悬沙净水平输运通量仅为-38.07 ～4.13 kg/m，观测期间净输运通量仅为-58.73 kg/m，即冬季潮周期内净水平输运通量显著大于夏季，并且冬季悬沙净向海输运，而夏季净向陆输运。

4 讨论

河口潮滩是由地质过程以及水动力过程、沉积物供应和生态过程的相互作用形成的，其地貌形态及沉积物组成受动力、物源及生物活动控制，具有高度的时空变化特征[4,50]。大量的河流入海泥沙进入河口后，在水动力的作用下发育了广泛的潮滩[6,51]，但由于我国河流泥沙入海主要集中在夏季，冬季入海泥沙通量很小[41]，因此，潮滩表层沉积物也表现出显著的季节差异[52-54]。闽江河流泥沙入海通量也主要集中在夏季[55]，本文断面表层沉积物分析结果（图3、图4）和冲淤观测结果（图6）也显示，夏季期间河口潮滩捕获了大量的河流入海沉积物，潮滩为淤积状态，潮滩表层沉积物以粉砂和黏土为主，平均粒径较小；冬季期间河流入海泥沙很少，河口潮滩很少能捕获到来自河流的泥沙，潮滩表层沉积物在水动力的作用下不断发生改造，最终细颗粒沉积物被搬运带走，潮滩发生侵蚀，表层沉积物以砂和粉砂为主。关于沉积物供应变化对闽江口潮滩沉积物组成的影响综合分析详见李海琪等[56]，本文不再赘述，但需要指出的是，随着闽江入海泥沙通量的显著减少[55]，河流输沙季节变化对琅岐岛潮滩季节性冲淤变化的影响将会减小。本文将主要从水动力方面阐述闽江河口潮滩冲淤的季节性变化机制。

影响潮滩沉积和地貌演化的动力过程主要包括潮汐、波浪、跨岸/沿岸流、风致环流、潮不对称等[50,57]，尤其是在极浅水环境下，潮流、波浪是影响潮滩沉积物分布的主要动力因素[27-28]。沉积物是否发生侵蚀和淤积，主要是看水动力引起的底部切应力是否会超过临界切应力，当水动力引起的底部切应力大于沉积物的临界侵蚀切应力，则沉积物发生再悬浮，滩面呈现侵蚀状态；当水动力引起的底部切应力小于泥沙临界沉降切应力，则泥沙发生沉降，滩面呈现淤积状态。根据本文观测期间采集的表层沉积物粒度分析结果，冬季表层沉积物砂含量达到了93.9%，中值粒径为2.85 Φ，属于典型的非粘性沉积物，而夏季表层沉积物砂、粉砂和黏土含量分别为24.4%、54.2%和21.4%，中值粒径为6.17 Φ，属于典型的粘性沉积物。因此，根据相应的计算公式可知，观测站位冬季和夏季表层沉积物临界侵蚀切应力分别为0.158 和0.107 N/m2。对比结果发现（图7d、8d），琅岐岛潮滩中上部地区近底层流速相对较小，由此产生的底部切应力也相对较小，对沉积物的活动性影响不明显，但冬季观测期间波浪引起的底部切应力显著大于表层沉积物的临界侵蚀切应力，表明冬季观测期间观测站位附近滩面发生强烈的侵蚀；夏季观测期间波浪引起的底部切应力的变化大部分时刻均小于表层沉积物临界侵蚀切应力，表明夏季观测期间观测站位附近滩面表层沉积物不易发生侵蚀。在长江口及江苏沿海潮间带观测结果也显示，在这种极浅水环境下波浪作用显著影响潮滩冲淤过程[28,58]。研究表明[59]，悬沙临界沉 降 切 应 力（τcrd）一 般 为0.06～0.1 N/m2，根 据Whitehouse 等的推荐，本文取τcrd=0.08 N/m2（其值适用于细颗粒泥沙），则对比结果显示（图8d），夏季观测期间大部分时刻由潮流和波浪引起的底部切应力小于临界沉降切应力，表明观测期间观测站位附近水体悬沙易发生沉降，导致滩面出现淤积。近底部悬沙净输运通量计算结果显示，冬季观测期间泥沙净向海输运，而夏季则净向岸输运。该结果表明，波浪是引起潮间带表层沉积物发生侵蚀的主要动力，而潮流则是搬运泥沙的主要动力。涨落潮不对称现象是河口海岸地区重要的水动力特征，是影响区域泥沙不对称输运的重要因素[60]，在河口海岸地区泥沙输运和地貌演变过程中扮演着十分重要的角色[57,61]。由图8 和图9 可知，冬季涨潮期间悬沙浓度总体大于落潮，但由于波浪引起的底部切应力显著大于悬沙临界沉降切应力，导致悬沙很少发生沉降，涨潮流带来的悬沙除部分被潮间带上部的互花米草盐沼植被捕获外，其他悬沙又在落潮流的作用下向海输运，而研究区无论冬季还是夏季，潮流均以落潮流占优，因此，落潮期间的底层悬沙输运通量总体大于涨潮，潮周期内表现为净向海输运；夏季，虽然是落潮流占优，但涨潮期间的底层悬沙浓度总体大于落潮，仅在大潮期间波浪作用较强时期出现落潮悬沙浓度大于涨潮的情况，因此涨潮期间的悬沙输运通量总体大于落潮，并且在多数时刻底部切应力小于临界沉降切应力，悬沙沿途发生沉降。该结论得到了本文现场冲淤观测结果的支持（图6），同时也与其他河口地区潮滩冲淤过程一致[52,62]。

图8 潮滩中上部光滩夏季近底部水动力时间序列变化过程a：淹没水深，b：近底部流速，c：有效波高，d：底部切应力，e：近底部悬沙浓度。Fig.8 Time-series of near-bottom hydrodynamic processes in summer on the mid-upper tidal flat a: Iinundation height, b: near-bottom mean current speed, c: significant wave height, d: bottom shear stress induced by tidal current and wave, e: near-bottom suspended sediment concentration.

图9 近底部瞬时悬沙水平输运率及潮周期内悬沙水平输运净通量a：冬季，b：夏季。正值为向海输运，负值为向岸输运；T1—T19 为观测期间的潮周期数。Fig.9 Seaward transport rate and net flux of near-bottom instantaneous suspended sediment in a tidal cycle in winter (a) and summer (b)Positive values indicate seaward transport, negative ones are landward transport; T1—T19: the number of tidal cycle during the observation.

本文观测结果显示，琅岐岛潮间带中部地区，冬季观测期间的潮流和波浪作用（图7b、c）均明显较夏季强（图8b、c），这主要一方面是由于研究区冬季主要受东北季风控制，而夏季主要受西南季风控制，冬季期间风速大于夏季，由此引起研究区冬季期间潮流和波浪作用较夏季强[43,63]；另一方面，则与潮间带地貌的季节性变化有关。现场调查发现，琅岐岛潮间带中部夏季被海三棱藨草覆盖，而冬季则变为光滩。现场观测和数值模拟结果均表明，盐沼植被具有显著的弱流、消波作用[64-66]。琅岐岛潮滩夏季由于生长了大面积的海三棱藨草，起到了一定的缓流作用，因此，夏季观测期间最大潮流流速和平均潮流流速均小于没有盐沼植被覆盖的冬季，波浪作用也是夏季显著小于冬季，表明琅岐岛潮间带海三棱藨草的生长具有一定弱流和显著的消波作用。根据不同站位冲淤变化对比分析（图6），在观测期间互花米草盐沼边缘滩面冲淤变化幅度很小，并且与其他两个站位的冲淤变化趋势没有明显的一致性，说明该地区滩面冲淤主要受水动力和互花米草的影响。位于海三棱藨草覆盖区域和光滩区域滩面冲淤变化总体表现为相反的变化趋势，这可能与海三棱藨草生长和消亡会引起潮间带中部地区冲淤变化有关，即在有海三棱藨草覆盖时，在动力-植被-泥沙的耦合作用下，涨潮带来的泥沙被植被捕获并沉积下来[30]，退潮水体悬沙浓度降低，导致在光滩地区退潮水流携沙能力有所增强，在一定程度上可能会引起潮间带下部光滩的侵蚀；当没有植被覆盖时，潮间带中部遭受显著侵蚀，侵蚀下来的泥沙在被落潮流携带至河口的过程中沿途发生沉降，引起潮间带下部滩面淤积。上述分析表明，盐沼植被在潮滩冲淤变化过程和潮滩发育演变中发挥着重要作用。

综上所述，琅岐岛潮间带中上部潮滩冲淤变化的动力学机制为：冬季，由于潮间带中部的海三棱藨草消失，潮间带被水淹没期间的波浪作用强，由此引起的底部切应力显著大于表层沉积物的临界侵蚀切应力，表层沉积物发生再悬浮，再悬浮起来的泥沙在潮流的作用下向海输运，导致滩面发生侵蚀；夏季潮间带中部被海三棱藨草覆盖，潮流和波浪均减弱，由此引起的底部切应力在大部分时刻小于表层沉积物的临界侵蚀切应力和临界沉降切应力，水体中的悬沙易发生沉降，导致滩面发生淤积。

5 结论

（1）琅岐岛潮滩表层沉积物粒度组分以粉砂和砂为主，黏土含量相对较少，呈现出显著的时空差异，夏季主要以粉砂和黏土组成的细颗粒沉积物为主，冬季则主要以砂和粉砂组成的粗颗粒沉积物为主。

（2）琅岐岛潮间带上部的互花米草盐沼总体以持续淤积为主，潮间带中部和下部表现出周期性的冲淤特征，总体以夏季淤积、冬季侵蚀为总特征，并且潮间带中部的冲淤变化幅度较潮间带下部大。

（3）沉积动力观测与计算结果表明，琅岐岛潮间带中上部近底部潮流流速相对较小，但波浪作用强，近底部悬沙浓度相对较高，其中冬季潮流流速、波浪作用及悬沙浓度均大于夏季，冬季潮周期内近底部悬沙以净向海输运为主，而夏季则以净向岸输运为主。

（4）琅岐岛潮滩冲淤变化主要受动力过程控制，冬季因潮间带海三棱藨草消失，波浪作用强，导致表层沉积物易发生再悬浮，再悬浮的泥沙在较强的潮流作用下净向海输运；夏季因潮间带有海三棱藨草生长，波浪作用较弱，水体中的悬沙易发生沉降，较弱的潮流将泥沙净向岸输运，导致潮间带发生淤积。

致谢：黄书仁、黄思添、余永泽、陈海煌、刘三善参与了野外采样工作，黄书仁参与了实验室样品分析，谨致谢忱！