Influence of a Background Shear Flow on Cyclone–Anticyclone Asymmetry in Ageostrophic Balanced Flows (2024)

1. Introduction

Large-scale flows in the ocean and atmosphere contain vortices. These vortices can evolve over long periods of time [1], having the potential to transport many physical properties, such as temperature and salinity, as well as biological properties over long distances [2,3]. These vortices are strongly influenced by the Earth’s rotation, referred to as the Coriolis effect, as well as the effects of density stratification [4]. As the sense of the Earth’s rotation is always the same (counterclockwise in the Northern Hemisphere, clockwise in the Southern Hemisphere), it affects cyclonic (anticlockwise) and anticyclonic (clockwise) vortices differently, leading to an asymmetry in their dynamical behavior. This asymmetry has been observed in large-scale observations [5,6] as well as in numerical simulations [7,8,9,10,11,12,13]. Despite these many studies, there is still much that is not understood about the factors that influence this asymmetry. In fact, a recent study [14] shows that there is a bias in the altimetry satellite measurements of large-scale eddies, with cyclonic eddies being under-detected; thus, there is still an important role for theoretical approaches in developing our understanding of cyclone–anticyclone asymmetry.

Here we consider an idealized model of a single vortex in the presence of a background shear flow in order to characterize the cyclonic–anticyclonic asymmetry induced by background shear flow. The model we consider is the nonlinear quasi-geostrophic model [15], a balanced model, i.e., one where all the dynamics are controlled by the evolution of the potential vorticity (PV) field, $ϖ$ , a materially conserved scalar field. The model assumes an f-plane with a fixed value for the planetary rotation rate, i.e., the Coriolis frequency, f, and the model is derived from an expansion of the full equations in the Rossby number up to the second order, where the Rossby number is defined as $R o \equiv {| ϖ |}_{m a x} / f$ . For the first order in the Rossby number, we obtain the quasi-geostrophic (QG) equations, whereas at the second order, we obtain the so-called QG + 1 equations. We use this model and study analytical solutions previously derived for an ellipsoidal vortex in the presence of a background shear flow [16,17]. We find equilibria and analyze their linear stability to determine the vortex characteristics at the margin of stability for a given background flow. As the QG solutions depend linearly on the PV, there is no dynamical difference in the behavior of cyclonic and anticyclonic vortices at this order, only the sign changes. However, the solutions up to the second order contain a linear and a quadratic term in the PV which, when combined, cause an asymmetry in dynamical behavior. Here we examine how this asymmetry changes as the background flow is changed in order to understand the impact of both horizontal strain and vertical shear.

In Section 2, we review the model of the ellipsoidal vortex and the method for obtaining equilibria and computing their linear stability. In Section 3, we present the results, showing how the vortex shape and orientation depend on the background flow and deriving an approximate model of the vortex shape characteristics that apply when the background flow is weak, in order to quantify the asymmetry between cyclonic and anticyclonic vortices. In Section 4, we discuss these results and draw our conclusions.

2. ProblemFormulation

2.1. Balanced Ellipsoidal Vortex Model

We consider the nonlinear quasi-geostrophic model [15], which is a balanced model derived from the full non-hydrostatic Oderdeck–Boussinesq equations [18]. The equations are expanded in terms of the PV-based Rossby number, $ϵ \equiv {| ϖ |}_{m a x}$ , where $ϖ$ is the PV anomaly given by:

$ϖ = \frac{ζ}{f} + \frac{1}{N^{2}} \frac{\partial b}{\partial z} + \frac{ω \cdot \nabla b}{f N^{2}},$

(1)

where $ω$ is the relative vorticity with $ζ$ being its vertical component, $b$ is the buoyancy anomaly, and $N$ and $f$ are the buoyancy and Coriolis frequencies, respectively. The other variables in the system are determined by the PV anomaly, which is materially conserved, i.e.,:

$\frac{\partial ϖ}{\partial t} + u \cdot \nabla ϖ = 0 .$

(2)

The velocity field can be expressed in terms of a vector potential, $φ \equiv (φ, ψ, ϕ)$ :

$u = - f \nabla \times φ,$

(3)

where the vector potential can be solved through an inversion problem at different orders in the Rossby number [15]. At the first order, we have the QG equations:

$\nabla^{2} ϕ^{(1)} = ϖ,$

(4)

where bracketed superscripts on field variables denote order in Rossby number (note $φ^{(1)} = ψ^{(1)} = 0$ ). Note here that the vertical coordinate z is stretched by $χ = N / f$ , whereas the buoyancy and Coriolis frequencies are considered constant here. At the next order $O (ϵ^{2})$ , we have the QG + 1 equations:

$\begin{matrix} \nabla^{2} φ^{(2)} & = & (2 / χ) (ϕ_{y y}^{(1)} ϕ_{z x}^{(1)} - ϕ_{y z}^{(1)} ϕ_{x y}^{(1)}), \end{matrix}$

(5a)

$\begin{matrix} \nabla^{2} ψ^{(2)} & = & (2 / χ) (ϕ_{x x}^{(1)} ϕ_{y z}^{(1)} - ϕ_{z x}^{(1)} ϕ_{x y}^{(1)}), \end{matrix}$

(5b)

$\begin{matrix} \nabla^{2} ϕ^{(2)} & = & {|\nabla ϕ_{z}^{(1)}|}^{2} - (\nabla^{2} ϕ^{(1)}) ϕ_{z z}^{(1)}, \end{matrix}$

(5c)

where subscripts x, y, and z on the fields denote partial differentiation.

These equations have been solved for the case of an isolated uniform ellipsoid of PV [16,19] and for the case where the vortex is embedded in a linear background shear flow [17]. We wish to revisit the latter solutions. The shape of the ellipsoid is specified by its axis half lengths a, b, and c, and by the unit vectors $\hat{a} = ({\hat{a}}_{1}, {\hat{a}}_{2}, {\hat{a}}_{3})$ , $\hat{b} = ({\hat{b}}_{1}, {\hat{b}}_{2}, {\hat{b}}_{3})$ , and $\hat{c} = ({\hat{c}}_{1}, {\hat{c}}_{2}, {\hat{c}}_{3})$ directed along these axes. When the flow field depends linearly on spatial coordinates, i.e., $u = S x$ , then the equation of motion for the ellipsoid can be written as [20]:

$\frac{d B}{d t} = S B + B S^{T},$

(6)

where $B$ and $S$ are $3 \times 3$ matrices. The matrix $B$ defines the shape and orientation of the vortex (the “shape” matrix) and is given by:

$B = M E M^{T},$

(7)

where $M$ is a rotation matrix defined by:

$M = (\begin{matrix} \hat{a} & \hat{b} & \hat{c} \end{matrix}),$

(8)

where the superscript T denotes transpose, and where matrix $E$ is a diagonal matrix:

$E = (\begin{matrix} a^{2} & 0 & 0 \\ 0 & b^{2} & 0 \\ 0 & 0 & c^{2} \end{matrix}) .$

(9)

We refer to $S$ as the “flow” matrix, and it is composed of two parts when $S = S_{v} + S_{b}$ : the self-induced motion, $S_{v}$ , and the motion induced by the external background shear flow, $S_{b}$ .

