海洋与湖沼  2023, Vol. 54 Issue (4): 1015-1027   PDF    
http://dx.doi.org/10.11693/hyhz20221200322
中国海洋湖沼学会主办。
0

文章信息

张德强, 庞重光, 王延平, 陈祥舰. 2023.
ZHANG De-Qiang, PANG Chong-Guang, WANG Yan-Ping, CHEN Xiang-Jian. 2023.
渤海表层悬沙分布季节变化特征与控制因子的数值研究
NUMERICAL STUDY ON SEASONAL DISTRIBUTION OF SURFACE SUSPENDED SEDIMENT IN THE BOHAI SEA AND THE CONTROLLING FACTORS
海洋与湖沼, 54(4): 1015-1027
Oceanologia et Limnologia Sinica, 54(4): 1015-1027.
http://dx.doi.org/10.11693/hyhz20221200322

文章历史

收稿日期:2022-12-08
收修改稿日期:2023-02-08
渤海表层悬沙分布季节变化特征与控制因子的数值研究
张德强1,2, 庞重光1,3,4, 王延平5,6, 陈祥舰1,7     
1. 中国科学院海洋研究所海洋环流与波浪重点实验室 山东青岛 266071;
2. 中国科学院大学 北京 100049;
3. 海洋动力过程与气候功能实验室 崂山实验室 山东青岛 266237;
4. 中国科学院海洋科学大数据研究中心 山东青岛 266071;
5. 北京应用气象研究所 北京 100029;
6. 地理信息工程国家重点实验室 陕西西安 710054;
7. 山东科技大学数学与系统科学学院 山东青岛 266590
摘要:针对整个渤海海域悬浮泥沙季节变化及其影响机制的数值研究相对缺乏且机制尚不清晰, 基于ROMS三维海洋模型对渤海海域水动力环境与悬沙分布开展数值模拟。模拟结果显示, 渤海海峡环流终年“北进南出”, 夏季环流明显强于冬季, 并呈现外围逆时针环、内部顺时针环的“双环”结构。渤海中部海域在夏季存在明显的温跃层现象, 其强度分布与等深线较为一致, 温跃层在4月开始形成, 7月最强。渤海表层悬沙分布具有显著的季节变化, 冬季悬沙浓度最大, 秋季次之, 春季再次之, 夏季最小。控制悬沙浓度的波流底切应力在秋冬季节较大, 春夏相对较小, 且流致切应力始终在波流切应力中占主导地位。秦皇岛海域悬沙浓度常年偏低的主要原因是位于M2无潮点附近, 属于弱潮流区, 底层流速相对较小, 底部沉积物发生再悬浮概率较小。夏季温跃层的存在在一定程度上减小了底边界层流速, 增大了流速的垂向梯度, 对底部悬浮泥沙的向上扩散有明显的抑制作用。因此, 温跃层的存在是造成夏季表层悬沙浓度最低的重要原因。
关键词ROMS模型    渤海    温跃层    表层悬沙浓度    季节变化    
NUMERICAL STUDY ON SEASONAL DISTRIBUTION OF SURFACE SUSPENDED SEDIMENT IN THE BOHAI SEA AND THE CONTROLLING FACTORS
ZHANG De-Qiang1,2, PANG Chong-Guang1,3,4, WANG Yan-Ping5,6, CHEN Xiang-Jian1,7     
1. Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China;
2. University of the Chinese Academy of Sciences, Beijing 100049, China;
3. Laboratory for Ocean and Climate Dynamics, Laoshan Laboratory, Qingdao 266237, China;
4. Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China;
5. Beijing Institute of Applied Meteorology, Beijing 100029, China;
6. State Key Laboratory of Geo-information Engineering, Xi'an 710054, China;
7. School of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
Abstract: Due to the lack of research on the numerical simulation of suspended sediment distribution in the whole Bohai Sea, the hydrodynamic environment and the distribution of suspended sediment concentration in the Bohai Sea was numerically simulated in the three-dimensional ocean model ROMS (regional ocean modeling system). Results show that the Bohai Strait circulation presents a "north in south out" pattern all year round. The circulation in summer is obviously stronger than that in winter, and presents a "double-ring" structure with an outer counterclockwise ring and an inner clockwise ring. There is a strong thermocline in the middle of the Bohai Sea in summer, and the distribution of its strength is close to the isobath. The thermocline begins to form in April and is the strongest in July. The distribution of surface suspended sediment concentration in the Bohai Sea has distinct seasonal change, with the highest concentration in winter, followed by autumn, spring, and summer. The wave-current bottom shear stress controlling suspended sediment concentration is larger in autumn and winter, but smaller in spring and summer, and the current-induced shear stress always plays a dominant role in wave-current shear stress. The main reason why suspended sediment concentration off the Qinhuangdao is low all the year round is that the current velocity at the bottom is relatively small, resulting in the weak sediment resuspension, as the location is adjacent to an amphidromic point of the M2 tidal constituent. The thermocline in spring and summer reduces the current velocity at the bottom boundary layer, and increases the vertical gradient of the velocity to a certain extent, and further suppresses the upward diffusion of the bottom sediment. Therefore, the existence of thermocline is one of the controlling factors of the lowest surface suspended sediment concentration in summer.
Key words: ROMS (regional ocean modeling system)    the Bohai Sea    thermocline    surface suspended sediment concentration    seasonal variation    

