The accuracy of state-of-the-art global barotropic tide models is assessed using bottom pressure data, coastal tide gauges, satellite altimetry, various geodetic data on Antarctic ice shelves, and independent tracked satellite orbit perturbations. Tide models under review include empirical, purely hydrodynamic (“forward”), and assimilative dynamical, i.e., constrained by observations. Ten dominant tidal constituents in the diurnal, semidiurnal, and quarter-diurnal bands are considered. Since the last major model comparison project in 1997, models have improved markedly, especially in shallow-water regions and also in the deep ocean. The root-sum-square differences between tide observations and the best models for eight major constituents are approximately 0.9, 5.0, and 6.5 cm for pelagic, shelf, and coastal conditions, respectively. Large intermodel discrepancies occur in high latitudes, but testing in those regions is impeded by the paucity of high-quality in situ tide records. Long-wavelength components of models tested by analyzing satellite laser ranging measurements suggest that several models are comparably accurate for use in precise orbit determination, but analyses of GRACE intersatellite ranging data show that all models are still imperfect on basin and subbasin scales, especially near Antarctica. For the M₂ constituent, errors in purely hydrodynamic models are now almost comparable to the 1980-era Schwiderski empirical solution, indicating marked advancement in dynamical modeling. Assessing model accuracy using tidal currents remains problematic owing to uncertainties in in situ current meter estimates and the inability to isolate the barotropic mode. Velocity tests against both acoustic tomography and current meters do confirm that assimilative models perform better than purely hydrodynamic models.

Two global ocean models ranging in horizontal resolution from 1/12° to 1/48° are used to study the space and time scales of sea surface height (SSH) signals associated with internal gravity waves (IGWs). Frequency-horizontal wavenumber SSH spectral densities are computed over seven regions of the world ocean from two simulations of the HYbrid Coordinate Ocean Model (HYCOM) and three simulations of the Massachusetts Institute of Technology general circulation model (MITgcm). High wavenumber, high-frequency SSH variance follows the predicted IGW linear dispersion curves. The realism of high-frequency motions (> 0.87??\text{cpd} ) in the models is tested through comparison of the frequency spectral density of dynamic height variance computed from the highest-resolution runs of each model (1/25° HYCOM and 1/48° MITgcm) with dynamic height variance frequency spectral density computed from nine in situ profiling instruments. These high-frequency motions are of particular interest because of their contributions to the small-scale SSH variability that will be observed on a global scale in the upcoming Surface Water and Ocean Topography (SWOT) satellite altimetry mission. The variance at supertidal frequencies can be comparable to the tidal and low-frequency variance for high wavenumbers (length scales smaller than ∼50 km), especially in the higher-resolution simulations. In the highest-resolution simulations, the high-frequency variance can be greater than the low-frequency variance at these scales.

The jets in the equatorial Pacific Ocean of a realistically forced global circulation model with a horizontal resolution of 1/12.5? cause a strong loss of phase coherence in semidiurnal internal tides that propagate equatorward from the French Polynesian Islands and Hawaii. This loss of coherence is quantified with a baroclinic energy analysis, in which the semidiurnal-band terms are separated into coherent, incoherent, and cross terms. For time scales longer than a year, the coherent energy flux approaches zero values at the equator, while the total flux is ∼500 W/m. The time variability of the incoherent energy flux is compared with the internal-tide travel-time variability, which is based on along-beam integrated phase speeds computed with the Taylor-Goldstein equation. The variability of monthly mean Taylor-Goldstein phase speeds agrees well with the phase speed variability inferred from steric sea surface height phases extracted with a plane-wave fit technique. On monthly time scales, the loss of phase coherence in the equatorward beams from the French Polynesian Islands is attributed to the time variability in the vertically sheared background flow associated with the jets and tropical instability waves. On an annual time scale, the effect of stratification variability is of equal or greater importance than the shear variability is to the loss of coherence. In the model simulations, low-frequency equatorial jets do not noticeably enhance the dissipation of the internal tide, but merely decohere and scatter it.

Motivated by the recent interest in ocean energetics, the widespread use of horizontal eddy viscosity in models, and the promise of high horizontal resolution data from the planned wide-swath satellite altimeter, this paper explores the impacts of horizontal eddy viscosity and horizontal grid resolution on geostrophic turbulence, with a particular focus on spectral kinetic energy fluxes Π(K) computed in the isotropic wavenumber (K) domain. The paper utilizes idealized two-layer quasigeostrophic (QG) models, realistic high-resolution ocean general circulation models, and present-generation gridded satellite altimeter data. Adding horizontal eddy viscosity to the QG model results in a forward cascade at smaller scales, in apparent agreement with results from present-generation altimetry. Eddy viscosity is taken to roughly represent coupling of mesoscale eddies to internal waves or to submesoscale eddies. Filtering the output of either the QG or realistic models before computing Π(K) also greatly increases the forward cascade. Such filtering mimics the smoothing inherent in the construction of present-generation gridded altimeter data. It is therefore difficult to say whether the forward cascades seen in present-generation altimeter data are due to real physics (represented here by eddy viscosity) or to insufficient horizontal resolution. The inverse cascade at larger scales remains in the models even after filtering, suggesting that its existence in the models and in altimeter data is robust. However, the magnitude of the inverse cascade is affected by filtering, suggesting that the wide-swath altimeter will allow a more accurate determination of the inverse cascade at larger scales as well as providing important constraints on smaller-scale dynamics.