On steady and travelling waves over bottom topography in the model of homogeneous flow in a beta-plane channel

Abstract

A spectral low-order model is proposed in order to investigate some effects of bottom corrugation on the dynamics of forced and free Rossby waves. The analysis of the interaction between the waves and the topographic modes in the linear version of the model shows that the natural frequencies lie between the corresponding Rossby wave frequencies for a flat bottom and those applying in the “topographic limit” when the beta-effect is zero. There is a possibility of standing or eastward-travelling free waves when the integrated topograhic effect exceeds the planetary beta-effect.

The nonlinear interactions between forced waves in the presence of topography and the beta-effect give rise to a steady dynamical mode correlated to the topographic mode. The periodic solution that includes this steady wave is stable when the forcing field moves to the West with relatively large phase speed. The energy of this solution may be transferred to the steady zonal shear flow if the spatial scale of this zonal mode exceeds the scale of the directly forced large-scale dynamical mode.     

Quasigeostrophic planetary waves in a two-layer ocean with one-dimensional periodic bottom topography

Although the study of topographic effects on the Rossby waves in a stratified ocean has a long history, the wave property over a periodic bottom topography whose lateral scale is comparable to the wavelength is still not clear. The present paper treats this problem in a two-layer ocean with one-dimensional periodic bottom topography by a simple numerical method, in which no restriction on the wavelength and/or the horizontal scale of the topography is required. The dispersion diagram is obtained for a wavenumber range of [?π/L _b , π/L _b ], where L _b is the periodic length of the topography. When the topographic?β?is not negligible compared to the planetary β, the Rossby wave solutions around the wavenumbers which satisfy the resonant condition among the waves and topography disappear and separate into an infinite number of discrete modes. For convenience, each mode is numbered in order of frequency. As topographic height is increased, the high frequency barotropic Rossby wave (mode 1) becomes a topographic mode which can exist even on the f plane, and the highfrequency baroclinic mode (mode 2) becomes a surface intensified mode. Behaviors of low frequency modes are somewhat complicated. When the topographic amplitude is small, the low frequency baroclinic modes tend to be bottom trapped and the low frequency barotropic modes tend to be surface intensified. As topographic amplitude further increases, the relation between the mode number and vertical structure changes. This change can be attributed to the increase of the frequency of the topographic mode with the topographic amplitude.     

Generation of baroclinic topographic waves by a tropical cyclone impacting a low-latitude continental shelf

Numerical model experiments have been performed to analyze the low-latitude baroclinic continental shelf response to a tropical cyclone. The theory of coastally trapped waves suggests that, provided appropriate slope, latitude, stratification and wind stress, bottom-intensified topographic Rossby waves can be generated by the storm. Based on a scale analysis, the Nicaragua Shelf is chosen to study propagating topographic waves excited by a storm, and a model domain is configured with simplified but similar geometry. The model is forced with wind stress representative of a hurricane translating slowly over the region at 6 km h⁻¹. Scale analysis leads to the assumption that baroclinic Kelvin wave modes have minimal effect on the low-frequency wave motions along the slope, and coastal-trapped waves are restricted to topographic Rossby waves. Analysis of the simulated motions suggests that the shallow part of the continental slope is under the influence of barotropic topographic wave motions and at the deeper part of the slope baroclinic topographic Rossby waves dominate the low-frequency motions. Numerical solutions are in a good agreement with theoretical scale analysis. Characteristics of the simulated baroclinic waves are calculated based on linear theory of bottom-intensified topographic Rossby waves. Simulated waves have periods ranging from 153 to 203 h. The length scale of the waves is from 59 to 87 km. Analysis of energy fluxes for a fixed volume on the slope reveals predominantly along-isobath energy propagation in the direction of the group velocity of a topographic Rossby wave. Another model experiment forced with a faster translating hurricane demonstrates that fast moving tropical cyclones do not excite energetic baroclinic topographic Rossby waves. Instead, robust inertial oscillations are identified over the slope.     

Instability of planetary waves in a high vertical resolution model

Abstract