The self-induced motion of the vortex is obtained through the solutions to the equations above, where the QG $O (ϵ)$ solution is obtained by inverting Equation (4), giving [21]:

$ϕ_{v}^{(1)} = \frac{1}{2} x^{T} Φ_{v}^{(1)} x,$

(10)

where $Φ_{v}^{(1)}$ is a $3 \times 3$ symmetric matrix:

$Φ_{v}^{(1)} = M D M^{T},$

(11)

where the matrix $D$ is a diagonal matrix with:

$\begin{matrix} D_{11} = ξ_{a} & = & κ_{v} R_{D} (b^{2}, c^{2}, a^{2}), \end{matrix}$

(12a)

$\begin{matrix} D_{22} = ξ_{b} & = & κ_{v} R_{D} (c^{2}, a^{2}, b^{2}), \end{matrix}$

(12b)

$\begin{matrix} D_{33} = ξ_{c} & = & κ_{v} R_{D} (a^{2}, b^{2}, c^{2}), \end{matrix}$

(12c)

where $κ_{v} = ϖ a b c / 3$ is the vortex strength, and where $R_{D}$ is the symmetric elliptic integral of the second kind [22] given by:

$R_{D} (f, g, h) \equiv \frac{3}{2} \int_{0}^{\infty} \frac{d t}{\sqrt{(t + f) (t + g) {(t + h)}^{3}}} .$

(13)

The solution at the next order, $O (ϵ^{2})$ , is obtained by solving the QG + 1 Equation (5), whose solutions are [16]:

$\begin{matrix} φ_{v}^{(2)} & = & \frac{1}{2} x^{T} Γ_{v}^{(2)} x, Γ_{v}^{(2)} \equiv \frac{1}{χ} M H^{1} M^{T}, \end{matrix}$

(14a)

$\begin{matrix} ψ_{v}^{(2)} & = & \frac{1}{2} x^{T} Ψ_{v}^{(2)} x, Ψ_{v}^{(2)} \equiv \frac{1}{χ} M H^{2} M^{T}, \end{matrix}$

(14b)

$\begin{matrix} ϕ_{v}^{(2)} & = & \frac{1}{2} x^{T} Φ_{v}^{(2)} x, Φ_{v}^{(2)} \equiv - M H^{3} M^{T}, \end{matrix}$

(14c)

where $H^{k}$ are $3 \times 3$ symmetric matrices whose elements are given by:

$\begin{matrix} H_{11}^{k} & = & (Ω_{a b} + Ω_{c a}) ξ_{a} {\hat{a}}^{T} J^{k} \hat{a} + (ξ_{a} - Ω_{a b}) ξ_{b} {\hat{b}}^{T} J^{k} \hat{b} + (ξ_{a} - Ω_{c a}) ξ_{c} {\hat{c}}^{T} J^{k} \hat{c}, \end{matrix}$

(15a)

$\begin{matrix} H_{22}^{k} & = & (Ω_{b c} + Ω_{a b}) ξ_{b} {\hat{b}}^{T} J^{k} \hat{b} + (ξ_{b} - Ω_{b c}) ξ_{c} {\hat{c}}^{T} J^{k} \hat{c} + (ξ_{b} - Ω_{a b}) ξ_{a} {\hat{a}}^{T} J^{k} \hat{a}, \end{matrix}$

(15b)

$\begin{matrix} H_{33}^{k} & = & (Ω_{c a} + Ω_{b c}) ξ_{c} {\hat{c}}^{T} J^{k} \hat{c} + (ξ_{c} - Ω_{c a}) ξ_{a} {\hat{a}}^{T} J^{k} \hat{a} + (ξ_{c} - Ω_{b c}) ξ_{b} {\hat{b}}^{T} J^{k} \hat{b}, \end{matrix}$

(15c)

$\begin{matrix} H_{12}^{k} & = & [ξ_{a} ξ_{b} - Ω_{a b} (ξ_{a} + ξ_{b})] {\hat{a}}^{T} J^{k} \hat{b}, \end{matrix}$

2.2. Vortex Equilibria and Linear Stability

From the equation governing the evolution of the vortex, Equation (6), we can determine steady states that satisfy the following:

$S B + B S^{T} = 0 .$

(23)

This equation can be solved numerically using an iterative linear method introduced by [23]. An initial guess for the matrix $B$ is used, and then the linearized form of Equation (23) is used to obtain the next iteration. Equation (23) cannot be inverted directly, but instead one equation must be removed and conservation of volume enforced to close the equations. For the lowest strain value, the initial guess is a vortex having a circular cross section aligned with the coordinate axes, while for the initial guess, for each further increment of the strain value, we use the equilibrium obtained for the previous strain values. We apply the iterative process until the root mean square difference between the elements of the matrix $B$ at successive iterations is less than $10^{- 10}$ . If the difference between iterations is greater than 1 or if the number of iterations exceeds 10,000, then the procedure is stopped, and we assume there is no steady state for the particular parameters considered. This procedure is applied for each set of specified parameter values $(h / r, β, θ)$ and for the QG, QG + 1 cyclonic, and QG + 1 anticyclonic cases. The procedure is first applied to the smallest value of the strain rate, $| γ | = 10^{- 5}$ and, once the equilibrium is found, the strain rate is incremented by $d γ = 10^{- 5}$ . The procedure is repeated until we reach a critical turning point strain, $γ_{c}$ , i.e., the strain value beyond which there are no more steady states. In general, the critical strain is different for the QG, QG + 1 cyclonic, and QG + 1 anticyclonic equilibria.

As well as determining the equilibria, we wish to analyze their linear stability to ellipsoidal modes, i.e., the $m = 2$ mode that corresponds to changes in the vortex shape that preserve its ellipsoidal form. This mode has been shown to be the marginal mode of instability for QG ellipsoids [24]. We compute the linear stability following the method of [24], where the shape matrix is decomposed into:

$B (t) = B_{e} + \hat{B} e^{σ t},$

(24)

where $B_{e}$ is the equilibrium part, and the second term is an infinitesimal perturbation. By Taylor expansion of the flow matrix $S$ and by substituting Equation (24) into Equation (6), one can obtain an eigenvalue problem where $σ$ and $\hat{B}$ are the eigenvalue and eigenvector, respectively, and where $σ = σ_{r} + i σ_{i}$ , where the real and imaginary parts are the growth rate and frequency of the mode, respectively (see [24] for details).

3. Results

Here we compute the equilibria for a range of values of the background flow parameters $β$ and $θ$ and for a number of values of the height-to-width aspect ratio $h / r$ . For each height-to-width aspect ratio considered, we examine both like-signed interactions with $1 \leq β \leq 6$ and opposite-signed interactions with $- 4 \leq β < 1$ , while we consider values of the angle of incidence $0^{\circ} \leq θ < 90^{\circ}$ . For every case, we will compute the QG, QG + 1 cyclonic, and QG + 1 anticyclonic equilibria up to the margin of stability, and we define the magnitude of the marginal strain rate for QG, QG + 1 cyclonic, and QG + 1 anticyclonic equilibria as $γ_{m}^{q g}$ , $γ_{m}^{c y}$ , and $γ_{m}^{a c}$ , respectively.

