南黄海西部日照近海悬浮泥沙分布、输运及控制因素

2022-07-27卢鹏飞岳英洁朱龙海胡日军尹砚军冷星

海洋地质与第四纪地质 2022年3期

关键词：盐度泥沙站位

卢鹏飞，岳英洁，朱龙海,3，胡日军,3，尹砚军，冷星

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

2.山东省海洋工程咨询协会，济南 250000

3.中国海洋大学海底科学与探测技术教育部重点实验室，青岛 266100

悬浮泥沙是沉积物由陆地向海洋输运的主要载体，同时悬浮泥沙在运移过程中会携带各类营养盐和污染物，因此悬浮泥沙的输运和沉积不仅能反映海洋沉积动力过程，而且也深刻影响着陆架海海洋生态环境和资源分布[1-5]。

针对南黄海海域悬浮泥沙分布、运移和物质来源的研究，前人利用实测资料分析[6-9]、数值模拟[10]、遥感资料反演[11-12]等方法开展了相关研究。如李文建等[13]基于现场激光粒度仪(LISST)观测数据研究了南黄海的悬浮物粒度分布特征；仲毅等[14]通过水样抽滤得到的悬浮体质量浓度数据和现场激光粒度仪(LISST)观测到的悬浮体体积浓度数据研究了南黄海35°断面的悬浮体浓度垂直分布及其季节变化。李建超等[15]基于现场激光粒度仪(LISST)观测数据证明南黄海底部存在多个悬浮体高浓度区域，沿岸流控制下的近岸悬浮体高浓度区，悬浮体多为细颗粒物质；黄海暖流流经海域浓度相对较低，悬浮体多为粗颗粒物质；海洋温跃层会抑制底层悬浮体向表层运移，促进其在底部富集。部分学者对悬浮泥沙输运机制的研究表明，不同海域悬浮泥沙输运机制存在差异，如鸭绿江河口地区、苏北近岸西洋水道、洋山港海域和长江口悬浮泥沙输运以平流输运为主，主要控制因素是欧拉余流[16-19]；天津港和福宁湾附近海域悬浮泥沙净输运以潮泵效应为主[20-21]；椒江河口悬浮泥沙输运以潮泵效应和平流输运为主[22]。

前人对南黄海海域的悬浮泥沙分布及输运取得了丰富的研究成果，但对日照近岸海域悬浮泥沙分布及输运的研究稍显匮乏。日照近岸海域位于南黄海西部，水体稳定，地表径流对其影响微弱，是研究波浪及海流作用对悬浮泥沙分布及输运影响的理想区域[23]。鉴于此，本文以日照近岸海域为研究区（图1），通过开展多站位、大比例尺的同步实测潮流、温盐和悬浮泥沙调查，运用通量机制分解法，研究日照近岸海域悬浮泥沙的平面和垂向分布特征，探究不同水动力因素对悬浮泥沙输运贡献率的差别。研究成果对完善南黄海近岸海域悬浮泥沙分布及输运机制，揭示日照近岸海域海底地貌演变机制，指导海岸工程建设具有一定的理论和现实意义。

1 研究区概况

研究区位于南黄海西部的日照海岸带地区，岸线以砂质海岸为主。研究区底质类型以砂和粉砂为主，东部、中部、西部海域多见砾质泥质砂和含砾泥质砂，沉积动力较强，以滚动和跳跃组分为主；岚山港和日照港附近海域多见含砾泥、砂质泥和泥等细粒物质，沉积动力较弱，以悬浮组分为主[24]。日照海域潮汐类型属规则半日潮，平均海平面3.10 m，历年最高潮位5.40 m，历年最低潮位0.27 m，平均高潮位4.22 m，平均低潮位1.24 m。波浪以风浪为主，常浪向、强浪向均为E向；次常浪向为ESE向和SE向[25]。

2 数据来源与研究方法

2.1 数据来源

2020年10月18 日10：00—19日11：00（大潮期），中国海洋大学在日照近岸海域进行了5个站位25 h的海流连续同步观测；2020年10月20—23日（秋季），中国海洋大学在日照东部海域进行了35个站位的大面站水文泥沙现场调查，调查内容包括水体浊度和水体温盐，观测期间海况良好，调查站位如图1所示。

图1 日照海域调查站位图Fig.1 Location map of survey stations in Rizhao Sea area

