黄海是一个半封闭式的陆架浅海^[1]，是世界上悬浮泥沙含量最高的海域之一，每年有大量来自中国和朝鲜半岛河流携带的陆源泥沙汇入黄海^[2]。山东半岛东部海域位于南、北黄海的交界区域，水动力作用较为强烈，悬浮泥沙含量高，是现代黄河沉积物向外海扩散的重要通道^[3]。该海域既是陆源沉积物的汇，也是外海沉积物的源^[4]，并且海域内沉积有记录环境变化信息的泥质体，因而引起了国内外学者的广泛关注。

泥质体主要发育在大陆架边缘，不仅记录着不同历史时期的海侵海退和气候变化^[5]等信息，同时也记录了人类活动的信号^[6]。二十世纪八十年代发现在山东半岛东部海区存在一向海延伸的泥质楔^[7]，山东半岛泥质楔保存着全新世以来黄河泥沙及入海口的变化，也保存了山东半岛东部海域的海流、气候、环境变化等信息^[8]。根据高分辨率浅地层剖面资料^[9-10]、同位素测年^{[9, 11-12]}、矿物和地球化学分析^[13-14]等多种分析手段，认为泥楔形成于11.6~4 kaBP^[8-10]，其表层沉积速率为6~12 mm/a^[15]。

黄河入海泥沙大部分在山东半岛水下三角洲区域沉积，剩余约有1%~30%的黄河泥沙通过渤海海峡向北黄海输运，而这其中有一部分的泥沙沉积在山东半岛北岸，其余的泥沙受沿岸流的影响绕过成山头进入山东半岛泥质区，最终可输运至黄海海域^{[10, 15-18]}。泥沙在渤黄海的扩散和沉积具有明显的季节性变化^[19]，泥沙多在夏季沉积入海，而冬季是泥沙向深海输运的主要季节，由于受冬季强东北风的影响，山东半岛沿岸流增强、温跃层消失^[20]，沿岸流会携带大量泥沙沿岸向外海输运，整个海域沉积物输运体系呈现典型的“夏储冬输”特征^[21]；但由于山东半岛沿岸流和黄海暖流形成的近岸强海流切变锋的阻碍，可抑制悬浮泥沙运移至山东半岛东部陆架末梢，有利于山东半岛Ω状泥质沉积体的形成^[22]。

众多学者采用大面站观测、遥感反演和数值模拟等方法，研究了山东半岛东部海域悬浮泥沙输运特征，但该区域多站位同步连续潮周期定点观测的研究成果较少，而且利用通量机制分解方法进行本海域悬浮泥沙输运动力机制的研究也很少。本文将在前人研究的基础上，根据4个站位25 h海流、悬沙同步观测和取样分析，运用通量机制分解方法，研究山东半岛东部海域冬季悬浮泥沙浓度时空变化特征，探讨山东半岛东部海域冬季悬浮泥沙的输运机制。

1 研究区概况

研究区位于山东半岛东部海域地处中纬度，属于温带季风气候，冬季主要是受山东半岛沿岸流、黄海暖流的影响^[23]。该区域位于南黄海与北黄海过渡区域，水深在20~40 m之间，海底凹凸不平，研究站位西侧与山东半岛之间有一道近NW向较深的海底冲刷潮槽，发育一环绕山东半岛的反C形楔状泥质沉积体^[8]，整体呈现中间厚、向海和陆减薄的Ω形状^[10]，最大厚度可达40 m，约在123°10′E处尖灭^[24]，最远可以延伸至80 m等深线处^[25]。

2 资料与方法 2.1 资料来源

2015年12月12~13日，中国海洋大学在山东半岛东部海域泥质区进行了4个站位的大潮期25 h同步连续海流观测和悬浮泥沙取样，观测站位见图 1。

采用声学多普勒流速剖面仪(ADCP，600 KHz，美国)及直读式电磁海流计(JFE AEM213-D，日本)观测各层海流的流向和流速数据；采用CTD(美国TRDI公司CTD-NV型CTD)获取剖面温度、盐度数据，每隔1 h采集一组剖面数据；利用采水器现场同步采集水样；采用蚌式抓斗获取底质样品；流速、流向和悬浮泥沙按表层、中层、底层三层进行数据处理分析。表层沉积物及悬浮泥沙粒度采用Mastersizer 2 000激光粒度仪测定，分析范围为0.02~2 000 μm，样品重复测量相对误差小于3%。悬浮泥沙浓度采用抽滤实验测定，测量准确度约2%，详细步骤如下：

