1. Introduction
A complete investigation of hurricane intensity requires as a vital component (to paraphrase Emanuel 1988) an understanding of the processes that determine the upper bound on intensity. This upper bound, for a given set of environmental conditions, is the maximum possible intensity, customarily referred to as maximum potential intensity, or MPI. If an MPI principle can be found, then the manifold of possible intensities of the hurricane can be constrained. Camp (1999) provides an in-depth review of all the known theoretical approaches for constructing various MPI models.
Camp and Montgomery (2001) argue that the MPI theory of Emanuel “is the closest to providing a useful calculation of maximum intensity.” This MPI theory was first introduced in Emanuel (1986, which shares some similarities to the approach of Kleinschmidt 1951) employing a frictional boundary layer beneath a conditionally neutral outflow layer. Within this view, the energy for the storm comes from enhanced sea-to-air exchanges of enthalpy in the hurricane boundary layer due to enhanced wind speeds. The Emanuel (1986, 1988) theory, which was presented as a rigorous upper bound, and subsequent modifications [improved eye parameterization (Emanuel 1995b); inclusion of dissipational heating (Bister and Emanuel 1998)] will be referred to here as E-MPI. E-MPI theory has gained much acceptance in the community (Bosart et al. 2000; AMS Council 2000; Bengtsson 2001) and is widely considered a paradigm.
E-MPI is formulated on two-dimensional, axisymmetric principles, and neglects certain processes such as vertical wind shear (Jones 1995; Wang and Holland 1996; Frank and Ritchie 2001), convective asymmetries, spiral arm bands, vortex Rossby waves (Guinn and Schubert 1993; Montgomery and Kallenbach 1997; Schubert et al. 1999; Möller and Montgomery 2000; Chen and Yau 2001; Wang 2002a,b), ocean spray (Fairall et al. 1994; Andreas and Emanuel 2001), and wind-induced ocean cooling (Shay et al. 1998; Bender and Ginis 2000; Jacob et al. 2000). Given these simplifications, one might anticipate that E-MPI can only serve as an approximate upper bound on real hurricanes that are three-dimensional and include all of these unmodeled processes. On the other hand, in an axisymmetric computer model that has a simplified one-way ocean–air interaction one might anticipate that E-MPI could be an exact upper-bound on the intensity of the resulting axisymmetric storms. Indeed, E-MPI theory was previously “validated” against results from such axisymmetric models (Rotunno and Emanuel 1987; Emanuel 1989, 1995a).
Using a nonhydrostatic, many-layered, axisymmetric hurricane model (based on Ooyama 2001) with standard tropical parameters, Hausman (2001) recently demonstrated a systematic increase of storm intensity with model resolution, with convergence found at approximately 1 km horizontal grid spacing [using 28°C sea surface temperature (SST)] with a maximum sustained tangential velocity of 140 m s–1. Needless to say, this greatly exceeds the E-MPI for this model storm with any reasonable estimation of an outflow temperature from a Jordan (1958) mean sounding. It is worth noting that, long ago, Yamasaki (1983) reported an axisymmetric simulation with >90 m s–1 with 28.85°C SST that also exceeds E-MPI. Because of these findings and the fundamental nature of the problem, we believe a reexamination and a more thorough testing of E-MPI theory is warranted.
This paper will demonstrate a hitherto unreported resolution sensitivity in the Rotunno and Emanuel (1987) model that results in a hurricane vortex that greatly exceeds E-MPI. To explain the simulated phenomenon, we examine each of the assumptions of E-MPI and find that the underlying philosophy of E-MPI, namely that there exists a point balance of entropy and angular momentum at the top of the subcloud layer at the base of the eyewall, cannot explain the observed “superintensity.” We will show that the addition of heat from the eye to the eyewall, which violates two key assumptions of E-MPI (that there is no such interaction and that the eyewall is neutral to moist ascent), is responsible for the enhanced wind speeds.
The plan of the paper is as follows. Section 2 briefly reviews the Rotunno and Emanuel 1987 model and our numerical experiments with it. Section 3 briefly reviews E-MPI. Section 4 documents the behavior of high-resolution simulations that produce very intense storms compared to E-MPI. Section 5 documents the evolution of air with high entropy in the simulated eye and demonstrates that this air interacts with the eyewall. Section 6 provides evidence that the pertinent element of the eye–eyewall interaction is as a source of latent heat rather than of local buoyancy. Section 7 provides supporting evidence of the superintensity mechanism from observations and recent three-dimensional cloud-resolving simulations. Section 8 provides a discussion. Section 9 presents our conclusions.
2. Rotunno and Emanuel model reviewed
The Rotunno and Emanuel (1987, hereafter RE87) hurricane model is an axisymmetric (two-dimensional), nonhydrostatic, cloud-resolving version of the Klemp and Wilhelmson (1978) formulation originally developed for the study of cumulus clouds and supercell storms. This model was kindly provided to us by R. Rotunno of the National Center for Atmospheric Research (NCAR). The model is integrated on a staggered C grid using a fixed radial and fixed vertical grid spacing. There is no ice in the model and explicit convection is employed with a fixed precipitation fallspeed of 7 m s–1. Subgrid-scale turbulence is parameterized using a standard Smagorinsky (1963) formulation modified to allow a variable mixing length, a function of the local stability using a Richardson-number-dependent closure following Lilly (1962). Radiation is represented crudely by a standard Newtonian cooling parameterization that tends to restore the local potential temperature to that found initially in the ambient environment. In the numerical experiments presented here the radiational heating/cooling rate is capped at 2 K day–1. Sponge layer damping occurs above the model tropopause (z > 19.4 km). The outer boundary at r = 1500 km is open; using a doubled radius for the outer boundary does not give significantly different results. Surface interaction is through a standard bulk aerodynamic parameterization. The surface drag coefficient CD is represented by a linear function of wind speed using Deacon’s formula (Roll 1965). The surface enthalpy exchange coefficient Ck is set equal to the drag coefficient (Ck/CD = 1). The radial and vertical wind fields (u and w, respectively) are computed using a “fast” time step to handle sound waves. All of these processes are unchanged from those described in RE87. The initial thermodynamic sounding of the environment is created by a procedure described in RE87 as a modification of the Jordan (1958) sounding so that it is neutral to convection as realized in the model.
Model parameters used in the default run are summarized in Table 1. These are the same as experiment A in RE87 except there the SST reported is 26.3°C1 and the radiational cooling is uncapped.2 For each doubling of model resolution, the time step is halved, the number of radial and vertical grid points is doubled, the number of integration steps is doubled, and the initial soundings of potential temperature (θ) and vapor mixing ratio (qυ) are interpolated to the new vertical grid. As a consequence of these changes, the radial and vertical grid spacings are halved, the fast time step is halved, and the horizontal mixing length is halved. Simulations will be referred to as the “default run,” “2x run,” “4x run,” and “8x run” with each doubling of radial, vertical, and temporal resolution. For the presentation of results, days (8.64 × 104 s) are used, with the simulations starting at day 0. Figure 1 shows that the 4x run produces a realistic hurricane vortex like Fig. 5 of RE87, except that it is more structured and much more intense, possessing stronger updrafts.
3. E-MPI theory reviewed
a. Foundations
The energetically based MPI theory considered here is that first presented by Emanuel (1986) and updated with an improved representation of the eye in Emanuel (1995b, hereafter E95b) yielding a minimum surface pressure that depends explicitly on the ratio of exchange coefficients of enthalpy and angular momentum. Camp and Montgomery (2001) provided a summary of E-MPI and its physical basis. Bister and Emanuel (1998) updated E-MPI to include the effects of dissipational heating in the boundary layer, but since the model results presented here do not include this process the theory of E95b will be considered.






b. Boundary layer balances


c. A model-appropriate θe






