A numerical study of a tide with four harmonic constituents at the open boundary

The tidal current in Kagoshima Bay is simulated by the two-dimensional subdomain finite-element method. At the open boundary, the sea level is forced to be oscillatory with a linear combination of the four harmonic constituents (M2, S2, K1, O1) of the tide. A calculation having only M2 sea level at the open boundary is also performed and the harmonic constants of M2 are compared with those calculated by the four components.In the calculated velocity vector, a large difference between the two maxima or the two minima in one day appears when the vector is directed to the open sea. On the other hand, when it is directed into the bay, the two extreme velocities are almost equal. In tide-killer residual flow, the flow pattern is almost independent of the tidal age; however, the absolute value of the current velocity depends on the tidal age. The pattern is nearly the same as that calculated by giving only M2 sea level at the open boundary.In relation to the nonlinear interaction, the kinetic energy ratio KER is defined by the ratio of kinetic energy of the induced harmonic components to that of harmonic constituents given at the open boundary. KER is large where the tidal current is disturbed by obstacles,e.g. along the west coast at the mouth of the bay and in the southern sea of Mt. Sakurajima. The distribution of KER seems to be independent of the number of harmonic constituents (one or four) given at the open boundary. The difference of kinetic energy of the M2 tidal constituent between the two calculations.i.e. with the open boundary conditions of four components (M2, S2, K1, O1) and of a single component (M2), is found to be large where the current velocity is large.     

南海东北部内潮数值研究

导出二维三层（准三维）潮流与起伏地形相互作用而产生内潮波的非线性数学模型，并分别以Ｍ２分潮、Ｋ１分潮为策动期波，数值求解南海东北部大陆斜坡对策动潮波的响应。结果表明：台湾浅滩南面斜坡对Ｍ２分潮响应最为强烈，东沙群岛东面斜坡次之，其余海域无显著响应；对Ｋ１分潮的响应则以一统暗沙东南面斜坡最甚，其次为东沙群岛东南的斜坡，其余海域响应微弱。     

Numerical modeling of tide and tidal current in the kangjin bay,south sea,Korea

This study is based on a series of numerical modeling experiments to understand the tidal circulation in the Kangjin Bay (KB). The tidal circulation in the KB is mostly controlled by the inflow from two channels, Noryang and Daebang which introduce the open ocean water into the northern part of the KB with relatively strong tidal current, while in the southern part of the KB, shallowest region of the entire study area, weak tidal current prevails. The model prediction of the sea level agrees with observed records at skill scores exceeding 90 % in terms of the four major tidal constituents (M2, S2, K1, O1). However, the skill scores for the tidal current show relatively lower values of 87, 99, 59, 23 for the semi-major axes of the constituents, respectively. The tidal ellipse parameters in the KB are such that the semi-major axes of the ellipse for M2 range from 1.7 to 38.5 cm/s and those for S2 range from 0.5 to 14.4 cm/s. The orientations of the major-axes show parallel with the local isobath. The eccentricity values at various grid points of ellipses for M2 and S2 are very low with 0.2 and 0.06 on the average, respectively illustrating that the tidal current in the KB is strongly recti-linear. The magnitude of the tidal residual current speed in the KB is on the order of a few cm/s and its distribution pattern is very complex. One of the most prominent features is found to be the counter-clockwise eddy recirculation cell at the mouth of the Daebang Channel.     

基于FVCOM的南海北部海域潮汐潮流数值模拟

基于非结构三角形网格的FVCOM(finite-volume coastal ocean model )数值模型, 对南海北部海域的潮汐、潮流进行了精细化数值模拟研究, 并根据模拟结果详细分析了M2, S2, K1, O1 分潮的潮汐和潮流特征。研究结果表明: 神泉港到甲子港海域表现为正规全日潮性质, 珠江口附近海区潮汐以不正规半日潮为主, 其他海域主要表现为不规则全日潮; 陆架海域和深水海域主要表现为往复流, 陆架坡折区存在较强的旋转流, 陆架坡折区为不规则半日潮流和不规则全日潮流的分界线; 东沙群岛附近海域以不规则全日潮流为主, 旋转方向为顺时针; 整个海域的最大流速分布与等深线基本平行, 东沙群岛附近速度明显变大, 最大值出现在台湾浅滩附近, 最大值超过70 cm/s; 南海潮波系统以巴士海峡传入的大洋潮波为主, 分为三支潮流, 以不同的形式进出南海北部海域; 余流在台湾浅滩附近达到最大, 超过6 cm/s, 自南向北进入台湾海峡, 近岸余流自东向西沿岸流动。本研究在东沙群岛周边的模拟结果与前人基于实测资料的分析吻合较好, 并且由于采用了高精度的三角网格, 本文对东沙群岛周边海域的潮汐潮流结构和性质的刻画和分析是迄今为止较为精细的, 同时本研究还提高了对沿岸验潮站调和常数的模拟精度。     