The steady states of this system have recently been solved in [17], although the focus there was on the magnitude of the marginal and critical turning point strain and how they change with the background flow parameters. Here we wish to revisit this, although we focus on how the vortex shape and orientation change with the background flow and quantify the asymmetry between the cyclonic and the anticyclonic equilibria. We characterize the shape of the vortex by defining the horizontal aspect ratio as $δ \equiv a / b$ and the vertical aspect as $τ \equiv c / \sqrt{a b}$ , while the orientation is specified by the tilt angle $η$ (Figure 1). Using these aspect ratios, the ellipsoid semi-axis lengths are given by $a = δ^{1 / 2} / τ^{1 / 3}$ , $b = 1 / (δ^{1 / 2} τ^{1 / 3})$ , and $c = τ^{2 / 3}$ , where $a b c = 1$ , thus fixing the volume of the ellipsoid without loss of generality. Note that the vertical aspect ratio is different than the height-to-width aspect ratio $h / r$ and that it can vary with different strain values.

3.1. Vortex Characteristics as a Function of Strain

In Figure 2, we plot the horizontal and vertical aspect ratios and the tilt angle as a function of strain for QG (black), QG + 1 cyclonic (blue), and QG + 1 anticyclonic (red) equilibria for three different values of $θ$ in the case of an oblate vortex with $h / r = 0.4$ and considering the case of like-signed interactions with $β = 2$ .

For the case when $θ = 0$ , Figure 2a, only the horizontal aspect ratio varies with the strain. That is, when there is no vertical shear, there is no tilting of the vortex, and the vertical aspect ratio is equal to the height-to-width aspect ratio. On the other hand, the horizontal aspect ratio decreases with increasing strain, with the cyclonic values being less than the anticyclonic values for a given strain value. However, as the marginal strain rate of the anticyclonic equilibria are higher than those of the cyclones, the marginal horizontal aspect ratio ends up being the same for anticyclones and cyclones. For the cases $θ = 30^{\circ}$ and $θ = 60^{\circ}$ in Figure 2b and Figure 2c, respectively, the equilibria tilt angle increases as the strain increases. As the equilibria increasingly tilt with the strain, the vertical aspect ratio $τ$ decreases. For a given strain rate, the deformation of the cyclonic vortex is greater than the anticyclonic case. However, the marginal strain is greater for the anticyclonic vortex, and so at the marginal strain, the anticyclone is more tilted and has a smaller vertical aspect ratio than the cyclone.

In Figure 3, we plot the parameters for the case of a prolate vortex with $h / r = 2.5$ . When there is no vertical shear, $θ = 0$ (Figure 3a), the behavior is similar to that which was seen for the oblate case, although in this case, at the marginal strain, there is stronger asymmetry of the horizontal aspect ratio, with the horizontal aspect ratio of the cyclonic equilibrium being much less than the anticyclonic equilibrium.

When there is vertical shear, Figure 3b,c, there is even stronger asymmetry in the horizontal aspect ratio, with the cyclonic case decreasing with strain, while the anticyclonic case initially decreases but eventually increases and becomes larger than its initial value. Unlike the oblate case, the vertical aspect ratio increases with the strain rate, while the vortex has a much larger degree of tilting. Additionally, conversely to the oblate case, for a given value of the strain rate, the anticyclonic vortex is more deformed and has a greater tilt angle than the cyclonic case.

Next, we consider the case of opposite-signed interactions with $β = 0$ and show the cases $h / r = 0.4$ and $h / r = 2.5$ in Figure 4 and Figure 5, respectively.

One significant difference between the opposite and like-signed cases is that the opposite-signed horizontal aspect ratio increases with the strain rate for both oblate and prolate cases, whereas this only occurs for anticyclonic prolate equilibria when there is vertical shear. Overall, the opposite-signed equilibria are more deformed than the like-signed equivalent cases, probably linked to the fact that the marginal strain values for opposite-signed interactions are much larger those for like-signed interactions. However, the overall asymmetry between the cyclonic and anticyclonic equilibria for the opposite-signed interactions is reduced in nearly all cases, except for the prolate case with $θ = 60^{\circ}$ (Figure 5c).

3.2. Asymmetry at the Marginal Strain

Now we examine how the shape characteristics of the cyclonic and anticyclonic equilibria at the marginal strain rate change with the background flow parameters. Similar to how we refer to the marginal strain rates for QG, QG + 1 cyclonic, and QG + 1 anticyclonic equilibria, we define the marginal shape characteristics for the QG, QG + 1 cyclonic, and QG + 1 anticyclonic respectively as: $δ_{m}^{q g}$ , $δ_{m}^{c y}$ , and $δ_{m}^{a c}$ for the marginal horizontal aspect ratio; $τ_{m}^{q g}$ , $τ_{m}^{c y}$ , and $τ_{m}^{a c}$ for the marginal vertical aspect ratio; and $η_{m}^{q g}$ , $η_{m}^{c y}$ , and $η_{m}^{a c}$ for the marginal tilt angle. First, we will consider the dependence on the vertical offset $θ$ , beginning with the case of like-signed interactions in Figure 6.

We plot the marginal strain rate $γ_{m}$ and the shape and orientation variables at the marginal strain rate as a function of $θ$ for the cases when $β = 2$ and for (a) $h / r = 0.4$ , (b) $h / r = 0.8$ , (c) $h / r = 1$ , (d) $h / r = 1.25$ , and (e) $h / r = 2.5$ . The marginal strain rate shows similar dependence on $θ$ , with the the largest values occurring at $θ = 90^{\circ}$ . However, for the oblate cases ( $h / r \leq 1$ ), $γ_{m}^{a c} > γ_{m}^{c y}$ , meaning that the anticyclonic oblate vortices are more stable, while the reverse is the case for the prolate cases ( $h / r > 1$ ). The horizontal aspect ratio $δ_{m}$ also tends towards larger values for large $θ$ , although the dependence on $θ$ changes with $h / r$ . In particular, the most prolate case $h / r = 2.5$ differs from the other cases, as it has the strongest cyclone–anticyclone asymmetry at large $θ$ , whereas for the other values of $h / r$ , the asymmetry is strongest for low $θ$ . For all aspect ratios, $δ_{m}^{a c} > δ_{m}^{c y}$ , and the value of $δ$ is closest to unity for $θ = 90^{\circ}$ , i.e., for the largest value of vertical shear, the equilibrium horizontal cross section becomes more circular. For the vertical aspect ratio, there is similar behavior for the cases where $h / r < 2.5$ , with $τ_{m}$ decreasing with $θ$ until it reaches $θ \approx 70^{\circ}$ , at which point it increases again, returning to near the initial value, with $τ_{m}^{a c} < τ_{m}^{c y}$ . However, for $h / r = 2.5$ (Figure 6e), $τ_{m}$ increases with $θ$ , reaching a maximum at about $θ \approx 70^{\circ}$ , before returning to the initial value. In this case, $τ_{m}^{a c} > τ_{m}^{c y}$ . The marginal tilting angle $η_{m}$ for all $h / r$ initially increases with $θ$ , although this increase is more abrupt for $h / r$ close to unity. For most values of $h / r$ , the tilting angle decreases at large $θ$ , except for the case $h / r = 1$ , where the value remains high for the largest $θ$ . Overall, looking at the asymmetry in the shape and orientation, it is clear that the anticyclonic equilibria undergo greater deformation than the cyclonic ones. The greatest difference between the anticyclonic and cyclonic marginal shape most often occurs at $θ \approx 70^{\circ}$ , with the exception of when $h / r$ is close to unity (Figure 6b,c,e), with $δ_{m}$ and $η_{m}$ having the largest asymmetry when $θ = 0^{\circ}$ .