4. Superintensity
The phenonmenon in which axisymmetric hurricane models tend to form more intense storms with increasing model resolution was reviewed by Hausman (2001) and was supplemented by his own simulations using a new thermodynamic formulation by Ooyama (2001). For a 28°C SST simulation at the highest resolution employed (Δr = 500 m) without ice physics, the intensity of the modeled hurricane reached 140 m s–1. Inclusion of ice physics was found to weaken the storm intensity somewhat (120 m s–1), but since this process is not included in E-MPI theory, the subject of our comparison, the simulations lacking ice are the most appropriate for comparison with E-MPI and the RE87 model simulations presented here. This sensitivity to resolution is confirmed here with the entirely different RE87 model (Fig. 3), where the modeled intensity (solid curves) greatly exceeds E-MPI [Eq. (2); dotted curves]. The overlays on Fig. 3 are estimated using daily averaged data. The dashed curves [using Eq. (6), evaluated at the lowest grid level at the radius of maximum winds (RMW)] show a slight increase with increasing intensity, but not nearly enough to explain the modeled superintensity. Use of (6) in this manner does justify in part the use of a fixed
The 2x, 4x, and 8x runs appear to reach a quasi-steady-state intensity by 10 days, with generally stronger steady-state intensities at higher resolution. The default case (Fig. 3a) also appears to complete its initial spinup by day 10, but slowly transitions to a second apparent steady state near 20 days.
Table 2 summarizes the output from each model run. From the first three rows we see evidence for convergence in the ultimate quasi-steady-state vortex intensity with increases in resolution by the 4x run. While MAX(V) and MAX[D(V)] continue to increase slightly for the 8x run, these can be attributed to short-lived events on the order of a day (Fig. 3). MEDIAN(V) for the 8x run is actually smaller than for the 4x run. The “typical” quasi-steady-state intensity can then be considered about the same for the 4x and 8x runs. Hausman (2001) proposed that, at higher resolution, a smaller RMW is resolvable and will result because the limiting factor is the ability of numerical diffusion to mix low angular momentum air outward from the center of the storm, a process whose spatial scale is on the order of the grid scale. Because the RMW is found to be the approximate limit for inward advection of environmental angular momentum, a smaller RMW at higher resolution was hypothesized to be the primary control for maximum intensity of the model hurricane in Hausman (2001). From our default and 2x runs, MEAN(R) is about 3Δr, but MEAN(R) does not become systematically smaller at higher resolution, and the observed difference in steady-state intensity between the 2x and 4x runs cannot be explained solely by a change in RMW. Although the 4x run appears most intense in terms of central pressure, MEDIAN(P) seems also to have largely converged by the 2x run at about 915 mb. Both extreme updrafts, MAX(W), and the “smoothed” peak updrafts MAX[D(W)] appear not to have converged with these simulations.






5. The eye as a latent heat reservoir
a. Phenomenological evidence
The hypothesis of this paper, namely that a beneficial interaction between the eyewall and warm, moist air in the low-level eye permits modeled storm intensities that greatly exceed E-MPI, is motivated by the behavior of the default run before and after 20 days, when the transition from one quasi-steady intensity to another occurs (Fig. 3a). Prior to day 18 there is a θ anomaly (θ – θinitial) in the boundary layer of the eye (RMW at the surface ∼40 km), but almost no temperature anomaly (T – Tinitial; Fig. 4). After this time, warming is sufficient to increase the temperature near the surface as well. The higher-resolution simulations show this fully developed eye structure during the initial contraction of the vortex (before 10 days), and thus do not pass through multiple intensity states. Associated with this eye structure is an enhancement of Θe in low levels. In the case of the default run, the low-level eye Θe is a relative maximum near r = 20 km, z = 3 km on day 18 (Fig. 5). Because this enhanced Θe (Liu et al. 1999), if introduced into the convectively modified environment of the eyewall, can represent a source of local heat, this feature can be referred to as a reservoir of potential heat.
The establishment of the reservoir in our default run appears to occur in two steps. The first step involves the elimination of the initial midlevel Θe minimum. By day 6, the minimum has been pinched in half near r = 40 km with the development of persistent convection in the eyewall. This persistent convection equilibrates the eyewall close to a moist adiabat by day 8 (Fig. 5). After the convection is established, a secondary circulation is created in the eye that forces descent and recirculation into the eyewall at low levels [Gray (1998) notes the angular momentum impact of such an interaction]. This eye circulation advects the original Θe minimum to lower levels where it can come into contact with the top of the boundary layer and can be modified by the ocean. Surface fluxes then eradicate the Θe minimum, but by day 18 the vertical Θe gradient is reversed. The descent also advects down the high Θe (interchangeably θ with small qυ) that defines the tropopause. This warms the eye (Gray 1998).
At higher resolution, a multiscale response of the eye can now be resolved, with strong descent concentrated in a narrow column just interior to (i.e., at smaller radii than) the eyewall, with much weaker descent throughout the rest of the eye (e.g., Willoughby 1998). Because of this, the original midlevel Θe minimum is never totally eliminated from the eye (Fig. 6). On the other hand, Θe is locally enhanced just interior to the eyewall, which results from a combination of dry descent and subgrid-scale inward mixing (which is locally nonconservative) of vapor from the saturated eyewall. The channel of warmest Θe at day 10 is also under weak saturated ascent, interior to the peak updraft. By day 10, the 360-K contour connected to the upper troposphere descends down to 3 km. After this time this contour ascends and the lower-level reservoir appears to expand with the 360-K contour moving up to 5-km altitude. This occurs as a result of both an inward radial diffusion of vapor and vertical advection of heat and moisture from the eye boundary layer, which is enriched by interaction with the ocean. The default run (Δr = 15 km), which resolves the eye by one or two radial grid points, cannot represent this detailed circulation.
The ultimate source of enhanced entropy for the low-level eye in the 4x run is an upward flux of moisture from the ocean surface at greatly reduced surface pressure (Fig. 7), in contrast to the default case, which has a much weaker latent heat reservoir. The warming of the eye has led to a reversal of the vertical potential temperature gradient such that the eye is actually warming the ocean. However, subsaturation in the eye leads to a much greater flux of qυ into the atmosphere. The net effect is to positively force Θe. A negative flux removing moist entropy from the low-level eye is therefore necessary to maintain a steady state. Smith (2003) presents an idealized boundary layer model for hurricane-like vortices and finds the boundary layer remains subsaturated. For his cases, the surface latent heat flux is maximized at the RMW, but is also significantly positive beneath the eye, like that found in Fig. 7.
b. Lagrangian trajectories
To confirm whether the eyewall updraft is being affected by warm, moist air from the eye, forward/backward Lagrangian trajectories bisecting the eyewall updraft are computed and displayed in Figs. 8 and 9. These are “seeded” at day 24.4 of the 4x run across the updraft at the 2-km height at equal radial spacing and are calculated forward for 0.1 day and backward for 0.4 day. Figure 8 shows that there are three subjectively identified source regions for air upon introduction to the eyewall updraft: 1) the boundary layer inflow (solid), 2) low-level inflow above the boundary layer (dashed), and 3) the eye (dotted). Selection of a higher seeding altitude can also expose a midlevel (z ≈ 5 km) inflow (not shown). Of the inflow parcels shown here (Fig. 8; solid and dashed), some originate from a downdraft of a convective feature near r = 45 km and others from the boundary layer inflow from further radii. A similar feature was found by Braun (2002) using trajectories in a three-dimensional simulation of Hurricane Bob (1991; reproduced here as Fig. 10).
Figure 9 shows that the entropy of the downdraft trajectories is indistinguishable from the other inflow trajectories (solid and dashed) by the time they encounter the eyewall updraft. The general nature of these trajectories (Fig. 9) are quite similar to trajectories shown by Braun (2002) in Fig. 10c. The Θe of the inflow trajectories is mostly a function of their originating height, and as they enter the eyewall, lower trajectories warm even more while penetrating to smaller radii (similar to Braun 2002). This is entirely consistent with a dominant ocean interaction. A handful of boundary layer inflow trajectories obtain much greater Θe values (355 K or more; lower-right part of Fig. 9, dotted) as they enter the base of the eye; typically the radial inflow weakly extends to the center of the storm due to the Ekman pumping in the weak positive vorticity regime of the eye. These penetrating inflow trajectories linger up to 2 h in the eye before being introduced into the eyewall updraft. The left-to-right traces of trajectories in Fig. 9 (dotted) indicate mixing of entropy with neighboring air. Boundary layer inflow trajectories increase in entropy as they ascend in the eyewall updraft to z = 2.5 km. Since much of this increase occurs above the boundary layer, an exchange of entropy must result from an additional source of warm Θe other than the ocean, namely the eye.4 Above z = 2.5 km, trajectories ascending in the outer portion of the eyewall continue to increase in Θe, while those in the inner portion first decrease in Θe below z = 4.5 km (indicative of an exchange of entropy within the eyewall) and increase again above this level. Braun (2002) reached a similar conclusion with three-dimensional trajectories, suggesting that the results presented here are not peculiar to axisymmetric geometry. Parcels located near the inner edge of the eyewall detrain into the eye, where they increase in Θe, and frequently are reintroduced into the eyewall at a later time.
Based on the foregoing results, we can suggest that two related assumptions of E-MPI are violated in this model: the entropy exchange from the eye to the eyewall is nontrivial (which is confirmed with an Eulerian Θe budget that is not shown here) and the eyewall updraft is not moist neutral. The violation of moist neutrality is a consequence of a beneficial turbulent entrainment into the eyewall of high latent heat air from the low-level eye.
c. Experimental suppression of the reservoir
To confirm that the low-level entropy source in the eye is crucial for maintaining the hurricane at superintensity, we might propose to eliminate the source through some artificial means. An experiment was performed where enhanced cooling by Newtonian relaxation (our surrogate for radiation; 8 K day–1) was applied in the lower eye after day 15 of the 4x run. This leads to a rapid reduction of storm intensity from 90 to 75 m s–1 (not shown), but the intensity stabilized above E-MPI (55 m s–1). This occurred because the RMW expanded from 20 to 30 km, thus placing 10 km of radial separation between the eyewall and our imposed forcing (Fig. 11). This experimental run produced a radial displacement of the latent heat reservoir centered on r = 25 km (compared to the control experiment where the reservoir is nominally centered at the origin). While E-MPI does not explicitly depend on the RMW, the expansion of the RMW could in part explain the weakening. On the other hand, the displaced entropy reservoir also experiences an enhanced heat sink due to its proximity to the imposed forcing, which also might in part explain the weaker storm.
6. Heat or buoyancy?
What remains is to establish the means by which the introduction of air with high entropy from the eye to the eyewall produces a stronger storm. The candidate phenomena come in two classes: those that necessarily invoke the nonhydrostatic properties of the RE87 model and those phenomena that can be well described with hydrostatic principles.
It is known that nonhydrostatic processes in tornadoes can drive wind speeds to more than twice the limit suggested by thermodynamic forcing of buoyancy in the environment of the tornado (Fiedler and Rotunno 1986; observationally confirmed by Bluestein et al. 1993). Fielder (1994) shows that this can come about through either one of two phenomena: 1) an end-wall vortex, or 2) a drowned vortex jump. An end-wall vortex has its updraft directly at the axis of rotation that does not resemble our simulated hurricane. The drowned vortex jump has its peak updraft at a nonzero radius and bears superficial similarities to a hurricane vortex. The drowned vortex jump is out of cyclostrophic balance and results from “nearly loss-free ‘overshoot’ of the radial inflow into nearly the lowest pressure of the core [of the tornado-like vortex]” (Fielder 1994). The drowned vortex jump of Fielder (1994) shows a local pressure minimum near the radius of maximum low-level winds at the point where the inflow is stopped and directed upward that is different from Fig. 1b. Figure 1b suggests that the pressure perturbation is downward-directed and in part argues against a dynamic pressure force.