渤海作为我国北部典型的半封闭陆架海, 是我国最早进行悬沙分布和输运研究的海域。前人首先进行了大量的悬沙浓度现场观测, 得到了渤海悬沙浓度的空间分布特征及其季节变化规律(秦蕴珊等, 1982; 中国科学院海洋研究所海洋地质研究室, 1985; Jiang et al, 2004)。通过与渤海水动力环境相结合, 能够大体明确悬沙输运的模式和机理。早期应用的数值模型不论是二维模型(Fang et al, 1985; Zhang et al, 1986; 董礼先等, 1989)还是三维模型(Feng, 1987; Graber et al, 1989; Zhao et al, 1993; 王凯等, 1999), 大都仅考虑潮流或风浪的作用。随着数值模拟技术的发展, 对渤海悬沙输运的数值模拟开始综合考虑浪-潮-流等水动力耦合过程对悬沙输运的影响(江文胜等, 2000; 庞重光等, 2004)。

渤海水动力环境的研究和模拟是悬沙输运数值模拟的基础。渤海水动力环境的研究自20世纪60年代以来开展至今, 其中基于各类数值模式的研究较多(Fang et al, 1985, 2000; Zhang et al, 1986; Feng, 1987; Zhao et al, 1993; Guan, 1994; 乔方利, 2012), 数值模拟结果揭示了渤海水动力环境的基本特征。Xia等(2019)利用高分辨率波浪-潮汐-环流耦合模式模拟渤海夏季的环流, 发现辽东湾西部至黄河口近岸海域沿温度锋面存在明显的密度流, 在40°N以南海域环流呈双环流结构。马伟伟等(2016)利用ROMS (regional ocean modeling system)海洋模式对渤海冬季流场进行高分辨率的数值模拟, 发现冬季深度平均流场在辽东湾为顺时针环流, 且在渤海中央表现为顺时针的环流结构。相比于研究较为成熟的渤海水动力特征, 渤海悬沙的研究大都集中在黄河口, 渤海海峡等热点区域(Bi et al, 2011; Yang et al, 2011; Wang et al, 2020)。受现场观测范围的限制, 整个渤海海域的悬沙分布研究一般通过卫星遥感数据的反演得到(庞重光等, 2013, 2014; 周舟等, 2017; Zhao et al, 2022)。庞重光等(2014)通过分析长时间序列悬沙浓度和风场遥感反演资料, 指出渤海不同海域表层悬沙浓度差别较大, 高浓度集中在辽东湾、渤海湾和莱州湾, 低浓度主要位于秦皇岛海域以及渤海中部等。庞重光等(2013)基于SeaWiFs反演的渤海表层悬沙浓度资料, 通过经验正交函数分解, 展示了渤海表层悬浮泥沙的空间分布特征及其随时间的变化; 其中第一空间模态显示了渤海表层悬沙的季节变化特征。周舟等(2017)利用2003~2014年的MODIS (Moderate Resolution Imaging Spectroradiometer)卫星遥感数据, 对渤海表层悬浮物浓度的空间分布形态进行反演, 结果显示表层悬浮物浓度在时间尺度上有明显的季节变化, 且渤海表层悬浮物从2003~2014年呈现下降的趋势。近年来, 相对于卫星遥感资料反演悬沙浓度, 利用数值模式对整个渤海海域悬沙分布的研究相对较少(Lu et al, 2011; Liu et al, 2020; Wang et al, 2020)。

