## 1. Introduction

It has been established that numerically simulated tropical cyclone intensity can exceed the maximum intensity that is predicted by analytic theories. In particular, the steady-state analytical maximum potential intensity (MPI) theory of Emanuel (1986, 1995)—hereafter referred to as E-MPI—has recently been evaluated against time-dependent numerical model simulations by Persing and Montgomery [2003 (hereafter PM03), 2005], Hausman et al. (2006), and Cram et al. (2007). Additionally, recent observational studies (e.g., Tonkin et al. 2000; Montgomery et al. 2006; Bell and Montgomery 2008) have found cases in which observed intensities are greater than E-MPI. The fact that E-MPI can be significantly weaker (by ∼10%–50%) than simulated and/or observed tropical cyclone intensity is not challenged here.

What remains to be established is *why* E-MPI can be significantly weaker than the intensities of some simulated and observed tropical cyclones. To this end, a hypothesis was put forth by PM03, which has further been evaluated using observations by Montgomery et al. (2006), Aberson et al. (2006), and Bell and Montgomery (2008) and using numerical simulations by Persing and Montgomery (2005) and Cram et al. (2007). The theory posits that the locally high-entropy air at low levels in the tropical cyclone’s eye can provide an additional source of energy that is not considered in E-MPI. This process has been referred to as the “superintensity” mechanism (PM03; Cram et al. 2007) and as the “turbo-boost” mechanism (Montgomery et al. 2006; Bell and Montgomery 2008). We will refer to it hereafter as the PM03 mechanism.

Evidence is clear in both observations and numerical simulations that there can be anomalously (and significantly; i.e., ∼10–20 K) higher equivalent potential temperature (*θ _{e}*) in the low levels of tropical cyclone eyes as compared to

*θ*in the eyewall. It is also clear that this air can be transported from the eye into the eyewall. It is not clear to us that this process has any significant effect on maximum azimuthally averaged hurricane intensity. Furthermore, to our knowledge, the

_{e}*total magnitude*of

*θ*transport from the eye to the eyewall has never been evaluated in comparison to

_{e}*θ*transport from other locations of the tropical cyclone; in other words, even though

_{e}*some*high-

*θ*air from the eye is transported into the eyewall, its effect on average

_{e}*θ*in the eyewall has never been quantified clearly.

_{e}The purpose of this article is to evaluate quantitatively the PM03 mechanism. To that end, we run experiments with a time-dependent axisymmetric numerical model and we conduct control volume analyses around the eye and eyewall using the model output. We demonstrate that the PM03 mechanism does happen as envisioned but that it is quantitatively small and has only a small effect on maximum axisymmetric intensity. Furthermore, we demonstrate that this mechanism cannot explain the subintensity of E-MPI compared to hurricane intensity in these numerical simulations.

## 2. Methodology

As did PM03, we use the nonhydrostatic axisymmetric numerical model developed by Rotunno and Emanuel (1987, hereafter RE87). Our default configuration, hereafter referred to as the Control simulation, is identical to the 4x-resolution simulations by PM03, except for the initial sounding. We use the same sounding that was used by RE87 (Fig. 1a), but we have interpolated their sounding to a higher-resolution grid; below their lowest model level and above their highest model level, we extrapolate potential temperature and mixing ratio to new grid points as needed. In contrast, the sounding used in the 4x simulations by PM03 is warmer in the boundary layer and cooler above the boundary layer (Fig. 1b); this different sounding was arrived at through interpolation error and by different assumptions near the bottom and top of the model (J. Persing 2007, personal communication). We have retained their sounding for simulations shown later in this article to allow for comparisons with results published by PM03. However, we use the RE87 sounding (Fig. 1a) for most of our analysis. The sea surface temperature (*T _{s}*) is 26.13°C in simulations with this sounding. All simulations reported here achieve a steady intensity after

*t*= 8 days, so we report the average of maximum tangential velocity (

*υ*

_{max}) over

*t*= 8–12 days as a measure of maximum tropical cyclone intensity.

For calculations of entropy (*s*) and equivalent potential temperature (*θ _{e}*), we use the approximate pseudoadiabatic formulations that were presented by Bryan (2008). These formulations use a constant latent heat of vaporization (

*L*

_{0}), which can be adjusted to increase accuracy for a specific numerical model. Using the code from Bryan (2008), and using the equations and constants from the RE87 model for the reference solution, we find that

*L*

_{0}= 2.604 × 10

^{6}J kg

^{−1}is appropriate for this model. The maximum error in

*θ*is 0.4 K.

_{e}For the sake of reference, we have calculated estimates of E-MPI for the Control sounding using the “a priori” method of Emanuel (1986), which uses only parameters from an environmental sounding. Specifically, we use (43) of Emanuel (1986), with *r*_{0} = 400 km, *f* = 5 × 10^{−5} s^{−1}, *T _{B}* = 292 K,

*p*= 1015.1 mb, and

_{a}*C*

_{θ}/

*C*= 1 [see Emanuel (1986) for definitions of these variables]. Given uncertainty in the outflow temperature (

_{D}*T*

_{0}) and the ambient boundary layer relative humidity (RH

*), we use a range of values for these two parameters. For*

_{as}*T*= 26.13°C, we find a minimum value of E-MPI = 37 m s

