渤海湾泥质区南部悬浮体夏季输运过程与机制

王燕，陈斌，李日辉，徐晓达，苏大鹏

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

2. 青岛海洋科学与技术试点国家实验室海洋矿产资源评价与探测技术功能实验室，青岛 266237

渤海湾毗邻现代黄河三角洲，沿岸主要为砂质淤泥质海岸，其沉积物主要来自沿岸河流携带的大量泥沙，受河流与沿岸细粒物质大量供给和潮流作用为主的海洋动力控制，呈不规则带状和斑块状分布，沿岸粒度较粗，多为粉砂和黏土粉砂，东北部沿岸多为砂质粉砂，中部海域粒度较细，多为黏土软泥和粉砂质软泥[1]。渤海湾沿岸有多条河流，如黄河、海河、滦河等，其中黄河和滦河是渤海湾沉积物的主要来源[2]，并以黄河入海泥沙量最大[3]。黄河每年入海泥沙约11亿t，大量入海泥沙在黄河口沉积下来，不仅使河流入海流路发生频繁变迁，而且导致河口三角洲快速淤积[4-7]。在海洋动力作用下，大量细颗粒泥沙自黄河口三角洲向北、东北、西北方向海域传输，可达到渤海湾湾顶[8-9]。自1976年以来其北部废弃的三角洲叶瓣遭受强烈的海岸侵蚀，在冬季强海洋动力作用下，大量再悬浮沉积物为渤海湾区域提供充足的物源供应[3,10-11]。随着人类活动的增强，河流入海径流量和泥沙量明显减少，粒度粗化[6,10,12]，大量泥沙快速沉积在河口附近，在波浪和潮流的控制下，泥沙再悬浮并在渤海海域传输，在沿岸悬浮体浓度明显高于渤海湾中部[13-14]，水体底层浓度明显高于表层浓度，悬浮体浓度的分布呈季节性分布特征[3,11,15-16]，冬季渤海湾在强风浪的再悬浮和沿岸流的搬运作用下，悬浮体浓度分布明显高于海洋动力较弱的夏季[3,17-18]，水体的层化混合过程对黄河入海悬浮体传输过程起到重要的控制作用[6,19]。除此之外，渤海湾海域悬浮体的分布还具有明显的大小潮差异，大潮时期悬浮体浓度高于小潮时期[20]。因此，不同时段、不同区域内，渤海湾海域悬浮体的分布特征与海域内海洋动力条件有着密不可分的关系，要了解该海域悬浮体分布的变化规律，就需要分不同时段对不同区域进行连续观测研究，才能更全面地掌握渤海湾悬浮体的运移特征。但是迄今为止对黄河三角洲北部和渤海湾泥质区毗邻海域再悬浮沉积物的输运过程及其沉积分布尚缺乏系统性观测研究，这方面的观测有待加强。大量的研究表明[21-23]，近海沉积物输运过程及其沉积分布与污染物和生源要素的扩散具有极为密切的关系，是影响区域生物地球化学循环和生态系统的重要因素。渤海湾周边城市人口稠密，大中城市密集，有重要能源生产基地，是目前我国重要的经济发展区之一。随着大规模的围海造地工程和周边临港工程建设，填海面积增加的同时，大量来自陆域的重金属和污染物较易附着在悬浮体上随径流入海，并在海洋动力的作用下，在近岸海域扩散和沉积，对近岸海洋生态系统产生了重要影响[24-27]。因此渤海湾及毗邻海域不同时间尺度的沉积动力过程及其环境效应，有利于推进对该海域沉积动力学和“源–汇”过程研究的进一步认识，对区域环境可持续发展具有重要的意义。

本文在前人的研究基础上，对渤海湾泥质区南部靠近黄河三角洲北部海域的夏季悬浮体浓度、水体温盐结构进行了25 h连续观测调查，对悬浮体输运过程展开了时空对比分析，并探讨了其水动力机制。

1 材料及方法

1.1 数据采集

