检索
E-mail
RSS

作者投稿 主编审稿 专家审稿 远程编辑

热带海洋学报 >

2021 , Vol. 40 >Issue 2: 7 - 16

DOI: https://doi.org/10.11978/2020057

海洋水文学

莫桑比克海峡及其邻近海区正压潮流数值模拟与能量收支分析*

  • 张华 , 1, 4 ,
  • 温茜茜 1 ,
  • 彭世球 , 1, 2, 3, 4
展开
  • 1.热带海洋环境国家重点实验室(中国科学院南海海洋研究所), 广东 广州 510301
  • 2.南方海洋科学与工程广东实验室, 广东 广州 511458
  • 3.中国科学院南海生态环境工程研究院, 广东 广州 510301
  • 4.中国科学院大学, 北京 100049
彭世球。email:

张华(1993—), 男, 广东省茂名市人, 博士研究生, 研究方向为潮致混合。email:

Copy editor: 殷波

收稿日期: 2020-05-27

  要求修回日期: 2020-07-20

  网络出版日期: 2020-07-21

基金资助

国家自然科学基金项目(41931182)

国家自然科学基金项目(41521005)

国家自然科学基金项目(41676016)

广东省重点工程(2019BT2H594)

南方海洋科学与工程广东实验室(广州)引进人才队伍重点专项(GML2019ZD0303)

南方海洋科学与工程广东实验室(广州)引进人才队伍重点专项(GML2019ZD0304)

中国科学院重点部署项目(ZDRW-XH-2019-2)

中国科学院南海生态环境工程创新研究院自主部署项目(ISEE2018PY05)

热带海洋学国家重点实验室独立研究项目(LTOZZ1902)

热带海洋学国家重点实验室独立研究项目(LTOZZ1802)

版权

版权所有，未经授权，不得转载、摘编本刊文章，不得使用本刊的版式设计。
收起

Numerical simulation of barotropic tides in Mozambique Strait and its adjacent coastal area and energy budget analysis*

  • ZHANG Hua , 1, 4 ,
  • WEN Xixi 1 ,
  • PENG Shiqiu , 1, 2, 3, 4
Expand
  • 1. State Key Laboratory of Tropical Oceanography (South China Sea Institute of Oceanology, Chinese Academy of Sciences), Guangzhou 510301, China
  • 2. Southern Marine Science and Engineering Guangdong Laboratory, Guangzhou 511458, China
  • 3. Institution of South China Sea Ecology and Environmental Engineering, Chinese Academy of Sciences, Guangzhou 510301, China
  • 4. University of Chinese Academy of Sciences, Beijing 100049, China
PENG Shiqiu. email:

Copy editor: YIN Bo

Received date: 2020-05-27

  Request revised date: 2020-07-20

  Online published: 2020-07-21

Supported by

National Natural Science Foundation of China(41931182)

National Natural Science Foundation of China(41521005)

National Natural Science Foundation of China(41676016)

Guangdong Key Project(2019BT2H594)

Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou)(GML2019ZD0303)

Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou)(GML2019ZD0304)

Key Deployment Projects of the Chinese Academy of Sciences(ZDRW-XH-2019-2)

Institution of South China Sea Ecology and Environmental Engineering, Chinese Academy of Sciences(ISEE2018PY05)

Independent Research Project Program of State Key Laboratory of Tropical Oceanography under contract(LTOZZ1902)

Independent Research Project Program of State Key Laboratory of Tropical Oceanography under contract(LTOZZ1802)

Copyright

Copyright reserved © 2021.
Fold

摘要

莫桑比克海峡及其邻近海区是全球海洋潮流和潮能耗散最强的海区之一。文章利用高分辨率通用环流模式对该海区的正压潮流进行模拟, 并对该海区潮能通量和潮能耗散特征进行分析。结果表明, 莫桑比克海峡及其邻近海区的潮波主要是半日分潮占主导地位, 全日分潮可忽略不计, M2分潮形成1个左旋潮波系统和1个右旋潮波系统, S2分潮形成1个左旋潮波系统。莫桑比克海峡和马达加斯加岛南部等绝大数区域的M2和S2半日潮流是逆时针旋转, 在马达加斯加岛顶部等局部区域是顺时针旋转, 而且在海峡通道等复杂地形处潮流流速量级较大。潮能通量矢量主要来自东边界, 大部分潮能通量沿马达加斯岛北部传入莫桑比克海峡区域, 其中经过马达加斯加岛北部和进入莫桑比克海峡的M2 (S2)分潮的潮能通量分别为156.86GW (40.53GW)和148.07GW (36.05GW), S2分潮潮能通量的量级大约为M2分潮的1/5~1/4。底摩擦耗散主要发生莫桑比克海峡和马达加斯加岛南北部, 其中莫桑比克海峡M2 (S2)分潮的底摩擦耗散为1.762GW (0.460GW), 占其底部总耗散的43.74% (39.72%)。

关键词： 莫桑比克海峡; 潮能通量; 潮能耗散; 数值模拟; 正压潮

本文引用格式

张华 , 温茜茜 , 彭世球 . 莫桑比克海峡及其邻近海区正压潮流数值模拟与能量收支分析*[J]. 热带海洋学报, 2021 , 40(2) : 7 -16 . DOI: 10.11978/2020057

Abstract

