An experimental investigation of blocking by partial barriers in a rotating baroclinic annulus

We present a series of experimental investigations in which a differentially-heated annulus was used to investigate the effects of topography on rotating, stratified flows with similarities to the Earth’s atmospheric or oceanic circulation. In particular, we compare and investigate blocking effects via partial mechanical barriers to previous experiments by the authors utilising azimuthally-periodic topography. The mechanical obstacle used was an isolated ridge, forming a partial barrier, employed to study the difference between partially blocked and fully unblocked flow. The topography was found to lead to the formation of bottom-trapped waves, as well as impacting the circulation at a level much higher than the top of the ridge. This produced a unique flow structure when the drifting flow and the topography interacted in the form of an “interference” regime at low Taylor number, but forming an erratic “irregular” regime at higher Taylor number. The results also showed evidence of resonant wave-triads, similar to those noted with periodic wavenumber-3 topography by Marshall and Read (Geophys. Astrophys. Fluid Dyn., 2015, 109), though the component wavenumbers of the wave-triads and their impact on the flow were found to depend on the topography in question. With periodic topography, wave-triads were found to occur between both the baroclinic and barotropic components of the zonal wavenumber-3 mode and the wavenumber-6 baroclinic component, whereas with the partial barrier two nonlinear resonant wave-triads were noted, each sharing a common wavenumber-1 mode.     

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.     

Topographically generated steady currents in barotropic turbulence

Abstract

Steady currents develop in oceanic turbulence above topography even in the absence of steady forcing. Mesoscale steady currents are correlated with mesoscale topography with anticyclonic eddies above topographic bumps, and large scale westward flows develop when β is non-zero. The relationship between those two kinds of steady currents, as well as their dependence on various parameters, is studied using a barotropic quasi-geostrophic channel model. The percentage of steady energy is found to depend on the forcing, friction and topography in a non-monotonic fashion. For example, the percentage of steady currents grows with the energy level in the linear regime (low energies) and decreases when the energy level increases in the nonlinear regime (high energies). Mesoscale steady currents are the energy source for the steady westward flow U, and therefore U is the maximum when large scale and mesoscale currents are of the same order of magnitude. This happens when the ratio S of the large scale slope βH/f ₀ and the mesoscale rms topographic slope α is of order one. U decreases for both small and large values of S.     

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.     

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.     

Rossby wave resonance in the presence of a nonlinear critical layer

Abstract

The behavior of Rossby waves on a shear flow in the presence of a nonlinear critical layer is studied, with particular emphasis on the role played by the critical layer in a Rossby wave resonance mechanism. Previous steady analyses are extended to the resonant case and it is found that the forced wave dominates the solution, provided the flow configuration is not resonant for the higher harmonics induced by the critical layer. Numerical simulations for the forced initial value problem show that the solution evolves towards the analysed steady state when conditions are resonant for the forced wave, and demonstrate some of the complications that arise when they are resonant for higher harmonics. In relating the initial value and steady problems, it is argued that the time dependent solution does not require the large mean flow distortion that Haberman (1972) found to be necessary outside the critical layer in the steady case.     

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.     

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     

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.     

Generation of waves by advected pressure fluctuations

Abstract

Because of the dispersion of surface gravity waves, the rate at which advected pressure fluctuations generate wave energy is independent of the coherence time of the pressure fluctuations provided that this time is long compared with the period of the waves. Phillips' (1957) analysis is correct insofar as it concerns wave components exactly “resonant” with the advected pressure fluctuations, but it does not deal adequately with the “bandwidth” of the resonance.     

The influence of mesoscale topography on the stability and growth rates of a two-layer model of the open ocean

Abstract

The influence of mesoscale topography on the baroclinic instability of a two-layer model of the open ocean is considered. For westward velocities in the top layer (U), and for a sinusoidal topography independent of x or longitude (a cross-stream topography), the critical value of U (U_c ) leading to instability is the same as when there is no topography. The wavelength of the unstable perturbation corresponding to U _c is shortened. For a given wavevector (k) of the perturbation the system becomes stable (as also in the absence of topography) for large values of |U|. The minimum value of the shear leading to stability is, however, significantly reduced by the topography.