利用数值模拟方法不仅可以再现渤海悬沙浓度分布的时空特征, 同时数值模拟也是揭示悬沙浓度分布的动力成因的有效手段。本文基于ROMS海洋模式, 对渤海的水动力环境和悬沙分布进行气候态数值模拟, 重点展示渤海海域的冬夏季流场特征、悬沙分布, 并对悬沙分布及其影响机制进行讨论。

1 研究方法

本文所采用的ROMS模式是基于三维非线性斜压原始方程的开源海洋模型, 模拟范围为36°~41°N, 117°~123°E, 如图 1所示。地形数据主要采用ETOPO1数据, 水平方向上采用正交经纬网格, 网格数为73×61, 水平分辨率为5′, 垂向网格采用σ坐标, 共分为20层, 最小水深设置为1 m。模型计算采用冷启动, 初始水位和流速均设为0, 初始场的温盐数据来自WOA13 (World Ocean Atlas 2013), 在初始场中将沉积物分为3类, 分别为黏土、粉砂和砂, 三种沉积物的占比分布来自Wang等(2020)的研究结果。其大气强迫场数据来自COADS (comprehensive ocean atmosphere data set)资料, 边界条件采用SODA (simple ocean data assimilation)气候态平均场, 潮强迫是在开边界上采用10个分潮(M2, S2, N2, K2, K1, O1, P1, Q1, Mf, Mm)的调和常数来确定水位变化, 其数据来源于全球潮波模式TPXO 7.0。波浪强迫是将SWAN (simulating waves nearshore)单独计算的有效波高、波向、波长以及波周期等波浪参数作为外在强迫添加到ROMS中。在边界强迫中, 开边界水位采用Chapman边界条件, 正压流速采用Flather边界条件, 三维斜压流速及温盐场均采用Rad边界条件。在本模型中, 由于在渤海入海的河流中, 辽河、滦河以及海河流量远小于黄河, 因此仅考虑黄河径流及悬浮物的输运作用。对于黄河径流量和悬浮物的数据来自于利津水文站近30 a的气候态数据, 且黄河中的黏土、粉砂和砂的占比分别为20%、70%和10% (Wang et al, 2020)。

图 1 模拟区域的水深和站位(a)及水平网格(b) Fig. 1 Water depth and stations (a) and horizontal grid (b) in the model simulated region

本文用气候态数据将模式运行了2 a, 取第2年的12个月份的月平均值来表征每个月的特征。本文采用了2018年7月在渤海进行观测的大面站温盐深(conductivity temperature depth, CTD)数据与模拟的气候态的月平均值进行对比, 同时采用2017年2~3月与2018年8月在秦皇岛海域的实测数据, 对其求10 d的平均值, 与模拟10 d的平均值进行对比。

2 模式结果验证 2.1 潮位对比

对于潮汐的准确模拟是正确模拟渤海水动力环境的前提。取第二年1月份的水位数据, 采用调和分析方法分离出潮位的调和常数。由于渤海的主要半日分潮和全日分潮分别为M2和K1, 因此本文仅对M2和K1进行分析。图 2c2d分别为M2和K1分潮的同潮图。

图 2 渤海M2 (a, c)和K1 (b, d)分潮的同潮图分布 Fig. 2 The amplitude and phase distribution of M2 (a, c) and K1 (b, d) in the Bohai Sea 注: a与b分别为渤海图集(Fang, 1986)中M2和K1分潮的同潮分布图, c和d为模式模拟的同潮图; 虚线表示振幅(单位: cm), 实线表示迟角(单位: °)

通过与Fang (1986)的渤海图集中同潮图进行对比, 可以发现两者之间分潮M2与K1的无潮点位置, 振幅以及相位的分布均较为一致。通过对模拟的调和常数与7个渤海验潮站实测的M2和K1分潮的调和常数进行定量对比, M2分潮振幅误差的绝对平均为8.2 cm, 相对误差的绝对平均为11.8%; 迟角误差的绝对平均为14.3°, 相对误差的绝对平均为8.7%。K1分潮振幅误差的绝对平均为0.7 cm, 相对误差的绝对平均为4.6%; 迟角误差的绝对平均为22.3°, 相对误差的绝对误差为19.4%。以上结果表明, 该模型对渤海潮汐的模拟是准确的。

2.2 流场对比

