## 1. Introduction

*D*/

*Dt*is the material derivative,

*α*the specific volume,

**= 2**

*ζ***Ω**+

**×**

**∇****u**the absolute vorticity vector,

**u**the three-dimensional velocity vector,

*θ*the potential temperature,

*θ̇*

**F**the frictional force per unit mass (for a succinct review of PV ideas and their relevance to tropical cyclone dynamics, see (McIntyre 1993). In a hurricane, the potential vorticity

*α*

**·**

*ζ*

**∇***θ*is not materially conserved because of the diabatic and frictional terms on the right-hand side of (1.1). The diabatic term is particularly important and can be written as

*α*|

**|**

*ζ***k**·

**∇***θ̇,*

**k**=

**/|**

*ζ***| is a unit vector pointing along the vorticity vector and**

*ζ***k**·

**∇***θ̇*

*θ̇*

*α*

**·**

*ζ*

**∇***θ̇*

*α*|

**| and**

*ζ***k**·

**∇***θ̇*

*α*

**·**

*ζ*

**∇***θ̇*

The two most complete observational studies of polygonal eyewalls are those of Lewis and Hawkins (1982) and Muramatsu (1986). Lewis and Hawkins examined polygonal eyewalls using time-lapse plan position indicator (PPI) film recordings of storms observed by land-based radars (Hurricanes David 1979; Anita 1977;Caroline 1975; Betsy 1965) and by airborne radars (Hurricanes Debbie 1969; Anita 1977). Their photographic analyses revealed that hurricane eyes can be in the shapes of hexagons, pentagons, squares, and triangles. Although the polygons were frequently incomplete, the straight line configurations were easily discernible. Lewis and Hawkins also noted that, while circular and elliptical eyes were also observed in the film recordings, they were seldom present for very long before straight lines and angles would appear. In addition, by showing that individual patterns could be simultaneously observed with two different radar systems, they disproved the possibility that polygonal shapes are artifacts of a particular radar system. Finally, Lewis and Hawkins also argued that polygonal eyes were not caused by topographic features near land-based radars, since polygonal features were also observed from airborne radars hundreds of kilometers from land.

At present, the most detailed radar observations of polygonal eyes have been obtained from Typhoon Wynne 1980, when it passed Japan’s Miyakojima radar (Muramatsu 1986). Observations from this small, remote island consisted of 15 h of PPI images showing polygonal eye features. During this period Wynne’s central surface pressure was between 920 and 935 mb. The polygonal features varied between square, pentagonal, and hexagonal, all of which rotated counterclockwise around the center. The pentagonal (which occurred most frequently) and hexagonal shapes had rotational periods of approximately 42 min, while the square shape had a rotational period of approximately 48 min. In addition, by presenting other radar and satellite images, Muramatsu offered further evidence that polygonal eyewalls are indeed real phenomena. Combining his observations with those of Lewis and Hawkins, Muramatsu defined some common characteristics of hurricanes exhibiting polygonal eyewalls. His results indicated that polygonal eyewalls were present in well-developed hurricanes having concentric eyewalls. Central pressures ranged from 920 to 950 mb, square to hexagonal eyewalls were most frequent, and the polygonal patterns tended to persist for tens of minutes. Contrary to Lewis and Hawkins, Muramatsu found no evidence of triangular patterns.

Lewis and Hawkins’ explanation for the existence of polygonal eyewalls was based largely on the internal gravity wave theories for spiral bands presented by Willoughby (1978) and Kurihara (1976). The basic idea is that gravity wave interference patterns due to the superposition of differing wavenumbers and periods would tend to produce polygonal bands. However, the existence of spiral bands can be explained without the use of transient gravity waves, but rather with PV dynamics and vortex Rossby waves (Guinn and Schubert 1993; Montgomery and Kallenbach 1997). According to these PV arguments, the inner spiral bands of hurricanes can form during the merger or axisymmetrization of asymmetric PV anomalies initiated internally by moist convection or externally through environmental forcing. It should also be mentioned that sharp bends in spiral bands, such as observed by Lewis and Hawkins for the case of Hurricane Caroline (1975), were observed in the numerical integrations of Guinn and Schubert (1993) whenever an existing vortex merged with another region of relatively high PV. In the present paper we shall explore the possibility that polygonal eyewalls can also be explained using PV dynamics.

Although Muramatsu (1986) did not develop a formal dynamical theory for the formation of polygonal eyes, he did suggest the phenomenon was related to an instability in the large shears of tangential wind near the inner edge of the eyewall (barotropic instability). Of special significance was the connection Muramatsu made between the formation of polygonal eyewalls and the formation of tornado suction vortices (Fujita et al. 1972). Although the two phenomena involve different scales, Muramatsu conjectured that both are related to barotropic instability. The barotropic instability of hurricane-scale annular PV rings was first explored by Guinn (1992) using a shallow water model.

With the exception of the recent work by Guinn (1992), the study of polygonal eyewalls has been given little attention since their discovery, perhaps for two reasons. First, some researchers simply do not believe polygonal features actually exist; they believe such “apparent features” are due to peculiarities in particular radar observing systems. Second, though they accept their existence, many researchers view polygonal eyewalls as a curiosity having little bearing on fundamental scientific questions such as hurricane motion, structure, and intensity change. Here we shall argue that polygonal eyes are more than just a curiosity and are of interest because they are symptoms of finescale PV mixing processes that are occurring continuously near the center of a hurricane.

Questions regarding the development of large horizontal shears, the possible barotropic instability of hurricane flows, and the formation of eyewall mesovortices are also of critical importance for the safe operation of reconnaissance aircraft and for the prediction of surface wind damage in landfalling storms (Wakimoto and Black 1994; Willoughby and Black 1996). Black and Marks (1991) have presented NOAA/P3 aircraft data from the eyewall region of Hurricane Hugo (1989) showing horizontal wind shears of 60 m s^{−1} over less than 1 km. They have also shown that such large horizontal shears can be associated with mesovortices having strong signatures in the vertical motion and radar reflectivity fields. More recently, Hasler et al. (1997) have used 1-min *GOES-9* images and concurrent NOAA/P3 radar observations to confirm the existence of an eyewall mesovortex in Hurricane Luis (1995). The mesovortex in Luis was apparently associated with hurricane track oscillations that could be interpreted as a series of eight cycloidal loops.

The outline of the present paper is as follows. In section 2 we argue, through use of linear stability analysis, that polygonal eyewalls can form simply as a result of barotropic instability. Although somewhat new to tropical cyclones, this idea has been used previously to examine the related subject of tornado suction vortices (e.g., Snow 1978; Staley and Gall 1979, 1984; Rotunno 1982, 1984; Gall 1982, 1983; Steffens 1988; Lin 1992;Finley 1997). Direct numerical simulations with the fully nonlinear equations, interpreted in terms of vorticity dynamics, are presented in sections 3 and 4 in support of our hypothesis. The method adopted here involves numerical simulations with the simplest model—an unforced barotropic nondivergent model. We demonstrate that the evolution of the PV field bears a striking resemblance to both radar and satellite images of eyewall features. However, since our simple model is not capable of predicting clouds and the observational data in hurricanes is not yet of sufficient quality and quantity to produce fine grain (Rossby–Ertel) PV maps, the connection between the evolution of PV fields and the evolution of convective fields is indirect but suggestive. Analytical predictions of the ultimate end-state of the mixing process based on minimum enstrophy and maximum entropy arguments are presented in sections 5 and 6, respectively. Section 7 concludes with a brief discussion of the relevance of the ideas developed here to real hurricanes.

## 2. Linear stability analysis of an annular region of vorticity

In a mature hurricane the frictional stress varies quadratically with wind speed and consequently the vertical velocity at the top of the boundary layer tends to be maximized just inside the radius of maximum tangential wind (Eliassen and Lystad 1977). Neglecting asymmetric processes, the frictionally forced vertical velocity organizes the convection into an annular ring near the radius of maximum winds. The cyclonic shear zone on the inner edge of the eyewall convection can be envisaged as an annular ring of uniformly high PV, with large radial PV gradients on its edges. On the inner edge of the annular ring the PV increases with radius, while on the outer edge the PV decreases with radius. In terms of PV (Rossby) wave theory, a PV wave on the inner edge of the annular ring will propagate counterclockwise relative to the flow there, while a PV wave on the outer edge will propagate clockwise relative to the flow there. Thus, it is possible for these two counterpropagating (relative to the tangential flow in their vicinity) PV waves to have the same angular velocity relative to the earth, that is, to be phase locked. If the locked phase is favorable, each PV wave will make the other grow, and exponential instability will result. This barotropic idealization is one view of the origin of polygonal eyewalls. In reality baroclinic and diabatic effects must also be important (after all, the basic-state PV field results from diabatic effects), but here we shall isolate the barotropic processes in order to investigate the extent to which they explain polygonal eyewall features.

*υ*

*r*) is a given function of radius

*r.*Using cylindrical coordinates (

*r, ϕ*), assume that the small-amplitude perturbations of the streamfunction,

*ψ*′(

*r, ϕ, t*), are governed by the linearized barotropic nondivergent vorticity (Rossby wave) equation

*ω*

*r*) =

*υ*

*r*)/

*r*is the basic-state angular velocity,

*ζ*

*r*) =

*d*(

*r*

*υ*

*rdr*the basic-state relative vorticity, (

*u*′,

*υ*′) = (−∂

*ψ*′/

*r*∂

*ϕ,*∂

*ψ*′/∂

*r*) the perturbation radial and tangential components of velocity, and ∂(

*rυ*′)/

*r*∂

*r*− ∂

*u*′/

*r*∂

*ϕ*= ∇

^{2}

*ψ*′ the perturbation vorticity.

^{1}Searching for modal solutions of the form

*ψ*′(

*r, ϕ, t*) =

*ψ*(

*r*)

*e*

^{i(mϕ−νt)}

*m*is the tangential wavenumber and

*ν*the complex frequency, we obtain from (2.1) the radial structure equation

*ξ*

_{1},

*ξ*

_{2},

*r*

_{1}, and

*r*

_{2}are constants. The constants

*ξ*

_{1}and

*ξ*

_{2}are the vorticity jumps (as one moves inward) at

*r*

_{1}and

*r*

_{2}. The case in which we are most interested has

*ξ*

_{1}< 0 and

*ξ*

_{2}> 0, that is, a ring of elevated vorticity.

^{2}

*ψ*′ = 0 for

*r*≠

*r*

_{1}and

*r*≠

*r*

_{2}, or equivalently, (2.2) reduces to

*r*

_{1}and

*r*

_{2}can be constructed from different linear combinations of

*r*

^{m}and

*r*

^{−m}in each region. This approach results in six undetermined constants. Requiring boundedness of

*ψ*(

*r*)

*r*→ 0 and

*r*→ ∞, and requiring continuity of

*ψ*(

*r*)

*r*=

*r*

_{1}and

*r*=

*r*

_{2}, reduces the number of undetermined constants to two [see section 4b of Michalke and Timme (1967) for details]. An alternative and more physically revealing approach is to write the general solution of (2.5), valid in any of the three regions, as a linear combination of the basis functions

*B*

^{(m)}

_{1}

*r*) and

*B*

^{(m)}

_{2}

*r*), defined by

*ψ*(

*r*)

*ψ*(

*r*)

_{1}

*B*

^{(m)}

_{1}

*r*

_{2}

*B*

^{(m)}

_{2}

*r*

_{1}and Ψ

_{2}are complex constants. Since

*dB*

^{(m)}

_{1}

*dr*is discontinuous at

*r*=

*r*

_{1}, the solution associated with the constant Ψ

_{1}has vorticity anomalies concentrated at

*r*=

*r*

_{1}and the corresponding streamfunction decays away from

*r*=

*r*

_{1}. Similarly, since

*dB*

^{(m)}

_{2}

*dr*is discontinuous at

*r*=

*r*

_{2}, the solution associated with Ψ

_{2}has vorticity anomalies concentrated at

*r*=

*r*

_{2}and the corresponding streamfunction decays away from

*r*=

*r*

_{2}.

_{1}and Ψ