The Mozambique Strait and its adjacent coastal area are one of the areas with the strongest tidal current and energy dissipation in the world. In this study, the high-resolution MIT general circulation model (MITgcm) is used to simulate the barotropic current in this area, and the characteristics of tidal energy flux and dissipation are analyzed. The results show that the tidal waves in the Mozambique Strait and its adjacent coastal area are dominated by the semidiurnal tide while the diurnal tide can be ignored. The M2 barotropic tide forms an anticlockwise rotary tidal system and a clockwise rotary system, and the S2 barotropic tide forms an anticlockwise rotary tidal system only. The M2 and S2 barotropic tidal currents in the Mozambique Strait and the southern part of Madagascar are counter-clockwise, and those north of Madagascar and other local areas are clockwise. The current velocity is large in the strait and over rough topography. The tidal energy mainly divergences at the east boundary while most of that convergences at the Mozambique Strait along the northern part of Madagascar Island. The tidal energy flux of M2 (S2) barotropic tide passing through the northern part of Madagascar Island and the Mozambique Strait is 156.86 GW (40.53 GW) and 148.07 GW (36.05 GW), respectively. The dissipation induced by bottom friction mainly occurs in the Mozambique Strait and the south and north of Madagascar Island. The bottom friction dissipation of M2 (S2) barotropic tide is 1.762 GW (0.460 GW) in the Mozambique Strait, accounting for 43.74% (39.72%) of the total dissipation in the water column.

Key words： Mozambique Strait; tidal energy flux; tidal energy dissipation; numerical simulation; barotropic tides

*感谢南海海洋研究所的高性能计算中心为本文数值模拟提供的技术支持。
风和潮汐是驱动海洋内部混合的主要机械能来源, 其中表层海水混合所需能量主要由风提供, 中深层海水的跨等密面混合所需能量主要由潮汐提供。据估计, 通过潮汐传入海洋的能量约为3.5TW (Munk et al, 1998), 绝大部分潮能耗散在边缘海的底边界层(Jeffreys, 1921), 对海水的温度、盐度和密度有重要影响作用。潮汐能通量是指单位时间内通过自海底至海面单位密度宽度断面的潮能(方国洪 等, 1994; 朱怀鑫 等, 2018), 是表征潮汐强弱的基本指标(佟景全 等, 2010; 朱学明 等, 2012, 2014)。对潮汐能通量和潮汐能量耗散的早期研究是自20世纪初开始的。Taylor (1919)首次通过潮能底边界层(Bottom Boundary Layer, BBL)耗散公式估算了爱尔兰海潮能底摩擦耗散, Jeffreys (1921)利用BBL耗散公式估算了全球的潮能耗散。Zhong等(2006)通过潮能低摩擦耗散公式计算了切萨皮克湾的地边界耗散。张学庆等(2016)由耗散公式估算了辽河口的低摩擦耗散。李蔷等(2018)通过耗散公式计算了白令海峡及其邻近海区底边界耗散。周相乾等(2019)利用底摩擦耗散公式估算了南极布兰斯菲尔德海峡的潮能耗散。
莫桑比克海峡位于西南印度洋边缘, 非洲大陆与马达加斯加岛之间, 海峡及周边海域的拥有复杂的海底地形和丰富的环流系统。印度南赤道流(SEC)在到达马达加斯加岛东北时, 其西北分支(NEMC)绕过琥珀角后, 部分海水以涡旋的形式南下进入莫桑比克海峡, 通过海峡后汇入阿古拉斯(Agulhas)海流系统(Beal et al, 2011; Peng et al, 2015); 其西南分支(SEMC)则沿着马达加斯加岛东部沿岸南下, 在岛屿南部汇入Agulhas海流系统(图1)。马达加斯加岛南北部海区以及莫桑比克海峡北部窄口处的地形十分复杂, 有大量海山、海岭和山脊, 使得莫桑比克海峡及其邻近海区是全球海洋潮能耗散最强的地区之一(Egbert et al, 2000, 2001)。研究该海区的潮流、潮能通量和潮能耗散特征, 对于深入认识莫桑比克海峡及其邻近海区的涡旋和局地环流以及印度洋和大西洋之间的水交换过程具有重要意义。
以往已有学者对莫桑比克海峡及邻近海区进行了一些基础的物理海洋学研究, 获得了对该海区的水团性质、环流特征等方面的基本认识(Swart et al, 2010; Ridderinkhof et al, 2010; Breitzke et al, 2017)。然而, 由于马达加斯加岛附近海区潮汐潮流观测资料仍然十分匮乏, 以往对莫桑比克海峡及其邻近海的潮流、潮能通量和潮能耗散的定量分析研究还比较少。为此, 本研究基于高分辨率通用环流模式(MITgcm)的模拟结果, 对莫桑比克海峡及邻近海区的正压潮流、潮能通量和潮能耗散特征进行系统的分析与研究。

1 模式设置与分析方法

1.1 模式简介

MITgcm (MIT General Circulation Model)是麻省理工学院开发的大气-海洋通用环流模式(Marshall et al, 1997)。该模式在海洋领域得到广泛的应用, 具有静力近似、准静力近似和非静力近似的模拟功能, 模拟尺度几乎涵盖了所有的海洋动力过程, 对从几米的对流扩散过程到上千米的大洋环流等过程都有良好的模拟效果。该模式采用有限体积法、正交曲线网格和网格切削等技术对模式方程进行离散, 能够在合理地处理不规则边界和地形的同时, 保持模式通量守恒, 从而保证模式结果不产生漂移(刘谊 等, 2018)。由于MITgcm模式具有上述优势功能, 目前已经被广泛应用于内波和内潮破碎导致混合的模拟中(Buijsman et al, 2012, 2014; Wang et al, 2016)。

1.2 模式设置

本研究选取的模拟区域为包括整个马达加斯加岛以及西南印度洋部分海区(25°—57°E, 34°S—0), 如图1所示。模式水平分辨率为1/24°×1/24°, 垂向上不等间距地分为 67 层, 层厚自上而下的变化范围是5~100m, 采用z垂直坐标系。研究区域地形采用美国国家地球物理数据中心(NGDC)提供的高分辨率全球水深地形数据(ETOPO1)(Amante et al, 2009)。设置均一分布的初始温盐场, 均一分布场分别为常数值20℃和35‰。边界采用正压潮潮流的法向流速驱动, 潮流的振幅和迟角数据来自俄勒冈州立大学的海洋潮汐反演模式系统(Oregon State University Tidal Data Inversion Software, OTIS) (Egbert et al, 2002), 并选取包含莫桑比克海峡的印度洋区域模式解, 水平分辨率为 1/12°×1/12°。此外, 为了减弱边界的虚假反射, 在每个开边界处添加了宽度为0.5°的海绵区。
图1 莫桑比克海峡及其邻近海区地形分布、验潮站站位分布和海流示意图