图 3a图 3b分别是渤海冬季(1月)和夏季(7月)的气候态垂向平均流场, 从模拟结果来看, 在渤海海峡终年呈现“北进南出”的结构; 冬季沿岸线流动的海流最强, 辽东湾以及渤海中部均呈现顺时针环流结构, 渤海湾和莱州湾海流逆时针流动, 整体海流流速较小, 一般都在0.05 m/s以下, 这与已有的数值模拟结果较为类似(毕聪聪等, 2015; 马伟伟等, 2016; Jiang et al, 2020; Liu et al, 2021)。夏季环流明显强于冬季, 流速很多在0.05~0.1 m/s。夏季, 20 m等深线以深的渤海中部, 环流呈现“双环”结构, 即外部的逆时针环包围着内部的顺时针环, 这与Xia等(2019), Liu等(2021)的研究结果一致。本文的气候态模拟结果可以再现渤海环流的基本结构及其季节变化。

图 3 渤海1月(a)和7月(b)气候态的深度平均流场 Fig. 3 The climatological depth-averaged current fields in the Bohai Sea in January (a) and July (b), respectively
3 结果与讨论 3.1 渤海温跃层的季节变化

为了研究渤海温跃层的季节性变化, 将模式模拟的气候态结果与实测数据进行对比。实测数据分为两部分, 第一部分是利用在S1站点(图 1a黑色空心圆所示)于2017年2月15日~3月9日(冬季), 2019年8月15~29日(夏季)进行的海床基平台观测数据, 并对其温度取10 d平均, 与模拟的气候态温度进行对比, 对比结果如图 4a所示; 第二部分是利用2018年7月对渤海进行的大面站观测数据与模式模拟结果(月平均)进行对比, 如图 5所示, 站点B29~B30, B53~B59 (红色空心圆)位置如图 1a所示。

图 4 S1站点模式模拟温度(10 d平均)底层(第一层, 第二层)、表层(二十层)与观测温度(10 d平均)的时间序列(a)以及S1站点8月垂向温度分布(b) Fig. 4 The time series of the bottom layers (Layer 1, Layer 2) and surface layer (Layer 20) in model simulation results (10-day average) and observed values (10-day average) at S1 (a), and vertical temperature distribution at S1 in August (b)

图 5 模式模拟温度与站位观测温度的垂向分布 Fig. 5 The vertical distribution of model simulated temperature and observed temperature

图 4a可以看出, 模式对于S1站点季节性变化的模拟准确。图 4b中的S1站点CTD实测数据在夏季的垂向分布显示, 表层温度为27 ℃, 底层温度为24.5 ℃, 没有明显的温跃层存在, 这也与图 4a中的模拟结果相对应。图 4a的模拟温度显示底层温度(Layer 1, Layer 2)与表层温度(Layer 20)相差无几。

图 5显示, 在渤海水深较深的海域, 夏季存在明显的温跃层。由于大面站观测获得的温度是瞬时值, 且这些站点的温度大都是白天所测, 故测得的温度普遍高于模式模拟的温度, 且温跃层以上温度更高, 垂向变化更大。而模式计算的温度是7月的平均值, 温度略低且垂向变化范围较小。此外, 模式模拟的盐度在31~32, 而B29~B30, B53~B59站点所测的盐度大都也在31~32, 进一步验证了模式模拟的准确性, 故可以用模式模拟的结果来分析渤海温度的季节性变化。

由于温跃层的存在使得表层温度与底层温度差距较大, 因此采用表层温度与底层温度的差值来体现温跃层的强度大小与位置分布, 其1~12月份的空间分布如图 6所示。图 6显示, 温跃层在4月开始形成, 7月达到最强, 在9月之后, 开始慢慢消退, 10月已经基本无温跃层的存在, 这与屠金钊(1992)的研究结果较为一致, 且在温跃层强盛期(6~8月)表底层温差的分布与等深线的分布较为一致。

图 6 表层温度与底层温度差值的空间分布 Fig. 6 Spatial distribution of temperature difference between the surface layer and the bottom layer
3.2 渤海悬沙分布的季节变化

