Ocean tides around New Zealand

The distribution of amplitude and phase for eight ocean tidal constituents (M₂, S₂, N₂, K₂, K₁, O₁, P₁, Q₁) is presented as tidal maps for the New Zealand area. The distribution was calculated using a barotropic tidal model driven by TOPEX/ Poseidon data on the outer ocean boundaries. The maps exhibit the known features of the tides in this area such as a complete rotation of the semi‐diurnal tides around New Zealand and the reduced spring‐neap variations on the east coast. They also point out several new features for which there are few or no observations, such as diurnal trapped waves and shelf waves. A comparison of the model results with observations shows that sea level errors are within 0.1 m in amplitude and 10° in phase for the largest constituents at all locations, including sites where the data are of low quality and where the geometry is not adequately resolved. For locations where the geometry is adequately represented and the observations are of high quality, sea level errors are within 0.02 m in amplitude and 7° in phase. These results represent the most accurate and highest resolution calculations of tides and currents yet attained for this area.     

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.     

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

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

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.     

基于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。此外,插值所得四个主要分潮的无潮点位置与前人研究结果相近。上述实验结果表明:利用双调和样条插值方法对卫星高度计所得调和常数进行插值能够获取较为准确的同潮图。     

渤海潮波系统的数值模拟

利用二维非线性潮波方程组,讨论了渤黄海主要分潮(全日潮、半日潮及浅水分潮) 数值模拟中的有关问题。数值模拟中同时考虑了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%左右。根据计算结果分析了印尼海域的潮汐特征及潮能传播规律,结果显示:爪哇海以外的印尼海域主要为不规则半日潮区;全日潮潮能主要由太平洋传入印尼海域,而半日潮潮能则是从印度洋传入印尼海域。     

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.     

Three-dimensional numerical model study of the residual current in the South China Sea

A three-dimensional baroclinic shelf sea model is employed to simulate the tidal and non-tidal residual current in the South China Sea. The four most significant constituents, M₂, S₂, K₁ and O₁, are included in the experiments with tidal effect. At most stations, the computed harmonic constants agree well with the observed ones. The circulations of the South China Sea in summer (August) and winter (December) are mainly discussed. It is shown that the barotropic tidal residual current is too weak to affect the South China Sea circulation, whilst the contribution of the baroclinic tidal residual current to the South China Sea circulation would be important in the continental shelf sea areas, especially in the Gulf of Thailand and Gulf of Tonkin. In the deep-sea areas, the upper barotropic or baroclinic tidal residual current is relatively very weak, however, the speed order of the deep baroclinic tidal residual current can be the same as that of the mean current without tidal effect. Moreover, the baroclinic tidal residual current seems to be related to the different seasonal stratification of ocean.     

Tidal Currents in the Tsushima Straits Estimated from ADCP Data by Ferryboat

Tidal currents in the Tsushima Straits have been analyzed using measurements obtained since February 1997 by an acoustic Doppler current profiler (ADCP) mounted on the ferryboat Camellia. Tidal current constituents (M ₂, S ₂, K ₁, O ₁) are dominant among the ten tidal current constituents (Q ₁, O ₁, P ₁, K ₁, N ₂, M ₂, S ₂, K ₂, MSf, Mf), and generally 1.4–2.1 times stronger at the western channel of the straits than those at the eastern channel. The ratio between amplitude of M ₂, S ₂, K ₁ and O ₁ averaged along the ferryboat track is 1:0.45:0.59:0.51. The major axis directions of tidal current ellipses are generally SW to NE, exceptionally in the vicinity of the Tsushima Islands. Approaching the Tsushima Islands from the Korean Peninsula side, the major axis gradually rotates clockwise. At the western channel, the M ₂ and K ₁ constituents change the rotation direction of current vectors from clockwise to counterclockwise at about 90–130 m depth. The contributions of the tidal currents to the mean kinetic energy and the mean eddy kinetic energy along the ferryboat track are, on average, 0.56 and 0.71, respectively. This suggests that tidal current activities are generally more dominant than the mean current activities and much more dominant than eddy activities. The only region where the eddy activities are comparable to the tidal current activities is located on the east side of the Tsushima Islands. This revised version was published online in July 2006 with corrections to the Cover Date.     

内潮耗散与自吸-负荷潮对南海潮波影响的数值研究

利用非结构三角形网格的FVCOM海洋数值模式,在其传统二维潮波方程中加入参数化的内潮耗散项和自吸-负荷潮项,计算了南海及其周边海域的M_2、S_2、K_1和O_1分潮的分布。与实测值的比较表明,引入这两项对模拟准确度的提高有明显效果。根据模式结果本文计算分析了研究海域的潮能输入和耗散。能量输入计算表明,能通量是潮能输入的最主要构成部分,通过吕宋海峡断面进入南海的M_2和K_1分潮能通量分别为38和29GW;半日周期的自吸-负荷潮能量输入以负值居多,而全日周期的自吸-负荷潮能量输入以正值居多,因而自吸-负荷潮减弱了南海的半日潮,并加强了南海的全日潮。引潮力的作用也减弱了半日潮而加强了全日潮,但其作用要小于自吸-负荷潮。潮能耗散的分析显示底摩擦耗散在沿岸浅水区域起主导作用,内潮耗散则主要发生在深水区域。内潮耗散的最大值出现在吕宋海峡,且位于南海之外的海峡东部的耗散量大于位于南海之内的海峡西部的耗散量。对M_2和K_1分潮吕宋海峡的内潮耗散总值分别达到16和23GW。     

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.     

