Skip to main content

ORIGINAL RESEARCH article

Front. Earth Sci., 08 November 2022
Sec. Atmospheric Science
This article is part of the Research Topic Tropical Cyclone Intensity and Structure Changes: Theories, Observations, Numerical Modeling and Forecasting View all 19 articles

Tropical cyclone rapid intensification and the excitation of inertia-gravity waves on the edges of an evolving potential vorticity structure

Wayne H. SchubertWayne H. SchubertRichard K. Taft
Richard K. Taft*
  • Department of Atmospheric Science, Colorado State University, Fort Collins, CO, United States

The problem of tropical cyclone rapid intensification is reduced to a potential vorticity (PV) equation and a second order, inhomogeneous, partial differential equation for the azimuthal wind. The latter equation has the form of a Klein-Gordon equation, the right-hand side of which involves the radial derivative of the evolving PV field. When the PV field evolves rapidly, inertia-gravity waves are excited at the edges of the evolving PV structure. In contrast, when the PV field evolves slowly, the second order time derivative term in the Klein-Gordon equation is negligible, inertia-gravity waves are not excited, and the equation reduces to an invertibility principle for the PV. The above concepts are presented in the context of an axisymmetric shallow water model, in both its linear and nonlinear forms. The nonlinear results show a remarkable sensitivity of vortex intensification to the percentage of mass that is diabatically removed from the region inside a given absolute angular momentum surface.

1 Introduction

While modern numerical weather prediction has produced significant skill in forecasting tropical cyclone tracks, the forecasting of rapid intensification remains an important challenge (DeMaria et al., 2014). The status of present observational, theoretical, and numerical modeling efforts to improve intensity forecasting can be found in the special collection of thirty articles that have appeared recently in five of the journals published by the American Meteorological Society (see the Doyle and Ferek, 2017, reference for the website that describes this special collection associated with the Tropical Cyclone Intensity Experiment). A concise and excellent review of much of this previous work can be found in Martinez et al. (2019). The present paper explores some of the theoretical aspects of the rapid intensification problem. In particular, this paper discusses how the problem of tropical cyclone intensification can be reduced to the solution of a coupled pair of partial differential equations—a first order equation for the time evolution of the potential vorticity (PV) field, and a second order, Klein-Gordon equation for the evolution of the rotational wind field. In many interesting cases the potential vorticity equation can be solved analytically and the result used as a forcing term in the Klein-Gordon equation. Examples of this approach, in the case when the Klein-Gordon equation reduces to an elliptic invertibility principle, are given in Schubert and Alworth (1987) and Möller and Smith (1994).

The forcing term in the Klein-Gordon equation can be especially large at the edges of an evolving PV structure. When this forcing effect varies slowly in time, the second order time derivative term in the Klein-Gordon equation can be neglected, which converts this hyperbolic equation into an elliptic, PV invertibility principle. When the forcing does not vary slowly in time, the second order time derivative term must be retained and the Klein-Gordon equation maintains its hyperbolic character, thereby describing both the balanced flow and the generation and propagation of inertia-gravity waves. These basic concepts are presented in the simple context of the axisymmetric shallow water model (Section 2), where curvature effects and gradient balance come into play. In the linear case (Section 3), the primary results are the PV solution, given in Eq. 8, and the linear Klein-Gordon equation for v (r, t), given in the middle entry of Eq. 11, or the linear invertibility principle for vb (r, t), given by the middle entry of Eq. 12. These linear results are generalized to the nonlinear case through the use of a Lagrangian description of the flow (Section 4). In the nonlinear case, the primary results are the PV solution, given in Eq. 21, and the nonlinear Klein-Gordon equation for V (R, τ), given by Eq. 24, or the nonlinear invertibility principle for Vb(R, τ), given by Eq. 26. Solutions of the nonlinear problem illustrate how rapid intensification can occur when a large percentage of the original lower tropospheric mass inside a given absolute angular momentum surface is diabatically removed via eyewall convection.

2 Shallow water model

Consider the inviscid, axisymmetric, nonlinear shallow water equations

DuDtf+vrv+ghr=0,DvDt+f+vru=0,DhDt+hrurr=hS,(1)

where u (r, t) is the radial velocity component, v (r, t) the azimuthal velocity component, h (r, t) the fluid depth, g the acceleration of gravity, f the constant Coriolis parameter, (D/Dt) = (/∂t) + u (/∂r) the material derivative, and S (r, t) the specified mass sink. The potential vorticity principle associated with Eq. 1 is

DPDt=PSwhereP=f+ζh̄handζ=rvrr,(2)

with the constant h̄ denoting the far-field fluid depth. Our goal in Section 4 is to solve Eq. 1 as a forced, nonlinear, transient, gradient adjustment problem. The approach is to solve analytically for P and then use this solution as a forcing in the second order partial differential equation for v. However, the derivation of a single nonlinear partial differential equation for v is not a straightforward task. For example, one approach might be as follows. First, take (/∂r) of hP=h̄(f+ζ), and use the result to eliminate (∂h/∂r) from the radial momentum equation. Then, write the azimuthal momentum equation in an angular momentum form. The momentum equations then take the form

DuDtf+vrv+gh̄Prrvrr=ghPPrand1rDrvDt+fu=0,(3)

which we can regard as a coupled pair of equations for u (r, t) and v (r, t) if P (r, t) is regarded as known from the solution of the PV principle (Eq. 2). Taking D/Dt of the second entry in Eq. 3 and using the result to eliminate Du/Dt from the first entry, we obtain

rrvrrf̂2ghvPfgh̄DDt1rDrvDt=hh̄Pr,wheref̂2gh=1ghf+vrf+rvrr=1gh̄f+vrP.(4)

It is tempting to regard Eq. 4 as a nonlinear Klein-Gordon equation for the single variable v (r, t), with a time dependent forcing term that is proportional to the radial derivative of the potential vorticity. However, it is important to note that the variables u (r, t) and h (r, t) have not been completely eliminated from Eq. 4, with u (r, t) appearing in the D/Dt operator. In Section 4 we discuss the transformation of Eqs 24 to a Lagrangian coordinate, which simplifies the material derivative D/Dt to a local derivative. This not only simplifies the (D/Dt)-terms in Eqs 24, but also allows for time integration of the potential vorticity equation (Eq. 2), thereby making the forcing term in Eq. 4 a known function. If the forcing term varies slowly enough in time, the term in Eq. 4 involving D/Dt becomes negligible, so that the diagnostic equation for the gradient balanced wind vb (r, t) is

rrvbrrf̂b2ghvb=hh̄Pr,wheref̂b2gh=1ghf+vbrf+rvbrr(5)

which is the shallow water version of the nonlinear invertibility relation discussed by Hoskins et al. (1985). However, before discussing the nonlinear problem, we consider in Section 3 the linearized version of the shallow water dynamics (Eq. 1). Although it is limited to weak vortices, the linear analysis of Section 3 provides a useful guide to the nonlinear analysis of Section 4.

Although the analysis presented here is in the context of the shallow water equations, a more general interpretation is that these equations approximately describe the dynamics of an atmospheric layer confined between two isentropic surfaces. The analogy is as follows. When there is a mass sink in the shallow water continuity equation, the two isentropic surfaces bounding the layer should be regarded as lower tropospheric surfaces (e.g., with the layer defined by 320 ≤ θ ≤ 325K). In contrast, when there is a mass source, the two isentropic surfaces should be regarded as upper tropospheric surfaces (e.g., with the layer defined by 355 ≤ θ ≤ 360K). This interpretation is justified by noting that the shallow water continuity equation in the independent variables (r, t) and the continuity equation for a stratified hydrostatic model in the independent variables (r, θ, t) can respectively be written as

ht+ruhrr=hSandσt+ruσrr=σσθ̇σθ,(6)