本文参考Bian等(2013)的泥沙参数方案, 对泥沙各级粒径参数的相关设置见表 1。数值模拟的渤海表层悬沙浓度的季节性分布如图 7所示, 渤海表层的悬沙分布具有明显的季节变化。整体上讲, 渤海的悬沙浓度在冬季(12~2月)最大, 秋季(9~11月)次之, 春季(3~5月)再次之, 夏季(6~8月)最小, 且不论在哪个季节, 悬沙浓度在辽东湾、渤海湾以及莱州湾均相对较高, 大都在70 mg/L以上, 悬沙的时空分布特征与庞重光等(2013)利用遥感资料对渤海表层悬沙的反演结果基本一致。图 7还显示, 秦皇岛海域悬沙浓度显著低于其他海域, 常年偏低, 其表层悬沙浓度基本都维持在20 mg/L以下, 这一现象也与前人对渤海悬沙浓度的反演结果类似(Jiang et al, 2004; Wang et al, 2014; 王震等, 2016; Zhao et al, 2022)。通过对模式模拟的表层悬沙浓度和卫星遥感反演的海表悬沙浓度进行对比, 如图 8所示(位置点A~F如图 1a所示), 发现在渤海的不同位置, 如典型的高、低悬沙浓度点, 模式均能够模拟出表层悬沙浓度的季节变化特征, 且大致与卫星遥感反演的悬沙浓度值一致, 冬季浓度最高, 夏季浓度最低。

表 1 ROMS泥沙模块中各级粒级泥沙参数 Tab. 1 Sediment parameters in the ROMS Sediment Model
泥沙类型 泥沙粒径/mm 沉降速度/(mm/s) 临界切应力/(N/m2) 侵蚀速率/[kg/(m2·s)] 密度/(kg/m3) 疏松性
细砂 0.200 15.3 0.16 2.5×10−5 2 650 0.4
粉砂 0.015 0.4 0.10 2.5×10−5 2 650 0.4
黏土 0.004 0.1 0.05 2.5×10−5 2 650 0.4

图 7 渤海表层悬沙浓度不同季节的空间分布 Fig. 7 Spatial distribution of surface suspended sediment concentration in the Bohai Sea in different seasons

图 8 模式模拟的表层悬沙浓度与多年平均的卫星遥感反演的海表悬沙浓度的对比 Fig. 8 Comparison of surface suspended sediment concentrations simulated by model and retrieved from multi-year average satellite remote sensing
3.3 渤海底部切应力的季节变化

泥沙模型分别计算了时间步长为3 min的流致切应力、波致切应力以及波流相互作用产生的最大切应力。通过对剪切应力求季节平均, 得到不同季节剪切应力的大小, 波流相互作用切应力、流致切应力和波致切应力分别如图 9~11所示。从整体来看, 由波流相互作用引起的最大剪切应力在冬季最大、秋季次之、春季和夏季相对较小。流致切应力与波流切应力几乎相当, 波致切应力仅仅在沿岸的浅水区域(20 m等深线以浅)发挥一定的作用, 其中在辽东湾和莱州湾较为显著, 而在深水区, 波致切应力与流致切应力相比, 几乎可以忽略不计, 这与Wang等(2020)的研究较为一致。从波流切应力的季节分布来看, 辽东湾东部的波流切应力在冬季明显高于夏季, 冬季波流切应力基本在0.2 N/m2以上, 夏季波流切应力主要在0.12~0.20 N/m2。但在辽东湾北部, 波流切应力的季节变化不大, 常年维持在0.3 N/m2以上, 这也是辽东湾海域常年悬沙浓度较高的原因。在渤海湾, 波流切应力在夏季甚至要大于冬季, 这主要是因为夏季海流相对冬季更强。渤海湾与辽东湾都是流致切应力占主导, 而莱州湾的流致剪切应力相对较小, 其悬沙浓度则是波流相互作用与黄河径流输入共同作用的结果。此外, 秦皇岛海域其底部的波流切应力常年保持在0.06 N/m2以下, 远远小于其他海域底部的波流切应力。因此, 之所以秦皇岛海域悬沙浓度常年低于其他海域, 最根本的原因在于秦皇岛海域位于M2无潮点附近, 属于弱潮流区, 底层流速相对较小, 流致切应力相应较小, 底部沉积物发生再悬浮的概率同步较小。

图 9 波流相互作用下季节平均的最大剪切应力的空间分布 Fig. 9 Spatial distribution of seasonal mean maximum shear stress under wave-current interaction

图 10 季节平均的流致切应力的空间分布 Fig. 10 Spatial distribution of seasonal mean current-induced shear stress

图 11 季节平均的波致切应力的空间分布 Fig. 11 Spatial distribution of seasonal mean wave-induced shear stress
3.4 渤海夏季表层悬沙分布的控制因子

