江苏海岸中部近岸沉积物重力流输运过程

彭雲，于谦，万洁颖，杜智运，王黎，王韫玮，高抒

1.南京大学地理与海洋科学学院，南京 210023

2.华东师范大学河口海岸学国家重点实验室，上海 200062

3.南京师范大学海洋科学与工程学院，南京 210023

0 引言

重力流是一种高密度流体，由于其流体密度大于周围环境水体密度，能在重力的驱动下离岸下坡输运[1-3]。因触发机制与支撑机理的不同，重力流的类型多种多样，如浊流、碎屑流、异重流等[3-6]。本研究中的重力流，指的是波浪和潮流支撑的浮泥层运动[1-2]。波浪和潮流支撑的浮泥层运动是一种比较特殊的重力流，特殊在维持其输运的机理不是自悬浮，而是波浪和潮流提供的底部切应力。因此，相比于依靠自悬浮维持，需要较大坡度才能发生的浊流，波浪和潮流支撑的浮泥层运动在坡度很缓的海岸和陆架区域也能发生[1-2]。

过去30 年间，在世界范围内的许多海岸和陆架区域都发现了波浪和潮流支撑的浮泥层运动，如中国黄河口渤海湾海域[1]、南美亚马逊陆架[7-8]、美国Eel河邻域陆架[9-10]、意大利Po 河水下三角洲[11]、新西兰Waipaoa河邻域陆架[12-13]、和荷兰莱茵河海岸[14]等。波浪和潮流支撑的浮泥层运动能输运大量的沉积物跨过陆架[1-2]，对地貌演化产生重要影响，如大尺度水下三角洲和倾斜地层的演化[15]，以及中部陆架泥质沉积的形成[16]。因此，研究波浪和潮流支撑的浮泥层运动这一重力流输运过程具有重要意义。

然而，波浪和潮流支撑的浮泥层（重力流）的发生具有偶然性，其持续时间也很短暂[1-2]。重力流的这两个固有特性，导致在野外现场捕捉其输运过程一直以来都是一个巨大挑战。已有研究表明，充足的细颗粒沉积物物源供给和较强的水动力是触发重力流的必要条件[2]，尽管最近在砂质海岸近岸海域也观测到了重力流[14]。江苏海岸中部近岸海域素以潮流强劲著称，最大潮差可达9.39 m[17]。自黄河1855年北归以后，江苏北部的废黄河三角洲遭受到严重侵蚀，侵蚀下来的沉积物在苏北沿岸流的作用下向南输运，为江苏海岸中部近岸海域提供了大量的沉积物供给[18-19]。因此，江苏海岸中部近岸海域是研究沉积物重力流输运过程的绝佳场所。

Penget al.[20]在江苏海岸中部近岸海域，于不同季节实施了三次野外观测（2019年春季、2019年秋季和2020 年秋季），均捕捉到多次（大于15 次）重力流输运过程，其研究表明重力流发生的频率可能更高[20]。然而，Penget al.[20]观测到的数次重力流的触发机制具有多样性，对重力流的维持和停止机理也没有进行详细的介绍。长久以来，触发、维持和停止重力流的准确机理，一直都是重力流研究的难点和热点之一[2,13]。因此，有必要依靠更多的观测数据，更加细致地研究重力流的触发、维持和停止机理。

Penget al.[20]的研究表明重力流可能有着更加复杂的动力行为。最普遍用于描述重力流动力行为的理论分析模型是由Wrightet al.[1]首次提出的浮力—阻力模型。过去二十年间，不论是对重力流野外观测数据的分析[11-14,21]，还是对重力流数值模拟[22-24]，以及水槽实验[25]研究，Wrightet al.[1]的浮力—阻力模型被多次运用，模型中的参数也被多次验证。Penget al.[26]继续分析了Penget al.[20]中2020年秋季的观测数据，运用在2020 年秋季观测到的23 次重力流事件，对Wrightet al.[1]的浮力—阻力模型中的两个关键参数（底部拖曳系数CD和整体理查德森数Rib）进行了验证。结果显示，CD明显大于前人研究报道的0.001～0.006 的区间范围，Rib接近于模型给出的临界值0.25。然而，也有关于重力流的研究表明，Rib不一定是0.25，而是维持在0.01 附近，比0.25 要小一个数量级[11,14]。因此，有必要基于重力流的野外观测数据，对Wrightet al.[1]的浮力—阻力模型中的重要参数（CD和Rib）进行验证。