_{2}, let us now integrate (2.2) over the narrow radial intervals between

*r*

_{1}−

*ϵ*and

*r*

_{1}+

*ϵ*and between

*r*

_{2}−

*ϵ*and

*r*

_{2}+

*ϵ*to obtain the jump (pressure continuity) conditions

*ω*

_{1}=

*ω*

*r*

_{1}) and

*ω*

_{2}=

*ω*

*r*

_{2}), and where we have assumed

*ν*≠

*m*

*ω*

_{1}and

*ν*≠

*m*

*ω*

_{2}. Substituting the solution (2.7) into the jump conditions (2.8) yields the matrix eigenvalue problem

*ν*

_{1}=

*m*

*ω*

_{1}− ½

*ξ*

_{1}and

*ν*

_{2}=

*m*

*ω*

_{2}− ½

*ξ*

_{2}are the pure (noninteracting) discrete vortex Rossby wave frequencies at the inner and outer interfaces. One can verify from (2.10) that

*ν*must be real for

*m*= 1, 2. This implies the vorticity field will remain stable to these disturbance patterns. The remaining wavenumbers can, however, produce frequencies with nonzero imaginary parts.

^{2}

If the basic-state vorticity jump at the outer interface were removed, the lower-left matrix element in (2.9) would disappear and the vortex Rossby wave on the inner interface would propagate with angular velocity *ν*_{1}/*m.* Similarly, if the basic-state vorticity jump at the inner interface were removed, the upper-right matrix element in (2.9) would disappear and the vortex Rossby wave on the outer interface would propagate with angular velocity *ν*_{2}/*m.* Thus, the system (2.9) can be regarded as a concise mathematical description of the interaction of two discrete counterpropagating vortex Rossby waves. The upper-right matrix element in (2.9) gives the effect of the outer vorticity anomaly pattern on the behavior of the inner interface, while the lower-left matrix element in (2.9) gives the effect of the inner vorticity anomaly pattern on the behavior of the outer interface. Note that the effect of these interactions decays with increasing wavenumber and decreasing values of the ratio *r*_{1}/*r*_{2}. For *ξ*_{1} < 0 and *ξ*_{2} > 0, the inner PV wave propagates counterclockwise relative to the basic-state tangential flow at *r* = *r*_{1}, while the outer PV wave propagates clockwise relative to the basic-state tangential flow at *r* = *r*_{2}.

*r*⩽

*r*

_{2}as

*ζ*

_{av}=

*ξ*

_{1}

*δ*

^{2}+

*ξ*

_{2}, where

*δ*=

*r*

_{1}/

*r*

_{2}. Then, defining

*γ*= (

*ξ*

_{1}+

*ξ*

_{2})/

*ζ*

_{av}as the ratio of the inner-region vorticity to the average vorticity, we can express

*ξ*

_{1}and

*ξ*

_{2}in terms of

*ζ*

_{av},

*δ,γ*as

*ξ*

_{1}= −

*ζ*

_{av}(1 −

*γ*)/(1 −

*δ*

^{2}) and

*ξ*

_{2}=

*ζ*

_{av}(1 −

*γδ*

^{2})/(1 −

*δ*

^{2}). Using these last two relations, and noting that

*ω*

_{1}= ½

*γζ*

_{av}and

*ω*

_{2}= ½

*ζ*

_{av}we can rewrite (2.10) as

*ν*/

*ζ*

_{av}as a function of the disturbance tangential wavenumber

*m*and the two basic-state flow parameters

*δ*and

*γ.*The imaginary part of

*ν*/

*ζ*

_{av}, denoted by

*ν*

_{i}/

*ζ*

_{av}, is a dimensionless measure of the growth rate. Isolines of

*ν*

_{i}/

*ζ*

_{av}as a function of

*δ*and

*γ*for

*m*= 3, 4, . . . , 8 are shown in Fig. 1. Note that all basic states with

*γ*< 1 satisfy the Rayleigh necessary condition for instability but that most of the region

*γ*< 1,

*δ*< ½ is in fact stable. Clearly, thinner annular regions (larger values of

*r*

_{1}/

*r*

_{2}) should produce the highest growth rates but at much higher tangential wavenumbers. Note also the overlap in the unstable regions of the

*γ–δ*plane for different tangential wavenumbers. For example, the lower right area of the

*γ–δ*plane is unstable to all the tangential wavenumbers

*m*= 3, 4, . . . , 8. We can collapse the six panels in Fig. 1 into a single diagram if, for each point in the

*γ–δ*plane, we choose the largest growth rate of the six wavenumbers

*m*= 3, 4, . . . , 8. This results in Fig. 2, which shows clearly the preference for higher wavenumbers as the annular ring becomes thinner.

Excluding values of *δ* larger than 0.9, the lower-right quarter of Fig. 2 contains isolines of the dimensionless growth rate *ν*_{i}/*ζ*_{av} between approximately 0.1 and 1.0, so that the dimensional growth rate *ν*_{i} satisfies 0.1*ζ*_{av} ⩽ *ν*_{i} ⩽ *ζ*_{av}. The average vorticity inside *r*_{2} can be expressed as the circulation around *r*_{2} (2*πr*_{2}*υ*_{2}) divided by the area inside *r*_{2} (*πr*^{2}_{2}*ζ*_{av} = 2*υ*_{2}/*r*_{2}. For a strong hurricane we can choose *υ*_{2} = 50 m s^{−1} at *r*_{2} = 50 km, in which case *ζ*_{av} = 2 × 10^{−3} s^{−1}. Then we have 2.0 × 10^{−4} ⩽ *ν*_{i} ⩽ 2.0 × 10^{−3} s^{−1}, or in terms of the *e*-folding time *τ*_{e} = *ν*^{−1}_{i}*τ*_{e} ⩽ 83 min. Although consideration of tangential wavenumbers larger than 8 and values of *δ* larger than 0.9 will yield even faster growth rates, such rapid growth rates are probably not realistic since they require very thin vorticity rings. However, basic-state structures yielding instabilities with 8–80 min *e*-folding times may be relevant to understanding hurricane structure.

## 3. Pseudospectral barotropic nondivergent model

Given a sign reversal in the radial PV gradient at lower-tropospheric levels in the hurricane, one expects the stage to be set for a convectively modified version of combined barotropic–baroclinic instability and PV redistribution. Since the PV field in the hurricane is induced by both boundary layer and moist processes, we generally expect these same processes, along with barotropic and baroclinic instability effects, to play a role in the evolution of asymmetric disturbances developing out of this background state. Regardless of the importance of boundary layer and moist processes, it is instructive to understand the conservative dynamics before nonconservative processes are included. With this viewpoint in mind we examine first the conservative (weakly dissipative) nonlinear vorticity dynamics using an unforced barotropic nondivergent model. Future work will investigate the implications of these preliminary results in a divergent barotropic and three-dimensional setting with and without nonconservative processes.

Insight into the early time nonlinear evolution of an annular vortex has been obtained by Dritschel (1986) and Lin (1992) using the method of contour dynamics/contour surgery (Zabusky et al. 1979; Zabusky and Overman 1983; Dritschel 1988, 1989; Ritchie and Holland 1993). This method is specifically designed for piecewise-constant vorticity distributions such as the one used in the three region model of section 2. In short, the method predicts the position of the contours separating the regions of constant vorticity. An initially unstable annulus of uniform vorticity with small undulations on its inner and outer contours can distort and evolve into a pattern in which the vorticity becomes “pooled” into rotating elliptical regions connected to each other by filaments or strands of high vorticity fluid. As the filaments become more and more intricately stretched and folded, contour dynamics/contour surgery adds more nodes to accurately follow the elongating contours, but also surgically removes finescale features. Surgical removal of vorticity during the simulation of a vorticity mixing process generally implies nonconservation of net circulation around the mixing region. As one of our goals is to characterize the ultimate end-state of this process, we prefer a numerical method that preserves exactly the net circulation. We consequently employ a pseudospectral model with ordinary (∇^{2}) diffusion.

*f*-plane barotropic nondivergent model are

*u*and

*υ*are the eastward and northward components of velocity,

*p*is the pressure, and

*ρ*is the constant density. Expressing the velocity components in terms of the streamfunction by

*u*= −∂

*ψ*/∂

*y*and

*υ*= ∂

*ψ*/∂

*x,*we can write the vorticity equation, derived from (3.1) to (3.3), as

*ζ*= ∇

^{2}

*ψ*is the relative vorticity and ∂( · , · )/∂(

*x, y*) is the Jacobian operator. Two quadratic integral properties associated with (3.4) on a closed or periodic domain are the energy and enstrophy relations

**∇***ψ*·

**∇***ψ dx dy*is the energy,

*ζ*

^{2}

*dx dy*is the enstrophy, and

**∇***ζ*·

**∇***ζ dx dy*is the palinstrophy. The diffusion term on the right-hand side of (3.4) controls the spectral blocking associated with the enstrophy cascade to higher wavenumbers.

*f*and

*ρ,*(3.7) can then be used to determine

*p*from

*ψ.*

Section 4 presents a representative numerical integration of (3.4) that demonstrates the formation of polygonal eyewalls, asymmetric eye contraction, and vorticity mixing. The solutions of (3.4) to be presented were obtained with a double Fourier pseudospectral code having 512 × 512 equally spaced collocation points on a doubly periodic domain of size 600 km × 600 km. The code was run with a dealiased calculation of the quadratic advection terms in (3.4). This results in 170 × 170 Fourier modes. Although the collocation points are only 1.17 km apart, a more realistic estimate of resolution is the wavelength of the highest Fourier mode, which is 3.53 km. While the gross features of the flow evolution presented below have been confirmed using lower spatial resolution and larger computational domains, the current configuration of numerical parameters yields a good compromise between computer speed and memory limitations and the desire for an adequately resolved inertial range. Other numerical details are as follows. Time differencing was accomplished with a standard fourth-order Runge–Kutta scheme using a 7.5-s time step. The chosen value of *ν* was 100 m^{2} s^{−1}, resulting in a 1/*e* damping time of 53 min for all modes having total wavenumber 170. This damping time lengthens to 3.5 h for modes having total wavenumber 85.

## 4. Redistribution of PV for initially hollow PV structures

### a. Initial condition

*r*

_{1},

*r*

_{2},

*d*

_{1},

*d*

_{2},

*ζ*

_{1}−

*ζ*

_{3},

*ζ*

_{2}−

*ζ*

_{3}, and

*ζ*

_{amp}are independently specified quantities, and the constant

*ζ*

_{3}is determined in order to make the domain average of

*ζ*(

*r, ϕ,*0) vanish. Here

*S*(

*s*) = 1 − 3

*s*

^{2}+ 2

*s*

^{3}is the basic cubic Hermite shape function satisfying

*S*(0) = 1,

*S*(1) = 0, and

*S*′(0) =

*S*′(1) = 0. Sensitivity tests using other broadbanded initial asymmetries yield qualitatively similar results provided the perturbation excites a nontrivial azimuthal wavenumber one component (either initially or through wave–wave interaction).

^{3}

For the representative numerical experiment we set the initial condition parameters to be *r*_{1} = 37.5 km, *d*_{1} = 7.5 km, *r*_{2} = 57.5 km, *d*_{2} = 7.5 km, *ζ*_{1} − *ζ*_{3} = 4.1825 × 10^{−4} s^{−1}, *ζ*_{2} − *ζ*_{3} = 3.3460 × 10^{−3} s^{−1}, and *ζ*_{amp} = 1.0 × 10^{−5} s^{−1} ≈ 0.003*ζ*_{2}. To make the domain average vorticity vanish we must choose *ζ*_{3} = −6.0653 × 10^{−5} s^{−1}. The symmetric part of the initial vorticity, tangential wind, and angular velocity fields are shown by the solid lines in Fig. 4 (the other lines in this figure will be discussed later). The vorticity of the ring is approximately nine times the vorticity of the central region. The initial tangential wind is weak inside 35 km but increases rapidly between 40 km and 50 km, where the vorticity is large. Note that the maximum tangential wind is approximately 54 m s^{−1} and lies on the outer edge of the vorticity ring near 60-km radius. The initial angular velocity *ω**π*/*ω**e*-folding time of 57.5 min. An approximate interpretation of this continuous vorticity distribution in terms of the three-region model yields a point in Fig. 2 near the transition between wavenumbers 4 and 3 along the line *γ* ≈ 0.20, which also predicts en *e*-folding time just under 60 min. A detailed analysis of the direct numerical integration (Fig. 3) at early times also yields a dominant mode-4 disturbance with a consistent *e*-folding time of approximately 60 min.