A high vertical resolution model is used to examine the instability of a baroclinic zonal flow and a finite amplitude topographically forced wave. Two families of unstable modes are found, consisting of zonally propagating most unstable modes, and stationary unstable modes. The former have time scale and spatial structure similar to baroclinic synoptic disturbances, but are localized in space due to interaction with the zonally asymmetric forcing. These modes transport heat efficiently in both the zonal and meridional directions. The second family of stationary unstable modes has characteristics of modes of low frequency variability of the atmosphere. They have time scales of 10 days and longer, and are of planetary scale with an equivalent barotropic vertical structure. The horizontal structure resembles blocking flows. They are maintained by available potential energy of the basic wave, and have large zonal heat fluxes. The results for both families of modes are interpreted in terms of an interaction between forcing and baroclinic instability to create favoured regions for eddy development. Applications to baroclinic planetary waves are also considered.     

Wave-induced topographic formstress in baroclinic channel flow

Large-scale zonal flow driven across submarine topography establishes standing Rossby waves. In the presence of stratification, the wave pattern can be represented by barotropic and baroclinic Rossby waves of mixed planetary topographic nature, which are locked to the topography. In the balance of momentum, the wave pattern manifests itself as topographic formstress. This wave-induced formstress has the net effect of braking the flow and reducing the zonal transport. Locally, it may lead to acceleration, and the parts induced by the barotropic and baroclinic waves may have opposing effects. This flow regime occurs in the circumpolar flow around Antarctica. The different roles that the wave-induced formstress plays in homogeneous and stratified flows through a zonal channel are analyzed with the BARBI (BARotropic-Baroclinic-Interaction ocean model, Olbers and Eden, J Phys Oceanogr 33:2719–2737, 2003) model. It is used in complete form and in a low-order version to clarify the different regimes. It is shown that the barotropic formstress arises by topographic locking due to viscous friction and the baroclinic one due to eddy-induced density advection. For the sinusoidal topography used in this study, the transport obeys a law in which friction and wave-induced formstress act as additive resistances, and windstress, the effect of Ekman pumping on the density stratification, and the buoyancy forcing (diapycnal mixing of the stratified water column) of the potential energy stored in the stratification act as additive forcing functions. The dependence of the resistance on the system parameters (lateral viscosity ε, lateral diffusivity κ of eddy density advection, Rossby radius λ, and topography height δ) as well as the dependence of transport on the forcing functions are determined. While the current intensity in a channel with homogeneous density decreases from the viscous flat bottom case in an inverse quadratic law ~δ ^–2 with increasing topography height and always depends on ε, a stratified system runs into a saturated state in which the transport becomes independent of δ and ε and is determined by the density diffusivity κ rather than the viscosity: κ/λ ² acts as a vertical eddy viscosity, and the transport is λ ²/κ times the applied forcing. Critical values for the topographic heights in these regimes are identified.     

Traveling planetary waves in the stratosphere

This paper reviews the theory and observations of some traveling planetary waves in the stratosphere. Two categories of waves which appear prominently in the literature are discussed: westward propagating waves of periods in the range 3–7 days (the 5-day wave) and in the range 10–20 days (the 16-day wave). Although the observations seem to indicate that these waves are waves of the Rossby type (planetary waves), the evidence is less clear regarding (1) the question of whether these waves are forced internal waves or free (resonant) external waves, and (2) the identification of the observed waves with specific theoretical waves of the Rossby type. When recent observations are compared with theory, the evidence seems to favor the notion that the 5-day and 16-day waves of longitudinal wave number 1 may be identified, respectively, with the gravest and next gravest symmetric free Rossby modes. However, the observational evidence seems to be less clear regarding the nature of the 16-day wave than the 5-day wave.     

Magnetic Rossby waves in a stably stratified layer near the surface of the Earth's outer core

Abstract