研究的主要目的：（1）展示2018年秋季在江苏海岸中部近岸海域观测到的重力流事件；（2）基于现场观测到的重力流事件，分析重力流的形成、维持和停止机理；（3）基于重力流观测数据，对重力流进行参数化，检验重力流浮力—阻力模型中重要参数（CD和Rib）的敏感性。

1 观测区域

研究区域位于江苏海岸中部近岸海域（图1）。1128—1855 年，黄河在700 多年的时间里从江苏北部入海，黄河携带的巨量沉积物造就了广阔的水下三角洲（向海延伸到40 m水深）和江苏沿岸的滨海平原[18]。自1855 年黄河北归以后，由于沉积物供应急剧减少，1128—1855 年由黄河输运而来的沉积物所形成的三角洲遭受到强烈的侵蚀，开始急剧后退，形成废黄河三角洲。不断受侵蚀的苏北废黄河三角洲提供了巨量的沉积物，在苏北沿岸流的作用下向南输运，为江苏中部海岸近岸海域提供了大量的沉积物供给[18-19]。

图1 观测区域示意图（a）区域概况，红色矩形表示研究区域；（b）野外观测站位概况；（c）使用多波束于2018年8月测量的跨岸地形剖面（子图b中的虚线）；子图（b）中的黑色箭头指示沿岸和跨岸流速正方向；子图（c）的横坐标S和纵坐标h分别表示向海距离和水深；蓝色三角形表示三脚架的放置位置Fig.1 Map of observation area(a) inset showing general setting of the Jiangsu coast,China;the red rectangle denotes the study area;(b) field measurement site;and (c) cross-shore profile measured by a multibeam survey in August 2018 (dashed line in (b));the black arrows in (b) denote the positive alongshore and crossshore current directions;the abscissa and ordinate in (c) denote distance seaward and water depth;the blue triangle denotes the location of the tripod

由于沉积物供给丰富，研究区是中国最浑浊的海岸之一[27-28]。研究区悬沙浓度在冬季最高，夏季最低。自潮滩向岸外浅水区，悬沙浓度不断减小[28]。研究区近岸10 m 水深范围内，底质类型主要是砂和砂质粉砂，平均粒径介于50.8～133.9 µm[29]。

研究区主要受两大潮波系统的影响，分别是来自南黄海的旋转潮波和来自东海的前进波[15,30]。研究区域的潮汐类型为正规半日潮，平均潮差介于3.9～5.5 m。潮流以涨潮流为主导，最大流速超过1.5 m/s，沿岸流速相对于跨岸流速占据主导地位。研究区余流较为复杂，大小潮变化较大，余流值介于0.007～0.370 m/s，余流方向变化趋势不明显。受岸外辐射沙脊群的阻挡，研究区波浪相对较弱，有效波高大于2 m的出现频率仅为5%[18]。

2 研究方法

2.1 野外观测

野外观测时间为2018 年11 月25 号7 点00 分至2018年11月29号8点30分，覆盖8个潮周期。观测站位位于江苏中部海岸岸外水下岸坡（33.415° N，120.748° E），该观测站位附近，沿岸方向上的地形起伏不大，而垂直于海岸的跨岸方向上的地形有明显坡度变化。近岸6 km，平均水深小于13 m的范围内，没有发育明显的沟槽。观测站位附近潮滩上的潮沟尺度均相对较小，并没有延伸到岸外水下岸坡[18]。2018年8月4日，在江苏海岸中部斗龙潮滩岸外水下岸坡区域，利用多波束随机测量了一条跨岸地形剖面（图1c），结果显示，近岸5 km，平均水深小于12 m的范围内，跨岸方向上的地形有明显坡度变化，但并没有发育明显的沟槽。