黑色三角形和数字表示站位和站位号; 绿色带箭头实线表示海流及其方向; 莫桑比克海峡的绿色虚线箭头显示了涡流传播的影响。该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1665的标准地图制作, 底图无修改

Fig. 1 Bathymetry and tidal station distributions used. Also in the diagram are currents in Mozambique Strait and its adjacent coastal area

The stations are represented by black triangles; The currents are represented by green solid lines; Dashed arrow in the Mozambique Strait shows the effect of eddy propagation

本研究采用MITgcm静力版本, 且为了简单起见, 不添加湍流封闭方案。在整个模拟的区域范围内, 模式的水平扩散系数Kh和水平黏性系数Ah设置为 1×10-2m2·s-1, 垂向扩散系数Kv和垂向黏性系数Av设置为 1×10-5m2·s-1, 底摩擦系数为Cd=2.5× 10-3。积分时间步长为90s, 从静止开始积分47d。首先在开边界分别加入8个主要分潮的混合正压潮法向流速进行驱动, 再使用单个主要正压分潮法向流速进行驱动, 基于模拟结果对潮能通量和潮能耗散特征进行分析研究。

1.3 分析方法

1.3.1 潮流椭圆要素
在研究潮流结构的基本特征时, 常常将分潮流速E(t)分解成向东和向北的分量, 分别记作:
${{u}_{\text{s}}}={{U}_{\text{s}}}\cos (\sigma t-{{\xi }_{\text{s}}})$
${{v}_{\text{s}}}={{V}_{\text{s}}}\cos (\sigma t-{{\eta }_{\text{s}}})$
式中: us, vs分别为表层潮流向东和向北的分量(单位: m·s-1); Us, Vs, ζs, ηs分别为表层潮流向东和向北分量振幅(单位: cm)和相位(单位: °); t为时间(单位: s)。这个分潮流可以用复数的形式表示:
$E(t)=a{{\text{e}}^{-i\sigma t}}+b{{\text{e}}^{i\sigma t}}$
$a=A{{\text{e}}^{i\alpha }},b=B{{\text{e}}^{i\beta }}$
式中: i为虚数单位; A和B是正实数; α和β是角度(单位: °)。潮流椭圆要素可写为如下的简洁形式(夏综万, 1993):
$W=A+B$
$w=kW=B-A$
$H=\frac{\alpha +\beta }{2}$
$\theta =\frac{\alpha -\beta }{2}$
式中: W为最大潮流(单位: m·s-1); w为最小潮流(单位: m·s-1); H为最大潮流方向(单位: °); θ以角度表示最大潮流发生时间(单位: °), 每个周期(T)对应360°。k是椭圆旋转率, 并规定分潮流矢逆时针旋转k为正, 顺时针旋转k为负。其中,
$A=\frac{1}{2}{{\left[ U_{\text{s}}^{2}+V_{\text{s}}^{2}+2{{U}_{\text{s}}}{{V}_{\text{s}}}\sin ({{\xi }_{\text{s}}}-{{\eta }_{\text{s}}}) \right]}^{\frac{1}{2}}}$
$B=\frac{1}{2}{{\left[ U_{\text{s}}^{2}+V_{\text{s}}^{2}-2{{U}_{\text{s}}}{{V}_{\text{s}}}\sin ({{\xi }_{\text{s}}}-{{\eta }_{\text{s}}}) \right]}^{\frac{1}{2}}}$
$\alpha ={{\tan }^{-1}}\frac{{{U}_{\operatorname{s}}}\sin {{\xi }_{\operatorname{s}}}+{{V}_{\operatorname{s}}}\cos {{\eta }_{\operatorname{s}}}}{{{U}_{\operatorname{s}}}\cos {{\xi }_{\operatorname{s}}}-{{V}_{\operatorname{s}}}\sin {{\eta }_{\operatorname{s}}}}$
$\beta ={{\tan }^{-1}}\frac{-{{U}_{\operatorname{s}}}\sin {{\xi }_{\operatorname{s}}}+{{V}_{\operatorname{s}}}\cos {{\eta }_{\operatorname{s}}}}{{{U}_{\operatorname{s}}}\cos {{\xi }_{\operatorname{s}}}+{{V}_{\operatorname{s}}}\sin {{\eta }_{\operatorname{s}}}}$
1.3.2 潮能通量
本文采用Gill (1982)方法估算潮能通量。根据潮能量收支守恒关系, 在一个潮周期内, 通过任何海区水平方向的潮能通量与该海区的潮能耗散总量相平衡。由于对流项对能量通量的贡献通常比正压项小两个数量级, 因此如果只考虑由正压潮汐运动所产生的压力做功, 单位时间内通过单位宽度水柱的潮能通量可以进一步简化为(Zhong et al, 2006; 佟景全 等, 2010; 朱学明 等, 2012; Wang et al, 2016; 孟云 等, 2019; 陈元杰, 2020):
${{F}_{x}}={{\rho }_{0}}gh\langle {{u}_{\text{bt}}}\eta \rangle =\frac{1}{2}{{\rho }_{0}}ghHU\ \cos ({{\varphi }_{u}}-{{\varphi }_{\eta }})$
${{F}_{y}}={{\rho }_{0}}gh\langle {{v}_{\text{bt}}}\eta \rangle =\frac{1}{2}{{\rho }_{0}}ghHV\ \cos ({{\varphi }_{u}}-{{\varphi }_{\eta }})$
式中: Fx, Fy分别为潮能通量在x, y方向的分量(单位: kW·m-1); g为重力加速度(取9.8m·s-2); h为水深(单位: m); η为海表面高度(单位: m); H, U, V分别为垂向积分平均分潮潮位振幅及分潮潮流东、北分量的振幅(单位: cm); ${{\varphi }_{\eta }}$、${{\varphi }_{u}}$、${{\varphi }_{v}}$分别为相应的分潮潮位及潮流东、北分量的相位(单位: °); $\langle \cdot \rangle $代表在一个潮周期内进行平均。垂向积分平均流速ubt公式为:
${{u}_{\text{bt}}}(t)=\frac{1}{D}\int_{-d}^{n}{u\text{d}z}=\frac{1}{D}\bar{u}$
式中: u为潮流向东的分量(单位: m·s-1); $D=\eta +d$为自海底到海表面的厚度(单位: m); d为海水深度(单位: m); z为深度坐标(单位: m)。
1.3.3 潮能底摩擦耗散
海洋中的潮能耗散有很多种来源, 包括潮汐过程导致的垂向水体混合、底摩擦耗散和海底地形阻力耗散等。本研究采用底边界潮能耗散公式(Kang et al, 2012; 李蔷 等, 2018):
$d=\rho {{C}_{\text{d}}}\left| {{{\vec{u}}}_{\text{b}}} \right|(u{{u}_{\text{bt}}}+v{{v}_{\text{bt}}})$
式中: d为底边界潮能耗散(单位: W·m-2); Cd为底摩擦系数, 取0.0025; ρ为海水密度(单位: kg·m-3); ${{\vec{u}}_{\text{b}}}$为底层速度(单位: m·s-1); v为潮流向北的分量 (单位: m·s-1); vbt为垂向积分平均流速向北的分量(单位: m·s-1)。