where σ = −(1/g) (∂p/∂θ) is the pseudodensity, p is the pressure, θ is the potential temperature, and σθ̇ is the diabatic mass flux. Because of the close correspondence of the two equations in Eq. 6, we can make the following analogies: hσ and S[(σθ̇)/σθ]. Note that in regions where the fractional variation of σ with respect to θ is much less than the fractional variation of θ̇ with respect to θ, this analogy simplifies to S(θ̇/θ). In the lower tropospheric core region of a tropical cyclone, [(σθ̇)/σθ]>0, so there should be an analogous mass sink in the shallow water equations to simulate a lower tropospheric layer. In the upper tropospheric core region of a tropical cyclone, [(σθ̇)/σθ]<0, so there should be an analogous mass source in the shallow water equations to simulate an upper tropospheric layer. In Section 4 and Supplementary Appendix C we develop an alternative analogy based on potential vorticity dynamics rather than simply the mass conservation principles in Eq. 6. Although these analogies help interpret the shallow water model results in terms of the dynamics of a continuously stratified fluid, the analysis of the detailed vertical structure of tropical cyclones requires a complete θ-coordinate model, in which case the Klein-Gordon equations and the invertibility principles discussed in Section 3 and Section 4 would include terms that involve second order partial derivatives in θ.

3 Linear dynamics

As a guide to the nonlinear analysis of Section 4, let us first linearize the system Eq. 1 about a state of rest with constant mean depth h̄. The linear equations are

utfv+ghr=0,vt+fu=0,ht+h̄rurr=h̄S,(7)

where h(r,t)=h(r,t)h̄ is the deviation depth. As initial conditions for Eq. 7 we assume there is no flow and no anomaly in the fluid depth, i.e., u (r, 0) = 0, v (r, 0) = 0, and h′(r, 0) = 0. With these assumptions, any flow that develops is due to the mass sink S (r, t). The equation for the potential vorticity anomaly P′(r, t), derived from the second and third entries in Eq. 7, and its solution are given by

Pt=fSwhereP=Pf=rvrrfh̄hPr,t=f0tSr,tdt.(8)

In order to better understand the stationary and oscillatory solutions of the unforced version of Eq. 7, assume the solutions are separable in r and t, with the time dependent part given by eiωt, where ω is the frequency of oscillation. It can then be shown that the radially dependent part of h′ satisfies an order zero Bessel equation and the radially dependent parts of u and v satisfy order one Bessel equations (for further discussion, see Schubert et al., 1980). In other words, we can simply begin with the assumed separable forms u(r,t)=cûJ1(kr)eiωt, v(r,t)=cv̂J1(kr)eiωt, and h(r,t)=(c2/g)ĥJ0(kr)eiωt, where û,v̂,ĥ are dimensionless complex constants, k is the radial wavenumber, J0 and J1 are Bessel functions of order zero and order one respectively, and c=(gh̄)1/2 is the pure gravity wave speed. Substituting these assumed solutions into the unforced version of Eq. 7 and using the Bessel function derivative relations dJ0 (kr)/dr = −kJ1 (kr) and d [rJ1 (kr)]/rdr = kJ0 (kr), we obtain

iωfckfiω0ck0iωûv̂ĥ=0ωω2f2c2k2=0ω=ω0=0ω1=f2+c2k21/2ω2=f2+c2k21/2.(9)

The three eigenvectors corresponding to the three eigenvalues ω0, ω1, ω2 are

ûv̂ĥ0=A0ckf,ûv̂ĥ1=A2iω1fck,ûv̂ĥ2=A2iω2fck,(10)

where A=f2+c2k21/2 is the normalization constant. Using the assumed form of solution and the derivative relation for rJ1 (kr), the potential vorticity anomaly can be written as P=(ckv̂fĥ)J0(kr)eiωt. For the eigenvector corresponding to the eigenvalue ω0, the potential vorticity anomaly is P=f2+c2k21/2J0(kr). In contrast, since the factor ckv̂fĥ vanishes for both the eigenvector corresponding to the eigenvalue ω1 and the eigenvector corresponding to the eigenvalue ω2, both of these inertia-gravity modes have P′ = 0. In other words, although the inertia-gravity waves contribute to the u (r, t), v (r, t), and h′(r, t) fields, these waves are invisible in the P′(r, t) field. The remainder of this section is devoted to understanding how inertia-gravity waves can be excited on the edge of a rapidly intensifying PV disk or the edges of a rapidly intensifying PV annulus.

Single, independent, second order partial differential equations for u (r, t), v (r, t), and h′(r, t) can be found. The easiest to derive is the equation for the divergent component u (r, t), which can be found by taking /∂t of the radial momentum equation and then using the other two equations to eliminate ∂v/∂t and ∂h′/∂t. Then, including boundary and initial conditions, the complete linear problem for the radial flow is given in the first entry of Eq. 11. Since the rotational wind component v (r, t) plays such a central role in tropical cyclone dynamics, the single partial differential equation for v (r, t) is of particular interest. This equation can be obtained by taking (/∂t) of the v-equation in Eq. 7, then using the u-equation to eliminate (∂u/∂t), and finally using the second entry in Eq. 8 to eliminate h′. Then, including boundary and initial conditions, the complete linear problem for the azimuthal flow is given in the second entry of Eq. 11. Similarly, and again making use of the second entry in Eq. 8, we can derive the single partial differential equation for h′(r, t), which is given in the third entry of Eq. 11. To summarize, the linear primitive equation problems for u (r, t), v (r, t), and h′(r, t) as independent equations are as follows:

Linear Klein-Gordon Equations forur,t,vr,t,hr,t̲rrurr1c2f2+2t2u=SrwithBCs:u=0atr=0,ru0asrICs:u=0andut=0att=0,rrvrr1c2f2+2t2v=PrwithBCs:v=0atr=0,rv0asrICs:v=0andvt=0att=0,rrrhr1c2f2+2t2h=fgP+1gStwithBCs:hr=0atr=0,rh0asrICs:h=0andht=h̄Satt=0.(11)

Equations of the type Eq. 11 are often referred to as linear, inhomogeneous, Klein-Gordon equations (Whitham, 1974), or, in the context of electromagnetic theory, as equations of telegraphy (see Cahn Jr, 1945, and sections 4.10 and 5.6 of Copson, 1975). See Supplementary Appendix A for historical notes on the Klein-Gordon equation. The initial conditions listed in Eq. 11 guarantee that the evolving flow results entirely from the forcing S (r, t). Note that the three equations in Eq. 11 contain time-dependent inhomogeneities that depend on the forcing. For the idealized forcing given below in the first entry of Eq. 13, P′(r, t) grows as shown in the second entry of Eq. 13, so that |∂P′/∂r| at the PV edge increases with time and inertia-gravity waves can be excited there. Thus, more generally, the inner and outer edges of a PV annulus can be regions for the generation of inertia-gravity waves. Since the three Klein-Gordon equations in Eq. 11 are linear, they can be solved by transform methods, e.g., Laplace transforms in time or Hankel transforms in radius. The Hankel transform solution of the Klein-Gordon equation for v (r, t) is discussed in Supplementary Appendix B.

When the forcing is slow enough, the pressure field and the azimuthal wind field remain close to a state of balance, so that the radial wind equation could be discarded and replaced by vb=(g/f)(hb/r). This is equivalent to saying that the (2/∂t2) terms in Eq. 11 all become negligible. Also, when the forcing is slow enough, the time derivative of the forcing becomes small and can be neglected. Thus, in this slow forcing case, each of the equations in Eq. 11 simplifies to a diagnostic equation, with the first equation in Eq. 11 reducing to a diagnostic equation for the slowly evolving divergent flow ub (r, t), the second equation in Eq. 11 reducing to an invertibility principle for the balanced rotational flow vb (r, t), and the third equation in Eq. 11 reducing to an invertibility principle for the balanced mass field hb(r,t). The only initial condition required is on P′(r, t), so that the hyperbolic problems in Eq. 11 reduce to

Linear Balanced Problems forubr,t,vbr,t,hbr,t̲rrubrrμ2ub=Srwith BCs:ub=0atr=0,rub0asr,rrvbrrμ2vb=Prwith BCs:vb=0atr=0,rvb0asr,rrrhbrμ2hb=fgPwith BCs:hbr=0atr=0,rhb0asr,(12)

where μ = f/c is the inverse of the Rossby length.