### b. Redistribution into a monopole

The results of the experiment are shown in the form of vorticity maps in Fig. 3 (every 2 h until *t* = 24 h and every 6 h thereafter). Near the beginning of the experiment the vorticity wave on the inner edge of the ring (30 ⩽ *r* ⩽ 45 km) is embedded in weak cyclonic flow of approximately 9 m s^{−1} (*ω*^{−3} s^{−1} or 7 h circuit time) and is propagating cyclonically relative to this flow. The vorticity wave on the outer edge of the ring (50 ⩽ *r* ⩽ 65 km) is embedded in strong cyclonic flow of approximately 54 m s^{−1} (*ω*^{−3} s^{−1} or 1.9 h turnaround time) and is propagating anticyclonically relative to this flow. As the two vorticity waves phase lock and help each other grow, the wave on the inner edge of the high vorticity ring is the first to reach the wave-breaking stage. During the breaking stage, interior particles with low and intermediate vorticity are drawn out into the high vorticity ring. This low and intermediate vorticity fluid spirals cyclonically into and carves up the high vorticity fluid in the original annular ring. This stretching and folding process turns the original high vorticity ring into an area of active enstrophy cascade. Two consequences of the withdrawal of intermediate and low PV fluid are a dramatic tightening of the PV gradient on the inside edge of the vorticity ring and a significant decrease in the area of the inner, low vorticity region. For example, the area of the inner, square-shaped, dark blue region at *t* = 8 h is approximately 60% of the dark blue inner area at *t* = 0 h. In other words, 40% of the low vorticity inner fluid has been mixed into the high vorticity ring. Thus, the annular ring shrinks inward and, through mixing, becomes broader and less intense. At *t* = 6 h the vorticity wave on the outer edge of the ring has increased enough in amplitude so that its four tips are in regions where the differential rotation is anticyclonic. In a period of 2 h, spiral bands form and, as vorticity is stripped off the outer edge of the high vorticity ring, the vorticity gradient is significantly tightened there also. Between 10 and 12 h the central region of low vorticity (dark blue) is knocked off-center and a long filament of it is withdrawn just to the south and east of the center of the domain. The remaining patch of low vorticity now begins to circle the vortex approximately every 2 h. For example, at 24 h it lies west-southwest of the center between *r* = 35 and 55 km. At these radii, the tangential wind is approximately 40 m s^{−1}, the orbital time is approximately 2 h, and there is weak anticyclonic differential rotation (see bottom two panels of Fig. 4).

In terms of the integral quantities enstrophy and palinstrophy, the period between 6 and 16 h is remarkable. During this period there is a large pulse of palinstrophy (5 times its initial value) and an associated [see Eq. (3.6)] rapid decay of enstrophy (as shown in Fig. 8). During these 10 h, approximately 32% of the original enstrophy is sent to small scales, where it is dissipated. In contrast, the angular momentum and kinetic energy are nearly conserved. For example, at *t* = 24 h, the angular momentum, kinetic energy, and enstropy are 99.9%, 99.2%, and 59.5% of their respective initial values.

Between 24 and 48 h the patch of low vorticity circles the vortex 10 times while the patch of high vorticity (located 30 km east-southeast of the domain center at *t* = 24 h) settles into the vortex center, accompanied by its associated trailing spiral bands of vorticity. At 48 h the nearly circular vorticity contour at the vortex center is 0.0025 s^{−1} and the central vorticity maximum is 0.0026 s^{−1}, which is 79% of the initial vorticity value in the annular ring. Thus, the fluid in the center of the vortex at 48 h should not be considered as a large patch that has been simply moved from the original vortex ring to the center under the constraint of material conservation of vorticity. Rather, the fluid in the vortex center at 48 h can be considered a mixture of some fluid which was initially in the high vorticity annular ring, some fluid which was inside this ring, and some fluid which was outside this ring. This interpretation will be further discussed in section 6, where the fractional contributions of fluid from each of the three regions will be determined via the maximum entropy hypothesis.

The symmetric (azimuthal mean) parts of the vorticity and tangential wind fields at selected times during the evolution are shown in Fig. 4. For this calculation, we take the vortex center as being the center of the domain.^{4} Early on (*t* = 4 h) a slight inward flux of mean vorticity on the interior side of the vorticity ring and a slight outward flux of vorticity on the exterior side of the ring is evident. Then there is a rather sudden transition to more extensive lateral mixing. Eventually, the lateral mixing process terminates and the mean vorticity profile assumes an almost monopolar structure except for a weak sign reversal in the radial vorticity gradient near *r* = 30 km. The mean tangential winds undergo a significant structure change as well. Corresponding to the vorticity transport to the center of the vortex, the mean tangential winds in the “eye” region increase substantially, while the overall tangential wind maximum decreases.

As one means of assessing whether the model resolution is adequate, we show the kinetic energy and enstrophy spectra in Fig. 5 at selected times through the violent mixing phase (0 ⩽ *t* ⩽ 16 h). An extensive enstrophy cascade is evident around *t* = 8 h. Subsequently, the small-scale energy increases somewhat but diffusion ultimately prevents the energy from piling up near the wavenumber cutoff.

It is interesting to speculate as to the possible physical effect of such a mesoscale mixing process on the thermodynamic structure of the hurricane eye region and vortex intensity. We generally think of the hurricane eye as being so dry that mixing of air from the eye into the eyewall cloud would be detrimental to convection. However, in addition to being dry, the air in the eye is also warm, so warm in fact that the air in the eye can often have an equivalent potential temperature higher than the air in the eyewall cloud (e.g., as shown below 700 mb in Fig. 16 of Hawkins and Imbembo 1976). In such a circumstance the mixing of air from the eye into the eyewall cloud as idealized here may actually be beneficial to convection. Emanuel (1997) goes a step further by arguing that the mechanical spinup of the eye region by vorticity mixing is critical for a hurricane to achieve its maximum intensity.

### c. Particle trajectories

*X*(

*t*),

*Y*(

*t*)] of a single particle is calculated from

*Ẋ*= −

*ψ*

_{y}(

*X, Y*) and

*Ẏ*=

*ψ*

_{x}(

*X, Y*). A small perturbation of the trajectory (

*δX, δY*) evolves according to the linear system

*e*

^{λt}, we find that

*λ*= ±

*Q*

*Q*=

*ψ*

^{2}

_{xy}

*ψ*

_{xx}

*ψ*

_{yy}. In regions where

*Q*< 0, two neighboring particles do not separate exponentially in time. Defining the two components of strain by

*S*

_{1}=

*u*

_{x}−

*υ*

_{y}= −2

*ψ*

_{xy}and

*S*

_{2}=

*υ*

_{x}+

*u*

_{y}=

*ψ*

_{xx}−

*ψ*

_{yy}, we can interpret

*Q*as a measure of the relative magnitudes of strain and vorticity, that is,

*Q*= ¼

*S*

^{2}

_{1}

*S*

^{2}

_{2}

*ζ*

^{2}). The coherent monopolar vortices that emerge during the evolution of decaying two-dimensional turbulence (McWilliams 1984; Benzi et al. 1988) have central regions with strongly negative

*Q*surrounded by regions of weakly positive

*Q.*

*Q*field can be obtained by taking the

*x*and

*y*derivatives of (3.4). Neglecting diffusion one obtains without any approximation

The time evolution of the *Q* field for the numerical experiment is shown in Fig. 6. For the first 4 h the *Q* field remains negative everywhere within the radius of maximum wind, with a thin ring of positive values just outside the radius of maximum wind. As the wave number 4 disturbance develops, four regions of positive *Q* appear near *r* = 40 km, where the intermediate and low vorticity fluid is being mixed into the high vorticity ring. At *t* = 8 h, the value of *Q* in these four regions is approximately 1 × 10^{−6} s^{−2}, which corresponds to an *e*-folding time for particle separation, *Q*^{−1/2}, of approximately 17 min.

There is an active enstrophy cascade occurring within the ring of high vorticity. Low vorticity fluid from the interior region is drawn out into the ring in four places and subjected to continual stretching and folding. This asymmetric eye contraction mechanism is physically quite different from the symmetric eye contraction mechanism proposed by Shapiro and Willoughby (1982) and Willoughby et al. (1982). The relative importance of these symmetric and asymmetric mechanisms in actual hurricanes remains an open question.

### d. A flawed heuristic argument

*r*

_{1}<

*r*<

*r*

_{2}. The vorticity is zero inside

*r*

_{1}and outside

*r*

_{2}. Now imagine that all of the vorticity becomes redistributed, without any mixing, into a circular patch of uniform vorticity centered at

*r*= 0. Let

*r*

_{1},

*r*

_{2},

*υ*

_{max}denote the inner radius, outer radius, and maximum tangential wind at the initial time. Let

*R*

_{2}, V

_{max}denote the corresponding values at the final time once the redistribution is complete. It is then easy to show that upon preserving the net circulation

## 5. Analytical prediction of the equilibrated end state via the minimum enstrophy hypothesis

Two approaches have been proposed in the literature that allow one to predict the equilibrated end state without explicitly simulating the details of the time-dependent nonlinear evolution of the flow. In this section we consider the minimum enstrophy approach and in section 6 we consider the maximum entropy approach. These approaches rest on different assumptions and lead to somewhat different solutions; we discuss these differences at the end of section 6.

### a. The selective decay hypothesis

An improvement to the flawed argument of section 4d that is more consistent with the nonlinear evolution shown in Fig. 3 invokes the selective decay hypothesis (e.g., Matthaeus and Montgomery 1980), which says that the barotropic instability process leads to an active cascade of the spectral distribution of enstrophy to high wavenumbers, where dissipation operates to decrease the enstrophy. By contrast, the area-integrated kinetic energy and angular momentum tend to be rugged integrals that are essentially unchanged.

Figure 8 shows a time series of the area-integrated kinetic energy, angular momentum, enstrophy, and palinstrophy for the numerical experiment shown in Fig. 3. Noteworthy features are the approximate invariance of the kinetic energy and angular momentum (1.4% reduction in kinetic energy and 0.1% reduction in angular momentum) over the duration of the mixing process. The enstrophy on the other hand is reduced substantially (47.4% overall), beginning near *t* = 8 h. The increase in the palinstrophy is consistent with the enstrophy decrease as must be the case from (3.6). A theory predicting the end-state of this mixing process in a dissipative (though nearly inviscid) evolution should be consistent with these properties.

Let *a* denote the outer edge of the vorticity mixing region. Following Leith (1984) we first hypothesize that, out of the family of vortices that have the same integrated angular momentum inside *r* = *a* and the same tangential wind at each radius outside *r* = *a,* the vortex with minimum integrated enstrophy inside *r* = *a* is the one toward which the flow actually evolves. This vortex is derived in section 5b and is called MinEV-*r* = *b* and the same tangential wind at each radius outside *r* = *b,* the vortex with minimum integrated enstrophy inside *r* = *b* is the one toward which the flow actually evolves. This vortex is called MinEV-

### b. Minimum enstrophy vortex with constrained circulation and angular momentum (MinEV-M )