If we appeal to the conceptualization of E-MPI as a Carnot heat engine (Emanuel 1988), then the hurricane is a mechanism for producing work by processing heat. The measure of work performed ultimately is the dissipation of the system, and the rate of dissipation increases with the intensity of the storm. A stronger storm would thus require more heat addition to produce more work. The explanation of superintensity then simply comes down to added heat in the eyewall, which in part motivates the section title. Heat addition is of course not inconsistent with the generation of local buoyancy as discussed in Braun (2002) and the stretching of vortex tubes, but the distinction here is that, if the storm intensity is a function of the addition of heat (in the spirit of E-MPI theory) the role of local buoyancy is secondary or simply consistent.


7. Three-dimensional evidence for superintensity
Our explanation of superintensity in the RE87 model is that the eyewall is able to utilize of a second source of heat available from the eye. Although, in reality, E-MPI is rarely exceeded (Emanuel 2000), this mechanism may play a role in boosting the intensity of observed hurricanes that are subject to adverse conditions. Figure 13 shows results from an MM5 simulation of Hurricane Bob employing an unprecedented horizontal grid spacing of 1.3 km on the finest grid (Braun 2002; J. Fulton 2001). The mean RMW (dotted circle) at this time was 37 km and the mean radius of peak updraft (not shown) is 41 km. Equivalent potential temperature (exponential form) θe at the z = 1135 m level is enhanced in the eye (pink) relative to the eyewall (blue). Positive anomalies in θe at the edge of the θe enhancement (the interface between the eye and the eyewall) are on the southeast (outward advecting) side of positive invertible moist potential vorticity (PV) (IMPV; Schubert et al. 2001)7 anomalies. Fulton (2001) provides anecdotal evidence that, when these θe anomalies encounter the eyewall, they are associated with localized enhancements in updraft speed. Braun (2002) shows quantitatively that these vortical “hot towers” carry the majority of the updraft mass flux of the eyewall. Braun (2002) has further shown that these vortical hot towers are associated in this simulation of Hurricane Bob with the end points of a wavenumber-2 asymmetry that has an angular velocity that is well described by that of a barotropic vortex Rossby wave. There is an obvious suggestion then for a possible local generation of buoyancy, but Braun (2002) also shows that, in contrast to Zhang et al. (2000), a dynamic pressure force does not drive hot tower ascent above the approximate inflow boundary layer. Zhang et al. (2002) show supporting evidence of this in an MM5 simulation of Hurricane Andrew (1992) with a pronounced maximum of θe (to 391 K) in the low-level eye.
Frank and Ritchie (2001, their Fig. 16) present another MM5 simulation with peak tangential wind speeds greater than 100 m s–1 for a simulated hurricane over 28.5°C SST for a case of low vertical shear. E-MPI, assuming
Willoughby (1998) provides independent observational evidence that the low-level eye at times (Fig. 14) may have higher entropy air than the eyewall (by more than 10 K). From Willoughby (1998) we can also infer that aircraft flight-level data is not always indicative of the low-level eye entropy structure, depending on the altitude of the inversion level (Kossin and Eastin 2001). The boundary layer in our RE87 model is given 20 days to respond (Fig. 15) to the intensity and surface fluxes of the hurricane, while real hurricanes are hardly allowed to attain a steady state for a period of many days. To illustrate this point the model eye sounding does not show a sharp inversion and Θe is minimized at midlevels (400 mb) rather than remaining constant with height as observations show (e.g., Willoughby 1998). After constructing a composite of dropsonde data obtained from the National Oceanic and Atmospheric Administration (NOAA) P-3 aircraft, LeeJoice (2000) finds that eye soundings at 0 ≤ z < 3000 m average 5 to 10 K warmer (in terms of θe) than in the eyewall. This attribute is confirmed in the RE87 model in Fig. 15. LeeJoice (2000) also finds that dropsondes show little vertical gradient in terms of θe in the eyewall, but a significant gradient (10 K over 3000 m) in the eye, generally supporting our interpretation of Willoughby (1998) and the observed moist entropy structure in our simulations.
8. Discussion
We have demonstrated that the RE87 model will produce a model storm of increasing intensity with increasing resolution. Convergence in the maximum tangential winds was found with sufficiently high resolution. Hausman (2001) suggested that the increase in storm intensity is due to decreases in the RMW at higher resolution; smaller RMWs simply were not resolvable at coarse resolution. Our finding with the RE87 model is that the RMW stabilizes by the 2x run, but increases in intensity further occur to the 4x run (Table 2). Because E-MPI has been presented previously as a rigorous upper bound on hurricane intensity, an investigation into its theoretical bases was warranted. E-MPI theory provides the context for extracting and summarizing the key phenomena for superintensity.
Figure 16 summarizes many of the important interactions associated with the superintensity phenomenon, where the abscissa is displayed as a function of absolute angular momentum (M = rυ + fr2/2). We have demonstrated a low-level entropy maximum in the eye (Fig. 16, Θe contours in red) that is associated with superintensity relative to E-MPI and an interaction between the warm eye and the eyewall (Fig. 16, blue streamlines near M = 1 × 106 m2 s–1). After showing that the storm evolves largely in gradient and hydrostatic balance, the related assumptions of E-MPI theory that the eyewall is in slantwise, moist adiabatic ascent without interaction with the eye are violated (Fig. 16, noting from the vertical alignment of streamlines in the green shaded, saturated eyewall region of the storm that angular momentum is better conserved than Θe). Through a simple addition of heat much of the superintensity of the modeled storm can be explained with an ad hoc modification of E-MPI.
The eddy features in the streamlines of Fig. 16 can be recognized as the breakdown of an azimuthal vortex sheet and as a primary advecting mechanism for high-entropy air between the eye and the eyewall. With a vertical scale of only about 1 km, these features are not resolvable at low vertical resolution. Associated with these eddy features is a tangential vorticity of ≈–1 × 10–2 s–1 (not shown). A low-level outflow jet is identifiable, perhaps, below z ;lt 1 km; however, the eddies extending up to z = 5 km are distinct from the low-level jet.
The internal redistribution of heat within the eyewall by diffusion raises the question whether the basic result presented here is a peculiarity of the diffusion scheme of the model. An experiment at 4x resolution using a greatly reduced diffusion (by a factor of 16 from the 4x run, mixing length lh by a factor of 4; see Table 3) suggests that the role of diffusion is replaced by temporal eddies, and the time series of Vmax (not shown) exhibits a shorter period variability of ≈±4 m s–1 from the long-term mean, but still is superintense (Vmax ≈ 85 m s–1). Although the energetic inconsistency of the RE87 model is a potential issue (K. Emanuel 2002, personal communication), the fact that superintensity persists as the mixing length is decreased substantially suggests that the subgrid-scale effects are not the reason for the superintensity phenomenon. Moreover, the availability of at least one source of high entropy (the ocean beneath the eye) and the resolved-scale breakdown of a transverse vortex sheet argue that some aspect of this superintensity mechanism should still operate in a variety of real situations.
A mathematical proof is not known to the authors of whether a nonlinear dynamical system prefers to exist in a state that maximizes the rate of entropy production, but Sawada (1981) provides a compelling argument for such a system bounded by two infinite thermal reservoirs. Under different conditions, different dynamical pathways are simply not available, which was used to describe, for example, different crystallization habits or convection modes between plates as states where pathways become available as geometric or thermodynamic constraints are relaxed. Different physical systems have been demonstrated to exist in or relax toward states of maximum entropy or maximum rate of entropy production (e.g., Malkus and Veronis 1958; Malkus and Riehl 1960; Stephens and O’Brien 1993; O’Brien 1997; Ozawa and Ohmura 1997; Jin and Dubin 1998; Schubert et al. 1999; Schecter et al. 1999; Paltridge 2001; Lorenz et al. 2001). The superintense behavior of the RE87 model at high resolutions may simply be the heat engine existing in a more effective dynamical pathway. What remains as future work is to explain why superintense behavior of hurricane models is a preferred pathway and to discover the constraints on accessing this pathway.
Sensitivity studies with the horizontal subgrid-scale parameterization in the RE87 model (not shown) indicate that a more diffusive simulation will produce a weaker (yet still superintense) storm. Since the eyewall warming (ΔΘe) observed in these experiments is likewise reduced, we believe the picture presented here is still consistent. Diffusion in the axisymmetric model not only models the effect of subgrid-scale axisymmetric mixing, but more prominently serves as a crude surrogate for the impact of otherwise resolvable asymmetries. Barotropic instability (Michalke and Timme 1967; Schubert et al. 1999) naturally results from a ring of convection (an ideal eyewall) due to vortex tube stretching, which would lead to the growth of vortex Rossby waves and eventually nonlinear, coherent vortex structures (eyewall mesovortices; Kossin and Schubert 2001; Montgomery et al. 2002). The radial structure of axisymmetric mean vorticity in a steady-state three-dimensional storm should be an approximate balance between sources due to convection and losses due to redistribution and friction. In a similar manner, the low-level entropy reservoir is a balance between an ocean source and loss to the eyewall (and to a lesser extent radiation). In a steady-state view, a slower mixing rate (less vigorous vortex Rossby waves and mesovortices) allows the entropy reservoir to build to greater levels and potentially allow for more eyewall warming (ΔΘe). In a more realistic, time-varying hurricane, the existence of more vigorous mixing may promote short-term intensification of the vortex and a temporary exhaustion of the heat reservoir. What remains as a basic question for future research is the role of vortex Rossby waves and eyewall mesovortices (Montgomery et al. 2002) in mixing thermodynamic properties and their impact on MPI.
Many mechanisms have been proposed that would provide a storm intensity that would be weaker than E-MPI. Adverse vertical wind shear is believed to disrupt the thermal maximum of the upper-level eye, providing for a weaker central pressure (Gray 1968; DeMaria 1996; Frank and Ritchie 1999, 2001). Secondary eyewalls are found in very intense hurricanes (Willoughby et al. 1982) that interrupt the radial inflow of environmental angular momentum (Camp and Montgomery 2001); intensity fluctuations are found with eyewall replacement cycles. Hausman (2001) found that axisymmetric hurricane simulations using ice microphysics are weaker than comparable rain-only simulations. DeCosmo et al. (1996) suggest that the ratio of Ck/CD becomes less than unity at high wind speeds, although Andreas and Emanuel (2001) present an argument that an accounting of the entropy flux due to reentrant sea spray could lead to a positive flux of entropy that nearly cancels the anticipated reduction of Vmax from changes in Ck/CD. Drag with the ocean surface generates mixing currents in the ocean that bring cool, deep ocean water to the surface, thus decreasing the effective surface temperature that the hurricane is in contact with (Shay et al. 1998; Bender and Ginis 2000; Jacob et al. 2000). In the face of these negative influences on intensity, the hurricane may take advantage of dissipational heating of the boundary layer (Bister and Emanuel 1998; Businger and Businger 2001) as a secondary heat source for the storm. Bister and Emanuel (1998) ran a version of the RE87 model that includes dissipational heating as a forcing term. This showed a 13 m s–1 increase in intensity with the inclusion of dissipation heating (Table 3). This compares to the 35 m s–1 increase associated with the eye entropy reservoir demonstrated here. Thus, eyewall heating may be as important or more important to storm intensity as dissipational heating.
9. Conclusions
While E-MPI has been previously presented as a rigorous upper bound on storm intensity, at high spatial and temporal resolution the RE87 hurricane model produces storms that greatly exceed E-MPI. Because the superintense state is well captured by Eliassen’s (1951) balanced vortex model, nonhydrostatic processes on the vortex scale are not viable candidates for explaining superintensity. The cause of superintensity is found to lie in the presence of high-entropy air in the low-level eye, which is entrained into the eyewall. The introduction of high-entropy air from the eye represents an additional source of heat to the eyewall of the storm and leads to a modified Carnot cycle that supports a stronger storm. Observational evidence points to the existence of this mechanism in real hurricanes and also in three-dimensional numerical models.
It is proposed that a new MPI should be formulated from first principles that incorporates this effect, whether as a novel approach to the problem or as a modification of E-MPI. It is proposed that further analyses of the thermodynamic structure of both idealized and real-case simulations of hurricanes be conducted in three dimensions. An interesting outcome of this work that bears further scrutiny is that the axisymmetric hypercane regime described by Emanuel (1988) may become more accessible with small perturbations of presently observed thermodynamic characteristics of the atmosphere (Knutson and Tuleya 1999). Whether hypercanes are indeed even possible in a three-dimensional context is an open question that will be the consideration of upcoming work.
Acknowledgments
This work was supported in part by Office of Naval Research Grants N00014-93-1-0456 P00006, N00014-02-1-0474, National Science Foundation ATM-0101781, National Aeronautics and Space Administration (NASA) Fourth Convective and Moisture Experiment (CAMEX) Grant NAG5-11140, and Colorado State University. We would like to thank R. Rotunno for kindly making his model available to us and for valuable discussion that helped clarify thought and presentation. Thanks also to H. Willoughby at National Oceanic and Atmospheric Administration Atlantic Oceanic and Meteorological Laboratory’s Hurricane Research Division and S. Braun of NASA Goddard Space Flight Center for the use of their figures. We would also like to thank W. Gray, W. Schubert, K. Emanuel, M. Eastin, S. Grandlienard, T. Cram, S. Hausman, M. Kirby, G. Stephens, and D. O’Brien, and reviewers J. Kossin and S. Gray for their helpful suggestions and input.
REFERENCES
AMS Council, 2000: Policy statement: Hurricane research and forecasting. Bull. Amer. Meteor. Soc., 81 , 1341–1346.
Andreas, E. L., and K. A. Emanuel, 2001: Effects of sea spray on tropical cyclone intensity. J. Atmos. Sci., 58 , 3741–3751.
Bender, M. A., and I. Ginis, 2000: Real-case simulations of hurricane–ocean interaction using a high-resolution coupled model: Effects on hurricane intensity. Mon. Wea. Rev., 128 , 917–946.
Bengtsson, L., 2001: Hurricane threats. Science, 293 , 440–441.
Bister, M., and K. A. Emanuel, 1998: Dissipative heating and hurricane intensity. Meteor. Atmos. Phys., 65 , 233–240.
Bluestein, H. B., J. G. Ladue, H. Stein, D. Speheger, and W. P. Unruh, 1993: Doppler radar wind spectra of supercell tornadoes. Mon. Wea. Rev., 121 , 2200–2221.
Bolton, D., 1980: The computation of equivalent potential temperature. Mon. Wea. Rev., 108 , 1046–1053.
Bosart, L. F., C. S. Velden, W. E. Bracken, J. Molinari, and P. G. Black, 2000: Environmental influences on the rapid intensification of Hurricane Opal (1995) over the Gulf of Mexico. Mon. Wea. Rev., 128 , 322–352.
Braun, S. A., 2002: A cloud-resolving simulation of Hurricane Bob (1991): Storm structure and eyewall buoyancy. Mon. Wea. Rev., 130 , 1573–1592.
Businger, S., and J. A. Businger, 2001: Viscous dissipation of turbulent kinetic energy in storms. J. Atmos. Sci., 58 , 3793–3796.
Camp, J. P., 1999: Hurricane maximum intensity: Past and present. M.S. thesis, Colorado State University, 147 pp.
Camp, J. P., and M. T. Montgomery, 2001: Hurricane maximum intensity: Past and present. Mon. Wea. Rev., 129 , 1704–1717.
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.
DeCosmo, J., K. B. Katsaros, S. D. Smith, R. J. Anderson, W. A. Oost, K. Bumke, and H. Chadwick, 1996: Air–sea exchange of water vapor and sensible heat: The Humidity Exchange Over the Sea (HEXOS) results. J. Geophys. Res., 101 , 12001–12016.
DeMaria, M., 1996: The effect of vertical shear on tropical cyclone intensity change. J. Atmos. Sci., 53 , 2076–2087.
Eliassen, A., 1951: Slow thermally or frictionally controlled meridional circulation in a circular vortex. Astrophys. Norv., 5 , 19–60.
Emanuel, K. A., 1986: An air–sea interaction theory for tropical cyclones. Part I: Steady-state maintenance. J. Atmos. Sci., 43 , 585–604.
Emanuel, K. A., 1988: The maximum intensity of hurricanes. J. Atmos. Sci., 45 , 1143–1155.
Emanuel, K. A., 1989: The finite-amplitude nature of tropical cyclogenesis. J. Atmos. Sci., 46 , 3431–3456.
Emanuel, K. A., 1995a: The behavior of a simple hurricane model using a convective scheme based on subcloud-layer entropy equilibrium. J. Atmos. Sci., 52 , 3960–3968.
Emanuel, K. A., 1995b: Sensitivity of tropical cyclones to surface exchange coefficients and a revised steady-state model incorporating eye dynamics. J. Atmos. Sci., 52 , 3969–3976.
Emanuel, K. A., 2000: A statistical analysis of tropical cyclone intensity. Mon. Wea. Rev., 128 , 1139–1152.
Fairall, C. W., J. D. Kepert, and G. J. Holland, 1994: The effect of sea spray on surface energy transports over the ocean. Global Atmos. Ocean Syst., 2 , 121–142.
Fiedler, B. H., 1994: The thermodynamic speed limit and its violation in axisymmetric numerical simulations of tornado-like vortices. Atmos.–Ocean, 32 , 335–359.
Fiedler, B. H., and R. Rotunno, 1986: A theory for the maximum windspeeds in tornado-like vortices. J. Atmos. Sci., 43 , 2328–2340.
Frank, W. M., and E. A. Ritchie, 1999: Effects of environmental flow upon tropical cyclone structure. Mon. Wea. Rev., 127 , 2044–2061.
Frank, W. M., and E. A. Ritchie, 2001: Effects of vertical wind shear on the intensity and structure of numerically simulation hurricanes. Mon. Wea. Rev., 129 , 2249–2269.
Fulton, J. D., 2001: Insights into the hurricane vortex using a high resolution numerical simulation of Hurricane Bob (1991). M.S. thesis, Colorado State University, 84 pp.
Gray, S. L., 1998: Analysis of the eyes formed in simulated tropical cyclones and polar lows. Quart. J. Roy. Meteor. Soc., 124 , 2357–2375.
Gray, W. M., 1968: Global view of tropical disturbances and storms. Mon. Wea. Rev., 96 , 669–700.
Guinn, T. A., and W. H. Schubert, 1993: Hurricane spiral bands. J. Atmos. Sci., 50 , 3380–3403.
Halliday, D., and R. Resnick, 1988: Fundamentals of Physics, 3d ed. John Wiley and Sons, 977 pp.
Hausman, S. A., 2001: Formulation and sensitivity analysis of a nonhydrostatic, axisymmetric tropical cyclone model. Atmospheric Science Paper 701, Dept. of Atmospheric Science, Colorado State University, 210 pp.
Heisler, H., 1995: Advanced Engine Technology. Arnold, 794 pp.
Holton, J. R., 1992: An Introduction to Dynamic Meteorology. 3d ed. Academic Press, 511 pp.
Hoskins, B. J., and F. P. Bretherton, 1972: Atmospheric frontogenesis models: Mathematical formulation and solution. J. Atmos. Sci., 29 , 11–37.
Jacob, S. D., L. K. Shay, A. J. Mariano, and P. G. Black, 2000: The 3D oceanic mixed layer response to Hurricane Gilbert. J. Phys. Oceanogr., 30 , 1407–1429.
Jin, D. Z., and D. H. E. Dubin, 1998: Regional maximum entropy theory for vortex crystal formation. Phys. Rev. Lett., 80 , 4434–4437.
Jones, S. C., 1995: The evolution of vortices in vertical shear. I: Initially barotropic vortices. Quart. J. Roy. Meteor. Soc., 121 , 821–851.
Jordan, C. L., 1958: Mean soundings for the West Indies area. J. Meteor., 15 , 91–97.
Kleinschmidt, E., 1951: Grundlagen einer Theorie der Tropischen Zyklonen (Basic principles for a theory of tropical cyclones). Arch. Meteor. Geophys. Bioklimatol., A4. 53–72.
Klemp, J. B., and R. B. Wilhelmson, 1978: The simulation of three-dimensional convective storm dynamics. J. Atmos. Sci., 35 , 1070–1096.
Knutson, T. R., and R. E. Tuleya, 1999: Increase hurricane intensities with CO2-induced warming as simulated using the GFDL hurricane prediction system. Climate Dyn., 15 , 503–519.
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.
LeeJoice, R. N., 2000: Hurricane inner-core structure as revealed by GPS dropwindsondes. M.S. thesis, Colorado State University, 56 pp.
Lilly, D. K., 1962: On the numerical simulation of buoyant convection. Tellus, 14 , 148–172.
Liu, Y., D-L. Zhang, and M. K. Yau, 1999: A multiscale numerical study of Hurricane Andrew (1992). Part II: Kinematics and inner-core structures. Mon. Wea. Rev., 127 , 2597–2616.
Lorenz, R. D., J. I. Lunine, P. G. Withers, and C. P. McKay, 2001: Titan, Mars and Earth: Entropy production by latitudinal heat transport. Geophys. Res. Lett., 28 , 415–418.
Malkus, J. S., and H. Riehl, 1960: On the dynamics and energy transformations in steady-state hurricanes. Tellus, 12 , 1–20.
Malkus, W. V. R., and G. Veronis, 1958: Finite amplitude cellular convection. J. Fluid Mech., 4 , 225–260.
McIntyre, M. E., 1992: Isentropic distributions of potential vorticity and their relevance to tropical cyclone dynamics. Tropical Cyclone Disasters: Proceedings of ISCU/WMO International Symposium, J. Lighthill et al., Eds., Peking University Press, 143–156.
Michalke, A., and A. Timme, 1967: On the inviscid instability of certain two-dimensional vortex-type flows. J. Fluid Mech., 29 , 647–666.
Möller, J. D., and M. T. Montgomery, 2000: Tropical cyclone evolution via potential vorticity anomalies in a three-dimensional balance model. J. Atmos. Sci., 57 , 3366–3387.
Möller, J. D., and L. J. Shapiro, 2002: Balance contributions to the intensification of Hurricane Opal as diagnosed from a GFDL model forecast. Mon. Wea. Rev., 130 , 1866–1881.
Montgomery, M. T., 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.
Montgomery, M. T., V. A. Vladimirov, and P. V. Denissenko, 2002: An experimental study on hurricane mesovortices. J. Fluid Mech., 471 , 1–32.
O’Brien, D. M., 1997: A yardstick for global entropy flux. Quart. J. Roy. Meteor. Soc., 123 , 243–260.
Ooyama, K. V., 2001: A dynamic and thermodynamic foundation for modeling the moist atmosphere with parameterized microphysics. J. Atmos. Sci., 58 , 2073–2102.
Ozawa, A., and H. Ohmura, 1997: Thermodynamics of a global-mean state of the atmosphere: A state of maximum entropy increase. J. Climate, 10 , 441–445.
Paltridge, G. W., 2001: A physical basis for a maximum thermodynamic dissipation of the climate system. Quart. J. Roy. Meteor. Soc., 127 , 305–313.
Persing, J., M. T. Montgomery, and R. E. Tuleya, 2002: Environmental interactions in the GFDL hurricane model for Hurricane Opal. Mon. Wea. Rev., 130 , 298–317.
Roll, H. V., 1965: Physics of the Marine Atmosphere. Academic Press, 426 pp.
Rotunno, R., and K. A. Emanuel, 1987: An air–sea interaction theory for tropical cyclones. Part II: Evolutionary study using a nonhydrostatic axisymmetric numerical model. J. Atmos. Sci., 44 , 542–561.
Sawada, Y., 1981: A thermodynamic variational principal in nonlinear non-equilibrium phenomena. Prog. Theor. Phys., 66 , 68–76.
Schecter, D. A., D. H. E. Dubin, K. S. Fine, and C. F. Driscoll, 1999: Vortex crystals from 2D Euler flow: Experiment and simulation. Phys. Fluids, 11 , 905–914.
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.
Schubert, W. H., S. A. Hausman, M. Garcia, K. V. Ooyama, and H-C. Kuo, 2001: Potential vorticity in a moist atmosphere. J. Atmos. Sci., 58 , 3148–3157.
Shapiro, L. J., 1983: The asymmetric boundary layer flow under a translating hurricane. J. Atmos. Sci., 40 , 1984–1998.
Shay, L. K., A. J. Mariano, S. D. Jacob, and E. H. Ryan, 1998: Mean and near-inertial ocean current response to Hurricane Gilbert. J. Phys. Oceanogr., 28 , 858–889.
Smagorinsky, J., 1963: General circulation experiments with the primitive equation. I: The basic experiment. Mon. Wea. Rev., 91 , 99–164.
Smith, R. K., 2003: A simple model for the hurricane boundary layer. Quart. J. Roy. Meteor. Soc., 129 , 1007–1027.
Stephens, G. L., and D. M. O’Brien, 1993: Entropy and climate. I: ERBE observations of the entropy production of the earth. J. Roy. Meteor. Soc., 119 , 121–152.
Wang, Y., 2002a: Vortex Rossby waves in a numerically simulated tropical cyclone. Part I: Overall structure, potential vorticity, and kinetic energy budgets. J. Atmos. Sci., 59 , 1213–1238.
Wang, Y., 2002b: Vortex Rossby waves in a numerically simulated tropical cyclone. Part II: The role in tropical cyclone structure and intensity changes. J. Atmos. Sci., 59 , 1239–1262.
Wang, Y., and G. J. Holland, 1996: The beta drift of baroclinic vortices. Part I: Adiabatic vortices. J. Atmos. Sci., 53 , 411–427.
Willoughby, H. E., 1990: Gradient balance in tropical cyclones. J. Atmos. Sci., 47 , 265–274.
Willoughby, H. E., 1998: Tropical cyclone eye thermodynamics. Mon Wea. Rev., 126 , 1653–1680.
Willoughby, H. E., J. A. Clos, and M. G. Shoreibah, 1982: Concentric eye walls, secondary wind maxima, and the evolution of the hurricane vortex. J. Atmos. Sci., 39 , 395–411.
Yamasaki, M., 1983: A further study of the tropical cyclone without parameterizing the effects of cumulus convection. Pap. Meteor. Geophys., 34 , 221–260.
Zhang, D-L., Y. Liu, and M. K. Yau, 2000: A multiscale numerical study of Hurricane Andrew (1992). Part III: Dynamically induced vertical motion. Mon. Wea. Rev., 128 , 3772–3788.
Zhang, D-L., Y. Liu, and M. K. Yau, 2002: A multiscale numerical study of Hurricane Andrew (1992). Part V: Inner-core thermodynamics. Mon. Wea. Rev., 130 , 2745–2763.
APPENDIX
An Assessment of the Assumptions and Approximations of E-MPI
This appendix discusses some sensitivities and approximations in E-MPI that can potentially impact the estimate of intensity. We find that all of the candidate sensitivities are too small to explain the superintensity shown in high-resolution simulations of the RE87 model (see Fig. 3).
Balance of entropy and angular momentum