In order to provide examples of the solutions of the linear balanced flow problems (Eq. 12), we now assume that the mass sink S (r, t) vanishes for r > r0 and is horizontally uniform for rr0, where r0 is a specified constant. The time dependence of the mass sink is assumed to be (t/ts2)et/ts, where the specified constant ts is the time at which the mass sink is a maximum. Small values of ts correspond to rapid forcing and large values of ts to slow forcing, but the total forcing is independent of ts since 0(t/ts2)et/tsdt=1. This forcing S (r, t) and the result of using it in Eq. 8 are given by

Sr,t=Vπr02h̄t/ts2et/ts0rr00r0<r<Pr,t=fVπr02h̄11+t/tset/ts0rr00r0<r<,(13)

where V is a specified constant. Thus, the P′ field is a disk of fixed radius r0, with the value of P′ in its interior increasing in time. The ultimate value of P′ in the disk is independent of the specified constant ts since P(r,t)f(V/πr02h̄) as t for rr0. With this particular forcing S (r, t), integration of the linearized continuity equation over all r and all t yields 2π0[h̄h(r,)]rdr=V, so the constant V can be interpreted as the total volume of fluid removed by the mass sink over its life cycle. We now assume that V=πr02h̄, i.e., the mass sink removes all the fluid initially inside r = r0. Because the linear dynamics have the limitation that |P′| < f, this value of V is essentially the maximum allowable value, since, according to the second entry in Eq. 13, it allows the magnitude of the growing PV jump at r = r0 to maximize at the value f. The time dependence of the mass sink S (r, t) for rr0 is shown in Figure 1A for the three choices ts = 3, 6, 12 h. The time dependence of P′(r, t)/f for rr0, as given by the second entry in Eq. 13, is shown in Figure 1B. Thus, without restating the boundary conditions, the linear balanced problems (Eq. 12) with this forcing (Eq. 13) take the forms

Linear Balanced Problems for the Specified Forcing(13)̲r22ubr2+rubrμ2r2+1ub=0forrr0withubr0r0+=0andrubrrr0r0+=t/ts2et/ts,r22vbr2+rvbrμ2r2+1vb=0forrr0withvbr0r0+=0andrvbrrr0r0+=f11+t/tset/ts,r22hbr2+rhbrμ2r2hb=f2r2g11+t/tset/tsif0r<r00ifr0<r<,withhbr0r0+=0andhbrr0r0+=0.(14)

FIGURE 1
www.frontiersin.org

FIGURE 1. Panel (A) shows the time dependence of the mass sink S (r, t) for rr0, as given by the first entry in Eq. 13, with V=πr02h̄, for the three choices ts =3,6,12 h, where the specified constant ts is the time at which the mass sink is a maximum. Since the areas under the three curves are identical, the total mass removed by the mass sink is independent of ts. Panel (B) shows the corresponding time dependence of P′(r, t)/f for rr0, as given by the second entry in Eq.13, again for the three choices ts =3,6,12 h. Although the rate of increase of P′(r, t) depends on ts, the ultimate value of the P′(r, t) is independent of ts, with P′(r, t)→ f as t for rr0.

Note that the relations Eq. 13 for S (r, t) and P′(r, t) result in Dirac delta functions on the right-hand sides of the ub and vb equations in Eq. 12. In Eq. 14, these become jump conditions on the radial derivatives of rub and rvb. The equations governing ub (r, t) and vb (r, t) for rr0 are homogeneous modified Bessel equations of order one, while the equation governing hb(r,t) is an inhomogeneous modified Bessel equation of order zero. This difference in orders arises from the difference in the radial derivative operators in Eq. 12. The jump conditions on the radial derivatives of rub and rvb can be derived via integration of the first two equations in Eq. 12 over a narrow radial interval surrounding r = r0.

As can be checked by direct substitution, the solutions of Eq. 14, along with the associated formula for the relative vorticity, are

Solutions of the Linear Balanced Problems(14)̲ubr,t=r0t/ts2et/tsI1μrK1μr0if0rr0I1μr0K1μrifr0r<,vbr,t=r0f11+t/tset/tsI1μrK1μr0if0rr0I1μr0K1μrifr0r<,hbr,t=h̄11+t/tset/ts1μr0I0μrK1μr0if0rr0μr0I1μr0K0μrifr0r<,ζbr,t=f11+t/tset/tsμr0I0μrK1μr0if0rr0μr0I1μr0K0μrifr0<r<,(15)

where I0 and K0 are the zero order modified Bessel functions, and I1 and K1 are the first order modified Bessel functions. All four of the solution fields in Eq. 15 are continuous across r = r0 except ζb (r, t). The fact that the solutions ub (r, t) and vb (r, t) given in Eq. 15 satisfy the jump conditions can be checked by using the derivative relations d [rI1 (μr)]/rdr = μI0 (μr) and d [rK1 (μr)]/rdr = −μK0 (μr), and the Wronskian I0 (μr0)K1 (μr0) + K0 (μr0)I1 (μr0) = 1/(μr0). Using this Wronskian, it is easily seen that the jump in ζb across r = r0 is the same as the jump in P′. The maximum inflow occurs at r = r0 when t = ts, while the maximum azimuthal flow occurs at r = r0 as t and the minimum fluid depth anomaly occurs at r = 0 as t.

We now examine actual solutions of the linear balanced model and the linear primitive equation model, using the parameters specified in Table 1 and the three forcing cases shown in Figure 1. The solutions of the linear primitive equation model (Eq. 7) are obtained numerically using a second-order, centered finite differencing scheme radially, a fourth-order Runge-Kutta scheme temporally, an outer radial boundary set at r = 5000 km, a 100 m radial resolution, and a 1 s temporal resolution. Select properties of the linear balanced solutions, specifically the values of ub (r0, ts), vb (r0, ), and hb(0,), are given in Table 2. For the rapid forcing case ts = 3 h, the linear balanced solutions ub(r,t),vb(r,t),hb(r,t) and the linear primitive equation solutions u (r, t), v (r, t), h′(r, t) are respectively plotted in the left and right columns of Figure 2. Similarly, Figure 3 shows linear balanced (left column) and linear primitive equation (right column) results for parameter settings identical to those in Figure 2, except that the time scale ts has been increased from 3 h to 6 h. Note from Figures 2, 3 that the differences between the balanced and primitive equation results become smaller as the forcing time scale ts increases from 3 h to 6 h. In fact, we have chosen not to show the ts = 12 h case because the balanced solutions and the primitive equation solutions are so similar for this case. To better understand the nature of the inertia-gravity wave motion excited in the primitive equation model, Figure 4 shows the time behavior of the u, v, h′ fields at r = 200, 400, 600, 800, 1000 km for the ts = 3 h case (left column) and for the ts = 6 h case (right column). The h′(r, t) field most clearly shows the outward-propagating inertia-gravity wave packet. For the ts = 3 h case, that wave packet causes the negative depth anomaly to overshoot before returning to an equilibrium value that is the same as that given by the balanced solution hb(r,). Although individual inertia-gravity wave components can have shorter time scales, the superposition of these individual components into an inertia-gravity wave packet results in an outward-propagating wave pattern that has a time scale (period) of approximately 11 h.

TABLE 1
www.frontiersin.org

TABLE 1. Parameters used for the linear balanced and linear primitive equation solutions.

TABLE 2
www.frontiersin.org

TABLE 2. The maximum radial inflow, maximum azimuthal wind, and final central depth anomaly for the three forcing cases ts =3,6,12 h, as computed from the linear balanced solutions (Eq. 15) using the specified and computed constants listed in Table 1. The maximum radial inflow is proportional to 1/ts, but the weaker inflows last progressively longer. Since the same volume of fluid is removed in each case, the final values vb (r0,) and hb(0,) are the same. Even though all the original mass inside r = r0 has been removed by the mass sink, almost all of this mass has been replaced by the time-integrated inflow at r = r0, so that the final central depth anomaly is only −11.43 m. Detailed plots of the ts =3 h and ts =6 h cases are presented in Figures 2, 3. The ts =12 h case is not shown since the balanced solutions and the primitive equation solutions are so similar.

FIGURE 2
www.frontiersin.org

