Baroclinic and topographic influences on the transport in western boundary currents

Abstract

One of the central unsolved theoretical problems of the large scale ocean circulation is concerned with explaining the very large transports measured in western boundary currents such as the Gulf Stream and the Kuroshio. The only theory up to now that can explain the size of these transports is that of non-linear recirculation in which the advective terms in the momentum equations became important near the western boundary. In this paper an alternative explanation is suggested. When bottom topography and baroclinic effects are included in a wind-driven ocean model it is shown that the western boundary current can have a transport larger than that predicted from the wind stress distribution even when the nonlinear advective terms are ignored. The explanation lies in the presence of pressure torques associated with bottom topography which can contribute to the vorticity balance in the same sense as the wind stress curl.

Three numerical experiments have been carried out to explore the nature of this process using a three dimensional numerical model. The first calculation is done for a baroclinic ocean of constant depth, the second for a homogeneous ocean with an idealized continental slope topography, and the third for a baroclinic ocean with the same continental slope topography. The nature of the vorticity balance and of the circulation around closed paths is examined in each case, and it is shown that bottom pressure torques lead to enhanced transport in the western boundary current only for the baroclinic case with variable depth.     

The nonlinear spin-up of a stratified ocean

Abstract

The adjustment of a nonlinear, quasigeostrophic, stratified ocean to an impulsively applied wind stress is investigated under the assumption that barotropic advection of vortex tube length is the most important nonlinearity. The present study complements the steady state theories which have recently appeared, and extends earlier, dissipationless, linear models.

In terms of Sverdrup transport, the equation for baroclinic evolution is a forced advection-diffusion equation. Solutions of this equation subject to a “tilted disk” Ekman divergence are obtained analytically for the case of no diffusion and numerically otherwise. The similarity between the present equation and that of a forced barotropic fluid with bottom topography is shown.

Barotropic flow, which is assumed to mature instantly, can reverse the tendency for westward propagation, and thus produce regions of closed geostrophic contours. Inside these regions, dissipation, or equivalently the eddy field, plays a central role. We assume that eddy mixing effects a lateral, down-gradient diffusion of potential vorticity; hence, within the closed geostrophic contours, our model approaches a state of uniform potential vorticity. The solutions also extend the steady-state theories, which require weak diffusion, by demonstrating that homogenization occurs for moderately strong diffusion.

The evoiution of potential vorticity and the thermocline are examined, and it is shown that the adjustment time of the model is governed by dissipation, rather than baroclinic wave propagation as in linear theories. If dissipation is weak, spin-up of a nonlinear ocean may take several times that predicted by linear models, which agrees with analyses of eddy-resolving general circulation models. The inclusion of a western boundary current may accelerate this process, although dissipation will still play a central role.     

Propagation of barotropic modons over topography

Abstract

This is a broad survey of the interaction of modons with topography in a one-layer, quasigeostrophic model. Numerical simulations of modons interacting with ridges, hills, random topography and other obstacles are presented. The behavior of the modon is compared to numerical simulations of a two-point-vortex model, which proves a useful guide to the basic trajectory deflection mechanism. Under sufficiently strong but quasigeostrophically valid topographic perturbations, the modon is shown to fission into two essentially independent, oppositely-signed vortices. In the breakup of a modon near a hill it is found that the positive vortex migrates to the top of the hill. The resulting correlation between the positive vorticity trapped over the hill and the topography is in sharp contrast with the theories of turbulent flow over topography and generation of flow over topography by large scale forcing, both of which describe the development of vorticity anticorrelated with topography. A heuristic explanation of this new behavior is provided in terms of the dynamics of β bT-plane vortices. Further, it is found that a modon travelling over rough topography homogenizes the field of potential vorticity in its vicinity. This behavior is explained in terms of the induced eddy activity near the modon.     

Topographic effects on the wind‐driven ocean circulation

Abstract

This is a study of the influence of bottom topography of an ocean basin on the wind‐driven, barotropic ocean circulation. A detailed investigation is made of the role of vorticity transfer to the ocean bottom in the presence of varying topography. It is shown that the wind‐driven gyre over the topography of the North Atlantic has a transport in the western part of the basin only half of that obtained in an ocean of constant depth.     