野外观测采用底基三脚架进行，三脚架放置在由多波束于2018年8月4日测量到的跨岸海底地形上，平均水深约7 m（图1b，c）。海底三脚架具有能够搭载多种水文仪器，长期、连续、稳定和自容进行几乎无干扰原位观测的优点。三脚架上搭载一系列水文仪器，用于长时间原位连续观测水体的流速和浊度等信息。

在距底1.25 m 处放置1 个上视声学多普勒流速剖面仪（ADP，1 MHz），用于测量上覆水体的流速剖面，采样分辨率为0.2 m 一层，采样频率为2 Hz。在距底0.3 m 处放置1 个声学多普勒流速仪（ADV），用于测量近底三维高频（16 Hz）流速，每15 min 的采样间隔内测量128 s。水深由ADV 的压力记录经过压力校正后得到。底床高度变化由ADV换能器与底床之间的相对高程得到[31-32]。

水体浊度由AQUAlogger 310 TY 和光学后向散射仪（OBS-3A）这两种浊度仪器测量，在距底0.1 m处放置1 个AQUAlogger 310 TY，在距底0.3 m 和1 m 处分别放置一个OBS-3A。三个浊度仪器的采样参数均设置为以1 Hz 频率每3 min 测量30 s。浊度仪器的详细标定流程见Wanget al.[33]。三脚架上搭载的所有仪器的采样参数设置详情见表1。观测时间段内风场信息选自欧洲中期天气预报中心（ECMWF）第五代全球大气再分析资料（ERA5）再分析近地面10 m 风速数据。ERA5 再分析风速数据具有时空分辨率高的特点，时间分辨率为1 h，空间分辨率达到0.25°×0.25°，能够较好地提供研究区域的背景风场信息。

表1 野外观测仪器的采样参数设置详情Table 1 Deployment scheme for the measurements

2.2 数据处理

基于ADV的压力记录，通过谱分析方法[34]，计算得到有效波高（Hs），波浪谱峰周期（Tp），按照线性波理论，计算得到波浪轨道速度（Uw）。

浪致底部切应力的计算方法如下[35]：

式中：τw是浪致底部切应力，ρw是海水密度，fw是波浪摩擦系数，计算方式为[35]：

式中：A=UwTp/2π 是波浪的半轨道偏移量（semiorbital excursion），Z0是摩阻高度。

将ADV 和ADP 测量的流速分解到沿岸和跨岸方向，正方向分别为沿岸向南和跨岸向海（图1b），沿岸正方向为从正北顺时针旋转157°。海上实测的水流，包括周期性潮流和余流两部分，余流是指从实际海流中去掉周期性潮流后的剩余部分，对水体各层的沿岸和跨岸流速的调和分析得到的非周期项即为余流，对垂向上各个层位的余流流速进行平均即得这些层位的垂向余流平均值。

潮致底部切应力可以由雷诺应力法得到[31,36]：

式中：τcx和τcy分别是平行和垂直于主流向方向的切应力，u′，v′和w′是距底某一高度处水平（平行和垂直于主流向方向）和竖直流速的瞬时湍流部分，可以表示为测量的瞬时流速与测量间隔平均流速之差（比如：′），上横线表示测量间隔的时间平均。由此：

式中：τc是潮致底部切应力。

3 结果

3.1 时间序列的观测

观测得到的时间序列数据如图2 所示。水深在4.43～8.86 m 范围变化，平均水深为6.94 m（图2a）。风速最大值小于5.00 m/s，平均风速为2.82 m/s，主要风向为西北（NW）（图2b）。距底0.1 m、0.3 m 和1 m处悬沙浓度的变化范围分别为1.17～15.75 kg/m3、0.73～7.88 kg/m3和0.67～3.78 kg/m3（图2c）。距 底0.3 m处的沿岸流速介于-0.47～0.87 m/s（图2e），而距底0.3 m处的跨岸流速介于-0.24～0.22 m/s（图2d），说明沿岸流速相对于跨岸流速占据主导地位。距底0.3 m 处的沿岸余流和跨岸余流流速分别为0.07 m/s和-0.01 m/s，上覆水体（＞1.25 m）垂向平均的沿岸余流和跨岸余流流速分别为0.13 m/s 和-0.04 m/s。在观测中发生了一次强度相对较小的波浪事件，持续时间约31 h（图2f 中的黑色矩形）。波浪事件前，有效波高小于0.25 m（图2f），波浪轨道速度小于0.05 m/s（图2d），浪致底部切应力小于0.15 Pa（图2g）。随着波浪事件的发生，有效波高、波浪轨道速度和浪致底部切应力均先增大后减小，它们的最大值分别为0.89 m、0.16 m/s和0.68 Pa（图2e～g）。波浪谱峰周期介于3.46～9.33 s，平均值为5.10 s（图2f）。潮致底部切应力变化范围为0.03～3.02 Pa，最小值大致出现在涨憩和落憩相位，最大值大致出现在涨急和落急相位，且涨急相位附近的最大值通常要比落急相位附近的最大值更大（图2g）。底床高度在0.03～0.13 m范围变化（图2h）。