FIGURE 2. Panels (A-C) show the linear balanced solutions ub(r,t),vb(r,t),hb(r,t) obtained from Eq. 15, while panels (D-F) show the corresponding linear primitive equation (PE) solutions u (r, t), v (r, t), h′(r, t) obtained from Eq. 7 or equivalently from Eq. 11, all for the rapid forcing case ts =3 h. The isoline intervals are 0.5 m s−1 for the radial and azimuthal flow and 2.0 m for the perturbation depth, and the gray regions indicate initial regions where the plotted field is zero. The peak radial inflow is −3.3 m s−1 for the balanced case and −3.8 m s−1 for the primitive equation case, which has a slight overshoot. As t, both vb (r0, t) and v (r0, t) approach 4.8 m s−1, while both hb(0,t) and h′(0, t) approach −11.4 m. Note that in the outer region 800< r <1000 km, the balanced model response to the inner core forcing is instantaneous (action at a distance), while in the primitive equation model there is a 4–5 h delay as the inertia-gravity wave packet propagates outward at 50 m s−1. The pink area in the u (r, t) panel is weak outward radial velocity and is associated with inertia-gravity waves having very small radial wavenumber k, so that the relevant frequency (f2+c2k2)1/2 is close to the inertia frequency ωf (i.e., a period of approximately 34.9 h).

FIGURE 3
www.frontiersin.org

FIGURE 3. As in Figure 2, panels (A-C) show the linear balanced solutions ub(r,t),vb(r,t),hb(r,t) while panels (D-F) show the corresponding linear primitive equation solutions u (r, t), v (r, t), h′(r, t). All parameter settings used in this figure are identical to those used in Figure 2, except that ts has been changed, i.e., Figure 2 corresponds to the rapid forcing case ts =3 h, while the present figure corresponds to the less rapid forcing ts =6 h. Note that, with less rapid forcing, the differences between the balanced and primitive equation solutions are smaller.

FIGURE 4
www.frontiersin.org

FIGURE 4. This figure shows the time dependence of the linear primitive equation solutions for the radial flow u (r, t), the azimuthal flow v (r, t), and the deviation depth h′(r, t) at the selected radii r =200,400,600,800,1000 km for the rapid forcing case ts =3 h (A-C) and for the less rapid forcing case ts =6 h (D-F). Note that these curves are simply vertical cross sections from Figures 2, 3. The outward propagating inertia-gravity wave packet is most readily apparent in the h′(r, t)-field, where it is immediately seen at r =200 km, and is seen, for example, with an approximately 4.5 h delay at r =1000 km. At each radius for the rapid forcing case (ts =3 h) these linear primitive equation solutions for h′(r, t) overshoot before coming into adjustment at a value that equals hb(r,). The less rapid forcing case (ts =6 h) generates a weaker inertia-gravity wave packet that does not result in such an overshoot in the h′(r, t) solutions.

The fundamental difference between the radial and azimuthal wind fields plotted in Figures 2, 3 can be understood as follows. The Dirac delta function on the right-hand side of the Klein-Gordon equation for u (r, t) and on the right-hand side of the diagnostic equation for ub (r, t) is characterized by a single pulse in time, so that both u (r, t) and ub (r, t) go to zero as t, similar to S (r, t). In contrast, the Dirac delta function on the right-hand side of the Klein-Gordon equation for v (r, t) and on the right-hand side of the invertibility principle for vb (r, t) is characterized by the same time behavior as the PV, so that both vr,17 and vb (r, t) approach nonzero values as t.

The solution of the invertibility problem for vb (r, t) or the solution of the invertibility problem for hb(r,t) can be viewed as essentially converting the radial distribution of P′ into the radial distributions of ζb and hb. This is a nonlocal process since it involves the solution of a second order differential equation. Since P=ζ(f/h̄)h and since P′ = 0 for r > r0, we obtain ζb=(f/h̄)hb<0 for r > r0. The situation is quite different in the region rr0, where P′ > 0. In this region, ζb > 0 and hb<0, so both the wind field and the mass field contribute to the positive PV anomaly. To better understand how the partitioning between ζb and hb in the core (rr0) depends on vortex size, note that the average vorticity in the core is given by 2vb (r0, t)/r0, which allows us to write

Average Core VorticityCore PV Anomaly=2vbr0,t/r0P0,t=2I1μr0K1μr0=0.960,0.891,0.817,0.746,0.680forμr0=0.2,0.4,0.6,0.8,1.0,(16)

where the second equality follows from the cancellation of the time dependence of vb (r0, t) with the time dependence of P′(0, t). Thus, as long as the radius of the vortex core is less than a Rossby length (i.e., μr0 < 1), most of the PV anomaly in the core is partitioned to the vorticity field rather than to the mass field.

Since P′(r, t) grows as shown in Eq. 8, then |∂P′/∂r| at the PV edge increases with time and inertia-gravity waves can be excited there. Thus, the outer edge of a PV tower can be a region for the spontaneous generation of inertia-gravity waves. The inertia-gravity wave packets initiated at r = 200 km propagate inward and outward at c=(gh̄)1/2=50 m s−1 = 180 km h−1, so that at r = 920 km the radial flow u, the azimuthal flow v, and the free surface height anomaly h′ remain zero for the first 4 hours, immediately after which the height anomaly quickly becomes negative, the radial flow quickly becomes inward, and the azimuthal flow becomes cyclonic. Although the inertia-gravity wave packets are visible in the u (r, t), v (r, t), and h′(r, t) fields, they are invisible in the P′(r, t) field. Even though the balanced solutions ub(r,t),vb(r,t),hb(r,t) filter inertia-gravity waves, they are useful approximations of the primitive equation solutions u (r, t), v (r, t), h′(r, t) even when the time scale of the mass sink is as short as ts = 3 h. Since the balanced dynamics have an elliptic rather than a hyperbolic character, the response at large radius to a mass sink at small radius is immediate, i.e., it can be described as “action at a distance.” In this regard, the following analogy is useful: The balanced model is to the primitive equation model as Newton’s gravitational theory is to general relativity, i.e., in the shallow water primitive equations information cannot radially propagate faster than c=(gh̄)1/2, while in general relativity information cannot propagate faster than the speed of light.

We may use the results of this section as the basis for the following physical picture. Each pulse of eyewall convective mass transport from the lower to the upper troposphere is accompanied by a corresponding pulse in the divergent circulation, which ratchets up the rotational circulation and ratchets down the fluid depth. According to the balanced model, the time dependence of the divergent circulation is in lock-step with the forcing S (r, t); this behavior is due to the mathematically elliptic character of the balanced divergent flow (i.e., action at a distance). Convection in the hurricane eyewall can fluctuate on a variety of time scales, so the secondary circulation is expected to fluctuate in a roughly similar fashion. If the convective fluctuations have a vertical structure similar to the first internal mode of the troposphere, the associated gravity wave speed is approximately 50 m s−1 (180 km h−1), so the information from eyewall convective fluctuations will propagate outward to 1000 km in approximately 4.5 h. For eyewall convective fluctuations having time scales larger than 6 h, the amplitude of the associated inertia-gravity wave packet is small, but for eyewall convective fluctuations with time scales less than 3 h, the amplitude of the associated wave packets may be considerably larger. In a continuously stratified fluid the inertia-gravity wave packets also propagate vertically, an important aspect that is further discussed in Section 5.

Although the linear analysis presented in this section provides a helpful guide to the nonlinear analysis of Section 4, the linearization used in Eq. 7 means that we are not analyzing true gradient adjustment and are restricting applications of the analytical results to situations where ζ < f, whereas real tropical cyclones involve situations in which ζ ∼ 100f. In Section 4 we examine the full nonlinear case.

4 Nonlinear dynamics in the Lagrangian coordinate

We now return to the nonlinear shallow water equations of Section 2, and perform numerical integrations using two types of forcing. The first type uses the mass sink defined by Eq. 13, i.e., the mass sink is confined within a disk of fixed radius (0 ≤ rr0). The second type, defined below in Eq. 27, uses a mass sink that is confined within the disk 0 ≤ RR0, where R is a Lagrangian coordinate that collapses as the vortex intensifies. As we shall see, there are some remarkable differences in the vortices that can be produced by these two types of forcing.