Relative humidity




Another difficulty with applying
Boundary layer height
According to Emanuel (1995a), the top of the boundary layer is the proper reference for evaluating χb. In this paper, when χb has been evaluated, it has been at the first grid level, the level of saturation at the eyewall. At low resolution (the default run), this level (z = 625 m) is also saturated in the near environment of the eyewall extending to r ∼ 80 km, except where downdrafts are occurring. As resolution is increased, the level of saturation in the near environment is generally found at z = 400 m. Selecting this level in computation of χb leads to only a small increase in estimating Vmax using (6) at higher resolutions (4x and 8x runs) of 5 to 10 m s–1, which is insufficient to explain the observed degree of superintensity.
Gradient wind balance




Thermal wind balance


We can relax the moist neutral assumption in one specific way such that (A6) would still be valid (subject to other assumptions of the derivation like gradient wind balance; hydrostatic balance; and persistent, slantwise convection, for example). Specifically we can relax the assumption that ∂s*/∂M)p is a constant along M surfaces, which would permit some moist static stability so long as that stabilization has more or less the same vertical structure on different M surfaces.










Thus, although the approximations used in deriving (A6) can be somewhat relaxed using values found in our model runs, they are incapable of explaining superintensity.
Outflow temperature
From (A8), it is clear that the meaning of Tout should be at the end point of the integral defined above. This would appear to permit different values of Tout as a function of radius. Also Tout would be defined at a point, rather than as an average as defined by (10), and would be generally colder than (10). There is no reason to anticipate that this definition of Tout would be greatly different than that defined by (11), used elsewhere in this paper. In any case, no reasonable value of Tout supplied into E-MPI theory can explain the degree of superintensity observed in our model runs.

