## 1. Introduction

Tropical cyclone modeling has now advanced to a stage where it is possible to predict, using three-dimensional nonhydrostatic nested-grid models, not only the synoptic-scale flow patterns that determine the cyclone track, but also the mesoscale inner core vorticity structures that determine the maximum wind and minimum central pressure (Chen and Yau 2001; Braun 2002). An observational basis for understanding the inner core vorticity structure was recently provided by Kossin and Eastin (2001), who used aircraft data to show that annular rings of high vorticity are generated in the hurricane eyewall during intensification but are typically eradicated when intensification stops. During prolonged periods of intensification, eyewall flows appear to evolve toward vortex sheets, in agreement with the theoretical notion that the dynamics of the eyewall are frontogenetic and that there is a tendency for the azimuthal flow to evolve toward a discontinuity (Emanuel 1989, 1997). Motivated by such observed flow fields in hurricanes, a theoretical and modeling study by Schubert et al. (1999) considered the barotropic stability and nonlinear evolution of unforced “hollow-tower” vorticity structures. Such structures are characterized by annular rings of enhanced vorticity (representing the eyewall) surrounding a region of relatively weak vorticity (representing the eye). Their study suggests that remarkable vorticity rearrangement processes are occurring in the inner core regions of tropical cyclones, and that these processes can have a significant effect on the structure, intensity, and evolution of the storms. Thus the simultaneous simulation of tropical cyclone track and inner core intensity remains a challenge, at least partly because of the very different, but comparably important, horizontal scales involved. In particular, it is quite possible to incorporate horizontal diffusion in such a way that it has little effect on the scales of motion involved in track prediction but has a dominant effect on the inner core vortex structure.

The purpose of this note is to illustrate the effects of diffusion on hollow-tower vorticity structures and to determine how small a value of viscosity we should specify in three-dimensional hurricane models in order that diffusion does not erroneously dominate the advective vorticity rearrangement process. We will also comment on some of the problems associated with the use of horizontal diffusion as a parameterization of the advective vorticity rearrangement process.

## 2. Diffusive evolution of a circular vortex sheet

*ζ*(

*r*,

*t*) iswhere

*ν*is the kinematic viscosity. As the initial condition, consider a vortex that is initially characterized by an infinitesimally thin annular ring of enhanced vorticity embedded in an irrotational flow. The initial azimuthal wind profile is given by

*υ*(

*r*, 0) = 0 for 0 ≤

*r*<

*a*, and

*υ*(

*r*, 0) =

*υ*

_{0}(

*a*/

*r*) for

*a*≤

*r*< ∞, where

*υ*

_{0}is the initial maximum wind. This describes a circular vortex “sheet” with radius

*r*=

*a*.

*ζ*(

*r*,

*t*) is given by

*ζ̂*

*k*,

*t*) =

*ζ*(

*r*,

*t*)

*J*

_{0}(

*kr*)

*r*

*dr*, and the inverse transform is given by

*ζ*(

*r*,

*t*) =

*ζ̂*

*k*,

*t*)

*J*

_{0}(

*kr*)

*k*

*dk*, where

*J*

_{0}is the zero-order Bessel function and

*k*is the radial wavenumber. Applying the transform to (1), integrating by parts twice, and using the zero-order Bessel equation

*rd*/

*dr*[

*rdJ*

_{0}(

*kr*)/

*dr*] = −

*k*

^{2}

*r*

^{2}

*J*

_{0}(

*kr*), we obtain

*d*

*ζ̂*

*dt*= −

*ν*

*k*

^{2}

*ζ̂*

*ζ̂*

*k*,

*t*) =

*ζ̂*

*k*, 0) exp(−

*ν*

*k*

^{2}

*t*). Substituting this last result into the inverse transform gives

*ζ*(

*r*,

*t*) =

*ζ̂*

*k*, 0) exp(−

*ν*

*k*

^{2}

*t*)

*J*

_{0}(

*kr*)

*k*

*dk*, where

*ζ̂*

*k*, 0) =

*ζ*(

*r*, 0)

*J*

_{0}(

*kr*)

*r*

*dr*=

*d*[

*r*

*υ*(

*r*, 0)]/

*dr*}

*J*

_{0}(

*kr*)

*dr*=

*a*

*υ*

_{0}

*J*

_{0}(

*ka*) in the limit

*ε*→ 0. Using these last two results we obtainwhich is an integral representation for the solution of our diffusion problem. The integral in (2) can be evaluated (see Gradshteyn and Ryzhik 1994, p. 739) to obtainwhere

*I*

_{0}is the modified Bessel function of order zero, and

*τ*=

*ν*

*t*/

*a*

^{2}is a dimensionless time.