渤海表层悬沙浓度(图 7)和波流底切应力(图 9)的季节分布显示, 在冬季、秋季以及春季, 波流底切应力的空间分布与表层的悬沙分布较为一致。据此推论, 冬季、秋季以及春季渤海表层悬沙大都来源于底部再悬浮泥沙的垂向扩散。然而, 在夏季, 其波流底切应力的空间分布与春季分布相当接近(图 9b9c), 但其表层悬沙浓度的分布却明显低于春季(图 7b7c)。本文数值模式的垂向分层为20层, 因此能够清楚地展示季节平均温度、悬沙浓度和流速大小的垂向分布, 分别如图 12~14所示。如图 12a所示, 冬季渤海海域海水温度在垂向几乎没有变化, 与其对应的悬沙浓度在垂向上也基本保持不变(图 13a), 流速大小自海表向下略有衰减(图 14a), 如流速最大的辽东湾北部位置, 流速大小从海表的0.41 m/s, 减小到海水下层第4层的0.35 m/s, 继续衰减到底层第1层的0.30 m/s。垂向基本均匀分布的水温和悬沙浓度表明在渤海, 冬季整个深度都是混合层, 垂向混合强度大、混合较为充分。然而, 在夏季, 辽东湾、渤海湾、莱州湾以及其他近岸海域由于水深较浅, 温度在垂向上依然没有明显变化, 因此浅水区的变化主要还是受动力条件的季节性波动控制, 但在渤海中部, 表层温度与底层温度差别显著, 存在明显的温跃层现象, 如图 12c所示。夏季, 渤海中部的表层温度一般在22~26 ℃, 而底部温度则大都在12~20 ℃。图 12c还显示, 从底层到中层Layer 1~Layer 10, 即温跃层之下, 温度的上升幅度较小, 但相对应的悬沙浓度在垂向上递减的幅度却较大(图 13c)。主要原因是温跃层的存在降低了风、海浪等外力对海底的扰动, 抑制了均温层水体在垂向的混合, 不利于底层再悬浮泥沙的向上扩散。夏季, 虽然水体在垂向上混合强度明显减弱, 但其平流项却更加发育, 使得中上层海水流速增大, 不论是季节平均的流速(图 14)还是冬夏季深度平均流场(图 3), 都是夏季流速大于冬季, 这与孙运佳等(2016)的研究一致。

图 12 季节平均温度的垂向分布 Fig. 12 Vertical distribution of seasonal mean temperature

图 13 季节平均悬沙浓度的垂向分布 Fig. 13 Vertical distribution of seasonal mean suspended sediment concentration

图 14 季节平均流速大小的垂向分布 Fig. 14 Vertical distribution of seasonal mean current velocity

图 12c展示了夏季温跃层中心所在层位(Layer 13), 该层温度位于上层(Layer 16)和下层(Layer 10)之间, 该层的悬沙浓度相对于下层(Layer 10)也有明显的降低, 除三湾外, 悬沙浓度基本上都低于10 mg/L, 如图 13c所示。由于温跃层具有层结强, 垂向扩散系数低的特征, 温跃层的存在对于底部悬浮颗粒物向表层的扩散起到明显的抑制作用(乔璐璐等, 2010; 雷栋等, 2019), 故温跃层上层的悬沙浓度(Layer 16和Layer 19)显著低于下层的悬沙浓度(Layer 1, Layer 4, Layer 7和Layer 10)。

水体层化现象对流速的垂向梯度具有明显的增强作用, 如图 14所示。在水体下层Layer 1~Layer 10, 夏季垂向流速梯度明显大于冬季, 体现在水体中层(Layer 10), 夏季流速明显大于冬季, 但在水体底层(Layer 1), 夏季流速要小于冬季。这一现象在辽东湾东部最为显著, 夏季中层水体流速约为0.48 m/s, 明显大于冬季中层水体的0.39 m/s, 而夏季的底层流速约为0.28 m/s, 小于冬季对应的底层流速0.32 m/s。水体的层化在一定程度上减弱了底层流速, 从而间接影响了底部沉积物的再悬浮, 这与齐昱恺等(2021) 在三沙湾的观测结果一致。与夏季相比, 春季在一定程度上也有温跃层现象(图 12b), 但强度较低, 这也使得春季的悬沙浓度(图 13b)比秋季(图 13d)更低, 因为秋季水体混合均匀, 几乎不存在温跃层的现象。