(1) 将双层微孔滤膜(直径为47 mm，孔径为45 μm)在40 ℃条件下烘干至恒重后，在10^-4 g/g分析天平上称滤膜重量，并将其置于已编号的空滤膜盒中。

(2) 将取得的水样经过滤膜抽滤，抽干后记录过滤海水体积，用去离子水冲洗滤膜2次。

(3) 将带有悬浮泥沙样品滤膜烘干至恒重，并称量样品和滤膜总重量。

(4) 将20%的样品加空白校正膜平行进行空白校正。

2.2 计算方法 2.2.1 再悬浮、沉降通量

研究区4个站位悬浮泥沙粒度粒径全部小于0.062 5 mm(4Φ)，属于粘性泥沙，故再悬浮通量计算采用Partheniades^[26]提出的公式，沉降通量采用Krone^[27]提出的公式：

$\frac{{{\rm{d}}{\mathit{M}_\mathit{E}}}}{{{\rm{d}}t}} = E\left( {\frac{{{\tau _0}}}{{{\tau _{cr}}}} - 1} \right), $

(1)

$\frac{{{\rm{d}}{\mathit{M}_D}}}{{{\rm{d}}\mathit{t}}} = {C_b}{W_s}\left( {1 - \frac{{{\tau _0}}}{{{\tau _{cr}}}}} \right)。$

(2)

式中：E是再悬浮常数；τ₀为底切应力；τ_cr为临界切应力；c_b为近底部悬浮泥沙浓度；w_s为悬浮泥沙的沉降速度。

2.2.2 理查森数

在分层的海洋中，理查森数(Ri_g)用来判断动能和密度层化的相对重要性：

$R{\mathit{i}_g} = {N^2}/{M^2}, $

(3)

${N^2} = - \left( {\mathit{g/\rho }} \right)\left( {{\rm{d}}\mathit{\rho /}{\rm{d}}\mathit{z}} \right)。$

(4)

其中N²为浮力频率，反映了密度层化时，流体质点由于浮力而产生垂直振荡的频率。Ri_g < 0.25时，流体质点的速度的作用大于密度层化的抑制，流体会出现剪切不稳定，海洋中产生湍流；Ri_g>0.25时，水体层化现象明显。

2.2.3 通量机制分解

根据Dyer^[28]的物质通量计算方法将悬浮泥沙通量分解为多个动力项，从而探讨各动力项所对应的主要动力因子对研究区悬浮泥沙输运的相对贡献大小。瞬时流速可以分解为垂向平均量及其偏差量，和又可分解为潮平均量和潮偏差量，则；同样含沙量c可以分解为；水深可表示为h=h₀+h_t，则潮周期T_t平均瞬时单宽悬浮泥沙输移通量T的计算式为：

$\begin{array}{l} T = \frac{1}{T}\int\limits_0^{{T_t}} {\int\limits_0^h {uc\;\;{\rm{d}}{\mathit{z}_0}{\rm{d}}\mathit{t}} } = \frac{1}{T}\int\limits_0^{{T_t}} {\int\limits_0^1 {uch\;\;{\rm{d}}\mathit{z}{\rm{d}}\mathit{t}} } = \\ {h_0}\overline {{u_0}} \overline {{c_0}} + \left\langle {{h_t}\overline {{u_t}} } \right\rangle \overline {{c_0}} + \left\langle {{h_t}\overline {{c_t}} } \right\rangle \overline {{u_0}} + \left\langle {{h_t}\overline {{u_t}} \overline {{c_t}} } \right\rangle + \\ {h_0}\overline {u{\mathit{'}_0}c{\mathit{'}_0}} + \left\langle {{h_t}u{\mathit{'}_0}c{\mathit{'}_t}} \right\rangle + \left\langle {{h_t}u{\mathit{'}_t}c{\mathit{'}_0}} \right\rangle + \left\langle {{h_t}u{\mathit{'}_t}c{\mathit{'}_t}} \right\rangle = \\ \left( {{T_1}} \right) + \left( {{T_2}} \right) + \left( {{T_3}} \right) + \left( {{T_4}} \right) + \left( {{T_5}} \right) + \\ {\rm{\;\;\;\;\;\;\;\;\;\;}}\left( {{T_6}} \right) + \left( {{T_7}} \right) + \left( {{T_8}} \right)。\end{array} $

