Persian gulf and Oman sea tide modeling using satellite altimetry and tide gauge data (TM-IR01)

Satellite altimetry observations and tide gauge data are invaluable tools to diagnose and resolve tidal constituents over the Oceans and Seas. The aim of this study is to introduce a new purely empirical tide model named TM-IR01 in the Persian Gulf, Oman Sea, and North Indian Ocean. The observations of three altimeter sensors including TOPEX/POSIDON, JASON1, and JASON2 and 13 coastal tide gauge (TG) stations are processed and analyzed in this research. First of all, the least square spectral analysis is utilized to recover the significant tide components and consequently the amplitude and phases of the constituents are found during the tide modeling. Finally, the analysis results are interpolated into a grid of 1/4° using the Kriging method. TM-IR01 model is validated by comparing with TG stations and global tide models. It is shown that for main tidal frequencies M2, S2, K1, and O1 the root mean square error (RMSE) between TM-IR01 and TG stations results are 0.372, 0.130, 0.141, and 0.084?m, respectively, and also the RMSE between TM-IR01 and FES2004 models are 0.231, 0.087, 0.027, and 0.042?m, respectively. Validating with FES2012 and Tpxo7.2, the results obtained are close to the above values.     

山东邻海长周期分潮对深度基准面的影响分析

利用潮汐模型NAO.99Jb和FES2014确定了山东邻海的深度基准面模型并对其精度进行了评估,结果表明,NAO.99Jb模型确定的深度基准值L10的中误差为23.28 cm,FES2014模型确定的深度基准值L13的中误差为34.37 cm,长周期分潮的相对误差过大导致加入长周期分潮改正项后深度基准值中误差分别增大了11.04 cm和12.38 cm,较其他分潮对深度基准值精度的影响更明显,所以基于潮汐模型构建深度基准面模型时,长周期分潮部分必须加入实测数据改正。进一步采用山东邻海13个长期验潮站实测数据,定量地分析了长周期分潮对深度基准面确定的影响,结果表明,长周期分潮改正项的量值介于13.89~22.39 cm,平均改正值为18.03 cm,在深度基准值中占比达到15.15%。因此,长周期分潮改正对深度基准面的精确确定研究贡献较大,准确的长周期分潮模型是构建高精度深度基准面模型的基础。     

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.     

东印度洋海域最优深度基准面模型构建

利用东印度洋海域周边长期验潮站实测数据、TOPEX/Poseidon等系列卫星测高反演结果,评估了DTU10,EOT11a,FES2014,GOT4.8,OSU12和TPXO8六种全球潮汐模型精度,根据卫星测高结果给出了浅水分潮改正量和长周期分潮改正量的经验模型,又在此基础上分析并构建了研究区域精度最优的深度基准面模型。考虑到全球潮汐模型在近岸的影响因素及验潮站位置,将13个验潮站分成开阔海域与近海海域两类,与潮汐模型的对比,结果表明,DTU10和FES2014模型分别在开阔海域和近海海域精度最优。根据潮汐模型在不同分潮处的精度,如EOT11a模型在O1和K1分潮处精度较高,DTU10在N2,M2,S2和K2分潮处精度较高等,分别构建了开阔海域与近海海域的组合深度基准面模型,计算得知误差分别为11.33和20.95 cm,其精度显著提高。     

Accuracy assessment of global ocean tide models in the South China Sea using satellite altimeter and tide gauge data

In this study, to meet the need for accurate tidal prediction, the accuracy of global ocean tide models was assessed in the South China Sea (0°–26°N, 99°–121°E). Seven tide models, namely, DTU10, EOT11a, FES2014, GOT4.8, HAMTIDE12, OSU12 and TPXO8, were considered. The accuracy of eight major tidal constituents (i.e., Q₁, O₁, P₁, K₁, N₂, M₂, S₂ and K₂) were assessed for the shallow water and coastal areas based on the tidal constants derived from multi-mission satellite altimetry (TOPEX and Jason series) and tide gauge observations. The root mean square values of each constituent between satellite-derived tidal constants and tide models were found in the range of 0.72–1.90 cm in the deep ocean (depth>200 m) and 1.18–5.63 cm in shallow water area (depth<200 m). Large inter-model discrepancies were noted in the Strait of Malacca and the Taiwan Strait, which could be attributable to the complicated hydrodynamic systems and the paucity of high-quality satellite altimetry data. In coastal regions, an accuracy performance was investigated using tidal results from 37 tide gauge stations. The root sum square values were in the range of 9.35–19.11 cm, with the FES2014 model exhibiting slightly superior performance.     