Figure 5 compares the nonlinear (left panels) and linear (right panels) numerical integrations using the identical parameter settings and forcing used to obtain the results shown in the right column of Figure 2. A major difference between the linear and nonlinear cases is the inward shift of the radius of maximum wind in the nonlinear case. In the linear case the azimuthal equation is (∂v/∂t) = −fu, so the radial distributions of (∂v/∂t) and − u are identical. In the nonlinear case the azimuthal equation is (∂v/∂t) = −(f + ζ)u, so the radial distribution of (∂v/∂t) is shifted inward when compared to the radial distribution of − u. During the first 24 h, the radius of maximum azimuthal wind shifts inward from 200 km to about 150 km and the maximum azimuthal wind approaches 5.3 m s−1. Thus, for 24 ≤ t ≤ 48 h, the average relative vorticity inside r = 150 km is approximately 1.4f, somewhat outside the range of validity for linear theory.

FIGURE 5
www.frontiersin.org

FIGURE 5. Panels (D-F) show the linear primitive equation solutions u (r, t), v (r, t), h′(r, t) for the rapid forcing case ts =3 h as shown previously in the right panels of Figure 2. Panels (A-C) show the nonlinear primitive equation solutions u(r,t),v(r,t),h(r,t)h̄ obtained from Eq. 1 for the rapid forcing case ts =3 h obtained using the same parameters and forcing as those used for the linear results. The isoline intervals are identical to those used in the linear case (Figure 2), i.e., 0.5 m s−1 for the radial and azimuthal flow and 2.0 m for the perturbation depth, and the gray regions again indicate initial regions where the plotted field is zero. For the nonlinear results, the peak radial inflow is −3.8 m s−1 and occurs at r =200 km, the central depth anomaly h(0,t)h̄ approaches −15.4 m as t, and the maximum azimuthal flow approaches 5.3 m s−1 at about r =150 km. It follows that, for 24≤ t ≤48 h, the average relative vorticity inside r =150 km is approximately 1.47 f. A major difference between the linear and nonlinear cases is the inward shift of the radius of maximum wind in the nonlinear case. Further discussion of this effect is given in Hack and Schubert (1986).

Before discussing results of the second type of forcing, we note that insight into nonlinear model results can be obtained by transforming from the original independent variables (r, t) to the new independent variables (R, τ), where τ = t and R is defined by 12fR2=12fr2+rv=12fr2+RV, with V = (r/R)v denoting the new dependent variable. Note that in general (/∂τ) ≠ (/∂t) since (/∂τ) implies fixed R while (/∂t) implies fixed r. We can think of R as an absolute angular momentum coordinate with the unit of length, or equivalently, as the potential radius, i. e., the radius to which a fluid particle must be moved, conserving its absolute angular momentum, in order for its azimuthal wind to vanish. Finally, in the present context of inviscid flow, R can also be viewed as a Lagrangian coordinate, i. e., as a fluid particle label. Noting that

1+2vfr12VfR=1,(17)

we can easily confirm that transformations between r, v and R, V can be written in the form

R=1+2vfr1/2rV=1+2vfr1/2vr=12VfR1/2Rv=12VfR1/2V.(18)

Since the relative circulation 2πrv divided by the area πr2 is 2v/r, the factor (1 + 2v/fr) can be interpreted as the absolute circulation per unit area, measured in units of f, or equivalently, as the dimensionless average absolute vorticity inside the radius r. In the core of a hurricane we typically find (1 + 2v/fr) ∼ 100 and (1 − 2V/fR) ∼ (1/100), so that r ∼ (1/10) R and v ∼ 10 V. Because of the conservation of absolute angular momentum, the material derivative is (D/Dt) = (/∂t) + u (/∂r) = (/∂τ). As discussed below, this simplification of the material derivative will allow the PV equation to be directly integrated in time, so that the PV part of the dynamics can be treated analytically. The fact that (D/Dt) = (/∂τ) means that the radial advective effects u (∂u/∂r), u (∂v/∂r), u (∂h/∂r), and u (∂P/∂r) become implicit in the coordinate transformation rR. If the solutions of the nonlinear shallow water equations are obtained in R-space, the effects of the radial advection terms do not appear until the final step of transforming the results back to r-space. Many of the characteristics of rapid intensification appear in this transformation Rr.

To further understand the transformation between v(r) and V(R), it is useful to consider the special case of a Rankine vortex. This example is simple because if V has a Rankine structure in R, then v has a Rankine structure in r, i. e.,

VR=V0R/R0if0RR0R0/RifR0R<vr=v0r/r0if0rr0r0/rifr0r<,(19)

where r0, v0 are related to R0, V0 by

r0=12V0fR01/2R0andv0=12V0fR01/2V0.(20)

Thus, if we specify R0 and V0, then compute r0 and v0 from Eq. 20, we can plot V(R) and v(r) using Eq. 19. Figure 6 shows a family of Rankine vortices with R0 = 300 km and 0 ≤ v0 ≤ 65 m s−1 where isolines of the physical azimuthal wind v in the (r, v0)-plane are shown in the top panel and isolines of the transformed azimuthal wind V in the (R, v0)-plane are shown in the bottom panel. For a particular value of v0, the associated horizontal cross section from the top panel gives v(r) for that particular Rankine vortex and the same horizontal cross section from the bottom panel gives the corresponding V(R) for that same Rankine vortex. Since all the effects of radial advection are implicit in the R-coordinate, these often striking radial advective effects appear in the transformation back from R to r, i.e., in going from the bottom panel to the top panel in Figure 6. Some of the information in Figure 6 can also be presented in tabular form, as shown in Table 3. In constructing Table 3 we have chosen R0 = 300 km and a convenient set of values for v0, and then used Eq. 20 to determine r0 and V0. The determination of r0 is easily accomplished by noting that the two equations in Eq. 20 can be combined into r0=[(α2+1)1/2α]R0, where α = v0/(fR0). Once r0 is thereby found, V0 can be computed from V0 = (r0/R0)v0. Table 3 shows properties for a sequence of Rankine vortices with R0 = 300 km, running from tropical depression (TD), to tropical storm (TS), and on upward to category five hurricane strength.

FIGURE 6
www.frontiersin.org

FIGURE 6. An illustration of the transformation relations Eq. 18 for the special case of a family of Rankine vortices, as defined in Eqs. 19, 20, with R0=300 km and f =5×10–5 s−1, so that 12fR0=7.5 m s−1. Panel (A) shows isolines of the physical azimuthal wind v in the (r, v0)-plane with an isoline interval of 5 m s−1, and panel (B) shows isolines of the transformed azimuthal wind V in the (R, v0)-plane with an isoline interval of 0.5 m s−1. In both (A,B) the range on the left ordinate is 0≤ v0≤65 m s−1. This corresponds to the range 0≤ V0≤7.403 m s−1, as shown on the right ordinate of panel (B), and to the range 0≤ ζ0/f ≤76.10, as shown on the right ordinate of panel (A), where ζ0/f =2v0/(fr0) is the dimensionless core vorticity. Note the striking difference in the behavior of V(R) and v(r). As V (R0) increases from rest to 7.403 m s−1, v (r0) increases from rest to 65 m s−1 and the radius of maximum wind decreases from r0=300 km to r0=34.17 km.

TABLE 3
www.frontiersin.org

TABLE 3. Rankine vortex properties for a variety of storm categories, ranging from tropical depression (TD), to tropical storm (TS), to the five hurricane categories (C1–C5), and computed using Eq. 20 with R0=300 km and f =5×10–5 s−1. A range of values of the maximum physical azimuthal wind v0 is specified in the second column followed by computed values in the next three columns of the radius of maximum wind r0, the maximum transformed azimuthal wind V0, and the quantity 2v0/(fr0), which can be interpreted as the core’s relative vorticity, measured in units of f. Since the validity of the linearized dynamics discussed in Section 3 requires that 2v0/(fr0)<1, seven of the eight storm categories shown here fall outside the scope of linear theory.

Because of the conservation of absolute angular momentum, the material derivative is (D/Dt) = (/∂t) + u (/∂r) = (/∂τ). This allows the potential vorticity equation (Eq. 2) and its solution to be written as