（1）海流和潮位

采用美国LinkQuest公司生产的FlowQuest 600 kHz多普勒海流剖面仪（ADCP）和美国Teledyne RD Instruments公司生产的WHS 600 多普勒海流剖面仪（ADCP）采集流向、流速、水深数据，ADCP于船右舷后侧位置固定，每隔10 min采集一次数据。观测剖面按0.5 m分层，根据实际水深按6点法（表层、0.2H、0.4H、0.6H、0.8H、底层）进行数据提取。

潮位观测采用美国Global Water生产的DCX-22-2潮位仪和中国中船海鹰加科海洋技术有限责任公司生产的海鹰HY1300全数字潮位仪，采样间隔1 s。

（2）悬浮泥沙

5个连续站悬浮泥沙水样的采集和海流观测同步进行，采样间隔1 h，根据水深在表（水面以下1 m）、中、底（海底以上1 m）3层各采集1 L水样，用于悬浮体质量浓度的测定。对现场采集的水样在实验室进行抽滤、烘干、称重。过滤膜为孔径0.45 μm、直径47 mm的醋酸纤维膜，在电热培养恒温箱中以50℃烘干，恒温6 h后放入硅胶干燥器，放置6 h。称重在精度1/10万的天平进行。为保证结果准确，悬浮体浓度数据利用空白校正膜进行了校正。

水体浊度（FTU）的现场测量使用英国Aquatec公司生产的Aqualogger 310TY型浊度仪和美国D&A公司生产的OBS-3A型浊度仪完成，采样频率为 1 Hz，量程（0～200 FTU），测量精度为±0.01 FTU。

对实测浊度数据与悬浮泥沙浓度实验结果进行相关性分析，结果表明二者相关系数达到0.85 以上（图2），证明水体浊度是悬浮泥沙浓度的良好替代性指标，本文已根据相关性关系将浊度换算成悬浮泥沙浓度。

图2 实测浊度与悬浮泥沙浓度相关曲线图Fig.2 Correlation curve of measured turbidity and sediment concentration

（3）温度和盐度

水体温度（℃）、盐度（PSU）、水深（m）的现场测量使用美国TRDI公司生产的CTD-NV型温盐深仪完成，采样频率5 Hz，温度和盐度测量精度为±0.005 ℃ 和±0.005 PSU。

2.2 研究方法

（1）通量机制分解法

本文采用Dyer[26]的物质通量计算方法将悬浮泥沙通量分解为多个动力项，从而研究海流作用下各动力项所对应的主要动力因子对研究区悬浮泥沙输运的作用。依据Ingram[27]和Uncles 等[28]提出的用相对水深进行分解瞬时的物质输运量的方式，潮周期Tt平均瞬时单宽悬浮泥沙输移通量T的计算式为：

式中：T1为欧拉余流控制下发生的悬浮泥沙输运，T2为斯托克斯漂移输运量，表示潮波与潮流相关项；T1与T2之和为拉格朗日余流平流输运贡献项；T3为潮汐与悬浮泥沙浓度的潮变化相关项；T4为悬浮泥沙浓度与潮流变化相关项，即水位、垂向流速和垂向含沙量的三阶相关项，该项反映了底部泥沙再悬浮的强度；T3与T4之和为潮泵效应，表示一个潮周期内涨落潮造成的上层水体和底部泥沙之间发生泥沙输运；T5为垂向流速变化和悬浮泥沙浓度变化的相关项，垂向上净环流产生的贡献；T6、T7为时均量和潮汐振动引起的剪切扩散，T8为垂向潮振荡引起的剪切扩散，T6、T7与T8之和为剪切扩散项。

（2）Richardson数 (Rig)

利用Richardson数(Rig)描述水体的混合和层化程度，其表达式为：

式中，N2为浮力频率，M2为流速切强度，z’为垂向距离，ρ为海水密度，u为各层海水流速，g为重力加速度。当Rig＜0.25（或ln(Rig/0.25)＜0）时，水体不稳定，混合增强；当Rig＞0.25（或ln(Rig/0.25)＞0）时，水体稳定，混合被抑制[29]。

（3）盐度层化系数

利用盐度层化系数可以判断水体的稳定程度，层化系数与水体稳定程度呈正相关性，层化系数越大，水体越稳定，水体各层的物质交换越不频繁，反之则相反。盐度层化系数[30]表达式如下：