To begin the MinEV-*r* = *a* is given by 2*π* ^{a}_{0}*rυr dr.* Because this angular momentum integral may be unbounded as *a* → ∞, it is convenient to work with the angular momentum deficit with respect to the initial vortex. Thus, for a vortex with initial and final tangential wind profiles *υ*_{0}(*r*) and *υ*(*r*), let us define the angular momentum deficit with respect to the angular momentum of the initial vortex as 2*π* ^{a}_{0}*r*(*υ*_{0} − *υ*)*r dr* = *π* ^{a}_{0}*r*^{2}(*ζ* − *ζ*_{0})*r dr,* where *ζ* = *d*(*rυ*)/*r dr* and where the second form of the angular momentum deficit follows from an integration by parts. For the special case where the initial tangential wind *υ*_{0}(*r*) is the corresponding tangential flow associated with a point vortex having the same circulation as the final flow, the angular momentum deficit is proportional to the total moment, about the origin, of the force impulse required to generate the difference motion (*υ* − *υ*_{0}) from rest (Batchelor 1967, section 7c). In a similar fashion, noting that all the enstrophy deficit is contained within *r* = *a,* we can write the enstrophy deficit of the hypothesized final axisymmetric flow as *π* ^{a}_{0}*ζ*^{2}_{0}*ζ*^{2})*r dr.*

*a,*the tangential wind profile

*υ*(

*r*), and the associated vorticity profile

*ζ*(

*r*) in search of that vortex that has maximum enstrophy deficit (i.e., minimum enstrophy) for fixed circulation and angular momentum.

^{5}The constancy of circulation at

*r*=

*a*requires

*υ*(

*a*) =

*υ*

_{0}(

*a*). Since the mixing radius

*a*is unknown, its first variation is related to the first variation in

*υ*at that point by

*δυ*(

*a*) = [

*υ*′

_{0}(

*a*) −

*υ*′(

*a*)]

*δa*= [

*ζ*

_{0}(

*a*) −

*ζ*(

*a*)]

*δa*(e.g., Fox 1987). Using these results, introducing the Lagrange multiplier

*γ,*and recalling Leibniz’s rule, the variational problem then becomes

*δζ*=

*d*(

*rδυ*)/

*rdr.*For the independent variation

*δa,*we obtain the transversality condition

*ζ*

*a*

*ζ*

_{0}

*a*

*δυ,*we obtain the Euler–Lagrange equation

*ζ*(

*r*) =

*d*(

*rυ*)/

*rdr*=

*ζ*

_{0}(

*a*) − (

^{1}/

_{2})

*γ*(

*a*

^{2}−

*r*

^{2}) for 0 ⩽

*r*⩽

*a,*where the constant of integration has been chosen such that (5.2) is satisfied. One further integration of this last relation for

*d*(

*rυ*)/

*rdr*yields

*υ*(

*r*) = (

^{1}/

_{2})

*rζ*

_{0}(

*a*) − (

^{1}/

_{8})

*γa*

^{2}

*r*[2 − (

*r*/

*a*)

^{2}] for 0 ⩽

*r*⩽

*a.*The condition that the final and initial tangential winds are equal at

*r*=

*a,*that is,

*υ*(

*a*) =

*υ*

_{0}(

*a*), yields

*ζ*

_{0}(

*r*) = 2

*υ*

_{0}(

*r*)/

*r*is the initial average vorticity inside

*r.*Remembering that the flow is unchanged for

*r*≥

*a,*we can write the solution for the tangential wind as

*ζ*

_{0}(

*a*) = 0. In that case it is readily shown that the radius of maximum wind is

*r*= (

^{2}/

_{3})

^{1/2}

*a*≈ 0.8165

*a,*the maximum wind is

*υ*= (4

*υ*

_{0}(

*a*) ≈ 1.089

*υ*

_{0}(

*a*), the peak vorticity is twice the average vorticity inside

*r*=

*a*[i.e.,

*ζ*(0) = 4

*υ*

_{0}(

*a*)/

*a*], and the peak angular velocity is twice the angular velocity at

*r*=

*a*[i.e.,

*ω*(0) = 2

*ω*

_{0}(

*a*)].

*a*remains to be determined. If (5.5) is substituted into the angular momentum constraint

^{a}

_{0}

*rυ*

_{0}

*r dr*=

^{a}

_{0}

*rυr dr,*and the integral on the right-hand side is evaluated, we obtain

*υ*

_{0}(

*r*) and

*ζ*

_{0}(

*r*) are given by the initial condition, (5.8) determines

*a.*We can summarize the predictions of the MinEV-

*υ*

_{0}(

*r*) and associated vorticity

*ζ*

_{0}(

*r*), first determine

*a*from (5.8). If multiple roots exist, choose that root that maximizes the enstrophy deficit.

^{6}The final adjusted tangential wind profile

*υ*(

*r*), vorticity profile

*ζ*(

*r*), and angular velocity

*ω*(

*r*) are then given by (5.5), (5.6), and (5.7).

^{7}

Using the initial tangential wind field given by the dotted curve in the second panel of Fig. 9 (the same profile as used in the direct numerical simulation), the solution of (5.8) yields *a* ≈ 82.13 km, so that the outer edge of the nonzero vorticity region shifts outward approximately 17.13 km. When this value of *a* is inserted into (5.5) and (5.6) we obtain the *υ*(*r*) and *ζ*(*r*) profiles given by the dash–dotted curves in the top two panels of Fig. 9. The associated angular velocity profile *ω*(*r*) = *υ*(*r*)/*r* is shown in the bottom panel of Fig. 9. Note that MinEV-^{−1} decrease in tangential winds near 60-km radius and a 26 m s^{−1} increase in tangential winds near 35-km radius, all in such a way that the total angular momentum (within any disk of radius ≥ 82.13 km) is invariant.

To determine the final pressure field we integrate the gradient wind relation *ρ*(*f* + *υ*/*r*)*υ* = *dp*/*dr* inward from *r* = 300 km, assuming *ρ* = 1.13 kg m^{−3}, *f* = 5 × 10^{−5} s^{−1}, and *p* = 1000 mb at *r* = 300 km. A plot of the resulting *p*(*r*) is shown by the dash–dotted curve in the third panel of Fig. 9. The dotted curve in the third panel of Fig. 9 gives the radial profile of pressure for the initial condition. Note that, even though the maximum tangential wind in MinEV-^{−1} weaker than the maximum tangential wind in the initial vortex, the central pressure in MinEV-

*r*=

*a*is invariant, that is,

^{a}

_{0}

*ζr dr*=

^{a}

_{0}

*ζ*

_{0}

*r dr,*and since the angular momentum inside

*r*=

*a*is also invariant, that is,

^{a}

_{0}

*r*

^{2}

*ζr dr*=

^{a}

_{0}

*r*

^{2}

*ζ*

_{0}

*r dr,*we can define an invariant mean radius

*r*

*r*

*r*= 0. For MinEV-

*r*

*r*

*a*{(

^{1}/

_{3})[1 +

*aζ*

_{0}(

*a*)/(4

*υ*

_{0}(

*a*))]}

^{1/2}≈ 46.7 km. Thus, in the top panel of Fig. 9,

*r*

*ζ*profile and the final

*ζ*profile predicted by the MinEV-

*r*

### c. Minimum enstrophy vortex with constrained circulation and energy (MinEV-E )

To begin the MinEV-*r* = *b* is given by *π* ^{b}_{0}*υ*^{2}_{0}*υ*^{2})*r dr.* We assume that, in the region *r* ≥ *b,* the final axisymmetric flow *υ*(*r*) is equal to the initial axisymmetric flow *υ*_{0}(*r*). In a similar fashion, noting that all the enstrophy deficit is contained within *r* = *b,* we can write the enstrophy deficit of the hypothesized final axisymmetric flow as *π* ^{b}_{0}*ζ*^{2}_{0}*ζ*^{2})*r dr.*

*b,*the tangential wind profile

*υ*(

*r*), and the associated vorticity profile

*ζ*(

*r*) in search of that vortex that has maximum enstrophy deficit for fixed energy. Because the derivation parallels the derivation of the MinEV-

*J*

_{0}and

*J*

_{1}denote Bessel functions of the first kind of order zero and one, respectively. The unknowns

*μ*and

*b*are determined from the requirements that

*ζ*be continuous at

*r*=

*b,*that is,

*ζ*(

*b*) =

*ζ*

_{0}(

*b*), and that energy be conserved, that is,

^{b}

_{0}

*υ*

^{2}

_{0}

*r*)

*r dr*=

^{b}

_{0}

*υ*

^{2}(

*r*)

*r dr.*For the initial tangential wind profile given by the dotted curve in Fig. 9, we find that

*μ*≈ 0.02912 km

^{−1}and

*b*≈ 84.42 km, so that the outer edge of the nonzero vorticity region shifts outward approximately 19.42 km. The final

*υ*(

*r*) and

*ζ*(

*r*) profiles are given by the dashed curves in the top two panels of Fig. 9. The associated angular velocity profile

*ω*(

*r*) =

*υ*(

*r*)/

*r*is shown in the bottom panel of Fig. 9. Note that the tangential winds predicted by MinEV-

^{−1}of the tangential winds predicted by MinEV-

To determine the final pressure field we again integrate the gradient wind relation inward from *r* = 300 km, assuming *ρ* = 1.13 kg m^{−3}, *f* = 5 × 10^{−5} s^{−1}, and *p* = 1000 mb at *r* = 300 km. The resulting pressure *p*(*r*) is plotted as the dashed line in the third panel of Fig. 9. Note that the central pressure in MinEV-

Also shown in Fig. 9 are the tangential mean *ζ*(*r*), *υ*(*r*), *p*(*r*), *ω*(*r*) for the direct numerical integration at 48 h (solid curves). A comparison of the MinEV-

## 6. Analytical prediction of the equilibrated end-state via the maximum entropy principle

The direct numerical simulation shown in Fig. 3 illustrates how the mixing process can produce vorticity patterns of increasing intricacy. An adaptive numerical method, such as contour dynamics, requires an increasing amount of computer time to advance one time step as the vorticity field becomes more complex. On the other hand, the pseudospectral method used to produce Fig. 3 is not adaptive and requires a fixed amount of computer time to advance one time step, no matter how complex the vorticity field. In spectral methods the production of finer and finer scales in vorticity is arrested by the model resolution and by the diffusion (or hyperdiffusion) processes operating near the resolution limit. While it is tempting to run spectral models at higher and higher resolution in order to follow vorticity structures to finer and finer scales, such costly pursuits do not necessarily yield fundamental advances in our understanding of the dynamics. A statistical mechanics approach may be more useful. Using the standard point-vortex model, Persing and Montgomery (1995) simulated the evolution of perturbed vortex rings as a way of determining the final radial vorticity profile in statistical equilibrium. A theoretical approach applicable to continuous vorticity distributions and based on the maximum entropy principle has recently been developed by Miller (1990), Robert (1991), Robert and Sommeria (1991, 1992), Sommeria et al. (1991), Miller et al. (1992), Whitaker and Turkington (1994), Chavanis and Sommeria (1996), and Turkington and Whitaker (1996). Although the maximum entropy theory suffers from the weakness of predicting vortex merger when the 2D Euler equations predict corotation without merger (Whitaker and Turkington 1994), one nevertheless hopes that the theory makes reliable predictions when mixing processes are dominant. As a test of the theory and as a foundation for future work we present here a simple version of this maximum entropy argument for the unforced vortex mixing problem.

### a. Tertiary mixing case

The maximum entropy argument rests on two views of the vorticity field after mixing: a macroscopic view, which sees a smooth distribution of vorticity, and a microscopic view, which reveals the intricate details produced by stretching and folding of vorticity in nonlinear flow. Rather than attempting to describe the finescale distribution (“microstate”) directly, we represent the macroscopic vorticity at a point statistically by averaging over all possible microstates in a small neighborhood of the point.

*r*

_{1},

*r*

_{2},

*ζ*

_{1},

*ζ*

_{2}, and

*ζ*

_{3}are specified constants. With this initial condition the problem reduces to one of tertiary mixing. Suppose we sample the vorticity at

*N*points within a small neighborhood of

*r.*Let

*n*

_{1}denote the number of points at which the vorticity value

*ζ*

_{1}is found,

*n*

_{2}the number of points at which the vorticity value

*ζ*

_{2}is found, and

*N*−

*n*

_{1}−

*n*

_{2}the number of points at which the vorticity value

*ζ*

_{3}is found. Then

*ρ*

_{1}(

*r*) =

*n*

_{1}/

*N*denotes the probability, at point

*r,*of finding the vorticity

*ζ*

_{1},

*ρ*

_{2}(

*r*) =

*n*

_{2}/

*N*the probability of finding the vorticity

*ζ*

_{2}, and

*ρ*

_{3}(

*r*) = 1 −

*ρ*

_{1}(

*r*) −

*ρ*

_{2}(

*r*) the probability of finding vorticity

*ζ*

_{3}. The number of possible arrangements having

*n*

_{1}points with vorticity

*ζ*

_{1},

*n*

_{2}points with vorticity

*ζ*

_{2}, and

*N*−

*n*

_{1}−

*n*

_{2}points with vorticity

*ζ*

_{3}is the multiplicity function

*W,*which is given by

*W*= ln

*N*! − ln

*n*

_{1}! − ln

*n*

_{2}! − ln(

*N*−

*n*

_{1}−

*n*

_{2})!. Using the Stirling approximation (e.g., ln

*N*! ≈

*N*ln

*N*−

*N*for large

*N*), we obtain

*ρ*

_{1},

*ρ*

_{2}) plane for

*ρ*

_{1}≥ 0 and

*ρ*

_{2}≥ 0 with

*ρ*

_{1}+

*ρ*

_{2}⩽ 1 are shown in Fig. 10. Note that the entropy density −

*ρ*

_{1}ln

*ρ*

_{1}−

*ρ*

_{2}ln

*ρ*

_{2}− (1 −

*ρ*

_{1}−

*ρ*

_{2}) ln(1 −

*ρ*

_{1}−

*ρ*

_{2}) approaches zero as (

*ρ*

_{1},

*ρ*

_{2}) → (0, 0), (0, 1), (1, 0), and that its maximum value of ln3 occurs at

*ρ*

_{1}=

*ρ*

_{2}=

^{1}/

_{3}. In other words the multiplicity of microstates is a maximum when a third of the sampled points in the neighborhood of

*r*have vorticity

*ζ*

_{1}, a third have vorticity

*ζ*

_{2}, and a third have vorticity

*ζ*

_{3}.

*S*[

*ρ*

_{1}(

*r*),

*ρ*

_{2}(

*r*)] as

*S*[

*ρ*

_{1}(

*r*),

*ρ*

_{2}(

*r*)] measures the loss of information in going from the fine grain (microscopic) view to the coarse grain (macroscopic) view. The macroscopic vorticity is given in terms of

*ρ*

_{1}(

*r*) and

*ρ*

_{2}(

*r*) by

*ζ*(

*r*) =

*ζ*

_{1}

*ρ*

_{1}(

*r*) +

*ζ*

_{2}

*ρ*

_{2}(

*r*) +

*ζ*

_{3}[1 −

*ρ*

_{1}(

*r*) −

*ρ*

_{2}(

*r*)]. To find the most probable macroscopic state, we must find the particular

*ρ*

_{1}(

*r*) and

*ρ*

_{2}(

*r*), which maximize

*S*[

*ρ*

_{1}(

*r*),

*ρ*

_{2}(

*r*)] subject to all the integral constraints associated with the inviscid vorticity dynamics. In other words, the variational problem is to find the expectation functions

*ρ*

_{1}(

*r*) and

*ρ*

_{2}(

*r*) by maximizing (6.1a) subject to the circulation constraints

*ζ̃*

_{1}

*ζ*

_{1}−

*ζ*

_{3}and

*ζ̃*

_{2}

*ζ*

_{2}−

*ζ*

_{3}, and introducing the Lagrange multipliers

*α*

_{1},

*α*

_{2},

*β, γ,*the variational problem is

*υ*=

*dψ*/