Pτ=PSPR,τ=fexp0τSR,τdτandhR,τ=h̄exp0τSR,τdτ,(21)

where the potential depth h* is related to the potential vorticity by Ph=fh̄. Comparing the second entry in Eq. 21 with the linear version (Eq. 8), note that the appearance of the exponential function in Eq. 21 will lead to much faster growth of the P field in the nonlinear case. In fact, it is the appearance of this exponential function that is the crux of the rapid intensification process. In order to make practical use of this and related theoretical results, Hendricks et al. (2021) have introduced and tested a minimal, axisymmetric, shallow water modeling system (SWAMI) as an aid to forecasting short-term tropical cyclone intensity and wind structure changes.

Now consider the transformation of the nonlinear equations (Eq. 3). Using Eqs. 17, 18, and the definitions for R and V, the following useful relationships can be derived

rr=f+ζfRR,f+rvrrfRVRR=f2=f+2vrf2VR,r2f+vr=R2fVR,andrrrvrr=f+ζf3RRRVRR.(22)

Using Eq. 22 and the definitions for P and h*, the nonlinear equations (Eq. 3) become

uτRr3fVRV+ghfh2rh2RRRVRR=gh2rh2RhRandRrVτ+fu=0.(23)

The single, nonlinear partial differential equation for V (R, τ) can now be obtained by eliminating u from the system in Eq. 23. Taking (r/R) (/∂τ) of the V-equation in Eq. 23, and then making use of the u-equation, we obtain

Nonlinear Klein-Gordon Equation forVR,τ̲RRVRR1ghhRhr2fRr2fVRV+rτ1rVτ=fhhRwithBCs:V0,τ=0,RVR,τ0asRICs:VR,0=0,VτR,0=0,(24)

which is the nonlinear generalization of the middle entry in Eq. 11. The initial conditions and boundary conditions in Eq. 24 result from the assumption that the evolving flow is due entirely to a localized forcing S (R, τ). The factors r and h appearing in Eq. 24 can be expressed in terms of V by

r=R12VfR1/2andh=hffRVRR1.(25)

When the forcing term on the right-hand side of Eq. 24 is slow enough, the azimuthal flow remains close to gradient balance. Then, the r (/∂τ)[(1/r) (∂V/∂τ)] term becomes negligible and Eq. 24 simplifies to the following invertibility principle for the balanced transformed azimuthal wind Vb(R, τ):

Nonlinear Invertibility Problem forVbR,τ̲R22VbR2+RVbRμb2R2+1Vb=fR2hhRwith BCs:Vb0,τ=0,RVbR,τ0asRwhereμbR,τ=fghfVbR1/2f2VbR1fRVbRR.(26)

Once again the h appearing in Eq. 26 can be expressed in terms of V by using Eq. 25.

In analogy with the argument given in Section 3 for the linear case, we now assume that the mass sink S (R, τ) vanishes for R > R0 and is horizontally uniform for RR0, where R0 is a specified constant. The time dependence of S (R, τ) is again assumed to be (τ/τs2)eτ/τs, where τs is a specified constant. Again, small values of τs correspond to fast forcing and large values of τs to slow forcing. This forcing S (R, τ) and the result of using it in Eq. 21 are given by

SR,τ=ln1ϵτ/τs2eτ/τsif0RR00ifR0<R<,(27)

and

hR,τ=h̄1ϵ11+τ/τseτ/τsif0RR01ifR0<R<,(28)

where ϵ=V/(πR02h̄) and, as confirmed below, the specified constant V is the total volume of fluid removed by the mass sink. Note that ϵ is the percentage of the initial fluid volume inside the potential radius R0 that is removed by the mass sink and that − ln (1 − ϵ) ≈ ϵ for ϵ ≪ 1, in which case the forcing for the nonlinear case (Eq. 27) is essentially the same as the forcing for the linear case (Eq. 13). The h* field is a disk of potential radius R0, with values of h* in its interior decreasing with time. To confirm that the constant V is the total volume of fluid removed, integrate the mass continuity equation over the entire spatial domain to obtain

t2π0h̄hr,trdr=2π0Shrdr=2π0R0ShRdR=πR02h0,ττ,(29)

where we have made use of hr dr = h*R dR, the assumption that S = 0 for R > R0, and the assumption that S (R, τ) and h*(R, τ) are independent of R for RR0. The last equality in Eq. 29 results from writing the first entry of Eq. 21 in the form (∂h*/∂τ) = −Sh*. The total volume of fluid removed is now obtained by integrating Eq. 29 over time, giving

2π0h̄hr,rdr=πR02h̄h0,=πR02h̄ϵ=V,(30)

where the second to last equality in Eq. 30 results from the use of Eq. 28 to show that h(0,)=(1ϵ)h̄. This confirms that V is the total volume of fluid removed by the mass sink over its life cycle. In other words, for a fixed value of V in Eq. 27, the total volume of fluid removed by the mass sink is independent of the choice of the constants τs and R0, although the choice of R0 should satisfy πR02h̄>V so that 1 − ϵ > 0. Note that since R is a Lagrangian coordinate, fluid does not flow across R surfaces, so ϵ is constrained to be less than unity. As we shall now discuss, a crucial aspect of the forcing (Eq. 27) is its confinement to RR0, i.e., in physical space the forcing collapses as the angular momentum surface R0 collapses. This confines the forcing to the high PV region (Musgrave et al., 2012) and leads to rapid increases in the core PV.

Since P/f=h̄/h, the solution in Eq. 28 can also be expressed in the potential vorticity form

PR,τ=f1ϵ1+τ/τseτ/τs1if0RR01ifR0<R<,P0,=f1ϵ=f1VπR02h̄1,(31)

A plot of P (0, )/f as a function of V/(πR02h̄) is given in Figure 7. As discussed in Supplementary Appendix D, the solution Eq. 31 can be used to obtain a simple upper bound on the rotational flow v (r, t). If we regard R0 as fixed, the abscissa of Figure 7 can be interpreted as a measure of the total volume of fluid removed by the mass sink. If the total volume of fluid removed is V=0.5(πR02h̄), then P (0, ) = 2f. If the total volume of fluid removed is V=0.90(πR02h̄), then P (0, ) = 10f. In the extreme case when the total volume of fluid removed is V=0.99(πR02h̄), then P (0, ) = 100f. In other words, it is the removal of the last vestiges of the original mass inside R0 that spins up the PV to such large values. This effect can also be seen in Figure 8, which shows plots of P (0, τ)/f, as given by Eq. 31, for the time range 0 ≤ τ ≤ 72 h and for the fixed volume parameter V=3.174×1013 m3. The nine curves correspond to combinations of the three choices τs = 3, 6, 12 h and the three choices R0 = 200, 300, 400 km, the latter of which give ϵ = 0.990, 0.440, 0.248. The nine curves illustrate the sensitivity of the core PV to the value of R0 for fixed V. They also illustrate the concepts of long incubation times and vortex preconditioning (see section 15 of Ooyama, 1969). For example, although all nine cases in Figure 8 have the same value of V, the R0 = 400 km case has the mass sink spread over the largest area, so the vortex experiences a long period of incubation, with the PV increasing to only 1.3f by the end of the mass sink. The R0 = 300 km case also experiences long incubation, with the PV increasing to 1.7f by the end of the mass sink. In contrast, the PV in the R0 = 200 km case increases to 100f in 72 h when developing from a resting initial state. However, if the R0 = 200 km case were initialized with P (0, 0) = 5f, it could reach a core PV of 100f in less than 24 h. Such a vortex might be said to be preconditioned for rapid development. The most detailed observations and analyses of a rapidly intensifying hurricane are those of Hurricane Patricia, which occurred in the northeast Pacific in October 2015. The potential vorticity analysis of Martinez et al. (2019) indicates that Hurricane Patricia could be roughly interpreted as the preconditioned R0 = 200 km case, although any detailed comparison would involve the important roles of frictional boundary layer pumping and moist physical processes.

FIGURE 7
www.frontiersin.org