The effect of flow at Maud Rise on the sea-ice cover – numerical experiments

 The role of seamounts in the formation and evolution of sea ice is investigated in a series of numerical experiments with a coupled sea ice–ocean model. Bottom topography, stratification and forcing are configured for the Maud Rise region in the Weddell Sea. The specific flow regime that develops at the seamount as the combined response to steady and tidal forcing consists of free and trapped waves and a vortex cap, which is caused by mean flow and tidal flow rectification. The enhanced variability through tidal motion in particular modifies the mixed layer above the seamount enough to delay and reduce sea-ice formation throughout the winter. The induced sea-ice anomaly spreads and moves westward and affects an area of several 100 000 km². Process studies reveal the complex interaction between wind, steady and periodic ocean currents: all three are required in the process of generation of the sea ice and mixed layer anomalies (mainly through tidal flow), their detachment from the topography (caused by steady oceanic flow) and the westward translation of the sea-ice anomaly (driven by the time-mean wind).     

A note on the free-surface effect on the topographically induced vorticity field in a homogeneous rotating flow

Abstract

The vorticity field induced by a topography in a homogeneous rotating fluid bounded by a free-surface is investigated. It is found that the integrated strength of the topographically bound vortex is zero and that the volume of the free-surface perturbations is equal to the volume of the topography.     

Spectral multigrid and collocation methods for barotropic nondivergent flow over irregular coastal topography

Abstract

We describe the high-resolution spectral modelling of nondivergent barotropic linearized flow over steep irregular topography. We use collocation to evaluate spatial derivatives in the barotropic vorticity equation, and a spectral multigrid technique to accelerate the iterative solution of the vorticity—stream function relation. The computational domain is a rectangular channel, which can be conformally mapped into more interesting shapes, as we also discuss. A Fourier-series representation is used in the (periodic) direction parallel to the walls of the channel, and a sine series in the cross-channel direction. For much of the paper we concentrate on the numerical techniques, though results are provided, including an application to the Bass Strait region of southeast Australia.     

A numerical study of the string function using a primitive equation ocean model

Abstract

We use results from a primitive-equation ocean numerical model (SCRUM) to test a theoretical 'string function' formulation put forward by Tyler and Käse in another article in this issue. The string function acts as a stream function for the large-scale potential energy flow under the combined beta and topographic effects. The model results verify that large-scale anomalies propagate along the string function contours with a speed correctly given by the cross-string gradient. For anomalies having a scale similar to the Rossby radius, material rates of change in the layer mass following the string velocity are balanced by material rates of change in relative vorticity following the flow velocity. It is shown that large-amplitude anomalies can be generated when wind stress is resonant with the string function configuration.     

Rapid formation of taylor columns: Obstacles against sidewalls

Abstract

The problem of topographic forcing by an obstacle against the boundary of a rotating flow is considered in various parameter regimes. The timescale for the motion is the topographic vortex-stretching time, which is inversely proportional to the background rotation rate and the fractional height of the obstacle. For slow flows this time is short compared with the advection time and the governing equation of conservation of potential vorticity is linear. The final state satisfies the non-linear equation for the advection of potential vorticity, however, and so time dependence has given a specific solution to a non-linear problem. The presence of the sidewall causes a stagnant Taylor column to be set up far more rapidly than cases with no sidewall. It is shown that viscosity and mixing arrests the inviscid evolution at some stage, thus some fluid still crosses the obstacle in the steady state. These solutions suggest that experimental results on separation obtained by Griffiths and Linden (1983) can tentatively be ascribed to entrainment and expulsion of fluid through vertical shear layers at the edge of the topography.     

Internal hydraulic control in rotating fluids—applications to oceans

Abstract

Laboratory experiments and analysis of shallow water equations in a rotating fluid show that channel flow is governed by the ratio of the width of the channel to the Rossby radius of deformation R= √[g&Delta;ρH/ρf ²]. Flows through narrow ocean openings exhibit blocking and clear evidence of hydraulic control. These imply that formulae can be derived for width, volume flux, and velocity scales of the currents. A new version of the constant potential vorticity problem is solved, and it is shown to predict volume flux within 22% of the zero potential vorticity results. Next a systematic method of predicting volume flux through ocean passages is described. Some examples are given from the Denmark Straits overflow and the flow of Antarctic Bottom Water into the western Atlantic Ocean. Two-layer flows and counter-flows with rotation in a narrow passage, the so-called lock exchange flow problem, duplicate flows at a number of important straits and openings to bays. A potential vorticity formulation is reviewed. The flows in the mouths of various bays such as Funka Bay in Hokkaido, Japan, Spencer Gulf in South Australia, and Chesapeake Bay in the United States has R < width of the mouth, and the two currents are separated by a front. The width of the front and the density difference can be predicted with good results.     