3.2 浮泥层与重力流事件

时间序列的悬沙浓度清楚地揭示了浮泥层事件的存在（图2c）。波浪事件前和波浪事件期间，均观测到近底高悬沙浓度（＞10 kg/m3）现象。在研究观测区域，已有研究报道过浮泥层的存在，并将区分近底高浓度浮泥层和上覆低浓度水体的悬沙浓度临界值取为5 kg/m3[20]。以5 kg/m3为临界值，观测期间一共出现了8 次短暂的浮泥层事件，持续时间为0.25～2.75 h，按出现的时间顺序，依次标记为E1～E8（图2）。其中，发生在波浪事件前的4 次浮泥层事件均出现在低水位时段（E1～E4）。而发生在波浪事件期间的4次浮泥层事件，有2次出现在低水位时段（E5，E7），2 次出现在高水位时段（E6，E8）。发生在低水位和高水位时段的浮泥层事件分别用黄色和绿色阴影区域表示（图2）。发生在低水位时段的浮泥层事件（E1～E5，E7）包含晚落潮和始涨潮相位，发生在高水位时段的浮泥层事件（E6，E8）仅包含始落潮相位。

距底0.3 m处的沿岸流速（UA0.3）具有明显的潮变化特征，涨潮时流向向南，落潮时流向向北（图2e）。距底0.3 m处的跨岸流速（UC0.3）和深度平均的上覆水体跨岸流速（UCU＞1.25 m）也具有明显的潮变化特征，涨潮时流向向陆，落潮时流向向海（图2d）。然而，对应于浮泥层事件发育期间的某些时段，UC0.3会出现异常（图2d），具体时间段为：E2 和E5 的始涨潮相位（early flood，缩写为ef），E6和E8的始落潮相位（early ebb，缩写为ee），E7 的晚落潮相位（late ebb，缩写为le），详细解释如图3所示。

图3 跨岸流速（UC）和悬沙浓度（SSC）剖面（a～e）上覆水体跨岸流速剖面（蓝色点）和距底0.3 m处的跨岸流速（UC0.3：红色点）；（f～j）近底层悬沙浓度剖面。从左至右的五列子图分别表示浮泥层事件E2的始涨潮（ef）、E5的始涨潮（ef）、E6的始落潮（ee）、E7的晚落潮（le）和E8的始落潮（ee）相位；流速正方向为跨岸向海Fig.3 Vertical profiles of cross-shore current velocity (UC) and suspended sediment concentration (SSC)(a-e) UC profiles of the upper water layers (blue dots) and cross-shore current velocity at 0.3 m asb (UC0.3: red dot);(f-j) SSC profiles of the near-bottom layers.From left to right,the first to fifth columns correspond to early flood (ef) tide of E2;early flood (ef) tide of E5;early ebb (ee) tide of E6;late ebb (le) tide of E7;and early ebb (ee) tide of E8,respectively.Positive current speed denotes the offshore direction