式中，Sd、Sb分别为底层、表层盐度，S¯为垂向平均盐度，当N＜0.01时，水体属于强混合型，盐度密度分层不影响悬浮泥沙的垂向扩散；0.01＜N＜1时，水体为缓混合型；N＞1时，水体为弱混合型，盐度密度分层会抑制悬浮泥沙垂向扩散。

（4）遥感影像反演悬浮泥沙浓度

考虑到大面站取样时间不同步对于分析研究区悬浮泥沙浓度平面分布特征可能存在一定影响，为了弥补此方面的缺陷，利用取样时间段内（2020年10月23日）的遥感影像进行反演，获得研究海域表层悬浮泥沙浓度。本文采用韩国海洋卫星中心提供的GOCI-TSS影像来反演悬浮泥沙浓度，GOCI-TSS是基于GOCI数据的二级产品，其中的TSS是利用GDPS软件基于Case-2算法得到的产品。

3 结果

3.1 海流特征

根据实测海流资料，将各个站位表层、中层、底层海流特征值列于表1，并绘制各站位垂向海流矢量图（图3）。各站位流向以NE-SW为主，海流运动形式为逆时针往复流，涨潮流主方向为SW，落潮流主方向为NE。总体上，3#站位流速最大，最大流速出现在涨潮期表层，为105.80 cm/s，5#站位流速最小，最小流速出现在落潮期底层，为11.30 cm/s。涨潮时，各站位最大流速为78.78～105.80 cm/s，落潮时，各站位最大流速为52.00 ～79.06 cm/s。在潮周期内，涨潮流速大于落潮流速，落潮历时更长，研究区各站位涨、落潮最大流速都出现在中潮位附近（图8）。垂向上，研究区各站位流速从表层到底层逐渐减小。

图3 各站位垂向潮流矢量Fig.3 Vertical power flow vector at each station

表1 各站位海流观测结果Table 1 Marine current observation results at each station

准调和分析结果表明研究区各测站潮流以M2分潮为主，潮流椭圆长轴比为0.07～0.3，小于0.5，说明潮流性质为规则半日潮流。研究区海域余流变化复杂，方向既有向海又有向陆，大体呈顺时针方向旋转，与前人研究结果一致[25]。离岸较近的1#、4#站位余流值明显大于离岸较远的2#、3#、5#站位，最大余流出现在1#站位表层为10.1 cm/s，对应流向为107.0°，最小余流出现在5#站位中层为2.5 cm/s，对应流向为 279.7°（表2）。

表2 各站位垂向余流特征值Table 2 Vertical residual power flow characteristic values at each station

3.2 悬浮泥沙时空分布及变化特征

3.2.1 悬浮泥沙浓度平面分布变化特征

研究区各海流观测站位涨潮阶段悬浮泥沙浓度稍微大于落潮阶段悬浮泥沙浓度（图4），涨潮阶段悬浮泥沙浓度为23.54～49.79 mg/L，落潮阶段悬浮泥沙浓度为20.75～46.09 mg/L。离岸较近的1#、4#站位平均悬浮泥沙浓度在涨潮、落潮阶段均明显大于离岸较远的2#、3#、5#站位。其中，平均悬浮泥沙浓度最大值出现在1#站位，其涨、落潮阶段平均悬浮泥沙浓度分别为 49.79、46.09 mg/L，其次是4#站位，其涨、落潮阶段平均悬浮泥沙浓度分别为 44.75、42.06 mg/L，最小值出现在3#站位，其涨、落潮阶段平均悬浮泥沙浓度分别为 23.40、20.17 mg/L。最大平均悬浮泥沙浓度出现在离岸较近的1#站位的涨潮阶段，为49.79 mg/L，最小平均悬浮泥沙浓度出现在离岸较远的3#站位的落潮阶段，为20.17 mg/L。

图4 平均悬浮泥沙浓度随涨落潮分布图Fig.4 Distribution of average suspended sediment concentration with fluctuating tide

根据大面站水文泥沙现场调查结果，研究区水体悬浮泥沙浓度呈现由陆向海逐渐降低的趋势（图5）。为了弥补大面站取样时间不同步对于分析研究区悬浮泥沙浓度平面分布特征存在的影响，本文利用取样时间段内的遥感影像进行反演，获得研究海域表层悬浮泥沙浓度（图6）。结果表明，遥感影像反演得到的悬浮泥沙浓度平面分布特征与实测结果大体相同，日照近海悬浮泥沙浓度呈现明显的近岸高、远岸低的分布特征。