In Figure 7, we consider the dependence on $θ$ for the opposite-signed case with $β = 0$ . Here, there is very different behavior for the most oblate case, $h / r = 0.4$ , and for the other cases. For the strongly oblate case in Figure 7a, we see that $γ_{c}^{a c} > γ_{c}^{c y}$ , although there are very few differences between the shape and orientation of the cyclonic and anticyclonic equilibria. The horizontal aspect ratio decreases with $θ$ , while the vertical aspect ratio initially decreases to a minimum at $θ \approx 60^{\circ}$ before increasing back to its original value. Conversely, the tilt angle increases to a maximum of $10^{\circ}$ before returning to zero. For $h / r > 0.4$ , the behavior of all the variables looks dramatically different from the most oblate case, with a kink appearing in the dependence at $θ \approx 60^{\circ}$ . This kink has been seen before in [17,24] and is associated with a swap in the horizontal axes a and b of the equilibria when the shape matrix component $B_{11} = B_{22}$ . It marks a transition in the behavior of the variables, with the vertical aspect ratio and the tilt angle reaching a maximum at this point before decreasing, while for the horizontal aspect ratio, it is close to the point where we go from $δ_{m}^{a c} < δ_{m}^{c y}$ to $δ_{m}^{a c} > δ_{m}^{c y}$ . For $δ_{m}$ , interestingly, when there is no vertical shear, $θ = 0^{\circ}$ , and the horizontal cross section is more deformed, with $δ_{m}$ much larger than 1; however, as the vertical shear increases, the horizontal cross section becomes more circular. In the case of the vertical aspect ratio and the tilt angle, the greatest asymmetry between the cyclonic and anticyclonic equilibria occurs where the kink occurs, whereas for the horizontal aspect ratio, the strongest asymmetry occurs for $θ = 0^{\circ}$ . Overall, the asymmetry is largest for the prolate cases ( $h / r > 1$ ).

Next, we look at how the shape characteristics of the equilibria at the marginal strain rate change with the parameter $β$ . In Figure 8, we plot the equilibria conditions at the marginal strain as a function of $β$ for the case of like-signed interactions.

For all values of $h / r$ , the marginal strain $γ_{m}$ decreases with increasing $β$ , although there is a difference between the oblate and prolate cases, with the oblate anticyclonic equilibria being more stable than the cyclonic equilibria, $γ_{m}^{a c} > γ_{m}^{c y}$ , while the reverse is true for the prolate cases. The marginal horizontal aspect ratios have similar dependence on $β$ for all values of $h / r$ , with $δ_{m}$ increasing with $β$ , moving towards unity; the only exception is $δ_{m}^{a c}$ for the most prolate case, where the value decreases as $β$ increases. In all cases, $δ_{m}^{a c} > δ_{m}^{c y}$ , with the strongest asymmetry occurring for $h / r = 2.5$ . Similarly, for the vertical aspect ratio, we see different behavior for the most prolate case, with $τ_{m}$ decreasing with $β$ , while for all the other values of $h / r$ , $τ_{m}$ increases with $β$ . In all cases, $τ_{m}$ approaches $h / r$ as $β$ increases. For the tilting angle $η_{m}$ , there seems to be a separation between the cases with $h / r$ close to unity and the extremely oblate and prolate cases, with the most oblate and prolate values decreasing with $β$ , while the other cases increase with $β$ . Overall, the asymmetry tends to increase with $β$ , with only the most oblate and prolate cases defying this trend for $δ_{m}$ .

In Figure 9, we look at the dependence on $β$ for opposite-signed interactions.

Here, the marginal strain behaves the same for all values of $h / r$ , with $γ_{m}$ increasing rapidly as $β$ goes to 1. These larger marginal strain values on the approach to $β = 1$ mark extreme deformation of the horizontal cross section of the equilibria, with $δ_{m}$ increasing dramatically for all values of $h / r$ . Instead, the behavior of the marginal vertical aspect ratio $τ_{m}$ changes with $h / r$ , with the most oblate cases (Figure 9a,b) decreasing with $β$ , whereas for the other values of $h / r$ , it increases with $β$ , except for $τ_{m}^{a c}$ for $h / r = 2.5$ (Figure 9e). Apart from the most oblate case $h / r = 0.4$ , both the vertical aspect ratio $τ_{m}$ and the tilt angle $η_{m}$ show increasing asymmetry as $β$ decreases, with the anticyclonic values tending to increase while the cyclonic values decrease. Generally, it appears that, as we move away from $β = 1$ for both like-signed and opposite-signed interactions, the asymmetry increases.

3.3. Approximate EquilibriumSolutions

In order to clarify how the vortex characteristics depend on the background flow, it is worth trying to find approximate solutions for the equilibria. Given the form of the background flow, we have found that the equilibria are always tilted about the x-axis, implying that the axis vectors are:

$\hat{a} = (\begin{matrix} 1 \\ 0 \\ 0 \end{matrix}), \hat{b} = (\begin{matrix} 0 \\ cos η \\ - sin η \end{matrix}), \hat{c} = (\begin{matrix} 0 \\ sin η \\ cos η \end{matrix}) .$

(25)

This means that the shape matrix always has a particular form, with $B_{12} = B_{31} = 0$ , and so the other components can be written in terms of the parameters $δ$ , $τ$ , and $η$ :

$\begin{matrix} B_{11} & = & \frac{δ}{τ^{2 / 3}}, \end{matrix}$

(26a)

$\begin{matrix} B_{22} & = & \frac{δ}{τ^{2 / 3}} [1 - (1 - δ τ^{2}) {sin}^{2} η], \end{matrix}$

(26b)

$\begin{matrix} B_{33} & = & \frac{δ}{τ^{2 / 3}} [δ τ^{2} + (1 - δ τ^{2}) {sin}^{2} η], \end{matrix}$

(26c)

$\begin{matrix} B_{23} & = & \frac{δ}{τ^{2 / 3}} (δ τ^{2} - 1) sin η cos η . \end{matrix}$

(26d)

Similarly the self-induced flow matrix has only certain non-zero components, which can be written as:

$\begin{matrix} S_{v 12} & = & q_{0} + q_{1} {sin}^{2} η + q_{2} {sin}^{4} η, \end{matrix}$

(27a)

$\begin{matrix} S_{v 13} & = & [r_{0} + r_{1} {sin}^{2} η] sin η cos η, \end{matrix}$

(27b)

$\begin{matrix} S_{v 21} & = & s_{0} + s_{1} {sin}^{2} η, \end{matrix}$

(27c)

$\begin{matrix} S_{v 31} & = & t_{0} sin η cos η, \end{matrix}$

(27d)

where the coefficients are:

$\begin{matrix} q_{0} & = & ξ_{b} (- 1 + 2 ξ_{c}) - Ω_{b c} (ξ_{b} + 2 ξ_{c}), \end{matrix}$

(28a)

