Data Assimilation Modeling of the Barotropic Tides in the Korea/Tsushima Strait

During 1999–2000, 13 bottom mounted acoustic Doppler current profilers (ADCPs) and 12 wave/tide gauges were deployed along two lines across the Korea/Tsushima Strait, providing long-term measurements of currents and bottom pressure. Tidally analyzed velocity and pressure data from the moorings are used in conjunction with other moored ADCPs, coastal tide gauge measurements, and altimeter measurements in a linear barotropic data assimilation model. The model fits the vertically averaged data to the linear shallow water equations in a least-squares sense by only adjusting the incoming gravity waves along the boundaries. Model predictions are made for the O₁, P₁, K₁, μ₂, N₂, M₂, S₂, and K₂ tides. An extensive analysis of the accuracy of the M₂ surface-height predictions suggests that for broad regions near the mooring lines and in the Jeju Strait the amplitude prediction errors are less than 0.5 cm. Elsewhere, the analysis suggests that errors range from 1 to 4 cm with the exception of small regions where the tides are not well determined by the dataset. The errors in the model predictions are primarily caused by bias error in the model’s physics, numerics, and/or parameterization as opposed to random errors in the observational data. In the model predictions, the highest ranges in sea level height occur for tidal constituents M₂, S₂, K₁, O₁, and N₂, with the highest magnitudes of tidal velocities occurring for M₂, K₁, S₂, and O₁. The tides exhibit a complex structure in which diurnal constituents have higher currents relative to their sea level height ranges than semi-diurnal constituents.     

Observed and modelled tidal currents in the New Zealand region

Tidal currents derived from current meter measurements are compared with the output from a barotropic tidal model of the New Zealand region. For the semi‐diurnal constituents there was very good agreement for the M₂ tide and good agreement for the S₂ tide. For the diurnal constituents (K_l, O_l) it was found that as the amplitude of the constituents decreased so did both the model/observation agreement and the accuracy of the observed tidal ellipse parameters. Consequently it was not possible to decide whether differences arose through shortcomings in the model or in the data. However, the overall performance of the model as a prognostic tool for ocean tidal current simulation appears to be good.     

Numerical study of tidal dynamics in the South China Sea with adjoint method

We adopt a parameterized internal tide dissipation term to the two-dimensional (2-D) shallow water equations, and develop the corresponding adjoint model to investigate tidal dynamics in the South China Sea (SCS). The harmonic constants derived from 63 tidal gauge stations and 24 TOPEX/Poseidon (T/P) satellite altimeter crossover points are assimilated into the adjoint model to minimize the deviations of the simulated results and observations by optimizing the bottom friction coefficient and the internal tide dissipation coefficient. Tidal constituents M₂, S₂, K₁ and O₁ are simulated simultaneously. The numerical results (assimilating only tidal gauge data) agree well with T/P data showing that the model results are reliable. The co-tidal charts of M₂, S₂, K₁ and O₁ are obtained, which reflect the characteristics of tides in the SCS. The tidal energy flux is analyzed based on numerical results. The strongest tidal energy flux appears in the Luzon Strait (LS) for both semi-diurnal and diurnal tidal constituents. The analysis of tidal energy dissipation indicates that the bottom friction dissipation occurs mainly in shallow water area, meanwhile the internal tide dissipation is mainly concentrated in the LS and the deep basin of the SCS. The tidal energetics in the LS is examined showing that the tidal energy input closely balances the tidal energy dissipation.     

泰国湾及邻近海域潮汐潮流的数值模拟

本文基于FVCOM(Finite-Volume Coastal Ocean Model)模式,模拟了泰国湾及其周边海域K1、O1、M2和S2四个主要分潮。采用47个验潮站实测调和常数与模拟结果进行比较,所得4个分潮的均方差分别为4.06cm、3.76cm、8.22cm和4.71cm,符合良好。根据计算结果分析了泰国湾及其周边海域的潮汐、潮流的分布特征和潮波的传播特征。数值试验表明,现有的数字水深资料(ETOPO1,ETOPO5,DBDB-V)的准确度不足以合理地模拟泰国湾潮波。     

Numerical modelling of the sea level under the actions of the tide and strong wind in the Bohai Sea