基于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°,结果显示各分潮模拟结果合理,较真实地反映了渤海海域四个分潮传播情况。     

粤西沿岸海域潮汐特征

对粤西海域水东、湛江、硇州岛、南渡和海安验潮站各1年水位资料进行了调和分析及统计。结果表明:粤西海域北部湛江、水东、硇洲岛、南渡4站主要分潮为M_2、K_1、S_2和O_1,是不规则半日潮;而南部海安站主要分潮为O_1和K_1分潮,是不规则日潮。粤西沿岸各站浅水分潮和平均水位从北到南有减小趋势。在日潮不等方面,粤西沿岸验潮站都存在明显的日潮高、日潮时不等现象。日潮时不等总体上从北到南有降低的趋势,北部湛江、水东和硇州岛涨潮时间比落潮长1~1.5 h,南部南渡和海安两站涨、落潮时相差不大。另外,调和分析和经验公式不再适用于南渡站,南渡河入海口处堤闸工程是主要原因。这为粤西海域环境资源开发、航运活动等提供环境支持和保障。     

胶州湾潮汐潮流动边界数值模拟

基于普林斯顿海洋模式,通过干湿网格判别法引入潮汐潮流的漫滩过程,考虑M₂,S₂,K₁,O₁,M₄和MS₄六个主要分潮,建立了胶州湾潮汐潮流数值模拟和预报模型,研究了该海域潮汐潮流特征,并讨论了漫滩对潮流模拟的影响。与实测资料的对比验证表明,该模式能够对胶州湾的潮汐和潮流做出较为合理的预测。给出了胶州湾潮汐、潮流、余流等分布特征,模拟的潮流场以及余流场涡旋等现象与观测符合良好;计算了潮波能通量,从能量角度探讨了潮波的传播特性;对潮位与潮流场演变规律,以及潮能通量的分析表明,胶州湾内的潮波以驻波为主。通过数值试验发现,漫滩过程的引入对胶州湾潮流速度的模拟至关重要,不考虑漫滩过程的模式会夸大或者低估潮流流速。对于滩涂面积广阔的海域来说,潮流数值模式中考虑漫滩的影响是必要的。     

Simulation of tides,residual flow and energy budget in the Gulf of California

With the application of a two-dimensional nonlinear hydrodynamical-numerical semi-implicit model, the principal tides M₂, S₂, K₂, N₂, K₁, P₁ and O₁ were studied. Energy budgets of the semi-diurnal M₂ and S₂ were calculated separately. The linear sum of these budgets was compared with the tidal energy budget obtained when these two tidal constituents interact. Since a quadratic form for the bottom friction was used, remarkable differences were found. The results show that in the area of the Colorado River delta, the dissipation of tidal energy is very strong. Intense tidal currents were observed in the same region and over the Salsipuedes Sill. Energy budgets calculated for forcing waves of different periods, but of the same amplitude, were used to estimate the principal periods of resonance. Although the topography of the Gulf is very complex, the model reproduced observed sea-surface elevation and current patterns. To study spring tide conditions, the above seven tidal constituents were simulated. Estimates of residual currents reveal the presence of several intense cyclonic and anticyclonic gyres. Over the Salsipuedes Sill, residual currents of the M₂ tide reach values of more than 15 cm s⁻¹. Horizontal distributions of dissipation rates of tidal energy and of kinetic energy were also obtained.     

潮汐非调和常数的计算方法——Ⅰ.原理

潮汐非调和常数是一组用来描述港口潮汐特征的数据,应当从长期观测的高、低潮时间及高度中统计得出。然而在只有短期记录的情况下,由于潮汐存在着长周期的变化,直接统计得出的结果会出现相当大的误差。因此,在只有短期验潮的港口,人们常常利用潮汐的调和常数间接计算非调和常数;而调和常数则可利用短期记录求得,不会带来严重的误差。各国使用的潮汐非调和常数计算方法不尽相同。在我国,郑交振同志编著的《实用潮汐学》一书中全面地介绍了     

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.     

Applications of the Mild-Slope Equation to Tidal Computations in the Taiwan Strait

Employing harmonic analysis of tidal data in the Taiwan Strait, the cross-strait tidal characteristics are completely illustrated. Based on the two dimensional mild-slope equation which can be reduced to the shallow-water wave equation, a finite element model (Tsay et al., 1989) is applied to investigate the characteristics of tides in the Taiwan Strait. The co-range and equi-phase charts of major tidal constituents, such as M₂, S₂, N₂, and K₁, are reproduced. Anomalous amplification of semidiurnal tides in the Taiwan Strait is verified. With rotation effects neglected and by applying a non-reflective condition on the open boundaries, the numerical results of phase-lag and co-range distributions show very good agreement with observed data for semidiurnal tides in the Taiwan Strait. Due to crude representation of the topography at two ends along the China coast, computed tidal distributions deviate from the observations. However, both computed amplitudes and phase-lags compare very well with observed data along the central half of the China coast.     

基于SCHISM模式的全球潮波模拟

基于非结构网格半隐式跨尺度海洋模式(semi-implicit cross-scale hydroscience integrated system model,SCHISM),作者采用非结构三角网格,对全球大洋潮波进行数值模拟.通过调和分析,将196个潮位站的实测数据与模拟结果进行比较验证,两者符合良好,M2、K1分...