4 结论

本文基于ROMS模式对渤海水动力环境与悬沙分布开展数值模拟, 得出如下结论:

在渤海海峡环流终年“北进南出”; 夏季环流明显强于冬季, 并呈现外围逆时针环、内部顺时针环的“双环”结构。在平均深度仅为18 m的渤海, 夏季存在明显的温跃层现象, 温跃层在4月开始形成, 7月达到最强, 9月开始慢慢消退, 温跃层强度分布与等深线较为一致。

渤海悬沙分布具有明显的季节变化, 冬季悬沙浓度最大, 秋季次之, 春季再次之, 夏季最小。控制悬沙浓度的波流底切应力, 秋冬季明显大于春夏季, 且不论哪个季节, 流致切应力在波流剪切应力中占主导地位, 波致切应力仅仅在浅水海域有一定的作用。秦皇岛海域由于位于M2无潮点附近, 属于弱潮流区, 底层流速相对较小, 相应的流致切应力较小, 故底部沉积物发生再悬浮的程度较小, 这是秦皇岛海域常年悬沙浓度较低的主要原因。

春夏季温跃层的存在对悬沙的垂向交换具有重要影响。温跃层在一定程度上减小了底边界层流速, 增大了流速的垂向梯度; 温跃层的形成对底部悬浮颗粒物向上层的扩散起到明显的抑制作用。夏季, 温跃层的存在是中上层悬沙浓度维持较低的重要原因。但在浅水区, 无论哪个季节, 温度在垂向上都没有明显的变化, 悬沙浓度的变化主要还是受动力条件的季节性波动控制。