在浮泥层事件E2和E5的始涨潮相位（ef），上覆水体的跨岸流速符合潮变化，流向向陆，但UC0.3的方向却是离岸向海的，大小约0.02 m/s（图3a，b，f，g）。在浮泥层事件E6的始落潮（ee）、E7的晚落潮（le）、和E8 的始落潮（ee）相位，上覆水体的跨岸流速和UC0.3均符合潮变化，流向向海，但是，UC0.3值却非常大，约0.2 m/s，明显大于上覆水体整个剖面的跨岸流速（图3c～e，h～j）。区别于上覆水体的异常离岸流速和浮泥层事件同时发生，是重力流发生最明显和最直接的证据[11,14,20,32]。因此，在浮泥层事件E2 的始涨潮（ef）、E5 的始涨潮（ef）、E6 的始落潮（ee）、E7 的晚落潮（le）、和E8的始落潮（ee）相位，基于现场观测的流速和悬沙浓度证据直接显示了重力流的存在。

4 讨论

4.1 重力流证据

研究观测到8次短暂的浮泥层事件（图2c），仅在其中5次发现了重力流输运过程（近底异常离岸流速和浮泥层同时发生）（图2c，d）。而且，近底异常离岸流速（UC0.3）并不是在浮泥层持续时间段内均有，而是发生在浮泥层发育时间段内的某个相位，如E2事件的始涨潮相位（图3）。

在E1、E3 和E4 事件发育期间，均未观测到类似E2 事件期间的异常离岸流速（图2c，d）。可能的解释是，E1 和E4 事件的持续时间太短，约半小时（图2c），重力流离岸流速不容易被观测到，或是浮泥层还无法形成重力流，就已经被潮流的作用停止[1,32]。而在E3 事件的始涨潮相位，虽然没有观测到如E2事件始涨潮相位UC0.3一般的离岸流速，但此阶段的UC0.3具有类似E2 事件始涨潮相位UC0.3的特征，即存在微小的凸起（图2c，d）。同样地，在E7 事件始涨潮相位也没有观测到异常的UC0.3，但此阶段的UC0.3也存在微小的凸起（图2c，d）。由此推测，在E3 和E7 事件的始涨潮相位，也可能发生了重力流，但由于重力流流速所处的垂向区间小于0.3 m，从而没有被捕捉到[1,32]。

在E6和E8事件的始落潮相位，均观测到异常离岸流速（图3c，e，h，j），且均发生在SSC0.3大于5 kg/m3期间（图2 中的绿色阴影区域）。由此推测，在E6 和E8事件发育期间，重力流持续的时间可能更久，但由于其厚度小于0.3 m而未被观测到[1,32]。

在研究观测区域，已有研究报道过发生在低水位和高水位时段的重力流事件[20,26]，近底高分辨率流速剖面（距底0.3 m，分辨率为0.01 m一层）显示，近底异常离岸流速并不是仅发生在浮泥层发育期间的某一个相位，而是在浮泥层发育的整个时间段均有（见Penget al.[20]中的图5和Penget al.[26]中的图3）。因此，在研究观测到的浮泥层事件发育期间，极有可能均发生了重力流输运过程，只不过由于重力流厚度小于0.3 m，导致在浮泥层事件发育期间的某些时段未观测到重力流流速。

4.2 浮泥层与重力流的形成、维持和停止机理

浮泥层形成最可能的两种机制是波浪导致的底床沉积物再悬浮[2,8-14]，和落憩或涨憩时段上覆水体沉积物沉降[2,37-38]。波浪事件前，对发生在低水位时段的浮泥层事件E2，波浪轨道速度（Uw）小于0.03 m/s（图4e），说明波浪对底床沉积物的再悬浮作用是可以忽略的[1]。因此，潮的作用是首要考虑因素。落憩相位前，当距底0.3 m 处的沿岸流速（UA0.3）减小到约-0.17 m/s 时（图4 左边列子图的第一条黑色竖直虚线），距底0.1 m处的SSC（SSC0.1）开始快速增大（大于5 kg/m3），而距底0.3 m（SSC0.3）和1 m（SSC1）处的SSC却不断减小（图4a，c）。落憩相位前后，即UA0.3从-0.17 m/s 变化到0.20 m/s 期间（图4左边列子图的第一条和第二条黑色竖直虚线之间），SSC0.1增大到15.75 kg/m3的最大值后几乎保持不变，SSC0.3先减小后增大，SSC1一直减小，且SSC0.3和SSC1均小于5 kg/m3，说明上覆水体的沉积物沉降形成了约0.1 m厚的浮泥层[37-38]。

