Calculating the ocean's mean dynamic topography from a mean sea surface and a geoid

In principle the global mean geostrophic surface circulation of the ocean can be diagnosed by subtracting a geoid from a mean sea surface (MSS). However, because the resulting mean dynamic topography (MDT) is approximately two orders of magnitude smaller than either of the constituent surfaces, and because the geoid is most naturally expressed as a spectral model while the MSS is a gridded product, in practice complications arise. Two algorithms for combining MSS and satellite-derived geoid data to determine the ocean’s mean dynamic topography (MDT) are considered in this paper: a pointwise approach, whereby the gridded geoid height field is subtracted from the gridded MSS; and a spectral approach, whereby the spherical harmonic coefficients of the geoid are subtracted from an equivalent set of coefficients representing the MSS, from which the gridded MDT is then obtained. The essential difference is that with the latter approach the MSS is truncated, a form of filtering, just as with the geoid. This ensures that errors of omission resulting from the truncation of the geoid, which are small in comparison to the geoid but large in comparison to the MDT, are matched, and therefore negated, by similar errors of omission in the MSS. The MDTs produced by both methods require additional filtering. However, the spectral MDT requires less filtering to remove noise, and therefore it retains more oceanographic information than its pointwise equivalent. The spectral method also results in a more realistic MDT at coastlines. 1. Introduction An important challenge in oceanography is the accurate determination of the ocean’s time-mean dynamic topography (MDT). If this can be achieved with sufficient accuracy for combination with the timedependent component of the dynamic topography, obtainable from altimetric data, then the resulting sum (i.e., the absolute dynamic topography) will give an accurate picture of surface geostrophic currents and ocean transports.

The influence of bottom topography on the decay of modeled Agulhas rings

The influence of a large meridional submarine ridge on the decay of Agulhas rings is investigated with a 1 and 2-layer setup of the isopycnic primitive-equation ocean model MICOM. In the single-layer case we show that the SSH decay of the ring is primarily governed by bottom friction and secondly by the radiation of Rossby waves. When a topographic ridge is present, the effect of the ridge on SSH decay and loss of tracer from the ring is negligible. However, the barotropic ring cannot pass the ridge due to energy and vorticity constraints. In the case of a two-layer ring the initial SSH decay is governed by a mixed barotropic–baroclinic instability of the ring. Again, radiation of barotropic Rossby waves is present. When the ring passes the topographic ridge, it shows a small but significant stagnation of SSH decay, agreeing with satellite altimetry observations. This is found to be due to a reduction of the growth rate of the m = 2 instability, to conversions of kinetic energy to the upper layer, and to a decrease in Rossby-wave radiation. The energy transfer is related to the fact that coherent structures in the lower layer cannot pass the steep ridge due to energy constraints. Furthermore, the loss of tracer from the ring through filamentation is less than for a ring moving over a flat bottom, related to a decrease in propagation speed of the ring. We conclude that ridges like the Walvis Ridge tend to stabilize a multi-layer ring and reduce its decay.

The spacing of orographic rainbands triggered by small-scale topography

A combination of idealized numerical simulations and analytical theory is used to investigate the spacing between convective orographic rainbands over the Coastal Range of western Oregon. The simulations, which are idealized from an observed banded precipitation event over the Coastal Range, indicate that the atmospheric response to conditionally unstable flow over the mountain ridge depends strongly on the subridge-scale topographic forcing on the windward side of the ridge. When this small-scale terrain contains only a single scale (l) of terrain variability, the band spacing is identical to l, but when a spectrum of terrain scales are simultaneously present, the band spacing ranges between 5 and 10 km, a value that is consistent with observations. Based on the simulations, an inviscid linear model is developed to provide a physical basis for understanding the scale selection of the rainbands. This analytical model, which captures the transition from lee waves upstream of the orographic cloud to moist convection within it, reveals that the spacing of orographic rainbands depends on both the projection of lee-wave energy onto the unstable cap cloud and the growth rate of unstable perturbations within the cloud. The linear model is used in tandem with numerical simulations to determine the sensitivity of the band spacing to a number of environmental and terrain-related parameters.