2012年8月20日—9月10日在渤海湾海域完成了对10个站位的海流连续25 h定点观测，并每隔1 h对垂向剖面进行温度、盐度和浊度的观测，采集表、中、底3层水位的悬浮体水样用于悬浮体浓度的实验室抽滤测试分析（图1）。其中，流速、流向及水深数据采集使用美国SonTek公司生产的声学多普勒流速剖面仪ADCP（工作频率为300 kHz，盲区为1 m，测流误差＜5 mm/s和＜1°）来进行观测，测层间隔设置为1 m，定点观测采样频率设为10 s，温度、盐度及浊度数据采集使用美国Seabird公司生产的SBE19 Plus多参数水质剖面仪来进行观测（温度、电导率和浊度的精度分别为＜0.1°C、0.001 mS/cm和0.01 NTU），用绞车以平均每2 s下放1 m的速度下放仪器，平均每0.2 s可获取一个数据，悬浮体水样采集使用美国Seabird公司生产的SBE ECO55自动采水器进行采集，该仪器与SBE19Plus集成在一起。悬浮体水样采集以ADCP测深为准，分层进行采取（表层样取自水深2 m以内，中层样取自0.5倍水深处，底层样取距底0.5～2 m处），与温度、盐度及浊度数据采集同步进行。本文采用本航次调查的10个站位中位于渤海湾泥质区南部的5个站位（A1—A5站）的观测结果进行悬浮体输运过程分析。

图 1 研究区域及定点连续观测站位泥质区位置出自文献[28]，遥感底图来自https://wvs.earthdata.nasa.gov/Fig.1 The study area and in-situ observation stationsThe location of mud area from reference [28], remote sensing map from https://wvs.earthdata.nasa.gov/

1.2 数据处理

悬浮体水样在实验室进行双滤膜抽滤测试，抽滤实验采用孔径0.45 μm的醋酸纤维滤膜进行双膜抽滤矫正，抽滤前将滤膜以45°C烘干24 h并称得上膜前质量（Wt1）和下膜前质量（Wb1），称重采用万分之一的电子天平，抽滤后再以45°C烘干24 h并称得上膜后质量（Wt2）和下膜后质量（Wb2），并记录抽滤水体体积（V）通过计算公式（1）计算悬浮体浓度（SSC）。

选取10个站位SSC测试结果，对采集悬浮体样品所在的垂向2 m范围内的悬浮体浊度进行平均计算，并与SSC进行拟合，拟合结果如图2所示，相关系数R2=0.83，拟合结果较好，可参考拟合结果，根据悬浮体浊度的观测结果，对SSC进行拟合判定。

图 2 悬浮体浓度（SSC）与浊度的相关性数据Fig.2 Correlation between suspended sediment concentration(SSC) and turbidity

图 3 各站位25 h流速矢量图潮汐和潮流类型出自文献[29]。Fig.3 Current vectors in 25 hours at each stationTidal and current types are from reference [29].

2 结果

2.1 潮流特征

渤海湾泥质区南部潮流以M2潮流为主[29-30]，从图3和图4的流速矢量图以及图3潮流类型分布[29]可以看出，A1、A2、A3、A5站的潮流类型为规则半日潮流，A4站潮流类型为不规则半日潮流，A1站为逆时针旋转流，A2站为E-W向往复流，A3、A4、A5站为顺时针旋转流。从图4水深变化以及图3潮汐类型分布[29]可以看出，A1、A2、A4站的潮汐类型为不规则半日潮，A3站为不规则全日潮，A5站处在不规则全日潮的位置，而在观测期间，A5站出现半日潮特征。各站位的流速东西分量u明显高于流速南北分量v（图4），u分量流速超过40 cm/s的历时超过观测历时的1/2，v分量流速低于40 cm/s的历时占观测历时的主要部分。A1站涨潮流为SW向，落潮流为NE向；A2站涨潮流为W向，落潮流为E向；A3站涨潮流和落潮流均有N向顺时针旋转为E向的现象；A4站涨潮流为SE向顺时针旋转为NW向，落潮流为NW向顺时针旋转为SE向；A5站涨潮流和落潮流均有NW向顺时针旋转一周回到NW向的现象。

图 4 各站位垂向平均流速矢量与流速各分量（u，v）时空分布图Fig.4 Vertical mean current vectors and temporal and spatial distributions of current components (u, v) at each station

2.2 温盐交换过程

从各站位温盐垂向剖面结构可以看出（图5），受各站位涨落潮流影响，水体的垂向温盐结构具有周期性分布的特点。整体来看（图6），水深最浅、最靠近岸边的A1站（水深5～15 m）温度最高，盐度最低，水深15～20 m的A2站和A5站温度和盐度都偏高，水深大于20 m的A3站和A4站温度偏低，盐度偏高。这说明渤海湾泥质区南部海域，15 m等深线以浅海域受周边径流高温淡水影响较大，15～20 m范围高温淡水与低温盐水交汇，盐度明显升高，20 m以深海域主要受渤海中部低温高盐水舌入侵影响，水体温度偏低。