*dr*and

*δζ*=

*d*(

*rδυ*)/

*rdr*=

*ζ̃*

_{1}

*δρ*

_{1}

*ζ̃*

_{2}

*δρ*

_{2}

*δρ*

_{1}and

*δρ*

_{2}, we obtain

*ρ*

_{1}(

*r*) and

*ρ*

_{2}(

*r*), we obtain

*ζ*(

*r*) =

*ζ*

_{1}

*ρ*

_{1}(

*r*) +

*ζ*

_{2}

*ρ*

_{2}(

*r*) +

*ζ*

_{3}[1 −

*ρ*

_{1}(

*r*) −

*ρ*

_{2}(

*r*)], we obtain

*ρ*

_{1}(

*r*) and

*ρ*

_{2}(

*r*) are given by (6.4a,b), Eq. (6.4c) is a nonlinear ordinary differential equation for

*ψ*(

*r*) with yet to be determined Lagrange multipliers

*α*

_{1},

*α*

_{2},

*β, γ.*The equations for

*α*

_{2},

*α*

_{2},

*β, γ*are obtained by enforcing the constraints (6.1b)–(6.1e). Thus, evaluating the integrals on the right-hand sides of (6.1b) and (6.1c), we obtain

^{∞}

_{0}

*ψ*(

*ζ*−

*ζ*

_{3})

*r dr*=

^{r2}

_{0}

*ψ*

_{0}(

*ζ*

_{0}−

*ζ*

_{3})

*r dr.*When the integral on the right-hand side of this form of the energy constraint is evaluated, we obtain

^{∞}

_{0}

*r*

^{2}(

*ζ*−

*ζ*

_{3})

*r dr*=

^{r2}

_{0}

*r*

^{2}(

*ζ*

_{0}−

*ζ*

_{3})

*r dr.*When the integral on the right-hand side of this form of the angular momentum constraint is evaluated, we obtain

*ρ*

_{1}(

*r*),

*ρ*

_{2}(

*r*),

*ψ*(

*r*),

*α*

_{1},

*α*

_{2},

*β, γ,*given the constants

*r*

_{1},

*r*

_{2},

*ζ*

_{1},

*ζ*

_{2},

*ζ*

_{3}.

Analytical solutions of the system (6.4) are not easily obtained, and numerical methods are generally required. To solve (6.4) we have used the iterative algorithm developed by Turkington and Whitaker (1996). For the constants defining the initial condition, we choose *r*_{1} = 37.5 km, *r*_{2} = 57.5 km, *ζ̃*_{1}^{−4} s^{−1}, *ζ̃*_{2}^{−3} s^{−1}, *ζ*_{3} = −6.0653 × 10^{−5} s^{−1}. This is the three region approximation to the initial condition used in the direct numerical simulation described in section 4. When these constants are used, the numerical solution of (6.4) yields the radial profiles of *ρ*_{1}(*r*), *ρ*_{2}(*r*), and *ρ*_{3}(*r*) = 1 − *ρ*_{1}(*r*) − *ρ*_{2}(*r*) shown in Fig. 11. These profiles can be interpreted as follows. In the MaxSV end-state, fluid particles from the initial annular ring have the highest probability (∼64%) of ending up in the central core, while fluid particles from the initial core have the highest probability (∼16%) of ending up just outside the initial annular ring near *r* = 67 km. In the central core of the MaxSV end-state there is a higher probability of finding air that was originally outside the annular ring (∼28%) than finding air that was originally inside the annular ring (∼9%). In this sense, the vortex has been “turned inside-out.” Such intense inside-out turning is the typical fate of highly unstable initial vortices with very low central vorticity. In contrast, for weakly unstable initial vortices whose central vorticity is only slightly lower than the vorticity of the annular ring, the amount of mixing is predicted to be much less, and typically *ρ*_{1} > *ρ*_{3} at *r* = 0.

The resulting vorticity profile, computed from *ζ*(*r*) = *ζ*_{1}*ρ*_{1}(*r*) + *ζ*_{2}*ρ*_{2}(*r*) + *ζ*_{3}[1 − *ρ*_{1}(*r*) − *ρ*_{2}(*r*)], is shown by the thin solid line in the top panel of Fig. 12. The MaxSV-predicted reduction (∼37%) in the maximum vorticity value from the initial condition to the final state is a consequence of the fundamental mixing processes occurring during the flow evolution. This emphasizes the fact that material conservation of vorticity on the macroscale is not a useful description of the flow evolution, even though material conservation of vorticity on the microscale is a useful description.

The integration of this vorticity profile yields the tangential wind profile shown by the thin solid line in the middle panel of Fig. 12. The corresponding pressure field, computed from the gradient balance relation as in sections 5b and 5c, is shown by the thin solid line in the bottom panel. For comparison, the MinEV-

### b. Maximum entropy vortex with constrained circulation and angular momentum only

The maximum entropy theory for tertiary mixing has four Lagrange multipliers (*α*_{1}, *α*_{2}, *β, γ*) and four constraints [the circulation constraints (6.4d), (6.4e), the energy constraint (6.4f), and the angular momentum constraint (6.4g)]. Can we obtain essentially equivalent solutions if we drop one of the constraints and, hence, one of the Lagrange multipliers? We have found that the answer to this question is, in many cases, no. To illustrate, consider the solution obtained when the energy constraint (6.4f) and the Lagrange multiplier *β* are dropped from the problem. Note that this significantly simplifies the problem (6.4) since the *βψ* terms disappear from (6.4a), (6.4b) and the determination of *α*_{1}, *α*_{2}, *γ* via (6.4d), (6.4e), (6.4g) decouples from the determination of *ψ* via (6.4c). Using the same values of *r*_{1}, *r*_{2}, *ζ̃*_{1}*ζ̃*_{2}*ζ*_{3} given in section 6a, we now obtain from this reduced MaxSV theory the radial profiles of *ζ*(*r*), *υ*(*r*), and *p*(*r*) shown by the dash–dotted curves in Fig. 12. As is evident from Fig. 12, this reduced version of MaxSV theory produces less mixing, and hence higher central vorticity, than the complete MaxSV theory. The differences are large enough to conclude that this particular reduced form of MaxSV theory is not a useful approximation to the complete MaxSV theory.

### c. Comparison of the MinEV and MaxSV theories

The MinEV and MaxSV approaches for predicting the final end-state of the barotropically unstable vorticity ring rest on different assumptions. For the initial basic state of Fig. 4, both produce solutions that agree well with the direct numerical solutions. However, for the tertiary mixing case, as the radial width of the initial annular vorticity patch gets larger, the MinEV-*r*_{1} = 0) where little or no adjustment is expected.^{8} The final vorticity profiles suggest that the MinEV-

In contrast, the maximum entropy adjustment process does not produce spuriously large vorticity, as can be seen as follows. Since *ζ*(*r*) = *ζ*_{1}*ρ*_{1}(*r*) + *ζ*_{2}*ρ*_{2}(*r*) + *ζ*_{3}[1 − *ρ*_{1}(*r*) − *ρ*_{2}(*r*)] with 0 ⩽ *ρ*_{1}(*r*) ⩽ 1, 0 ⩽ *ρ*_{2}(*r*) ⩽ 1, *ρ*_{1}(*r*) + *ρ*_{2}(*r*) ⩽ 1, the final vorticity *ζ*(*r*) is a mixture of *ζ*_{1}, *ζ*_{2}, *ζ*_{3}. In such a mixture it is impossible for the final vorticity to lie outside the range set by the maximum and minimum values of *ζ*_{1}, *ζ*_{2}, *ζ*_{3}. In this sense the predictions of the maximum entropy theory are more reliable than those of the minimum enstrophy theory.

While the MaxSV approach does not suffer from the same weaknesses of the MinEV approach for describing the redistribution of an initial vorticity ring, the MaxSV approach, when applied to atmospheric vortex merger problems, does not appear to be entirely defect-free either, since it predicts vortex merger when the 2D Euler equations predict corotation without merger (Whitaker and Turkington 1994).^{9}

## 7. Concluding remarks