图5 秋季大面站悬浮泥沙浓度平面分布图Fig.5 Plane distribution map of suspended sediment concentration at large surface station in autumn

图6 研究区表层悬浮泥沙浓度平面分布图Fig.6 Plane distribution of surface suspended sediment concentration in the study area

3.2.2 悬浮泥沙浓度垂向分布变化特征

垂向上，研究区悬浮泥沙浓度呈现从表层向底层逐渐升高的分布特征（图7）。根据各站位平均悬浮泥沙浓度垂向变化统计表（表3），4#站位的表、底层平均悬浮泥沙浓度差值最大，为29.06 mg/L，5#站位的表、底层平均悬浮泥沙浓度差值最小，为22.52 mg/L。研究区大面站调查结果表明，表层水体悬浮泥沙浓度为14.79～45.22 mg/L，平均值为23.97 mg/L；中层水体悬浮泥沙浓度为19.43～53.85 mg/L，平均值为30.22 mg/L；底层水体悬浮泥沙浓度为33.25～85.53 mg/L，平均值为50.19 mg/L。

表3 各海流观测站位平均悬浮泥沙浓度垂向变化统计Table 3 Statistical table of vertical variation of average suspended sediment concentration at each current observation station mg/L

图7 秋季典型断面悬沙浓度垂向分布图Fig.7 Vertical distribution of suspended sediment concentration at typical sections in autumn

3.2.3 悬浮泥沙浓度在潮周期内变化特征

研究区各站位悬浮泥沙浓度在潮周期内大小变化非常复杂（表4）。悬浮泥沙浓度变化周期与流速变化周期类似，各站位悬浮泥沙浓度在单日内出现4次峰值（图8）。由于当天海上风况良好，除1#站位外，其余站位涨潮阶段的悬浮泥沙浓度峰值都比落潮阶段更大。研究区海域悬浮泥沙浓度的变化滞后于流速的变化，从各站位流速和悬浮泥沙浓度随时间变化分布图（图8）中可以发现，悬浮泥沙浓度的峰值较流速峰值存在2 h的滞后。

图8 流速和悬浮泥沙浓度随时间变化分布图Fig.8 Distribution of velocity and suspended sediment concentration over time

表4 悬浮泥沙浓度特征值Table 4 List of characteristic values of suspended sediment concentration mg/L

3.3 温度、盐度特征

根据CTD实测数据，绘制 1#、3#、4#站位温度、盐度时间序列变化图（图9）。研究区3个站位温度为 18.66～21.28 ℃，盐度为 26.85～29.31 PSU，平均温度为20.31 ℃，平均盐度为28.39 PSU。3个站位观测期间表层平均温度为20.28 ℃，底层平均温度为20.34 ℃；表层平均盐度为28.34 PSU，底层平均盐度为28.42 PSU（表5）。由图9可以发现，一个潮周期内温度和盐度都存在两个峰值，涨潮时温度和盐度升高，落潮时温度和盐度降低，温度峰值对应盐度高点。对比研究区海水温度、盐度周期性波动曲线与潮位波动曲线，可以发现三者呈正相关，这表明当海况良好时，研究区水体温度和盐度的变化主要是由潮流驱动下的水体水平运动造成的[23]。

表5 1#、3#、4#站位温度、盐度特征值Table 5 The characteristic values of temperature and salinity at stations 1, 3 and 4

图9 1#、3#、4#站位温度、盐度随时间变化分布图Fig.9 Distribution of temperature and salinity at stations 1, 3 and 4 over time

研究区大面站调查结果表明，各层水体温度为18.30～21.23 ℃，表层、中层、底层温度平均值分别为 20.07、20.08、20.09 ℃；各层水体盐度为 27.38～29.70 PSU，表层、中层、底层盐度平均值分别为28.60、28.61、28.63 PSU。水体表层和底层的温度、盐度差值极小（表6），温度、盐度垂向梯度接近于零，平面上总体呈近岸低、远岸高的特征（图10）。

图10 大面站温度、盐度平面分布图Fig.10 Planar distribution of temperature and salinity of the large surface station