图 5 各站位温盐结构时空分布图Fig.5 Spatial and temporal distribution of thermohaline structure at each station

图 6 各站位温度和盐度范围比较Fig.6 Comparison in the range of temperature and salinity at each station

A1站SW向涨潮流从湾外带来低温高盐的水体，NE向落潮流向湾外带走近岸的高温低盐的水体；湾口处A2站W向涨潮流从湾外带来低温高盐水体，E向的落潮流从湾内带出高温低盐水体，主要在表层水深10 m范围内传输，10 m以深水体为盐度较高、温度较低的水体；A3站的温盐结构在涨潮和涨平阶段，在表层10 m范围内层化明显，出现高温低盐水体，10 m以深水体垂向较均匀，在落潮和落平阶段，水体温盐垂向结构整体均匀，主要为温度和盐度都相对较高的水体；A4站和A5站的温盐结构较其他站位更稳定，温度自上而下逐渐降低，盐度自上而下逐渐升高，A4站（站位水深约25 m）在水深15 m（自海面向下）处温盐梯度较大，A5站（站位水深约20 m）则是在水深10 m处温盐梯度较大，两站连线与黄河口附近岸线和等深线接近垂直，呈NE向展布，说明夏季黄河口高温低盐的淡水向20 m以深的渤海中部海域传输的方式，主要为表层传输的羽状流形式。

2.3 SSC时空变化特征

从图7可以看出，渤海湾内（A1站）和湾口（A2站）处高浊度悬浮体多出现在涨急和落急时期，较高的流速对底质掀沙作用加强（图8a），其中A1站涨平后到落急期间SSC明显降低，可见，在渤海湾南部湾内的近岸海域，大部分底质再悬浮的泥沙还是随着涨潮流向湾内汇聚，方向多为SW向。黄河口SE向的A3、A4、A5三站SSC分布的潮周期特点较湾内和湾口处不显著，且水体较湾内和湾口清澈，SSC值偏低。A3站SSC高值出现在N向和E向潮流时期，A4和A5站SSC高值出现在NW向和SE向潮流时期。

整体来看（图8a），夏季渤海湾南部海域高于20 mg/L的SSC值出现在A2、A3站表层5 m水深范围内，A3、A4、A5站底层5 m水深范围内，A2站底层12 m范围内，以及A1站的整个垂向水体范围内。靠近底层的SSC值均高于表层SSC值（图8a），各站位靠近底层最高SSC值对比结果显示，湾内（A1站）及湾口（A2站）处最高，分别可达约130和80 mg/L，黄河口SE向20 m及更深（A4、A5站）处最小，均低于40 mg/L。SSC值高于20 mg/L的流速值基本位于50～70 cm/s区间（图8b），各站位流速基本低于100 cm/s，结合图8a，可知该区域底层高浊度的悬浮细颗粒物质的输运流速在50～100 cm/s区间，且具备这种输运条件的高浊度悬浮体多出现在渤海湾南部湾内和湾口海域，其向湾内汇沙的效应也符合中国东部陆架边缘海沉积物“夏储冬输”的季节性特征[31]。

图 7 各站位悬浮体浓度的时空分布特征与垂向平均流速矢量对比Fig.7 Comparison of spatial and temporal distribution characteristics of SSC at each station and vertical mean current vector

图 8 各站位悬浮体浓度与距底高度值（a）和流速值（b）的分布范围Fig.8 Distribution range of SSC and height from bottom (a) and velocity (b) at each station

3 讨论

3.1 悬浮体输沙率的潮周期变化特征

采用相对水深对瞬时物质输移量进行分解[32]，设x轴为纵向坐标，t为时间，z为相对水深（0≤z≤1），不计流速脉动项，将瞬时流速u(x,z,t)分解成垂向平均量项及其偏项之和：

ū和u′均分解为潮平均量项和潮变化项之和：

瞬时流速即为：

水深可分解为：

单宽潮周期平均输水量为：

由式（8）计算出来的即为一维垂向平均拉格朗日余流，其中分别为一维垂向平均欧拉余流和斯托克斯余流。

与余流的计算原理相似（式（2）—（5）），悬浮泥沙浓度SSC的函数c(x,z,t)可分解为：

单宽输沙率E为：

潮平均单宽瞬时输沙率T为：