2 模式结果验证

由于缺乏正压潮流流速的现场观测资料, 本文通过验证潮能通量来间接评判模式模拟的正压流流速的合理性。基于潮能通量公式(12), 估算研究区域M2、S2、K1和O1 4大分潮正压潮能通量, 分别为186.83GW、29.09GW、1.49GW和1.43GW, 其中两个半日分潮(M2和S2)正压潮能通量之和占了4大分潮(M2、S2、K1和O1)总正压潮能通量的98.66%。与利用海洋潮汐反演模式系统(OSU Tidal Data Inversion Software, OTIS)数据(Egbert et al, 2002)计算所得的正压潮能通量相比, 基于本文模式结果计算的半日分潮M2和S2的正压潮能通量相对均方根误差分别为3.29%和5.21%, 表明模式能较好的模拟出正压潮流流速。
本文选取研究区域14个验潮站的水位观测资料与相应位置的模式模拟水位, 通过对两个占绝对优势的半日分潮(M2、S2)的调和分析来验证模式对该区域正压潮模拟的准确性。验潮站资料来源于夏威夷大学海平面中心(Caldwell et al, 2015), 站位分布如图1所示。表1给出各验潮站两个半日分潮(M2、S2)调和常数的观测与模拟对比。需要说明的是, 13、14和17号验潮站观测时长较短而且数据缺测过多, 所以在此忽略不计。从表1和图2可以看出, M2 (S2)分潮的模式计算结果与夏威夷大学海平面中心(University of Hawaii Sea Level Center, UHSLC)验潮站观测结果比较一致, 其振幅和迟角的均方根误差分别为10.32cm (7.45cm)和25.31° (26.63°), 相关系数分别为0.97 (0.99)和0.99 (0.98)。少数几个验潮站数据与模式计算结果相比偏差较大, 其主要原因可能是验潮站位于港口内部, 数值模式不能很好考虑局地效应。总体而言, 模式结果和观测结果之间具有较好的一致性, 为接下来进行潮能通量和潮能耗散特征分析提供可信的模式数据。
表1 模式结果与验潮站分析结果的对比

Tab. 1 Comparison between model results and tide station analysis results

编号 站名 南纬 东经 M2 S2
振幅/cm 迟角/(°) 振幅/cm 迟角/(°)
观测 模拟 观测 模拟 观测 模拟 观测 模拟
1 Mobasa 4°04′ 39°39′ 104.46 86.53 25.99 35.70 51.39 41.27 64.88 39.10
2 Dar Es Salaam 6°49′ 39°17′ 107.07 96.25 29.29 38.86 52.74 46.12 70.04 43.09
3 Praslin 4°21′ 55°46′ 35.71 40.92 11.61 18.10 15.77 17.27 52.36 24.53
4 Port Victoria 4°37′ 55°28′ 40.22 43.66 12.21 18.90 17.77 18.24 52.37 25.12
5 Pt. La Rue 4°40′ 55°32′ 39.68 43.67 14.18 18.60 17.68 18.14 51.78 24.79
6 Lamu 2°16′ 40°54′ 98.47 81.86 32.28 35.25 48.12 38.98 73.79 38.80
7 Nosy Be 13°24′ 48°18′ 108.42 101.50 22.30 37.28 56.41 48.02 62.77 40.21
8 Zanzibar 6°09′ 39°11′ 120.59 127.61 25.43 82.74 60.90 55.43 64.65 93.20
9 Dzaoudzi 12°47′ 45°16′ 104.82 91.62 27.16 36.37 53.27 43.47 67.52 38.62
10 Reunion 20°56′ 55°17′ 17.84 22.98 282.83 341.79 9.11 7.43 297.89 292.63
11 Durban 29°52′ 31°03′ 55.81 52.86 44.36 65.91 31.30 24.20 76.07 49.64
12 Port Elizabeth 33°58′ 25°38′ 52.41 38.86 43.72 71.91 27.17 18.11 69.86 40.47
15 Pemba 12°58′ 40°29′ 113.43 98.36 29.12 40.90 57.49 47.78 71.11 42.71
16 Nacala 14°28′ 40°41′ 102.26 99.85 28.55 42.81 59.67 48.82 66.39 43.84
r 0.97 0.99 0.99 0.98
RMSE 10.32 25.31 7.45 26.63