In this paper, the numerical modelling of the tidal level and current in the Bohai Sea was carried out with ADI method, by taking the sum of four main tidal components M2,S2K2,O1 as the open boundary condition. The calculated values were consistent with the predicted ones (the observed values in the case of calm) in the Tidal Table. On the basis of the modelling of the tide, the sea level and current fields under the effects of strong wind were simulated. The calculated results were also quite satisfactory.     

海南岛莺歌海近岸的潮汐不对称与潮致余流研究

潮汐不对称与潮致余流在河口海岸区的物质输运中扮演着重要角色。已往的研究表明,在驻波占主导的河口海湾中,涨落潮的历时不对称与流速不对称有较为良好的对应关系。而潮致余流主要由地形与潮波的非线性作用所致。本研究以海南岛莺歌海附近为代表,结合实际观测与数值模型,研究复杂地形的开阔近岸区的潮汐不对称与潮致余流。结果表明,在莺歌海近岸区,涨落潮历时不对称皆表现为涨潮历时短于落潮历时,而流速不对称则出现复杂的空间变化。对流速不对称的机制分解表明,研究区的流速不对称主要由K₁、O₁与M₂的相互作用,以及潮余流与各潮汐分潮的相互作用所控制。其中前者产生涨潮流速大于落潮流速的涨潮优势,而后者则与余流的方向相对应,出现多个涨潮优势与落潮优势的区域。总体而言,研究区的流速不对称由余流与各潮汐分潮的相互作用所决定。这表明,采用涨落潮历时的不对称来确定潮汐不对称的方法在开阔近海区可能并不适用。对潮致余流的研究表明,研究区的欧拉余流远大于斯托克斯余流。欧拉余流表现为多个顺时针与逆时针的涡流。涡流分布与地形具有较好的对应关系,潮流沙脊区多发育顺时针涡流,而深槽区则以发育逆时针涡流为主。摩擦力在涡流的发育中起着重要作用。     

基于T/P和Jason-1/2卫星高度计海面高度数据获取南海同潮图的研究

随着卫星高度计资料的不断丰富,通过对卫星高度计所得潮汐调和常数进行插值或拟合得到潮汐同潮图成为可能。本文拟对T/P(TOPEX/POSEIDON)、Jason-1和Jason-2卫星高度计数据进行分析,得到南海区域星下观测点处四个主要分潮(M2、S2、K1和O1分潮)的调和常数,进而利用双调和样条插值方法对其进行插值,获取南海同潮图。首先,以1992～2016年T/P和Jason卫星高度计所得海面高度数据为基础,利用调和分析方法计算了南海星下观测点处M2、S2、K1和O1四个主要分潮的调和常数,并与40个验潮站数据进行了对比,最大矢量均差为4.99cm,说明分析所得调和常数与利用验潮站资料提取的调和常数的误差较小。进而采用双调和样条插值方法对星下点调和常数进行插值,得到了南海四个主要分潮的同潮图,所得结果与全球潮汐模型TPXO7.2模式结果的矢量均差分别为4.69、2.46、3.13和2.42 cm,与141个验潮站处观测结果的矢量均差分别为22.59、10.26、10.24和8.51 cm。此外,插值所得四个主要分潮的无潮点位置与前人研究结果相近。上述实验结果表明:利用双调和样条插值方法对卫星高度计所得调和常数进行插值能够获取较为准确的同潮图。     

Three-dimensional structure of tidal current in the East China Sea and the Yellow Sea

A three-dimensional tidal current model is developed and applied to the East China Sea (ECS), the Yellow Sea and the Bohai Sea. The model well reproduces the major four tides, namely M₂, S₂, K₁ and O₁ tides, and their currents. The horizontal distributions of the major four tidal currents are the same as those calculated by the horizontal two-dimensional models. With its high resolutions in the horizontal (12.5 km) and the vertical (20 layers), the model is used to investigate the vertical distributions of tidal current. Four vertical eddy viscosity models are used in the numerical experiments. As the tidal current becomes strong, its vertical shear becomes large and its vertical profile becomes sensitive to the vertical eddy viscosity. As a conclusion, the HU (a) model (Davieset al., 1997), which relates the vertical eddy viscosity to the water depth and depth mean velocity, gives the closest results to the observed data. The reproduction of the amphidromic point of M₂ tide in Liaodong Bay is discussed and it is concluded that it depends on the bottom friction stress. The model reproduces a unique vertical profile of tidal current in the Yellow Sea, which is also found in the observed data. The reason for the reproduction of such a unique profile in the model is investigated.     