湾内A1站（图9a）涨平期间的单宽输沙率为68～143 g·m-1·s-1，方向为偏S向，向岸输沙。落平期间单宽输沙率为41～93 g·m-1·s-1，方向为偏N向，向湾外输沙，单宽输沙率明显小于涨平期间。落潮时期输沙以E、SE向为主，涨潮时期输沙率方向以W、SW向为主，大小相当。整体来看（图10a），渤海湾南部湾内海域以向近岸输沙为主，潮平均单宽输沙率为7.8 g·m-1·s-1，方向为280°（以E向为0°的逆时针旋转角度），其中，T1+T2所代表的拉格朗日输沙率为7.6 g·m-1·s-1，方向为274°；T4和T5数量级相当，分别为0.6和0.3 g·m-1·s-1，方向为偏E向和偏SE向；其他分量的数量级较小，T3、T6数量级相当，T7和T8数量级相当，T3、T6、T7方向偏E向和偏SE向，T8方向偏N向。

湾口A2站（图9b）落平和落潮时期输沙方向偏E向，以朝湾外输沙为主，在涨平和涨潮期间输沙方向偏W向，以向湾内输沙为主，前者的单宽输沙率范围为24～281 g·m-1·s-1，后者的单宽输沙率范围为28～258 g·m-1·s-1，后 者 历 时 更 长。整 体 上（图10b），湾口处海域拉格朗日输沙率（T1+T2）为12.7 g·m-1·s-1，方 向 为329°；垂 向 净 环 流 的 贡 献项T5为5.5 g·m-1·s-1，方向为149°，方向与拉格朗日输沙率方向相反；整体来看A2站潮平均输沙率为7.2 g·m-1·s-1，方向为328°，输沙以向湾外SE向输沙为主；其他分量数量级较小，T4和T7数量级相当，T3、T6、T8数量级最小，T3、T6、T7方向偏SE向，T4和T8方向偏NW向。

A3站（图9c）落平时输沙以偏E向为主（18～46 g·m-1·s-1），涨 平 时 输 沙 以 偏SE向 为 主（10～54 g·m-1·s-1），涨 潮 时 输 沙 以 偏W向 为 主（12～59 g·m-1·s-1），落潮时输沙以偏N向和偏NE向为主（5～44 g·m-1·s-1）。整体来看（图10c），A3站潮平均单宽输沙率为4.7 g·m-1·s-1，方向为77°，其中，拉格朗日输沙率（T1+T2）为5.2 g·m-1·s-1，方向为50°；垂向净环流的贡献项T5为2.4 g·m-1·s-1，方向为166°；T3、T4数量级相当，T6、T7、T8数量级最小，T4、T7方向偏E向，T6偏SE向，T3、T8偏SW向。

A4站（图9d）涨潮时输沙偏SE向，单宽输沙率为8～49 g·m-1·s-1；涨平时输沙偏NW向，单宽输沙率为19～52 g·m-1·s-1；落潮时输沙偏NW向，单宽输沙率为29～56 g·m-1·s-1；落平时输沙偏N向和偏NE向，单宽输沙率为11～51 g·m-1·s-1。A4站潮平均单宽输沙率为5.2 g·m-1·s-1，方向为94°，其中，拉格朗日输沙率（T1+T2）为5.2 g·m-1·s-1，方向为91°，其他各分量数量级较小，T3和T7方向相反，数量级相当；T4数量级与T3和T7相当，方向为偏NW向；T5数量级比T3大，方向为偏W向；T6数量级与T3相当，方向为偏NE向；T8数量级比T3小，方向为偏SE向。

A5站（图9e）落潮时输沙偏NW向，单宽输沙率 为50～88 g·m-1·s-1；落 平 时 单 宽 输 沙 率 为32～78 g·m-1·s-1，方 向 自 偏NW向 顺 时 针 旋 转 至 偏SE向；涨潮时输沙偏SE向，单宽输沙率为37～105 g·m-1·s-1；涨平时单宽输沙率 为7～55 g·m-1·s-1，输沙偏NW向。A5站潮平均单宽输沙率为7.7 g·m-1·s-1，方 向 为102°，其 中，拉 格 朗 日 输 沙 率（T1+T2）为8.8 g·m-1·s-1，方向86°；垂向净环流的贡献项T5为2.5 g·m-1·s-1，方向208°；其他各分量数量级较小，其中T4数量级最大，方向偏E向，T3、T6、T7、T8数量级相当，T3、T7方向偏NE向，T6、T8方向偏SW向。

图 9 各站位潮平均单宽输沙率的变化Fig.9 Variation of tidal averaged suspended sediment flux per unit meter width at each station

图 10 各站位潮平均输沙率的各分量值与输沙角度Fig.10 Values of each component of tidal-averaged sediment transport rate and transport angles at each station