表6 大面站温度、盐度特征值Table 6 Table of characteristic values of temperature and salinity of the main surface station

3.4 悬浮泥沙输运通量

根据各站位悬浮泥沙输运项及单宽净输沙率计算结果（表7），比较各站位单宽净输沙率（T总）可以发现，离岸较近的4#站位单宽净输沙率（T总）最大，为 24.68 g/(s·m)，对应输沙方向为 5.61°；其次为离岸较近的 1#站位，单宽净输沙率（T总）为19.12 g/(s·m)，对应输沙方向为 103.02°；最小单宽净输沙率（T总）出现在离岸较远的 5#站位，为7.62 g/(s·m)，对应输沙方向为 247.69°。T1最大值出现在 4#站位，为 38.96 g/(s·m)，对应输沙方向为25.47°；T2最大值出现在 1#站位，为 9.82 g/(s·m)，对应输沙方向为214.97°；T3+T4最大值出现在4#站位，为 7.34 g/(s·m)，对应输沙方向为 248.82°；T5最大值出现在2#站位，为3.75 g/(s·m)，对应输沙方向为234.00°；剪切扩散项（T6+T7+T8）净输沙率较低，所有站位均低于 1.3 g/(s·m)。

4 讨论

4.1 悬浮泥沙输运机制

结合各站位悬浮泥沙输运项及单宽净输沙率计算结果（表7）和各悬浮泥沙输运项在单宽净输沙率中占比（表8），研究区域悬浮泥沙输运整体以平流输运项（T1+T2）为主，其次为潮泵效应（T3+T4）或垂向净环流贡献项（T5）。平流输运项（T1+T2）在5个站位单宽净输沙率中平均占比为95.24%，其输沙方向与净输沙方向基本一致；潮泵效应（T3+T4）在1#、4#站位中占据次要地位，在单宽净输沙率（T总）中分别占比27.77%、29.73%；垂向净环流贡献项（T5）在 2#、3#、5#站位中占据次要地位，在单宽净输沙率（T总）中分别占比27.98%、23.39%、25.86%。剪切扩散项（T6+T7+T8）净输沙率较低，在单宽净输沙率（T总）中平均占比为6.33%，对研究区域悬浮泥沙输运贡献有限。

表7 悬浮泥沙输运项及单宽净输沙率Table 7 Suspended sediment transport items and net sediment transport rate per width

表8 各悬浮泥沙输运项在单宽净输沙率中占比Table 8 Proportion of suspended sediment transport terms in net sediment transport rate per width %

研究区单宽输沙量变化复杂，1#、2#站位单宽净输沙量向海，而3#、4#、5#站位的单宽净输沙量向陆。对比研究海域各站位垂向平均余流和悬浮泥沙输运项可以发现，垂向平均余流方向与拉格朗日平流输运项（T1+T2）方向和单宽净输沙项（T总）方向基本相同，所有站位余流方向与T1+T2、T总的输沙方向相差均不超过25°（图11）。可以初步判定，潮余流是主导悬浮泥沙输运的重要因素。比较各站位单宽净输沙率（T总）可以发现，离岸较近站位的单宽净输沙率（T总）明显大于离岸较远站位的单宽净输沙率（T总）。这是由于平流输运项（T1+T2）方向受余流控制，量值由流速和悬浮泥沙浓度控制，由于近岸1#、4#站位余流流速和悬浮泥沙浓度数值均高于远岸2#、3#、5#站位，所以近岸站位的拉格朗日平流输运项（T1+T2）数值更大，进而导致离岸较近站位的单宽净输沙率（T总）明显大于离岸较远站位的单宽净输沙率（T总）。

图11 悬浮泥沙输运通量与余流矢量叠置图注：1#、4#站位悬沙输运项缩小至原来的1/2，5#站位悬沙输运项扩大至原来的2倍。Fig.11 Overlay diagram of suspended sediment transport volume and residual flow vector