On the generation of baroclinic rossby waves in the presence of a semi-circular boundary

Abstract

An analytical model is constructed for the generation of baroclinic Rossby waves by a vorticity source in the presence of a semi-circular boundary. The vorticity source is used to represent the effect of the Agulhas retroflection to the south of Southern Africa. The displacement of the interface between the two layers of the model ocean consists of quantized waves near the coast and a train of Rossby waves drifting westward further offshore.     

On reduced-gravity flow through sills

Abstract

The hydraulic flow of a reduced-gravity fluid of non-negative potential vorticity through a sill is considered. It is shown that for any flow with a reversal of current, another, physically realisable, flow exists which is unidirectional and/or resting, and carries more flux than the original flow. Thus only non-negative flows need be considered when examining maximal hydraulic fluxes. Then, for a simple sill (one which slopes downward on the left and upward on the right, looking downstream), it is shown that zero potential vorticity flow, possibly modified by having a region of motionless fluid at its right, carries the maximum flux possible for that sill shape. This makes the calculation of maximal fluxes for a given sill considerably simpler, and examples of parabolic and V-shaped sills are computed.     

Displacement and rectification of planetary fluids

Abstract

Laboratory experiments on the decay (spin-up) of fluid motion on the β-plane are compared with theory. Under weakly dissipative conditions, some particles conserve potential vorticity during the entire decay. We also study the rectified mean flow which is produced by the lateral Reynolds Stress when a low frequency force is applied to the planetary fluid. The possible connection of effects with oceanic phenomena is briefly discussed.     

Taylor columns in horizontally sheared flow

Abstract

Solutions of the steady, inviscid, non-linear equations for the conservation of potential vorticity are presented for linearly sheared geostrophic flow over a right circular cylinder. The indeterminancy introduced by the presence of closed streamline regions is removed by requiring that the steady flow retains above topography a given fraction of that fluid initially present there, assuming the flow to have been started from rest. Those solutions which retain the largest fraction in uniform and negatively sheared streams satisfy the Ingersoll (1969) criterion (that, in the limit of vanishingly small viscosity, closed streamline regions are stagnant) and so are unaffected by Ekman pumping. These flows are set up on the advection time scale. In positively sheared flows the maximum retention solutions do not satisfy the Ingersoll criterion and thus would be slowly spun down on the far longer viscous spin-up time.

For arbitrary isolated topography, both the partial retention and Ingersoll problems are reduced to a one-dimensional non-linear integral equation and the solution of the Ingersoll problem obtained in the limit of strong positive shear. The stagnant region is symmetric about the zero velocity line and extends to infinity in the streamwise direction. Its cross-stream width is proportional to the rotation rate and fractional height occupied by the obstacle and inversely proportional to the strength of the shear, decreasing inversely as the square of distance upstream and downstream.     

Tidally generated island wakes and surface water cooling over Izu Ridge

The Kuroshio current is well known for generating cold wakes behind islands over Izu Ridge in Northwestern Pacific. Observational data from the geostationary Himawari-8 satellite for 2015–2017 revealed the occurrence of cold waters during the period when the Kuroshio current flows away from the islands. With a focus on tidal currents, this study presents an investigation of dynamical processes responsible for the formation of areas with low sea surface temperature (SST) through the adoption of a high-resolution numerical ocean model for an event that happened in July 2017. Areas with cold water emerged only when tidal currents are included in the numerical model. The model results indicate the cold surface waters are formed in the vicinity of the islands because of upwelling and vertical mixing. Qualitative features of the cold water formation for each island are found to depend on its size, topography, and ambient currents. Near Kozu Island, the tidal excursion is large enough to cause eddy shedding. These shed eddies are stirred by tidal currents to extend the surface cooling effect to wider areas. Near Hachijo Island, a persistent wake is formed by the ambient northward current. Inclusion of tidal currents destabilizes the wake, and consequently leads to the formation of a low SST area, although no clear eddy shedding is detected. The flow patterns around the islands are classified using an additional non-dimensional parameter, defined as the ratio between tidal excursion and island diameter.

     

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.     