The stratification profile of the Earth's magnetofluid outer core is unknown, but there have been suggestions that its upper part may be stably stratified. Braginsky (1984) suggested that the magnetic analog of Rossby (planetary) waves in this stable layer (the ‘H’ layer) may be responsible for a portion of the short-period secular variation. In this study, we adopt a thin shell model to examine the dynamics of the H layer. The stable stratification justifies the thin-layer approximations, which greatly simplify the analysis. The governing equations are then the Laplace's tidal equations modified by the Lorentz force terms, and the magnetic induction equation. We linearize the Lorentz force in the Laplace's tidal equations and the advection term in the magnetic induction equation, assuming a zeroth order dipole field as representative of the magnetic field near the insulating core-mantle boundary. An analytical β-plane solution shows that a magnetic field can release the equatorial trapping that non-magnetic Rossby waves exhibit. A numerical solution to the full spherical equations confirms that a sufficiently strong magnetic field can break the equatorial waveguide. Both solutions are highly dissipative, which is a consequence of our necessary neglect of the induction term in comparison with the advection and diffusion terms in the magnetic induction equation in the thin-layer limit. However, were one to relax the thin-layer approximations and allow a radial dependence of the solutions, one would find magnetic Rossby waves less damped (through the inclusion of the induction term). For the magnetic field strength appropriate for the H layer, the real parts of the eigenfrequencies do not change appreciably from their non-magnetic values. We estimate a phase velocity of the lowest modes that is rather rapid compared with the core fluid speed typically presumed from the secular variation.     

On the exact shape of the horizontal profile of a topographically rectified tidal flow

Abstract

Starting from the nonlinear shallow water equations of a homogeneous rotating fluid we derive the equation describing the evolution of vorticity by a fluctuating bottom topography of small amplitude, using a multiple scale expansion in a small parameter, which is the topographic length scale relative to the tidal wave length. The exact response functions of residual vorticity for a sinusoidal bottom topography are compared with those obtained by a primitive perturbation series and by harmonic truncation, showing the former to be invalid for small topographic length scales and the latter to be only a fair approximation for vorticity produced by planetary vortex stretching. In deriving the exact shape of the horizontal residual velocity profile at a step-like break in the bottom topography, it is shown that the Lagrangian profile only exists in a strip having the width of the amplitude of the tidal excursion at both sides of the break, and that it vanishes outside that interval. Moreover, in the limit of small amplitude topography at least, it vanishes altogether for the generation mechanism by means of planetary vortex stretching. The Eulerian profile is shown to extend over twice the interval of the Lagrangian profile both for production by vortex stretching and by differential bottom friction. These finite intervals over which the residual velocity profiles exist for a step-like topography are not reproduced by harmonic truncation of the basic equation. This method gives exponentially decaying profiles, indicating spurious horizontal diffusion of vorticity. In terms of orders of magnitude, the method of harmonic truncation is reliable for residual velocity produced by vortex stretching but it overestimates the residual velocity produced by differential bottom friction by a factor 2.     

On isallobaric effect in a quasi-geostrophic topographic-drag problem

Abstract

Isallobaric effect of a slowly varying quasi-geostrophic flow represented by propagating waves may give rise to a mean steady topographic drag component which turns out to be the principal one when viscous effects are negligibly small. This drag component decreases, in contrast to the quasi-geostrophic component, when statistical properties of the topography become isotropic. When the phase velocity of the incident wave is much larger (smaller) than the phase speed of Rossby waves, the isallobaric drag becomes independent (dependent) on the sign of that velocity.     

The computation of the free barotropic oscillations of a global ocean model including friction and loading effects

A set of 61 normal modes with periods between 7.8 and 133.1 h has been calculated, using a 1° model of the global ocean, including the Arctic Ocean. The model explicitly considers frictional forces and ocean self-attraction and loading effects. The latter effects have generally been taken into account by parameterization, but for some modes the effects have also been considered fully. Due to friction, the computed eigenfrequencies are complex, exhibiting also the varying dissipative properties of the modes and their dependence on the distribution of potential and kinetic energies over the oceanic regions. In detail, gravity modes having periods less then 80 h and dominating the semi-diurnal and the diurnal tides, topographically controlled vorticity modes with periods longer than diurnal, and two planetary vorticity modes with periods of 96.8 and 119.4 h have been identified. These planetary vorticity modes have their energies distributed over Pacific, Atlantic, and Indian Oceans, while the other modes with periods longer than 80 h, as vorticity modes, have their energies concentrated on topographic structures of restricted extension. The modes are discussed with respect to their wave properties, e.g., concerning quasi-standing-wave resonances and to the appearance of Kelvin waves of different orders and trapped by different coastlines. In particular, the relevance of specific modes for the development of the fields of the most important semi-diurnal and diurnal tidal constituents is investigated.     

