The quasi-geostrophic ellipsoidal vortex model

We present a simple approximate model for studying general aspects of vortex interactions in a rotating stably-stratified fluid. The model idealizes vortices by ellipsoidal volumes of uniform potential vorticity, a materially conserved quantity in an inviscid, adiabatic fluid. Each vortex thus possesses 9 degrees of freedom, 3 for the centroid and 6 for the shape and orientation. Here, we develop equations for the time evolution of these quantities for a general system of interacting vortices. An isolated ellipsoidal vortex is well known to remain ellipsoidal in a fluid with constant background rotation and uniform stratification, as considered here. However, the interaction between any two ellipsoids in general induces weak non-ellipsoidal perturbations. We develop a unique projection method, which follows directly from the Hamiltonian structure of the system, that effectively retains just the part of the interaction which preserves ellipsoidal shapes. This method does not use a moment expansion, e.g. local expansions of the flow in a Taylor series. It is in fact more general, and consequently more accurate. Comparisons of the new model with the full equations of motion prove remarkably close.Publisher PDFPeer reviewe

Balanced solutions for an ellipsoidal vortex in a rotating stratified flow

Support for this research has come from the UK Engineering and Physical Sciences Research Council (grant number EP/H001794/1).We consider the motion of a single ellipsoidal vortex with uniform potential vorticity in a rotating stratified fluid at finite Rossby number . Building on previous solutions obtained under the quasi-geostrophic approximation (at first order in ), we obtain analytical solutions for the balanced part of the flow at . These solutions capture important ageostrophic effects giving rise to an asymmetry in the evolution of cyclonic and anticyclonic vortices. Previous work has shown that, if the velocity field induced by an ellipsoidal vortex only depends linearly on spatial coordinates inside the vortex, i.e. , then the dynamics reduces markedly to a simple matrix equation. The instantaneous vortex shape and orientation are encapsulated in a symmetric matrix , which is acted upon by the flow matrix to provide the vortex evolution. Under the quasi-geostrophic approximation, the flow matrix is determined by inverting the potential vorticity to obtain the streamfunction via Poisson's equation, which has a known analytical solution depending on elliptic integrals. Here we show that higher-order balanced solutions, up to second order in the Rossby number, can also be calculated analytically. However, in this case there is a vector potential that requires the solution of three Poisson equations for each of its components. The source terms for these equations are independent of spatial coordinates within the ellipsoid, depending only on the elliptic integrals solved at the leading, quasi-geostrophic order. Unlike the quasi-geostrophic case, these source terms do not in general vanish outside the ellipsoid and have an inordinately complicated dependence on spatial coordinates. In the special case of an ellipsoid whose axes are aligned with the coordinate axes, we are able to derive these source terms and obtain the full analytical solution to the three Poisson equations. However, if one considers the homogeneous case, whereby the outer source terms are neglected, one can obtain an approximate solution having a compact matrix form analogous to the leading-order quasi-geostrophic case. This approximate solution proves to be highly accurate for the general case of an arbitrarily oriented ellipsoid, as verified through comparisons of the solutions with solutions obtained from numerical simulations of an ellipsoid using an accurate nonlinear balance model, even at moderate Rossby numbers.PostprintPeer reviewe

Influence of turbulent advection on a phytoplankton ecosystem with nonuniform carrying capacity

In this work we study a plankton ecosystem model in a turbulent flow. The plankton model we consider contains logistic growth with a spatially varying background carrying capacity and the flow dynamics are generated using the two-dimensional (2D) Navier-Stokes equations. We characterize the system in terms of a dimensionless parameter, γ TB / TF, which is the ratio of the ecosystem biological time scales TB and the flow time scales TF. We integrate this system numerically for different values of γ until the mean plankton reaches a statistically stationary state and examine how the steady-state mean and variance of plankton depends on γ. Overall we find that advection in the presence of a nonuniform background carrying capacity can lead to very different plankton distributions depending on the time scale ratio γ. For small γ the plankton distribution is very similar to the background carrying capacity field and has a mean concentration close to the mean carrying capacity. As γ increases the plankton concentration is more influenced by the advection processes. In the largest γ cases there is a homogenization of the plankton concentration and the mean plankton concentration approaches the harmonic mean, 1/K -1. We derive asymptotic approximations for the cases of small and large γ. We also look at the dependence of the power spectra exponent, β, on γ where the power spectrum of plankton is k-β. We find that the power spectra exponent closely obeys β=1+2/γ as predicted by earlier studies using simple models of chaotic advection

The Ellipsoidal Vortex: A Novel Approach to Geophysical Turbulence

We review the development of the ellipsoidal vortex model within the field of geophysical fluid dynamics. This vortex model is built on the classical potential theory of ellipsoids and applies to large-scale fluid flows, such as those found in the atmosphere and oceans, where the dynamics are strongly affected by the Earth's rotation. In this large-scale limit the governing equations reduce to the quasi-geostrophic system, where all the dynamics depends on a single scalar field, the potential vorticity, which is a dynamical marker for vortices. The solution of this system is achieved by the inversion of a Poisson equation, that in the case of an ellipsoidal vortex can be solved exactly. From this ellipsoidal solution equilibria have been determined and their stability properties have been studied. Many studies have shown that this ellipsoidal vortex model, while being conceptually simple, is an extremely powerful tool in eliciting some of the fundamental characteristics of turbulent geophysical flows