基于FVCOM的廉州湾及周边海域三维潮汐潮流数值模拟(英文)

基于采用不规则三角网格和有限体积方法的FVCOM模式,建立廉州湾三维潮流数值模型来重现廉州湾及周边海域的潮位和潮流变化状况。根据模拟结果计算得到了较以往更为精细的廉州湾及周边海域K1、O1和M2分潮的同潮图,并计算了由此三个分潮引起的潮汐不对称的变化情况。K1和O1分潮在廉州湾外主要以驻波的形式存在,进入廉州湾后转化为前进波;M2分潮在廉州湾外主要以前进波的形式存在,进入廉州湾后前进波特征更为明显。K1和O1分潮流在廉州湾外以旋转流为主,在廉州湾内以往复流为主;M2分潮流在整个研究海域以往复流为主。由潮余流场的分布特点可以看出自南向北由外海进入廉州湾的潮余流,在冠头岭处分为两支,一支逆时针转向西,另一支被冠头岭阻挡在其南侧形成顺时针封闭环流。在廉州湾内部同时存在两个环流系统,湾顶的气旋式环流和口门处的反气旋式环流。     

渤黄海潮汐潮流精细化数值模拟及可视化预报系统

利用有限元海洋模式ADCIRC (Advanced Circulation Model),建立了高分辨率的渤黄海二维潮汐潮流数值模型,该模型以M2,S2,K1,O1等8个分潮的水位作为驱动,模拟出了该8个分潮的潮汐潮流调和常数;利用该调和常数预报的潮位和二维平均潮流与实测资料相比,符合较好;利用模拟得到的潮汐潮流调和常...     

黄海、渤海TOPEX/Poseidon高度计资料潮汐伴随同化

首先将大约10a的TOPEX/Poseidon(T/P)高度计资料沿星下轨迹点做潮汐调和分析,提取得到各分潮的调和常数,利用伴随同化方法,同化到二维非线性潮汐数值模式中,模拟了黄海、渤海区域M₂,S₂,O₁,K₁等4个潮汐分潮,并根据计算结果给出了各分潮的同潮图.将计算值与观测值的进行偏差统计,结果表明计算值与验潮站资料符合良好.研究过程中做了两类试验:一类试验是针对不同的参数进行优化,一类试验是针对不同的资料进行同化.第一类试验表明:将开边界条件和底摩擦系数同时作为模型优化的控制参数,其结果明显优于单独优化开边界条件;第二类试验表明:同时同化高度计资料与验潮站资料,比单独同化其中任一种资料,对模式计算结果都有较好的改进.研究结果表明,采用伴随同化方法,利用T/P高度计资料和验潮站资料作为同化数据能有效改进模拟结果,用来反演黄海、渤海的潮波系统是可行的.     

Distribution of 9 Principal Tidal Constituents and Tide Prediction for the East China Sea

By use of the hydrodynamic model,the harmonic constants of 8 principal tidal constituents(Q_1,O_1,P_1,K_1,N_2,M_2,S_2andK_2)are obtained for the East China Sea,and the harmonic constant ofS_a is calculated by two-dimensional interpolation.The calculated results agree well with the observed dataaround the sea.The harmonic constants can be used to predict the tide in the East China Sea.The cotidalcharts of the 9 tidal constituents reveal their distribution.     

全球大洋潮汐模式在北印度洋潮汐预报准确性的评估