Open-ocean response and normal mode excitation in an eddy-resolving general circulation model

Abstract

Analysis of a two-layer, flat-bottom, steady-wind driven, eddy-resolving general circulation model reveals a distinct separation in frequency of baroclinic and barotropic motion in the region distant from the model Gulf Stream. The far-field motions at periods less (greater) than about 100 days are predominantly barotropic (baroclinic), unlike the near-field, eddy-generating, free-jet region which contains barotropic and baroclinic energy throughout the modei frequency range. The far-field barotropic energy produces a peak in the model sea-level spectra between 25 and 50 days with a magnitude comparable to energy levels observed in spectra of sea level from oceanic island tide gauges. The far-field barotropic motion is clearly composed of large-scale, resonant, barotropic normal modes drive by mesoscale activity of the turbulent, free-jet region. Oceanic mesoscale turbulence may therefore provide for planetary normal modes an excitation mechanism distinct from atmospheric forcing. The open-ocean, barotropic, model response is very similar to that of a fluctuating-wind driven model, which suggests that atmospheric and intrinsic forcing of mid-ocean eddies may be of comparable importance.     

Trapped internal waves over undular topography in a partially mixed estuary

The flow of a stratified fluid over small-scale topographic features in an estuary may generate significant internal wave activity. Lee waves and upstream influence generated at isolated topographic features have received considerable attention during the past few decades. Field surveys of a partially mixed estuary, the Rotterdam Waterway, in 1987, also showed a plethora of internal wave activity generated by isolated topography, banks and groynes. Additionally it revealed a spectacular series of resonant internal waves trapped above low-amplitude bed waves. The internal waves reached amplitudes of 3–4 m in an estuary with a mean depth of 16 m. The waves were observed during the decreasing flood tide and are thought to make a significant contribution to turbulence production and mixing. However, while stationary linear and finite amplitude theories can be used to explain the presence of these waves, it is important to further investigate their time-dependent and non-linear behaviour. With the development of advanced non-hydrostatic models it now becomes possible to further investigate these waves through numerical experimentation. This is the focus of the work presented here. The non-hydrostatic finite element numerical model FINEL3D developed by Labeur was used in the experiments presented here. The model has been shown to work well in a number of stratified flow investigations. Here, we first show that the model reproduces the field data and for idealised stationary flow scenarios that the results are in agreement with the resonant response predicted by linear theory. Then we explore the effects of non-linearity and time dependence and consider the importance of resonant internal waves for turbulence production in stratified coastal environments.Responsible Editior: Hans Burchard     

A critical analysis of the concept of planetary wave breaking

The concept of planetary wave breaking (McIntyre andPalmer, 1983; 1984) is critically reviewed. It is concluded that the wave breaking signature is not unique to any particular dynamic event in the stratosphere. Therefore, the classification of stratospheric transport events, such as wave breaking, groups fundamentally different events together. Better qualification of the wave breaking signature and a more solid theoretical basis of planetary wave breaking must be presented if the concept is to be of significant utility in describing stratospheric tracer transport.     

Topographic instability and multiple equilibria on an f-plane

Abstract

The flow of a two-layer flow in a rotating channel on an f-plane over topography with sinusoidal variation of height in a direction parallel to the flow is investigated. When the two layers flow in opposite directions a resonance is found when the topographic scale matches the free mode of the system. We examine the stability of the forced mode in the vicinity of this resonance by means of a perturbation expansion of the topographic height. Both subresonant and super-resonant instabilities are found and their equilibration is examined. For small values of the dissipation multiple equilibria are found. The topographic drag releases potential energy even when the flow is baroclinically stable.     

Resonance induced by Rossby Lee waves