For sufficiently large values of the height of the topographic features, instabilities appear which are localized within a narrow range of the shear. These instabilities are studied for a topography that depends both on x and y.

For a cross-stream topography the growth rates are somewhat smaller than those without topography and they depend only weakly on k_y . For the topographies considered here which depend both on x and y, perturbations with different values of k_y can again have roughly the same growth rate.

In the case of stable oscillations, variations in the eddy energy with very long periods are made possible by the coexistence of topographic modes with closely lying periods.     

Topographic effects on nonlinear edge waves

Abstract

Edge waves are known to give rise to beach cusps. This paper investigates the topographic feed-back upon the waves. For edge waves generated by subharmonic resonance with incident waves, the topography acts to decrease the edge wave response. As well as causing frequency detuning (Guza and Bowen, 1981) the topography can cause the scattering of edge wave energy. For synchronous waves the topographic irregularities have the opposite effect, and there can be a feed of energy into the edge waves by scattering from the incident waves.     

Nonlinearity and unsteadiness in river meandering: a review of progress in theory and modelling

River meandering has been extensively investigated. Two fundamental features to be explored in order to make further progress are nonlinearity and unsteadiness. Linear steady models have played an important role in the development of the subject but suffer from a number of limits. Moreover, rivers are not steady systems; rather their states respond to hydrologic forcing subject to seasonal oscillations, punctuated by the occurrence of flood events. We first derive a classification of river bends based on a systematic assessment of the various physical mechanisms affecting their morphodynamic equilibrium and their evolution in response to variations of hydrodynamic forcing. Using the database by Lagasse et al. ( 2004 ) we also show that natural meanders are typically mildly curved and long, i.e. such that both the centrifugal and the topographic secondary flows are weak, but they are almost invariably nonlinear. We then review some recent developments which allow us to treat analytically the flow and bed topography of mildly curved and long nonlinear bends subject to steady forcing, taking advantage of the fact that flow and bed topography in mildly curved long bends are slowly varying. Results show that nonlinearity has a number of consequences: most notably damping of the morphodynamic response and upstream shifting of the location of the nonlinear peak of the flow speed. Next we extend the latter model to the case of unsteady forcing. Results are found to depend crucially on the ratio between the flood duration and a morphodynamic timescale. It turns out that, in a channel subject to a repeated sequence of floods, the system reaches a dynamic equilibrium. We conclude the paper discussing how the present assessment relates to the debate on meander modelling of the late 1980s and suggesting what we see as promising lines of future developments.     

A nonlinear dynamo driven by rapidly rotating convection

In Kim et al. (Kim, E., Hughes, D.W. and Soward, A.M., “An investigation into high conductivity dynamo action driven by rotating convection”, Geophys. Astrophys. Fluid Dynam. 91, 303–332 ().) we investigated kinematic dynamo action driven by rapidly rotating convection in a cylindrical annulus. Here we extend this work to consider self-consistent nonlinear dynamo action in which the back-reaction of the Lorentz force on the flow is taken into account. In particular, we investigate, as a function of magnetic Prandtl number, the evolution of an initially weak magnetic field in two different types of convective flow – one chaotic and the other integrable. On saturation, the latter shows a systematic dependence on the magnetic Prandtl number whereas the former appears not to. In addition, we show how, in keeping with the findings of Cattaneo et al. (Cattaneo, F., Hughes, D.W. and Kim, E., “Suppression of chaos in a simplified nonlinear dynamo model”, Phys. Rev. Lett. 76, 2057–2060 ().), saturation of the growth of the magnetic field is brought about, for the originally chaotic flow, by a strong suppression of chaos.     

Calculation of diffusivities of passive fields in turbulent media

Abstract