(5)

式中：T₁为欧拉余流引起的悬浮泥沙输运；T₂为斯托克斯漂移输运量，T₁+T₂为拉格朗日平流输运贡献项；T₃为潮汐与悬浮泥沙浓度的潮变化相关项；T₄为悬浮泥沙与潮流场变化相关项，T₃、T₄是一个潮周期内涨落潮引起的水体和底部泥沙双向交换而产生的泥沙输运，称为潮泵效应；T₅为垂向流速变化和悬浮泥沙浓度变化的相关项，垂向上净环流产生的贡献；T₆，T₇为时均量和潮汐震动引起的剪切扩散；T₈为垂向潮振荡引起的剪切扩散。

3 结果 3.1 温度、盐度特征

根据温度、盐度调查结果，绘制1、2、4号站位温度、盐度垂向变化图(见图 2)和时间序列剖面图(见图 3)。

研究区3个站位温度、盐度垂向变化较小，混合较为均匀，水平梯度也较小，温度变化范围在8.78~10.30 ℃之间，盐度变化范围在31.21~31.60之间，整体呈低温低盐特征，主要受低温低盐的山东半岛沿岸流影响。

一个潮周期内，温度、盐度变化量较小且呈明显的周期性变化规律。1、2号站位涨潮时温度升高、盐度减小，落潮时温度降低、盐度增大；4号站变化规律正好相反，涨潮时温度降低、盐度增大，落潮时温度升高、盐度减小。

3.2 海流特征

海流调和分析表明，研究区潮流性质为正规半日潮，涨潮流速略大于落潮流速(见图 4)；潮流运动形式以往复流为主，表、中、底层流向基本一致(见图 5)，涨落潮方向随山东半岛沿岸走向变化而变化，1号站流向主要为WSW~ENE向，2、3号站主流向为SSW~NNE向，山东半岛北部东向海流流经成山头时，由于惯性作用，4号站表层海流流向呈SSE向。

涨落潮最大流速及平均流速波动范围较大，涨潮流时，研究区各站位潮流最大流速介于40.67~106.00 cm/s之间(见表 1)；落潮流时，最大流速介于42.67~83.33 cm/s之间。总体上，离岸最近的2号站流速最大，距岸最远的3号站流速整体较小；垂向上，各站位表层最大流速和平均流速均大于底层，流速自表层至底层逐渐减小。

3.3 悬浮泥沙浓度特征

研究区涨潮段悬浮泥沙浓度基本大于落潮段，但两者相差并不大。涨潮段含沙量介于15.4~64.8 mg/L之间(见表 2)，落潮段含沙量介于11.6~88.0 mg/L之间。平面上，悬浮泥沙浓度规律与流速相似，流速最大的2号站悬沙浓度也相对最大，3号站仍最小。垂向上，除1号站落潮段外，各站位平均含沙量自表层向底层递增，并且3号站和4号站垂向浓度变化较小，混合较好。

研究区海域潮周期内涨、落潮悬沙浓度大小变化过程非常复杂，悬沙浓度变化与流速的变化存在一定的相关性，含沙量对流速响应较好，流速高值区对应着悬浮泥沙浓度的峰值，但有轻微的滞后现象(见图 4)。

3.4 悬浮泥沙及海底表层沉积物粒度特征

根据谢帕德分类方法，研究区各站位、各层悬浮泥沙类型均为粘土质粉砂(见表 3)，粉砂、粘土含量高，平均值介于91.58%~97.77%之间；砂含量相对较低，平均值介于2.03%~8.42%之间。悬浮泥沙粒径较小，平均粒径介于6.66~6.87 Ф之间，各站位粒径相差不大，相对来说，3号站粒径较小，4号站位粒径最大；分选系数介于1.60~1.99之间，分选较差；偏态介于-1.51~1.10之间，从极负偏到正偏；峰态介于2.06~2.86之间，呈宽平分布。

4个站位海底表层沉积物粉砂含量最高，平均值介于65.65%~70.81%之间；1、2号站为砂质粉砂，3、4号站为粘土质粉砂。平均粒径较悬浮泥沙要大，介于5.81~6.30之间，3号站粒径仍最小，2号站粒径最大；分选系数介于2.23~2.31之间，分选差；偏态介于2.56~2.65之间，极正偏；峰态介于3.40~3.51之间，呈很宽峰态。