_{s}^{−1}, with

*T*

_{0}= 220 K and RH

*= 0.9, and a maximum value of E-MPI = 60 m s*

_{as}^{−1}, with

*T*

_{0}= 200 K and RH

*= 0.8.*

_{as}Given the large inherent uncertainty in the calculation of a priori E-MPI, hereafter we use instead the “local PBL balance” method from PM03 (their p. 2352). Evaluation is performed at the radius of maximum winds using output from the numerical simulation. To obtain a value for *T*_{0}, we compute a trajectory using the mean flow fields from the model (averaged from *t* = 8–12 days); for the parcel that passes through the maximum tangential winds, we calculate the average temperature along the parcel trajectory as it passes between *r* = 200–300 km. This calculation of *T*_{0} is consistent with the definition by Emanuel (1986, p. 587) as “the temperature the air has along *M* surfaces as they flare out to very large radii.”

## 3. Results with the RE87 environment

After about 5 simulated days, the tropical cyclone in the simulation with the RE87 sounding (Control) reaches an approximately steady intensity of 86 m s^{−1} (Table 1). We estimate E-MPI to be 56 m s^{−1} (Table 1). Clearly, this theoretical value is significantly less than the simulated storm’s intensity. We also note that this magnitude of subintensity is comparable to that found by PM03.

Near the surface in the eye, *θ _{e}* is 12 K higher than

*θ*at the radius of maximum winds (Fig. 2a). The magnitude of this

_{e}*θ*difference is comparable both to that in the simulation by PM03 and to that observed in Hurricane Isabel (Bell and Montgomery 2008).

_{e}To identify the processes that allow this *θ _{e}* pattern to exist, we have calculated a

*θ*budget using the same technique as RE87 (see p. 555), except we use the formulation of

_{e}*θ*from Bryan (2008). We confirm (following PM03, p. 2355) that the ultimate source of the high-

_{e}*θ*anomaly is the upward flux of entropy from the sea surface in the eye. All other processes (i.e., advection, horizontal turbulent diffusion, and “radiation”) act to lower

_{e}*θ*in the low-level eye.

_{e}Given that surface entropy fluxes in the eye ultimately create the high-*θ _{e}* anomaly, we conduct a series of simulations in which the surface fluxes are set to zero in the eye. This eliminates the high-

*θ*anomaly and thus allows us to evaluate the role it plays in maximum hurricane intensity. To this end, we set the surface exchange coefficient for entropy (

_{e}*C*; see RE87, p. 547) to zero for

_{E}*r*< 17 km (the average radius of maximum winds from Control). We note that using the radius of maximum winds to delineate the eye is a very generous definition; that is, this choice probably makes “the eye” too large. In fact, some grid points inside the eye, thus defined, are actually cloudy (and precipitating). We make this choice primarily because of similar simulations conducted by J. Persing and M. Montgomery (2007, personal communication) and because it biases the results toward the largest effect one could expect to find from fluxes in the eye. In other words, we give the eye fluxes the best possible chance to have a significant effect.

We first analyze a simulation in which surface fluxes are zero for *r* < 17 km for the entire simulation (referred to as NoEyeFlux). As expected, the high-*θ _{e}* anomaly in the eye does not form in this simulation (Fig. 2b). However, there is very little difference in maximum intensity compared to Control (Fig. 2c). The difference in

*υ*

_{max}is only 3 m s

^{−1}(Table 1). In Fig. 2c, the maximum difference in

*υ*at any grid point is actually 9 m s

^{−1}, but this occurs inside the eye itself; the maximum intensity in both simulations is located at

*r*≈ 17 km, where the difference is 3 m s

^{−1}. This difference in intensity is much less than the ∼30 m s

^{−1}subintensity of E-MPI for this case. This result contradicts the statement that numerically simulated intensity is much greater than E-MPI “only in the presence of an enhancement of low-level eye entropy” (PM03, p. 2349).

We also conduct a simulation in which surface fluxes are set to zero for *r* < 17 km after 7 days only (referred to as NoEyeFlux-7d). This simulation is identical to Control for the first 7 days, which is long enough into the simulation for the high-*θ _{e}* anomaly to form and reach a steady state. After the surface fluxes in the eye are turned off, the high-

*θ*anomaly is quickly eradicated (in about 8 h) (Fig. 3c). Notably, the maximum sustained winds are hardly affected by this change;

_{e}*υ*

_{max}decreases by only 2 m s

^{−1}(Table 1).

In addition, we conduct a simulation in which we track the radius of maximum winds (*r _{m}*) at every time step during the simulation and set

*C*to zero inward of this point after some nominal intensity (30 m s

_{E}^{−1}) is achieved (referred to as NoEyeFlux-rm). In this case, the zero surface fluxes are not confined to a specified region. Furthermore, when the intensity threshold of 30 m s

^{−1}is first exceeded (at

*t*= 3 days), the region of zero fluxes extends initially to

*r*≈ 50 km. Nevertheless, the same overall result is ultimately obtained. That is,

*r*eventually decreases to 17 km (Fig. 3d), and the steady maximum intensity is only slightly less than that from the Control simulation (Table 1).