整体来看，渤海湾南部各站位潮平均单宽输沙率以拉格朗日输沙贡献最显著，在湾口泥质区南部和东南部水深15～25 m的海域，垂向净环流的影响较大，有抵消一部分拉格朗日输沙率的作用，且对潮平均单宽输沙率的影响比湾内和25 m以深海域的大，其中对A2站的拉格朗日输沙率的抵消作用最大，这跟A2站垂向净环流输沙率较大，且方向与拉格朗日输沙率相反有关。其他分量数量级较小，对潮平均单宽输沙率贡献较小。由于本次观测站位非同步连续观测，时间跨度较大，部分站位观测的是小潮时期的输沙率（A1、A2、A3站），另一部分站位观测的是大潮时期的输沙率（A4、A5站），因此对各站位之间输沙率的比较还需要未来进一步开展同步观测调查。

3.2 水体层化混合过程与输沙率之间的相关性

水体的垂向混合状态可以用Richardson数来体现[6,33-34]，其计算公式如下：

其中，分子表示水体密度梯度导致的层化强度，分母表示水体剪切力所导致的湍动混合程度，Ri＞0.25，代表水体以层化为主，Ri＜0.25，代表水体以混合为主。结果如图11所示，结合图5的温盐结构来看，A1、A2、A3、A4站受近岸高温低盐水体影响，这4个站位的密度结构不稳定，层化强度和Ri值多数时刻为负，相应时刻对应的混合程度较高。Ri值整体以小于0.25为主，说明渤海湾南部海域水体以混合为主，单宽输沙率降低的时刻，与水体层化程度和Ri值增加的时刻相吻合，如A1站的0时、14时、20时，A2站的5时、18时、23时，A3站的4时、14时，A4站的5时、17时、23时，A5站的6时、11时、12时、17时、18时、22时，可见水体的层化程度加强对各站位悬浮体输运均有一定的抑制作用，这与黄河口、加利福尼亚北部的大陆架和福宁湾海域前人研究结论相符[6,35-36]。

图 11 各站位水体Ri、层化程度、混合程度与E的变化Fig.11 Variation of Ri, stratification stability, turbulent intension, and E at each station

4 结论

（1）受周边径流高温淡水影响，夏季渤海湾南部及黄河三角洲北部15 m以浅海域水体温度较高，盐度较低；15～20 m范围内高温淡水与低温盐水交汇，水体盐度明显升高；20 m以深海域主要受渤海中部低温高盐水舌入侵影响，水体温度偏低。高温淡水以羽状流的方式自近岸向渤海中部传输。渤海湾内和湾口处高浊度悬浮体多出现在涨急和落急时期，较高的流速对底质的再悬浮作用加强，距底5 m水深范围内SSC值较高，其中湾内和湾口的底层SSC值最高，黄河口外NE向剖面20 m以深海域SSC值最小。

（2）渤海湾泥质区南部海域夏季单宽输沙率具有潮周期性特点。湾内海域涨潮和涨平期以向湾内近岸输沙为主，落潮和落平期以向湾外输沙为主，整体以向湾内近岸输沙为主，潮平均单宽输沙率为7.8 g·m-1·s-1，方向为280°。湾口A2站涨落潮流输沙方向相反，整体以SE向朝湾外近岸输沙为主，潮平均单宽输沙率为7.2 g·m-1·s-1，方向为328°。湾外A3站涨潮和涨平时期输沙方向基本相反，大小相当，落潮和落平时期输沙方向以偏E和偏NE向为主，潮平均单宽输沙率为4.7 g·m-1·s-1，方向为77°。黄河口SE向的2个站位涨潮时输沙偏SE向，涨平时输沙偏NW向，落潮时输沙偏NW向，落平时输沙偏N向和偏NE向，两站潮平均单宽输沙率方向以偏N向为主，A4站潮平均单宽输沙率为5.2 g·m-1·s-1，方向为94°，A5站潮平均单宽输沙率为7.7 g·m-1·s-1，方向为102°。

（3）渤海湾泥质区南部各站位夏季潮平均单宽输沙率以拉格朗日输沙贡献最显著，在湾口泥质区南部和东南部水深15～25 m海域，垂向净环流的影响较大，有抵消一部分拉格朗日输沙率的作用，且对潮平均单宽输沙率的影响比湾内和25 m以深海域的大，其他分量数量级较小，对潮平均单宽输沙率贡献较小。

（4）渤海湾泥质区南部海域夏季水体以混合为主，各站位水体层化程度加强对各站位悬浮体输运均有一定的抑制作用。