Although considerable insight into the physics of tropical cyclones has been acquired using axisymmetric theory and models (e.g., Ooyama 1969), fundamental questions remain concerning the role of asymmetric processes in the cyclone life cycle. Questions associated with asymmetric potential vorticity redistribution in tropical cyclones can be studied with a hierarchy of dynamical models. These include the barotropic nondivergent model, the barotropic divergent model (shallow water equations), the quasigeostrophic model, the asymmetric balance model, the quasi-static primitive equation model, and the full nonhydrostatic primitive equation model. The latter two models may include parameterized or explicit moist physical processes. The present study complements recent work examining tropical cyclogenesis as a PV redistribution problem governed by vortex Rossby waves on PV monopoles (Montgomery and Enagonio 1998) and has focused on PV redistribution in or near the eyewall region of mature hurricanes possessing an elevated ring of PV near the eyewall. For simplicity, the study has been limited to the barotropic nondivergent model. We thus have bypassed all questions associated with vertical structure and moist physical processes. Despite these limitations, the results seem sufficiently interesting to warrant further study with more complete physical models.

In the numerical simulations presented here, a ring of elevated vorticity was perturbed with azimuthally broadbanded initial conditions as a means of demonstrating a simple dynamical mechanism for the formation of polygonal eyewalls, asymmetric eye contraction, and vorticity redistribution. The intent was not to replicate the evolution of a hurricane eyewall, but rather to isolate the fundamental dynamics believed responsible for the formation of polygonal eyewalls and mesovortices. The simulations indicate that the barotropic instability associated with annular regions of relatively high vorticity results in polygonal shapes before the ultimate rearrangement into a nearly monopolar circular vortex. Although the minimum enstrophy description tends to overestimate the degradation of the maximum vorticity, the numerically simulated end-state nevertheless suggests a significant reduction of the initial vorticity maximum. These results suggest a limitation of the idea that the vorticity within the “core” of vortices is a rugged invariant (Carnevale et al. 1992). Analytical predictions based on maximum entropy arguments support this view.

In the eyewall of a real hurricane, frictional convergence, and moist convection continually act to concentrate high vorticity there, thus satisfying the necessary condition for barotropic exponential instability and the sufficient condition for a mode-1 algebraic instablitiy (see footnote 2). Natural convective asymmetries near the eyewall provide the perturbations that allow these instabilities to grow. As the instabilities grow the vorticity pools into a small number of “pockets,” creating the appearance of a polygon on the inner edge of the original annular region. These pools or pockets of PV are also likely to be responsible for mesovortices. At later times in the numerical simulations most of the vorticity in these pockets gets strained by the mean shear and is ultimately diffused by viscosity at small scales. The remaining high vorticity gets advected into the central region, thereby spinning up the center or “eye.” Our numerical and analytical results therefore suggest an explicit quasi-two-dimensional mechanism for the spinup and maintenance of the hurricane’s eye circulation and thermal structure (via thermal wind balance), which is thought to be a crucial process for the attainment of maximum intensity (Emanuel 1997, and references).

There are a variety of ways the present work could be extended. Within the context of the barotropic nondivergent model, the theoretical predictions of sections 5 and 6 could be extended from the symmetric to the asymmetric case. Of particular interest would be the determination of the bifurcations between the symmetric and asymmetric structures predicted by the maximum entropy vortex arguments. Outside the context of the nondivergent model, both the direct numerical integrations and the theoretical predictions could be generalized to the shallow water equations. This generalization would allow a distinction between vorticity and potential vorticity and a distinction between enstrophy and potential enstrophy, and would allow the energy constraint to include both kinetic and potential energy. The generalization to the quasi-static primitive equations in isentropic coordinates would not be significantly more difficult than the shallow water case. Finally, there remains the question of the relevance of the concepts presented here to the moist convective environment in real hurricanes (e.g., Schade 1994; Guinn and Schubert 1994). This difficult question might partially be answered through the use of PV diagnostics on the output of “full-physics models.”

In closing we would like to point out that, in addition to the connections with tornado suction vortices and with the dynamic instability of the stratospheric polar vortex (e.g., Ishioka and Yoden 1994, 1995; Bowman and Chen 1994), there is a remarkable connection between the PV redistribution problem in hurricanes and the electron density redistribution problem recently studied in experimental plasma physics. The plasma experiment involves the creation of a confined pure electron plasma in a cylindrical trap, with a uniform axial magnetic field providing radial confinement, and bias voltages applied to the cylinder ends providing axial confinement. Once produced, the confined plasma evolves due to the drift flow, which is determined by the cross product of the electric and magnetic fields (i.e., **E** × **B** drift). The drift flow is two-dimensional in a plane perpendicular to the axis of the cylinder. Since the drift velocity changes the electron density according to an equation identical to the inviscid form of the vorticity equation (3.4), and since electron density is related to the Laplacian of the electrostatic potential, the two-dimensional equations for the evolution of the electron plasma are isomorphic with the two-dimensional Euler equations for inviscid incompressible flow (electron density ↔ vorticity and electrostatic potential ↔ streamfunction). Compared to laboratory experiments with ordinary fluids (e.g., Chomaz et al. 1988), the electron plasma experiments have the advantage of being able to simulate very high Reynolds number flow; this strengthens the analogy with geophysical fluid dynamics. For further details on these extraordinary experiments readers are referred to the recent papers by Peurrung and Fajans (1993), Huang and Driscoll (1994), and Fine et al. (1995).

## Acknowledgments

This work was supported in part by the National Science Foundation, under Grants ATM-9422525 and ATM-9529295, and by the Office of Naval Research under Grant ONR-N00014-93-1-0456. We also acknowledge the hospitality and support provided by the Isaac Newton Institute for Mathematical Sciences during the 1996 Program in The Mathematics of Atmosphere and Ocean Dynamics, where we benefited greatly from discussions with Chuck Leith, Rick Salmon, and Phil Morrison. We would also like to express our thanks to Paul Reasor and Mark Hidalgo for their assistance with the stability analysis of continuous profiles, and to Lloyd Shapiro and to an anonymous reviewer for their helpful comments.

## REFERENCES

Batchelor, G. K., 1967:

*An Introduction to Fluid Dynamics.*Cambridge University Press, 615 pp.Benzi, R., S. Patarnello, and P. Santangelo, 1988: Self-similar coherent structures in two-dimensional decaying turbulence.

*J. Phys. A: Math. Gen.,***21,**1221–1237.Black, P. G., and F. D. Marks, 1991: The structure of an eyewall meso-vortex in Hurricane Hugo (1989). Preprints,

*19*Miami, FL, Amer. Meteor. Soc., 579–582.^{th}Conf. on Hurricanes and Tropical Meteorology,Bowman, K. P., and P. Chen, 1994: Mixing by barotropic instability in a nonlinear model.

*J. Atmos. Sci.,***51,**3692–3705.Butler, K. M., 1991: Equilibration and optimal excitation in viscous shear flow. Ph.D. dissertation, Harvard University, 133 pp. [Available from Dept. of Earth and Planetary Sciences, Harvard University, Cambridge, MA 02138.].

Carnevale, G. F., J. C. McWilliams, Y. Pomeau, J. B. Weiss, and W. R. Young, 1992: Rates, pathways and end states of nonlinear evolution in decaying two-dimensional turbulence: Scaling theory versus selective decay.

*Phys. Fluids A,***4,**1314–1316.Chavanis, P. H., and J. Sommeria, 1996: Classification of self-organized vortices in two-dimensional turbulence: The case of a bounded domain.

*J. Fluid Mech.,***314,**267–297.——, and ——, 1988: Classification of robust isolated vortices in two-dimensional hydrodynamics.

*J. Fluid Mech.,***356,**259–296.Chomaz, J. M., M. Rabaud, C. Basdevant, and Y. Couder, 1988: Experimental and numerical investigation of a forced circular shear layer.

*J. Fluid Mech.,***187,**115–140.Dritschel, D. G., 1986: The nonlinear evolution of rotating configurations of uniform vorticity.

*J. Fluid Mech.,***172,**157–182.——, 1988: Contour surgery: A topological reconnection scheme for extended integrations using contour dynamics.

*J. Comput. Phys.,***77,**240–266.——, 1989: On the stabilization of a two-dimensional vortex strip by adverse shear.

*J. Fluid Mech.,***206,**193–221.Eliassen, A., and M. Lystad, 1977: The Ekman layer of a circular vortex. A numerical and theoretical study.

*Geophys. Norv.,***31,**1–16.Emanuel, K. A., 1997: Some aspects of hurricane inner-core dynamics and energetics.

*J. Atmos. Sci.,***54,**1014–1026.Fine, K. S., A. C. Cass, W. G. Flynn, and C. F. Driscoll, 1995: Relaxation of 2D turbulence to vortex crystals.

*Phys. Rev. Lett.,***75,**3277–3280.Finley, C., 1997: Numerical simulation of intense multi-scale vortices generated by supercell thunderstorms. Ph.D. dissertation, Dept. of Atmospheric Science, Colorado State University, 297 pp. [Available from Dept. of Atmospheric Science, Colorado State University, Fort Collins, CO 80523.].

Fox, C., 1987:

*An Introduction to the Calculus of Variations.*Dover, 271 pp.Fujita, T. T., K. Watanabe, K. Tsuchiya, and M. Shimada, 1972: Typhoon-associated tornadoes in Japan and new evidence of suction vortices in a tornado near Tokyo.

*J. Meteor. Soc. Japan,***50,**431–453.Gall, R. L., 1982: Internal dynamics of tornado-like vortices.

*J. Atmos. Sci.,***39,**2721–2736.——, 1983: A linear analysis of the multiple vortex phenomenon in simulated tornadoes.

*J. Atmos. Sci.,***40,**2010–2024.Guinn, T. A., 1992: A dynamical theory for hurricane spiral bands. Ph.D. dissertation, Dept. of Atmospheric Science, Colorado State University, 178 pp. [Available from Dept. of Atmospheric Science, Colorado State University, Fort Collins, CO 80523.].

——, and W. H. Schubert, 1993: Hurricane spiral bands.

*J. Atmos. Sci.,***50,**3380–3403.——, and W. H. Schubert, 1994: Reply.

*J. Atmos. Sci.,***51,**3545–3546.Hasler, A. F., P. G. Black, V. M. Karyampudi, M. Jentoft-Nilsen, K. Palaniappan, and D. Chesters, 1997: Synthesis of eyewall mesovortex and supercell convective structures in Hurricane Luis with

*GOES-8/9*stereo, concurrent 1-min*GOES-9*and NOAA airborne radar observations. Preprints,*22d Conf. on Hurricanes and Tropical Meteorology,*Ft. Collins, CO, Amer. Meteor. Soc., 201–202.Hawkins, H. F., and S. M. Imbembo, 1976: The structure of a small, intense hurricane—Inez 1966.

*Mon. Wea. Rev.,***104,**418–442.Huang, X.-P., and C. F. Driscoll, 1994: Relaxation of 2D turbulence to a metaequilibrium near the minimum enstrophy state.

*Phys. Rev. Lett.,***72,**2187–2190.Ishioka, K., and S. Yoden, 1994: Nonlinear evolution of a barotropically unstable circumpolar vortex.

*J. Meteor. Soc. Japan,***72,**63–80.——, and ——, 1995: Nonlinear aspects of a barotropically unstable polar vortex in a forced-dissipative system: Flow regimes and tracer transport.

*J. Meteor. Soc. Japan,***73,**201–212.Kurihara, Y., 1976: On the development of spiral bands in a tropical cyclone.

*J. Atmos. Sci.,***33,**940–958.Leith, C. E., 1984: Minimum enstrophy vortices.

*Phys. Fluids,***27,**1388–1395.Lewis, B. M., and H. F. Hawkins, 1982: Polygonal eye walls and rainbands in hurricanes.

*Bull. Amer. Meteor. Soc.,***63,**1294–1300.Lin, S.-J., 1992: Contour dynamics of tornado-like vortices.

*J. Atmos. Sci.,***49,**1745–1756.Matthaeus, W. H., and D. Montgomery, 1980: Selective decay hypothesis at high mechanical and magnetic Reynolds numbers.

*Ann. N. Y. Acad. Sci.,***357,**203–222.McInytre, M. E., 1993: Isentropic distributions of potential vorticity and their relevance to tropical cyclone dynamics.

*ICSU/WMO International Symposium on Tropical Cyclone Disasters,*Z. Zheng, G. Holland, and K. Emanuel, Eds., Peking University Press, 143–156.McWilliams, J. C., 1984: The emergence of isolated coherent vortices in turbulent flow.