_{m}As mentioned earlier, our definition for the eye in these simulations is quite generous and includes grid points that have precipitation. Consequently, these simulations have zero surface fluxes over part of the eyewall, conventionally defined. To investigate further, we have also conducted simulations that retain the fluxes over all of the eyewall (defined either as precipitating grid points, or as grid points with vertical velocity exceeding 1 m s^{−1} at *z* = 1 km). Results (not shown) show that the high-*θ _{e}* anomaly in the eye is eliminated in these simulations but that

*υ*

_{max}is unchanged from Control. Thus, the small (∼3 m s

^{−1}) change noted in Table 1 is an overestimate of the effect on maximum hurricane intensity due to the high-

*θ*anomaly in the low-level eye.

_{e}## 4. Simulations in other environments

To evaluate whether our results are sensitive to the environment chosen for the initial state, we also conduct simulations with the sounding used for the 4x simulations by PM03 (Fig. 1b), with *T _{s}* = 26.13°C, and also with the mean hurricane sounding by Jordan (1958) (Fig. 1c), with

*T*= 28.0°C. The latter environment is identical to that used by Hausman et al. (2006), and is nearly the same as the “Mid CAPE” simulation by Persing and Montgomery (2005). As in the previous section, we conduct additional simulations in which the high-

_{s}*θ*anomaly is eliminated by setting the surface fluxes to zero inside the radius of maximum winds. We use all three of the methodologies described in the previous section, although in these two environments the nominal value of

_{e}*r*is 24 km. As before, we bias the results toward the greatest possible effect of the high-

_{m}*θ*anomaly by including all grid points for

_{e}*r*<

*r*.

_{m}Consistent with the simulations in the previous section, the high-*θ _{e}* anomaly is completely eliminated when there are no surface fluxes in the eye. In the simulations without the eye fluxes,

*υ*

_{max}is again slightly weaker than in simulations that retain these fluxes (Table 1). Of all our experiments, the maximum decrease in

*υ*

_{max}is 7.5 m s

^{−1}and the average decrease is 4.5 m s

^{−1}. Nevertheless,

*υ*

_{max}in all simulations is 25–30 m s

^{−1}higher than E-MPI (Table 1).

Out of all simulations in which eye fluxes are included, the simulation JordanControl has the largest *θ _{e}* anomaly; the maximum

*θ*difference is 25 K (in which

_{e}*θ*in the eye is compared to

_{e}*θ*at the radius of maximum winds). This

_{e}*θ*difference is comparable to the maximum value documented in observations (e.g., Eastin et al. 2005b; Montgomery et al. 2006; Marks et al. 2008; Bell and Montgomery 2008), which suggests that our simulations are producing features that are consistent with observations. Although some of these observational studies suggested that this feature has a significant impact on the intensity of tropical cyclones, our simulations demonstrate that this feature has a negligible impact on maximum axisymmetric intensity.

_{e}We have repeated these simulations using an axisymmetric version of the cloud model developed by Bryan and Fritsch (2002). This newer model uses more accurate numerical techniques, conserves mass and energy, and includes dissipative heating; all of these are significant improvements over the RE87 model. We also use smaller grid spacing (Δ*r* = 1000 m; Δ*z* = 250 m) for these simulations. Our overall conclusions from these additional simulations are the same. The maximum decrease of *υ*_{max} when we remove the eye fluxes in these simulations is 4 m s^{−1}, and this occurs in the simulation JordanNoEyeFlux-rm.

## 5. Interpretation

We conclude, based on these simulations, that the PM03 mechanism has a negligible effect on maximum axisymmetric intensity and that it does not explain the significant subintensity of E-MPI for these simulated hurricanes. We acknowledge that the process identified by PM03 does happen [as shown, for example, by Cram et al. (2007)], but these numerical simulations demonstrate that it has a very small impact on azimuthally averaged winds. [The PM03 mechanism clearly affects the intensity of individual, unsteady, nonaxisymmetric features as measured by local buoyancy and/or updraft intensity, as shown by Braun (2002), Eastin et al. (2005a, b), and Cram et al. (2007). The focus here is on the *steady, axisymmetric intensity*, which these simulations show to be barely affected by the PM03 mechanism.] It also follows, based on these simulations, that the PM03 mechanism can be reasonably neglected in analytic MPI theories.

We reiterate that by shutting off the surface fluxes in the eye, we are *directly* evaluating the PM03 mechanism (as explained in section 3). In contrast, PM03 attempted to evaluate this mechanism by inserting a heat sink (via Newtonian relaxation) in the low-level eye (see their section 5c). By doing so, they not only countered the effect of the surface fluxes in the eye (which is the ultimate goal of such tests), but they also countered all other processes that contribute toward the warm-core structure of the low-level eye. Not surprisingly, they found a substantial decrease in intensity in their simulations (∼15 m s^{−1} decrease in *υ*_{max}) and an adjustment of the hurricane’s structure (i.e., an increase in radius of maximum winds) as the hurricane responded to this large heat sink. The weak (∼4 m s^{−1}, on average) decrease in intensity reported here is a more direct measure of the PM03 mechanism.

The small impact of the high-*θ _{e}* anomaly in these simulations is ultimately attributable to several factors, which we investigate in the next several subsections.

### a. Magnitude of surface fluxes

*ρ*is density,

*u*