The exact numerical and approximate analytical solutions of the simplest nonlinear integral equation with second order nonlinearity for the averaged Green function are presented. It is assumed that the turbulence is stationary, homogeneous, isotropic and incompressible. Numerous examples of turbulent spectra are considered (peak-like spectrum, spectra of Kolmogorov's type with different forms of “pumping” regions, stepwise spectra etc.). Special emphasis is given to investigating the case of so called “frozen” turbulence when the parameter ξ =u ₀τ/R→∞ where uτ₀,R ₀ are characteristic velocity, lifetime and space scale of turbulent pulsations, respectively. It is shown that these solutions allow us to calculate the turbulent diffusivities accurately for arbitrary spectra with any values of the parameter ξ. The results take into account the possible helicity of turbulence concerned only with scalar passive fields (number density and temperature).     

On the inadequacies of Long's model for steady two-dimensional stratified flows

Nonlinear analysis of two-dimensional steady flows with density stratification in the presence of gravity is considered. Inadequacies of Long's model for steady stratified flow over topography are explored. These include occurrence of closed streamline regions and waves propagating upstream. The usual requirements in Long's model of constant dynamic pressure and constant vertical density gradient in the upstream condition are believed to be the cause of these inadequacies. In this article, we consider a relaxation of these requirements, and also provide a systematic framework to accomplish this. As illustrations of this generalized formulation, exact solutions are given for the following two special flow configurations: the stratified flow over a barrier in an infinite channel; the stratified flow due to a line sink in an infinite channel. These solutions exhibit again closed-streamline regions as well as waves propagating upstream. The persistence of these inadequacies in the generalized Long's model appears to indicate that they are not quite consequences of the assumptions of constant dynamic pressure and constant vertical density gradient in Long's model, contrary to previous belief.

On the other hand, solutions admitted by the generalized Long's model show that departures from Long's model become small as the flow becomes more and more supercritical. They provide a nonlinear mechanism for the generation of columnar disturbances upstream of the obstacle and lead in subcritical flows to qualitatively different streamline topological patterns involving saddle points, which may describe the lee-wave-breaking process in subcritical flows and could serve as seats of turbulence in real flows. The occurrences of upstream disturbances in the presence of lee-wave-breaking activity described by the present solution are in accord with the experiments of Long (Long, R.R., “Some aspects of the flow of stratified fluids, Part 3. Continuous density gradients”, Tellus 7, 341--357 (1955)) and Davis (Davis, R.E., “The two-dimensional flow of a stratified fluid over an obstacle”, J. Fluid Mech. 36, 127–143 ()).     

On the incipient formation of bars and channels on alluvial fans

The morphodynamics of topographic expansion has been recently investigated both experimentally, by Sittoni et al., (2014) Shaw et al., (2018), and numerically Sittoni et al., 2014. Here, we study the basic mechanism that governs the evolution of topographic and expansions and explore the instability of the bottom topography under conditions of steady but spatially expanding flow. We model the expanding flow via a by configuration where water and sediments are supplied from a central hole and flow on a cone shaped surface confined by lateral walls. The governing equations are the shallow-water equations coupled with the Exner equation, written in cylindrical coordinates. We initially approach the problem analytically by considering the conditions required for the basic state, consisting of a pure radial flow and bottom profile, to lose stability to small amplitude perturbations. This analysis suggests that more than one mode may be unstable, encouraging us to extend the analysis to the nonlinear regime. We do this through numerical modeling of the full governing equations, which allows us to predict the establishment of a bar pattern whose features are similar to those experimentally observed. Two prominent features of the finite-amplitude bar pattern are (1) bar apices are distributed at a radial distance from the inflow consistent with work of Shaw et al. (2018); and (2) that the flow aspect ratio of the interbar areas remain high without provoking further instability. Both features imply that in general expansion acts to reduce bar development relative to an equivalent rectilinear flow. © 2019 John Wiley & Sons,  Ltd.     

Experiments in ocean circulation modelling

Abstract

In a laboratory model ocean, fluid in a rotating tank of varying depth is subjected to “wind-stress”, For a certain range of the parameters, Ekman number E and Rossby number R, a homogeneous fluid displays steady, westward intensified flow. For the same range of E and R, a two-layer fluid can have baroclinic instabilities. The parameter range for the various kinds of instabilities is mapped in a regime diagram. The northward transport in the western boundary current is measured as it varies with Rossby number for both homogeneous and two-layer fluid.     

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.