注: r为相关系数, RMSE为均方根误差

图2 正压潮模式模拟结果与验潮站观测数据的散点图

a. 分潮M2振幅; b. 分潮S2振幅; c. 分潮M2迟角; d. 分潮S2迟角(迟角以格林尼治时间为时间基准)。r为相关系数, RMSE为均方根误差

Fig. 2 Scatter diagrams of barotropic model results and tidal station data

The amplitudes of tides (a) M2 and (b) S2 (units: cm); phase of tides (c) M2 and (d) S2 (phase based on Greenwich mean time)

3 正压潮特征

3.1 潮汐类型

任何海区的潮汐变化中的包含有半日周期和日周期变化的震荡, 这两部分震荡的相对大小决定潮汐类型, 而潮汐类型主要是根据全日潮和半日潮振幅的比值来判别(王斌, 2010; 罗丹, 2015; 周相乾 等, 2019):
$A=\frac{({{H}_{\text{K}}}_{_{1}}+{{H}_{\text{O}}}_{_{1}})}{{{H}_{\text{M}}}_{_{2}}}$
式中: A为判别系数; ${{H}_{\text{K}}}_{_{1}}$、${{H}_{\text{O}}}_{_{1}}$和${{H}_{\text{M}}}_{_{2}}$分别为K1、O1和M2分潮的振幅(单位: cm)。判别依据为: 比值A在0~0.5为规则半日潮, 0.5~2.0为不规则半日潮, 2.0~4.0为不规则全日潮, A大于4.0为规则全日潮。
根据以上潮汐类型判别方法, 莫桑比克海峡及其邻近海区的潮汐类型主要是规则半日分潮(图3), 只在模式区域东北部附近海区分布着较小范围的不规则半日分潮。由于研究海区里半日分潮占绝对主导地位, 因此本文后面只对半日分潮的正压潮能及耗散特征进行分析。
图3 潮汐类型分布

图中色标值表示: 0~0.5为规则半日潮, 0.5~2.0为不规则半日潮, 2.0~4.0为不规则全日潮, > 4.0为规则全日潮。该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1665的标准地图制作, 底图无修改

Fig. 3 Distribution of tidal constituents

The color scale: 0 ~ 0.5 for regular semidiurnal tide, 0.5 ~ 2.0 for irregular semidiurnal tide, 2.0 ~ 4.0 for irregular diurnal tide, and > 4.0 for regular diurnal tide

3.2 潮汐特征

图4显示模拟得到的2个主要半日分潮M2和S2同潮图, 可以看到M2和S2分潮等振幅线分布大体一致, 其最大振幅都出现在莫桑比克海峡通道, 分别为156.84cm和79.74cm, 最小振幅在马达加斯加岛西南角。从迟角的分布可以看出M2分潮波从马达加斯加岛东南部进入研究区域, 根据Taylor (1922)的矩形海湾旋转潮波系统形成机制, M2分潮在马达加斯加岛东部沿岸呈现对称分布, 入射潮波经海岸反射分为两支向南北传播: 北支向北马达加斯加岛沿岸传入莫桑比克海峡, 在马达加斯加岛北端形成一个左旋潮波系统; 南支沿马达加斯加岛南岸, 在马达加斯加岛南端形成一个退化的右旋潮波系统, 并在莫桑比克海峡形成潮波共振现象, 其中海峡通道中部及岛屿南端地形变浅, 摩擦作用使得入射潮波波长变短、反射潮波减弱(图4a)。S2分潮波从马达加斯加岛西北部的开阔海域进入, 在科氏力作用下向西南方向偏转传入莫桑比克海峡通道, 沿着马达加斯加岛形成一个左旋潮波系统, 同样由于海峡通道中部地形变浅导致潮波逐渐减弱(图4b)。
图4 莫桑比克海峡及其附近海区的同潮图

a. M2分潮的振幅和迟角; b. S2分潮的振幅和迟角。黑色实线表示振幅(单位: cm); 色柱表示迟角(以格林尼治时间为时间基准, 单位: °)。该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1665的标准地图制作, 底图无修改

Fig. 4 Co-tidal charts in the Mozambique Strait and its adjacent coastal area

(a) M2 and (b) S2 constituents. The black solid line indicates the amplitudes (units: cm) and the color indicates phase (based on Greenwich mean time)

3.3 潮流性质

根据公式(1)~(12), 将模式输出的流场数据进行调和分析, 得到2个主要半日分潮M2和S2的表层潮流椭圆分布, 如图5所示。其中潮流椭圆越扁表明旋转率越小, 往复式潮流越强; 旋转率越大, 旋转式潮流越强。可以看到M2分潮在莫桑比克海峡非洲沿岸、马达加斯加岛南部和东部等绝大部分海区以逆时针旋转为主, 仅在马达加斯加岛西部和顶部等海区存在顺时针旋转, 其潮流流速在莫桑比克海峡和马达加斯加岛南北端较大, 大约为0.2~1.0m·s-1 (图5a)。S2分潮潮流椭圆分布与M2大体一致, 其潮流流速整体上小于M2分潮的潮流流速, 大约为0.1~0.4m·s-1 (图5a)。从图5中还可以看出, 强的往复式潮流主要分布在海峡、水道等特殊地形。
图5 表层潮流椭圆分布

a. M2分潮; b. S2分潮。色标表示椭圆旋转率, 其中正值(红色)表示潮流椭圆逆时针旋转, 负值(蓝色)表示潮流椭圆顺时针旋转。该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1665的标准地图制作, 底图无修改

Fig. 5 Elliptical maps of surface tidal current

(a) M2 and (b) S2 constituents. The color indicates the rotation rate of ellipse. Positive values (red) indicate counterclockwise rotation of the tidal current ellipse, and negative values (blue) indicate clockwise rotation of the tidal current ellipse