为评估DTU10、TPXO8、GOT00.2和NAO.99b 4个全球大洋潮汐模式对北印度洋潮汐的预报能力,采用英国海洋资料中心提供的海区中部和沿岸站潮汐调和常数资料,检验了这些模式4个主要分潮(M_2、S_2、K_1、O_1)的准确度。它们的各分潮调和常数资料准确度都比较高,振幅绝均差的最大值仅5.61 cm,迟角绝均差的最大值仅9.13°。这些模式的调和常数给出潮波传播特征差别不大。基于这些模式提供的调和常数,分别建立了北印度洋4、8和16分潮潮汐预报模型,将预报结果与中国海事服务网提供的沿岸24个站潮汐表资料进行对比。各模式的8分潮(M_2、S_2、N_2、K_2、K_1、O_1、P_1、Q_1)潮汐预报模型均优于4分潮(M_2、S_2、K_1、O_1)潮汐预报模型,NAO.99b模式可以提供16分潮(M_2、S_2、N_2、K_2、K_1、O_1、P_1、Q_1、MU_2、NU_2、T_2、L_2、2N_2、J_1、M1、OO_1)潮汐预报模型,但是对预报结果改善不明显;在各模式中,GOT00.2模式的8分潮潮汐预报模型对北印度洋沿岸的预报效果最好,平均绝均差为14.97 cm。     

Tidal current fluctuations in the Soya Current

Long-term current measurements were carried out near the Soya Strait in the Okhotsk Sea during a period from February 1980 to September 1982. The data were divided into five segments, each being 150 days long, and the tidal ellipse parameters of major axis, minor axis, orientation, and phase for the four major constituents (M2, S2, K1 and O1 tides) were calculated at each segment. The major axis of the mean tidal ellipse averaged over five segments was 29.9 cm sec^–1 for O1 tide, 28.3 cm sec^–1 for K1 tide, 10.4 cm sec^–1 for M2 tide, and 3.7 cm sec^–1 for S2 tide. The phase and orientation of the tidal ellipse were much stable. But, the root mean square deviations of the major axis reached 20% of the mean values for all four constituents. The tidal currents estimated from the sea level records at Wakkanai and Esashi along the Hokkaido coast in the Okhotsk Sea show that their amplitudes and phases are in good agreement with the observed ones for all four constituents.     

基于FVCOM的厦门湾及其周边海域三维潮流数值模拟

基于非结构三角形网格、干-湿判别技术和有限体积法的FVCOM（finite volume coastal oceanmodel）海洋数值模式,建立了厦门湾及其周边海域高分辨率（30 m）的三维潮汐、潮流数值模型.模拟结果同该海域2个验潮站和4个连续海流站的观测资料符合良好,较好地反映了厦门湾及其周边海域潮汐、潮流运动的变化状况和分布特征,并给出了M2、S2、K1、O1共4个主要分潮的同潮图、表层潮流椭圆、最大可能潮流流速及表、底层潮余流分布.厦门湾及其周边海域属正规半日潮类型,4个分潮的最大潮汐振幅分别为200、65、36、29 cm,厦门湾内外迟角差分别为20°、25°、18°、10°;镇海角至围头角连线东南侧湾口区为逆时针旋转的驻波,西北侧湾内为前进潮波.湾内潮流属正规半日潮流,湾口区潮流以逆时针方向的旋转流运动为主,湾内各水道为往复潮流,椭圆长轴与水道走向一致,4个分潮流表层最大流速分别为201、51、34、25 cm/s;九龙江口区3条港道内的流速以南港最大;表层余流大于底层余流,二者水平分布形态基本一致,都为北进南出.     

渤海湾水位、流场的一个数值预报模式

本文利用普林斯顿海洋模式(POM),建立了渤、黄海风增水预报模式和渤海湾水位海流预报模式;利用正交曲线网格提高重点区域的分辨率;采用美国海洋大气局(NOAA)全球预报风场和气压场作为模式表面强迫场,将计算域边界上的天文潮预报值与风增水模式预报的余水位相叠加构建模式的边界条件,在正压条件下,模拟了渤海湾2002年的水位流场过程。结果表明,模式能够较好地再现计算域内天文潮和综合水位的预报,域内10个潮位站模式与实测分析的m1和M2分潮的振幅与迟角差均不超过5.1cm和6.3°,15个潮流站模式与实测分析的m1和M2分潮流的振幅与迟角差均不超过7.5cm/s和15.8°,模式预报的水位值与塘沽站实测值非常接近,预报精度较单纯的天文潮预报有明显提高。     

风暴潮预报模式在渤海海域中的应用研究