The Ellipsoidal Vortex: A Novel Approach to Geophysical Turbulence

We review the development of the ellipsoidal vortex model within the field of geophysical fluid dynamics. This vortex model is built on the classical potential theory of ellipsoids and applies to large-scale fluid flows, such as those found in the atmosphere and oceans, where the dynamics are strongly affected by the Earth's rotation. In this large-scale limit the governing equations reduce to the quasi-geostrophic system, where all the dynamics depends on a single scalar field, the potential vorticity, which is a dynamical marker for vortices. The solution of this system is achieved by the inversion of a Poisson equation, that in the case of an ellipsoidal vortex can be solved exactly. From this ellipsoidal solution equilibria have been determined and their stability properties have been studied. Many studies have shown that this ellipsoidal vortex model, while being conceptually simple, is an extremely powerful tool in eliciting some of the fundamental characteristics of turbulent geophysical flows

Resonant plankton patchiness induced by large-scale turbulent flow

Here we study how large-scale variability of oceanic plankton is affected by mesoscale turbulence in a spatially heterogeneous environment. We consider a phytoplankton-zooplankton (PZ) ecosystem model, with different types of zooplankton grazing functions, coupled to a turbulent flow described by the two-dimensional Navier-Stokes equations, representing large-scale horizontal transport in the ocean. We characterize the system using a dimensionless parameter, γ=T/T, which is the ratio of the ecosystem biological time scale T and the flow time scale T . Through numerical simulations, we examine how the PZ system depends on the time-scale ratio γ and find that the variance of both species changes significantly, with maximum phytoplankton variability at intermediate mixing rates. Through an analysis of the linearized population dynamics, we find an analytical solution based on the forced harmonic oscillator, which explains the behavior of the ecosystem, where there is resonance between the advection and the ecosystem predator-prey dynamics when the forcing time scales match the ecosystem time scales. We also examine the dependence of the power spectra on γ and find that the resonance behavior leads to different spectral slopes for phytoplankton and zooplankton, in agreement with observations

The stability of an ellipsoidal vortex in a background shear flow

We consider the motion of a single quasi-geostrophic ellipsoid of uniform potential vorticity in equilibrium with a linear background shear flow. This motion depends on four parameters: the height-to-width aspect ratio of the vortex, h/r, and three parameters characterizing the background shear flow, namely the strain rate, y, the ratio of the background rotation rate to the strain, beta, and the angle from which the shear is applied, theta. We generate the equilibria over a large range of these parameters and analyse their linear stability. For the second-order (m = 2) modes which preserve the ellipsoidal form, we are able to derive equations for the eigenmodes and growth rates. For the higher-order modes we use a numerical method to determine the full linear stability to general disturbances (m &gt; 2).Overall we find that the equilibria are stable over most of the parameter space considered, and where instability does occur the marginal instability is usually ellipsoidal. From these results, we determine the parameter values for which the vortex is most stable, and conjecture that these are the vortex characteristics which would be the most commonly observed in turbulent flows.</p

The motion of a fluid ellipsoid in a general linear background flow

The study of the motion of a fluid ellipsoid has a long and fascinating history stretching back originally to Laplace in the late 18th century. Recently, this subject has been revived in the context of geophysical fluid dynamics, where it has been shown that an ellipsoid of uniform potential vorticity remains an ellipsoid in a background flow consisting of horizontal strain, vertical shear, and uniform rotation. The object of the present work is to present a simple, appealing, and practical way of investigating the motion of an ellipsoid not just in geophysical fluid dynamics but in general. The main result is that the motion of an ellipsoid may be reduced to the evolution of a symmetric, 3 x 3 matrix, under the action of an arbitrary 3 x 3 'flow' matrix. The latter involves both the background flow, which must be linear in the Cartesian coordinates at the surface of the ellipsoid, and the self-induced flow, which was given by Laplace.The resulting simple dynamical system lends itself ideally to both numerical and analytical study. We illustrate a few examples and then present a theory for the evolution of a vortex within a slowly varying background flow. We show that a vortex may evolve quasi-adiabatically, that is, it stays close to an equilibrium form associated with the instantaneous background flow. The departure from equilibrium, on the other hand, is proportional to the rate of change of the background flow.</p

A new twist to rotating stratified turbulence

The turbulent motion of the Earth's atmosphere and oceans is hugely influenced by the effects of rotation and stratification. These effects alter the nature of turbulence profoundly from that of a homogeneous fluid. In particular, motions are dominantly horizontal, with vertical motions some three to four orders of magnitude smaller than horizontal motions. Moreover, coherent structures - vortices - are highly anisotropic, with vertical scales one to two orders of magnitude smaller than horizontal scales. And, fluid particle motions are doubly constrained: they must remain on (nearly flat) density surfaces and must retain their scalar value of 'potential vorticity'. These constraints are shown to be powerful, even in flow regimes for which rotation and stratification are not dominant effects.</p