3.4 潮能通量

图6给出了基于模式输出结果和OTIS数据计算的M2 和S2分潮的正压潮能通量矢量分布, 由于这两个半日分潮的正压潮能通量占了4大分潮总潮能通量的98.66%, 因此基本能代表潮波的传播情况。从图6a和6c中可见, 来自西印度洋的M2半日潮波由东边界进入模式研究区域, 在马达加斯加岛东部经海岸阻挡分为南北两支, 大部分能量沿马达加斯加岛东岸北上, 到达岛屿东北角海域后转向西传, 由于非洲大陆的阻挡, 绝大部分能量经过莫桑比克海峡向南传播, 在马达加斯加岛南端与沿着马达加斯加岛东部沿岸南下的南支潮能通量汇集。S2 半日潮波潮能通量矢量分布与M2分潮波潮能通量矢量分布大体一致, 主要差异在于马达加斯加岛东部沿岸南支潮能通量较弱, 经过莫桑比克海峡的北支潮能部分东传。以上基于模式输出结果计算的M2和S2半日分潮正压潮能通量分布与基于OTIS数据计算的M2和S2半日分潮正压潮能通量分布吻合得很好(图6b、6d), 整个研究区域的相对均方根误差分别为3.29%和5.21%。
图6 莫桑比克海峡及其附近海区M2分潮(a, b)和S2分潮(c, d)的潮能通量矢量分布

a, c: 模式结果; b, d: OTIS计算结果。箭头和断面附近的数值表示通过该断面的潮能通量, 正值表示纬向向东或经向向北, 负值则相反。A1、A2、A3和A4断面的位置分别为16°12′S、46°E、14°S、49°18′E

Fig. 6 Tidal energy flux in the Mozambique Strait and its adjacent coastal area

Model results (a) and OTIS (b) calculated tidal energy flux distributions of M2 constituent. Model results (c) and OTIS (d) calculated tidal energy flux distribution of S2 constituent. The value near the section indicates the tidal energy flux through the section. Positive value indicates eastward or northward transport, and negative value is the opposite. Sections A1, A2, A3, and A4 are located along 16°12′S, 46°E, 14°S, and 49°18′E, respectively

为了进一步分析研究区域的潮能耗散情况, 本研究选取了A1 (16°12′S)、A2 (46°E)、A3 (14°S)和A4 (49°18′E) 4个断面(图6a、6c), 计算两个主要半日分潮M2和S2通过它们的潮能通量值。结果显示, M2分潮通过A1、A2、A3和A4断面的潮能通量依次为148.07GW、36.77GW、69.1GW和156.86GW; S2分潮通过A1、A2、A3和A4断面的潮能通量分别为36.05GW、6.91GW、21.23GW2和40.53GW, 其量级大约为M2 分潮的1/5~1/4。M2 和S2分潮通过A1断面的潮能通量分别约占通过A4断面的94% 和89%, 表明经过A4断面的潮能通量绝大部分经过了A1进入了莫桑比克海峡, 而两者通过A2断面的潮能通量在方向上相反, 即M2 分潮潮能通量向西传, S2分潮潮能通量向东传。

3.5 潮能底摩擦耗散

图7为莫桑比克海峡及其邻近海区由底摩擦耗散公式(15)估算的M2 和S2分潮底摩擦耗散分布图。图7a和7b的空间分布基本一致, 底摩擦耗散主要发生莫桑比克海峡和马达加斯加岛南北部。为了研究计算海区各个海区的底摩擦耗散情况, 本研究选取了3个代表区块, 依次为马达加斯加岛北部区域(MADN)、莫桑比克海峡区域(MC)和马达加斯加岛南部区域(MADS), 分别计算整个研究区域和这3个子区域的底摩擦耗散量级大小(图7)。结果表明, 整个研究海区的M2和S2分潮的底摩擦耗散分别为4.028GW和1.158GW。3个子区域MADN、MC和MADS的M2 (S2)分潮底摩擦耗散分别为0.592GW (0.151GW)、1.762GW (0.460GW)和0.235GW (0.051GW), 约占整个研究区域底摩擦耗散的14.70% (13.04%)、43.74% (39.72%)和5.83% (4.40%), 其中由于莫桑比克海峡是狭长的海水通道且地形崎岖, 潮汐在复杂的水道中造成的急流与较浅较窄的地形相互作用(Donohue et al, 2003; 周相乾 等, 2019), 形成了强底摩擦耗散区域, 其底摩擦耗散接近整个海区的底摩擦耗散的一半。
图7 莫桑比克海峡及其附近海区潮能底摩擦耗散分布(采用对数log10色标)

a. M2分潮; b. S2分潮。红色方框为马达加斯加岛北部, 黑色方框为莫桑比克海峡, 黄色方框为马达加斯加岛南部。方框内数值表示该区域的潮能耗散值, 图内最上面的数字是表示整个研究区域的耗散值(单位: GW)。

Fig. 7 Tidal energy dissipation induced by bottom friction in the Mozambique Strait and its adjacent coastal area (logarithmic (log10) scale is used here)

(a) M2 and (b) S2 constituents. The red box is the north of Madagascar, the black box is the Mozambique Strait, and the yellow box is the south of Madagascar. The value in the box represents the tidal energy dissipation value in this area, the number in the top of figure represents the dissipation value of the whole study area (units: GW)

4 结论