*J. Fluid Mech.,***146,**21–43.Michalke, A., and A. Timme, 1967: On the inviscid instability of certain two-dimensional vortex-type flows.

*J. Fluid Mech.,***29,**647–666.Miller, J., 1990: Statistical mechanics of Euler equations in two dimensions.

*Phys. Rev. Lett.,***65,**2137–2140.——, P. B. Weichman, and M. C. Cross, 1992: Statistical mechanics, Euler’s equations, and Jupiter’s red spot.

*Phys. Rev. A,***45,**2328–2359.Möller, J. D., and R. K. Smith, 1994: The development of potential vorticity in a hurricane-like vortex.

*Quart. J. Roy. Meteor. Soc.,***120,**1255–1265.Montgomery, M. T., and L. J. Shapiro, 1995: Generalized Charney-Stern and Fjørtoft theorems for rapidly rotating vortices.

*J. Atmos. Sci.,***52,**1829–1833.——, and R. J. Kallenbach, 1997: A theory for vortex Rossby waves and its application to spiral bands and intensity changes in hurricanes.

*Quart. J. Roy. Meteor. Soc.,***123,**435–465.——, and J. Enagonio, 1998: Tropical cyclogenesis via convectively forced vortex Rossby waves in a three-dimensional quasigeostrophic model.

*J. Atmos. Sci.,***55,**3176–3207.Muramatsu, T., 1986: The structure of polygonal eye of a typhoon.

*J. Meteor. Soc. Japan,***64,**913–921.Ooyama, K., 1969: Numerical simulation of the life cycle of tropical cyclones.

*J. Atmos. Sci.,***26,**3–40.Persing, J., and M. T. Montgomery, 1995: A point-vortex application to vortex stability, evolution and statistical equilibrium. Atmospheric Science Paper 587, Dept. of Atmospheric Science, Colorado State University, 86 pp. [Available from Dept. of Atmospheric Science, Colorado State University, Fort Collins, CO 80523.].

Peurrung, A. J., and J. Fajans, 1993: Experimental dynamics of an annulus of vorticity in a pure electron plasma.

*Phys. Fluids A,***5,**493–499.Ritchie, E., and G. J. Holland, 1993: On the interaction of tropical cyclone scale vortices. Part II: Discrete vortex patches.

*Quart. J. Roy. Meteor. Soc.,***119,**1363–1379.Robert, R., 1991: A maximum-entropy principle for two-dimensional perfect fluid dynamics.

*J. Stat. Phys.,***65,**531–551.——, and J. Sommeria, 1991: Statistical equilibrium states for two-dimensional flows.

*J. Fluid Mech.,***229,**291–310.——, and ——, 1992: Relaxation towards a statistical equilibrium state in two-dimensional perfect fluid dynamics.

*Phys. Rev. Lett.,***69,**2776–2779.Rotunno, R., 1982: A numerical simulation of multiple vortices.

*Intense Atmospheric Vortices,*L. Bengtsson and M. J. Lighthill, Eds., Springer-Verlag, 215–228.——, 1984: An investigation of a three-dimensional asymmetric vortex.

*J. Atmos. Sci.,***41,**283–298.Schade, L. R., 1994: Comments on “Hurricane spiral bands.”

*J. Atmos. Sci.,***51,**3543–3544.Shapiro, L. J., 1996: The motion of Hurricane Gloria: A potential vorticity diagnosis.

*Mon. Wea. Rev.,***124,**2497–2508.——, and H. E. Willoughby, 1982: The response of balanced hurricanes to local sources of heat and momentum.

*J. Atmos. Sci.,***39,**378–394.——, and J. L. Franklin, 1995: Potential vorticity in Hurricane Gloria.

*Mon. Wea. Rev.,***123,**1465–1475.Smith, R. A., and M. N. Rosenbluth, 1990: Algebraic instability of hollow electron columns and cylindrical vortices.

*Phys. Rev. Lett.,***64,**649–652.Snow, J. T., 1978: On inertial instability as related to the multiple-vortex phenomenon.

*J. Atmos. Sci.,***35,**1660–1677.Sommeria, J., C. Staquet, and R. Robert, 1991: Final equilibrium state of a two-dimensional shear layer.

*J. Fluid Mech.,***233,**661–689.Staley, D. O., and R. L. Gall, 1979: Barotropic instability in a tornado vortex.

*J. Atmos. Sci.,***36,**973–981.——, and ——, 1984: Hydrodynamic instability of small eddies in a tornado vortex.

*J. Atmos. Sci.,***41,**422–429.Steffens, J. L., 1988: The effect of vorticity profile shape on the instability of a two-dimensional vortex.

*J. Atmos. Sci.,***45,**254–259.Turkington, B., and N. Whitaker, 1996: Statistical equilibrium computations of coherent structures in turbulent shear layers.

*SIAM J. Sci. Comput.,***17,**1414–1433.Vladimirov, V. A., and V. F. Tarasov, 1980: Formation of a system of vortex filaments in a rotating liquid.

*Izv. Akad. Az. Nauk SSSR, Mekhanika Zhidkosti i Gaza,***1,**44–51.Wakimoto, R. M., and P. G. Black, 1994: Damage survey of hurricane Andrew and its relationship to the eyewall.

*Bull. Amer. Meteor. Soc.,***75,**189–200.Whitaker, N., and B. Turkington, 1994: Maximum entropy states for rotating vortex patches.

*Phys. Fluids,***6,**3963–3973.Willoughby, H. E., 1978: A possible mechanism for the formation of hurricane rainbands.

*J. Atmos. Sci.,***35,**838–848.——, and P. G. Black, 1996: Hurricane Andrew in Florida: Dynamics of a disaster.

*Bull. Amer. Meteor. Soc.,***77,**543–549.——, J. A. Clos, and M. G. Shoreibah, 1982: Concentric eyewalls, secondary wind maxima, and the evolution of the hurricane vortex.

*J. Atmos. Sci.,***39,**395–411.Zabusky, N. J., and E. A. Overman, 1983: Regularization of contour dynamical algorithms.

*J. Comput. Phys.,***52,**351–373.——, M. H. Hughes, and K. V. Roberts, 1979: Contour dynamics for the Euler equations in two dimensions.

*J. Comput. Phys.,***30,**96–106.

## APPENDIX

### Derivation of the MinEV-E Vortex

*λ,*the variational problem is

*υ*(

*b*) =

*υ*

_{0}(

*b*), the third line of (A.1) results from an integration by parts, along with the relation

*δζ*=

*d*(

*rδυ*)/

*rdr,*and the last line makes use of

*δυ*(

*b*) = [

*ζ*

_{0}(

*b*) −

*ζ*(

*b*)]

*δb.*For the independent variation

*δb,*we obtain the transversality condition

*ζ*

*b*

*ζ*

_{0}

*b*

*δυ,*we obtain

*dζ*/

*dr*=

*λυ*for 0 ⩽

*r*⩽

*b,*or

*μ*

^{2}= −

*λ.*The solution of (A.3) is a constant times

*J*

_{1}(

*μr*), where

*J*

_{1}is the first-order Bessel function. The constant is determined by requiring that

*υ*(

*b*) =

*υ*

_{0}(

*b*). Remembering that the flow is unchanged for

*r*≥

*b,*we can write the solution for the tangential wind as

*ζ*(

*r*) =

*d*(

*rυ*)/

*rdr*and

*d*[

*rJ*

_{1}(

*μr*)]/

*rdr*=

*μJ*

_{0}(

*μr*), we obtain from (A.4)

*ζ*(

*b*) =

*ζ*

_{0}(

*b*) for the expression in the upper half of (A.5) yields

*μ*and

*b.*Another relation between

*μ*and

*b*is obtained as follows. Substituting (A.4) into the energy constraint

^{b}

_{0}

*υ*

^{2}

_{0}

*r*)

*r dr*=

^{b}

_{0}

*υ*

^{2}(

*r*)

*r dr*yields

*xJ*

^{2}

_{1}

*x*) =

*d*{½

*x*

^{2}[

*J*

^{2}

_{0}

*x*) +

*J*

^{2}

_{1}

*x*)] −

*xJ*

_{0}(

*x*)

*J*

_{1}(

*x*)}/

*dx*and (A.6), we can perform the integration in (A.7) to obtain

*υ*

_{0}(

*r*) and the associated initial vorticity

*ζ*

_{0}(

*r*), (A.6) and (A.8) determine

*μ*and

*b.*