_{0}and

*υ*

_{0}are the radial and tangential velocities at the lowest model level,

*s*

_{surf}is surface entropy at saturation at temperature

*T*, and

_{s}*s*

_{0}is entropy at the lowest model level. Because

*F*just inside the eye is comparable to

^{s}*F*in the eyewall (Fig. 4a), we might conclude that surface fluxes in the eye are very important; in comparison, fluxes far outside of the eyewall (

^{s}*r*> 40 km) seem unimportant. However, energetically, fluxes integrated over area (as opposed to fluxes at individual grid points) are ultimately important for tropical cyclone intensity. We therefore need to calculate the area-integrated flux, which is given by

*r*is the radial grid spacing. This analysis, shown in Fig. 4b, suggests a much larger potential impact from the surface fluxes in the region well beyond the radius of maximum winds (i.e., for 50 km <

*r*< 100 km in Fig. 4b).

*θ*increases with decreasing radius in the boundary layer; this is the region where upward flux of

_{e}*θ*from the ocean is approximately balanced by inward radial advection of

_{e}*θ*toward the eyewall [referred to as region II by Emanuel (1986, p. 594)]. Based on plots of

_{e}*θ*, as well as our

_{e}*θ*budget analysis (not shown), we find the outer radius

_{e}*R*of this region to be at

*r*= 200 km. (In an analysis using the same numerical model and the same sounding, but different model resolution, RE87 also concluded that there is a region, extending beyond the eyewall, where surface fluxes balance radial advection of

*θ*; in their case, they found

_{e}*R*≈ 240 km.) Based on this analysis, we normalize the storm-integrated surface fluxes as a function of

*r*by the total value as follows:

The analysis in this subsection has focused solely on the energetics of tropical cyclones, but it is not clear how these surface fluxes affect the *intensity* of tropical cyclones (other than the fact that surface fluxes in the eye have already been shown to have a negligible impact on intensity, as shown in Table 1). To investigate further, we conduct a series of simulations in which the surface entropy fluxes are set to zero from *r* = 0 to some point *r* = *r _{f}* , where

*r*is specified at the beginning of each simulation. In Fig. 5, we show results for a range of values between

_{f}*r*= 3.75 km (i.e., one grid point) to

_{f}*r*= 100 km for all three environments studied here. Each data point on this figure is

_{f}*υ*

_{max}from a different simulation, and we note that

*r*= 0 is the control simulation for each case. Regardless of the environment used,

_{f}*υ*

_{max}does not decrease until surface fluxes are set to zero inside the nominal eyewall region (shaded region in Fig. 5). With further increases in

*r*, there is a gradual reduction in

_{f}*υ*

_{max}. For simulations in which

*r*>

_{f}*R*(not shown), no steady hurricane is produced in the simulations. These simulations support the conclusion that surface fluxes over a broad region (specifically, from the eyewall to

*r*=

*R*) are important to the energetics, and thus to the intensity, of hurricanes.

### b. Entropy transport

*θ*anomaly. We begin with a governing equation for

_{e}*θ*assuming steady flow:

_{e}*u*and

*w*are the radial and vertical velocities,

*F*

^{θe}

_{r}and

*F*

^{θe}

_{z}are the parameterized turbulent fluxes of

*θ*in the radial and vertical directions,

_{e}*ρ*

*θ*

_{υ}is the base-state virtual potential temperature. We neglect the small contribution from the “radiative cooling” term in the model. To write a flux-form version of this equation, we utilize the model’s equation for mass-conservation in steady flow [following RE87’s Eq. (4)]:

*θ*can be expressed as

_{e}*θ*anomaly in our set of simulations. The steady flow from this simulation suggests that a convenient control volume (dashed box in Fig. 6) can be defined from

_{e}*r*= 0 to

*r*=

*r*, where

_{m}*r*is the radius of maximum winds, and from

_{m}*z*= 0 to

*z*=

*H*, where

*H*is the top of the control volume (at

*z*= 3.5 km). We integrate (6) over this control volume:

*u*|

_{r=0}=

*w*|

_{z=0}=

*F*

^{θe}

_{r}|

_{r=0}= 0. Further, we find from analysis of model output that terms with

*F*

^{θe}

_{r}|

_{r=rm},

*F*

^{θe}

_{z}|

_{z=H}, and

*w*|

_{z=H}are small and therefore negligible compared to other terms on the boundaries. The control volume budget thus reduces to

*θ*through the right side of the control volume [left side of (8)] balances the total flux of

_{e}*θ*from the sea surface [right side of (8)]. With (8), it is straightforward to demonstrate how we can remove the PM03 mechanism from our simulations by simply setting

_{e}*F*

^{θe}

_{z}|

_{z=0}= 0 for

*r*<

*r*; this means that the total flux through the right side of the control volume [left side of (8)] must be zero, and thus there is no net enhancement of

_{m}*θ*for parcels that pass through the eye and then back into the eyewall.

_{e}*F*

^{eye}

_{surf}. We use output from simulation JordanControl, which we then average over

*t*= 8–10 days (the time period shown in Fig. 6). We find the value of

*F*

^{eye}

_{surf}to be 2 × 10

^{11}kg K

^{−2}s

^{−1}. This is a quantitative measure of the PM03 mechanism, and it incorporates the

