Abstract
During the 1990s, a large number of new tidal atlases were developed, primarily to provide accurate tidal corrections for satellite altimetry applications. During this decade, the French tidal group (FTG), led by C. Le Provost, produced a series of finite element solutions (FES) tidal atlases, among which FES2004 is the latest release, computed from the tidal hydrodynamic equations and data assimilation. The aim of this paper is to review the state of the art of tidal modelling and the progress achieved during this past decade. The first sections summarise the general FTG approach to modelling the global tides. In the following sections, we introduce the FES2004 tidal atlas and validate the model against in situ and satellite data. We demonstrate the higher accuracy of the FES2004 release compared to earlier FES tidal atlases, and we recommend its use in tidal applications. The final section focuses on the new dissipation term added to the equations, which aims to account for the conversion of barotropic energy into internal tidal energy. There is a huge improvement in the hydrodynamic tidal solution and energy budget obtained when this term is taken into account.
Similar content being viewed by others
References
Accad Y, Pekeris CL (1978) Solution of the tidal equations for the M2 and S2 tides in the world oceans from a knowledge of the tidal potential alone. Philos Trans R Soc Lond Ser A Math Phys Sci 290:235–266
Andersen OB, Egbert GD, Erofeeva SY, Ray RD (2006) Mapping nonlinear shallow-water tides: a look at the past and future (in press). DOI 10.007/s10236-006-0060-7
Arbic BK (2005) Atmospheric forcing of the oceanic semidiurnal tide. Geophys Res Lett 32:L02610, DOI 10.1029/2004GL021668
Bell TH (1975) Topographically generated internal waves in the open ocean. J Geophys Res 80:320–327
Bennett AF (1990) Inverse methods for assessing ship-of-opportunity networks and estimating circulation and winds from tropical expendable bathythermograph data. J Geophys Res 95(C9):16111–16148
Bennett AF, McIntosh PC (1982) Open ocean modelling as an inverse problem. J Phys Oceanogr 12:1004–1018
Carrère L (2003) Etude et modélisation de la réponse haute fréquence de l’Océan global aux forçages météorologiques. Ph.D. thesis, Université III Paul Sabatier
Carrère L, Lyard F (2003) Modeling the barotropic response of the global ocean to atmospheric wind and pressure forcing—comparisons with observations. Geophys Res Lett 30(6):1275, DOI 10.1029/2002GL016473
Carrère L, Le Provost C, Lyard F (2004) On the statistical stability of the M2 barotropic and baroclinic tidal characteristics from along-track TOPEX/Poseidon satellite altimetry analysis. J Geophys Res 109:C03033
Cartwright DE, Ray RD (1991) Energetics of global ocean tides from Geosat altimetry. J Geophys Res 96(C9):16897–16912
Cartwright DE, Spencer R (1988) The tides of the Atlantic Ocean. Philos Trans R Soc Lond A 324:513–563
Chabert d’Hières G, Le Provost C (1978) Atlas des composantes harmoniques de la marée dans la Manche. Ann Hydrogr 46:5–35
Egbert G, Ray RD (2000) Significant dissipation of tidal energy in the deep ocean inferred from satellite altimeter data. Nature 405:775–778
Egbert GD, Ray RD (2001) Estimates of M2 tidal energy dissipation from Topex/Poseidon altimeter data. J Geophys Res 106:22475–22502
Egbert GD, Ray RD (2003) Deviation of long period tides from equilibrium: kinematics and geostrophy. J Phys Oceanogr 33:822–839
Egbert GD, Bennett AF, Foreman GG (1994) TOPEX/POSEIDON tides estimated using a global inverse model. J Geophys Res 99:24,821–24,852
Egbert GD, Ray RD, Bills BG (2004) Numerical modeling of the global semidiurnal tide in the present day and in the last glacial maximum. J Geophys Res 109:C03003, DOI 10.1029/2003JC001973
Farrell WE (1972) Deformation of the earth by surface loads. Rev Geophys Space Phys 10:761–797
Francis O, Mazzega P (1990) Global charts of ocean tide loading effects. J Geophys Res 95(C7):11411–11424
Genco ML, Lyard F, Le Provost C (1994) The Oceanic tides in the South Atlantic Ocean. Ann Geophys 12:868–886
Greenberg D, Dupont F, Lyard F, Lynch DR, Werner FE (2006) Resolution issues in numerical models of Oceanic and coastal circulation. Cont Shelf Res (in press)
Groves GW, Reynolds RW (1975) An orthogonalized convolution method of tide prediction. J Geophys Res 80:4131–4138
Jayne SR, St. Laurent LC (2001) Parameterizing tidal dissipation over rough topography. Geophys Res Lett 28:811–814
Le Provost C, Fornerino M (1985) Tidal spectroscopy of the English Channel with a numerical model. J Phys Oceanogr 15(8):1009–1031
Le Provost C, Lyard F (1998) Energetics of the M 2 barotropic ocean tides: an estimate of bottom friction dissipation from a hydrodynamic model. Prog Oceanogr 40:37–52 (special issue on tidal science in honour of David E. Cartwright)
Le Provost C, Vincent P (1986) Some tests of precision for a finite element model of ocean tides. J Comput Phys 65:273–291
Le Provost C, Poncet A, Rougier G (1981) Numerical modelling of the harmonic constituents of the tides—application to the English Channel. J Phys Oceanogr 11(8):1123–1138
Le Provost C, Lyard F, Molines JM (1991) Improving Ocean tides predictions by using additional semidiurnal constituents from spline approximation in the frequency domain. Geophys Res Lett 18(5):845–848
Le Provost C, Bennett AF, Cartwright DE (1995) Ocean tides for and from TOPEX/POSEIDON. Science 267:639–642
Le Provost C, Lyard F, Genco ML, Rabilloud F (1998) A hydrodynamic ocean tide model improved by assimilation of a satellite altimeter-derived data set. J Geophys Res 103(C3):5513–5529
Lefevre F (2004) Modélisation des marées océaniques à l’échelle globale : assimilation de données in situ et altimétriques. Ph.D. thesis, Université III Paul Sabatier
Lefevre F, Lyard F, Le Provost C, Schrama EJO (2002) FES99: a global tide finite element solution assimilating tide gauge and altimetric information. Atmos Ocean Tech 19:1345–1356
Letellier T (2004) Etude des ondes de marée sur les plateaux continentaux. Ph.D. thesis, Université III Paul Sabatier
Lewellyn Smith SG, Young WR (2002) Conversion of the barotropic tide. J Phys Oceanogr 32:1554–1566
Lyard F (1999) Data assimilation in a wave equation: a variational representer approach for the Grenoble tidal model. J Comput Phys 149:1–31
Morozov E (1995) Semidiurnal internal wave global field. Deep-Sea Res Part 1 Oceanogr Res Pap 42:135–148
Munk WH, Cartwright DE (1966) Tidal spectroscopy and prediction. Philos Trans R Soc Lond A 259:533–581
Munk W, Wunsch C (1998) Abyssal recipes II: energetics of tidal and wind mixing. Deep-Sea Res Part 1 Oceanogr Res Pap 45:1977–2010
Ponchaut F, Lyard F, Le Provost C (2001) An analysis of the tidal signal in the WOCE Sea level dataset. J Atmos Ocean Technol 18:77–91
Ray RD (1999) A global ocean tide model from TOPEX/POSEIDON altimetry: GOT99. NASA technical memorandum 209478. Goddard space flight centre, Greenbelt
Ray RD, Egbert GD (2004) The global S1 tide. J Phys Oceanogr 34:1922–1935
Ray RD, Mitchum GT (1996) Surface manifestation of internal tides generated near Hawaii. Geophys Res Lett 23(16):2101–2104
Ray RD, Ponte RM (2003) Barometric tides from ECMWF operational analyses. Ann Geophys 21:1897–1910
Schrama E, Ray RD (1994) A preliminary tidal analysis of TOPEX/POSEIDON altimetry. J Geophys Res 99:24799–24808
Schwiderski EW (1980) Ocean tides, II, a hydrodynamic interpolation model. Mar Geod 3:219–255
Shum CK, Woodworth PL, Andersen OB, Egbert G, Francis O, King C, Klosko S, Le Provost C, Li X, Molines JM, Parke M, Ray R, Schlax M, Stammer D, Temey C, Vincent P, Wunsch C (1997) Accuracy assessment of recent ocean tide models. J Geophys Res 102(C11):25173–25194
Simmons H, St Laurent L, Jayne S, Weaver A (2004a) Tidally driven mixing in a numerical model of the ocean general circulation. Ocean Model 6:245–263
Simmons H, Hallberg R, Arbic B (2004b) Internal wave generation in a global baroclinic tide model. Deep-Sea Res Part 2 Top Stud Oceanogr 51:3043–3068, DOI 10.1016/j.dsr2.2004.09.015
Acknowledgements
Numerical simulations of the hydrodynamic and assimilation codes were performed at the Institut du Développement et des Ressources en Informatique Scientifique computational centre (Paris, France). The authors wish to thank all the people who have collaborated in the development of the FES atlases during the past 15 years, including the tidal group of the T/P mission, and the Centre National d’Etudes Spatiales for its constant support. They also thank R. Morrow for her encouragement and help in correcting earlier drafts.
Author information
Authors and Affiliations
Corresponding author
Additional information
Responsible editor: Bernard Barnier
Appendix
Appendix
1.1 Block resolution technique
As mentioned earlier, the tidal solutions are computed separately for each of the main oceanic basins. In the FES atlases prior to FES99, the continuity of the solutions along the basins’ open limits is obtained by adding some constraints in the assimilation procedures, which are performed on a global basis. This rather simple approach proved to be not totally satisfying; in particular it was responsible for some local inconsistencies (unbalanced mass or energy budget) and some undesirable side effects on the solution along the shorelines and the open ocean boundaries. To overtake this difficulty, we have developed an additional step in the tidal numerical model which allows us to retrieve the global solutions by merging together the basin-wide simulations. It is based on an approach similar to the block resolution technique, and the so-obtained solution is exactly similar to the one that would be computed on a global FE mesh (as far as the hydrodynamic model is linear).
Let us consider the global problem, obtained by merging the N-1 separate basins, with unknowns sorted by basins, except for the unknowns belonging to the open limits which are shared by at least two different basins, which are amalgamated in a Nth vector. Under this formalism, the global hydrodynamic problem of N-1 basins is equivalent to the following system:
where X i is the interior unknown vector for oceanic basin i, X N the shared open limits unknown vector for basin 1 to N-1, Y i the tidal forcing for oceanic basin i and Y N the tidal forcing along the shared open limits. The zero blocks in Eq. 30 denote the disconnection between inner nodes belonging to different basins.
The matrix Eq. 30 then yields:
The block solution technique consists in eliminating the X i in the Eq. 31.
The difficulty here is that Eq. 32 involves the computation of the inverse of A i,i , which would normally prevent us from using this formulation in a direct method. The solution to this problem would be rather to use an iterative solver instead. However, the inverse matrices are multiplied by matrices such that not all coefficients in the inverse matrices are needed. Actually, the interior unknowns of basin i (i.e. all unknowns except open boundary ones) can be separated into two groups. Group 1 contains unknowns with no direct interactions with the open boundary unknowns. Group 2 contains the “neighbours” of the open boundary unknowns. Sorting the basin interior unknowns into groups 1 and 2 allow us to redefine a simpler problem by rewriting \( A_{{N,i}} A^{{ - 1}}_{{i,i}} ,A_{{i,N}} \) as follow:
The C2,2 term is independently defined for each basin i and has the dimensions of the square of the number of computational nodes which are neighbours of at least one shared open boundary node. It can be efficiently computed using an impulse response technique from the basin, boundary condition constrained models. The basin problem with boundary conditions is given by:
where the superscript ∼ indicates that the vector/matrix are reduced to their blocks related with the unknowns of basin i. \( \widetilde{X}^{0}_{N} \) is the restriction of the \( X^{0}_{N} \) vector and represents the open boundary conditions for basin i. In our procedure, this system is solved for each basin using prior boundary conditions. We can notice that:
Rewriting Eq. 34 yields:
The most right terms of the right hand-side of Eq. 32 can be computed from the N-1 solutions of the basin problems, computed in the preliminary step:
Due to the dynamical disconnection of the interior nodes from different basins, we can infer the following equality:
The solution vector of the shared open boundary nodes is given by:
Due to the integral formulation of the wave equation, the A −, N blocks of the dynamic equations of the shared open boundary nodes can be easily obtained by adding the partial dynamic equations formed for the basin-wide problems. After Eq. 39 has been solved, the global solution is obtained by applying the following formula to each basin:
Again \( A^{{ - 1}}_{{i,i}} A_{{i,N}} \) involves only a limited number of coefficients of the inverse matrix, i.e. the coefficients related to open boundary nodes. Those necessary coefficients can be efficiently computed using an impulse response technique from the basin, boundary condition constrained models. In practice, the tidal problem is solved for each basin, and impulse response is simultaneously computed for the shared open boundary nodes and their immediate neighbours. Then the block resolution is carried out. The data assimilation uses the same approach for the backward and forward systems, except that no specific nor additional impulse response computations are needed (the impulse responses of the adjoint system are the complex conjugate of the impulse responses of the direct system). Note that in a linear problem, the prior boundary conditions are of no influence at all, and could be set to zero. In practice, specifying realistic boundary conditions is necessary for the (non-linear) dominant wave case and for computing realistic friction coefficients. This step also allows a quality control of the basin’s computation for the other (linearized) waves before passing to the computation of the global solution.
Rights and permissions
About this article
Cite this article
Lyard, F., Lefevre, F., Letellier, T. et al. Modelling the global ocean tides: modern insights from FES2004. Ocean Dynamics 56, 394–415 (2006). https://doi.org/10.1007/s10236-006-0086-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10236-006-0086-x