本研究基于MITgcm模式建立了莫桑比克海峡及其邻近海区的高分辨率潮汐数值模式, 并在模式模拟结果的基础上, 对莫桑比克海峡及其邻近海区潮流、潮能通量和潮能耗散进行分析, 取得以下结论:
1) 莫桑比克海峡及其邻近海区的潮汐类型以规则半日潮为主, 且M2分潮占主导, S2分潮次之, 全日分潮不显著。
2) M2和S2半日分潮振幅的最大值出现在莫桑比克海峡通道, 分别是156.84cm和79.74cm, 最小值在马达加斯加岛西南角。M2分潮在马达加斯加岛南北两端分别形成右旋潮波系统和左旋潮波系统, S2分潮在马达加斯加岛出现了一个左旋潮波系统。
3) M2和S2分潮潮流椭圆分布相似, 在莫桑比克海峡和马达加斯加岛南部等绝大部分区域是逆时针旋转, 仅在马达加斯加岛顶部等局部区域是顺时针旋转, 而且潮流流速在莫桑比克海峡和马达加斯加岛南北部等复杂地形处最大, 达到1.0m·s-1。
4) 在潮能通量分析中, M2与S2半日分潮的正压潮能通量占4大分潮总正压潮能通量的98.66%。半日分潮正压潮能通量矢量分布情况表明莫桑比克海峡及其邻近海区的正压潮能通量主要从东部传入, 绝大部分潮能通量沿马达加斯岛北部传入莫桑比克海峡区域, 其中通过马达加斯加岛北部断面的M2 (S2)分潮的正压潮能通量为156.86GW (40.53GW), 通过莫桑比克海峡断面的M2 (S2)分潮的正压潮能通量为148.07GW (36.05GW), S2分潮潮能通量值大约为M2分潮的1/5~1/4。
5) 研究区域的底摩擦耗散主要发生莫桑比克海峡和马达加斯加岛南北部, 其中莫桑比克海峡的底摩擦耗散最强, M2 (S2)分潮的底摩擦耗散为1.762GW (0.460GW), 占整个研究区域的43.74% (39.72%)。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序
[1]
陈元杰, 程鹏, 2020. 东中国海潮能通量与耗散的数值模拟分析[J]. 厦门大学学报(自然科学版), 59(1):61-70.

CHEN YUANJIE, CHENG PENG, 2020. Numerical modeling study of tidal energy and dissipation in the East China Seas[J]. Journal of Xiamen University (Natural Science), 59(1):61-70 (in Chinese with English abstract).

[2]
方国洪, 曹德明, 黄企洲, 1994. 南海潮汐潮流的数值模拟[J]. 海洋学报, 16(4):1-12 (in Chinese).

[3]
李蔷, 高郭平, 安佰超, 等, 2018. 白令海峡及其邻近海域潮汐潮能数值模拟[J]. 极地研究, 30(1):1-13.

LI QIANG, GAO GUOPING, AN BAICHAO, et al, 2018. Numerical study of the Barotropic tide and tidal energy distribution in the Bering strait and adjacent sea areas[J]. Chinese Journal of Polar Research, 30(1):1-13 (in Chinese with English abstract).

[4]
刘谊, 王晓玮, 彭世球, 2018. 印度尼西亚海内潮生成及传播过程研究[J]. 热带海洋学报, 37(2):1-9.

LIU YI, WANG XIAOWEI, PENG SHIQIU, 2018. On the generation and propagation of internal tides in the Indonesian Seas[J]. Journal of Tropical Oceanography, 37(2):1-9 (in Chinese with English abstract).

[5]
罗丹, 2015. 渤、黄、东海潮汐潮流数值模拟研究[D]. 上海: 上海海洋大学.

LUO DAN, 2015. Numerical simulations of tide and tidal current in the Bohai Sea, Yellow Sea and East China Sea[D]. Shanghai: Shanghai Ocean University (in Chinese with English abstract).

[6]
孟云, 谢蓉, 2019. 渤海潮流、潮能通量和耗散的数值模拟[J]. 上海船舶运输科学研究所学报, 42(1):70-77.

MENG YUN, XIE RONG, 2019. Numerical simulation of tidal currents, energy fluxes and dissipation in the Bohai Sea[J]. Journal of Shanghai Ship and Shipping Research Institute, 42(1):70-77 (in Chinese with English abstract).

[7]
佟景全, 雷方辉, 毛庆文, 等, 2010. 不考虑局地引潮势的南海正压潮能通量与潮能耗散[J]. 热带海洋学报, 29(3):1-9.

TONG JINGQUAN, LEI FANGHUI, MAO QINGWEN et al, 2010. Tidal energy fluxes and dissipation in the South China Sea without considering tide-generated potential energy[J]. Journal of Tropical Oceanography, 29(3):1-9 (in Chinese with English abstract).

[8]
王斌, 2010. 乐清湾潮能和潮能通量的数值研究[D]. 杭州: 浙江大学.

WANG BIN, 2010. Numerical study of Tidal energy and Tidal energy flux in Yueqing Bay[D]. Hangzhou: Zhejiang University (in Chinese with English abstract).

[9]
夏综万, 1993. 介绍一种潮流椭圆要素的计算方法[J]. 海洋信息, (11):14 (in Chinese).

[10]
张学庆, 王兴, 刘睿, 等, 2016. 辽河口潮能通量与潮能耗散的数值研究[J]. 海洋环境科学, 35(1):20-26, 67.

ZHANG XUEQING, WANG XING, LIU RUI, et al, 2016. Numerical simulation of the tidal energy flux and tidal energy dissipation in Liaohe River Estuary[J]. Marine Environmental Science, 35(1):20-26, 67 (in Chinese with English abstract).

[11]
周相乾, 胡松, 张瑜, 2019. 南极布兰斯菲尔德海峡正压潮数值模拟[J]. 极地研究, 31(1):56-68.

ZHOU XIANGQIAN, HU SONG, ZHANG YU, 2019. Numerical simulation of barotropic tide in the Bransfield Strait, Antarctica[J]. Chinese Journal of Polar Research, 31(1):56-68 (in Chinese with English abstract).

[12]
朱怀鑫, 俎婷婷, 李健, 等, 2018. 基于高频地波雷达观测的粤西近海潮流潮能分析[J]. 热带海洋学报, 37(5):25-32.