Numerical computation of tidal currents in the eastern and southern parts of the East China Sea

Numerical simulation on diurnal and semidiurnal tidal currents is made for the eastern and southern parts of the East China Sea bounded by the Tsushima, islands of Kyushu and Ryukyu in the east and extending to the central Taiwan Strait in the south. The numerical computation shows that the rapid variation in the seabed topography near the Okinawa Trough may cause computational instability and this instability can be removed by setting a suitable lateral eddy viscosity. The computed results are in good agreement with the observations and indicate that the tidal currents in the greater part of the area studied rotate clockwise. An M2 current-amphidromic point exists in the southern part of the Korea Strait and another one appears in the west of Okinawa with less certainty. K1 and M2 cotidal charts based on the observed and computed tidal currents in the sea and the observed tides on the coast are also given.     

Phase distribution of tidal constituents around New Zealand

The phase distributions of the M₂, S₂, K₁, and O₁ tidal constituents around New Zealand are plotted from existing harmonic analyses of tidal heights. Both semidiurnal constituents exhibit a complete 360° range of phase around New Zealand, with complex areas of rapid phase change through or near the strait separating the two main islands. The K₁ amphidrome and that for O₁, which previously were thought to be centred on New Zealand, are shown to be located cast of New Zealand. The distributions plotted highlight areas where tidal observations are lacking.     

渤海潮波系统的数值模拟

利用二维非线性潮波方程组,讨论了渤黄海主要分潮(全日潮、半日潮及浅水分潮) 数值模拟中的有关问题。数值模拟中同时考虑了4个主要分潮(M2,S2,K1,O1)和两个浅水分潮(M4,MS4)。分析表明,在渤黄海潮波系统数值模拟中,稳定后选取14 d的数值模拟结果进行调和分析能够取得最佳(最合理)的调和分析结果。计算出调和常数的模拟值与实测值之差的绝对平均值:M2分潮的振幅差为4cm,迟角差为3.3°,S2分潮的振幅差为2cm,迟角差为4.2°,K1 分潮的振幅差为1cm,迟角差为3.7°,O1分潮的振幅差为2 cm,迟角差为5.5°。实验结果较好地体现了渤黄海潮波系统的特征。     

印度尼西亚海域潮波的数值研究

基于ROMS模式构建了模拟区域为(15.52°S-7.13°N,110.39°~134.15°E)水平分辨率为2′的潮波数值模式,分别模拟了印尼海域M2、S2、K1、O1四个主要分潮。模拟结果与29个卫星高度计交叠点上的调和常数进行比较,符合较好。M2分潮的振幅均方根差为3.4cm,迟角均方根差为5.9°;S2分潮的振幅均方根差为1.7cm,迟角均方根差为6.3°;K1分潮振幅均方根差为1.1cm,迟角均方根差为5.8°;O1分潮振幅均方根差为1.2cm,迟角均方根差为4.4°。M2、S2、K1、O1分潮向量均方根差分别为3.8cm、2.4cm、1.9cm和1.3cm,模拟结果的相对偏差在10%左右。根据计算结果分析了印尼海域的潮汐特征及潮能传播规律,结果显示:爪哇海以外的印尼海域主要为不规则半日潮区;全日潮潮能主要由太平洋传入印尼海域,而半日潮潮能则是从印度洋传入印尼海域。     

Cotidal charts and tidal power input atlases of the global ocean from TOPEX/Poseidon and JASON-1 altimetry

The global distributions of eight principal tidal constituents, M₂ , S₂ , K₁ , O₁ , N₂ , K₂ , P₁ , and Q₁ , are derived using TOPEX/Poseidon and JASON-1（T/P-J） satellite altimeter data for 16 a. The intercomparison of the derived harmonics at 7000 subsatellite track crossover points shows that the root mean square （RMS） values of the tidal height differences of the above eight constituents range from 1.19 cm to 2.67 cm, with an average of about 2 cm. The RMS values of the tidal height differences between T/P-J solutions and the harmonics from ground measurements at 152 tidal gauge stations for the above constituents range from 0.34 cm to 1.08 cm, and the relative deviations range from 0.031 to 0.211. The root sum square of the RMS differences of these eight constituents is 2.12 cm, showing the improvement of the present model over the existing global ocean tidal models. Based on the obtained tidal model the global ocean tidal energetics is studied and the global distribution of the tidal power input density by tide-generating force of each constituent is calculated, showing that the power input source regions of semidiurnal tides are mainly concentrated in the tropical belt between 30 S and 30 N, while the power input source regions of diurnal tides are mainly concentrated off the tropic oceans. The global energy dissipation rates of the M₂ , S₂ , K₁ , O₁ , N₂ , P₁ , K₂ and Q₁ tides are 2.424, 0.401, 0.334, 0.160, 0.113, 0.035, 0.030 and 0.006 TW, respectively. The total global tidal dissipation rate of these eight constituents amounts to 3.5 TW.     