$\begin{matrix} q_{1} & = & ξ_{b} - ξ_{c} - 6 ξ_{b} ξ_{c} + 3 (ξ_{b} Ω_{a b} + ξ_{c} Ω_{c a}) + Ω_{b c} (7 ξ_{b} + 9 ξ_{c}), \end{matrix}$

(28b)

$\begin{matrix} q_{2} & = & 8 [ξ_{b} ξ_{c} - Ω_{b c} (ξ_{b} + ξ_{c})] - 3 (ξ_{b} Ω_{a b} + ξ_{c} Ω_{c a}), \end{matrix}$

(28c)

$\begin{matrix} r_{0} & = & ξ_{b} - ξ_{c} - 6 ξ_{b} ξ_{c} + 3 ξ_{c} Ω_{c a} + Ω_{b c} (5 ξ_{b} + 7 ξ_{c}), \end{matrix}$

(28d)

$\begin{matrix} r_{1} & = & 12 [ξ_{b} ξ_{c} - Ω_{b c} (ξ_{b} + ξ_{c})] - 3 (ξ_{b} Ω_{a b} + ξ_{c} Ω_{c a}), \end{matrix}$

(28e)

$\begin{matrix} s_{0} & = & ξ_{a} (1 - 2 ξ_{c}) + Ω_{c a} (2 ξ_{c} + ξ_{a}), \end{matrix}$

(28f)

$\begin{matrix} s_{1} & = & 2 ξ_{a} (ξ_{c} - ξ_{b}) + Ω_{a b} (ξ_{a} + 2 ξ_{b}) - Ω_{c a} (2 ξ_{c} + ξ_{a}), \end{matrix}$

(28g)

$\begin{matrix} t_{0} & = & \frac{f}{N} [Ω_{c a} (ξ_{c} - ξ_{a}) + Ω_{a b} (ξ_{a} - ξ_{b}) + ξ_{a} (ξ_{b} - ξ_{c})] . \end{matrix}$

(28h)

As a result of the form of the shape and flow matrices, equilibrium Equation (23) reduces significantly to just two equations:

$\begin{matrix} S_{31} B_{11} + S_{13} B_{33} + S_{12} B_{23} & = & 0, \end{matrix}$

(29a)

$\begin{matrix} S_{21} B_{11} + S_{12} B_{22} + S_{13} B_{23} & = & 0 . \end{matrix}$

(29b)

Given these two equations, as well as the fact that the $B_{33}$ component is constant, we have three equations in the three unknown variables $δ$ , $τ$ , and $η$ . However, since the self-induced flow matrix $S_{v}$ depends on these variables through the elliptical integrals, these equations cannot be directly solved. However, we can obtain approximate solutions by expanding these variables in the strain rate, i.e.,:

$\begin{matrix} δ & = & δ_{0} + γ δ_{1} + γ^{2} δ_{2} + O (γ^{3}), \end{matrix}$

(30a)

$\begin{matrix} τ & = & τ_{0} + γ τ_{1} + γ^{2} τ_{2} + O (γ^{3}), \end{matrix}$

(30b)

$\begin{matrix} η & = & η_{0} + γ η_{1} + γ^{2} η_{2} + O (γ^{3}) . \end{matrix}$

(30c)

We assume that the unperturbed solution when no strain is applied has a circular cross section, $δ_{0} = 1$ , is upright, $η_{0} = 0$ , and that the initial vertical aspect ratio is equal to the height-to-width aspect ratio, i.e., $τ_{0} = h / r$ . We can Taylor expand the self-induced flow matrix components about this unperturbed configuration as follows:

$S_{v} (δ, τ, η) = S_{v} (δ_{0}, τ_{0}, η_{0}) + {\frac{\partial S_{v}}{\partial δ}|}_{δ_{0}} (δ - δ_{0}) + {\frac{\partial S_{v}}{\partial τ}|}_{τ_{0}} (τ - τ_{0}) + {\frac{\partial S_{v}}{\partial η}|}_{η_{0}} (η - η_{0}) + \dots,$

(31)

where derivatives with respect to the aspect ratios can be related to derivatives with respect to the axis lengths by:

$\begin{matrix} \frac{\partial S_{v}}{\partial δ} & = & \frac{\partial S_{v}}{\partial a^{2}} \frac{\partial a^{2}}{\partial δ} + \frac{\partial S_{v}}{\partial b^{2}} \frac{\partial b^{2}}{\partial δ} + \frac{\partial S_{v}}{\partial c^{2}} \frac{\partial c^{2}}{\partial δ} = \frac{1}{δ^{2} τ^{2 / 3}} (δ^{2} \frac{\partial S_{v}}{\partial a^{2}} - \frac{\partial S_{v}}{\partial b^{2}}), \end{matrix}$

(32a)

$\begin{matrix} \frac{\partial S_{v}}{\partial τ} & = & \frac{\partial S_{v}}{\partial a^{2}} \frac{\partial a^{2}}{\partial τ} + \frac{\partial S_{v}}{\partial b^{2}} \frac{\partial b^{2}}{\partial τ} + \frac{\partial S_{v}}{\partial c^{2}} \frac{\partial c^{2}}{\partial τ} = - \frac{2}{3 δ τ^{5 / 3}} (δ^{2} \frac{\partial S_{v}}{\partial a^{2}} + \frac{\partial S_{v}}{\partial b^{2}} - 2 δ τ^{2} \frac{\partial S_{v}}{\partial c^{2}}), \end{matrix}$

(32b)

and where we use formulas for the derivatives of the elliptical integrals that have been derived in previous works [17,24]. By applying these formulas in the equilibrium relations in Equation (29) and matching the coefficients at different orders, we get, at the first order:

$δ_{1} = - \frac{3 (1 + cos 2 θ)}{A_{1}}; τ_{1} = 0; η_{1} = - \frac{3 τ_{0}^{2} sin 2 θ}{A_{2}},$

(33)

and at the second order, we find:

$\begin{matrix} δ_{2} & = & \frac{A_{3} η_{1}^{2} + τ_{0}^{2} [8 (β - 1) - A_{1} δ_{1}] δ_{1}}{2 τ_{0}^{2} A_{1}}, \end{matrix}$

(34a)

$\begin{matrix} τ_{2} & = & \frac{3 (τ_{0}^{2} - 1) η_{1}^{2}}{4 τ_{0}}, \end{matrix}$

(34b)

$\begin{matrix} η_{2} & = & \frac{[A_{4} δ_{1} + 2 (τ_{0}^{2} - 1) (1 - β)] η_{1}}{A_{2}}, \end{matrix}$

(34c)

where

$\begin{matrix} A_{1} & = & 4 (s_{0} + \frac{\partial s_{0}}{\partial δ}), \end{matrix}$

(35a)

$\begin{matrix} A_{2} & = & 2 [s_{0} + t_{0} + τ_{0}^{2} (r_{0} - s_{0})], \end{matrix}$

(35b)

$\begin{matrix} A_{3} & = & 4 [(τ_{0}^{2} - 1) (s_{0} + t_{0}) - τ_{0}^{2} (q_{1} + s_{1})], \end{matrix}$

(35c)

$\begin{matrix} A_{4} & = & 2 [τ_{0}^{2} (2 s_{0} - \frac{\partial r_{0}}{\partial δ}) - t_{0} - \frac{\partial t_{0}}{\partial δ}] + \frac{(τ_{0}^{2} - 1) A_{1}}{2}, \end{matrix}$