ZHU HUAIXIN, ZU TINGTING, LI JIAN, et al, 2018. Analysis of tidal current and tidal energy on the west coast of Guangdong Province based on High Frequency Ground Wave Radar observation[J]. Journal of Tropical Oceanography, 37(5):25-32 (in Chinese with English abstract).

[13]
朱学明, 刘桂梅, 2012. 渤海、黄海、东海潮流、潮能通量与耗散的数值模拟研究[J]. 海洋与湖沼, 43(3):669-677.

ZHU XUEMING, LIU GUIMEI, 2012. Numerical study on the tidal currents, tidal energy fluxes and dissipation in the China Seas[J]. Oceanologia Et Limnologia Sinica, 43(3):669-677 (in Chinese with English abstract).

[14]
朱学明, 宋德海, 鲍献文, 等, 2014. 西北太平洋的潮能通量与耗散[J]. 热带海洋学报, 33(1):1-9.

ZHU XUEMING, SONG DEHAI, BAO XIANWEN et al, 2014. Tidal energy flux and dissipation in the Northwest Pacific[J]. Journal of Tropical Oceanography, 33(1):1-9 (in Chinese with English abstract).

[15]
AMANTE C, EAKINS B W, 2009. ETOPO1 1 Arc-minute global relief model: procedures, data sources and analysis[Z]. NOAA Technical Memorandum NESDIS NGDC-24. National Geophysical Data Center, NOAA.

[16]
BEAL L M, DE RUIJTER W P M, BIASTOCH A, et al, 2011. On the role of the Agulhas system in ocean circulation and climate[J]. Nature, 472(7344):429-436.

DOI PMID

[17]
BREITZKE M, WILES E, KROCKER R, et al, 2017. Seafloor morphology in the Mozambique Channel: evidence for long-term persistent bottom-current flow and deep-reaching eddy activity[J]. Marine Geophysical Researches, 38(3):241-269.

[18]
BUIJSMAN M C, LEGG S, KLYMAK J M, 2012. Double-ridge internal tide interference and its effect on dissipation in Luzon strait[J]. Journal of Physical Oceanography, 42(8):1337-1356.

[19]
BUIJSMAN M C, KLYMAK J M, LEGG S, et al, 2014. Three-dimensional double-ridge internal tide resonance in Luzon strait[J]. Journal of Physical Oceanography, 44(3):850-869.

[20]
CALDWELL P C, MERRIFIELD M A, THOMPSON P R, 2015. Sea level measured by tide gauges from global oceans - the Joint Archive for Sea Level holdings (NCEI Accession 0019568), Version 5.5[Z]. NOAA National Centers for Environmental Information, Dataset.

[21]
DONOHUE K A, TOOLE J M, 2003. A near-synoptic survey of the Southwest Indian Ocean[J]. Deep Sea Research Part Ⅱ: Topical Studies in Oceanography, 50(12-13):1893-1931.

[22]
EGBERT G D, RAY R D, 2000. Significant dissipation of tidal energy in the deep ocean inferred from satellite altimeter data[J]. Nature, 405(6788):775-778.

PMID

[23]
EGBERT G D, RAY R D, 2001. Estimates of M2 tidal energy dissipation from TOPEX/Poseidon altimeter data[J]. Journal of Geophysical Research: Oceans, 106(C10):22475-22502.

[24]
EGBERT G D, EROFEEVA S Y, 2002. Efficient inverse modeling of barotropic ocean tides[J]. Journal of Atmospheric & Oceanic Technology, 19(2):183-204.

[25]
GILL A E, 1982. Atmosphere-ocean dynamics[M]. London: Academic Press.

[26]
JEFFREYS H, SHAW W N, 1920. Tidal friction in shallow seas[J]. Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character 221: 239-264.

[27]
KANG DUJUAN, FRINGER O, 2012. Energetics of Barotropic and Baroclinic tides in the Monterey bay area[J]. Journal of Physical Oceanography, 42(2):272-290.

[28]
MARSHALL J, ADCROFT A, HILL C, et al, 1997. A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers[J]. Journal of Geophysical Research, 102(C3):5753-5766.

[29]
MUNK W, WUNSCH C, 1998. Abyssal recipes Ⅱ: Energetics of tidal and wind mixing[J]. Deep Sea Research Part Ⅰ: Oceanographic Research Papers, 45(12):1977-2010.

[30]
PENG SHIQIU, QIAN YUKUN, LUMPKIN R, et al, 2015. Characteristics of the near-surface currents in the Indian Ocean as deduced from satellite-tracked surface drifters. Part Ⅰ: Pseudo-Eulerian statistics[J]. Journal of Physical Oceanography, 45(2):441-458.

[31]
RIDDERINKHOF H, VAN DER WERF P M, ULLGREN J E, et al, 2010. Seasonal and interannual variability in the Mozambique Channel from moored current observations[J]. Journal of Geophysical Research, 115(C6):C06010.

[32]
SWART N C, LUTJEHARMS J R E, RIDDERINKHOF H, et al, 2010. Observed characteristics of Mozambique Channel eddies[J]. Journal of Geophysical Research: Oceans, 115(C9):C09006.

[33]
TAYLOR G I, 1919. Tidal Friction in the Irish Sea[J]. Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character, 96(678):330.

[34]
TAYLOR G I, 1922. Tidal oscillations in gulfs and rectangular basins[J]. Proceedings of the London Mathematical Society, 20(1):148-181.

[35]
WANG XIAOWEI, PENG SHIQIU, LIU ZHIYU, et al, 2016. Tidal mixing in the South China Sea: An estimate based on the internal tide energetics[J]. Journal of Physical Oceanography, 46(1):107-124.

[36]
ZHONG LIEJUN, LI MING, 2006. Tidal energy fluxes and dissipation in the Chesapeake Bay[J]. Continental Shelf Research, 26(6):752-770.

Options
文章导航

/

〈 〉