图4 浮泥事件E2 和E5 的时间序列（a，b）水深（h）与距底0.1、0.3和1 m处的悬沙浓度（SSC）；（c，d）距底0.3 m处的沿岸（UA0.3）和跨岸（UC0.3）流速；（e，f）波浪轨道速度（Uw）；流速正方向为沿岸向南和跨岸向海Fig.4 Observed time-series of fluid-mud events E2 and E5(a,b) water depth (h) and SSCs at 0.1,0.3 and 1 m above sea bed (asb);(c,d) along-shore (UA0.3) and cross-shore (UC0.3) velocities at 0.3 m asb;and (e,f) wave orbital velocity (Uw).Positive current speeds denote offshore and southward alongshore directions

当UA0.3＞0 m/s时，SSC0.3开始增大，SSC1仍在减小。当UA0.3＞0.20 m/s 时，SSC0.3和SSC1均不断增大（图4a，c），说明潮流不断悬浮近底浮泥层内的沉积物向上扩散[12]。SSC0.3增大到大于5 kg/m3且一直大于5 kg/m3的时间范围内，UC0.3发生异常，具有约0.02 m/s的离岸流速（图4左边列子图的黄色阴影区域），说明潮流导致的再悬浮触发了重力流，并维持其运动[11-14,20]。当UA0.3增加且大于0.40 m/s 后（图4 左边列子图的第三条黑色竖直虚线），SSC0.1、SSC0.3和SSC1均快速降低（图4a，c），指示了在潮流的作用下，浮泥层内沉积物不断向上扩散，直至浮泥层消亡[1-2,12]。

需要指出的是，针对E2事件期间SSC0.1不随时间变化的区间（图4a），最可能的原因是真实浊度值超过了浊度仪的量程，从而导致浊度仪测量到的浊度一直维持在量程最大值附近[39]。因此，真实的SSC0.1可能更高。

浮泥层事件E1、E3和E4也发生在波浪明显变大前的低水位时段，且这些事件发育期间的近底SSC的变化情况同E2事件发育期间的近底SSC的变化情况类似（图2c）。因此，可以认为，对发生于波浪事件前低水位时段的浮泥层事件E1到E4，形成机理均为上覆水体的沉积物沉降[37-38]。

与发生在波浪事件前的E1 到E4 事件相比，E5和E7 事件也发生在低水位前后，但E5 和E7 发生在波浪事件期间（图2）。与E2事件相比，在E5事件落憩相位前后，即UA0.3从-0.17 m/s变化到0.20 m/s期间（图4右边列子图的第一条和第二条竖直虚线之间），SSC0.1同样是先急剧增大到最大值15.75 kg/m3后几乎保持不变，SSC1不断减小，然而，SSC0.3却是先增大后减小（图4b），且E5 事件期间的Uw大于0.04 m/s（图4f）。因此，除了上覆水体的沉积物沉降这一触发机理外，波浪再悬浮对E5 事件的形成也有影响[9-12,16]。在E5 事件的始涨潮相位，即UA0.3从0.20 m/s 变化到0.40 m/s期间，近底SSC的变化情况同E2事件始涨潮相位近底SSC 的变化情况类似，且观测到了约0.02 m/s 的异常离岸流速（图4b，d）。而且，在E5 事件始涨潮相位，Uw从0.04 m/s增大到大于0.10 m/s（图4f）。因此，E5事件始涨潮相位的重力流是由潮流和波浪共同导致的底部切应力维持[1-2,9-14]。

虽然均是发生在波浪事件期间的低水位前后，但与E5 事件相比，在E7 事件的UA0.3从-0.15 m/s 变化到0.15 m/s 期间（图5的第一条和第二条竖直虚线之间），SSC0.1不断增大，SSC0.3和SSC1均是先增大后减小（图5a），显现出很强的再悬浮特征。此阶段的Uw大于0.05 m/s（图5c），且SSC0.1、SSC0.3和SSC1增大期间，Uw大于0.08 m/s（图5）。因此，对E7 事件，波浪导致的底床沉积物再悬浮是其形成的主要机理，上覆水体的沉积物沉降是次要机理[9-13]。E7 事件始涨潮相位的近底SSC 变化情况与E2 和E5 事件始涨潮相位的近底SSC 变化情况大致相似，即随着UA0.3不断增加，潮流不断地悬浮近底约0.1 m厚的浮泥层内沉积物，并致其向上扩散，浮泥层因而消亡（图5a，b）[32]。