*total*contribution to

*θ*in the eyewall due to the eye/eyewall exchange mechanism identified by PM03. Apparently, this is a negligible contribution to

_{e}*θ*in the eyewall because setting this to zero (i.e., setting the surface fluxes in the eye to zero) has negligible effect on

_{e}*υ*

_{max}(Table 1) and a negligible impact on the

*θ*pattern in the eyewall (e.g., Fig. 2). Theoretical studies have concluded that the

_{e}*θ*distribution in the eyewall of hurricanes is important for maximum hurricane intensity (e.g., Emanuel 1986, 1997). Consequently, we compare

_{e}*F*

^{eye}

_{surf}to other terms in the

*θ*budget around the eyewall (dotted box in Fig. 6). We find that

_{e}*F*

^{eye}

_{surf}is one order of magnitude smaller than surface fluxes in the eyewall:

*r*

_{1}and

*r*

_{2}are radii that encompass the eyewall updraft, as shown in Fig. 6. We also find that

*F*

^{eye}

_{surf}is three orders of magnitude smaller than the

*θ*flux on the right side of the control volume (which incorporates all surface entropy fluxes between

_{e}*r*

_{2}and

*R*). In summary, consistent with the numerical experiments (Table 1 and Fig. 2), the PM03 mechanism is a negligible component of the entropy budget of the eyewall in these simulated hurricanes.

Using a Lagrangian-based analysis, Cram et al. (2007) found that parcels that pass from the low-level eye to the eyewall tended to have *θ _{e}* that was several K higher (on average) than parcels that did not enter the eye (e.g., their Fig. 9). In contrast to our findings, they concluded that the eye/eyewall exchange mechanism can “benefit the intensity of the storm.” One might wonder whether the simulated hurricane used by Cram et al. (2007) is significantly different from the simulated hurricanes used here and whether this could lead to their different conclusion. We note from JordanControl that average

*θ*exiting the eye (between

_{e}*z*=

*z*and

_{m}*z*=

*H*) is 4 K higher than average

*θ*entering the eye (between

_{e}*z*= 0 and

*z*=

*z*) (Fig. 6). This average value of

_{m}*θ*difference is comparable to that found by Cram et al. (2007) using Lagrangian parcel analysis. Thus, it seems likely that the magnitude of the PM03 mechanism is similar in these simulations.

_{e}Even though air parcels *are* exchanged between the eye and eyewall [e.g., Braun (2002), PM03, Cram et al. (2007)], and these parcels *do have* anomalously high *θ _{e}* such that they could contribute to higher mean

*θ*in the eyewall, it turns out that this process is just not quantitatively significant in these simulations (e.g., Fig. 2; the eyewall spans roughly 15 km <

_{e}*r*< 25 km). It follows that there must be a relatively small amount of mass that passes between the eye and the eyewall as compared to the total mass in the eyewall. Using a control-volume analysis of mass continuity around the eyewall (i.e., dotted box in Fig. 6), we find that the mass flux from the eye to the eyewall is only 8% of the total upward mass flux at

*z*=

*H*in the eyewall. Considering that

*θ*is elevated by 4 K (on average) for the parcels that pass through the eye, this implies that average

_{e}*θ*in the eyewall may be elevated by 0.3 K (at most) by the PM03 mechanism (assuming linear mixing of

_{e}*θ*). This analysis is consistent with results from numerical experiments (e.g., Fig. 2); in fact, average

_{e}*θ*in the eyewall is only a few tenths of a degree Kelvin higher in simulations with the PM03 mechanism as compared to simulations without the PM03 mechanism.

_{e}### c. A sensitivity test

We conclude, based on the preceding analyses, that the PM03 mechanism—and, hence, surface fluxes in the eye—would have to be much larger in magnitude to have a significant effect on the energetics of hurricanes. To estimate how much larger, we note that the total surface entropy flux in the nominal eyewall region (*F*^{eyewall}_{surf}) is one order of magnitude larger than total surface flux in the eye (*F*^{eye}_{surf}). From the sensitivity simulations shown in Fig. 5, we find that setting (*F*^{eyewall}_{surf}) to zero results in approximately a 15% decrease in *υ*_{max}. Therefore, we surmise that *F*^{eye}_{surf} would probably have to be one order of magnitude larger to have a significant effect on *υ*_{max}.

To test this hypothesis, we conduct simulations in which *C _{E}* in the eye is increased by a factor of 10. For these simulations, we track the radius of maximum winds every time step (as in the NoEyeFlux-rm simulations) and we multiply

*C*by 10, as compared to the formulation from the default model. Results are shown in Table 2. We find that maximum intensity increases by 5–14 m s

_{E}^{−1}, with the greatest difference occurring in the Jordan mean hurricane environment. These results support our hypothesis that surface fluxes in the eye would have to be at least one order of magnitude higher than typical values to have a significant positive effect on maximum axisymmetric hurricane intensity. However, we note that the

*θ*anomalies in these simulations are much larger (in amplitude) than values that have been reported in observations (e.g., Eastin et al. 2005b; Montgomery et al. 2006; Marks et al. 2008; Bell and Montgomery 2008). In the experiment with the Jordan environment,

_{e}*θ*in the eye is 46 K higher than