通过研究区各站位涨落潮时悬浮泥沙输运通量表（表9）可以发现，该海域的悬浮泥沙涨潮时沿SW向输运，落潮时沿NE向输运，与涨落潮的方向基本一致。比较涨落潮输沙率和净输沙率可以发现，研究区净输沙率远小于涨落潮输沙率，这表明研究区水体中的悬浮泥沙只是在涨落潮流携带下做往复搬运，实际产生的净输沙量很小。但潮余流是非周期性的，可以携带泥沙做远距离运移。研究区在潮周期内，涨潮流（SW向）占优势。2#、5#站位涨潮流方向和潮余流方向大体相同，受其影响2#、5#站位呈向南的净输沙趋势；与之不同的是，3#站位由于其周围及其西部海域存在一狭长地带，水深较深，受水深、地形影响，3#站位呈向西的净输沙趋势；1#、4#站位受到近岸岬角、防波堤、岚山港池以及航道等阻挡，呈现向东、向北的净输沙趋势。整体而言，研究区悬浮泥沙在远岸开阔海域呈向南的净输运趋势。

表9 各站位涨落潮悬浮泥沙输运通量Table 9 Suspended sediment transport at each station at ebb and flow tide

南黄海西部近岸海域影响悬浮泥沙输运的因素十分复杂，前人的成果表明废黄河口海域可能同时存在向北（海洲湾方向）和向南（苏北辐射沙洲方向）的输运趋势[31]。针对日照南侧废黄河口海域悬浮泥沙的输运观测，鲁号号等[32]认为悬浮泥沙输运方向和潮余流方向基本相同，涨潮期间悬浮泥沙沿岸向南输运(涨潮流方向)，落潮期间悬浮泥沙沿岸向北输运(落潮流方向)；周良勇等[33]认为悬浮泥沙涨潮时向南运移，落潮时向北运移，与涨落潮的方向基本相同；杨林等[34]认为悬浮泥沙输运的主要影响因素是T1和T2项，各测站净输沙方向与潮余流方向较为一致。徐粲等[35]对苏北辐射沙脊群的研究表明，该区域T1和T4项对研究区的悬沙输运贡献最大。

综合前人的成果并结合本文的研究，可以推断南黄海西部近岸海域悬浮泥沙在涨潮、落潮期间输运方向与涨潮、落潮方向基本一致，悬浮泥沙净输运方向和潮余流方向大体相同，受余流、悬浮泥沙浓度、水深地形等多因子影响，在不同潮周期内各输沙项对净输沙率的贡献不尽相同。

4.2 悬浮泥沙输运的影响因素

影响悬浮泥沙输运的水文动力因素包括海洋环流、潮汐、潮流和波浪等，在外陆架和中陆架地区，悬浮泥沙的长期输运受环流控制[36]，黄海陆架深海区悬浮泥沙输运受苏北沿岸流、黄海西部沿岸流、黄海暖流、黄海冷水团等影响，存在“夏储冬输”的季节性输运模式；在内陆架浅海区，潮流和波浪对悬浮泥沙的分布、输运起决定性作用，其中潮流引起悬浮泥沙浓度的潮周期变化，波浪控制悬浮泥沙浓度季节性变化[37-38]。前人对南黄海近海悬浮体无机物和有机物组分含量的研究证实，该区域的悬浮泥沙主要为陆源物质，其主要来源为近岸地区海底沉积物的再悬浮和河流输沙[39-41]。日照及海州湾海域周边沿岸无大型供沙河流，河流输沙对该海域悬浮泥沙分布及输运影响有限[24，42]。研究区岚山港和日照港附近海域表层沉积物以含砾泥、砂质泥和泥等细粒物质为主[24]，细粒物质被侵蚀、搬运、再悬浮之后，在波浪、潮流和水体状态等因素的控制下发生输运和再沉积。

4.2.1 水体垂向混合对悬浮泥沙垂向运移的影响

在近岸浅海区，水体的混合和层化程度会影响悬浮泥沙垂向分布[21]。水体垂向混合较弱时会出现分层从而抑制悬浮泥沙向上的运移扩散。根据1#、3#、4#站位的实测数据，计算Richardson数来判断水体混合和层化程度，并分析其对研究区悬浮泥沙浓度垂向分布的影响。当Rig＜0.25（或ln(Rig/0.25)＜0）时，水体不稳定，混合增强；当Rig＞0.25（或ln(Rig/0.25)＞0）时，水体稳定，混合被抑制。计算结果显示，研究区水体在观测期间Richardson数大部分时段小于0.25，表明水体混合良好，仅在少部分时段出现微弱分层现象（图12）。

图12 各站位垂向Richardson数时间序列图Fig.12 Vertical Richardson number time series diagram of each station