FIGURE 7. The final, dimensionless core potential vorticity P (0,)/f as a function of ϵ=V/(πR02h̄), as given by Eq. 31 and shown by the displayed formula. Note that the scale on the ordinate is logarithmic. If we regard R0 as given, the abscissa is a measure of the total volume of fluid removed by the mass sink. This graph implies that it is the removal of the last 10% of the initial volume of fluid inside R0 that can rapidly spin up the PV from 10f to more than 100f. Since P=(f+ζ)(h̄/h), large values of P are associated with large values of the absolute vorticity (f + ζ) and small values of the fluid depth h.

FIGURE 8
www.frontiersin.org

FIGURE 8. Plots of P (0, τ)/f, as given by Eq. 31 for the three choices τs =3,6,12 h, and for the three choices R0=200,300,400 km. The time evolution of P (0, τ)/f is plotted for 0≤ τ ≤72 h, and the volume parameter is chosen to be V=3.174×1013 m3. This specified value of V is 99.0% of the initial volume inside the potential radius R0=200 km, 44.0% of the initial volume inside the potential radius R0=300 km, and 24.8% of the initial volume inside the potential radius R0=400 km.

As a final illustration of the remarkable intensification that can result with the second type of forcing, Figure 9 shows the results of two nonlinear primitive equation model integrations using the forcing in Eq. 27 with ts = 3 h. The right three panels show the time evolution of u (r, t), v (r, t), and h(r,t)h̄ for the parameter settings R0 = 300 km and ϵ = 0.44, while the left three panels show the corresponding results for R0 = 200 km and ϵ = 0.99. Since V=πR02h̄ϵ, the same total volume of fluid is removed in each case. However, the vortex intensification is quite different in the two cases, with the ϵ = 0.99 case intensifying to hurricane strength and the ϵ = 0.44 case intensifying to only about 4 m s−1. This is consistent with the sensitivity of PV evolution illustrated in Figures 7, 8.

FIGURE 9
www.frontiersin.org

FIGURE 9. Panels (A–F) show the nonlinear primitive equation solutions obtained from Eq. 1 for the radial flow u (r, t), the azimuthal flow v (r, t), and the deviation depth h(r,t)h̄ for the forcing in Eq. 27 with ts =3 h. Panels (A-C) use the parameters R0=200 km and ϵ =0.99, while panels (D-F) use R0=300 km and ϵ =0.44. Since V=πR02h̄ϵ, the same volume of fluid is removed in each case. For the ϵ =0.44 case shown in (D-F), the maximum azimuthal winds are only about 4 m s−1, while for the ϵ =0.99 case shown in (A-C), the maximum azimuthal winds increase to hurricane strength. The gray regions indicate initial regions where the plotted field is zero.

5 Concluding remarks

It is important to note that the simple axisymmetric shallow water model framework of this work does not include environmental effects and the frictional boundary layer, both of which certainly play a role in the rapid intensification process. Nevertheless, the nonlinear results presented here provide what we believe is an important part of the puzzle and are consistent with the notion that intensity forecasting is a difficult problem. For example, suppose (quite hypothetically) that we have a tropical cyclone with its core diabatic forcing (convection) confined by the absolute angular momentum surface R = R0. As can be see from Figures 79, if that diabatic forcing is able to remove 99% of the original mass inside that R0 surface over the next 24 h, then that tropical cyclone will undergo rapid intensification, but if the forcing only removes 50% of the original mass, then it will only slightly intensify. This does not mean that larger storms, i.e., storms with larger R0 values, cannot undergo rapid intensification, just that such a larger storm may require a longer incubation period before rapid intensification can occur.

Concerning the frictional boundary layer, consider a region bounded by two R-surfaces and two lower tropospheric θ-surfaces just above the boundary layer, i.e., an (R, θ)-pseudovolume. The removal of mass from such a pseudovolume involves both diabatic processes and boundary layer pumping. Boundary layer pumping may play a dual role, helping to keep the diabatic heating radially confined so that it collapses with the R-surfaces, but also limiting the percentage of mass that can be removed from that (R, θ)-pseudovolume. We might conclude that accurate determination of the location and magnitude of both frictional and diabatic mass fluxes is necessary for forecasting rapid intensification. The inclusion of environmental effects, such as vertical shear, would require fully three-dimensional arguments and thus a more detailed PV analysis or a generalized omega equation analysis (DeMaria, 1996; Wang and Zhang, 2003), which is an important topic for future research. While the present paper certainly does not solve the rapid intensification problem, perhaps it better defines the nature of the problem.

In closing we note that it can be argued that inertia-gravity waves generated by tropical cyclones are simply “the zero-PV icing on the cake” of what is fundamentally a balanced, potential vorticity phenomenon. Although there is merit to this argument, the observational monitoring of inertia-gravity waves as they propagate upward into the stratosphere and mesosphere can reveal information about eyewall convection. Some of this information might be of practical use in the prediction of rapid intensification. Inertia-gravity waves excited by inner core tropical cyclone convection can sometimes be observed at mesosphere levels by airglow observations from the Earth’s surface or from low Earth orbit (Miller et al., 2012, 2015). For example, Suzuki et al. (2013) described inertia-gravity wave patterns in the mesopause region caused by tropospheric convection in Typhoon Pongsona on 10 December 2002. Concentric rings of inertia-gravity waves in the OH airglow were observed simultaneously by all-sky imagers in a Japanese network of surface stations located between 31.0N and 43.5N. Such ground-based optical measurements require clear-sky conditions, so must be obtained at some distance from the tropical cyclone core in order to avoid persistent cirrus overcast. Such concentric mesospheric airglow waves can also be observed from space (Yue et al., 2014). As a recent example, Figure 10 shows a nocturnal Suomi-NPP VIIRS Day/Night Band image that captured concentric mesospheric airglow waves propagating upward from Super Typhoon Surigae near the time of its peak intensity on 17 April 2021. Concerning evidence from data at stratospheric levels, signatures of stratospheric inertia-gravity waves that have large vertical wavelengths can be found in Atmospheric Infrared Sounder (AIRS) radiance measurements in the 4.3 μm CO2 waveband. Such measurements detect the wave-induced perturbations in stratospheric temperature at levels between 30 and 40 km. Using such data, Hoffmann et al. (2018) searched for an association of stratospheric inertia-gravity waves with the intensification of tropical cyclones. They found a statistical correlation between stratospheric inertia-gravity wave activity and tropical cyclone intensification, lending support to the idea that variability in inertia-gravity wave signatures as a result of changing storm behavior may constitute a useful diagnostic tool for rapid intensification. This interesting possibility led Tratt et al. (2018) to propose a satellite mission concept for monitoring stratospheric inertia-gravity waves from geostationary orbit, which would open up the possibility of continuous monitoring. Although the theoretical results presented here are limited to shallow water dynamics, they do highlight some of the subtleties that are involved in linking inertia-gravity wave activity to rapid intensification. Thus, as research on this topic evolves with the use of full physics numerical models, it would not be surprising to encounter subtleties in the relationship of stratospheric inertia-gravity wave activity to rapid intensification.

FIGURE 10
www.frontiersin.org

FIGURE 10. A nocturnal (several hours before sunrise) Suomi-NPP VIIRS Day/Night Band image (0.7 μm) showing concentric mesospheric airglow waves that have propagated upward from the convective core of Super Typhoon Surigae near the time of its peak intensity on 17 April 2021, when it was just east of the Philippines. From Bachmeier (2021).

An important result presented here is the nonlinear, inhomogeneous, Klein-Gordon equation (Eq. 24), which describes essential aspects of the rapid intensification of the potential vorticity and azimuthal wind fields in a tropical cyclone. Although this equation has an unfamiliar form, it may be of interest to note a relevant point made by Weinberg (1977, page 131) concerning results derived from sound physical and mathematical principles: “This is often the way it is in physics—our mistake is not that we take our theories too seriously, but that we do not take them seriously enough. It is always hard to realize that these numbers and equations we play with at our desks have something to do with the real world.”

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author contributions

WS researched and derived the scientific aspects of this work and wrote the majority of this paper. RT developed and ran the numerical models, created all of the figures, and provided editing and LaTeX support.

Funding

Our research has been supported by the National Science Foundation under grant AGS-1841326.