图5 浮泥事件E7 的时间序列（a）水深（h）与距底0.1 m、0.3 m和1 m处的悬沙浓度；（b）距底0.3 m处的沿岸（UA0.3）和跨岸（UC0.3）流速；（c）波浪轨道速度（Uw）；流速正方向为沿岸向南和跨岸向海Fig.5 Observed time-series of fluid-mud event E7(a) water depth (h) and SSCs at 0.1 m,0.3 m and 1 m above sea bed (asb);(b) along-shore (UA0.3) and cross-shore (UC0.3) velocities at 0.3 m asb;and (c) wave orbital velocity (Uw).Positive current speeds denote offshore and southward alongshore directions

波浪事件期间，浮泥层事件不仅发生在低水位时段，也发生在高水位时段，即E6 和E8（图2）。在E6 和E8 发育期间，SSC0.1和SSC0.3均是先增大后减小（图6a，b），呈明显的再悬浮特征。E6事件期间的Uw均大于0.09 m/s（图6e），E7 事件期间的Uw均大于0.05 m/s（图6f）。因此，E6 和E8 均是由波浪再悬浮底床沉积物而形成[9-11,13]。在E6 和E8 事件的始落潮相位，近底SSC的变化情况与E2、E5 和E7 事件始涨潮相位的近底SSC变化情况大致相似，即随着UA0.3不断增加，潮流不断地悬浮近底0.1～0.2 m 厚的浮泥层内的沉积物，并致其向上扩散，浮泥层因而消亡（图6）[32]。

图6 浮泥事件E6 和E8 的时间序列（a，b）水深（h）与距底0.1 m、0.3 m和1 m处的悬沙浓度（SSC）；（c，d）距底0.3 m处的沿岸（UA0.3）和跨岸（UC0.3）流速；（e，f）波浪轨道速度（Uw）；流速正方向为沿岸向南和跨岸向海Fig.6 Observed time-series of fluid-mud events E6 and E8(a,b) water depth (h) and SSC at 0.1,0.3 and 1 m above sea bed (asb);(c,d) along-shore (UA0.3) and cross-shore (UC0.3) velocities at 0.3 m asb;and (e,f) wave orbital velocity (Uw).Positive current speeds denote offshore and southward alongshore directions

综上，对观测到的浮泥层与重力流的形成、维持和停止机理总结如下：低水位时段上覆水体的沉积物沉降会形成约0.1 m 厚的高浓度浮泥层（E1～E4），波浪再悬浮和上覆水体沉积物沉降的共同作用也会导致浮泥层事件在低水位时段形成（E5，E7），厚约0.1 m，波浪的再悬浮作用会导致浮泥层事件在高水位时段形成（E6，E8），厚度为0.1～0.2 m。浮泥层事件一旦形成后，能在重力的作用下离岸下坡输运，依靠潮致的底部切应力（E2），或波浪和潮流共同导致的底部切应力维持（E5～E8）。随着沿岸流速的增大，近底0.1～0.2 m 厚的浮泥层会被冲淡，直至消亡。以上基于研究观测数据的推论，与重力流数值模拟研究得到的结果类似[40]。

4.3 理论模型运用

对于波浪和潮流支撑浮泥层这一比较特殊的重力流类型，Wrightet al.[1]基于观测数据，首次提出了理论分析模型。近20 年来，该理论分析模型被广泛用于重力流参数化和动力行为研究[11-14,21-25]。Wrightet al.[1]的理论分析提出，重力流所受的重力驱动的离岸下坡梯度力和底部阻力处于平衡态，可表示为：