Determination of Ocean Tide Loading Displacement by GPS PPP with Priori Information Constraint of NAO99b Global Ocean Tide Model

The published global ocean tide models show good agreement in deep oceans and exhibit differences in complex coastal areas, along with subsequent Ocean Tide Loading Displacement (OTLD) modeling differences. Meanwhile, OTLD parameters (amplitudes and phase lags) derived by Global Positioning System (GPS) Precise Point Positioning (PPP) approach need long time to converge to a stable state and show poor precision of S2, K1, and K2 constituents. Based on the fact that no constraint is imposed in the current kinematic solution, a new method is put forward, in which global ocean tide model predictions are taken as the priori information constraints to speed up the convergence rate and improve the accuracy of the GPS-derived OTLD parameters. First, the data of tide gauge from 01 January 2014 to 31 December 2016 are used to generate the harmonic parameters to evaluate the accuracy of six global ocean tide models and a regional ocean tide model (osu.chinesea.2010). Osu.chinesea.2010 model shows good agreement with the tide gauge results, while NAO99b model presents relatively large difference. The predictions from osu.chinesea.2010 and NAO99b model are employed as reference and the prior information, respectively. Second, continuous observations of 12 GPS sites during 2006–2013 in Hong Kong are collected to generate three dimensional OTLD amplitudes and phase lags of eight constituents using PPP with prior information constraints approach and harmonic analysis. Third, comparing the convergence time of eight constituents from PPP without and with priori information constraints approaches, the results show that the new method can significantly improve the convergence rate of OTLD amplitude estimates which obtain a certain level of stability seven years earlier than that derived by the PPP without priori information constraints. Precision of OTLD parameters derived by the new method is about 1 mm. By comparing with the precision of single PPP approach, the accuracy of eight constituents has been improved, especially for S2, K1, and K2 constituents. Finally, through comparing the different correction effects of OTLD estimates on the coordinates and their time series of the ground GPS stations, the results show that OTLD estimates derived by the new approach have similar influence as the osu.chinasea.2010 ocean tide model. The new method provides an effective means to improve the convergence and precision of the GPS-derived OTLD parameters, and achieve a similar correction as the high precision ocean tide model.     

远海航渡式水深测量水位改正方法研究

针对远海航渡式水深测量作业中的潮汐改正难题,基于全球潮汐场DTU10模型及GPS无验潮测深两种改正模式,通过潮汐场预报精度评估、验潮站实测数据比对分析以及GPS大地高计算潮汐值等多种手段,开展了大范围、长时段、单测线情况下水深测量水位改正研究,形成了一套适用性强的航渡水深测量水位改正方法与流程,为面向全球的海洋水深测量资料处理提供了潮汐、垂直基准和水位归算的方法和技术支持。     

Hydrodynamic Tidal Model of the Persian Gulf Based on Spatially Variable Bed Friction Coefficient

The aim of this study is to develop a two-dimensional hydrodynamic tidal model for the Persian Gulf (PG2017) using 2D-MIKE21 software. The advantages of present study is accounting for the spatial variation of bed friction coefficient besides a precise bathymetry together with a 23-year of combined records of satellite altimetry data. We found that the bed friction coefficient has a significant effect on sea level changes in the region under our modeling consideration. Since the tidal behavior in the northern part of the Qeshm Island is significantly different from the other parts of the Persian Gulf, to present a more accurate hydrodynamic tidal model, the Gulf is divided into two regions where the bed friction coefficient is modeled separately for each region. The root mean square value of the differences between the amplitude of dominant constituents; M2, S2, K1, and O1 derived from the PG2017 model and that of 98 altimetry and coastal tide gauge stations are respectively equal to 1.6, 1.9, 2.8, and 1.3?cm. Moreover, comparing the PG2017 model efficiency with the FES2014, OSU12, EOT11a, DTU10, and Admiralty models shows that the PG2017 model has an improvement of 22.1%, 47.2%, 43.2%, 44.2%, and 57.6% in terms of relative error, respectively.     

Absolute sea level variability of Arctic Ocean in 1993–2018 from satellite altimetry and tide gauge observations

