Ocean Tide Models Developed by Assimilating TOPEX/POSEIDON Altimeter Data into Hydrodynamical Model: A Global Model and a Regional Model around Japan

A global ocean tide model (NAO.99b model) representing major 16 constituents with a spatial resolution of 0.5° has been estimated by assimilating about 5 years of TOPEX/POSEIDON altimeter data into barotropic hydrodynamical model. The new solution is characterized by reduced errors in shallow waters compared to the other two models recently developed; CSR4.0 model (improved version of Eanes and Bettadpur, 1994) and GOT99.2b model (Ray, 1999), which are demonstrated in comparison with tide gauge data and collinear residual reduction test. This property mainly benefits from fine-scale along-track tidal analysis of TOPEX/POSEIDON data. A high-resolution (1/12°) regional ocean tide model around Japan (NAO.99Jb model) by assimilating both TOPEX/POSEIDON data and 219 coastal tide gauge data is also developed. A comparison with 80 independent coastal tide gauge data shows the better performance of NAO.99Jb model in the coastal region compared with the other global models. Tidal dissipation around Japan has been investigated for M₂ and K₁ constituents by using NAO.99Jb model. The result suggests that the tidal energy is mainly dissipated by bottom friction in localized area in shallow seas; the M₂ ocean tidal energy is mainly dissipated in the Yellow Sea and the East China Sea at the mean rate of 155 GW, while the K₁ energy is mainly dissipated in the Sea of Okhotsk at the mean rate of 89 GW. TOPEX/POSEIDON data, however, detects broadly distributed surface manifestation of M₂ internal tide, which observationally suggests that the tidal energy is also dissipated by the energy conversion into baroclinic tide.     

Modeling Local Tidal Constituents Using TOPEX/POSEIDON SSH for the Malaysian Marine Region

A new Local Ocean Tide Model, has been produced for the Exclusive Economic Zone (EEZ) of Malaysia, which incorporates some of the latest TOPEX/POSEIDON data for the years 1992 to 1998. Local tide gauge data are used as a comparison, along with another leading Global Ocean Tide Model, Ori96. The leading diurnal and semidiurnal constituents M₂, S₂, N₂, K₁, O₁, P₁ and Q₁ are reproduced using TOPEX/POSEIDON Sea Surface Heights (SSH) in a response analysis type least squares derivation following Munk and Cartwright (1966).     

Comparison of TOPEX/POSEIDON Altimeter Derived Wind Speed and Wave Parameters with Ocean Buoy Data in the North Indian Ocean

Results of comparison exercises carried out between the state-of-the-art TOPEX/POSEIDON altimeter-derived ocean surface wind speed and ocean wave parameters (significant wave height and wave period) and those measured by a set of ocean data buoys in the North Indian Ocean are presented in this article. Altimeter-derived significant wave height values exhibited rms deviation as small as ±0.3 m, and surface wind speed of ±1.6 m/s. These results are found consistent with those found for the Pacific Ocean. For estimation of ocean wave period, the spectral moments-based semiempirical approach, earlier applied on GEOSAT data, was extended to TOPEX/POSEIDON. For this purpose, distributions of first four years of TOPEX/POSEIDON altimeter data and climatology over the North Indian Ocean were analyzed and a new set of coefficients generated for estimation of wave period. It is shown that wave periods thus estimated from TOPEX/POSEIDON data (for the subsequent two years), when compared with independent data set of ocean data buoys deployed in the North Indian Ocean, exhibit improved accuracy (rms ~ ±1.4 nos) over those determined earlier with GEOSAT data.     

Errors in recent ocean tide models: possible origin and cause

Recently, the TOPEX/POSEIDON Science Working Team has recommended the FES95.2.1 and CSR3.0 ocean tide models for reprocessing the TOPEX/POSEIDON Geophysical Data Records. Without doubt, the performance of these models, especially in the deep oceans, is excellent. However, from a comparison of these hydrodynamically consistent models with the purely empirical DW3.2 and DEOS96.1 models, it appears that FES95.2.1 and CSR3.0 are affected by basin boundary related errors which are caused by the basin-wise solution procedure of the FES ocean tide model series. In their turn, the empirical DW3.2 and DEOS96.1 models seem to suffer from significant errors in the Antarctic seas due to the seasonal growth and decay of Antarctic sea ice. Also, bathymetry-induced differences were found between the hydrodynamically consistent models and the empirical models. Concerning these differences, TOPEX/POSEIDON and ERS-1 crossover statistics unfortunately do not provide conclusive results on which models are in error.     

Tidal Correction of Altimetric Data in the Japan Sea

Satellite altimetric data have been very useful in the study of variation in the eddy field of the ocean. In order to investigate the variation in the eddy field, we have to remove tidal signals from altimetric data. However, global tidal models do not have sufficient accuracy in marginal seas such as the Japan Sea. In this study, we carried out harmonic analysis of temporal fluctuations of sea surface height data in the Japan Sea measured by TOPEX/POSEIDON. We could eliminate the tidal signals from altimetric data of TOPEX/POSEIDON and also from ERS-2 altimetric data with use of the harmonic constants derived from TOPEX/POSEIDON and tide gauge data along the coast. We draw co-tidal and co-range charts in the Japan Sea using the result of the harmonic analysis of TOPEX/POSEIDON altimetric data and tide gauge data along the coast. The results obtained turn out to be very useful for the tidal correction of altimetric data from satellite in the Japan Sea.     

Comparison of TOPEX/POSEIDON Altimeter Derived Wind Speed and Wave Parameters with Ocean Buoy Data in the North Indian Ocean