The triggering of orographic rainbands by small-scale topography

The triggering of convective orographic rainbands by small-scale topographic features is investigated through observations of a banded precipitation event over the Oregon Coastal Range and simulations using a cloud-resolving numerical model. A quasi-idealized simulation of the observed event reproduces the bands in the radar observations, indicating the model’s ability to capture the physics of the band-formation process. Additional idealized simulations reinforce that the bands are triggered by lee waves past small-scale topographic obstacles just upstream of the nominal leading edge of the orographic cloud. Whether a topographic obstacle in this region is able to trigger a strong rainband depends on the phase of its lee wave at cloud entry. Convective growth only occurs downstream of obstacles that give rise to lee-wave-induced displacements that create positive vertical velocity anomalies w_c and nearly zero buoyancy anomalies b_c as air parcels undergo saturation. This relationship is quantified through a simple analytic condition involving w_c, b_c, and the static stability N_m^2 of the cloud mass. Once convection is triggered, horizontal buoyancy gradients in the cross-flow direction generate circulations that align the bands parallel to the flow direction.

Approximations to the scattering of water waves by steep topography

A new method is developed for approximating the scattering of linear surface gravity waves on water of varying quiescent depth in two dimensions. A conformal mapping of the fluid domain onto a uniform rectangular strip transforms steep and discontinuous bed profiles into relatively slowly varying, smooth functions in the transformed free-surface condition. By analogy with the mild-slope approach used extensively in unmapped domains, an approximate solution of the transformed problem is sought in the form of a modulated propagating wave which is determined by solving a second-order ordinary differential equation. This can be achieved numerically, but an analytic solution in the form of a rapidly convergent infinite series is also derived and provides simple explicit formulae for the scattered wave amplitudes. Small-amplitude and slow variations in the bedform that are excluded from the mapping procedure are incorporated in the approximation by a straightforward extension of the theory. The error incurred in using the method is established by means of a rigorous numerical investigation and it is found that remarkably accurate estimates of the scattered wave amplitudes are given for a wide range of bedforms and frequencies.

The quasi-nondispersive regimes of long extratropical baroclinic rossby waves over (slowly varying) topography

Actual energy paths of long, extratropical baroclinic Rossby waves in the ocean are difficult to describe simply because they depend on the meridional-wavenumber-to-zonal-wavenumber ratio tau, a quantity that is difficult to estimate both observationally and theoretically. This paper shows, however, that this dependence is actually weak over any interval in which the zonal phase speed varies approximately linearly with tau, in which case the propagation becomes quasi-nondispersive (QND) and describable at leading order in terms of environmental conditions (i.e., topography and stratification) alone. As an example, the purely topographic case is shown to possess three main kinds of QND ray paths. The first is a topographic regime in which the rays follow approximately the contours f/h(alpha c) = a constant (alpha(c) is a near constant fixed by the strength of the stratification, f is the Coriolis parameter, and h is the ocean depth). The second and third are, respectively, "fast" and "slow" westward regimes little affected by topography and associated with the first and second bottom-pressure-compensated normal modes studied in previous work by Tailleux and McWilliams. Idealized examples show that actual rays can often be reproduced with reasonable accuracy by replacing the actual dispersion relation by its QND approximation. The topographic regime provides an upper bound ( in general a large overestimate) of the maximum latitudinal excursions of actual rays. The method presented in this paper is interesting for enabling an optimal classification of purely azimuthally dispersive wave systems into simpler idealized QND wave regimes, which helps to rationalize previous empirical findings that the ray paths of long Rossby waves in the presence of mean flow and topography often seem to be independent of the wavenumber orientation. Two important side results are to establish that the baroclinic string function regime of Tyler and K se is only valid over a tiny range of the topographic parameter and that long baroclinic Rossby waves propagating over topography do not obey any two-dimensional potential vorticity conservation principle. Given the importance of the latter principle in geophysical fluid dynamics, the lack of it in this case makes the concept of the QND regimes all the more important, for they are probably the only alternative to provide a simple and economical description of general purely azimuthally dispersive wave systems.