(35d)

where all of these coefficients are evaluated at the unperturbed state with $δ = δ_{0}$ , $τ = τ_{0}$ , and $η = η_{0}$ . From these equations, we see that, when there is no vertical shear ( $θ = 0$ ), the $τ$ and $η$ disturbance terms vanish, with only the horizontal aspect ratio $δ$ varying. In fact, the vertical aspect ratio has no linear dependence on the strain-rate, i.e., $τ_{1} = 0$ . Additionally, there is no dependence on $β$ at the first order. The dependence of $δ_{1}$ on $cos 2 θ$ implies that the values with the largest magnitude will tend to occur for $θ = 0^{\circ}$ ; therefore, as the cyclonic and anticyclonic strain rates are differently signed, we would expect the strongest asymmetry to occur at this value, i.e., when there is no vertical shear. This is overall what we have seen when examining the dependence on $θ$ in Figure 6 and Figure 7 for like-signed and opposite-signed interactions, respectively. On the other hand, the dependence of $τ_{2}$ on $sin θ$ means that the largest magnitudes occur for $θ = 45^{\circ}$ . This explains to some extent why the strongest asymmetries in $τ$ seen in Figure 6 and Figure 7 occur for larger values of $θ$ , although it tends to occur at even higher values, $θ \approx 60^{\circ}$ , implying that higher-order terms play a role.

The dependence on $β$ first appears at the second order for $δ$ and $η$ , but not for $τ_{2}$ . The appearance of $β - 1$ in $δ_{2}$ and $η_{2}$ means there will be a change in sign for these terms passing from the like-signed ( $β \geq 1$ ) and opposite-signed cases ( $β < 1$ ). In order to see the dependence of $τ$ on $β$ , we have to go to the next order, which can be easily done using the condition of constant $B_{33}$ and gives:

$τ_{3} = \frac{3 [δ_{1} η_{1} + 2 (τ_{0}^{2} - 1) η_{2}] η_{1}}{4 τ_{0}} .$

(36)

while $τ_{3}$ has a relatively simple form, the third order of $δ$ and $η$ is significantly more complicated depending on second- and third-order derivatives of the elliptical integrals, and so we do not compute them here. However, the third-order coefficient for $τ$ reveals the dependence on $β$ through the appearance of $η_{2}$ in Equation (36). For $η_{2}$ , the $β - 1$ term is multiplied by $τ_{0}^{2} - 1$ , meaning that the sign of this term also depends on whether the vortex is oblate or prolate, and so both $η$ and $τ$ will have a strong dependence on $h / r$ , as we have seen in Figure 8 and Figure 9.

In Figure 10a,b we plot $δ$ , $τ$ , and $η$ as a function of $γ$ for the cyclonic and anticyclonic cases, respectively, where $h / r = 2.5$ , $β = 2$ , and $θ = 30^{\circ}$ .

We plot the perturbation solutions at the first order (dashed blue curve) and the second order (dashed red curve), as well as the full solution (solid black curve). For the vertical aspect ratio, we also plot the third-order solution (green curve). For low values of the strain, both the first- and second-order perturbation solutions remain close to the full solution. At larger strain rates, the variables depart from the linear solution, with the direction of this deformation being captured by the higher-order perturbation solutions. In particular, the large asymmetry in $δ$ , with the cyclonic value decreasing as the anticyclonic value increases, is initially captured by the second-order perturbation. The change in the anticyclonic trend is due to the fact that, for like-signed interactions, the anticyclonic strain is negative, and so there is a sign change between the first- and second-order terms, i.e., $δ^{a c} \approx δ_{0} - γ δ_{1} + γ^{2} δ_{2}$ . However, as the strain-rate approaches the marginal value, even the second-order and third-order corrections cannot capture the change in the shape and orientation, implying that even higher-order terms are needed to fully capture the deformation.

In Figure 11, we plot the perturbation values for the opposite-signed case of an oblate vortex, with $h / r = 0.8$ , $β = - 4$ , and $θ = 30^{\circ}$ .

This case is one in which we have previously seen (Figure 9b) that there is strong asymmetry in the tilt angle. This asymmetry in $η$ is fairly well captured up to high strain values. For $δ$ and $τ$ , the perturbation approximation captures the equilibria for low strain rates. However, as was seen for the like-signed case, as we approach the marginal strain rate, the error grows, particularly for $δ^{c y}$ , implying the need for higher-order corrections.

4. Discussion andConclusions

In this work, we have examined how the shape and orientation of a vortex subjected to a background shear flow change as the background flow parameters are changed. In particular, we looked at how differences between cyclonic and anticyclonic equilibria can be induced by changing the level of strain and shear. It was shown that, in the absence of vertical shear, horizontal strain had the greatest impact on the horizontal cross section of the ellipsoidal vortices. Horizontal strain produced the strongest asymmetry in the horizontal aspect ratio of cyclonic and anticyclonic vortices, with cyclonic vortices becoming more deformed. This has been seen in shallow water turbulence, where anticyclones tend to be more axisymmetric than cyclones [7,11,25]. Instead, vertical shear tends to tilt the vortex from an upright position and changes the vertical aspect ratio. Vertical shear affects anticyclonic vortices more strongly than cyclonic vortices, with the strongest asymmetry occurring for large values of vertical shear. A mechanism that can damp the tilt angle produced by vertical shear has been seen in QG columnar vortices [26,27], although it is not known if this produces cyclonic–anticyclonic asymmetry in such vortices at the next order to QG. However, more recently, a study of a vortex merger at the finite Rossby number [28] showed that anticyclonic vortices are more susceptible to vertical shear, with vertically offset anticyclones merging from further apart than the equivalent cyclonic cases.

We also found that the shape and orientation of the equilibria also depend on the parameter $β \equiv 1 + κ_{v} / κ_{b}$ , i.e., the ratio of the target and background vortex strengths, with the vortex shape being most deformed as $β \to 1$ . However, the asymmetry between the cyclonic and anticyclonic equilibria increases as $β$ moves away from $β = 1$ for both like-signed ( $β \geq 1$ ) and opposite-signed interactions ( $β < 1$ ). The importance of both the PV ratio and the vortex volume have recently been investigated in the case of counter-rotating vortex pairs at the finite Rossby number [29] and for vertically aligned hetons [30]. Both these studies point out that, as a result of stronger motions induced by the anticyclonic vortex relative to the cyclone, this leads to the anticyclone dominating the evolution of the two vortices. Here we have seen that, over the range of values of the parameter $β$ , the horizontal aspect ratio of the cyclonic vortex is more deformed than the anticyclone, while the vertical aspect ratio of the anticyclone is generally more deformed. This deformation of the vertical aspect ratio when there is vertical shear is compensated by the increased tilting of the anticyclone with respect to the cyclone, an effect that is seen in the two vortex merger experiments of [28,29] for co-rotating and counter rotating vortices, respectively.

In order to understand the dependence of the equilibria on the background flow parameters, we derived an approximate model based on an expansion in the strain rate $γ$ . While the formulas derived capture the equilibria shape and orientation up to intermediate values of the strain rate, they cannot fully capture the dependence at high strain values up to the marginal strain, where strong interactions like mergers between vortices would be expected to take place. However, the higher-order terms do provide some insights into how the asymmetry between cyclonic and anticyclonic cases arises and its dependence on the background flow parameters. Future work extending the QG + 1 ellipsoidal model to the multi-vortex case would help provide better insights into this asymmetry by first providing a better understanding of ageostrophic corrections in two vortex interactions.