_{e}*θ*at the radius of maximum winds, which is a factor of 2 larger than the maximum observed amplitude of the

_{e}*θ*anomaly.

_{e}### d. Comparison to observations

*F*

^{eye}

_{surf}can be estimated using either side of (8). The left side of (8) is especially convenient because it contains terms that are easily observed. To make the following interpretation easier, we simplify this expression. First, we divide the control volume into two parts: a bottom part where

*u*< 0, from

*z*= 0 to

*z*=

*z*; and an upper part where

_{m}*u*> 0, from

*z*=

*z*to

_{m}*z*=

*H*. Next, for simplicity, we further assume that

*u*and

*θ*in the bottom part can be represented by average values

_{e}*u*

_{1}and

*θ*

_{e}_{1}, and the same variables in the upper part can be represented by average values

*u*

_{2}and

*θ*

_{e}_{2}. By further assuming that

*ρ*

*θ*

_{υ}≈ constant, we arrive at

*H*−

*z*)

_{m}*u*

_{2}= −

*z*

_{m}u_{1}. Thus, we arrive at an approximate relation:

*F*

^{eye}

_{surf}needs to be

*at least*one order of magnitude larger to have a significant impact. The right side of (11) suggests that this would be difficult, if not impossible, to accomplish. The depth of the inflowing layer into the eye (

*z*) seems to have little variation in observations (e.g., Bell and Montgomery 2008, their Fig. 4). The radius of maximum winds

_{m}*r*certainly would not be expected to be substantially larger for very intense storms (i.e., those having intensities that greatly exceed E-MPI). Furthermore, values of

_{m}*u*

_{1}and (

*θ*

_{e}_{2}−

*θ*

_{e}_{1}) from available observations (e.g., Bell and Montgomery 2008) are comparable to values from the JordanControl simulation. In summary, to raise the right side of (11) by an order of magnitude would require a drastic change in hurricane structure that seems highly unlikely. We conclude from this analysis that surface fluxes in the eye are unlikely to have a significant effect (on maximum axisymmetric intensity) in observed tropical cyclones. Further analysis with observations and/or numerical simulations in different environments, and for hurricanes with different

*r*, and with three-dimensional effects, would be needed to confirm this conclusion, especially because indirect effects are not considered in this relatively simple analysis.

_{m}## 6. Summary

We use a nonhydrostatic axisymmetric numerical model to evaluate the “superintensity” (or “turbo-boost”) mechanism put forth by PM03. The theory posits that anomalously high-*θ _{e}* air from near the surface in hurricane eyes is transported into the eyewall and that this mechanism (which is not considered in analytic theories for maximum intensity) can substantially increase the maximum intensity of hurricanes. We find that although the PM03 mechanism does occur in the numerical simulations, it is too small in magnitude to significantly affect the maximum axisymmetric intensity of hurricanes. In simulations in which we set surface fluxes in the eye to zero, which effectively and quickly eliminates the high-

*θ*anomaly, the intensity of the cyclone is only slightly affected (the typical decrease in intensity is ∼4%, despite a very generous definition for “the eye”). We also find that all of our simulations retain an intensity that is 25–30 m s

_{e}^{−1}greater than the maximum theoretical intensity of Emanuel (1986) (E-MPI; see Table 1), regardless of whether there is a high-entropy anomaly in the low-level eye.

From further analysis, we find that total surface entropy fluxes in the eye are negligible compared to total surface fluxes near and outside the eyewall, and therefore the total magnitude of *θ _{e}* flux between the eye and eyewall is a negligible component of the

*θ*budget in the eyewall. A cursory comparison to observations suggests that no combination of relevant parameters (e.g., radius of maximum winds,

_{e}*θ*distribution in the eye) is likely to occur that would suggest that the PM03 mechanism is important in natural hurricanes. Additional sensitivity tests with the numerical model support these conclusions.

_{e}With regard to future work on this topic, additional numerical simulations could be undertaken in different environments (i.e., with different soundings and/or *T _{s}*) to check whether our results are sufficiently general. Furthermore, a similar set of numerical simulations could be undertaken in three dimensions to investigate whether asymmetric features play any significant role that is not considered here. To that end, a recent study by Yang et al. (2007) compared results from three-dimensional simulations to comparable axisymmetric simulations. Based on figures in their article, there do not appear to be any significant changes in the relevant proxy parameters [e.g.,

*r*,

_{m}*z*,

_{m}*u*

_{1}; see (11)] that would suggest any significant changes in surface fluxes from the eye in their simulations.

Finally, although the mechanism proposed by PM03 might seem like a reasonable qualitative explanation for the subintensity of E-MPI, our quantitative analysis shows that it is far too weak to explain the discrepancy. We do not propose an alternative explanation for the subintensity of E-MPI here, although we are exploring this issue for a future article. Nevertheless, based on the analyses here, we conclude that the neglect of surface fluxes in the eye is not an inherent limitation of E-MPI.

## Acknowledgments

We thank John Persing and Mike Montgomery for providing reviews of our manuscript, for sharing their code, and for our many discussions on this topic. We also thank Chris Davis, Greg Holland, and Kerry Emanuel for their reviews of earlier versions of this manuscript.

## REFERENCES