The 20.0–21.0-day average fields in the 4x run near the eyewall. (a)–(e) Shaded regions indicate negative values of the field: (a) tangential velocity (contour interval is 5 m s–1); (b) dimensional pressure deviation from the outer boundary (at 1500-km radius; contour interval is 5 mb); (c) radial velocity (contour interval is 4 m s–1); (d) potential temperature deviation from the initial state (contour interval is 2 K); (e) vertical velocity (contour interval is 0.5 m s–1); (f) liquid water mixing ratio (contour interval is 1 g kg–1); dashed line denote the 0.1 g kg–1 contour and water vapor mixing ratio; lightest shading indicates 3 < qυ < 8 g kg–1; darker shading indicates 8 < qυ < 13 g kg –1; darkest shading indicates qυ > 13 g kg–1
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

The 20.0–21.0-day average fields in the 4x run near the eyewall. (a)–(e) Shaded regions indicate negative values of the field: (a) tangential velocity (contour interval is 5 m s–1); (b) dimensional pressure deviation from the outer boundary (at 1500-km radius; contour interval is 5 mb); (c) radial velocity (contour interval is 4 m s–1); (d) potential temperature deviation from the initial state (contour interval is 2 K); (e) vertical velocity (contour interval is 0.5 m s–1); (f) liquid water mixing ratio (contour interval is 1 g kg–1); dashed line denote the 0.1 g kg–1 contour and water vapor mixing ratio; lightest shading indicates 3 < qυ < 8 g kg–1; darker shading indicates 8 < qυ < 13 g kg –1; darkest shading indicates qυ > 13 g kg–1
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
The 20.0–21.0-day average fields in the 4x run near the eyewall. (a)–(e) Shaded regions indicate negative values of the field: (a) tangential velocity (contour interval is 5 m s–1); (b) dimensional pressure deviation from the outer boundary (at 1500-km radius; contour interval is 5 mb); (c) radial velocity (contour interval is 4 m s–1); (d) potential temperature deviation from the initial state (contour interval is 2 K); (e) vertical velocity (contour interval is 0.5 m s–1); (f) liquid water mixing ratio (contour interval is 1 g kg–1); dashed line denote the 0.1 g kg–1 contour and water vapor mixing ratio; lightest shading indicates 3 < qυ < 8 g kg–1; darker shading indicates 8 < qυ < 13 g kg –1; darkest shading indicates qυ > 13 g kg–1
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Vmax from E-MPI as a function of Ts and Tout. The boxed area denotes the range of values found in the present paper.
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Vmax from E-MPI as a function of Ts and Tout. The boxed area denotes the range of values found in the present paper.
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Vmax from E-MPI as a function of Ts and Tout. The boxed area denotes the range of values found in the present paper.
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Time series of maximum tangential winds for the (a) default, (b) 2x, (c) 4x, and (d) 8x runs. Overlaid are theoretical predictions of maximum intensity using E-MPI (dotted) from Eq. (2) (with a fixed 80% RH) and using a boundary layer balance (dashed) from Eq. (6) evaluated at the RMW. The overlays are computed once daily using daily averaged model output
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Time series of maximum tangential winds for the (a) default, (b) 2x, (c) 4x, and (d) 8x runs. Overlaid are theoretical predictions of maximum intensity using E-MPI (dotted) from Eq. (2) (with a fixed 80% RH) and using a boundary layer balance (dashed) from Eq. (6) evaluated at the RMW. The overlays are computed once daily using daily averaged model output
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Time series of maximum tangential winds for the (a) default, (b) 2x, (c) 4x, and (d) 8x runs. Overlaid are theoretical predictions of maximum intensity using E-MPI (dotted) from Eq. (2) (with a fixed 80% RH) and using a boundary layer balance (dashed) from Eq. (6) evaluated at the RMW. The overlays are computed once daily using daily averaged model output
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Temperature T and potential temperature θ (K) from the default run displayed as departures from the initial state for days 12 and 28
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Temperature T and potential temperature θ (K) from the default run displayed as departures from the initial state for days 12 and 28
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Temperature T and potential temperature θ (K) from the default run displayed as departures from the initial state for days 12 and 28
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Equivalent potential temperature Θe(K) at various times in the default run
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Equivalent potential temperature Θe(K) at various times in the default run
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Equivalent potential temperature Θe(K) at various times in the default run
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Equivalent potential temperature Θe (K) at various times in the 4x run
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Equivalent potential temperature Θe (K) at various times in the 4x run
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Equivalent potential temperature Θe (K) at various times in the 4x run
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Surface forcing terms from the 4x run
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Surface forcing terms from the 4x run
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Surface forcing terms from the 4x run
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Half-day trajectories (when contained within the display domain) from the 4x run. These are forward/back trajectories seeded at 24.40 days at a height of 2 km, bisecting the eyewall draft at equal spatial intervals. Trajectories are computed forward for 0.1 days and backward 0.4 days using the same seed points. Dotted trajectories enter the eyewall from the eye; solid trajectories enter the eyewall from the boundary layer inflow; dashed trajectories enter the eyewall from above the boundary layer; Xs mark the seeding locations
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Half-day trajectories (when contained within the display domain) from the 4x run. These are forward/back trajectories seeded at 24.40 days at a height of 2 km, bisecting the eyewall draft at equal spatial intervals. Trajectories are computed forward for 0.1 days and backward 0.4 days using the same seed points. Dotted trajectories enter the eyewall from the eye; solid trajectories enter the eyewall from the boundary layer inflow; dashed trajectories enter the eyewall from above the boundary layer; Xs mark the seeding locations
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Half-day trajectories (when contained within the display domain) from the 4x run. These are forward/back trajectories seeded at 24.40 days at a height of 2 km, bisecting the eyewall draft at equal spatial intervals. Trajectories are computed forward for 0.1 days and backward 0.4 days using the same seed points. Dotted trajectories enter the eyewall from the eye; solid trajectories enter the eyewall from the boundary layer inflow; dashed trajectories enter the eyewall from above the boundary layer; Xs mark the seeding locations
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