Figure 1a, which has been constructed using (3), shows plots of the dimensionless vorticity *ζ*(*r*, *t*)/(*υ*_{0}/*a*) as a function of *r*/*a* for seven different values of the dimensionless time *τ*. To construct the corresponding wind plots shown in Fig. 1b we have used numerical quadrature on *υ*(*r*, *t*) = *r*^{–1} *ζ*(*r*′, *t*)*r*′ *dr*′ with the integrand *ζ*(*r*′, *t*)*r*′ computed using (3). Figure 1 shows how the diffusive flux of vorticity smooths the profiles and drives the vorticity from a hollow-tower to a monotonic profile. The tangential wind profile becomes increasingly less “U shaped” as the flow increases in the eye and decreases in the eyewall.

For *r* = 0 and *t* > 0, (3) reduces to *ζ*(0, *τ*) = (*υ*_{0}/*a*)[1/(2*τ*)] exp[−1/(4*τ*)], from which we can plot the dimensionless central vorticity *ζ*(0, *τ*)/(*υ*_{0}/*a*) as a function of *τ*. This is shown in Fig. 2. The central flow remains approximately irrotational until *τ* ≈ 0.03, after which the central vorticity increases until it achieves its maximum at *τ* = 0.25. When *τ* = 0.25, the Laplacian operator on the right-hand side of (1), evaluated at *r* = 0, vanishes, and the central vorticity is (2*υ*_{0}/*a*)*e*^{–1}, that is, *e*^{–1} times the initial average vorticity inside *r* = *a*. When *τ* > 0.25 the central vorticity decreases as all the diffusive flux is then directed outward. The increase of the central vorticity due to the inward flux when *τ* < 0.25 is much more rapid than the “slow diffusive spindown” that occurs later. In other words, the diffusion is more efficient at eradicating the annular ring than the subsequent monopole. Since *τ* = 0.25 is the dimensionless diffusive time to monopole formation, we can define the dimensional diffusive time to monopole formation as *t*_{mf} = *a*^{2}/(4*ν*). For *a* = 10 km, *t*_{mf} is, respectively, 694, 69.4, and 6.94 h for *ν* = 10, 100, and 1000 m^{2} s^{–1}, a range of eddy viscosity values typically used in hurricane models.

## 3. Advective–diffusive evolution of a circular vortex sheet

*a*= 10 km, the vorticity and wind profiles at

*τ*= 0.008 in Fig. 1 are very similar to the initial condition used by Kossin and Schubert (2001) in experiment 7 of their vorticity mixing study. In their results, reproduced here in Fig. 3, the initial annular ring of vorticity advectively rearranges itself into a vorticity monopole in less than 8 h, according to the nondivergent barotropic dynamics:

*ζ*

_{t}

*u*

*ζ*

_{x}

*υ*

*ζ*

_{y}

*ν*

^{2}

*ζ*

*u*= −

*ψ*

_{y},

*υ*=

*ψ*

_{x}, and

*ζ*= ∇

^{2}

*ψ*. Their results are nearly inviscid since their specified viscosity is

*ν*= 5 m

^{2}s

^{–1}, which gives a dimensional diffusive time to monopole formation of

*t*

_{mf}=

*a*

^{2}/(4

*ν*) ≈ 1389 h, much longer than the roughly 7 h required for advective rearrangement to a monopole. Values of

*ν*somewhat larger than

*ν*= 5 m

^{2}s

^{–1}would also seem acceptable. For example,

*ν*= 100 m

^{2}s

^{–1}gives

*t*

_{mf}= 69.4 h, which is about 10 times the time required for advective rearrangement.

Now suppose we want to parameterize the asymmetric advective rearrangement process in an axisymmetric model using simple diffusion. What value of *ν* should be used? If we want the timescale of monopole formation to be correctly parameterized in the axisymmetric model, we should choose *ν* to satisfy *ν* = *a*^{2}/(4*t*_{mf}) = (10 km)^{2}/(4 × 7 h) ≈ 1000 m^{2} s^{–1}. However, this simple parameterization has serious deficiencies since, during the advective process, the peak value of vorticity remains unchanged (see Fig. 3), while during the pure diffusive process illustrated in Fig. 1, the peak vorticity at *τ* = 0.25 is only 24% of the peak value at *τ* = 0.008. It is in this regard that the results of advective rearrangement and the results of diffusion are quite different. We can summarize by saying that the purely diffusive evolution of an annular ring does not accurately parameterize the process of vorticity rearrangement, particularly in the hurricane inner core.

*ζ*

^{2}

*dx*

*dy*is selectively decayed over the global kinetic energy

*u*

^{2}+

*υ*

^{2})

*dx*

*dy*. This is a consequence of the fact that, according to the governing dynamics (4),

*t*) and

*t*) obey

*d*

*dt*= −2

*ν*

*d*

*dt*= −2

*ν*

*ζ*· ∇

*ζ*

*dx*

*dy*is the global palinstrophy, a measure of the overall vorticity gradient. Global palinstrophy evolution is described byAlthough the last term in (5) always tends to reduce

*ν*,

*d*

*dt*= −

*ν*2

*π*∫ {∂[

*r*(∂

*ζ*/∂

*r*)]/(

*r*∂

*r*)}