Here we have considered an idealized case of an isolated vortex in the presence of a background shear flow. This case can provide insights into cyclone–anticyclone asymmetry when there is no forcing, such as subsurface eddies in the deep ocean. However, in the upper ocean boundary, forcing and Ekman effects can also drive cyclone–anticyclone asymmetry [9,31,32,33]. Further work is needed to understand how the approach considered here can be applied to the case where there is such forcing.

Funding

This research received no external funding.

Data Availability Statement

The data in this study are available on request from the corresponding author.

Conflicts of Interest

The author declares no conflicts of interest.

References

Armi, L.; Herbert, D.; Oakey, N.; Price, J.; Richardson, P.; Rossby, T.; Ruddick, B. Two years in the life of a Mediterranean salt lens. J. Phys. Oceanogr. 1989, 19, 354–370. [Google Scholar] [CrossRef]
Dong, C.; McWilliams, J.C.; Liu, Y.; Chen, D. Global heat and salt transports by eddy movement. Nat. Commun. 2014, 5, 3294. [Google Scholar] [CrossRef] [PubMed]
Dai, J.; Wang, H.; Zhang, W.; Yuzhu, A.; Ren, Z. Observed spatiotemporal variation of three-dimensional structure and heat/salt transport of anticyclonic mesoscale eddy in Northwest Pacific. J. Ocean. Limnol. 2020, 38, 1654–1675. [Google Scholar] [CrossRef]
Griffiths, R.W.; Linden, P.F. The stability of vortices in a rotating, stratified fluid. J. Fluid Mech. 1981, 105, 283–316. [Google Scholar] [CrossRef]
Bell, G.D.; Bosart, L.F.A. A 15-Year climatology of Northern Hemisphere 500 mb closed cyclone and anticyclone centers. Mon. Weather Rev. 1989, 117, 2142–2163. [Google Scholar] [CrossRef]
Hakim, G.J.; Canavan, A.K. Observed cyclone-anticyclone tropopause vortex asymmetries. J. Atmos. Sci. 2005, 62, 231–240. [Google Scholar] [CrossRef]
Polvani, L.M.; McWilliams, J.C.; Spall, M.A.; Ford, R. The coherent structures of shallow-water turbulence: Deformation-radius effects, cyclone/anticyclone asymmetry and gravity-wave generation. Chaos 1994, 4, 177–186. [Google Scholar] [CrossRef]
Pallàs-Sanz, E.; Viúdez, Á. Three-dimensional ageostrophic motion in mesoscale vortex dipoles. J. Phys. Oceanogr. 2007, 37, 84–105. [Google Scholar] [CrossRef]
Roullet, G.; Klein, P. Cyclone-anticyclone asymmetry in geophysical turbulence. Phys. Rev. Lett. 2010, 104, 218501. [Google Scholar] [CrossRef]
Moisy, F.; Morize, M.; Rabaud, M.; Sommeria, J. Decay laws, anisotropy and cyclone–anticyclone asymmetry in decaying rotating turbulence. J. Fluid Mech. 2011, 666, 5–35. [Google Scholar] [CrossRef]
Perret, G.; Dubos, T.; Stegner, A. How large-scale and cyclogeostrophic barotropic instabilities favor the formation of anticyclonic vortices in the ocean. J. Phys. Oceanog. 2011, 41, 303–328. [Google Scholar] [CrossRef]
Naso, A. Cyclone-anticyclone asymmetry and alignment statistics in hom*ogeneous rotating turbulence. Phys. Fluids 2015, 27, 035108. [Google Scholar] [CrossRef]
Boffetta, G.; Toselli, F.; Manfrin, M.; Musacchio, S. Cyclone–anticyclone asymmetry in rotating thin fluid layers. J. Turbul. 2021, 22, 242–253. [Google Scholar] [CrossRef]
Stegner, A.; Le Vu, B.; Dumas, F.; Ali Ghannami, M.; Nicolle, A.; Durand, C.; Faugere, Y. Cyclone-anticyclone asymmetry of eddy detection on gridded altimetry product in the Mediterranean Sea. J. Geophys. Res. Ocean. 2021, 126, e2021JC017475. [Google Scholar] [CrossRef]
McKiver, W.J.; Dritschel, D.G. Balance in non-hydrostatic rotating stratified turbulence. J. Fluid Mech. 2008, 596, 201–219. [Google Scholar] [CrossRef]
McKiver, W.J. Balanced ellipsoidal vortex at finite Rossby number. Geophys. Astrophys. Fluid Dyn. 2020, 114, 453–480. [Google Scholar] [CrossRef]
McKiver, W.J. Balanced ellipsoidal vortex equilibria in a background shear flow at finite Rossby number. J. Fluid Mech. 2021, 926, A38. [Google Scholar] [CrossRef]
Dritschel, D.G.; Viúdez, Á. A balanced approach to modelling rotating stably stratified geophysical flows. J. Fluid Mech. 2003, 488, 123–150. [Google Scholar] [CrossRef]
McKiver, W.J.; Dritschel, D.G. Balanced solutions for an ellipsoidal vortex in a rotating stratified flow. J. Fluid Mech. 2016, 114, 453–480. [Google Scholar] [CrossRef]
McKiver, W.J.; Dritschel, D.G. The motion of a fluid ellipsoid in a general linear background flow. J. Fluid Mech. 2003, 474, 147–173. [Google Scholar] [CrossRef]
Meacham, S.P. Quasigeostrophic, ellipsoidal vortices in a stratified fluid. Dyn. Atmos. Ocean. 1992, 16, 189–223. [Google Scholar] [CrossRef]
Carlson, B.C. On computing elliptic integrals and functions. J. Math. Phys. 1965, 44, 35–51. [Google Scholar] [CrossRef]
Reinaud, J.N.; Dritschel, D.G.; Koudella, C.R. The shape of vortices in quasi-geostrophic turbulence. J. Fluid Mech. 2003, 474, 175–191. [Google Scholar] [CrossRef]
McKiver, W.J.; Dritschel, D.G. The stability of a quasi-geostrophic ellipsoidal vortex in a background shear flow. J. Fluid Mech. 2006, 560, 1–17. [Google Scholar] [CrossRef]
Graves, L.P.; McWilliams, J.C.; Montgomery, M.T. Vortex evolution due to straining: A mechanism for dominance of strong, interior anticyclones. Geophys. Astrophys. Fluid Dyn. 2006, 100, 151–183. [Google Scholar] [CrossRef]
Schecter, D.A.; Montgomery, M.T.; Reasor, P.D. A theory for the vertical alignment of a quasigeostrophic vortex. J. Atmos. Sci. 2001, 59, 150–168. [Google Scholar] [CrossRef]
Reasor, P.D.; Montgomery, M.T. Three-dimensional alignment and corotation of weak, TC-like vortices via linear vortex Rossby waves. J. Atmos. Sci. 2001, 58, 2306–2330. [Google Scholar] [CrossRef]
Reinaud, J.N.; Dritschel, D.G. The merger of geophysical vortices at finite Rossby and Froude number. J. Fluid Mech. 2018, 848, 388–410. [Google Scholar] [CrossRef]
Reinaud, J.N. Finite Froude and Rossby numbers counter-rotating vortex pairs. J. Fluid Mech. 2023, 971, A41. [Google Scholar] [CrossRef]
Reinaud, J.N.; McKiver, W.J.; Dritschel, D.G. Dynamical asymmetries in finite Rossby number baroclinic vortex pairs. Phys. Fluids, 2024; submitted. [Google Scholar]
Zavala Sansón, L. The asymmetric Ekman decay of cyclonic and anticyclonic vortices. Eur. J. Mech. B-Fluids 2001, 20, 541–556. [Google Scholar] [CrossRef]
Benthuysen, J.A.; Thomas, L.N. Asymmetries in vertical vorticity and vertical velocity arising during nonlinear hom*ogeneous spindown. Phys. Fluids 2012, 24, 076601. [Google Scholar] [CrossRef]
Zavala Sansón, L.; García-Martínez, I.M.; Sheinbaum, J. Nonlinear surface Ekman effects on cyclonic and anticyclonic vortices. J. Fluid Mech. 2023, 971, A35. [Google Scholar] [CrossRef]