An important special case of (A.6) and (A.8) occurs when the initial vorticity vanishes at those radii to which mixing extends. In that case the term in braces in (A.8) becomes unity, and *b* is entirely determined by (A.8). In addition, (A.6) simplifies to *J*_{0}(*μb*) = 0, so that *μb* must be one of the zeroes of the *J*_{0} Bessel function. The minimum enstrophy will occur for the smallest possible value of *μ.* Thus, to minimize the enstrophy, we must choose *μb* to be the first zero of *J*_{0}, so that *μb* ≈ 2.4048. From (A.4) it can then be shown that the radius of maximum wind is *r* ≈ 0.7656*b* and that the maximum wind is *υ* ≈ 1.121*υ*_{0}(*b*). Evaluating (A.5) at *r* = 0, we also obtain *ζ*(0) = [*μb*/2*J*_{1}(*μb*)](2*υ*_{0}(*b*)/*b*) ≈ 2.3161(2*υ*_{0}(*b*)/*b*), so that the peak vorticity is approximately 2.3161 times the average vorticity inside *r* = *b.*

We can summarize the MinEV-*υ*_{0}(*r*), first determine *μ* and *b* from (A.6) and (A.8). The final adjusted tangential wind profile *υ*(*r*) and vorticity profile *ζ*(*r*) are then given by (A.4) and (A.5).

Isolines of the maximum dimensionless growth rate *ν*_{i}/*ζ*_{av} among the tangential wavenumbers *m* = 3, 4, . . . , 8. The displayed isolines are the same as in Fig. 1, and again the largest growth rates occur in the lower-right corner of the figure. Shading indicates the wavenumber associated with the maximum dimensionless growth rate at each point.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Isolines of the maximum dimensionless growth rate *ν*_{i}/*ζ*_{av} among the tangential wavenumbers *m* = 3, 4, . . . , 8. The displayed isolines are the same as in Fig. 1, and again the largest growth rates occur in the lower-right corner of the figure. Shading indicates the wavenumber associated with the maximum dimensionless growth rate at each point.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Isolines of the maximum dimensionless growth rate *ν*_{i}/*ζ*_{av} among the tangential wavenumbers *m* = 3, 4, . . . , 8. The displayed isolines are the same as in Fig. 1, and again the largest growth rates occur in the lower-right corner of the figure. Shading indicates the wavenumber associated with the maximum dimensionless growth rate at each point.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Vorticity contour plots for the representative numerical experiment. The model domain is 600 km × 600 km, but only the inner 200 km × 200 km is shown. The contours begin at 0.0005 s^{−1} and are incremented by 0.0005 s^{−1}. Low vorticity values are shaded blue and high vorticity values are shaded red. (a) Vorticity from *t* = 0 h to 8 h.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Vorticity contour plots for the representative numerical experiment. The model domain is 600 km × 600 km, but only the inner 200 km × 200 km is shown. The contours begin at 0.0005 s^{−1} and are incremented by 0.0005 s^{−1}. Low vorticity values are shaded blue and high vorticity values are shaded red. (a) Vorticity from *t* = 0 h to 8 h.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Vorticity contour plots for the representative numerical experiment. The model domain is 600 km × 600 km, but only the inner 200 km × 200 km is shown. The contours begin at 0.0005 s^{−1} and are incremented by 0.0005 s^{−1}. Low vorticity values are shaded blue and high vorticity values are shaded red. (a) Vorticity from *t* = 0 h to 8 h.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

(*Continued*) (b) Vorticity from *t* = 10 h to 20 h.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

(*Continued*) (b) Vorticity from *t* = 10 h to 20 h.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

(*Continued*) (b) Vorticity from *t* = 10 h to 20 h.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

(*Continued*) (c) Vorticity from *t* = 22 h to 48 h with the time interval switched to 6 h after *t* = 24 h.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

(*Continued*) (c) Vorticity from *t* = 22 h to 48 h with the time interval switched to 6 h after *t* = 24 h.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

(*Continued*) (c) Vorticity from *t* = 22 h to 48 h with the time interval switched to 6 h after *t* = 24 h.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Azimuthal mean vorticity (*ζ**υ**ω**t* = 0 (solid), 4 h (dotted), 8 h (dashed), 12 h (dash–dot), 24 h (dash–dot–dot–dot), 48 h (long dashes). The scale on the right of the bottom panel is for 2*π*/*ω*

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Azimuthal mean vorticity (*ζ**υ**ω**t* = 0 (solid), 4 h (dotted), 8 h (dashed), 12 h (dash–dot), 24 h (dash–dot–dot–dot), 48 h (long dashes). The scale on the right of the bottom panel is for 2*π*/*ω*

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Azimuthal mean vorticity (*ζ**υ**ω**t* = 0 (solid), 4 h (dotted), 8 h (dashed), 12 h (dash–dot), 24 h (dash–dot–dot–dot), 48 h (long dashes). The scale on the right of the bottom panel is for 2*π*/*ω*

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Kinetic energy (*E*(*k*)) and enstrophy (*Z*(*k*)) spectra during the extensive mixing phase (0 ⩽ *t* ⩽ 16 h) for the experiment shown in Fig. 3 at the selected times *t* = 0 (solid), 4 h (dotted), 8 h (dashed), 12 h (dash–dot), and 16 h (dash–dot–dot–dot). Spectra are obtained in the usual way by binning into rings of radius *k* centered at the origin in Cartesian wavenumber space. The *k*^{−3} (energy) and *k*^{−1} (enstrophy) spectra, as expected for the enstrophy cascade from two-dimensional turbulence theory (neglecting logarithmic corrections), are shown for comparison.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Kinetic energy (*E*(*k*)) and enstrophy (*Z*(*k*)) spectra during the extensive mixing phase (0 ⩽ *t* ⩽ 16 h) for the experiment shown in Fig. 3 at the selected times *t* = 0 (solid), 4 h (dotted), 8 h (dashed), 12 h (dash–dot), and 16 h (dash–dot–dot–dot). Spectra are obtained in the usual way by binning into rings of radius *k* centered at the origin in Cartesian wavenumber space. The *k*^{−3} (energy) and *k*^{−1} (enstrophy) spectra, as expected for the enstrophy cascade from two-dimensional turbulence theory (neglecting logarithmic corrections), are shown for comparison.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Kinetic energy (*E*(*k*)) and enstrophy (*Z*(*k*)) spectra during the extensive mixing phase (0 ⩽ *t* ⩽ 16 h) for the experiment shown in Fig. 3 at the selected times *t* = 0 (solid), 4 h (dotted), 8 h (dashed), 12 h (dash–dot), and 16 h (dash–dot–dot–dot). Spectra are obtained in the usual way by binning into rings of radius *k* centered at the origin in Cartesian wavenumber space. The *k*^{−3} (energy) and *k*^{−1} (enstrophy) spectra, as expected for the enstrophy cascade from two-dimensional turbulence theory (neglecting logarithmic corrections), are shown for comparison.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Contour plots of the *Q* field (*Q* = *ψ*^{2}_{xy}*ψ*_{xx}*ψ*_{yy}) at selected times for the experiment shown in Fig. 3. The contour interval is 0.5 × 10^{−6} s^{−2}, with the zero isoline omitted for clarity (to suppress the effect of small oscillations about zero). The shading (where *Q* ≥ 0.5 × 10^{−6} s^{−2}) corresponds to regions of particle separation.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Contour plots of the *Q* field (*Q* = *ψ*^{2}_{xy}*ψ*_{xx}*ψ*_{yy}) at selected times for the experiment shown in Fig. 3. The contour interval is 0.5 × 10^{−6} s^{−2}, with the zero isoline omitted for clarity (to suppress the effect of small oscillations about zero). The shading (where *Q* ≥ 0.5 × 10^{−6} s^{−2}) corresponds to regions of particle separation.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Contour plots of the *Q* field (*Q* = *ψ*^{2}_{xy}*ψ*_{xx}*ψ*_{yy}) at selected times for the experiment shown in Fig. 3. The contour interval is 0.5 × 10^{−6} s^{−2}, with the zero isoline omitted for clarity (to suppress the effect of small oscillations about zero). The shading (where *Q* ≥ 0.5 × 10^{−6} s^{−2}) corresponds to regions of particle separation.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Schematic of an initial vorticity distribution and the corresponding hypothesized end-state after a redistribution without mixing. The redistribution argument illustrated here is flawed because it violates kinetic energy and angular momentum invariance. This points out the necessity of vorticity mixing during the redistribution process.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Schematic of an initial vorticity distribution and the corresponding hypothesized end-state after a redistribution without mixing. The redistribution argument illustrated here is flawed because it violates kinetic energy and angular momentum invariance. This points out the necessity of vorticity mixing during the redistribution process.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Schematic of an initial vorticity distribution and the corresponding hypothesized end-state after a redistribution without mixing. The redistribution argument illustrated here is flawed because it violates kinetic energy and angular momentum invariance. This points out the necessity of vorticity mixing during the redistribution process.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Time series of kinetic energy (

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Time series of kinetic energy (

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Time series of kinetic energy (

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Plots of the azimuthal mean vorticity *ζ*(*r*), tangential wind *υ*(*r*), pressure *p*(*r*), as determined from gradient balance, and angular velocity *ω*(*r*) for MinEV-

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Plots of the azimuthal mean vorticity *ζ*(*r*), tangential wind *υ*(*r*), pressure *p*(*r*), as determined from gradient balance, and angular velocity *ω*(*r*) for MinEV-

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Plots of the azimuthal mean vorticity *ζ*(*r*), tangential wind *υ*(*r*), pressure *p*(*r*), as determined from gradient balance, and angular velocity *ω*(*r*) for MinEV-

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Isolines of the entropy density −*ρ*_{1} ln*ρ*_{1} − *ρ*_{2} ln*ρ*_{2} − (1 − *ρ*_{1} − *ρ*_{2}) ln(1 − *ρ*_{1} − *ρ*_{2}) in the (*ρ*_{1}, *ρ*_{2})-plane for *ρ*_{1} ≥ 0 and *ρ*_{2} ≥ 0 with *ρ*_{1} + *ρ*_{2} ⩽ 1. The isolines are 0.1, 0.2, . . . , 1.0 and the asterisk indicates the maximum value of entropy density (ln 3), which occurs at *ρ*_{1} = *ρ*_{2} = ^{1}/_{3}.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Isolines of the entropy density −*ρ*_{1} ln*ρ*_{1} − *ρ*_{2} ln*ρ*_{2} − (1 − *ρ*_{1} − *ρ*_{2}) ln(1 − *ρ*_{1} − *ρ*_{2}) in the (*ρ*_{1}, *ρ*_{2})-plane for *ρ*_{1} ≥ 0 and *ρ*_{2} ≥ 0 with *ρ*_{1} + *ρ*_{2} ⩽ 1. The isolines are 0.1, 0.2, . . . , 1.0 and the asterisk indicates the maximum value of entropy density (ln 3), which occurs at *ρ*_{1} = *ρ*_{2} = ^{1}/_{3}.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

Isolines of the entropy density −*ρ*_{1} ln*ρ*_{1} − *ρ*_{2} ln*ρ*_{2} − (1 − *ρ*_{1} − *ρ*_{2}) ln(1 − *ρ*_{1} − *ρ*_{2}) in the (*ρ*_{1}, *ρ*_{2})-plane for *ρ*_{1} ≥ 0 and *ρ*_{2} ≥ 0 with *ρ*_{1} + *ρ*_{2} ⩽ 1. The isolines are 0.1, 0.2, . . . , 1.0 and the asterisk indicates the maximum value of entropy density (ln 3), which occurs at *ρ*_{1} = *ρ*_{2} = ^{1}/_{3}.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

The MaxSV expectation functions *ρ*_{1}(*r*) (dashed), *ρ*_{2}(*r*) (dotted), and *ρ*_{3}(*r*) = 1 − *ρ*_{1}(*r*) − *ρ*_{2}(*r*) (solid), as computed from (6.4a)–(6.4g) for the initial condition displayed by the thick solid curves in Fig. 12.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

The MaxSV expectation functions *ρ*_{1}(*r*) (dashed), *ρ*_{2}(*r*) (dotted), and *ρ*_{3}(*r*) = 1 − *ρ*_{1}(*r*) − *ρ*_{2}(*r*) (solid), as computed from (6.4a)–(6.4g) for the initial condition displayed by the thick solid curves in Fig. 12.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

The MaxSV expectation functions *ρ*_{1}(*r*) (dashed), *ρ*_{2}(*r*) (dotted), and *ρ*_{3}(*r*) = 1 − *ρ*_{1}(*r*) − *ρ*_{2}(*r*) (solid), as computed from (6.4a)–(6.4g) for the initial condition displayed by the thick solid curves in Fig. 12.

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

A comparison of the predictions of MinEV-

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

A comparison of the predictions of MinEV-

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

A comparison of the predictions of MinEV-

Citation: Journal of the Atmospheric Sciences 56, 9; 10.1175/1520-0469(1999)056<1197:PEAECA>2.0.CO;2

^{1}

Throughout this paper the symbols *u, υ* are used to denote radial and tangential components of velocity when working in cylindrical coordinates and eastward and northward components of velocity when working in Cartesian coordinates.

^{2}

In the analogous problem possessing a continuous basic-state vorticity profile, a) *m* = 2 can become exponentially unstable and b) an algebraic continuous-spectrum instability for *m* = 1 also exists and becomes nonlinear at large times even if the profile supports no exponential instabilities for *m* = 2, 3, . . . (Smith and Rosenbluth 1990). Whether a similar algebraic instability exists in the piecewise constant model when sheared disturbances are allowed is an interesting, unresolved question.

^{3}

For exceptional initial conditions (i.e., those that do not contain an azimuthal wavenumber 1 contribution nor excite wavenumber 1 through nonlinear interaction) the macroscopic results (not shown) are observed to be qualitatively much different. In such cases the vorticity ring never mixes completely into a monopole whose largest vorticity resides at the center. For example, if the initial perturbation is in mode 4 only, then the nonlinear interactions excite modes 8, 12, 16, . . . . The high vorticity in the ring thus never gets advected to the storm center due to the absence of a wavenumber 1 wind component. [The relative contribution from the mode 1 algebraic instability (see footnote 1) to the vorticity mixing process reported here for the general case of a broadbanded initial disturbance has not yet been determined.] For such exceptional initial conditions the mathematical theory developed in sections 5 and 6 is no longer applicable as it implicitly assumes no mixing barriers.

^{4}

Since the model is on an *f* plane, the centroid of vorticity is invariant and always remains at the domain center.

^{5}

The variational formulation developed here and in the sequel is a slight generalization of Leith’s (1984) formulation since it does not require the flow to be irrotational outside the mixing region and obtains the edge of the mixing region in terms of the initial vortex.

^{6}

In general, conservation of angular momentum and circulation are not sufficient to guarantee a unique solution to (5.8). For the initial vortex discussed above one finds an infinite number of roots inside the high vorticity region (*a* ⩽ 46 km) and a distinct root outside (*a* ≈ 82 km). The root yielding the maximum enstrophy deficit (the distinct root) is deemed the solution to the MinEV-

^{7}

Upon taking the second variation it can be shown that the MinEV-

^{8}

This fact is readily demonstrated for both the MinEV-

^{9}

Note added in proof: A recent paper by Chavanis and Sommeria (1998) has generalized Leith’s MinEV theory to simultaneously enforce the circulation, angular momentum, and energy constraints. They have shown that the resulting relation between vorticity and streamfunction is the same as the linearized version of the corresponding MaxSV relation (e.g., the linearized version of (6.4a)–(6.4c) for the three region case), thus establishing a connection between MinEV and MaxSV theories in the limit of strong mixing.