^{2}

*r*

*dr*, which is the nonadvective, axisymmetric, special case of (5). Thus, in the pure diffusive case,

Although pure diffusion is not consistent with selective decay of enstrophy over energy, there are a number of parameterizations that are consistent. The anticipated vorticity method for the nondivergent barotropic model (Leith 1985), and its generalization to the shallow water equations and to the quasi-static primitive equations [referred to as the anticipated potential vorticity method (Sadourny and Basdevant 1985)] are formed by recasting the equations of motion using a modified diffusion that conserves energy exactly and allows dissipation of enstrophy or potential enstrophy. In addition, certain spatial discretization schemes (Arakawa and Hsu 1990) have been designed to mimic the anticipated potential vorticity method and thus to be consistent with selective decay. The maximum entropy production principle (Kazantsev et al. 1998) is a statistical method that relies on solution of the variational problem of globally maximizing the rate of production of Boltzmann mixing entropy under the physical constraints of the system (e.g., global energy conservation), while allowing dissipation of enstrophy. These parameterizations could be considered an improvement over pure diffusion, but they can still suffer from deficiencies when modeling the dramatic vorticity rearrangements that can occur in the hurricane inner core. For example, the maximum entropy solutions to barotropically unstable initial flows do not capture the process whereby coherent vorticity structures are transported via advection while remaining relatively unmixed in the process (e.g., Jin and Dubin 1998) and thus the method tends to underpredict (overly mix) vorticity values.

## Acknowledgments

We thank Kerry Emanuel for his comments on the original manuscript. This work was supported by NASA/CAMEX Contract NAG5-11010, NSF Grant ATM-0087072, and NOAA Grant NA67RJ0152.

## REFERENCES

Arakawa, A., , and Y-J. G. Hsu, 1990: Energy conserving and potential-enstrophy dissipating schemes for the shallow water equations.

,*Mon. Wea. Rev.***118****,**1960–1969.Braun, S. A., 2002: A cloud-resolving simulation of Hurricane Bob (1991): Storm structure and eyewall buoyancy.

,*Mon. Wea. Rev.***130****,**1573–1592.Chen, Y., , and M. K. Yau, 2001: Spiral bands in a simulated hurricane. Part I: Vortex Rossby wave verification.

,*J. Atmos. Sci.***58****,**2128–2145.Emanuel, K. A., 1989: The finite-amplitude nature of tropical cyclogenesis.

,*J. Atmos. Sci.***46****,**3431–3456.Emanuel, K. A., 1997: Some aspects of hurricane inner-core dynamics and energetics.

,*J. Atmos. Sci.***54****,**1014–1026.Gradshteyn, I. S., , and I. M. Ryzhik, 1994:

*Tables of Integrals, Series, and Products*. 5th ed. Academic Press, 1204 pp.Jin, D. Z., , and D. H. E. Dubin, 1998: Regional maximum entropy theory of vortex crystal formation.

,*Phys. Rev. Lett.***80****,**4434–4437.Kazantsev, E., , J. Sommeria, , and J. Verron, 1998: Subgrid-scale eddy parameterization by statistical mechanics in a barotropic ocean model.

,*J. Phys. Oceanogr.***28****,**1017–1042.Kossin, J. P., , and M. D. Eastin, 2001: Two distinct regimes in the kinematic and thermodynamic structure of the hurricane eye and eyewall.

,*J. Atmos. Sci.***58****,**1079–1090.Kossin, J. P., , and W. H. Schubert, 2001: Mesovortices, polygonal flow patterns, and rapid pressure falls in hurricane-like vortices.

,*J. Atmos. Sci.***58****,**2196–2209.Leith, C. E., 1985: Two-dimensional coherent structures.

*Proceedings of the International School of Physics “Enrico Fermi,” Course 88: Turbulence and Predictability in Geophysical Fluid Dynamics and Climate Dynamics*, M. Ghil, R. Benzi, and G. Parisi, Eds., North-Holland Physics Publishing, 266–280.Sadourny, R., , and C. Basdevant, 1985: Parameterization of subgrid scale barotropic and baroclinic eddies in quasi-geostrophic models: Anticipated potential vorticity method.

,*J. Atmos. Sci.***42****,**1353–1363.Schubert, W. H., , M. T. Montgomery, , R. K. Taft, , T. A. Guinn, , S. R. Fulton, , J. P. Kossin, , and J. P. Edwards, 1999: Polygonal eyewalls, asymmetric eye contraction, and potential vorticity mixing in hurricanes.

,*J. Atmos. Sci.***56****,**1197–1223.