Arctic absolute sea level variations were analyzed based on multi-mission satellite altimetry data and tide gauge observations for the period of 1993–2018. The range of linear absolute sea level trends were found ?2.00 mm/a to 6.88 mm/a excluding the central Arctic, positive trend rates were predominantly located in shallow water and coastal areas, and negative rates were located in high-latitude areas and Baffin Bay. Satellite-derived results show that the average secular absolute sea level trend was (2.53±0.42) mm/a in the Arctic region. Large differences were presented between satellite-derived and tide gauge results, which are mainly due to low satellite data coverage, uncertainties in tidal height processing and vertical land movement (VLM). The VLM rates at 11 global navigation satellite system stations around the Arctic Ocean were analyzed, among which 6 stations were tide gauge co-located, the results indicate that the absolute sea level trends after VLM corrected were of the same magnitude as satellite altimetry results. Accurately calculating VLM is the primary uncertainty in interpreting tide gauge measurements such that differences between tide gauge and satellite altimetry data are attributable generally to VLM.     

渤海风暴减水特征及其对深水航路影响的数值模拟

利用环渤海9个沿岸站近10a潮位资料分析渤海海域的风暴减水特征,结果表明:渤海年均出现50cm和100cm风暴减水分别超过30d和6d,每年的9月至翌年4月份风暴减水最为频繁;建立了一套精细化天文潮-风暴潮耦合模型用于渤海深水航路的潮位预报,各站天文潮模拟验证的平均均方根误差为18.5cm,由此计算得到航路代表点的潮汐特征值并作潮汐预报;后报模拟了近10a重大风暴减水过程,模拟与实测吻合较好,说明该耦合模型可为该航路的潮位预报提供有益参考。     

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).     

Improving the Coastal Mean Dynamic Topography by Geodetic Combination of Tide Gauge and Satellite Altimetry

Abstract

The ocean mean dynamic topography (MDT) is the surface representation of the ocean circulation. The MDT may be determined by the ocean approach, which involves temporal averaging of numerical ocean circulation model information, or by the geodetic approach, wherein the MDT is derived using the ellipsoidal height of the mean sea surface (MSS), or mean sea level (MSL) minus the geoid as the geoid. The ellipsoidal height of the MSS might be estimated either by satellite or coastal tide gauges by connecting the tide gauge datum to the Earth-centred reference frame. In this article we present a novel approach to improve the coastal MDT, where the solution is based on both satellite altimetry and tide gauge data using new set of 302 tide gauges with ellipsoidal heights through the SONEL network. The approach was evaluated for the Northeast Atlantic coast where a dense network of GNSS-surveyed tide gauges is available. The typical misfit between tide gauge and satellite or oceanographic MDT was found to be around 9?cm. This misfit was found to be mainly due to small scale geoid errors. Similarly, we found, that a single tide gauge places only weak constraints on the coastal dynamic topography.     

潮汐对多波束测深的影响及改正

讨论了潮汐对多波束测深的影响；通过分析验潮站控制及潮汐场解算的潮汐改正传统模式的特点，设计了该模式的工作流程；比较了传统改正模式与GPS无验潮模式的不同，研究了无验潮模式及陆海图拼接所涉及的垂直基准转换技术。     

Comparison of sea surface heights observed by TOPEX altimeter with sea level data at Chichijima

The sea surface heights (SSHs) observed by the TOPEX altimeter are compared with tide gauge data at Chichijima in Ogasawara (Bonin) Islands and hydrographic data taken around the islands, in order to quantitatively verify the altimeter observations and oceanic tide corrections by three tide models proposed by Cartwright and Ray (1991), Rayet al. (1994), and Maet al. (1994). First, performance of the new tide models is assessed by comparing tidal variations consisting of diurnal and semi-diurnal constituents with the tide gauge data at Chichijima. The tide model proposed by Rayet al. gives the smallest root-mean-squared (rms) difference of 2.61 cm. Errors in amplitude and phase in each tide model are evaluated by spectral analysis. The TOPEX SSHs corrected by the tide models are compared with sea level data at Chichijima. A long-term variation of a period of about 1 year is found in the residual between the SSHs and the Chichijima sea levels. This variation is also found in the difference between the dynamic height anomalies calculated from hydrographic data around the island and the Chichijima sea levels. By subtracting the long-term variation, the rms difference between the TOPEX SSHs and the Chichijima sea levels is reduced to about 4 cm and the slope of the regression line is improved to unity. The residual shows variations related to aliasing caused by incompleteness of the ocean tide correction with the repeat cycle of the altimeter observation.     

Tides in the Tongan Region of the Pacific Ocean