Aberson, S. D., M. T. Montgomery, M. Bell, and M. Black, 2006: Hurricane Isabel (2003): New insights into the physics of intense storms. Part II: Extreme localized wind.

,*Bull. Amer. Meteor. Soc.***87****,**1349–1354.Bell, M. M., and M. T. Montgomery, 2008: Observed structure, evolution, and potential intensity of category 5 Hurricane Isabel from 12 to 14 September.

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

,*Mon. Wea. Rev.***130****,**1573–1592.Bryan, G. H., 2008: On the computation of pseudoadiabatic entropy and equivalent potential temperature.

, in press.*Mon. Wea. Rev.*Bryan, G. H., and J. M. Fritsch, 2002: A benchmark simulation for moist nonhydrostatic numerical models.

,*Mon. Wea. Rev.***130****,**2917–2928.Cram, T. A., J. Persing, M. T. Montgomery, and S. A. Braun, 2007: A Lagrangian trajectory view on transport and mixing processes between the eye, eyewall, and environment using a high-resolution simulation of Hurricane Bonnie (1998).

,*J. Atmos. Sci.***64****,**1835–1856.Eastin, M. D., W. M. Gray, and P. G. Black, 2005a: Buoyancy of convective vertical motions in the inner core of intense hurricanes. Part I: General statistics.

,*Mon. Wea. Rev.***133****,**188–208.Eastin, M. D., W. M. Gray, and P. G. Black, 2005b: Buoyancy of convective vertical motions in the inner core of intense hurricanes. Part II: Case studies.

,*Mon. Wea. Rev.***133****,**209–227.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., 1995: Sensitivity of tropical cyclones to surface exchange coefficients and a revised steady-state model incorporating eye dynamics.

,*J. Atmos. Sci.***52****,**3969–3876.Emanuel, K. A., 1997: Some aspects of hurricane inner-core dynamics and energetics.

,*J. Atmos. Sci.***54****,**1014–1026.Hausman, S. A., K. V. Ooyama, and W. H. Schubert, 2006: Potential vorticity structure of simulated hurricanes.

,*J. Atmos. Sci.***63****,**87–108.Jordan, C. L., 1958: Mean soundings for the West Indies area.

,*J. Meteor.***15****,**91–97.Marks, F. D., P. G. Black, M. T. Montgomery, and R. W. Burpee, 2008: Structure of the eye and eyewall of Hurricane Hugo (1989).

,*Mon. Wea. Rev.***136****,**1237–1259.Montgomery, M. T., M. M. Bell, S. D. Aberson, and M. L. Black, 2006: Hurricane Isabel (2003): New insights into the physics of intense storms. Part I: Mean vortex structure and maximum intensity estimates.

,*Bull. Amer. Meteor. Soc.***87****,**1335–1347.Persing, J., and M. T. Montgomery, 2003: Hurricane superintensity.

,*J. Atmos. Sci.***60****,**2349–2371.Persing, J., and M. T. Montgomery, 2005: Is environmental CAPE important in the determination of maximum possible hurricane intensity?

,*J. Atmos. Sci.***62****,**542–550.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.Tonkin, H., G. J. Holland, N. Holbrook, and A. Henderson-Sellers, 2000: An evaluation of thermodynamic estimates of climatological maximum potential tropical cyclone intensity.

,*Mon. Wea. Rev.***128****,**746–762.Yang, B., Y. Wang, and B. Wang, 2007: The effect of internally generated inner-core asymmetries on tropical cyclone potential intensity.

,*J. Atmos. Sci.***64****,**1165–1188.

Output averaged from *t* = 8–12 days: (a) from Control, (b) from NoEyeFlux, and (c) the difference between the two simulations (NoEyeFlux − Control). In (a) and (b), thin contours are *θ _{e}* every 1 K, shading is

*θ*≥ 360 K, and thick contours are

_{e}*υ*every 20 m s

^{−1}. In (c), shading is

*θ*difference (as indicated by color bar on the right), and contours are

_{e}*υ*difference every 2 m s

^{−1}with the zero contour excluded.

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Output averaged from *t* = 8–12 days: (a) from Control, (b) from NoEyeFlux, and (c) the difference between the two simulations (NoEyeFlux − Control). In (a) and (b), thin contours are *θ _{e}* every 1 K, shading is

*θ*≥ 360 K, and thick contours are

_{e}*υ*every 20 m s

^{−1}. In (c), shading is

*θ*difference (as indicated by color bar on the right), and contours are

_{e}*υ*difference every 2 m s

^{−1}with the zero contour excluded.

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Output averaged from *t* = 8–12 days: (a) from Control, (b) from NoEyeFlux, and (c) the difference between the two simulations (NoEyeFlux − Control). In (a) and (b), thin contours are *θ _{e}* every 1 K, shading is

*θ*≥ 360 K, and thick contours are

_{e}*υ*every 20 m s

^{−1}. In (c), shading is

*θ*difference (as indicated by color bar on the right), and contours are

_{e}*υ*difference every 2 m s

^{−1}with the zero contour excluded.

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Time series of model output from: (a) Control, (b) NoEyeFlux, (c) NoEyeFlux-7d, and (d) NoEyeFlux-rm. Shading is *θ _{e}* (K) at the lowest model level, the thick black contour is