首先运用COHEREN S三维多功能大陆架水动力学模型,对渤海四个主要分潮进行了调和分析,建立了渤海天文潮预报模型。结合9712号热带风暴的路径记录资料,运用风场、气压场估算公式,估算了风暴所经地带风场、气压场的变化过程。以此为输入条件,运用COHEREN S模型模拟了9712号热带风暴作用下渤海的增水过程,还模拟讨论了底摩阻系数对水位的影响。     

台湾海峡及其邻近海域潮汐数值计算

建立二维潮波模式,模拟了台湾海峡及其邻近海域（18～30°N,110～130°E）八个主要分潮（M2、S2、K1、O1、P1、Q1、K2、N2）,并利用中国大陆及环台湾岛20多个潮位站的实洲资料进行验证,计算结果与实测值吻合良好。此外,给出了八个主要分潮的同潮图,并逐个讨论了潮汐特征。结果艟示：（1）台湾海峡中的潮波运动是北部蜕化了的旋转潮波系统和南部的前进潮波系统共同作用的结果。（2）半日分潮南、北两支潮波在台湾海峡中部汇合,而今日分潮则在台湾海峡南部海域汇合后继续朝西南方向传播。（3）半日分潮振幅最高值发生在福建省湄洲湾－兴化湾一带,全日分湖最高值则出现在雷州半岛以东一带近岸海域。（4）N2、K2和O1、P1、Q1分湖的振幅、迟角分布分别同M2与K1分潮的整体分布趋势相似。     

北部湾三维潮流场数值研究(英文)

本文对北部湾三维潮波进行了数值研究,并将所得结果与现有观测资料作了比较。通过比较发现两者相当一致。文中还绘制了M_2和K_1分潮的同潮图,并对其作了讨论。可以断定,北部湾的日潮优势是由日潮共振所引起的。M_2和K_1分潮流的水平分布表明,强流区位于它们各自的蜕化无潮点附近。还讨论了分潮流的铅直结构,可以看出这种结构十分接近实际分潮流场的特征。     

基于FVCOM的廉州湾及周边海域三维潮汐潮流数值模拟

基于采用不规则三角网格和有限体积方法的FVCOM模式,建立廉州湾三维潮流数值模型来重现廉州湾及周边海域的潮位和潮流变化状况。根据模拟结果计算得到了较以往更为精细的廉州湾及周边海域K₁、O₁和M₂分潮的同潮图,并计算了由此3个分潮引起的潮汐不对称的变化情况。K₁和O₁分潮在廉州湾外主要以驻波的形式存在,进入廉州湾后转化为前进波;M₂分潮在廉州湾外主要以前进波的形式存在,进入廉州湾后前进波特征更为明显。K₁和O₁分潮流在廉州湾外以旋转流为主,在廉州湾内以往复流为主;M₂分潮流在整个研究海域以往复流为主。由潮余流场的分布特点可以看出自南向北由外海进入廉州湾的潮余流,在冠头岭处分为两支,一支逆时针转向西,另一支被冠头岭阻挡在其南侧形成顺时针封闭环流。在廉州湾内部同时存在2个环流系统,湾顶的气旋式环流和口门处的反气旋式环流。     

Surface tidal currents in the open sea area to the east of the Zhoushan Islands measured with high frequency surface wave radar

Based on the quasi-harmonic analysis of 11 d vector ocean currents obtained from two high frequency sur- face wave radars located at Zhujiajian Island and Shengshan Island, the spatial distribution characteristics of surface tidal currents in the open sea area to the east of the Zhoushan Islands of Zhejiang Province, China are studied. The following conclusions are drawn from the analysis： the tidal current pattern in the open sea area to the east of Zhoushan Islands is primarily regular semidiurnal, which is significantly affected by the shallow water constituents. The directions of the major axes of tidal current ellipses of M2 lie approx- imately in the NW-SE direction. With the increasing of distance away from the coast, the directions of the tidal current ellipses gradually shift toward the E-W direction. The tidal currents are mainly reversing cur- rents. The spatial distribution of probable maximum current velocities decreases gradually from northeast to southwest which is basically in accordance with the spatial distribution of measured maximum current velocities. The residual currents near the coast are larger than those far away from the coast. The directions of the residual currents are basically north by east, and the angle to the due north increases gradually with the increasing distance away from the coast. The topography shows a certain impact on the spatial distri- bution of shallow water constituents, the rotation of tidal currents, the probable maximum currents and the residual currents.     