Rectification of the wind-driven ocean circulation on the beta plane

Abstract

A nonlinear Stommel model of the ocean circulation on the beta plane, driven by a time periodic wind stress, is investigated in order to study symmetry properties of the observed time-mean ocean gyres. Due to the presence of vorticity advection terms the model will have a steady or rectified response to fluctuating wind fields. In this paper a small inverse Ekman number, “the small beta regime”, is considered. It is demonstrated that for this case all qualitative features of the residual circulation, obtained numerically by Veronis (1970). are reproduced in an analytical way. They include the dipole character of the gyre, its maximum symmetry breaking around the north-south axis for intermediate Reynolds numbers, measuring the ratio of vorticity forcing and dissipation, and the maximum residual response for intermediate forcing frequencies.     

Rotating channel flow: Control and upstream currents

Abstract

Theory and experiments are presented for critically controlled flow of a layer of inviscid rotating fluid. Flow is controlled by a level passage. For a wide upstream channel of fixed depth (i.e. constant potential vorticity) the volume flux on the right-hand wall is unaffected by passage flow. This suggests that specifying Bernoulli potential on the right-hand passage wall produces a physically well-posed condition. The specification results in one less dimensionless number than was required by previous formulations to specify flow in the controlled passage. The upstream flow needs the same number as before, so that a range of upstream conditions produce exactly the same passage flow. A laboratory study is conducted using a thin layer of water under air. This is pumped in steadily at various locations in a deep rotating upstream basin, with fluid leaving through a level passage. All currents in the upstream basin cross to the left-hand wall as the current approaches the passage over a sloping bottom. The current crosses back to the right-hand wall within the passage. Velocity profiles of currents agree reasonably well with constant potential vorticity theory. To the right of the detached upstream current is a closed gyre that connects the upstream flows (that have different patterns depending on source location) with the unique passage flows. The results suggest that gyres upstream of critically controlling passages in the ocean might serve as adjustment regions between the relatively unconstrained upstream flows and the tightly controlled passage flows.     

Inertial effects in an oceanic circulation model

Abstract

The flow in a mechanically driven thin barotropic rotating fluid system is analysed. The linear theory of Baker and Robinson (1969) is modified and extended into the non-linear regime.

An internal parameter, the “local Rossby number”, is indicative of the onset of nonlinear effects. If this parameter is 0(1) then inertial effects are as important as Coriolis accelerations in the interior of the transport-turning western boundary layer and both of its Ekman layers. The inertial effects in the Ekman layers, ignored in previous explorations of non-linear wind driven oceanic circulation, are retained here and calculated using an approximation of the Oseen type. The circulation problem is reduced to a system of scalar equations in only two independent variables; the system is valid for non-small local Rossby number provided only that the approximate total vorticity is positive.

To complete the solution for small Rossby number a boundary condition for the inertially induced transport is needed. It is found by examining the dynamics controlling this additional transport from the western boundary layer as the transport recirculates through the rest of the ocean basin. The strong constraint of total recirculation within the western boundary layer (zero net inertial transport) is derived.

The calculated primary inertial effects are in agreement with the observations of the laboratory model of Baker and Robinson (1969).

The analysis indicates the extent to which three-dimensional non-linear circulation can be reduced to a two dimensional problem.     

Instability of a geostrophic front and its energetics

Abstract

The instability of a current with a geostrophic surface density front is investigated by means of a reduced gravity model having a velocity profile with nearly uniform potential vorticity. It is shown that currents are unstable when the mean potential vorticity decreases toward the surface front at the critical point of the frontal trapped waves investigated by Paldor (1983). This instability is identical with that demonstrated by Killworth (1983) in the longwave limit.

The cross-stream component of mass flux and the rates of energy conversions among the five energy forms defined by Orlanski (1968) are also calculated. The main results are as follows, (a) The mass flux toward the surface front is positive near the front and negative around the critical point. The positive mass flux near the front does not vanish at the position of the undisturbed surface front, so that the mean position of the front moves outward and the region of the strong current spreads. (b) The potential energy of the mean flow integrated over the fluid is released through the work done by the force of the pressure gradient of the mean flow on the fluid, and is converted into the kinetic energy of the mean flow. (c) In the critical layer, the mean flow is rapidly accelerated with the growth of the unstable wave. This acceleration is caused by the rapid phase shift of the unstable wave in the critical layer.