*υ*= 75 m s

^{−1}at

*z*= 1 km, and the dots denote the location of

*υ*

_{max}. E-MPI is approximately 55 m s

^{−1}for all cases.

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Time series of model output from: (a) Control, (b) NoEyeFlux, (c) NoEyeFlux-7d, and (d) NoEyeFlux-rm. Shading is *θ _{e}* (K) at the lowest model level, the thick black contour is

*υ*= 75 m s

^{−1}at

*z*= 1 km, and the dots denote the location of

*υ*

_{max}. E-MPI is approximately 55 m s

^{−1}for all cases.

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Time series of model output from: (a) Control, (b) NoEyeFlux, (c) NoEyeFlux-7d, and (d) NoEyeFlux-rm. Shading is *θ _{e}* (K) at the lowest model level, the thick black contour is

*υ*= 75 m s

^{−1}at

*z*= 1 km, and the dots denote the location of

*υ*

_{max}. E-MPI is approximately 55 m s

^{−1}for all cases.

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Analysis of surface entropy fluxes from the Control simulation. Open dots denote grid points that are inside the hurricane’s eye (as defined in the experiments in which surface entropy fluxes in the eye are set to zero). (a) Magnitude of flux at individual grid points, (b) area-integrated flux, and (c) total-storm-integrated flux normalized by the total value.

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Analysis of surface entropy fluxes from the Control simulation. Open dots denote grid points that are inside the hurricane’s eye (as defined in the experiments in which surface entropy fluxes in the eye are set to zero). (a) Magnitude of flux at individual grid points, (b) area-integrated flux, and (c) total-storm-integrated flux normalized by the total value.

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Analysis of surface entropy fluxes from the Control simulation. Open dots denote grid points that are inside the hurricane’s eye (as defined in the experiments in which surface entropy fluxes in the eye are set to zero). (a) Magnitude of flux at individual grid points, (b) area-integrated flux, and (c) total-storm-integrated flux normalized by the total value.

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Values of *υ*_{max} (m s^{−1}) from simulations in which surface entropy fluxes are set to zero from *r* = 0 to *r* = *r _{f}* using: (a) the RE87 environment, (b) the PM03 environment, and (c) the Jordan mean hurricane environment. The dashed line marks the nominal location of the radius of maximum winds from control simulations (i.e., from the

*r*= 0 simulations), and the gray shaded region denotes the nominal location of the eyewall from the control simulations.

_{f}Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Values of *υ*_{max} (m s^{−1}) from simulations in which surface entropy fluxes are set to zero from *r* = 0 to *r* = *r _{f}* using: (a) the RE87 environment, (b) the PM03 environment, and (c) the Jordan mean hurricane environment. The dashed line marks the nominal location of the radius of maximum winds from control simulations (i.e., from the

*r*= 0 simulations), and the gray shaded region denotes the nominal location of the eyewall from the control simulations.

_{f}Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Values of *υ*_{max} (m s^{−1}) from simulations in which surface entropy fluxes are set to zero from *r* = 0 to *r* = *r _{f}* using: (a) the RE87 environment, (b) the PM03 environment, and (c) the Jordan mean hurricane environment. The dashed line marks the nominal location of the radius of maximum winds from control simulations (i.e., from the

*r*= 0 simulations), and the gray shaded region denotes the nominal location of the eyewall from the control simulations.

_{f}Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Output from JordanControl averaged from *t* = 8–10 days. Shading is *θ _{e}* (K), contours are

*υ*every 20 m s

^{−1}, and vectors illustrate radial and vertical flow at every grid point; vector length corresponding to 1 km on this figure is equivalent to a vector magnitude of 10 m s

^{−1}. The dashed box denotes a control volume for the low-level eye, and the dotted box denotes a control volume for the low-level eyewall (see analysis in section 5).

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Output from JordanControl averaged from *t* = 8–10 days. Shading is *θ _{e}* (K), contours are

*υ*every 20 m s

^{−1}, and vectors illustrate radial and vertical flow at every grid point; vector length corresponding to 1 km on this figure is equivalent to a vector magnitude of 10 m s

^{−1}. The dashed box denotes a control volume for the low-level eye, and the dotted box denotes a control volume for the low-level eyewall (see analysis in section 5).

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Output from JordanControl averaged from *t* = 8–10 days. Shading is *θ _{e}* (K), contours are

*υ*every 20 m s

^{−1}, and vectors illustrate radial and vertical flow at every grid point; vector length corresponding to 1 km on this figure is equivalent to a vector magnitude of 10 m s

^{−1}. The dashed box denotes a control volume for the low-level eye, and the dotted box denotes a control volume for the low-level eyewall (see analysis in section 5).

Citation: Journal of the Atmospheric Sciences 66, 1; 10.1175/2008JAS2707.1

Summary of simulations in this study. For *υ*_{max}, the value listed is an average of the maximum instantaneous value of *υ* every time step for *t* = 8–12 days. For E-MPI, average conditions from *t* = 8–12 days are used for the calculation.

Results from simulations in which *C _{E}* is increased by a factor of 10 in the eye; Δ

*υ*

_{max}is the change in

*υ*

_{max}relative to the control simulation.

^{}

* The National Center for Atmospheric Research is sponsored by the National Science Foundation.