3.5 余流特征

余流是指从实测海流中分离出周期性潮流后所余下的非周期的流。通过准调和分析得到的欧拉余流，其表征经过同一点流体微团的平均速度；斯托克斯漂流表示的是水体的净漂移量；拉格朗日余流则是欧拉余流与斯托克斯漂流之和，描述了水体中悬浮物潮周期的净输移速度，能较好指示悬浮泥沙的输运。各站位欧拉余流、斯托克斯漂流及拉格朗日余流的计算结果见表 4。

欧拉余流在拉格朗日余流中占主导地位。2、3、4号站拉格朗日余流数值相近，并且三个站位的斯托克斯漂流较弱，欧拉余流与拉格朗日余流相差较小；1号站斯托克斯漂流是4个站位中最强的，对拉格朗日余流影响相对较大，但其欧拉余流较其他三个站位小得多，故1号站拉格朗日余流最弱。

4 讨论 4.1 悬浮泥沙垂向运移特征 4.1.1 悬浮泥沙和表层沉积物的交换

涨停、落停时悬浮泥沙粒径较小(见图 6)，但涨急、落急时粒径相对较大，频率曲线向粗颗粒方向(φ值减小)偏移，表层沉积物平均粒径较悬浮泥沙要大，说明涨急、落急时刻有底部颗粒沉积物加入水体，发生了再悬浮。

将实测沉积物粒度、近底悬浮泥沙浓度、流速数据带入式(1)、(2)计算表明，潮周期内各站位再悬浮通量、沉降通量值较小(见图 7)，但再悬浮通量略大于沉降通量，且再悬浮通量变化与潮流流速有关。研究区海底表层沉积物主要为粘土质粉砂和粉砂，颗粒较细，在研究区较强水动力作用下，泥沙较易起动，流速增大，再悬浮通量随之增大，反之，再悬浮通量逐渐减小；当流速降至最低时，因底部切应力小于临界值，悬浮泥沙发生沉降。其中，2号站虽表层沉积物粒度最粗，但其流速最大，近底剪切应力最大使得再悬浮通量值较高；而1号站表层沉积物粒度较粗，流速较小，因而再悬浮通量值最小、沉降通量值最大。

悬浮泥沙和底沙的交换直接影响底层悬浮泥沙浓度的变化，进而影响研究区含沙量的变化。但是再悬浮通量峰值与含沙量浓度峰(见图 4)并不完全一致，说明研究区悬浮泥沙浓度不仅会受到再悬浮的影响，在再悬浮通量峰值发生后，还受平流输沙等其他因素影响，悬浮泥沙浓度仍会继续升高。

4.1.2 悬浮泥沙的垂向混合

研究区海域悬浮泥沙受湍流作用垂向混合剧烈，只有个别时段理查森数(Ri_g)略大于0.25，水体有轻微的层化现象(见图 8)。1、2号站Ri_g略高于3、4号站位，但研究区4个站位Ri_g大多都小于0.25，基本无层化现象，水体中悬浮泥沙在强湍流混合作用下扩散到整个水柱中，并且研究区4个站位表、中、底各层悬浮泥沙粒径相差较小(见表 3)，底层和表层的平均粒径平均仅差0.08 Ф，表底悬浮泥沙混合均匀；只在1号站19点、2号站18点，Ri_g大于0.25，此时水体处于层化状态，湍流被抑制，悬浮泥沙的垂向混合减弱，致使底层再悬浮的泥沙无法进入上层水体，底层含沙量明显大于中、上层(见图 4)。

4.2 悬浮泥沙输运机制

根据悬浮泥沙输运通量机制分解结果(见表 5)，拉格朗日平流输移项(T₁+T₂)占据主导地位，其次为垂向净环流(T₅)和潮泵效应(T₃+T₄)，剪切扩散项(T₆+T₇+T₈)数值极小。拉格朗日平流输移项(T₁+T₂)不仅取决于拉格朗日余流的大小，还与悬浮泥沙浓度有关，1号站余流最小、3号站悬浮泥沙浓度最小，故1、3号站T₁+T₂较小。四个站位欧拉余流输运项(T₁)对净输沙的贡献最大；2、3、4号斯托克斯余流输运量(T₂)较T₁小得多，对整体悬浮泥沙输运影响较小，而1号站斯托克斯余流相对较大，对平流输运及整体净输沙影响较大；T₃+T₄主要取决于泥沙再悬浮和沉降作用，T₃、T₄、T₅对悬浮泥沙输运也有影响，但是数值都不大，影响十分小；T₆+T₇+T₈贡献极小，对研究区整体悬浮泥沙输运几乎无影响。