Abstract

The mutual interaction of fields induced by spatially separated potential vorticity sources in a quasi-geostrophic barotropic flow is investigated using the weakly nonlinear approach. It is found that a powerful nonlinear response can be triggered by Rossby lee waves. This resonance phenomenon which dominates all other nonlinear corrections depends on certain global resonance conditions and on the change in the phase of the Rossby lee wave across the distance separating the sources. The response is particularly strong for topographic forcing possessing δ-function characterisitics.     

Relative role of subinertial and superinertial modes in the coastal long wave response forced by the landfall of a tropical cyclone

A set of numerical experiments has been performed in order to analyze the long-wave response of the coastal ocean to a translating mesoscale atmospheric cyclone approaching the coastline at a normal angle. An idealized two-slope shelf topography is chosen. The model is forced by a radially symmetric atmospheric pressure perturbation with a corresponding gradient wind field. The cyclone's translation speed, radius, and the continental shelf width are considered as parameters whose impact on the long wave period, modal structure, and amplitude is studied. Subinertial continental shelf waves (CSW) dominate the response under typical forcing conditions and on the narrower shelves. They propagate in the downstream (in the sense of Kelvin wave propagation) direction. Superinertial edge wave modes have higher free surface amplitudes and faster phase speeds than the CSW modes. While potentially more dangerous, edge waves are not as common as subinertial shelf waves because their generation requires a wide, gently sloping shelf and a storm system translating at a relatively high (∼10 m s⁻¹ or faster) speed. A relatively smaller size of an atmospheric cyclone also favors edge wave generation. Edge waves with the highest amplitude (up to 60% of the forced storm surge) propagate upstream. They are produced by a storm system with an Eulerian time scale equal to the period of a zero-mode edge wave with the wavelength of the storm spatial scale. Large amplitude edge waves were generated during Hurricane Wilma's landfall (2005) on the West Florida shelf with particularly severe flooding occurring upstream of the landfall site.     

The effects of vertical shear and stratification on stationary rossby waves

Abstract

The generation of stationary Rossby waves by sources of potential vorticity in a westerly flow is examined here in the context of a two-layer, quasi-geostrophic, β-plane model. The response in each layer consists of a combination of a barotropic Rossby wave disturbance that extends far downstream of the source, and a baroclinic disturbance which is evanescent or wave-like in character, depending on the shear and degree of stratification. Contributions from each of these modes in each layer are strongly dependent on the basic flows in each layer; the degree of stratification; and the depths of the two layers. The lower layer response is dominated by an evanescent baroclinic mode when the upper layer westerlies are much larger than those in the lower layer. In this case, weak stationary Rossby waves of large wavelengths are confined to the upper layer and the disturbance in the lower layer is confined to the source region.

Increasing the upper layer flow (with the lower layer flow fixed) increases the Rossby wavelength and decreases the amplitude. Decreasing the lower layer flow (with the upper layer flow fixed) decreases the wavelength and increases the amplitude. Stratification increases the contribution from the barotropic wave-like mode and causes the response to be confined to the lower layer.

The finite amplitude response to westerly flow over two sources of potential vorticity is also considered. In this case stationary Rossby waves induced by both sources interact to reinforce or diminish the downstream wave pattern depending on the separation distance of the sources relative to the Rossby wavelength. For fixed separation distance, enhancement of the downstreatm Rossby waves will only occur for a narrow range of flow variables and stratification.     

Effects of wave-wave and wave-mean flow interactions on the evolution of a baroclinic wave

Abstract

The weakly nonlinear evolution of a free baroclinic wave in the presence of slightly supercritical, vertically sheared zonal flow and a forced stationary wave field that consists of a single zonal scale and an arbitrary number of meridional harmonics is examined within the context of the conventional two-layer model. The presence of the (planetary-scale) stationary wave introduces zonal variations in the supercriticality and is shown to alter the growth rate and asymptotic equilibrium of the (synoptic-scale) baroclinic wave via two distinct mechanisms: The first is due to the direct interaction of the stationary wave with the shorter synoptic wave (wave-wave mechanism), and the second is due to the interaction of the synoptic wave with that portion of the mean field that is corrected by the zonally rectified stationary wave fluxes (wave-mean mechanism). These mechanisms can oppose or augment each other depending on the amplitude and spatial structure of the stationary wave field. If the stationary wave field is confined primarily to the upper (lower) layer and consists of only the gravest cross-stream mode, conditions are favorable (unfavorable) for nonzero equilibrium of the free wave.