The same trajectories as Fig. 8, displayed as a function of Θe and height. Dotted trajectories enter the eyewall from the eye; solid trajectories enter the eyewall from the boundary layer inflow dashed trajectories enter the eyewall from above the boundary layer; Xs mark the seeding locations
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

The same trajectories as Fig. 8, displayed as a function of Θe and height. Dotted trajectories enter the eyewall from the eye; solid trajectories enter the eyewall from the boundary layer inflow dashed trajectories enter the eyewall from above the boundary layer; Xs mark the seeding locations
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
The same trajectories as Fig. 8, displayed as a function of Θe and height. Dotted trajectories enter the eyewall from the eye; solid trajectories enter the eyewall from the boundary layer inflow dashed trajectories enter the eyewall from above the boundary layer; Xs mark the seeding locations
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Figure 13 of Braun (2002), displaying trajectories from a 3D simulation of Hurricane Bob using the MM5 model. The trajectories are seeded at a point of locally enhanced ascent within the eyewall. (a) Horizontal projection of the trajectory paths in a reference frame moving with the storm. Shading indicates simulated reflectivities greater than 40 dBZ at 66 h from initialization. Filled circles and arrowheads indicate the beginnings and ends, respectively, of the trajectories. The small rectangular box indicates where the trajectories were seeded. (b) Profiles of trajectory radius vs height. (c) Profiles of trajectory θe vs height. In (b) and (c), small arrows indicate the general direction of movement of the air parcels
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Figure 13 of Braun (2002), displaying trajectories from a 3D simulation of Hurricane Bob using the MM5 model. The trajectories are seeded at a point of locally enhanced ascent within the eyewall. (a) Horizontal projection of the trajectory paths in a reference frame moving with the storm. Shading indicates simulated reflectivities greater than 40 dBZ at 66 h from initialization. Filled circles and arrowheads indicate the beginnings and ends, respectively, of the trajectories. The small rectangular box indicates where the trajectories were seeded. (b) Profiles of trajectory radius vs height. (c) Profiles of trajectory θe vs height. In (b) and (c), small arrows indicate the general direction of movement of the air parcels
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Figure 13 of Braun (2002), displaying trajectories from a 3D simulation of Hurricane Bob using the MM5 model. The trajectories are seeded at a point of locally enhanced ascent within the eyewall. (a) Horizontal projection of the trajectory paths in a reference frame moving with the storm. Shading indicates simulated reflectivities greater than 40 dBZ at 66 h from initialization. Filled circles and arrowheads indicate the beginnings and ends, respectively, of the trajectories. The small rectangular box indicates where the trajectories were seeded. (b) Profiles of trajectory radius vs height. (c) Profiles of trajectory θe vs height. In (b) and (c), small arrows indicate the general direction of movement of the air parcels
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Cross section of Θe (K) from the 4x run before and after imposition of an artificial heat sink in the lower eye
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Cross section of Θe (K) from the 4x run before and after imposition of an artificial heat sink in the lower eye
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Cross section of Θe (K) from the 4x run before and after imposition of an artificial heat sink in the lower eye
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