基于FVCOM 的渤海潮波数值模拟

基于有限体积法海洋数值模型(FVCOM),对渤海当前水深岸线状况下的潮汐潮流进行了数值计算。模式采用不规则三角形网格,较好地提高了黄河口处网格分辨率,模拟了渤海海域K1,O1,M2和S2四个主要分潮。利用渤海沿岸19个验潮站的资料对模拟结果进行了验证,K1分潮振幅绝均差2.39 cm,迟角绝均差4.36°,O1分潮振幅绝均差1.40 cm,迟角绝均差4.29°,M2分潮振幅绝均差为3.55 cm,迟角绝均差为5.69°,S2分潮振幅绝均差1.72 cm,迟角绝均差8.86°,结果显示各分潮模拟结果合理,较真实地反映了渤海海域四个分潮传播情况。     

烟墩角海域二测站周日水温及潮流变化特征

为加深对山东半岛东部养殖区的海流、海温变化特征,及其与该地区生态特征分布之间关系的认识,应用短期资料的潮流准调和分析方法,计算了烟墩角海域两个测站的O1,K1,M2,S2,M4,MS46个主要分潮的北、东分量潮流调和常数,并给出了各测站在各层的潮流椭圆要素。结果表明:该海区水温呈现层化结构,在1天内呈现两个波动,与N-S向潮流变化关系密切。海区除湾内表层外属日潮性质外,潮流属于正规半日潮流性质,半日分潮流的北分量大于东分量,呈现典型往复流特征,最大流流向多为南北向。分析结果也反映出表层、中层的浅水分潮在观测海流中所占的份额高于底层的浅水分潮。     

南海北部陆架陆坡区海流观测研究

针对2006-2009年期间,南海北部陆架陆坡区3个站ADCP海流连续观测资料,采用功率谱分析、潮流调和分析方法,重点分析了陆架陆坡区100 m,200 m和1 200 m水深海域海流的垂向结构,探讨了环流的季节变化和空间分布特征,特别讨论了南海暖流和北陆坡流的时空变化特征。结果表明,陆架陆坡区潮流类型属于不规则日潮,深水站点中层表现为正规全日潮类型,垂向为"三层结构",甚至更加复杂。O₁,K₁,M₂,S₂等分潮总体上为顺时针旋转,在深水站点,基本表现为西北－东南走向的往复流形态。从能量角度看,表层和底层海流中,潮流所占份额较大,分别占30%~40%和40%~50%,中层较小,约为20%。对东沙群岛西南陆架陆坡区环流,观测计算结果证实了西向强流的存在,且垂向结构具有显著的季节变化,在200 m水深处没有明显的南海暖流,只是10~30 m以上层次存在逆风海流。海南岛以东海域连续15个月表层环流的结果表明,冬季明显受到南海暖流的影响,存在东北向的逆风海流,夏秋季的环流表现为西南向,流速较强,夏季也存在逆风情况,造成上述情形的原因可能是该地南海暖流的流轴具有季节性变化——冬季偏南,夏季偏北。     

Tidal transport through the Tsugaru Strait — part I: Characteristics of the major tidal flow and its residual current

The potential role of the tide-induced time-mean flow (the tidal residual current) in determining transport through the Tsugaru Strait (located between the East/Japan Sea and the North Pacific) is investigated using a high-resolution numerical barotropic model. The calculated K₁, O₁, M₂, and S₂ tidal fields agree well with available observational records derived from both tide gauge and current meter measurements in the strait and the adjacent seas. The tidal residual current speed reaches 0.3 ms⁻¹ in two narrow “neck” areas where topographic sills are located. This result suggests that tides should be taken into account in estimating the long-term water mass and nutrient transport through narrow regions between the East/Japan Sea and the North Pacific. An interesting aspect of the tidal residual current field is the prediction of several active eddy zones in which sequences of eddy triplets develop in the vicinity of capes. Our vorticity analysis reveals that the interplay of topographic effects arising from both the headland and the sill around capes plays a critical role in the formation of these triple eddy patterns.     