Acknowledgments

We thank Chris Slocum and Mark DeMaria for their helpful comments and Eric Hendricks, Jonathan Vigh, and Chris Rozoff for sharing with us their ideas on how PV concepts associated with a forced, balanced, axisymmetric shallow water model can be used for understanding short-term tropical cyclone intensification and wind structure changes. We would also like to thank the two reviewers for their useful comments and suggestions.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2022.1038351/full#supplementary-material

References

Bachmeier, S. (2021). Rapid intensification of super Typhoon Surigae. Wisconsin: University of Wisconsin. CIMSS Satellite Blog. Available at: https://cimss.ssec.wisc.edu/satellite-blog/archives/date/2021/04/16 (Accessed 04 25, 2022).

Google Scholar

Cahn, A. (1945). An investigation of the free oscillations of a simple current system. J. Meteor. 2, 113–119. doi:10.1175/1520-0469(1945)002<0113:AIOTFO>2.0.CO;2

CrossRef Full Text | Google Scholar

Copson, E. T. (1975). Partial differential equations. Cambridge: Cambridge University Press, 280.

Google Scholar

DeMaria, M., Sampson, C. R., Knaff, J. A., and Musgrave, K. D. (2014). Is tropical cyclone intensity guidance improving? Bull. Am. Meteorol. Soc. 95, 387–398. doi:10.1175/BAMS-D-12-00240.1

CrossRef Full Text | Google Scholar

DeMaria, M. (1996). The effect of vertical shear on tropical cyclone intensity change. J. Atmos. Sci. 53, 2076–2088. doi:10.1175/1520-0469(1996)053<2076:TEOVSO>2.0.CO;2

CrossRef Full Text | Google Scholar

Doyle, J. D., and Ferek, R. (2017). Tropical cyclone intensity experiment (TCI). AMS Special Collection. Boston: American Meteorological Society. Available at: https://journals.ametsoc.org/collection/TCI (Accessed 10 07, 2022).

Google Scholar

Hack, J. J., and Schubert, W. H. (1986). Nonlinear response of atmospheric vortices to heating by organized cumulus convection. J. Atmos. Sci. 43, 1559–1573. doi:10.1175/1520-0469(1986)043<1559:NROAVT>2.0.CO;2

CrossRef Full Text | Google Scholar

Hendricks, E. A., Vigh, J. L., and Rozoff, C. M. (2021). Forced, balanced, axisymmetric shallow water model for understanding short-term tropical cyclone intensity and wind structure changes. Atmosphere 12, 1308. doi:10.3390/atmos12101308

CrossRef Full Text | Google Scholar

Hoffmann, L., Wu, X., and Alexander, M. J. (2018). Satellite observations of stratospheric gravity waves associated with the intensification of tropical cyclones. Geophys. Res. Lett. 45, 1692–1700. doi:10.1002/2017GL076123

CrossRef Full Text | Google Scholar

Hoskins, B. J., McIntyre, M. E., and Robertson, A. W. (1985). On the use and significance of isentropic potential vorticity maps. Q. J. R. Meteorol. Soc. 111, 877–946. doi:10.1002/qj.49711147002

CrossRef Full Text | Google Scholar

Martinez, J., Bell, M. M., Rogers, R. F., and Doyle, J. D. (2019). Axisymmetric potential vorticity evolution of Hurricane Patricia (2015). J. Atmos. Sci. 76, 2043–2063. doi:10.1175/JAS-D-18-0373.1

CrossRef Full Text | Google Scholar

Miller, S. D., Mills, S. P., Elvidge, C. D., Lindsey, D. T., Lee, T. F., and Hawkins, J. D. (2012). Suomi satellite brings to light a unique frontier of nighttime environmental sensing capabilities. Proc. Natl. Acad. Sci. U. S. A. 109, 15706–15711. doi:10.1073/pnas.1207034109

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, S. D., Straka, W. C., Yue, J., Smith, S. M., Alexander, M. J., Hoffmann, L., et al. (2015). Upper atmospheric gravity wave details revealed in nightglow satellite imagery. Proc. Natl. Acad. Sci. U. S. A. 112, E6728–E6735. doi:10.1073/pnas.1508084112

PubMed Abstract | CrossRef Full Text | Google Scholar

Möller, J. D., and Smith, R. K. (1994). The development of potential vorticity in a hurricane-like vortex. Q. J. R. Meteorol. Soc. 120, 1255–1265. doi:10.1002/qj.49712051907

CrossRef Full Text | Google Scholar

Musgrave, K. D., Taft, R. K., Vigh, J. L., McNoldy, B. D., and Schubert, W. H. (2012). Time evolution of the intensity and size of tropical cyclones. J. Adv. Model. Earth Syst. 4, 15. doi:10.1029/2011MS000104

CrossRef Full Text | Google Scholar

Ooyama, K. (1969). Numerical simulation of the life cycle of tropical cyclones. J. Atmos. Sci. 26, 32–40. doi:10.1175/1520-0469(1969)026<0003:NSOTLC>2.0

CrossRef Full Text | Google Scholar

Schubert, W. H., and Alworth, B. T. (1987). Evolution of potential vorticity in tropical cyclones. Q. J. R. Meteorol. Soc. 113, 147–162. doi:10.1002/qj.49711347509

CrossRef Full Text | Google Scholar

Schubert, W. H., Hack, J. J., Silva Dias, P. L., and Fulton, S. R. (1980). Geostrophic adjustment in an axisymmetric vortex. J. Atmos. Sci. 37, 1464–1484. doi:10.1175/1520-0469(1980)037<1464:GAIAAV>2.0.CO;2

CrossRef Full Text | Google Scholar

Suzuki, S., Vadas, S. L., Shiokawa, K., Otsuka, Y., Kawamura, S., and Murayama, Y. (2013). Typhoon-induced concentric airglow structures in the mesopause region. Geophys. Res. Lett. 40, 5983–5987. doi:10.1002/2013GL058087

CrossRef Full Text | Google Scholar

Tratt, D. M., Hackwell, J. A., Valant-Spaight, B. L., Walterscheid, R. L., Gelinas, L. J., Hecht, J. H., et al. (2018). Ghost: A satellite mission concept for persistent monitoring of stratospheric gravity waves induced by severe storms. Bull. Am. Meteorol. Soc. 99, 1813–1828. doi:10.1175/BAMS-D-17-0064.1

CrossRef Full Text | Google Scholar

Wang, X., and Zhang, D.-L. (2003). Potential vorticity diagnosis of a simulated hurricane. Part I: Formulation and quasi-balanced flow. J. Atmos. Sci. 60, 1593–1607. doi:10.1175/2999.1

CrossRef Full Text | Google Scholar

Weinberg, S. (1977). The first three minutes: A modern view of the origin of the universe. New York: Basic Books, 188.

Google Scholar

Whitham, G. B. (1974). Linear and nonlinear waves. New York: John Wiley & Sons, 636.

Google Scholar

Yue, J., Miller, S. D., Hoffmann, L., and Straka, W. C. (2014). Stratospheric and mesospheric concentric gravity waves over tropical cyclone Mahasen: Joint AIRS and VIIRS satellite observations. J. Atmos. Sol. Terr. Phys. 119, 83–90. doi:10.1016/j.jastp.2014.07.003

CrossRef Full Text | Google Scholar

Keywords: tropical cyclone, rapid intensification, inertia-gravity waves, potential vorticity, axisymmetric shallow water model, Klein-Gordon equation

Citation: Schubert WH and Taft RK (2022) Tropical cyclone rapid intensification and the excitation of inertia-gravity waves on the edges of an evolving potential vorticity structure. Front. Earth Sci. 10:1038351. doi: 10.3389/feart.2022.1038351

Received: 06 September 2022; Accepted: 11 October 2022;
Published: 08 November 2022.

Edited by:

Eric Hendricks, National Center for Atmospheric Research (UCAR), United States

Reviewed by:

Zifeng Yu, China Meteorological Administration, China
Binod Sreenivasan, Indian Institute of Science (IISc), India

Copyright © 2022 Schubert and Taft. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Richard K. Taft, rick.taft@colostate.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.