The daily averaged secondary circulation observed from the 4x run at day 25, and the same inferred from the Eliassen (1951) balanced vortex model using the observed axisymmetric structure and forcings at the same time
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

The daily averaged secondary circulation observed from the 4x run at day 25, and the same inferred from the Eliassen (1951) balanced vortex model using the observed axisymmetric structure and forcings at the same time
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
The daily averaged secondary circulation observed from the 4x run at day 25, and the same inferred from the Eliassen (1951) balanced vortex model using the observed axisymmetric structure and forcings at the same time
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

IMPV at z = 1135 m at (a) 268, (b) 272, and (c) 276 min into an MM5 simulation of Hurricane Bob, kindly provided by S. A. Braun, NASA/GSFC (see Braun 2002 for details). Data is reported at the 1.3-km grid spacing of the innermost mesh. (d)–(f) Equivalent potential temperature θe displayed at respective times. The hurricane symbol denotes the storm center. Levels of shading in (a)–(c) denote no shading: IMPV < 0 PVUs; grey shading: 0 < IMPV < 15 PVUs; blue shading: 15 < IMPV < 30 PVUs; pink shading: 30 < IMPV < 60 PVUs; red shading: IMPV > 60 PVUs. Levels of shading in (d)–(f) denote no shading: θe < 345 K; grey shading: 345 < θe < 351 K; blue shading: 351 < θe < 357 K; pink shading: 357 < θe < 363 K; red shading: θe > 363 K. The dotted circle is the radius of maximum winds, r = 37 km. The “+” is to mark the position of a typical enhancement of IMPV interior to the RMW. Cyclonic motion around this IMPV enhancement appears to position high-θe air radially outward to the RMW around the southeast side of the feature. Figure adapted from Fulton (2001)
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

