Ocean tides from HAMTIDE model

HAMTIDE global Phase of m2 constituent.

HAMTIDE global amplitude and phase of m₂ constituent. Copyright: UHH/CEN/AJB

Topics

Access

UNRESTRICTED:

RESTRICTED: This link to the data set is only available for a restricted user group. The data set is only accessible in CEN/MPI net or accessible from external nets with a customer account. Please contact ICDC if you would like to access this data from outside the network.

Data access via file system: /data/icdc/ocean/hamtide

to top

Description

All numerical models are subjected to errors and results are hence incomplete. Specifically with respect to tide models there may be many causes for concern, including unknown dissipation parameters [Taguchi et al., 2010] inexact bathymetry, uncertainty in initial and boundary conditions, missing physics, especially in shallow seas [Taguchi, 2004; Setiawan, 2007, see references]. Constraining tide models through geophysical data reduce this deficiency and can help to determine missing physics or uncertain parameters. However conventional data assimilation methods are afflicted with difficulty in handling the rapid growing number of altimeter data and above all in reproducing the non-linear tides in shallow water regions. The tide model HAMTIDE (Hamburg direct data Assimilation Methods for TIDEs) of the Institut für Meereskunde of the University Hamburg was developed to overcome those issues.

Fig 1: Schematic variational techniques with strong (left) and weak (right) formulation. Effective weight matrices α and β must be defined.

It is based on the generalized inverse methods for tides developed at the University of Hamburg (Zahel, 1995). The principle of the HAMTIDE is the direct minimization of the model deficiency and the inaccuracy of the recordings in a least square sense, resulting to solve the over determined algebraic equations and so called normal equation respectively. The equations are then solved by a memory saving iterative method for the given sparse matrix and the model is corrected simultaneously by inferring the physics from data. The dynamic residuals are then used for the detection of possible model errors such as bathymetry, parameterization of dissipation, loading and self-attraction (LSA) and so on [Taguchi et al., 2010, see references].

The direct method differs from usual initialization problems like quasi-geostrophic ocean model, where the initial conditions are varied as free parameters (controls) while the model is utilized as a strong constraint, and the variational technique is usually used for the solution approach (Fig. I). On the other hand the direct method puts the truth of both dynamics (weak constraint) and data in to question (Fig. I). Hence all variables x(for the whole tidal period) of the model are taken into account and adjusted till the distances of the measurements and the model-counterparts reach the minimum and the exremum of the weighted sum of dynamic and data residuals is achieved.

HAMTIDE is run with spatial resolution 7.5’. Solutions were constrained by data from the past on ongoing satellite altimeter missions as compiled by the DGFI altimeter data bank obtained from 15 years time series observed by TOPEX and Jason-1 (Bosch et al., 2009). Altimeter data were provided in form of an empirical analysis of cross-calibrated multi-mission satellite altimeter data providing constituents such as M₂, S₂, N₂, K₂, K₁, O₁, Q₁, P₁ und 2N₂.

Amplitude (AMPL) and phase (PHAS) of water elevation ζ(x,y) are denoted by a complex number like:

ζ(x,y) = RE + i * IM

with RE = AMPL * cos(PHAS) and IM = AMPL * sin(PHAS) and they are given in the following order:

Order of numbers
LON	LAT	RE	IM	AMPL	PHAS
324.000	83.625	0.4718	3.1347	3.17	81.44

Available constituents:

Astronomical semidiurnal tides:

M₂: Principal lunar,
S₂: Principal solar,
N₂: Larger elliptical lunar,
K₂: Declinational solar/lunar,
2N₂: Second-order elliptical lunar

Astronomical diurnal tides:

K₁: Principal solar/lunar,
O₁: Principal lunar,
P₁: Principal solar,
Q₁: Larger elliptical lunar

This data set is an ICDC product. Acknowledgement: This project was funded by Deutsche Forschungsgemeinschaft (SPP1257).

Last update of data set at ICDC:

to top

Parameters

Parameters of the HAMTIDE data set
Name	Unit
Amplitude	cm
Phase (Greenwich)	degree (º)

to top

Coverage, spatial and temporal resolution

Period and temporal resolution:

No period, Result from Model / Assimilationexperiment using data from a total of 15 years
No temporal resolution

Coverage and spatial resolution:

Global
Spatial resolution: 7.5' x 7.5', cartesian grid
Geographic longitude: 0.0°E to 360.0°E
Geographic latitude: -90.0°N to 90.0°N
Dimension: 1441 rows x 2881 columns
Altitude: 0.0m

Format:

NetCDF

to top

Data quality

HAMTIDE RMS differences of M2 — Fig. 2: RMS differences of M₂ tidal constants between ST102p data set and HAMTIDE with 10’ grid resolution. Diameter of blue circles is the RMS difference and proportioned to the size. Thin red lines are contour lines of 2000 and 4000 meters ocean depth [Taguchi et al., 2010].

Comparisons are made between HAMTIDE and reference pelagic data set (Fig. 2), ST102p that are provided by R. Ray [2007], also against the FES tide model [Lyard et al., 2006], taking M₂, S₂, N₂, K₂, K₁, O₁, Q₁, P₁ and 2N₂ major constituents as an example.

Note that the FES model assimilated the ST102p data set, even though, RMS differences between HAMTIDE and FES2004 for 9 constituents are within 1 mm at the ST102p.

.

RMS difference (in cm) between ST102p pelagic tidal constituents and HAMTIDE/FES2004 solutions, respectively.
Model	Resolution	M₂	S₂	N₂	K₂	K₁	O₁	P₁	Q₁	2N₂
HAMTIDE	7.5’	1.5	1.0	0.7	0.4	1.1	0.8	0.4	0.3	0.3
FES2004	-	1.4	0.9	0.7	0.5	1.0	0.7	0.4	0.3	0.3

to top

Contact

Remon Sadikni
ICDC / CEN / University of Hamburg
E-Mail: remon.sadikni"AT"uni-hamburg.de

to top

References

Literature:

Bosch B., R. Savcenko, F. Flechtner, C. Dahle, T. Mayer-Gürr, D. Stammer, E. Taguchi, and K-H. Ilk (2009), Residual ocean tide signals from satellite altimetry, GRACE gravity fields, and hydrodynamic modeling, GJI, 178(3), 1185-1192. doi :10.1111/j.1365-246X.2009.04281.x. https://doi.org/10.1111/j.1365-246X.2009.04281.x
Lyard, F., F. Lefevre, T. Letellier, and O. Francis (2006), Modelling the global ocean tides: a modern insight from FES2004, Ocean Dynamics, 56(5-6), 394-415, doi:10.1007/s10236-006-0086-x. https://doi.org/10.1007/s10236-006-0086-x
Ray R. G. (2007), ST102p data set, personnel communication
Setiawan, A. (2007), Modelling Over- and Compound Tides of the Irish and Celtic Seas using Variational Data Assimilation Methods, Thesis, University of Hamburg, Hamburg, Germany, https://ediss.sub.uni-hamburg.de/handle/ediss/1837
Taguchi, E. (2004), Inverse Modellierung nichtlinearer Flachwassergezeiten und ihre Anwendung auf ein Randmeer, Berichte aus dem Zentrum fuer Meeres- und Klimaforschung,47. https://ediss.sub.uni-hamburg.de/handle/ediss/419
Taguchi, E., Stammer, D., and Zahel, W. (2014), Inferring deep ocean tidal energy dissipation from the global high-resolution data-assimilative HAMTIDE model, J. Geophys. Res. Oceans, 119, 4573– 4592, https://doi.org/10.1002/2013JC009766
Zahel, W. (1995), Assimilating ocean tide determined data into global tidal models. J. Mar. Syst. 6, 3-13, doi:10.1016/0924-7963(94)00014-3. https://doi.org/10.1016/0924-7963(94)00014-3

to top

Data citation

Please cite the data as follows:

Taguchi, E., Stammer, D., and Zahel, W. (2014), Inferring deep ocean tidal energy dissipation from the global high-resolution data-assimilative HAMTIDE model, J. Geophys. Res. Oceans, 119, 4573– 4592, https://doi.org/10.1002/2013JC009766

Distributed in netCDF format by the Integrated Climate Data Center (ICDC), CEN, University of Hamburg, Hamburg, Germany.

to top