Figure 1. Schematic of the ellipsoidal vortex of strength $κ_{v}$ whose axes of length a, b, and c are directed along the unit vectors $\hat{a}$ , $\hat{b}$ , and $\hat{c}$ , respectively, and where the vertical axis c is tilted about the x-axis by an angle $η$ . The vortex is in the presence of background flow induced by a vortex of strength $κ_{b}$ at a distance R and at an angle $θ$ with respect to the vortex. Note that $\hat{a}$ is parallel to the x-axis, which is directed out of the page in this perspective.

Figure 2. The vortex equilibrium values of $δ$ (left column), $τ$ (middle column), and $η$ (right column) as a function of the strain rate for QG (black line), QG + 1 cyclonic (blue), and QG + 1 anticyclonic and for (a) $θ = 0^{\circ}$ , (b) $θ = 30^{\circ}$ , and (c) $θ = 60^{\circ}$ for an oblate vortex with $h / r = 0.4$ and for like-signed interactions with $β = 2$ .

Figure 3. The vortex equilibrium values of $δ$ (left column), $τ$ (middle column), and $η$ (right column) as a function of the strain rate for QG (black line), QG + 1 cyclonic (blue), and QG + 1 anticyclonic and for (a) $θ = 0^{\circ}$ , (b) $θ = 30^{\circ}$ , and (c) $θ = 60^{\circ}$ for a prolate vortex with $h / r = 2.5$ and for like-signed interactions with $β = 2$ .

Figure 4. The vortex equilibrium values of $δ$ (left column), $τ$ (middle column), and $η$ (right column) as a function of the strain rate for QG (black line), QG + 1 cyclonic (blue), and QG + 1 anticyclonic and for (a) $θ = 0^{\circ}$ , (b) $θ = 30^{\circ}$ , and (c) $θ = 60^{\circ}$ for an oblate vortex with $h / r = 0.4$ and for opposite-signed interactions with $β = 0$ .

Figure 5. The vortex equilibrium values of $δ$ (left column), $τ$ (middle column), and $η$ (right column) as a function of the strain rate for QG (black line), QG + 1 cyclonic (blue), and QG + 1 anticyclonic and for (a) $θ = 0^{\circ}$ , (b) $θ = 30^{\circ}$ , and (c) $θ = 60^{\circ}$ for a prolate vortex with $h / r = 2.5$ and for opposite-signed interactions with $β = 0$ .

Figure 6. Plot of the marginal values of the strain rate $γ_{m}$ (left column), the horizontal aspect ratio $δ_{m}$ (second column), the vertical aspect ratio $τ_{m}$ (third column), and the tilt angle $η_{m}$ (right column) as a function of $θ$ for like-signed interactions with $β = 2$ and for (a) $h / r = 0.4$ , (b) $h / r = 0.8$ , (c) $h / r = 1$ , (d) $h / r = 1.25$ , and (e) $h / r = 2.5$ . The black, blue, and red curves are the QG, QG + 1 cyclonic, and QG + 1 anticyclonic cases, respectively.

Figure 7. Plot of the marginal values of the strain rate $γ_{m}$ (left column), the horizontal aspect ratio $δ_{m}$ (second column), the vertical aspect ratio $τ_{m}$ (third column), and the tilt angle $η_{m}$ (right column) as a function of $θ$ for opposite-signed interactions with $β = 0$ and for (a) $h / r = 0.4$ , (b) $h / r = 0.8$ , (c) $h / r = 1$ , (d) $h / r = 1.25$ , and (e) $h / r = 2.5$ . The black, blue, and red curves are the QG, QG + 1 cyclonic, and QG + 1 anticyclonic cases, respectively.

Figure 8. Plot of the marginal values of the strain rate $γ_{m}$ (left column), the horizontal aspect ratio $δ_{m}$ (second column), the vertical aspect ratio $τ_{m}$ (third column), and the tilt angle $η_{m}$ (right column) as a function of $β$ for like-signed interactions with $θ = 30^{\circ}$ and for (a) $h / r = 0.4$ , (b) $h / r = 0.8$ , (c) $h / r = 1$ , (d) $h / r = 1.25$ , and (e) $h / r = 2.5$ . The black, blue, and red curves are the QG, QG + 1 cyclonic, and QG + 1 anticyclonic cases, respectively.

Figure 9. Plot of the marginal values of the strain rate $γ_{m}$ (left column), the horizontal aspect ratio $δ_{m}$ (second column), the vertical aspect ratio $τ_{m}$ (third column), and the tilt angle $η_{m}$ (right column) as a function of $β$ for opposite-signed interactions with $θ = 30^{\circ}$ and for (a) $h / r = 0.4$ , (b) $h / r = 0.8$ , (c) $h / r = 1$ , (d) $h / r = 1.25$ , and (e) $h / r = 2.5$ . The black, blue, and red curves are the QG, QG + 1 cyclonic, and QG + 1 anticyclonic cases, respectively.

Figure 10. The vortex equilibrium quantities $δ$ (left column), $τ$ (middle column), and $η$ (right column) as a function of the strain rate for (a) QG + 1 cyclonic, and (b) QG + 1 anticyclonic cases where $h / r = 2.5$ , $β = 2$ , and $θ = 30^{\circ}$ . The solid black curve is the full solution, while the dashed blue, red, and green curves are for the approximate solutions up to the first ( $n = 1$ ), second ( $n = 2$ ), and third orders ( $n = 3$ for $τ$ only) in the strain rate, respectively.

Figure 11. The vortex equilibrium quantities $δ$ (left column), $τ$ (middle column), and $η$ (right column) as a function of the strain rate for (a) QG + 1 cyclonic, and (b) QG + 1 anticyclonic cases where $h / r = 0.8$ , $β = - 4$ , and $θ = 30^{\circ}$ . The solid black curve is the full solution, while the dashed blue, red, and green curves are for the approximate solutions up to the first ( $n = 1$ ), second ( $n = 2$ ), and third orders ( $n = 3$ for $τ$ only) in the strain rate, respectively.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).