In addition to the time dependent heat flux generated by baroclinic growth of the free wave, its interaction with a stationary wave field consisting of two or more meridional harmonics generates time dependent heat fluxes that vary with period of the free wave. However, if the stationary wave field contains several meridional harmonics of sufficiently large amplitude, the free baroclinic wave is destroyed.     

Hydromagnetic waves in rapidly rotating spherical shells generated by poloidal decay modes

Abstract

This paper presents the first attempt to examine the stability of a poloidal magnetic field in a rapidly rotating spherical shell of electrically conducting fluid. We find that a steady axisymmetric poloidal magnetic field loses its stability to a non-axisymmetric perturbation when the Elsasser number A based on the maximum strength of the field exceeds a value about 20. Comparing this with observed fields, we find that, for any reasonable estimates of the appropriate parameters in planetary interiors, our theory predicts that all planetary poloidal fields are stable, with the possible exception of Jupiter. The present study therefore provides strong support for the physical relevance of magnetic stability analysis to planetary dynamos. We find that the fluid motions driven by magnetic instabilities are characterized by a nearly two-dimensional columnar structure attempting to satisfy the Proudman-Taylor theorm. This suggests that the most rapidly growing perturbation arranges itself in such a way that the geostrophic condition is satisfied to leading order. A particularly interesting feature is that, for the most unstable mode, contours of the non-axisymmetric azimuthal flow are closely aligned with the basic axisymmetric poloidal magnetic field lines. As a result, the amplitude of the azimuthal component of the instability is smaller than or comparable with that of the poloidal component, in contrast with the instabilities generated by toroidal decay modes (Zhang and Fearn, 1994). It is shown, by examining the same system with and without fluid inertia, that fluid inertia plays a secondary role when the magnetic Taylor number T_m ? 10⁵. We find that the direction of propagation of hydromagnetic waves driven by the instability is influenced strongly by the size of the inner core.     

Baroclinic,Kelvin and inertia-gravity waves in the barostrat instability experiment

The differentially heated rotating annulus is a laboratory experiment historically designed for modelling large-scale features of the mid-latitude atmosphere. In the present study, we investigate a modified version of the classic baroclinic experiment in which a juxtaposition of convective and motionless stratified layers is created by introducing a vertical salt stratification. The thermal convective motions are suppressed in a central region at mid-depth of the rotating tank, therefore double-diffusive convection rolls can develop only in thin layers located at top and bottom, where the salt stratification is weakest. For high enough rotation rates, the baroclinic instability destabilises the flow in the top and the bottom shallow convective layers, generating cyclonic and anticyclonic eddies separated by the stable stratified layer. Thanks to this alternation of layers resembling the convective and radiative layers of stars, the planetary’s atmospheric troposphere and stratosphere or turbulent layers at the sea surface above stratified waters, this new laboratory setup is of interest for both astrophysics and geophysical sciences. More specifically, it allows to study the exchange of momentum and energy between the layers, primarily by the propagation of internal gravity waves (IGW). PIV velocity maps are used to describe the wavy flow pattern at different heights. Using a co-rotating laser and camera, the wave field is well resolved and different wave types can be found: baroclinic waves, Kelvin and Poincaré type waves. The signature of small-scale IGW can also be observed attached to the baroclinic jet. The baroclinic waves occur at the thin convectively active layer at the surface and the bottom of the tank, though decoupled they show different manifestation of nonlinear interactions. The inertial Kelvin and Poincaré waves seem to be mechanically forced. The small-scale wave trains attached to the meandering jet point to an imbalance of the large-scale flow. For the first time, the simultaneous occurrence of different wave types is reported in detail for a differentially heated rotating annulus experiment.