根据上述计算可知，研究区涨落潮水体的平流输运对泥沙运移起主要作用，而再悬浮、垂向交换是相对次要的。这说明研究区涨落潮流虽然会携带搬运大量悬浮泥沙，但其中大部分泥沙只是随潮流作用往复搬运，并未产生净输运；泥沙再悬浮、沉降与湍流也都会引起研究区悬浮泥沙浓度变化，但对悬浮泥沙净输运的影响较小。而余流是非周期性的，可以携带泥沙进行长距离搬运，研究区4个站位净输沙方向与拉格朗日余流方向(见图 9)大致相同，平流输运才是影响研究区海域悬浮泥沙运移的最主要因素。

悬浮泥沙净输运主要向泥质沉积体中心或沿岸输运。北侧4号站悬浮泥沙主要向东南方向输运(见图 9)，3号站向正南方向运向泥质沉积体的中心(见图 1)，西侧2号站和南侧1号站分别沿岸向西南和西北方向输沙。3号站与2、4号站的悬浮泥沙输运通量值差距大，而2、3、4号站余流大小相差较小，这是因为3号站的悬浮泥沙浓度明显小于2、4号站。3号站位于泥质沉积体顶积层之上，顶积层附近海域存在温度锋面^[29](见图 10)，起到“水障”的作用，故此处流速较低，悬浮泥沙浓度也相应较小。

黄海是一半封闭的陆架浅海，几乎不受外海物质的影响^[30]。山东半岛陆域沉积物通过黄河及沿岸短源河流搬运作用入海，山东半岛沿岸流再携带入海悬浮泥沙及侵蚀冲刷黄河三角洲、滨岸带再悬浮的泥沙通过渤海海峡，搬运出渤海^[31]，之后沿山东半岛沿岸向东输运，绕过成山头后继续沿岸向南运移到山东半岛东部海域。在强冬季风作用下，沿岸流增强，山东半岛东部海域会发生再悬浮现象，并且垂向混合剧烈。南下的沿岸流与北上的黄海暖流在山东半岛东部海域发生强海流切变^[1]，形成的锋面抑制悬浮泥沙的离岸输运，使悬浮泥沙只能扩散至泥质体顶积层附近，或者沿岸向山东半岛南侧海域输运，并不能穿越锋面继续向外海输运^[1]，因此山东半岛东部海域冬季悬浮泥沙浓度高。而夏季悬浮泥沙不易再悬浮，且温跃层也阻止底层悬浮泥沙的垂向扩散，有助于平流输运携带的泥沙及之前水体中的泥沙在此沉积，堆积形成山东半岛东部海域特有的Ω形泥质沉积体。

5 结论

(1) 研究区潮流性质为正规半日潮, 潮流运动形式为往复流，各站位涨落潮方向大致与山东半岛海岸线平行，潮流流速自表层至底层、由岸向海逐渐减小，涨潮流速略大于落潮流速。

(2) 悬浮泥沙浓度涨潮时段大于落潮时段，且自表层至底层逐渐增大，但变化梯度较小，垂向混合较好。悬浮泥沙浓度变化与潮流流速大致呈正相关关系，但相对于潮流流速存在轻微的滞后现象。

(3) 悬浮泥沙类型均为粘土质粉砂，粒径较小，分选较差；海底表层沉积物平均粒径大于悬浮泥沙，距岸近的站位为砂质粉砂，离岸相对较远的站位为粘土质粉砂，分选差。

(4) 悬浮泥沙浓度变化一定程度上受泥沙再悬浮、沉降与湍流的影响。研究区各站位的再悬浮通量略大于沉降通量，但量值均较小；悬浮泥沙受湍流作用垂向混合剧烈(理查森数 < 0.25)，个别时段有轻微的层化现象。

(5) 悬浮泥沙通量机制分解表明，平流输运是山东半岛东部海域悬浮泥沙净输运中最主要的输沙机制、占据主导地位，悬浮泥沙净输运方向主要向泥质沉积体中心或沿岸输运。