参考文献
马伟伟, 万修全, 万凯, 2016. 渤海冬季风生环流的年际变化特征及机制分析. 海洋与湖沼, 47(2): 295-302
王凯, 方国洪, 冯士, 1999. 渤海、黄海、东海M2潮汐潮流的三维数值模拟. 海洋学报, 21(4): 1-13
王震, 乔璐璐, 王云飞, 2016. 东中国海表层悬浮体浓度卫星遥感反演研究进展. 沉积学报, 34(2): 292-307
中国科学院海洋研究所海洋地质研究室, 1985. 渤海地质. 北京: 科学出版社, 50-56
毕聪聪, 鲍献文, 万凯, 2015. 渤海盐度年代际变异对环流结构的影响. 中国海洋大学学报, 45(1): 1-8
乔方利, 2012. 中国区域海洋学: 物理海洋学. 北京: 海洋出版社
乔璐璐, 李广雪, 邓声贵, 等, 2010. 夏季渤海湾中北部悬浮体分布. 海洋地质与第四纪地质, 30(3): 23-30
齐昱恺, 涂俊彪, 于俊杰, 等, 2021. 三沙湾夏季底边界层动力过程及悬沙输运特征. 海洋科学, 45(6): 1-12
江文胜, 孙文心, 2000. 渤海悬浮颗粒物的三维输运模式Ⅰ. 模式. 海洋与湖沼, 31(6): 682-688
孙运佳, 童朝锋, 孟艳秋, 2016. 水体层化对流速结构及悬沙分布影响. 中国港湾建设, 36(2): 1-4
周舟, 张万磊, 江文胜, 等, 2017. 渤海表层悬浮物浓度长期变化(2003—2014年)的卫星反演研究. 中国海洋大学学报, 47(3): 10-18
庞重光, 于炜, 2013. 渤海表层悬浮泥沙的空间模态及其时间变化. 水科学进展, 24(5): 722-727
庞重光, 白学志, 胡敦欣, 2004. 渤、黄、东海海流和潮汐共同作用下的悬浮物输运、沉积及其季节变化. 海洋科学集刊, 46: 32-41
庞重光, 李坤, 于炜, 2014. 渤海表层悬沙的时空分布特征及其动力成因. 海洋科学进展, 32(4): 450-458
秦蕴珊, 李凡, 1982. 渤海海水中悬浮体的研究. 海洋学报, 4(2): 191-200
屠金钊, 1992. 渤、黄海温度跃层的分布与变化. 海洋通报, 11(4): 27-32
董礼先, 苏纪兰, 王康墡, 1989. 黄渤海潮流场及其与沉积物搬运的关系. 海洋学报, 11(1): 102-114
雷栋, 王珍岩, 李文建, 2019. 夏季东海陆架不同水团对悬浮体物质组成及粒度分布的影响. 海洋与湖沼, 50(4): 740-751
BI N S, YANG Z S, WANG H J, et al, 2011. Seasonal variation of suspended-sediment transport through the southern Bohai Strait. Estuarine, Coastal and Shelf Science, 93(3): 239-247
BIAN C W, JIANG W S, GREATBATCH R J, 2013. An exploratory model study of sediment transport sources and deposits in the Bohai Sea, Yellow Sea, and East China Sea. Journal of Geophysical Research: Oceans, 118(11): 5908-5923
FANG G H, 1986. Tide and tidal current charts for the marginal seas adjacent to China. Chinese Journal of Oceanology and Limnology, 4(1): 1-16
FANG Y, FANG G H, ZHANG Q H, 2000. Numerical simulation and dynamic study of the wintertime circulation of the Bohai Sea. Chinese Journal of Oceanology and Limnology, 18(1): 1-9
FANG G H, YANG J F, 1985. A two-dimensional numerical model of the tidal motions in the Bohai Sea. Chinese Journal of Oceanology and Limnology, 3(2): 135-152
FENG S Z, 1987. A three-dimensional weakly nonlinear model of tide-induced Lagrangian residual current and mass-transport, with an application to the Bohai Sea. Elsevier Oceanography Series, 45: 471-488
GRABER H C, BEARDSLEY R C, GRANT W D, 1989. Storm-generated surface waves and sediment resuspension in the East China and Yellow Seas. Journal of Physical Oceanography, 19(8): 1039-1059
GUAN B X, 1994. Patterns and structures of the currents in Bohai, Huanghai and East China Seas [M] // ZHOU D, LIANG Y B, ZENG C K. Oceanology of China Seas. Dordrecht: Springer: 17-26.
JIANG M, PANG C G, LIU Z L, et al, 2020. Sediment resuspension in winter in an exceptional low suspended sediment concentration area off Qinhuangdao in the Bohai Sea. Estuarine, Coastal and Shelf Science, 245: 106859
JIANG W S, POHLMANN T, SUN J, et al, 2004. SPM transport in the Bohai Sea: field experiments and numerical modeling. Journal of Marine Systems, 44(3/4): 175-188
LIU H W, PANG C G, YANG D Z, et al, 2021. Seasonal variation in material exchange through the Bohai Strait. Continental Shelf Research, 231: 104599
LIU X M, QIAO L L, ZHONG Y, et al, 2020. Pathways of suspended sediments transported from the Yellow River mouth to the Bohai Sea and Yellow Sea. Estuarine, Coastal and Shelf Science, 236: 106639
LU J, QIAO F L, WANG X H, et al, 2011. A numerical study of transport dynamics and seasonal variability of the Yellow River sediment in the Bohai and Yellow seas. Estuarine, Coastal and Shelf Science, 95(1): 39-51
WANG C H, LIU Z Q, HARRIS C K, et al, 2020. The impact of winter storms on sediment transport through a narrow strait, Bohai, China. Journal of Geophysical Research: Oceans, 125(6): e2020JC016069
WANG H J, WANG A M, BI N S, et al, 2014. Seasonal distribution of suspended sediment in the Bohai Sea, China. Continental Shelf Research, 90: 17-32
XIA C S, GUO J S, WANG G S, et al, 2019. The summertime circulation of the Bohai Sea simulated from a high-resolution wave-tide-circulation coupled model. Acta Oceanologica Sinica, 38(1): 32-37
YANG Z S, JI Y J, BI N S, et al, 2011. Sediment transport off the Huanghe (Yellow River) delta and in the adjacent Bohai Sea in winter and seasonal comparison. Estuarine, Coastal and Shelf Science, 93(3): 173-181
ZHANG S Z, WANG H T, FENG S Z, et al, 1986. A simulation of the residual flow in the Bohai Sea [M] // VAN DE KREEKE J. Physics of Shallow Estuaries and Bays. New York: Springer-Verlag: 114-119.
ZHAO G B, JIANG W S, WANG T, et al, 2022. Decadal variation and regulation mechanisms of the suspended sediment concentration in the Bohai Sea, China. Journal of Geophysical Research: Oceans, 127(3): e2021JC017699
ZHAO J P, SHI M C, 1993. Numerical modelling of three-dimension characteristics of wind-driven current in the Bohai Sea. Chinese Journal of Oceanology and Limnology, 11(1): 70-79