北部湾潮波数值研究

利用普林斯顿海洋模式(POM08)建立了北部湾及其临近海区潮汐潮流数值模式,模拟了K1,O1,M2和S2这4个主要分潮,分析了模拟的潮汐和潮流分布特征,从潮波能量的角度讨论了琼州海峡对北部湾潮波系统的影响,并给出北部湾潮能的耗散情况。研究表明,北部湾是典型的全日潮海区,K1和O1分潮在南部湾口形成半个旋转潮波系统,无潮点位于越南顺安附近岸边。琼州海峡中的欧拉潮汐余流为西向流,潮余流造成的水通量约为0.034×106m3/s;余流出海峡西口后,先折向北,然后转向南流出湾外。研究海区中两个强潮流区分别位于琼州海峡和海南岛的西侧,同时这也是两个潮能的高耗散区。北部湾的潮能自南部湾口由外海传入,通过西口涌入琼州海峡,到达海峡东口时日潮波的能量已基本耗散殆尽,在海峡内耗散的4个分潮的潮能约为3.33 GW,相当于北部湾潮能耗散量的35%左右。数值试验表明,琼州海峡作为潮能耗散的重要海区,其存在对于北部湾潮波系统的形成具有较大影响。计算了底边界潮能耗散,结果表明在北部湾和琼州海峡,底边界耗散的潮能分别占该海区总耗散的83%和80%。     

Characteristics of Tides in the Bay of Bengal

A vertically integrated 2D numerical model was developed for the simulation of major tidal constituents (M₂, S₂, N₂, K₁ and O₁) in the Bay of Bengal. The bathymetry for the model domain was derived from an improved ETOPO5 dataset prepared in our earlier work. The simulated tidal elevations showed good agreement with the hourly tide gauge observations at Paradip, Visakhapatnam, and Chennai. The amplitudes and phases of M₂, S₂, K₁, and O₁ at the coastal stations, obtained from harmonic analysis of simulated tides, were found to agree well with those obtained from Admiralty Tide Tables with the RMS misfit 9.2, 5.6, 2.9 and 3.1 cm, respectively. In the Bay of Bengal, semi-diurnal tides (M₂, S₂, and N₂) attain highest amplitudes (180, 80, 30 cm, respectively) in the Gulf of Martaban while amplitudes of diurnal tides (K₁, O₁) reaches maximum (20, 12 cm, respectively) in the Malacca Strait. The continental shelf in the head bay and along the southern coast of Myanmar is about 200 km wide and the amplitudes of semi-diurnal tides are doubled in these regions while the diurnal tides amplify only marginally, which is consistent with Clarke and Battisti theory. In the north eastern end of the head bay and the Gulf of Martaban, the geometrical configuration of the coastline, in addition to the wide continental shelf, could contribute to the amplification of both semi-diurnal and diurnal constituents. In the Malacca Strait, the amplitudes of both semi-diurnal and diurnal tides are found to increase gradually from the northern end to the 2.5°N and decreases towards southern boundary. The co-tidal and co-range charts of M₂ and S₂ tidal constituents also show the presence of two degenerate amphidromic points in the head bay. A virtual amphidromic point for M₂ is identified in the Malacca Strait.     

Clockwise phase propagation of semi-diurnal tides in the Gulf of Thailand

The phase of semi-diurnal tides (M₂ and S₂) propagates clockwise in the central part of the Gulf of Thailand, although that of the diurnal tides (K₁, O₁ and P₁) is counterclock-wise. The mechanism of clockwise phase propagation of semi-diurnal tides at the Gulf of Thailand in the northern hemisphere is examined using a simple numerical model. The natural oscillation period of the whole Gulf of Thailand is near the semi-diurnal period and the direction of its phase propagation is clockwise, mainly due to the propagation direction of the large amplitude part of the incoming semi-diurnal tidal wave from the South China Sea. A simplified basin model with bottom slope and Coriolis force well reproduces the co-tidal and co-range charts of M₂ tide in the Gulf of Thailand.