在近岸海域，由盐度产生的水体密度分层会抑制悬浮泥沙的垂向扩散，进而影响悬浮泥沙浓度垂向分布[43]。盐度层化系数可以判断水体的稳定程度，层化系数越大，水体越稳定，水体各层的物质交换越不频繁，反之则相反。为了探究本次观测期内水体密度分层对悬浮泥沙浓度垂向分布的影响，本文利用实测数据计算了研究区3个海流观测站位和35个大面站的盐度分层系数。数据计算结果表明，观测期间所有站位的盐度分层系数均小于0.005，表明水体处于强混合状态，盐度密度分层极小，不影响悬浮泥沙浓度垂向扩散（图13）。

图13 大面站盐度分层系数平面分布图Fig.13 Planar distribution of salinity stratification coefficient of Dameen Station

整体上，研究区水体垂向混合良好，对底层悬浮泥沙的向上扩散运移无抑制作用，推断可能受大潮期间强潮流作用影响，潮流的混合搅拌作用增强，导致海水垂向混合作用强烈。

4.2.2 潮流在潮周期内对悬浮泥沙分布、输运的影响

再悬浮作用是南黄海近岸海域悬浮泥沙浓度变化的主导因素[44]，研究区海域悬浮泥沙浓度在潮周期内的变化与潮流流速的变化存在高度相关性，悬浮泥沙浓度的高值区在流速高值区时刻附近，悬浮泥沙浓度的低值区也对应着流速低值区附近（图8）。以往研究表明，同样的温盐条件下，潮流强度对悬浮泥沙的垂向分布有较大影响[45]，潮流强度越高，悬浮泥沙的垂向扩散作用和再悬浮作用就越强，反之则弱，同时潮流的加速和减速也深刻影响着悬浮泥沙垂向分布。

在涨、落潮过程中，初期随着潮流流速的增大，海底的表层沉积物受到侵蚀和冲刷，当底层潮流流速大于泥沙起动流速时，大量海底泥沙会被起动，悬浮于水中，导致近底层悬浮泥沙浓度快速增大，随后底层悬浮泥沙在剧烈的紊动作用下向上运移、扩散，致使中上层悬浮泥沙浓度增大；在涨潮落潮相互转换时期，由于流向的转变，潮流流速较低，再悬浮作用会变弱，悬浮泥沙会发生沉降，导致悬浮泥沙浓度变低。总体而言，由于涨潮流速略大于落潮流速，所以涨潮阶段的平均悬浮泥沙浓度略大于落潮阶段。然而，海底表层沉积物的起动、再悬浮均需要一定时间，因此悬浮泥沙浓度的变化滞后于流速的变化，研究区海域悬浮泥沙浓度的峰值较流速峰值存在2 h的滞后。由于研究区海域潮流为典型的往复流，悬浮泥沙在涨、落潮流挟带下沿涨、落潮方向进行物质运移和交换，涨潮时沿涨潮流方向（SW）净输运，落潮时沿落潮流方向（NE）净输运。

根据宋红瑛等[24]人对日照近海表层沉积物粒度特征的研究，利用窦国仁泥沙起动公式，大致计算出研究区潮流流速达到18 cm/s时便可起动粒径小于0.16 mm（大于2.6 Φ）的细粒物质，日照近海表层沉积物以细砂（2～4 Φ）为主，海流调查站位附近沉积物平均粒径为2～5 Φ，近岸1#、4#站位沉积物粒径较远岸2#、3#、5#站位更小[24]。研究区大潮期潮流流速为11.30～105.80 cm/s，因此绝大部分海域海底泥沙都能在潮流作用下发生起动（图14），NESW向潮流既能输沙又能掀沙。表层沉积物粒径大小是影响沉积物再悬浮的重要因素，细粒物质沉降速度慢，发生再悬浮时需要的水动力条件更弱，相反粗粒沉积物更难于向水体上部运移、扩散[46]。相较于远岸的2#、3#、5#站位，岚山港附近的1#、4#站位沉积物粒径更小，起动所需要的潮流流速更小，且1#、4#站位水动力条件较2#、3#、5#站位更强，潮流施加于海床沉积物的底切应力作用更强，使得1#、4#站位沉积物更易发生起动和再悬浮。因此，受底质类型和潮流流速的共同影响，1#、4#站位悬浮泥沙浓度明显高于2#、3#、5#站位。