Results of comparison exercises carried out between the state-of-the-art TOPEX/POSEIDON altimeter-derived ocean surface wind speed and ocean wave parameters (significant wave height and wave period) and those measured by a set of ocean data buoys in the North Indian Ocean are presented in this article. Altimeter-derived significant wave height values exhibited rms deviation as small as - 0.3 m, and surface wind speed of - 1.6 m/s. These results are found consistent with those found for the Pacific Ocean. For estimation of ocean wave period, the spectral moments-based semiempirical approach, earlier applied on GEOSAT data, was extended to TOPEX/POSEIDON. For this purpose, distributions of first four years of TOPEX/POSEIDON altimeter data and climatology over the North Indian Ocean were analyzed and a new set of coefficients generated for estimation of wave period. It is shown that wave periods thus estimated from TOPEX/POSEIDON data (for the subsequent two years), when compared with independent data set of ocean data buoys deployed in the North Indian Ocean, exhibit improved accuracy (rms ~ - 1.4 nos) over those determined earlier with GEOSAT data.     

引入差比关系法分析西北太平洋TOPEX/POSEIDON卫星高度计测高数据

TOPEX/POSEIDON(T/P)卫星高度计数据信息中存在周期成分混淆问题.对其中的一类混淆引入差比关系方法对混淆的分潮进行分离.卫星轨道交叉点资料包括升轨和降轨资料,资料量比沿轨点资料多1倍,经分析发现:在已有为期6a多的观测资料时间序列中,在沿轨处混淆的分潮如K₁和SSA在交叉点处不再混淆,可以直接分离.因此首先对交叉点资料进行调和分析.然后由交叉点的分析结果得到分潮间的差比关系,处理到相近的沿轨点处,从而得到沿轨点的调和常数.用引入差比关系方法,对西北太平洋海区6a多的T/P卫星高度计资料进行了潮汐分析,并与沿岸及岛屿验潮站资料进行了比较,所得结果较满意.     

Tidal harmonic analysis of TOPEX/POSEIDON data in the Northwest Pacific by introducing difference-ratio relations

The difference-ratio relations are introduced to separate tidal constituents that are aliaseddue to the sampling interval and sampling span of the TOPEX/POSEIDON altimeter. It is found that some tidal constituents such as K_1 and SSA, though aliased at along track points, are not aliased at crossover points Ixrausf the data at crossover points are double those at along track points. So the harmonic analysis method can be employed directly for the analysis of time series at crossover points. Then the difference-ratio relations from crossover points are introduced to analyze the time series at along track points. The TOPEX/POSEIDON data in the Northwest Pacific are analyzed with this method. The results from this method agree well with tidal constants from tidal gauges.     

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

Seafloor mapping from high-density satellite altimetry

In this paper, we publish the results of a bathymetry survey based on the processing of satellite altimetry data. Data gathered from GEOSAT (Geodetic Mission), SEASAT, ERS-1 and TOPEX/POSEIDON satellites were processed to recover the seafloor topography over new seamounts in a test area located in the south central Pacific. We show that by processing high-density satellite altimetry data, alone or in combination with shiptrack bathymetric data, it is possible to produce full coverage bathymetric maps.     

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

本文基于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.     

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

基于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%左右。根据计算结果分析了印尼海域的潮汐特征及潮能传播规律,结果显示:爪哇海以外的印尼海域主要为不规则半日潮区;全日潮潮能主要由太平洋传入印尼海域,而半日潮潮能则是从印度洋传入印尼海域。     

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.     

Estimation of baroclinic tide energy available for deep ocean mixing based on three-dimensional global numerical simulations

The global distributions of the major semidiurnal (M₂ and S₂) and diurnal (K₁ and O₁) baroclinic tide energy are investigated using a hydrostatic sigma-coordinate numerical model. A series of numerical simulations using various horizontal grid spacings of 1/15–1/5° shows that generation of energetic baroclinic tides is restricted over representative prominent topographic features. For example, nearly half of the diurnal (K₁ and O₁) baroclinic tide energy is excited along the western boundary of the North Pacific from the Aleutian Islands down to the Indonesian Archipelago. It is also found that the rate of energy conversion from the barotropic to baroclinic tides is very sensitive to the horizontal grid spacing as well as the resolution of the model bottom topography; the conversion rate integrated over the global ocean increases exponentially as the model grid spacing is reduced. Extrapolating the calculated results in the limit of zero grid spacing yields the estimate of the global conversion rate to be 1105 GW (821, 145, 102, 53 GW for M₂, S₂, K₁, and O₁ tidal constituents, respectively). The amount of baroclinic tide energy dissipated in the open ocean below a depth of 1000 m, in particular, is estimated to be 500–600 GW, which is comparable to the mixing energy estimated by Webb and Suginohara (Nature 409:37, 2001) as needed to sustain the global overturning circulation.     

渤海潮波系统的数值模拟

利用二维非线性潮波方程组,讨论了渤黄海主要分潮(全日潮、半日潮及浅水分潮) 数值模拟中的有关问题。数值模拟中同时考虑了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°。实验结果较好地体现了渤黄海潮波系统的特征。     

2008年黄海绿潮路径的数值模拟

利用FVCOM模式拉格朗日粒子跟踪模块模拟了2008年5月到7月黄海绿潮漂移路径,根据遥感图像选取粒子跟踪的初始位置和初始时刻,模式中加入M2,S2,K1,01四个主要分潮和Quickscat每日风场资料作为驱动,模拟的粒子运动的主要路径和到达青岛近岸的时间与遥感图像对比都比较吻合。该结论进一步验证了前人关于青岛绿潮来...     

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.     

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.