式中：左边和右边分别代表浮力异常项和阻力项；B=gsCH是浮力异常；g 是重力加速度；s是沉积物相对于海水的比重，取值为1.6；C=Cs/ρs是基于浮泥层厚度平均的体积浓度；Cs是基于浮泥层厚度平均的质量浓度；ρs=2 650 kg/m3是沉积物密度；H是浮泥层的厚度；sina是底床坡度。CD是底部拖曳系数，取值0.003，，其中，Uw、Vc和Ug分别为波浪轨道速度、浮泥层顶部的沿岸流速和重力流流速。

整体理查德森数Rib是Wrightet al.[1]理论分析模型的另一种描述方式：

一般来说，B和Ug难以在现场观测中直接获得，尤其是在早期的野外观测中。因此，通常假定沉积物物源供给是无限的，以Rib等于0.25 的临界值作为重力流处于平衡态的边界。反馈机制如下：如果Rib＜0.25，增加的湍流将会悬浮更多的沉积物，Rib增大到0.25。反之，如果Rib＞0.25，悬沙浓度梯度导致的沉积物分层将会抑制湍流，沉积物沉降，从而Rib降低到0.25[1,12,20-23]。

研究中，近底可获得的最低层流速是距底0.3 m处的，可获得的近底悬沙浓度是距底0.1 m、0.3 m 和1 m处的。因此，对观测到的重力流事件（距底0.3 m处的异常离岸流速和浮泥层同时发生），根据图3～6，取浮泥层厚度H=0.3 m，Ug=0.2 m/s，Vc=0.24 m/s，Uw=0.09 m/s，观测站位的坡度为0.007 2（图1c）。

重写公式（6）得到：

将上述基于观测数据选取好的参数代入公式（8），得到Cs=15.26 kg/m3。重力流发育期间，实测距底0.1 m处的悬沙浓度约为15 kg/m3（图3）。因此，由理论模型得到的距底0.3 m 内深度平均的质量悬沙浓度与实测结果高度吻合。

或者，对观测到的重力流事件，假设湍流混合和沉积物分层之间的反馈达到平衡，即Rib=0.25。重写公式（7）得到：

将Rib=0.25 的临界值和其他上述基于观测数据选取好的参数代入公式（9），得到Cs=14.88 kg/m3，同样与观测结果高度吻合。

拖曳系数（CD）和整体理查德森数（Rib）是重力流参数化研究中的关键参数[11-14,22-25]。近底高浓度浮泥层内的高分辨率悬沙浓度和流速剖面的野外观测，一直以来都是重力流研究的难点之一，尤其是在早期的野外观测中，从而浮力异常B和重力流流速Ug难以从野外观测中直接获得[1-2]。因此，许多重力流参数化的研究，假设沉积物的供给是无限的，以保证重力流始终处于临界分层状态，从而将Rib等于0.25的临界值用于实现模型闭合[9-10,19-20]。模型给的CD等于0.003 的参考值，也是当Wrightet al.[1]的理论浮力—阻力模型被用于重力流野外观测数据分析[12-13,21]和数值模拟研究[22-24]时，最常运用的CD值。

基于江苏海岸重力流观测数据，对重力流进行参数化后，不仅提供了又一个野外观测实例支持理论模型给的CD等于0.003和Rib等于0.25 的假设，也说明当Wrightet al.[1]的理论浮力—阻力模型被用于野外重力流参数化[1,12-13,21]和重力流数值模拟[22-24]研究时，Rib取值0.25和CD取值0.003的合理性。

5 结论

（1）2018 年秋季，在江苏海岸中部近岸海域实施的野外观测中，捕捉到了8 次高浓度浮泥层事件，其中5次发生了短暂的重力流输运过程。

（2）重力流由上覆水体沉积物沉降或波浪再悬浮底床沉积物形成，由潮致底部切应力或波浪与潮流共同导致的底部切应力维持，近底高浓度浮泥层内的沉积物在潮流的作用下，会不断向上扩散，重力流从而消亡。

（3）研究区域观测到的重力流的动力行为支持经典重力流浮力—阻力模型的使用。

致谢 南京大学樊一阳和林杭杰参加了野外工作和室内标定，广东海洋大学祁雅莉帮助绘图。谨致谢忱！