Tidal levels and currents in the Tongan region of the Pacific were simulated using a two-dimensional frequency-domain finite element model. The eight major diurnal and semidiurnal tidal constituents were modeled successfully, using open boundary conditions taken from a global tidal model based on the Topex/Poseidon satellite altimeter. Comparison of model results with observations from the single tide gauge site in the area were later used to adjust the boundary conditions. The validity of omitting horizontal eddy viscosity from the finite element model was checked by running an equivalent finite difference model. The results show that although the submarine Tongan ridge does not appear to trap tidal energy, there are residual tidal currents and possible recirculations which are capable of influencing biological productivity around Tonga. The model results are reduced to a simple method for predicting tidal heights in outlying areas, based only on the tidal calendar for the capital, Nuku'alofa.     

Book Reviews

The UK Hydrographic Office (UKHO)-sponsored Vertical Offshore Reference Frames (VORF) project aims to develop tidal level transformation models that are referenced to the GRS80 ellipsoid and thus compatible with GNSS positioning; in particular, heighting. Benefits include increasing the efficiency of hydrographic surveying, providing a stable consistent reference frame and enabling integration with land data in the coastal zone. Seven contemporary global ocean tide models are used to derive Lowest Astronomical Tide (LAT) surfaces which are each assessed by comparison with LAT values from the 7,389-strong UKHO tide gauge database, with the results correlated with distance from land. The proportion of truly offshore and pelagic gauges is relatively limited; however, the transition zone whereby the global ocean tide models commence to deteriorate in accuracy is evident at approximately 30km from the coast. The DTU10 model was selected as the strongest candidate overall. Subsequently, a thin plate spline method is used with the tide gauge dataset to enhance the DTU10 LAT surface in the coastal zone, creating a high resolution global LAT surface with respect to mean sea level. It is seen by cross-validation that the method may be used to predict LAT in near-shore locations with a standard error of 0.23 m.     

Tides in the Tongan Region of the Pacific Ocean

Tidal levels and currents in the Tongan region of the Pacific were simulated using a two-dimensional frequency-domain finite element model. The eight major diurnal and semidiurnal tidal constituents were modeled successfully, using open boundary conditions taken from a global tidal model based on the Topex/Poseidon satellite altimeter. Comparison of model results with observations from the single tide gauge site in the area were later used to adjust the boundary conditions. The validity of omitting horizontal eddy viscosity from the finite element model was checked by running an equivalent finite difference model. The results show that although the submarine Tongan ridge does not appear to trap tidal energy, there are residual tidal currents and possible recirculations which are capable of influencing biological productivity around Tonga. The model results are reduced to a simple method for predicting tidal heights in outlying areas, based only on the tidal calendar for the capital, Nuku'alofa.     

基于POM模式与blending同化法建立中国近海潮汐模型

利用POM海洋数值模式建立了中国近海（2°N-41°N，99°E～132°E）分辨率为5′×5′的潮汐模型，模式采用blending同化法同化了由10年TOPEX／Poseidon测高数据反演的潮汐参数与沿岸52个验潮站观测。精度分析表明建立的潮汐模型的8分潮RSS为12．5cm。     

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

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

A numerical study of sea level and current responses to Hurricane Frederic using a coastal ocean model for the Gulf of Mexico

A three-dimensional, nonlinear, primitive equation ocean general circulation model is used to study the response of the Gulf of Mexico to Hurricane Frederic. The model has free surface dynamics and a second order turbulence closure scheme for the mixed layer. Realistic coastlines, bottom topography and open boundary conditions are used in the study. The model has a vertical sigma coordinate with 18 levels, and a horizontal resolution of 0.2°×0.2° for the entire Gulf. The study focuses on hurricane generated sea level, current, and coastally trapped wave (CTW) responses of the Gulf. Time series of sea levels from U.S. coastal tide gauge stations and the numerical model simulation of sea levels and currents on the shelf are used to study sea level, current and CTW responses. Both model sea levels and observations from tide gauge stations show a westward progression of the surge as a CTW response. The results of the study of sea levels and currents indicate that CTW propagate to the west with phase speeds of 7–10 m s^–1. There is also a strong nonlinear interaction between the Loop Current and hurricane induced currents. The surface current attains a maximum of 200 cm s^–1 in the eastern Gulf. The model surface elevation at several locations is compared with tide gauge data. The current meter data at three moorings are also compared with the model currents. The model simulations show good agreement with observed data for the hurricane induced coastally trapped wave, storm surge, and current distribution in the Gulf.