IMPV at z = 1135 m at (a) 268, (b) 272, and (c) 276 min into an MM5 simulation of Hurricane Bob, kindly provided by S. A. Braun, NASA/GSFC (see Braun 2002 for details). Data is reported at the 1.3-km grid spacing of the innermost mesh. (d)–(f) Equivalent potential temperature θe displayed at respective times. The hurricane symbol denotes the storm center. Levels of shading in (a)–(c) denote no shading: IMPV < 0 PVUs; grey shading: 0 < IMPV < 15 PVUs; blue shading: 15 < IMPV < 30 PVUs; pink shading: 30 < IMPV < 60 PVUs; red shading: IMPV > 60 PVUs. Levels of shading in (d)–(f) denote no shading: θe < 345 K; grey shading: 345 < θe < 351 K; blue shading: 351 < θe < 357 K; pink shading: 357 < θe < 363 K; red shading: θe > 363 K. The dotted circle is the radius of maximum winds, r = 37 km. The “+” is to mark the position of a typical enhancement of IMPV interior to the RMW. Cyclonic motion around this IMPV enhancement appears to position high-θe air radially outward to the RMW around the southeast side of the feature. Figure adapted from Fulton (2001)
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
IMPV at z = 1135 m at (a) 268, (b) 272, and (c) 276 min into an MM5 simulation of Hurricane Bob, kindly provided by S. A. Braun, NASA/GSFC (see Braun 2002 for details). Data is reported at the 1.3-km grid spacing of the innermost mesh. (d)–(f) Equivalent potential temperature θe displayed at respective times. The hurricane symbol denotes the storm center. Levels of shading in (a)–(c) denote no shading: IMPV < 0 PVUs; grey shading: 0 < IMPV < 15 PVUs; blue shading: 15 < IMPV < 30 PVUs; pink shading: 30 < IMPV < 60 PVUs; red shading: IMPV > 60 PVUs. Levels of shading in (d)–(f) denote no shading: θe < 345 K; grey shading: 345 < θe < 351 K; blue shading: 351 < θe < 357 K; pink shading: 357 < θe < 363 K; red shading: θe > 363 K. The dotted circle is the radius of maximum winds, r = 37 km. The “+” is to mark the position of a typical enhancement of IMPV interior to the RMW. Cyclonic motion around this IMPV enhancement appears to position high-θe air radially outward to the RMW around the southeast side of the feature. Figure adapted from Fulton (2001)
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Dropsonde retrieval of θe (solid curve) from the eye of Hurricane Jimena (eastern Pacific, 2058 UTC 23 Sep 1991) compared to θe measured on board aircraft at the moment of eye penetration (large dots), defined by the last measurement of saturation. Figure adapted with kind permission from Willoughby (1998)
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Dropsonde retrieval of θe (solid curve) from the eye of Hurricane Jimena (eastern Pacific, 2058 UTC 23 Sep 1991) compared to θe measured on board aircraft at the moment of eye penetration (large dots), defined by the last measurement of saturation. Figure adapted with kind permission from Willoughby (1998)
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Dropsonde retrieval of θe (solid curve) from the eye of Hurricane Jimena (eastern Pacific, 2058 UTC 23 Sep 1991) compared to θe measured on board aircraft at the moment of eye penetration (large dots), defined by the last measurement of saturation. Figure adapted with kind permission from Willoughby (1998)
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Soundings for the eye (solid; computed at r = 0 km), the eyewall (dotted; computed at the radius of peak updraft), and the environment (dashed; averaged over r > 800 km). The eyewall sounding shows jaggedness because of the discrete outward jumps in radial grid steps of the peak updraft
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Soundings for the eye (solid; computed at r = 0 km), the eyewall (dotted; computed at the radius of peak updraft), and the environment (dashed; averaged over r > 800 km). The eyewall sounding shows jaggedness because of the discrete outward jumps in radial grid steps of the peak updraft
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2
Soundings for the eye (solid; computed at r = 0 km), the eyewall (dotted; computed at the radius of peak updraft), and the environment (dashed; averaged over r > 800 km). The eyewall sounding shows jaggedness because of the discrete outward jumps in radial grid steps of the peak updraft
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Equivalent potential temperature Θe (K; red, 1-day average) displayed as a function of absolute angular momentum and height at day 20 in the 4x run. Overlaid are 2D streamlines of the daily averaged secondary circulation (blue). Light green shading denoted ql > 0.3 g kg–1. Dark green shading denotes ql > 1 g kg–1. The inset figure repeats the shading for ql, displayed as a function of physical radius. The contours of absolute angular momentum (black, inset) are provided for physical reference for the main figure
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2

Equivalent potential temperature Θe (K; red, 1-day average) displayed as a function of absolute angular momentum and height at day 20 in the 4x run. Overlaid are 2D streamlines of the daily averaged secondary circulation (blue). Light green shading denoted ql > 0.3 g kg–1. Dark green shading denotes ql > 1 g kg–1. The inset figure repeats the shading for ql, displayed as a function of physical radius. The contours of absolute angular momentum (black, inset) are provided for physical reference for the main figure
Citation: Journal of the Atmospheric Sciences 60, 19; 10.1175/1520-0469(2003)060<2349:HS>2.0.CO;2