不考虑局地引潮势的南海正压潮能通量与潮能耗散

利用ECOM模式模拟南海正压M2、S2、K1、O1分潮, 对南海潮能通量及潮能耗散进行研究.结果显示, M2、S2、K1和O1分潮分别有38.93、5.77、29.73和28.97GW的能通量经吕宋海峡传入南海, 并有2.42、0.36、8.67和7.86GW的能通量由南海经卡里马塔海峡传入爪哇海.由东海及吕宋海峡西北部传入台湾海峡的M2分潮能通量为25.28GW.半日潮进入北部湾和泰国湾的能通量较少(6.52GW), 全日潮则较大(24.74GW).通过民都洛和巴拉巴克海峡断面, 全日潮由南海向苏禄海共输送12.28GW的能通量, 而半日潮则由苏禄海向南海输送1.92GW的能通量.由模式输出结果估计得到的南海各局部海域的底摩擦耗散与净潮能通量存在差异, 为使二者平衡, 可对南海不同海域的底摩擦系数进行调整.依净潮能通量与底摩擦耗散平衡关系计算得到台湾海峡、北部湾、泰国湾及南海深水海域的底摩擦系数分别为0.0023、0.0024、0.0023和0.0021.     

渤、黄、东海8个主要分潮的数值模拟研究

应用MIKE数值模拟软件,采用无结构三角形网格,建立一套计算区域包括整个渤海、黄海、东海以及东海大陆架和琉球群岛的高分辨率数值模型,考虑了实际水深和岸线,外海开边界采用西北太平洋大模型结果的潮位提供,模拟了东中国海潮波的波动过程,对潮波垂直运动过程进行调和分析,得到了渤海、黄海、东海的M2,S2,K1,O1以及N2,K2,P1,Q1八个主要分潮的传播和分布特征。利用中国沿海14个潮位站的调和常数对模型结果进行了验证,验证结果显示模型较为准确可靠。研究结果表明:4个主要半日潮(全日潮)在渤、黄、东海的传播情形基本相似,即潮波在渤海、黄海、东海沿岸的传播性质上类似沿岸开尔文波的传播形态,并且成功再现了计算海域的4个半日分潮无潮点和2个全日分潮无潮点。全日潮振幅各无潮点附近振幅最小,而海湾的波腹区振幅最大,东海潮差呈现近岸方向振幅大、离岸方向振幅小,浙闽沿海振幅也较大,黄海振幅相对较小,渤海振幅在辽东湾和渤海湾顶最大,两个无潮点周边振幅较小。     

Modelling tidal currents on the coast of Portugal

A tide circulation model of the Atlantic coast of the Iberian Peninsula has been constructed. This regional numerical model covers the whole continental shelf. The finite element computational grid is made of some 16,300 triangular elements with sizes ranging between 13 km (on the offshore boundary) and 1 km (near the coast), with local refinements on the continental shelf and in the area of Figueira da Foz. This site was selected as experimental site for the study of waves and currents in the frame of the MAST2/WAVEMOD research project.Boundary conditions along the three oceanic limits of this widely open domain are obtained from the North-Atlantic component of a World Ocean tidal numerical model known as FES94 [Le Provost, C., Genco, M.L., Lyard, F., Vincent, P. and Canceil, P. (1994) Spectroscopy of the world ocean tides from a finite element hydrodynamic model, Journal of Geophysical Research, Topex-Poseidon Special Issue]. A new radiation-like boundary condition has been introduced in the modelling system used (Telemac-2D), which solves the Shallow Water Equations (SWE), in order to interface the two models and to allow for the tidal wave to leave the northern limit without reflection. Model calibration has been performed on the dominating M₂ constituent. The introduction of the astral static potential generating the tide in the SWE improved this regional model.A long duration run (1 month) has been performed, the model being forced by the eight major tide constituents. Harmonic analysis of results has been performed on 17 tide constituents, due to non-linear interactions of constituents on the continental shelf. Comparisons with the FES94 model on one hand, and with a set of coastal tide gauges on the other hand, are good. A database of tidal harmonics is now available for forecasting sea levels and currents in this area.This work has shown that diurnal shelf trapped tide waves exist in some places along the Portuguese continental shelf, which induce diurnal dominant tidal currents in these places (North of Figueira da Foz).