渤海表层悬沙分布季节变化特征与控制因子的数值研究.pdf
《渤海表层悬沙分布季节变化特征与控制因子的数值研究.pdf》由会员分享,可在线阅读,更多相关《渤海表层悬沙分布季节变化特征与控制因子的数值研究.pdf(13页珍藏版)》请在咨信网上搜索。
1、第 54 卷 第 4 期 海 洋 与 湖 沼 Vol.54,No.4 2 0 2 3 年7 月 OCEANOLOGIA ET LIMNOLOGIA SINICA Jul.,2023 *中国科学院科研仪器设备研制项目,YJKYYQ20190047 号。张德强,博士研究生,E-mail: 通信作者:庞重光,研究员,E-mail: 收稿日期:2022-12-08,收修改稿日期:2023-02-08 渤海表层悬沙分布季节变化特征与 控制因子的数值研究*张德强1,2 庞重光1,3,4 王延平5,6 陈祥舰1,7(1.中国科学院海洋研究所海洋环流与波浪重点实验室 山东青岛 266071;2.中国科学院大学
2、 北京 100049;3.海洋动力过程与气候功能实验室 崂山实验室 山东青岛 266237;4.中国科学院海洋科学大数据研究中心 山东青岛 266071;5.北京应用气象研究所 北京 100029;6.地理信息工程国家重点实验室 陕西西安 710054;7.山东科技大学数学与系统科学学院 山东青岛 266590)摘要 针对整个渤海海域悬浮泥沙季节变化及其影响机制的数值研究相对缺乏且机制尚不清晰,基于 ROMS 三维海洋模型对渤海海域水动力环境与悬沙分布开展数值模拟。模拟结果显示,渤海海峡环流终年“北进南出”,夏季环流明显强于冬季,并呈现外围逆时针环、内部顺时针环的“双环”结构。渤海中部海域在夏
3、季存在明显的温跃层现象,其强度分布与等深线较为一致,温跃层在 4 月开始形成,7 月最强。渤海表层悬沙分布具有显著的季节变化,冬季悬沙浓度最大,秋季次之,春季再次之,夏季最小。控制悬沙浓度的波流底切应力在秋冬季节较大,春夏相对较小,且流致切应力始终在波流切应力中占主导地位。秦皇岛海域悬沙浓度常年偏低的主要原因是位于 M2无潮点附近,属于弱潮流区,底层流速相对较小,底部沉积物发生再悬浮概率较小。夏季温跃层的存在在一定程度上减小了底边界层流速,增大了流速的垂向梯度,对底部悬浮泥沙的向上扩散有明显的抑制作用。因此,温跃层的存在是造成夏季表层悬沙浓度最低的重要原因。关键词 ROMS 模型;渤海;温跃层
4、;表层悬沙浓度;季节变化 中图分类号 P734.2+3 doi:10.11693/hyhz20221200322 渤海作为我国北部典型的半封闭陆架海,是我国最早进行悬沙分布和输运研究的海域。前人首先进行了大量的悬沙浓度现场观测,得到了渤海悬沙浓度的空间分布特征及其季节变化规律(秦蕴珊等,1982;中国科学院海洋研究所海洋地质研究室,1985;Jiang et al,2004)。通过与渤海水动力环境相结合,能够大体明确悬沙输运的模式和机理。早期应用的数值模型不论是二维模型(Fang et al,1985;Zhang et al,1986;董礼先等,1989)还是三维模型(Feng,1987;Gr
5、aber 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),数值模拟结果揭示了渤海水动力环
6、境的基本特征。Xia 等(2019)利用高分辨率波浪-潮汐-环流耦合模式模拟渤海夏季的环流,发现辽东湾西部至黄河口近岸海域沿温度锋面存在明显的密度流,在 40N 以南海域环流呈双环流结构。马伟伟等(2016)利用 ROMS(regional ocean modeling system)海洋模式对渤海冬季流场进行高分辨率的数值模拟,发现冬季深度平均1016 海 洋 与 湖 沼 54卷 流场在辽东湾为顺时针环流,且在渤海中央表现为顺时针的环流结构。相比于研究较为成熟的渤海水动力特征,渤海悬沙的研究大都集中在黄河口,渤海海峡等热点区域(Bi et al,2011;Yang et al,2011;Wa
7、ng et al,2020)。受现场观测范围的限制,整个渤海海域的悬沙分布研究一般通过卫星遥感数据的反演得到(庞重光等,2013,2014;周舟等,2017;Zhao et al,2022)。庞重光等(2014)通过分析长时间序列悬沙浓度和风场遥感反演资料,指出渤海不同海域表层悬沙浓度差别较大,高浓度集中在辽东湾、渤海湾和莱州湾,低浓度主要位于秦皇岛海域以及渤海中部等。庞重光等(2013)基于 SeaWiFs 反演的渤海表层悬沙浓度资料,通过经验正交函数分解,展示了渤海表层悬浮泥沙的空间分布特征及其随时间的变化;其中第一空间模态显示了渤海表层悬沙的季节变化特征。周舟等(2017)利用20032
8、014 年的 MODIS(Moderate Resolution Imaging Spectroradiometer)卫星遥感数据,对渤海表层悬浮物浓度的空间分布形态进行反演,结果显示表层悬浮物浓度在时间尺度上有明显的季节变化,且渤海表层悬浮物从 20032014 年呈现下降的趋势。近年来,相对于卫星遥感资料反演悬沙浓度,利用数值模式对整个渤海海域悬沙分布的研究相对较少(Lu et al,2011;Liu et al,2020;Wang et al,2020)。利用数值模拟方法不仅可以再现渤海悬沙浓度分布的时空特征,同时数值模拟也是揭示悬沙浓度分布的动力成因的有效手段。本文基于 ROMS 海洋
9、模式,对渤海的水动力环境和悬沙分布进行气候态数值模拟,重点展示渤海海域的冬夏季流场特征、悬沙分布,并对悬沙分布及其影响机制进行讨论。1 研究方法 本文所采用的 ROMS 模式是基于三维非线性斜压原始方程的开源海洋模型,模拟范围为 3641N,117123E,如图 1 所示。地形数据主要采用 ETOPO1数据,水平方向上采用正交经纬网格,网格数为7361,水平分辨率为 5,垂向网格采用 坐标,共分为 20 层,最小水深设置为 1 m。模型计算采用冷启动,初始水位和流速均设为 0,初始场的温盐数据来自WOA13(World Ocean Atlas 2013),在初始场中将沉积物分为 3 类,分别为
10、黏土、粉砂和砂,三种沉积物的占比分布来自 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)单独计算的有效波高、波向、波长以及波周期等波浪参数作为外在强
11、迫添加到 ROMS 中。在边界强迫中,开边界水位采用 Chapman 边界条件,正压流速采用Flather 边界条件,三维斜压流速及温盐场均采用 Rad边界条件。在本模型中,由于在渤海入海的河流中,图 1 模拟区域的水深和站位(a)及水平网格(b)Fig.1 Water depth and stations(a)and horizontal grid(b)in the model simulated region 4 期 张德强等:渤海表层悬沙分布季节变化特征与控制因子的数值研究 1017 辽河、滦河以及海河流量远小于黄河,因此仅考虑黄河径流及悬浮物的输运作用。对于黄河径流量和悬浮物的数据来自
12、于利津水文站近 30 a 的气候态数据,且黄河中的黏土、粉砂和砂的占比分别为 20%、70%和 10%(Wang et al,2020)。本文用气候态数据将模式运行了 2 a,取第 2 年的 12 个月份的月平均值来表征每个月的特征。本文采用了 2018 年 7 月在渤海进行观测的大面站温盐深(conductivity temperature depth,CTD)数据与模拟的气候态的月平均值进行对比,同时采用 2017 年 23 月与2018年8月在秦皇岛海域的实测数据,对其求10 d的平均值,与模拟 10 d 的平均值进行对比。2 模式结果验证 2.1 潮位对比 对于潮汐的准确模拟是正确模拟
13、渤海水动力环境的前提。取第二年 1 月份的水位数据,采用调和分析方法分离出潮位的调和常数。由于渤海的主要半日分潮和全日分潮分别为 M2和 K1,因此本文仅对 M2和 K1进行分析。图 2c 和 2d 分别为 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 为模式模拟的同潮图;虚线表示振幅(单位:c
14、m),实线表示迟角(单位:)通过与 Fang(1986)的渤海图集中同潮图进行对比,可以发现两者之间分潮 M2与 K1的无潮点位置,振幅以及相位的分布均较为一致。通过对模拟的调和常数与 7 个渤海验潮站实测的 M2和 K1分潮的调和常数进行定量对比,M2分潮振幅误差的绝对平均为8.2 cm,相对误差的绝对平均为 11.8%;迟角误差的绝对平均为 14.3,相对误差的绝对平均为 8.7%。K1分潮振幅误差的绝对平均为 0.7 cm,相对误差的绝对1018 海 洋 与 湖 沼 54卷 平均为 4.6%;迟角误差的绝对平均为 22.3,相对误差的绝对误差为 19.4%。以上结果表明,该模型对渤海潮汐
15、的模拟是准确的。2.2 流场对比 图 3a 和图 3b 分别是渤海冬季(1 月)和夏季(7 月)的气候态垂向平均流场,从模拟结果来看,在渤海海峡终年呈现“北进南出”的结构;冬季沿岸线流动的海流最强,辽东湾以及渤海中部均呈现顺时针环流结构,渤海湾和莱州湾海流逆时针流动,整体海流流速较小,一般都在 0.05 m/s 以下,这与已有的数值模拟结果较为类似(毕聪聪等,2015;马伟伟等,2016;Jiang et al,2020;Liu et al,2021)。夏季环流明显强于冬季,流速很多在 0.050.1 m/s。夏季,20 m 等深线以深的渤海中部,环流呈现“双环”结构,即外部的逆时针环包围着内
16、部的顺时针环,这与 Xia 等(2019),Liu等(2021)的研究结果一致。本文的气候态模拟结果可以再现渤海环流的基本结构及其季节变化。3 结果与讨论 3.1 渤海温跃层的季节变化 为了研究渤海温跃层的季节性变化,将模式模拟的气候态结果与实测数据进行对比。实测数据分为两部分,第一部分是利用在S1站点(图1a黑色空心圆所示)于 2017 年 2 月 15 日3 月 9 日(冬季),2019 年 8月 1529 日(夏季)进行的海床基平台观测数据,并对其温度取 10 d 平均,与模拟的气候态温度进行对比,对比结果如图 4a 所示;第二部分是利用 2018 年 7 月对渤海进行的大面站观测数据与
17、模式模拟结果(月平均)进行对比,如图 5 所示,站点 B29B30,B53B59(红色空心圆)位置如图 1a 所示。图 3 渤海 1 月(a)和 7 月(b)气候态的深度平均流场 Fig.3 The climatological depth-averaged current fields in the Bohai Sea in January(a)and July(b),respectively 从图 4a 可以看出,模式对于 S1 站点季节性变化的模拟准确。图 4b 中的 S1 站点 CTD 实测数据在夏季的垂向分布显示,表层温度为 27 C,底层温度为24.5 C,没有明显的温跃层存在,这
18、也与图 4a 中的模拟结果相对应。图 4a 的模拟温度显示底层温度(Layer 1,Layer 2)与表层温度(Layer 20)相差无几。图 5 显示,在渤海水深较深的海域,夏季存在明显的温跃层。由于大面站观测获得的温度是瞬时值,且这些站点的温度大都是白天所测,故测得的温度普遍高于模式模拟的温度,且温跃层以上温度更高,垂向变化更大。而模式计算的温度是 7 月的平均值,温度略低且垂向变化范围较小。此外,模式模拟的盐度在 3132,而 B29B30,B53B59 站点所测的盐度大都也在 3132,进一步验证了模式模拟的准确性,故可以用模式模拟的结果来分析渤海温度的季节性变化。由于温跃层的存在使得
19、表层温度与底层温度差距较大,因此采用表层温度与底层温度的差值来体现温跃层的强度大小与位置分布,其 112 月份的空间分布如图 6 所示。图 6 显示,温跃层在 4 月开始形成,7 月达到最强,在 9 月之后,开始慢慢消退,10 月已经基本无温跃层的存在,这与屠金钊(1992)的研究结果较为一致,且在温跃层强盛期(68 月)表底层温差的分布与等深线的分布较为一致。4 期 张德强等:渤海表层悬沙分布季节变化特征与控制因子的数值研究 1019 图 4 S1 站点模式模拟温度(10 d 平均)底层(第一层,第二层)、表层(二十层)与观测温度(10 d 平均)的时间序列(a)以及 S1站点 8 月垂向温
20、度分布(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 s
21、imulated temperature and observed temperature 1020 海 洋 与 湖 沼 54卷 图 6 表层温度与底层温度差值的空间分布 Fig.6 Spatial distribution of temperature difference between the surface layer and the bottom layer 3.2 渤海悬沙分布的季节变化 本文参考 Bian等(2013)的泥沙参数方案,对泥沙各级粒径参数的相关设置见表 1。数值模拟的渤海表层悬沙浓度的季节性分布如图 7 所示,渤海表层的悬沙分布具有明显的季节变化。整体上讲,渤海的悬
22、沙浓度在冬季(122 月)最大,秋季(911 月)次之,春季(35 月)再次之,夏季(68 月)最小,且不论在哪个季节,悬沙浓度在辽东湾、渤海湾以及莱州湾均相对较高,大都在 70 mg/L 以上,悬沙的时空分布特征与庞重光等(2013)利用遥感资料对渤海表层悬沙的反演结果基本一致。图 7 还显示,秦皇岛海域悬沙浓度显著低于其他海域,常年偏低,其表层悬沙浓度基本都维持在 20 mg/L 以下,这一现象也与前人对渤海悬沙浓度的反演结果类似(Jiang et al,2004;Wang et al,2014;王震等,2016;Zhao et al,2022)。通过对模式模拟的表层悬沙浓度和卫星遥感反演
23、的海表悬沙浓度进行对 4 期 张德强等:渤海表层悬沙分布季节变化特征与控制因子的数值研究 1021 图 7 渤海表层悬沙浓度不同季节的空间分布 Fig.7 Spatial distribution of surface suspended sediment concentration in the Bohai Sea in different seasons 表 1 ROMS 泥沙模块中各级粒级泥沙参数 Tab.1 Sediment parameters in the ROMS Sediment Model 泥沙类型 泥沙粒径/mm 沉降速度/(mm/s)临界切应力/(N/m2)侵蚀速率/kg
24、/(m2s)密度/(kg/m3)疏松性 细砂 0.200 15.3 0.16 2.5105 2 650 0.4 粉砂 0.015 0.4 0.10 2.5105 2 650 0.4 黏土 0.004 0.1 0.05 2.5105 2 650 0.4 比,如图 8 所示(位置点 AF 如图 1a 所示),发现在渤海的不同位置,如典型的高、低悬沙浓度点,模式均能够模拟出表层悬沙浓度的季节变化特征,且大致与卫星遥感反演的悬沙浓度值一致,冬季浓度最高,夏季浓度最低。3.3 渤海底部切应力的季节变化 泥沙模型分别计算了时间步长为 3 min 的流致切应力、波致切应力以及波流相互作用产生的最大切应力。通
25、过对剪切应力求季节平均,得到不同季节剪切应力的大小,波流相互作用切应力、流致切应力和波致切应力分别如图 911 所示。从整体来看,由波流相互作用引起的最大剪切应力在冬季最大、秋季次之、春季和夏季相对较小。流致切应力与波流切应力几乎相当,波致切应力仅仅在沿岸的浅水区域(20 m等深线以浅)发挥一定的作用,其中在辽东湾和莱州湾较为显著,而在深水区,波致切应力与流致切应力相比,几乎可以忽略不计,这与 Wang 等(2020)的研究较为一致。从波流切应力的季节分布来看,辽东湾东部的波流切应力在冬季明显高于夏季,冬季波流切应力基本在 0.2 N/m2以上,夏季波流切应力主要在0.120.20 N/m2。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 渤海 表层 分布 季节 变化 特征 控制 因子 数值 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。